MATLAB处理向前映射和向后映射

clc;
clear;
close all;

%向后映射法
I1=imread（‘3.2.bmp‘）;
[a b]=size（I1）;
I2=I1;
A2=[989898*981;9815898*1581;158158158*1581;15898158*981];
B2=[11888;90140;139168;168116];
X2=A2\B2;
U=zeros（1a）;
Y=zeros（1b）;
for x=1:a
    for y=1:b
        U（x）=[xyx*y1]*X2（:1）;
        V（y）=[xyx*y1]*X2（:2）;
        v=U（x）-floor（U（x））;
        u=V（y）-floor（V（y））;
        if U（x）<1
            U（x）=1;
        end
        if U（x）>a
            U（x）=a;
        end
        if V（y）<1
            V（y）=1;
        end
        if V（y）>b
            V（y）=b;
        end
        f1=I1（floor（U（x））floor（V（y）））;%f（00）
        f2=I1（floor（U（x））ceil（V（y）））;%f（10）
        f3=I1（ceil（U（x））floor（V（y）））;%f（01）
        f4=I1（ceil（U（x））ceil（V（y）））;%f（11）     
        I2（xy）=（f2-f1）*u+（f3-f1）*v+（f4+f1-f2-f3）*u*v+f1;      
    end
end
I2=uint8（I2）;
subplot（121）;
imshow（I2）;
title（‘向后映射的矫正图像‘）;


%向前映射法
I3=zer