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
转载于:https://www.cnblogs.com/denny402/p/4684770.html
相关资源：经纬度坐标转换为平面坐标的matlab实现
