phi1 = double(rgb2gray(imread('peppers.png')));
xsize = min(大小(phi1));
fftsize = 2。^(nextpow2(xsize));
phi1m = phi1-mean2(phi1);
IMFFT = FFTSHIFT(FFT2(PHI1M,FFTSIZE,FFTSIZE)/ FFTSIZE);
imfftnorm = sqrt(imfft。*结合(imfft))/ max(max(abs(imfft)));
XMP =尺寸(IMFFT,1)/ 2;
YMP =尺寸(IMFFT,2)/ 2;
半径=圆形(FFTSIZE / 3);
imrot = imrootate(imfftnorm,-90);
AMPL =零(1800,1);
为了Ang = 1:1800
xend = xmp round(cosd(-ang / 10)*半径);
YEND = YMP-圆形(SIND( - ang / 10)*半径);
l1 = Improfile(imrot,[xmp xend],[ymp yend]);
Xend2 = XMP +圆形(COSD(-ang / 10)*半径);
YEND2 = YMP +圆形(SIND( - - ang / 10)*半径);
L2 = Improfile(IMROT,[XMP XEND2],[YMP YEND2]);
剧(Ang)= NANSUM(L1(:))+ NANSUM(L2(:));
结尾