A和B我定义空间参考坐标系
RA = imref3d(大小(A), [min (xA), max (xA)], [min(是的),max(丫)],[min (zA), max (zA)]);
RB = imref3d(大小(B)、(min (xB), max (xB)], [min (yB), max (yB)], [min (zB), max (zB)]);
然后我定义了一个仿射变换和变换矩阵TM
在= affine3d (TM);
现在我可以变换B的坐标系统
(C, RC) = imwarp (B, RB,,“OutputView”RA);
[rp, Bfit, ErrorStats] = absor (PA、PB、“DoScale”1);
TM = rp.M。”;