圆柱投影就是将一张二维的图像投影到三维的圆柱体上,不过在显示图像的时候依然是以二维的形式给出。

投影最重要的步骤就是计算投影变换公式,和图像旋转类似,只要得到变换公式,再依照公式进行代码编写就很容易了。

这里就不写投影变换公式的推导过程了,直接给出变换公式。公式分为正变换和反变换,编程时,反变换公式通常更有用。

正变换公式如下:

gif.latex?%0D%0Ax%E2%80%98=f%5Ccdot%5Carctan%5Cleft(%7B%5Cfrac%7Bx-%5Cfrac%7BW%7D%7B2%7D%7D%7Bf%7D%7D%5Cright)+f%5Ccdot%5Carctan%5Cleft(%7B%5Cfrac%7BW%7D%7B2f%7D%7D%5Cright)

gif.latex?%0D%0Ay%E2%80%98=%5Cfrac%7Bf%5Ccdot%5Cleft(y-%5Cfrac%7BH%7D%7B2%7D%5Cright)%7D%7B%5Csqrt%7B%5Cleft(x-%5Cfrac%7BW%7D%7B2%7D%5Cright)%5E2+f%5E2%7D%7D+%5Cfrac%7BH%7D%7B2%7D

其中,x,y为原图的坐标,x‘,y‘为变换后图像的坐标,W,H为原图的宽和高,f=W/(2*tan(hfOV/2)),这里hfOV为相机水平视角,我们通过设置这个参数来实现不同的投影效果,hfOV取值为(0,pi)之间。

反变换公式可以通过正变换公式推导出来,其实在程序中已经写出来了。

下面是hfOV为pi/2时得到的变换结果:

原图:

c5967bad3e0c50e6a93ee1a522dadc2f.png

处理后结果:

ecdeba26a9bb5adb40b6cabb35a061e6.png

matlab代码如下:

clear all; close all;clc;

img=imread(‘lena.jpg‘);

[h,w]=size(img);

hfOV=pi/2; %可取区间为(0,pi)

f=w/(2*tan(hfOV/2));

newh=h; %变换后高度不变

neww=floor(f*atan((w-w/2)/f)+f*atan(w/(2*f))); %通过正变换公式计算最大宽度

imgn=zeros(newh,neww);for i=1:newhfor j=1:neww

%反变换公式的应用

x=floor(f*tan(j/f-atan(w/(2*f)))+w/2);

y=floor((i-h/2)*sqrt((x-w/2)^2+f^2)/f+h/2);if x>=1 && x<=w && y>=1 && y<=h

imgn(i,j)=img(y,x);end

end

endimshow(img);

figure;

imshow(imgn,[]);

原文:http://www.cnblogs.com/tiandsp/p/3815119.html

更多推荐