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