《单幅圆形鱼眼图像的校正》matlab程序

声明：本文转载至http://blog.csdn.net/wd1603926823/article/details/48374369。转载至此，方便以后查阅

理论部分来自《单幅圆形鱼眼图像的校正》 其实他是在《应用经纬映射的鱼眼图像校正设计方法》得到经纬度后  再采用柱面压缩 把经纬度变换到一个图像平面即可

这篇其实我早就看了  只是之前经纬度校正程序编出来是一个乱七八糟的 现在既然把经纬度校正程序弄正确了 那就可以编这篇了  其实就最后一步不一样而已

function C=zongxiangyasuo(A,R)
%《单幅圆形鱼眼图像的校正》
[m,n,k]=size(A);
C=zeros(m,n);
f=2*R/pi;
dx=pi/m;
dy=dx;
for i=1:m
      for j=1:n
         u=i-R;
         v=R-j;
         ro=sqrt(u^2+v^2);
         if(ro>R)
            continue;
         end
         if(ro==0)
            phi=0;
         elseif(u>0)
             phi=asin(v/ro);
         else
             phi=pi-asin(v/ro);
         end
         theta=ro/f;
         x=f*sin(theta)*cos(phi);
         y=f*sin(theta)*sin(phi);
         z=f*cos(theta);
         rr=sqrt(x^2+z^2);
         sita=pi/2-atan(y/rr);
         if(z>=0)
            fai=acos(x/rr);
         else
            fai=pi-acos(x/rr);
         end
         vv=round(sita/dy);
         if(fai>pi/2)
            uu=2+sqrt(2-2*sin(fai))/sin(pi/2-atan((1-sin(fai))/abs(cos(fai))));
         else
            uu=2-sqrt(2-2*sin(fai))/sin(pi/2-atan((1-sin(fai))/abs(cos(fai))));
         end
         uu=2*R-round(uu/dx);
         if(uu<1||vv<1||uu>m||vv>n)
              continue;
         end
         C(vv,uu,1)=A(j,i,1);
         C(vv,uu,2)=A(j,i,2);
         C(vv,uu,3)=A(j,i,3);
     end
end
C=uint8(C);
T=15;
Y=[];
Z1=[];
Z2=[];
Z3=[];
for i=1:m
    for j=1:n
         I(i,j)=0.59*C(i,j,1)+0.11*C(i,j,2)+0.3*C(i,j,3);
         if(I(i,j)>=T)
            Y=[Y j];
            Z1=[Z1 C(i,j,1)];
            Z2=[Z2 C(i,j,2)];
            Z3=[Z3 C(i,j,3)];
         end
    end
    if(size(Y)<3)
          Y=[];
          Z1=[];
          Z2=[];
          Z3=[];
          continue;
    else
       for j=1:n
          I(i,j)=0.59*C(i,j,1)+0.11*C(i,j,2)+0.3*C(i,j,3);
           if(I(i,j)<T)
             C(i,j,1)=interp1(Y,double(Z1),j,'cubic');
             C(i,j,2)=interp1(Y,double(Z2),j,'cubic');
             C(i,j,3)=interp1(Y,double(Z3),j,'cubic');
           end
       end
    end
    Y=[];
    Z1=[];
    Z2=[];
    Z3=[];
end

结果：

可以看到最后一步处理与经纬度校正不一样  导致了出来结果和上一篇有点不一样 即 图片好像横向拉伸了  之前经纬度我用原图大小可以正好装下 不会丧失原图信息   可是这个柱面压缩校正不行  原图大小不能够容下它  我再改大点试试  也不能设定不能超过原图大小的条件看来      结果如下：

可以看到柱面压缩校正的确像拉伸了一样  校正图比原图大  右边这些黑的 是因为我刚刚只改了校正部分的大小  插值那里没改  所以只把原图大小的地方插值了   我再改下  结果：

这就是王志周童鞋这篇论文中用柱面压缩校正的效果  

经纬度校正和柱面压缩校正对比下：