图像畸变校正

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)&&(u1)&&(v