傅里叶变换对干涉图进行研磨

  clear 
   T=imread（‘52与标准.png‘）;
  
  T=im2double（T）;
  h=fspecial（‘gaussian‘[7 7]）;
T=imfilter（Th）;
T=medfilt2（T[7 7]）;
T=adapthisteq（T）;
 T0=imresize（T[256 256]‘bilinear‘）;
 T0=im2double（T0）;
 l=1;
f=256;
 x=linspace（-llf）;
y=linspace（-llf）;
for i=1:f
    for j=1:f
      if （x（i）^2+y（j）^2）<（0.8*l）^2
          BW（ij）=1;
      else
          BW（ij）=0;
      end
    end
end
T=T0.*BW;figure;imshow（T）;
edge1=edge（BW‘canny‘）;figure;imshow（edge1）;
figure;imshow（BW）;
a=pi/8;