已知点的经纬度坐标计算/读取DEM高程信息/海拔高度

szZack

已于 2022-08-10 22:08:31 修改

阅读量3.2k

点赞数 1

分类专栏：常见问题文章标签：经纬度坐标 DEM 高程信息海拔

于 2022-08-10 11:14:22 首次发布

本文链接：https://blog.csdn.net/zengNLP/article/details/126262286

版权

常见问题专栏收录该内容

34 篇文章 1 订阅

订阅专栏

已知点的经纬度坐标计算/读取DEM高程信息/海拔高度

思路
环境
依赖
安装GDAL
安装依赖包
测试GDAL
数据依赖
python代码
参考

思路

实现条件
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
```

方法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

szZack的博客

方法3：
【推荐方法】机器需要可以连外网
```
apt-get update
apt-get install -y python3-gdal
```

安装依赖包

方法1：
```
pip install GDAL
```
方法2：
先下载 GDAL-3.4.1.tar.gz
然后解压后执行下面的指令：
```
python setup.py build
python setup.py install
```

测试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

不报错即说明安装成功！

数据依赖

下载高程 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)