%??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個文件信息