%**************************************************************************
%該程序為運動模糊參數識別的頻域方法
%模糊方向識別的MRT算法
%**************************************************************************
clear?all
clc

%xrgb=imread（‘0121.jpg‘）;
%x=0.11*xrgb（::1）+0.59*xrgb（::2）+0.3*xrgb（::3）;

x=imread（‘cameraman.tif‘）;
[rowcol]=size（x）;
Blurlength=10;
Blurtheta=10;
for?Bn=1:length（Blurtheta）
??
PSF=fspecial（‘motion‘BlurlengthBlurtheta（Bn））;
f=imfilter（xPSF‘circular‘‘conv‘）;
figure（1）
imshow（f）
FFTf=fft2（im2double（f））;
FrequencyF=log（1+abs（fftshift（FFTf）））;
%FrequencyF=im2bw（FrequencyF0.9999）;
%FrequencyF=abs（fftshift（FFTf））;
figure（2）
imshow（FrequencyF）
%__________________________________________________________________________
%Radon?變換求模糊方向
%________________________

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????3646??2013-05-24?10:14??AngleExtentTest.m