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