如何获取特定经纬度在SRTM3中的高程值

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. 读出相应栅格点高程

可以利用不同的工具，将相应位置的数据读出。具体读出方法，可以参考其它博文。