iterative_blind_deconvolution 实现了迭代盲反卷积的图像去模糊

实现了迭代盲反卷积的图像去模糊，自己写的代码，没有调用matlab已有函数

clear all; 
star_rgb = imread（‘C:\Users\freeboy\Desktop\涂扬—去模糊\RobustMotionDeblur\RobustMotionDeblur_CPU\Sample Image\Large\roma.png‘）; 
star_gray = rgb2gray（star_rgb）; 
figure; imshow（star_gray）; imwrite（star_gray‘oldstar.bmp‘）; title（‘原图‘）; 
f = double（star_gray）; 
size_f = size（f）; 
 
m1=min（f（:））; 
m2=max（f（:））; 
l=m2-m1+1; 

for i=m1:m2 
    His（i-m1+1）=length（find（f==i））;%统计处各灰度的像素点数 
end 

A=sum（His）;%一阶动量矩 
Temp=double（m1:m2）; 
M=His.*Temp; 
B=sum（M）;%二阶动量矩 
t=round（（m1+m2）/2）;%初值 
oldt=t-1; 

%用平均值平均法求出阈值 
while（oldt~=t） 
    oldt=t; 
    B1=sum（M（1:t））; 
    A1=sum（His（1:t））; 
    u1=B1/A1; 
    u2=（B-B1）/（A-A1）; 
    t=round（（u1+u2）/2）; 
end 
t=t-1+m1;%输出阈值 
%二值化分割 
k=find（f%以上找到支撑域 
size_psf = 15; %psf大小 
h_z = zeros（size_f）; 
for u= 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1945  2012-07-13 09:14  iterative_blind_deconvolution\BD.m

     文件     249610  2012-07-13 09:28  iterative_blind_deconvolution\best.bmp

     文件      42061  2011-08-20 09:09  iterative_blind_deconvolution\cameraman.jpg

     文件     405574  2007-12-13 21:27  iterative_blind_deconvolution\gray_star.bmp

     文件       2608  2012-07-13 10:27  iterative_blind_deconvolution\IDB.asv

     文件       2651  2012-07-13 10:32  iterative_blind_deconvolution\IDB.m

     文件       4602  2012-07-12 10:44  iterative_blind_deconvolution\IDBCopy.m

     文件       9966  2012-07-12 15:21  iterative_blind_deconvolution\ipexblind.m

     文件     249610  2012-07-13 09:14  iterative_blind_deconvolution\newstar.bmp

     文件     249610  2012-07-13 09:14  iterative_blind_deconvolution\oldstar.bmp

     文件     249610  2012-07-13 09:14  iterative_blind_deconvolution\psf.bmp

    ..A.SH.      5120  2008-09-04 22:23  iterative_blind_deconvolution\Thumbs.db

     目录          0  2012-07-13 09:13  iterative_blind_deconvolution

----------- ---------  ---------- -----  ----

              1473435                    15