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

 像素坐标转地理坐标：pix2map

地理坐标转像素坐标：map2pix