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/); %转换为地理坐标