鱼眼镜头的distortion校正【matlab】

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%         作者：WWC        %%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  功能：畸变矫正
clc;
clear;
close all;
%% 读取图像
A=imread('D:\文件及下载相关\图片\distortion2.jpg');
[m,n,k]=size(A);
R=n/;
w=*R;
h=*R;
xo=m/;
yo=n/;
f=*R/pi;
for u=:w
    for v=:h
        ho=sqrt((u-xo)^+(v-yo)^)+;
        h1=f*atan2(ho,f);
        x=h1*(u-xo)/ho+xo;
        y=h1*(v-yo)/ho+yo;
        x=round(x);
        y=round(y);
        if(x>*R || y>*R || x< || y<)
           continue;
        end
        C(u,v,)=A(x,y,);
        C(u,v,)=A(x,y,);
        C(u,v,)=A(x,y,);
    end
end
C=uint8(C);
imshow(C);

鱼眼镜头畸变校正。（直接处理彩色图）

代码存在错误：貌似只能处理宽比长大的，竖着的。和方的。