图像配准程序matlab

图像畸变配准的matlab程序，不需在做编译，直接就可以用

%  Basic demon registration code. （To easy understand the algorithm）

% Clean
clc; clear all; close all;

% Compile the mex files
compile_c_files

% Read two images
I1=im2double（imread（‘images/lenag1.png‘））;  
I2=im2double（imread（‘images/lenag2.png‘））; 

% Set static and moving image
S=I2; M=I1;

% Alpha （noise） constant
alpha=2.5;

% Velocity field smoothing kernel
Hsmooth=fspecial（‘gaussian‘[60 60]10）;

% The transformation fields
Tx=zeros（size（M））; Ty=zeros（size（M））;

[SySx] = gradient（S）;
for itt=1:200
	    % Difference image between moving and static image
        Idiff=M-S;

        % Default demon force （Thirion 1998）
        %Ux = -（Idiff.*Sx）./（（Sx.^2+Sy.^2）+Idiff.^2）;
        %Uy = -（Idiff.*Sy）./（（Sx.^2+Sy.^2）+Idiff.^2）;

        % Extended demon force. With forces from the gradients from both
        % moving as static image. （Cachier 1999 He Wang 2005）
        [MyMx] = gradient（M）;
        Ux = -Idiff.*  （（Sx./（（Sx.^2+Sy.^2）+alpha^2*Idiff.^2））+（Mx./（（Mx.^2+My.^2）+alpha^2*Idiff.^2）））;
        Uy = -Idiff.*  （（Sy./（（Sx.^2+Sy.^2）+alpha^2*Idiff.^2））+（My./（（Mx.^2+My.^2）+alpha^2*Idiff.^2）））;
 
        % When divided by zero
        Ux（isnan（Ux））=0; Uy（isnan（Uy））=0;

        % Smooth the transformation field
        Uxs=3*imfilter（UxHsmooth）;
        Uys=3*imfilter（UyHsmooth）;

        % Add the new transformation field to the total transformation field.
        Tx=Tx+Uxs;
        Ty=Ty+Uys;
        M=movepixels（I1TxTy）; 
end

subplot（131） imshow（I1[]）; title（‘image 1‘）;
subplot（132） imshow（I2[]）; title（‘image 2‘）;
subplot（133） imshow（M[]）; title（‘Registered image 1‘）;
	

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       14719  2010-06-03 22:32  demon_registration_version_8f\register_images.m
     目录           0  2016-03-30 11:02  functions_affine\
     目录           0  2016-03-30 10:03  demon_registration_version_8f\
     文件        1690  2016-03-30 10:04  demon_registration_version_8f\basic_demon_example.m
     文件         830  2009-08-14 21:31  demon_registration_version_8f\compile_c_files.m
     目录           0  2016-03-30 10:01  demon_registration_version_8f\functions_nonrigid\
     文件       49611  2009-08-14 22:25  demon_registration_version_8f\functions_nonrigid\image_interpolation.c
     文件        2518  2009-07-15 01:13  demon_registration_version_8f\functions_nonrigid\image_interpolation.h
     文件        1677  2009-07-14 03:35  demon_registration_version_8f\functions_nonrigid\movepixels.m
     文件        8658  2009-07-28 02:40  demon_registration_version_8f\functions_nonrigid\movepixels_2d_double.c
     文件        2500  2009-07-14 02:07  demon_registration_version_8f\functions_nonrigid\movepixels_2d_double.m
     文件       16896  2016-03-30 10:01  demon_registration_version_8f\functions_nonrigid\movepixels_2d_double.mexw64
     文件        6225  2009-07-28 02:40  demon_registration_version_8f\functions_nonrigid\movepixels_3d_double.c
     文件       14848  2016-03-30 10:01  demon_registration_version_8f\functions_nonrigid\movepixels_3d_double.mexw64
     文件        6282  2009-07-28 02:40  demon_registration_version_8f\functions_nonrigid\movepixels_3d_single.c
     文件       14848  2016-03-30 10:01  demon_registration_version_8f\functions_nonrigid\movepixels_3d_single.mexw64
     目录           0  2016-03-30 10:00  demon_registration_version_8f\images\
     文件        1848  2009-08-14 21:31  demon_registration_version_8f\images\checkboard2.png
     文件        7053  2009-08-14 21:31  demon_registration_version_8f\images\checkboard3.png
     文件        7107  2009-08-14 21:31  demon_registration_version_8f\images\checkboard4.png
     文件       34678  2009-08-14 21:31  demon_registration_version_8f\images\lenag1.png
     文件       36184  2009-08-14 21:31  demon_registration_version_8f\images\lenag2.png
     文件       41466  2009-08-14 21:31  demon_registration_version_8f\images\lenag3.png
     文件        3792  2009-08-14 21:31  demon_registration_version_8f\images\modtest1.png
     文件        3481  2009-08-14 21:31  demon_registration_version_8f\images\modtest2.png
     文件        3793  2009-08-14 21:31  demon_registration_version_8f\images\modtest3.png
     文件        1314  2014-02-12 13:08  demon_registration_version_8f\license.txt
     文件        1677  2009-07-14 03:35  demon_registration_version_8f\movepixels.m
     文件        2500  2009-07-14 02:07  demon_registration_version_8f\movepixels_2d_double.m
     目录           0  2016-03-30 11:02  functions_nonrigid\
     文件       13139  2009-10-02 04:00  demon_registration_version_8f\register_volumes.m
............此处省略75个文件信息