matlab demons算法用来实现图像的配准

demons算法用来实现图像的配准，特别是医学图像的配准，已经达到了很高的精度

function [fx fy ft] = computeDerivatives（im1 im2）

if size（im21）==0
    im2=zeros（size（im1））;
end

% Horn-Schunck original method
fx = conv2（im10.25* [-1 1; -1 1]‘same‘） + conv2（im2 0.25*[-1 1; -1 1]‘same‘）;
fy = conv2（im1 0.25*[-1 -1; 1 1] ‘same‘） + conv2（im2 0.25*[-1 -1; 1 1] ‘same‘）;
ft = conv2（im1 0.25*ones（2）‘same‘） + conv2（im2 -0.25*ones（2）‘same‘）;

% derivatives as in Barron
% fx= conv2（im1（1/12）*[-1 8 0 -8 1]‘same‘）;
% fy= conv2（im1（1/12）*[-1 8 0 -8 1]‘‘same‘）;
% ft = conv2（im1 0.25*ones（2）‘same‘） + conv2（im2 -0.25*ones（2）‘same‘）;
% fx=-fx;fy=-fy;

% An alternative way to compute the spatiotemporal derivatives is to use simple finite difference masks.
% fx = conv2（im1[1 -1]）;
% fy = conv2（im1[1; -1]）;
% ft= im2-im1;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-11-11 14:50  光流法matlab程序\
     文件         759  2011-11-11 09:15  光流法matlab程序\computeDerivatives.m
     文件      526106  2011-10-12 09:48  光流法matlab程序\CT0.dcm
     文件      526110  2011-10-12 09:48  光流法matlab程序\CT5.dcm
     文件         921  2011-11-11 09:15  光流法matlab程序\ffd.m
     文件        1080  2011-11-11 09:15  光流法matlab程序\gaussFilter.m
     文件        2818  2011-11-11 09:15  光流法matlab程序\HS.asv
     文件        2788  2011-11-11 14:40  光流法matlab程序\HS.m
     文件      236079  2011-11-11 09:15  光流法matlab程序\HSsnapshot.JPG
     文件       34678  2011-11-11 09:15  光流法matlab程序\lenag1.png
     文件       36184  2011-11-11 09:15  光流法matlab程序\lenag2.png
     文件        1382  2011-11-11 09:15  光流法matlab程序\plotFlow.asv
     文件        1382  2011-11-11 09:15  光流法matlab程序\plotFlow.m
     文件         944  2011-11-11 09:15  光流法matlab程序\smoothImg.m
     文件       78886  2011-11-11 09:15  光流法matlab程序\yos10.tif
     文件       78908  2011-11-11 09:15  光流法matlab程序\yos9.tif