Matlab像素得到经纬度/经纬度转像素坐标

我想要利用matlab读取geotif后,输入经纬度得到该图像的坐标以及DN值。

首先在网上能找到的是利用图像的坐标值得到经纬度坐标,就是这个博客。

那么首先是

[A,refmat,bbox] = geotiffread('H:\georef.tif');
info = geotiffinfo('H:\georef.tif');

A获取的是图像的矩阵本身,info好像获取的是各种图像的信息吧。

之后

[x,y] = pix2map(info.RefMatrix,3303, 3278);
[lat,lon] = projinv(info,x,y);

pix2map函数是将3303,3278这个行列号转化为直角坐标系的坐标,官网说的是Convert pixel coordinates to map coordinates,(话说这个to变成2也是很搞笑。。。我还以为是版本的意思)

projinv函数是直角坐标投影到地理坐标(直角坐标系坐标不是tif图像的行列号!!!)

就可以得到lat和lon了

那么逆过程如下,就很清晰了

[x1,y1] =projfwd(info,lat,lon);
[row, col ] =map2pix(info.RefMatrix,x1,y1);
z=A(int32(col),int32(row));

注意的是这个int32是为了将row和col是double类型转换成int类型

已标记关键词 清除标记
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页