360度全景图像展开MATLAB程序

         本程序主要是将360度全景相机拍摄的图片展开为一般图像。MATLAB代码及实验结果如下

function Expand2()
Src=imread('100.jpg');
Src=im2double(Src);
x_dot=1010;      % 全景图像中心点横坐标
y_dot=562;        % 全景图像中心点纵坐标
InnerR=98;        % 全景相机中间黑色圆圈半径
OuterR=497;      % 注意这个外径大小的选取必须都在图像内部  
Width=2*pi*OuterR;   % 展开后的图像的宽度
Height=OuterR-InnerR;  % 展开后的图像的高度
ImgDst=zeros(Height,Height,3);

for j=1:Height
    for i=1:Width
        Angle=GetAngle(i,j);
        Radius=GetRadius(i,j);
        [x,y,TempSrc]=GetRGB(Angle,Radius,Src,InnerR,x_dot,y_dot);
        ImgDst(j,i,1)=TempSrc(1);
        ImgDst(j,i,2)=TempSrc(2);
        ImgDst(j,i,3)=TempSrc(3);
    end
end


figure(2)
imshow(Src);
[m,n]=ginput(1);close(gcf);
Radius2=Radius+InnerR;
Angle2=acos((m-x_dot)/Radius2);
ExpandW2=497*Angle2;
ExpandH2=sqrt((m-x_dot)^2 + (n-y_dot)^2 );
figure(1);
imshow(ImgDst);
% hold on;plot(ExpandW2,ExpandH2,'r*','MarkerSize',15)


function dw_Angle=GetAngle(ExpandWidth,~)
dw_Angle=ExpandWidth/497;          % 若主函数中ExpandWidth的大小变了 497也得变  变得准则是ExpandWidth/497的pi值的2倍即可



function Radius=GetRadius(~,ExpandHeight)
Radius=ExpandHeight;

function [x,y,TempSrc]=GetRGB(dw_Angle,Radius,Src,InnerR,x_dot,y_dot)
Radius=Radius+InnerR;
x=round(Radius*cos(dw_Angle)+x_dot);
y=round(Radius*sin(dw_Angle)+y_dot);

if x<=0 x=1;end
if y<=0 y=1;end
if y>=size(Src,1) y=size(Src,1); end
if x>=size(Src,2) x=size(Src,2); end

TempSrc(1)=Src(y,x,1);
TempSrc(2)=Src(y,x,2);
TempSrc(3)=Src(y,x,3);

图示说明：r1为全景中的黑色圆圈的半径，对应程序中的InnerR，r2或者r3为第二个黑圈的半径长度，在此图示中，由于r2对应的黑圈出了图像区域，因此以r2作为程序中的OuterR，若r2和r3都没有出边界，随便选就行。

结果如下：

说明：上面为原始全景图像，下方为对其展开后的图像。