MATLAB傅里叶变换轮廓术三维重建

clear all;
clc;




%读取全息图像
fr1=imread（‘0.copy-filt.bmp‘）;
fr2=imread（‘4.5.copy-clear.bmp‘）;

info=imfinfo（‘0.copy-filt‘‘bmp‘）;

w=info.Width; 
h=info.Height; 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5555
%采用傅立叶变换法
fr1=im2double（fr1）;
% fr2=double（fr2）;
FFTfr1=fftshift（fft2（fr1））;
% FFTfr2=fftshift（fft2（fr2））;
%显示log频谱
 tmp1=abs（log（FFTfr1））;
Frequen1=uint8（（tmp1-min（tmp1（:）））/（max（tmp1（:））-min（tmp1（:）））*255）;
%  Frequen1=log（1+abs（FFTfr1））;
 figure（）;
 imshow（Frequen1[]）;

 
%  tmp2=abs（log（FFTfr2））;
%  Frequen2=uint8（（tmp2-min（tmp2（:）））/（max（tmp2（:））-min（tmp2（:）））*255）;
%  figure（2）;
%  imshow（Frequen2[]）;

%%
%零级坐标
Xmax1=1141;
Ymax1=541;
Xmax2=1141;
Ymax2=541;
%%正一级谱坐标
XoneOrder1=1213;
YoneOrder1=541;
Xone