读取Geotiff栅格文件赋值给点

最新推荐文章于 2023-10-23 16:39:47 发布

sunshinei928

最新推荐文章于 2023-10-23 16:39:47 发布

阅读量497

点赞数

文章标签： python

本文链接：https://blog.csdn.net/sunshinei928/article/details/125785430

版权

本文讲述了在Python 3.8环境下安装Rasterio库时遇到的问题，包括依赖GDAL库的编译挑战和对numpy版本的要求。作者提供了下载预编译包的建议，并展示了如何通过这些库处理DEM数据并解决实际应用中的高程赋值问题。

摘要由CSDN通过智能技术生成

python3.8安装rasterio库遇到几个问题，

1、rasterio库依赖GDAL库，

2、rasterio库对numpy库版本有要求。

直接安装GDAL库会编译要安装vc 2005那个，4G大小，太大了，需要直接下载wheel，

在https://www.lfd.uci.edu/~gohlke/pythonlibs/下载

rasterio‑1.2.10‑cp38‑cp38‑win_amd64.whl
GDAL‑3.4.3‑cp38‑cp38‑win_amd64.whl

另外numpy版本升级到1.22.3，否则会提示np错误。

import rasterio as rio##要安装GDAL库
import pandas as pd
import os

shotfile=r'E:\prc\python\多因素量化难度\炮点.xlsx'
filesave=os.path.splitext(shotfile)[0]
rasterfile=r'E:\prc\python\多因素量化难度\dem.tif'
shotdata=pd.read_excel(shotfile,usecols=[0,1,2,3])###只读取线号，点号，x，y坐标
shotdata['高程']=5000
rast=rio.open(rasterfile)##读取栅格图像
temp_r=rast.read(1)##读取1波段
for i in range(le