坐标转换（四参数）

[row col]=size(source)

for i=1:row

u(i)=source(i,2)

v(i)=source(i,3)

x(i)=source(i,4)

y(i)=source(i,5)



B(2*i-1,1)=1

B(2*i-1,2)=0

B(2*i-1,3)=u(i)

B(2*i-1,4)=-v(i)



B(2*i,1)=0

B(2*i,2)=1

B(2*i,3)=v(i)

B(2*i,4)=u(i)



L(2*i-1,1)=x(i)

L(2*i,1)=y(i)

end

X=inv(B *B)*(B *L)



Txy=[X(1) X(2)]

xzjz(1,1)=X(3)

xzjz(1,2)=-X(4)

xzjz(2,1)=X(4)

xzjz(2,2)=X(3)



[row1 col1]=size(target)

for i=1:row1

u1(i)=target(i,2)

v1(i)=target(i,3)

d=[u1(i) v1(i)]

md=Txy+xzjz*d

end

end

-Coordinate conversion：

[row col]=size(source)

for i=1:row

　　u(i)=source(i,2)

　　v(i)=source(i,3)

　　x(i)=source(i,4)

　　y(i)=source(i,5)

　　

　　B(2*i-1,1)=1

　　B(2*i-1,2)=0

　　B(2*i-1,3)=u(i)

　　B(2*i-1,4)=-v(i)

　　

　　B(2*i,1)=0

　　B(2*i,2)=1

　　B(2*i,3)=v(i)

　　B(2*i,4)=u(i)

　　

　　L(2*i-1,1)=x(i)

　　L(2*i,1)=y(i)

end

X=inv(B *B)*(B *L)



Txy=[X(1) X(2)]

xzjz(1,1)=X(3)

xzjz(1,2)=-X(4)

xzjz(2,1)=X(4)

xzjz(2,2)=X(3)



[row1 col1]=size(target)

for i=1:row1

　　u1(i)=target(i,2)

　　v1(i)=target(i,3)

　　d=[u1(i) v1(i)]

　　md=Txy+xzjz*d

end

end