Matlab影像像素坐标得到经纬度/经纬度转影像坐标（已知经纬度获取影像DN值）

我想要利用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了

那么已知经纬度获取DN值就是上述的逆过程，如下

[A,refmat,bbox] = geotiffread('H:\georef.tif');
info = geotiffinfo('H:\georef.tif');
[x1,y1] =projfwd(info,lat,lon);
[row, col ] =map2pix(info.RefMatrix,x1,y1);
z=A(int32(row),int32(col));

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

行和列千万别搞混，可以自己用数据验证一下!