function [refimgout] = warpimgrigid2d(refimg, rigidmatrix, movimg, interptype)
% 将movimg图像根据刚性配准矩阵warp到refimg, 计算ref img上个像素点的值
if nargin < 3
error('input params are too small');
end
if nargin == 3
interptype = 'cubic';
end
size2dref = [refimg.dim(1), refimg.dim(2)];
% 插值的点 movimg变换到refimg
xi = zeros(size2dref); %height方向
yi = zeros(size2dref); %width方向
size2dmov = [movimg.dim(1), movimg.dim(2)];
% xmov = zeros(size2dmov); % 原始的点 movimg
% ymov = z