博文

gdal_translate 计算每个像元的经纬度解决

已有 7665 次阅读 2016-3-10 11:22 |个人分类:遥感地信|系统分类:科研笔记|关键词:学者

问题：

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文件读取。

转载本文请联系原作者获取授权，同时请注明本文来自陈兴峰科学网博客。
链接地址：https://m.sciencenet.cn/blog-474887-961746.html

上一篇：C# String.Replace 无效解决办法
下一篇：清明遇到的熟脸人

数据加载中...

返回顶部

扫一扫，分享此博文