gamma转stamps中经纬度文件的制作

（1）使用mkgrid命令，建立一个和主影像大小一致的pointlist文件.

生成1400*5088的plist文件。

（2）使用data2pt命令，把栅格的dem转换为点格式的，这一步是给第三步用的。

生成pdata文件。

（3）用pt2geo对第一步生成的pointlist，就是主影像的每个像元，生成其对应的空间经纬度数据

其中，生成的plat_lon文件，就是我们需要的文件。这个文件是二进制存储的，里面的数据是这样的：[104.01,30.02,0.1.02,30.03……]，经纬度交替分布，而我们想要的lon文件是[104.01,104.02……]，lat文件是[30.01,30.02……]。

首先需要对这个文件swap以下，使用swap_bytes命令。

然后用Python进行拆分。

注意，拆分后的lon和lat文件，我这了还要做一次swap_bytes 

最后，20180115.lat和20180115.lon文件就是生成的经纬度文件！

转载本文请联系原作者获取授权，同时请注明本文来自吴远昆科学网博客。
链接地址：https://wap.sciencenet.cn/blog-3473451-1328271.html