matlab图像畸变校正

你自己看这一段

for i=1:h %从理想图像矩阵出发处理

for j=1:w

x=[1,j-og(1),i-og(2),(j-og(1))^2,(i-og(2))*(j-og(1)),(i-og(2))^2]

u=x*a0+og(2) % 逆向映射（j，i）到畸变图像矩阵（v，u）

v=x*b0+og(1)

if (u>1)&&(u<w)&&(v>1)&&(v<h) %处理在图像大小范围内的像素野搜弯点

uu=floor(u) %对u取整

vv=floor(v) %对v取整

arf=u-uu%计算上面提到的

bta=v-vv%计算上面提到的

for k=1:3 %进行灰度双线性插值

ft1=(1-bta)*b(vv,uu,k)+bta*b(vv+1,uu,k)

ft2=(1-bta)*b(vv,uu+1,k)+bta*b(vv+1,uu+1,k)

sp(i,j,k)=(1-arf)*ft1+arf*ft2

end

end

imshow(uint8(sp)) %显示校正图像

3个for加一个if，可只有2个end，程序都没写完当然出错。

——颂闷——————————————————————漏中——————————

然后还有这一段

for k=1:n%转换到以对称点为原点的空间关系并构造矩阵A

A(k)=[1,gm(k,1)-og(1),gm(k,2)-og(2),(gm(k,1)-og(1)^2), (gm(k,1)-og(1))*(gm(k,2)-og(2)),(gm(k,2)-og(2) ^2)]

end

A(k)是一个元素，可你却把它定义为一个数组，肯定也要出错。

————————————————————————————————

你要我改，我只能保证程序能运行，但对不对我无法保证。

function gmodify(pic,uv,gm,og) %pic表示要处理的图像的路径文件名

%uv是一个二维矩阵,uv(:,1)代表上面提到的,uv(:,2)表示

%gm是一个二维矩阵,gm(j,:)代表在校正图空间上与uv(j,:)一一应的点

%og 代表对称中心，它是一个二维向量

a=imread(pic)

b=double(a)

n=size(gm(:,1))

for k=1:n%转换到以对称点为原点的空间关系并构造矩阵A

A(k,:)=[1,gm(k,1)-og(1),gm(k,2)-og(2),(gm(k,1)-og(1)^2), (gm(k,1)-og(1))*(gm(k,2)-og(2)),(gm(k,2)-og(2) ^2)]

end

[h,w]=size(b(:,:,1))

sp=zeros(h,w,3)+255

a0=pinv(A)* uv(:,2) %计算上面提到的地址映射的系数估计a

b0=pinv(A)* uv(:,1) %计算上面中提到的地址映射的系数估计b

for i=1:h %从理想图像矩阵出发处理

for j=1:w

x=[1,j-og(1),i-og(2),(j-og(1))^2,(i-og(2))*(j-og(1)),(i-og(2))^2]

u=x*a0+og(2) % 逆向映射（j，i）到畸变图像矩阵（v，u）

v=x*b0+og(1)

if (u>1)&&(u<w)&&(v>1)&&(v<h) %处理在图像大小范围内的像素点

uu=floor(u) %对u取整

vv=floor(v) %对v取整

arf=u-uu%计算上面提到的

bta=v-vv%计算上面提到的

for k=1:3 %进行灰度双线性插值

ft1=(1-bta)*b(vv,uu,k)+bta*b(vv+1,uu,k)

ft2=(1-bta)*b(vv,uu+1,k)+bta*b(vv+1,uu+1,k)

sp(i,j,k)=(1-arf)*ft1+arf*ft2

end

end

end

end

imshow(uint8(sp)) %显示校正图像

可以用RODAN法校正， 1. 先灰度化,求出春亩图像尺寸[m n] = size(I)2. theta = 0:180 R = radon(I, theta) 3. 求出图像中心点至边界的距离L = round((m/2).^2+(n/2).^2) [C,angle] = max(R(L,:)) angle=angle-1%angle为图像倾斜角度 4. A = imrotate(I,angle,'bilinear') imshow(A)%校正扒改森歼团后图片

---