matlab仿真的图像序列的柱面全景拼接算法 （Panorama stitching algorithm）

matlab仿真的图像序列的柱面全景拼接算法，其中有柱面投影程序可对其它投影公式有一定的启示。

function [T1coor_shift02]=coortransf（T0fcoor_shift）
%%%%图像坐标系到柱面坐标的转换
%对输入图像序列T0按焦距f进行变换
coor_shift02=coor_shift;%%%%第一维（行值）不变，第二维（列）在映射后更新
[HWrN]=size（T0）;
w2f=W/2/f;
h2=H/2;
constant2=f*atan（W/（2*f））;
constant1=h2;

for y=1:W      %%%%%%   列
    angle=atan（y/f-w2f）;%%%%atan（（y-W/2）/f）;
     y1=uint16（f*angle+constant2）;
     if  y1==0  
         y1=1; 
     end
    for x=1:H  %%%%%%%%%%%%   行
     x1=uint16（（x-h2）*cos（angle）+constant1）;
     if  x1==0 
         x1=1;
     end
     if r==3   %%%%%%%%%%%%%彩色图像
         for n=1:N  %%%%%%%%%
             if （y==coor_shift（n2）） 
                 coor_shift02（n2）=y1; 
             end   %%%对应偏移量
             T1（x1y1:n）=T0（xy:n）;%%%%点的映射
         end
     elseif r==1
%          T1（x1y1）=T0（xy）;
     end
    end
end
[hwaN]=size（T1）;
for i=1:60%%%边缘没有映射点的地方填白色
    for j=1:w
        if （T1（ij::）==0）
            T1（ij::）=255;
        end
        if （T1（h-ij::）==0）
            T1（ij::）=255;
        end
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-07-05 20:56  Panorama stitching algorithm\
     文件        1125  2014-07-03 11:47  Panorama stitching algorithm\coortransf.m
     文件        1481  2006-10-24 15:59  Panorama stitching algorithm\cylinpano2images.m
     文件         331  2006-06-27 16:55  Panorama stitching algorithm\imagefusion02.m
     文件         672  2006-11-14 16:24  Panorama stitching algorithm\inorm.m
     文件        2149  2014-07-03 11:27  Panorama stitching algorithm\main.m
     文件         619  2006-11-14 16:56  Panorama stitching algorithm\medfilt.m
     文件        1065  2006-10-23 21:42  Panorama stitching algorithm\mosaic.m
     文件         660  2006-10-24 15:22  Panorama stitching algorithm\multi_resolution.m
     文件        1537  2006-06-30 17:48  Panorama stitching algorithm\phase_correlation.m
     文件        2612  2006-10-24 19:53  Panorama stitching algorithm\poc_2pow.m