1. 问题描述
进行地理信息处理时,通常需要获得已知经纬度处在SRTM3文件中的高程值,这需要首先根据该经纬度来确定相应SRTM3的文件名,其次,确定该经纬度在相应文件中的行与列,即相应的栅格点。这里需要注意的是,由于给定的经纬度的值通常不会恰巧在栅格点上,因此,可以考虑距离最近的栅格点的值作为该经纬度的高程。最后,采用相应的软件将相应点的数据读出即可。
2. 根据经纬度确定SRTM3文件名称
SRTM3的数据组织方式为:每5度经纬度方格划分为一个文件,文件内包含6000x6000个采样点的高程数据。文件共分为24行(-60至60度)和72列(-180至180度)。文件命名规则为SRTM_XX_YY.zip(或SRTM_XX_YY.tiff),XX表示列数(01-72),YY表示行数(01-24)。示意图如下:
由上图可见,对于文件SRTM_XX_YY而言,其左上角的点的经纬度分别为
和
因此,假设已知经纬度分别为J和W,则该点所在文件的左上角经纬度与该点经纬度之间的关系为
和
将和
代入不等式得
和
即
式中,表示向下取整运算。
3. 确定已知点在SRTM3文件中的位置
每个数据文件中的数据组织细节可参考下图所示的两个文件 SRTM_38_03和SRTM_38_02相邻处的数据关系。
约定与
分别表示SRTM3文件中数据的列与行的编号。因为SRTM3文件数据精度为3 arc-seconds,所以距离已知点
最近的栅格点的列与行分别为
和
4. 读出相应栅格点高程
可以利用不同的工具,将相应位置的数据读出。具体读出方法,可以参考其它博文。