% http://maps.live.com/default.aspx?v=2&cp=45.40833~11.89368&style=h&lvl=19
% &tilt=-90&dir=0&alt=143.849469123408&scene=10539604&sp=Point.rk17gjhzr0sc
% _Untitled%20item___&encType=1

%%
close all, clear all
im = imread('mappa_fotografica.jpg');
imshow(im,[])

%% TIRA UNA LINEA SUL LATO LUNGO, DAL BASSO ALL'ALTO DELLA FIGURA
figure
im2 = imread('mappa_fotografica2.jpg');
imshow(im2,[]);
[x,y] = getline(gcf)

%% 
v = [x(2)-x(1) y(2)-y(1)];
phi = pi/2 + atan2(v(2),v(1));
phi_degree = phi*180/pi

%%
R = [cos(phi) sin(phi) 0; -sin(phi) cos(phi) 0; 0 0 1]

%%
figure
imshow(imtransform(im,maketform('projective',R')),[])

%%
save R R



