在matlab中進行地理坐標和像素坐標的相互轉換

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