python快速实现shp点提取栅格值 (两种方法对比)

最新推荐文章于 2024-05-21 09:38:19 发布
最新推荐文章于 2024-05-21 09:38:19 发布
阅读量5.5k 收藏 50
点赞数 5
文章标签： python gis
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41106093/article/details/122525017
版权

大量快速值提取至点的方法

在遥感数据的处理中，常常需要实现的一个功能就是“提取指定shp点对应栅格图像的value值”。实现该功能可以使用ArcGIS中的提取工具-值提取至点。但遇到大量数据，众多实验区的数据时，手动操作着实会带来很大的麻烦。这里我个人改进总结了两种方法，可实现python快速提取shp点对应的栅格值方法。

rasterio+GeoPandas

先上图说明：
底图为CHM，其上为对应的采样点.shp

import matplotlib.pyplot as plt
import geopandas as gpd
import rasterio
from rasterio.plot import show
pointData=gpd.read_file(r'\ABBY.shp')		
raster=rasterio.open(r'\ABBY_CHM.TIF')		#CRS为UTM, EPSG:32610
pointData=pointData.to_crs(
最低0.47元/天 解锁文章
白泽渌
关注
  • 5
    点赞
  • 踩
  • 50
    收藏
    觉得还不错? 一键收藏
  • 10
    评论
GIS数据处理：用Python提取站点处的栅格值
qq_41387508的博客
08-06 3090
在进行地理数据处理的时候，我们经常需要用点矢量来提取栅格数据中的值，在Arcgis或者R里面，这个很好实现，用extract工具或者函数就行，但是我看了下貌似Python里面没有现成的。那就自己写一个吧。 本文示例用一个点矢量的shp文件来批量提取一个栅格数据上各点处的值。 公众号《GIS与Climate》，欢迎关注和交流。 读取数据 分别用rasterio和geopandas包读取栅格和矢量数据： rs = rasterio.open('../input/rs_min_1km.tif'
根据点提取栅格值
weixin_33768481的博客
05-01 4537
栅格数据中每一个栅格都有一个值，比如DEM，每一格都有一个高程值，值都是内插得出的。现在有一点集，需要根据栅格将每一个赋上相应的栅格值。 用到的工具是提取分析中的值提取至点（ExtractValuesToPoints） 当然，如果栅格是多个栅格组成，必须还先把它们镶嵌，使用到的工具是镶嵌到新栅格，需要注意的是参数中"像素类型"的填写方法，如果需要结果的栅格值如原栅格一样，保留多位小数，必...
探索地理空间数据的利器：Rasterstats
最新发布
gitblog_00039的博客
05-21 411
探索地理空间数据的利器：Rasterstats 项目介绍 项目地址:https://gitcode.com/perrygeo/python-rasterstats 在地理信息科学领域，对矢量和栅格数据的融合分析是一项至关重要的任务。rasterstats 是一个强大的 Python 模块，专门设计用于基于矢量几何体对地理空间栅格数据进行总结统计。这个库提供了 分区统计 和 插值点查询 的功能，并且...
通过点获取栅格值
qq_29768651的博客
03-29 249
/// <summary> /// 通过点获取栅格值 /// </summary> /// <param name="raster">栅格对象</param> /// <param name="point">点对象</param> /// <returns></returns> public static doubl.
提取点坐标数据_利用ArcGIS根据点坐标位置提取栅格值
热门推荐
weixin_29019241的博客
01-14 1万+
在工作中经常需要根据给定点的坐标位置，提取栅格影像中的值。在ArcGIS中通过以下几个步骤实现。第一步：在ArcGIS中加载栅格数据用例数据是2015年1月1日，全球潜在蒸散发数据第二步 在Arctoolbox---Spatial Analyst Tools---Extraction中启动ExtractValues to Points第一行参数是坐标文件(在这里我们给大家一些气象站点的...
ArcGIS按点提取栅格
wlh2067的博客
11-19 4580
这里我想提取这些站点的土壤湿度，栅格图是土壤湿度值Spatial Analyst Tools——Extraction——Extract values to Points(注意不要点到下面的。很相似的两个工具)输入点数据，栅格，和输出的矢量文件即可。这样就能完成批量提取了，如图所示，Pro支持一个点图层同时提取多个栅格的内容。如图所示，Pro支持一个点图层同时提取多个栅格的内容。
利用python进行矢量转栅格操作
ljt_217的博客
07-09 1978
利用python进行矢量转栅格操作 基于rasterio和geopandas两个库包对矢量进行栅格化操作 #!usr/bin/env python # -*- coding: utf-8 -*- """ date：2021/7/8 author：甲戌_Tr email: liu_xxxi@163.com """ import sys,os import geopandas as gpd import rasterio from rasterio import features from rasterio.
ArcGIS根据SHP裁剪栅格的Python代码
12-10
编写了一个Python文件，用户可自行修改输入栅格、SHP、字段名、保存路径。代码会根据选择的字段名自动对栅格数据进行裁剪，生成多个栅格文件。
利用Python实现Shp格式向GeoJSON的转换方法
09-19
JSON（JavaScript Object Nonation）是利用键值对+嵌套来表示数据的一种格式，以其轻量、易解析的优点，这篇文章主要介绍了利用Python实现Shp格式向GeoJSON的转换,需要的朋友可以参考下
用python提取shp点文件的坐标信息
01-04
个人觉得比arcGIS方便
python脚本实现shp文件转geojson文件
01-29
针对QGIS中shp转geojson数据容易造成中文乱码问题 提出针对python脚本的一键转换，在代码中只需要更改自己的文件地址，运行即可实现代码的转换。 脚本依靠GIS常用GDAL库实现，安装简单，熟悉基本的python语法即可...
Arcgis使用Python代码将栅格数据批量转成矢量面（shp）
03-15
Arcgis使用Python代码将栅格数据批量转成矢量面（shp）,代码注释详细，只需更改栅格数据及输出文件夹所在路径即可
python 矢量数据转栅格数据代码实例
09-18
主要介绍了python 矢量数据转栅格数据代码实例,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
点转栅格批量处理arcpy
08-10
本程序基于arcpy编写，用于点数据批量转为栅格，可根据需要修改
关于Python 的简单栅格图像边界提取方法
09-19
今天小编就为大家分享一篇关于Python 的简单栅格图像边界提取方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
Python-提取矢量边界
04-23
Python-提取矢量边界
利用shp矢量数据裁剪遥感影像的Python实现
02-25
本文将深入探讨如何使用Python语言，结合Shp矢量数据对遥感影像进行裁剪操作，以便提取特定区域的信息。Shp文件是ESRI（Environmental Systems Research Institute）公司开发的地理信息系统中广泛使用的矢量数据格式...
python 矢量转栅格
weixin_41035425的博客
04-14 1718
虽然利用arcpy.PolygonToRaster_conversion可以将矢量文件转为栅格，但是arcpy的速度令人堪忧，有没有更快更好的办法呢？对于经常用python处理数据的人来说，会优先考虑用python来解决。百度一下“python 将矢量转栅格”，会有别人分享的一些代码，但是有些代码我自己使用后，发现一些bug。于是在他们的基础上，我对代码进行了优化和修改，封装成一个函数，经过自己数据的测试没问题，分享给大家。 import gdal import ogr import gdalc...
解决Python利用shp进行栅格tif数据计算时出现的问题：TypeError: can‘t multiply sequence by non-int of type ‘Affine‘
andyfradar的博客
04-29 393
利用到了rasterio这个库，进行shp的几何平面提取它的dem的均值，问题在于rasterio预计接收得是两个值，但我们shp给它传递了3个值，也就是说shp层不只是包含经纬度，它还包含了一个高度值，因此，我们的解决办法就是去掉shp的高度值，使它只包含几何的经纬度。首先查看shp的属性表，看它是否是包含高度值的，包含高度值时，它的【shape *】属性这一列几何要素后面会加一个【ZM】，而没有包含高度值的几何要素后面不会有【ZM】打开Toolbox的【数据管理工具】-【要素】-【复制要素】，
如何在python中利用矢量文件提取栅格值
05-31
在中，可以使用GDAL库来读取和处理矢量和栅格数据。以下是一些简单的步骤，可以利用GDAL库从矢量文件中提取栅格值： 1. 导入gdal包 ```python import gdal ``` 2. 打开栅格文件和矢量文件 ```python raster_file = gdal.Open('path/to/raster.tif') vector_file = ogr.Open('path/to/vector.shp') ``` 3. 获取栅格文件的坐标系 ```python proj = raster_file.GetProjection() ``` 4. 获取栅格文件的地理变换信息 ```python transform = raster_file.GetGeoTransform() ``` 5. 获取矢量文件的图层 ```python layer = vector_file.GetLayer() ``` 6. 遍历矢量文件的要素，并提取栅格值 ```python for feature in layer: geom = feature.GetGeometryRef() mx, my = geom.GetX(), geom.GetY() # 获取要素的中心点坐标 # 将地理坐标转换为栅格坐标 px = int((mx - transform[0]) / transform[1]) py = int((my - transform[3]) / transform[5]) # 获取栅格值 band = raster_file.GetRasterBand(1) value = band.ReadAsArray(px, py, 1, 1)[0][0] print(value) ``` 这样就可以从矢量文件中提取栅格值了。需要注意的是，在使用GDAL库时，需要先安装GDAL库和相关的依赖库。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 10
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值