在matlab中进行地理坐标和像素坐标的相互转换

最新推荐文章于 2024-04-25 00:54:06 发布

lingllllove

最新推荐文章于 2024-04-25 00:54:06 发布

阅读量740

点赞数 9

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/lingllllove/article/details/135144682

版权

本文介绍如何在MATLAB中进行地理坐标到像素坐标以及像素坐标到地理坐标的转换，包括使用pix2map函数将像素坐标转化为地理坐标，以及使用map2pix函数将地理坐标转化为像素坐标。此外，还探讨了大地坐标与三维直角坐标系之间的转换方法。

摘要由CSDN通过智能技术生成

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(