orthophoto = imread('cropped_image.tif');
figure, imshow(orthophoto)
unregistered = imread('cropped_image2.tif');
figure, imshow(unregistered)
pause
[movingPoints, fixedPoints] = cpselect(unregistered,orthophoto,'Wait',true);
cpselect(unregistered, orthophoto)
mytform = cp2tform(movingPoints, fixedPoints,'projective');
t = fitgeotrans(movingPoints, fixedPoints,'projective');
registered = imtransform(unregistered, mytform);
Rfixed = imref2d(size(orthophoto));
registered = imwarp(unregistered,t,'OutputView',Rfixed);
figure1 = imshow(registered);
saveas(figure1,'undistorted.tif');
0 Comments
Sign in to comment.