我想要利用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类型
行和列千万别搞混,可以自己用数据验证一下!