clear all;close all;clc;

img1=imread('TSD.png');

img1=rgb2gray(img1);

%[h1 w1]=size(img1);

imshow(img1);

in_points =1.0e+03 *[

? ? 0.55? ? 0.43

? ? 0.74? ? 0.43

? ? 0.25? ? 0.88

? ? 1.26? 0.88 ];

%ginput();

w1=100;h1=100;

out_points = [1,1;w1,1;1,h1;w1,h1];

%p1=[1,1;w1,1;1,h1;w1,h1];

%p2=ginput();? ? ? ? %依次點(diǎn)擊待變換區(qū)域，左上、右上、左下、右下

%%方法一

%T = maketform('projective', in_points, out_points);

%[imgn X Y]=imtransform(img1,T);? ? %投影

%%方法二

T=fitgeotrans(in_points,out_points,'projective');

imgn=imwarp(img1,T);?

figure;imshow(imgn)