傅里叶描述算子

利用matlab对图像的轮廓提取其傅里叶描述算子，其中还有将选取得到的描述子傅里叶反变换进行可视化的验证程序。

function B = bound2im（b M N x0 y0）

[np nc] = size（b）;
if np < nc
   b = b‘; % To convert to size np-by-2.
   [np nc] = size（b）;
end
 
% Make sure the coordinates are integers.
x = round（b（: 1））; 
y = round（b（: 2））;
 
% Set up the default size parameters. 
x = x - min（x） + 1;
y = y - min（y） + 1;
B = false（max（x） max（y））;
C = max（x） - min（x） + 1;
D = max（y） - min（y） + 1;
 
if nargin == 1  
   % Use the preceding default values.
elseif nargin == 3
   if C > M | D > N
      error（‘The boundary is outside the M-by-N region.‘）
   end
   % The image size will be M-by-N. Set up the parameters for this.
   B = false（M N）;
   % Distribute extra rows approx. even between top and bottom.
   NR = round（（M - C）/2）; 
   NC = round（（N - D）/2）; % The same for columns.
   x = x + NR; % Offset the boundary to new position.
   y = y + NC;  
elseif nargin == 5
   if x0 < 0 | y0 < 0
      error（‘x0 and y0 must be positive integers.‘）
   end
   x = x + round（x0） - 1;
   y = y + round（y0） - 1;
   C = C + x0 - 1;
   D = D + y0 - 1;
   if C > M | D > N
      error（‘The shifted boundary is outside the M-by-N region.‘）
   end
   B = false（M N）;
else
   error（‘Incorrect number of inputs.‘） 
end
 
B（sub2ind（size（B） x y）） = true;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4025  2017-12-03 18:19  源程序\a0.png

     文件       2465  2017-12-03 18:19  源程序\b0.png

     文件       1286  2017-11-18 19:07  源程序\bound2im.m

     文件       2563  2017-12-03 18:19  源程序\c0.png

     文件       4999  2017-12-03 10:21  源程序\ceshi1.png

     文件      13954  2017-12-03 10:39  源程序\ceshi2.png

     文件        195  2017-12-02 20:13  源程序\dist.m

     文件        437  2017-12-03 00:55  源程序\frdescp.m

     文件        492  2017-11-19 21:10  源程序\ifrdescp.m

     文件        622  2017-12-03 01:12  源程序\normalfrd.m

     文件       5446  2017-12-03 12:32  源程序\pingyi.png

     文件       8910  2017-12-03 12:32  源程序\shensuo.png

     文件        745  2017-12-03 18:20  源程序\test1.m

     文件        953  2017-12-03 18:25  源程序\test2.m

     文件        345  2017-12-03 18:24  源程序\test3.m

     文件      14451  2017-12-03 12:32  源程序\xuanzhuan.png

     文件       4025  2017-12-03 12:32  源程序\yuantu.png

     目录          0  2017-12-03 18:27  源程序

----------- ---------  ---------- -----  ----

                65913                    18