% 试验结果参考 相册中的图片
%只是用于镜头视角<pi/2的摄像头采集的图像的畸变校正。
%function DistortionCorrect();
clear all;
clc;
tic
pi=3.141592653589793;
I0=imread('pic.bmp');
I1=rgb2gray(I0);
[height0,width0]=size(I1);
height=2*height0;
width =2*width0;
I2=zeros(height,width);
a1=height/4;
a2=height*3/4;
b1=width/4;
b2=width*3/4;
for i=1:height
for j=1:width
if i<=a1 || i>a2 || j<=b1 || j>b2;
I2(i,j)=255;
else
I2(i,j)=I1(i-a1,j-b1);
end
end
end
Inew=DistortionRate(I2);
% I2=uint8(I2);
Inew=uint8(Inew);
%figure,imshow(I2);
figure,imshow(Inew);
toc
time=toc
%%=======================&