clc;close all;clear;
%地理坐标和像素坐标的相互转换
[pic,R]=geotiffread('boston.tif'); %读取带地理坐标信息的tif影像
[m,n,~]=size(pic); %像素坐标转换为地理坐标
figure(),imshow(pic),title('不带地理坐标的静态图片');
hold on; scatter(n/,m/,,'r.'); %选择1/4处像素坐标，并在图上标示
[lon,lat]=pix2map(R,m/,n/); %转换为地理坐标
figure(),mapshow(pic,R);
mapshow(lon,lat,'Marker','.','MarkerEdgeColor','r');
title('带地理坐标的tif影像');
axis off;
%打印结果
disp(['(',num2str(m/),',',num2str(n/),') -> (',num2str(lon),',',num2str(lat),')']); %地理坐标转换为像素坐标。注意地理坐标从左下角开始，而像素坐标从左上角开始
%左下角开始点地理坐标为[R.XLimWorld(),R.YLimWorld()]
%R.RasterWidthInWorld表示图像的地理宽度
x=R.XLimWorld()+(/)*R.RasterWidthInWorld; %获取图像3/4处的地理坐标
y=R.YLimWorld()+(/)*R.RasterHeightInWorld;
figure(),mapshow(pic,R),axis off;
mapshow(x,y,'Marker','*','MarkerEdgeColor','r');
[row,col]=map2pix(R,x,y);
figure(),imshow(pic);
hold on;
scatter(col,row,,'r*');
%命令行中打印结果
disp(['(',num2str(x),',',num2str(y),') -> (',num2str(row),',',num2str(col),')']);
像素坐标转地理坐标：pix2map
地理坐标转像素坐标：map2pix
