《基于几何成像模型的鱼眼镜头图像校正算法和技术研究》实现
这是根据2013一个期刊上的小论文《基于几何成像模型的鱼眼镜头图像校正算法和技术研究》中的校正原理式(11)和式(12)编写的，其实这两个式子给出的是二维的校正方法，就跟之前的经度坐标校正差不多都是平面校正，所以我不知道这篇论文中给出式(6)干嘛？有什么用？还有这论文里说校正后的图像宽为w，高为h，这两个参数怎么确定呢？在没校正之前我怎么知道校正后的图像是怎样大小的？有谁知道吗，如果有，请告诉我。这里，根据这两个式子编程其实很简单，我令校正和畸变鱼眼图一样大小，即2Rx2R。
function C=jihemoxing(A,R)
%基于几何成像模型的鱼眼镜头图像校正算法和技术研究
w=2*R;
h=2*R;
xo=w/2;
yo=h/2;
f=2*R/pi;
for u=1:w
for v=1:h
ho=sqrt((u-xo)^2+(v-yo)^2);
h1=f*atan2(ho,f);
x=h1*(u-xo)/ho+xo;
y=h1*(v-yo)/ho+yo;
x=round(x);
y=round(y);
C(u,v,1)=A(x,y,1);
C(u,v,2)=A(x,y,2);
C(u,v,3)=A(x,y,3);
end
end
C=uint8(C);
这是M文件 其中A是鱼眼图，R是鱼眼图半径
实验结果如下：
A=imread('F:\orl_zhifangtu\s3.jpg');
[A,R]=kuaisusaomiao(A,40);
C=jihemoxing(A,R);
>> imshow(C)
而原畸变鱼眼图是这样的：
可以看到 用这篇论文的方法校正其实不怎么好 丢失了原图的信息 当然 这可能和我将校正后的图规定为2Rx2R有关 于是 我将M文件里的校正后的图的大小改大 改成了3Rx3R这样 重新试了下
A=imread('F:\orl_zhifangtu\s3.jpg');
[A,R]=kuaisusaomiao(A,40);
>> w=3*R;
>> h=3*R;
>> xo=w/2;
yo=h/2;
f=2*R/pi;
for u=1:w
for v=1:h
ho=sqrt((u-xo)^2+(v-yo)^2);
h1=f*atan2(ho,f);
x=h1*(u-xo)/ho+xo;
y=h1*(v-yo)/ho+yo;
x=round(x);
y=round(y);
if(x>2*R || y>2*R || x<1 || y<1)
continue;
end
C(u,v,1)=A(x,y,1);
C(u,v,2)=A(x,y,2);
C(u,v,3)=A(x,y,3);
end
end
C=uint8(C);
>> imshow(C)
结果：
所以应该不是改预设图像大小的问题 有谁知道那个wxh怎么确定的 告诉我啊
我比较了一下经度坐标校正的效果 经度坐标校正的效果如下：
这样一比较 很明显看到 这篇论文的校正效果不够好 和经度坐标校正比起来差了很多 都是尚未插值的比较的 明显看得出来 感觉是不是这个算法不好
还有一个期刊上发表的《鱼眼图像校正和配准算法研究》这上面的校正原理我觉得没交代清楚 而是直接给出了公式，从目标图像到鱼眼图像反向映射的公式(3)(4)(5)，既然给出来了 那就按照这个直接写哦：
A=imread('F:\orl_zhifangtu\s3.jpg');
[A,R]=kuaisusaomiao(A,40);
[m,n,k]=size(A);
for i=1:m
for j=1:n
theta=i/R;
fi=j/R;
x=R*cos(fi)*sin(theta);
y=R*sin(fi);
z=R*cos(fi)*cos(theta);
u=R*cos(1/tan(y/x))/tan(sqrt(x^2+y^2)/z);
v=R*sin(1/tan(y/x))/tan(sqrt(x^2+y^2)/z);
u=round(u);
v=round(v);
if(u<1||v<1||u>m||v>n)
continue;
end
C(i,j,1)=A(u,v,1);
C(i,j,2)=A(u,v,2);
C(i,j,3)=A(u,v,3);
end
end
>> imshow(C)
结果原图和校正图如下：
这是什么东西 明明按照公式写的 很简单啊 是期刊上那三个公式错了 没原理？？
我觉得我没写错 这个乱七八糟的图 是公式的原因吧