傅里叶描述子代码

傅里叶描述子有效描述物体轮廓形状，代码简洁易懂

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  Fourier Descriptor\find_close_indeces.m
     文件        1340  2010-03-01 12:27  Fourier Descriptor\license.txt
     文件        1215  2010-02-28 18:03  Fourier Descriptor\make_fft_sec.m
     文件         814  2016-01-03 17:40  Fourier Descriptor\test_fourier.m