matlab练习程序(全景图到穹顶图)
这个程序我最初是用FreeImage写的,这两天改成了matlab,再不贴上来,我就要忘了。
看到一篇文章有这样的变换,挺有意思的,就拿来试了一下,文章点此。
全景图到穹顶图变换,通俗的说就是将全景图首尾相接做成一个圆环的样子。
先看下面这张图:
下面的矩形就是我们要处理的全景图,上面的矩形是变换后的图像。下面图像的底边对应穹顶图的内圆,顶边对应穹顶图的外圆,当然,反过来也是可以的。
程序流程:
1.定义穹顶图内圆和外圆的半径,变换后的像素就填充在这个内外半径的圆环中。
2.遍历穹顶图,当所处理当前像素位于圆环内,则通过极坐标反变换去全景图中寻找相应位置的像素进行填充。
3.遍历完图像就行了。
用的技巧和图像旋转或放大缩小都是类似的。
处理结果:
原图:
结果:
matlab代码如下:
clear all; close all; clc; img=imread('pan.jpg'); imshow(img); [m,n]=size(img); r1=100; %内环半径 r2=r1+m; %外环半径 imgn=zeros(2*r2,2*r2); [re_m,re_n]=size(imgn); for y=1:re_m for x=1:re_n dis_x=x-re_n/2; dis_y=y-re_m/2; l=sqrt(dis_x^2+dis_y^2); if l<=r2 && l>=r1 theta=0; if y>re_m/2 theta=atan2(dis_y,dis_x); end if y<re_m/2 theta=pi+atan2(-dis_y,-dis_x); end if y==re_m/2 theta=atan2(dis_y,dis_x)+0.0001; end xx=ceil(n*theta/(2*pi)); yy=ceil(l-r1); if yy>=1 && yy<=m && xx>=1 && xx<=n imgn(y,x)=img(yy,xx); end end end end figure; imshow(imgn,[])
最后要说的是,一般我们要是有一张全景图,通常会用cubic映射,将图像变换为立方体的六个面,然后通过图形学方法贴到立方体上,就能做出类似谷歌街景的样子。cubic映射应该才是全景图最常用的处理方法,不过那又是另一类变换了。