t_concord = cp2tform(input_points,base_points,'linear conformal');
info=imfinfo('cameraman.tif');
registered = imtransform(input,t_concord, 'XData',[1 info.Width], 'YData',[1 info.Height]);
figure,imshow(base);
figure, imshow(registered)
X=0;
Y=0;
A=double(base)-double(registered);
B=A.*A;
for i=1:256*256;
if registered(i)~=0
X=X+B(i);
else
Y=Y+1;
end
end
MSE=X/(256*256-Y)
评论0