（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文件就是生成的经纬度文件！

参考：www.jianshu.com/p/a34bc5c54cc5