matlab练习程序(全景图到穹顶图)

  • Post author:
  • Post category:其他


这个程序我最初是用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映射应该才是全景图最常用的处理方法,不过那又是另一类变换了。

转载于:https://www.cnblogs.com/tiandsp/p/3670201.html