1、modis的mod02产品下载后,hdf文件中直接把经纬度存为两个图层,每个像元都有对应值,但是经过HEG软件裁件拼接重采样重新投影以后,这俩图层没有了,只是在结果文件中存储了投影和角点位置信息。
2、大多国产卫星GeoTiff文件,都不会直接提供每个像元的经纬度。
在进行各种定量化计算的时候,逐个像元的经纬度往往直接或者间接需要的,如何得到?
解决:
使用gdal编译完后的exe程序,在cmd中运行处理。gdal编译完的诸多exe中,有一个gdal_translate.exe,其作用是实现文件格式的翻译转换,比如hdf变为tif。
如何通过这条命令来得到每个像元的经纬度值呢?使用下面的命令:
F:gdal_translate -of XYZ mod02resam.hdf mod02resamLatLon.txt
其中,-of意为指定输出格式,我们指定了ASCII格式(简写为XYZ,详见gdal官网format)。
输出的mod02resamLatLon.txt中每个像元为一行,分别为经纬度、像元值,由此,每个像元的经纬度可以从该txt文件读取。