fid=fopen('P02560.7','r+');
status=fseek(fid,61465,'bof');
A=fread(fid,[256,256],'int16');
A=double(A);
I1=A/(max(max((abs(A)))));
figure,imshow(abs(I1));
data = 1/sqrt(length(I1(:)))*ifftshift(ifft2(fftshift(I1)));
figure,imshow(abs(data),[]);
IG=ifftshift(fft2(ifftshift(data)));
figure,imshow(abs(IG),[]);
这个程序出现的第2副图是这样的
问下怎样让图显示完整
status=fseek(fid,61465,'bof');
A=fread(fid,[256,256],'int16');
A=double(A);
I1=A/(max(max((abs(A)))));
figure,imshow(abs(I1));
data = 1/sqrt(length(I1(:)))*ifftshift(ifft2(fftshift(I1)));
figure,imshow(abs(data),[]);
IG=ifftshift(fft2(ifftshift(data)));
figure,imshow(abs(IG),[]);
这个程序出现的第2副图是这样的
