PIIFD描述符matlab源代码

《A Partial Intensity Invariant Feature Descriptor for multimodal retinal image registration》论文的matlab源代码

% Auther: Jian Chen
% April 2008
% Department of Biomedical imaging Columbia University New York USA
% Institute of Automation Chinese Academy of Sciences Beijing China
% email: jc3129@columbia.edu  jian.chen@ia.ac.cn
% All rights reserved

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% read images
[p1p2] = rr_getpath;
if ~strcmp（p2pwd）
    path（pathpwd）;
end
 

t = zeros（181）;
for i = 1:1
    close all
    f1 = [num2str（i）‘_1.jpg‘];
    f2 = [num2str（i）‘_2.jpg‘];
[I1 I2 f1 f2 pathname] = rr_readimage;%（f1f2）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fprintf（‘It may take half a minute to registering those two images\n‘）;
fprintf（‘Please be patient...\n‘）;

Iter = 2;
tic;
if Iter == 1
    t_fundus = rr_registration1（I1 I2 ‘linear conformal‘）;
    I1_c = imtransform（I1t_fundus‘XData‘[1 size（I22）] ‘YData‘[1 size（I21）]）;
    [I1_cI2_c] = rr_imagesize（I1_cI2）;
    clear I1
    clear I2
end

k = 1;
while Iter>1
    if k == 1
        t_fundus = rr_registration1（I1 I2 ‘linear conformal‘）;
        I1_c = imtransform（I1t_fundus‘XData‘[1 size（I22）] ‘YData‘[1 size（I21）]）;
        [I1_cI2_c] = rr_imagesize（I1_cI2）;
        clear I1
        clear I2
%         figureimshow（I1_c+I2_c[]）
        Iter = Iter+1;
    else
        t_fundus = rr_registration2（I1_c I2_c ‘affine‘）;
%         t_fundus = registration1（I1_c I2_c  ‘affine‘）;
        I1_c = imtransform（I1_ct_fundus‘XData‘[1 size（I2_c2）] ‘YData‘[1 size（I2_c1）]）;
        [I1_cI2_c] = rr_imagesize（I1_cI2_c）;
        
    end
    Iter = Iter-1;   
    k = 2;
end

t（i） = toc;
figureimshow（I1_c+I2_c[]）
title（‘fusion of results‘）
fprintf（‘Runtime = %6.2f seconds.\n\n‘t（i））;
fprintf（‘The results have been saved in the same directory as your original images\n\n‘）;

%%%%%%%%%%%%%%%%%%%%%%%%%% save images
I3 = I1_c+I2_c;
I3 = I3/max（I3（:））;
f11 = [pathname（1:end-3）‘moved_‘f1];
f22 = [pathname（1:end-3）‘moved_‘f2];
f33 = [pathname（1:end-3）‘fusion_‘f2（1:end-4）‘-‘f1];
% imwrite（I1_cf11）;
% imwrite（I2_cf22）;
% imwrite（I3f33）;
end

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

     文件        548  2008-04-14 14:38  rr_translation.m

     文件        492  2008-04-14 14:46  read me.txt

     文件       2146  2008-11-30 13:52  registration.m

     文件        447  2008-04-14 14:36  rr_appendimages.m

     文件       1449  2008-04-14 14:36  rr_checkerboardshow.m

     文件        941  2008-04-15 14:06  rr_corr.m

     文件        946  2008-04-14 14:36  rr_d2gauss.m

     文件       4209  2008-04-14 14:36  rr_descriptor.m

     文件        565  2008-04-14 14:36  rr_getpath.m

     文件       2161  2008-03-23 01:48  rr_harris.m

     文件        667  2008-04-14 14:36  rr_image_pat.m

     文件        465  2008-04-14 14:37  rr_imagesize.m

     文件       1384  2008-04-14 14:37  rr_match.m

     文件        804  2008-04-14 14:37  rr_msk.m

     文件       2155  2008-04-14 14:37  rr_MTransformation.m

     文件        886  2008-04-14 14:37  rr_orientation.m

     文件       1278  2008-04-15 14:20  rr_para_transform.m

     文件       1374  2008-04-14 14:37  rr_readimage.m

     文件       3200  2008-04-15 13:41  rr_refine.m

     文件       2284  2008-11-30 14:17  rr_reg.m

     文件       2531  2008-11-30 13:55  rr_registration1.m

     文件       2230  2008-04-18 12:16  rr_registration2.m

     文件        988  2008-04-14 14:38  rr_removeboundarypoint.m

     文件       2594  2008-04-14 14:38  rr_removeincorrectmatch.m

     文件       3369  2008-11-30 14:09  rr_removeincorrectmatch_new.asv

     文件       3346  2008-11-30 14:12  rr_removeincorrectmatch_new.m

     文件        935  2008-04-14 14:38  rr_rigidtransformation.m

     文件        780  2008-04-14 14:38  rr_showmatch.m

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

                45174                    28

............此处省略1个文件信息