思路

• 实现条件
  1.GDAL库
  2.DEM数据

• 开发语言
  python3

• 开发环境
  Linux

• 难度等级
  中

环境

• Ubuntu18.04 LTS
• python 3.6.9

依赖

• GDAL
• numpy
• python
  szZack的博客

安装GDAL

• 方法1：
  inux下python3的GDAL链接如下：
  https://sourceforge.net/projects/gdal-wheels-for-linux/files/
  下载对应python版本的GDAL
  再pip安装：

  pip install GDAL-3.4.1-cp36-cp36m-manylinux_2_5_x86_64.manylinux1_x86_64.whl
  1

• 方法2：

  # 安装的版本需要和python版本一致
  python3 -m pip install https://sourceforge.net/projects/gdal-wheels-for-linux/files/GDAL-3.4.1-cp36-cp36m-manylinux_2_5_x86_64.manylinux1_x86_64.whl/download
  1
  2

  szZack的博客

• 方法3：
  【推荐方法】机器需要可以连外网

  apt-get update
  apt-get install -y python3-gdal
  1
  2

安装依赖包

• 方法1：

  pip install GDAL
  1

• 方法2：
  先下载 GDAL-3.4.1.tar.gz
  然后解压后执行下面的指令：

  python setup.py build
  python setup.py install
  1
  2

测试GDAL

~# python3
Python 3.6.9 (default, Jan 26 2021, 15:33:00)
[GCC 8.4.0] on linux
Type "help", "copyright", "credits" or "license" for more information.
>>> import gdal
1
2
3
4
5

不报错即说明安装成功！

数据依赖

下载高程 tif 格式数据，自行百度一下

python代码

from osgeo import gdal

def get_elevation(tif_path, latitude, longitude):
    
    gdal.UseExceptions()
    ds = gdal.Open(tif_path)
    
    band = ds.GetRasterBand(1)
    elevation = band.ReadAsArray()
    nrows, ncols = elevation.shape
    
    x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()
    
    new_ncols, new_nrows = int((y0-latitude)/dx), int((longitude-x0)/dx)
    
    return elevation[new_ncols][new_nrows]


#经纬度点坐标
latitude, longitude = 31.15, 111.24
h = get_elevation('./xxx.tif', latitude, longitude)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

szZack的博客

参考

• https://pypi.org/project/GDAL/
• https://www.yisu.com/zixun/689106.html
  szZack的博客