PSF参数估计

模糊长度和尺度估计

I=imread（‘tire.tif‘）;
I=im2double（I）;
len=20;
theta=75;
PSF=fspecial（‘motion‘lentheta）;
I1=imfilter（IPSF‘circular‘）;
I2=fft2（I1）;
I2=abs（I2）;
figure
imshow（I2）;
%对图像做平滑滤波
I3=fftshift（I2）;
[mn]=size（I3）;
d0=80;
m1=fix（m/2）;
n1=fix（n/2）;
for i=1:m
    for j=1:n
        d=sqrt（（i-m1）^2+（j-n1）^2）;
        h（ij）=exp（-d^2/2/d0^2）;
    end
end
I4=I3.*h;
figure
imshow（I4）;
I4=histeq（I4512）; %直方图均衡
I4=im2bw（I4）;      %二值化
I4=edge（double（I4）‘sobel‘）; %边缘提取

%进行Hough变换
[HTR] = hough（I4）;
peaks=houghpeaks（H5）;
r_theta=[];
r_theta=peaks（:2）;
maxmum=max（r_theta）;
R_theta=90-maxmum;
ipsf=fspecial（‘motion‘lenR_theta）;
I5=deconvwnr（I1ipsf）;
figure
imshow（I5）;

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

     文件        746  2009-03-22 10:26  tai_frequency_domain.m

     文件       3480  2009-03-23 22:38  tai_time_domain.m

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

                 4226                    2