matlab 代码实现的傅里叶描述子

matlab 代码实现的傅里叶描述子，程序简单明了.

function [new_indeces]=find_close_indeces（im_in）

[rc]=find（im_in）;

p=pdist（[r c]）;
psqr=squareform（p）;

 nl=length（r）;
 
new_indeces=[r（1） c（1）];
ind_ind_data=1;
lenindata=1;
newind=1;
 for ind=2:nl
  
    mcur_dist=psqr（newind:）;
    [dist_mindist_min_ind]=sort（mcur_dist）;
    [dmin inddmin]=setdiff（ dist_min_indind_ind_data）;
    dist_min_ind =  dist_min_ind（sort（inddmin））;
    newind=dist_min_ind（1）;
    new_indeces  =[new_indeces; [r（newind） c（newind）]];
    
     ind_ind_data=[ind_ind_data; newind];
    lenindata=length（ind_ind_data）;
    
 end
 

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

     文件        595  2010-02-27 13:01  foirier\find_close_indeces.m

     文件       8780  2013-01-31 10:45  foirier\img_1.jpg

     文件       7624  2013-01-31 10:45  foirier\img_2.jpg

     文件       1340  2010-03-01 12:27  foirier\license.txt

     文件       1337  2013-01-31 14:18  foirier\make_fft_sec.m

     文件        532  2013-01-31 14:21  foirier\test_fourier.asv

     文件        491  2013-01-31 14:21  foirier\test_fourier.m

     目录          0  2013-01-31 11:52  foirier

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

                20699                    8