应用经纬映射的鱼眼图像校正设计方法

来源：互联网发布：微信请检查网络设置编辑：程序博客网时间：2024/05/16 11:28

理论是根据《应用经纬映射的鱼眼图像校正设计方法》作者由畸变鱼眼图上像素点(i，j)得到了经纬度坐标

然后我根据这个编出来效果不对出来是个全黑屏的figure 然后我又看了另几篇关于经纬度校正的论文不只是像那个论文作者写的那样有了经纬度坐标就行了她省略了后面的一步可能她以为我们都知道原来还要将得到的经纬度坐标再转换成图像平面坐标这次是按照正向映射编的只是为了看有没校正效果如果有校正效果再反一下采用反向映射编就好了

原图：

提取有效区域后的图：

经纬度校正法效果图：

可以看到我正向映射而且没有插值的有空隙点所以需要改成反向映射然后图左右反了所以应该程序哪里弄反了不过这已经比我前一段时间弄出的全黑屏好多了至少我感觉还是有希望的

A=imread('F:\fisheye\s2.jpg');
[A,R]=kuaisusaomiao(A,40);
[m,n,k]=size(A);
for i=1:m
for j=1:n
u=j-R;
v=R-i;
r=sqrt(u^2+v^2);
if(r==0)
fi=0;
elseif(u>=0)
fi=asin(v/r);
else
fi=pi-asin(v/r);
end
f=R*2/pi;
theta=r/f;
x=f*sin(theta)*cos(fi);
y=f*sin(theta)*sin(fi);
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
xx=round(f*sita);
yy=round(f*fai);
if(xx<1||yy<1||xx>m||yy>n)
continue;
end
C(xx,yy,1)=A(i,j,1);
C(xx,yy,2)=A(i,j,2);
C(xx,yy,3)=A(i,j,3);
end
end

刚刚在这基础上改了下结果：

的确是将左右重新反过来了可是可以看到正向映射是多对一比如下面这一块本来是校正成黑色的物体结果原图房子的像素以及某块天空的像素也映射在这黑色物体上了造成了叠加不行

刚刚又改了改上个改的不行这次改的才对结果：

这就比较好了只有待插值的空隙点没有之前的叠加现象

然后我用的是cubic插值结果如下：

初步插值的效果还算行因为我是一行一行插值的所以像被一个东西横向碾压过一样我再修改修改

调整了一个阈值的大小现在这样好点了下面黑色物体这里插值痕迹比较明显

又修改了下结果：

这个已经比较完美了吧

经纬度算法校正还是比较快的一下就出来图了但是插值相对较慢要等一下才弹出图来

好了全部编成M文件现在再拿几幅图试试经纬度校正的效果：

A=imread('F:\fisheye\s3.jpg');
>> [A,R]=kuaisusaomiao(A,40);
>> figure
>> imshow(A)
>> C=jingwei(A,R);
>> figure
>> imshow(C)

原图和校正图：

再试一个：原图和效果图：

这个的阈值我其实应该再调节小一点下面插值的痕迹太明显了

这个经纬映射校正总体还可以但并不是对所有的都非常适合就像一种插值方法不可能适用于所有的图像一样这个经纬度校正以及cubic插值也不是对所有图像都有很好的校正效果比如下面的：

原图和效果图

可以看到还是有变形这可能是不适用于这幅图的原因也可能是因为这幅图采用的鱼眼相机的视场角、焦距和我假设的不一样我编程时是假设视场角为180 焦距f=r/theta 还假设相机采用等距投影如果这幅图的相机不是按照我假设的这样那么较正效果成这个样子那我也能理解毕竟我们学校太穷了买不起鱼眼相机不可能用来做实验即使买得起也不是给我用。。。