［求助］如何将傅立叶变换提取的形状参数还原回原来的图像？急～～～
clear all;close all;
pic=imread('111.JPG');
pic=imresize(pic,0.1);
g=pic(:,:,1);
g=wiener2(g,[9 9]);
T=graythresh(g);
InputImageBin=im2bw(g,T);
%figure(1),imshow(InputImageBin);
pic4=edge(InputImageBin,'sobel');
%figure(2),imshow(pic4),title('sobel 算子');
se90=strel('line',4,90);
se0=strel('line',4,0);
BWsdil=imdilate(pic4,[se90 se0]);
%figure(3),imshow(BWsdil),title('膨胀操作');
bw=imfill(BWsdil,'holes');
[m,n]=size(bw);
bound = boundaries(bw); %生成边界
bound = bound{1};
z=frdescp(bound); %计算傅利叶描述子
zc=fftshift(z);
plot(abs(zc)) %频谱图
feature = abs(zc([1:13]))'
c=ifft(feature);
figure(4),imshow(log(abs(c)),[]);
title('傅利叶描述子图');
请问各位高手们，如何将
提取的形状参数还原回原来的图像．红色的程序哪里错了，怎么改？小弟在此先谢过了．
