光流场配准matlab实现

该文件提供了双线性插值、求解光流场方程、图像的梯度等相关M文件，实现了运动图像的光流场配准，配准效果精确。 同时，提供了两幅医学CT图像用于配准实验。

function G=chazhi2（im1uv）  %双线性插值计算灰度值
G=zeros（256256）;
for x=1:256
  for y=1:256
xx = x+u（xy）;
yy = y+v（xy）;
if （xx/double（uint8（xx）） == 1.0） &（yy/double（uint8（yy）） == 1.0）
       G（xy） = im1（uint8（xx）uint8（yy））;
else if （xx<1.0）|（yy<1.0）|（xx>=255）|（yy>=255）
       G（xy）=im1（xy）;
     else                                             % a or b is not integer
     a = double（uint8（xx））;                           
     b = double（uint8（yy））;
     x11 = double（im1（ab））;                            % x11 <- T（ab）
     x12 = double（im1（ab+1））;                          % x12 <- T（ab+1）
     x21 = double（im1（a+1b））;                          % x21 <- T（a+1b）
     x22 = double（im1（a+1b+1））;                        % x22 <- T（a+1b+1）
     G（xy） = uint8（ （b+1-yy） * （（xx-a）*x21 + （a+1-xx）*x11） + （yy-b） * （（xx-a）*x22 +（a+1-xx） * x12） ）; % calculate G（xy）
    end
end
  end
end 

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

     文件      66614  2010-09-12 09:03  GJC\1.bmp

     文件      66614  2010-09-12 09:03  GJC\2.bmp

     文件      66614  2011-05-20 15:14  GJC\3.bmp

     文件      66614  2011-05-20 15:14  GJC\4.bmp

     文件        923  2011-05-12 10:35  GJC\chazhi2.m

     文件        862  2011-05-11 21:33  GJC\ComputeDerivatives.m

     文件        482  2011-05-12 10:41  GJC\guangliu2.m

     文件        196  2011-05-10 21:56  GJC\msefunction.m

     文件        419  2011-05-12 11:10  GJC\tuxiangpeizhun.m

     目录          0  2011-05-20 15:14  GJC

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

               269338                    10