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