matlab练习程序（弧形投影）

这个其实也算是圆柱体投影了，不过上一篇文章是从正面看，得到的是凸形的结果，而这个是从反面看，得到的是凹形的结果。

计算公式就不写了，大致介绍一下，计算公式中关于x坐标求法和上篇一样，y坐标则正好是上篇公式的反变换，结合上篇公式代码和本篇的代码，应该都不是很难理解的。

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

原图：

处理后结果：

matlab代码如下：

clear all; close all;clc;

img=imread('lena.jpg');
[h,w]=size(img);

hfOV=pi/;     %可取区间为(,pi)
f=w/(*tan(hfOV/));

x1=;
x2=floor(*f*atan(w/(*f)));
y1=floor(h/-h*(sqrt((w/)^+f^))/(*f));
y2=floor(h/+h*(sqrt((w/)^+f^))/(*f));

newh=y2-y1;
neww=x2-x1;
imgn=zeros(newh,neww);
for i=+y1:newh+y1
    for j=:neww

       %反变换公式的应用
       x=floor(f*tan(j/f-atan(w/(*f)))+w/);
       y=floor(h/+f*(i-h/)/sqrt(f^+(w/-x)^)); 

        if x>= && x<=w && y>= && y<=h
            imgn(i-y1,j)=img(y,x);
        end       

    end
end

imshow(img);
figure;
imshow(imgn,[]);