python如何利用fiona判定那些point在那些polygon内

2023-11-05 05:10

本文主要是介绍python如何利用fiona判定那些point在那些polygon内,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

如何利用fiona判定那些point在那些polygon内)

  • 什么是foina
  • 几何类型
  • 什么是shapely
  • 代码展示
  • 利用了foina的那些函数

问题来源,我有两个数据数据:一个是矢量面数据、一个是矢量point数据,那么如何判定那些point数据在包含在矢量面范围内,并把point在对应面的 属性表的字段,写入到point内;

什么是foina

fiona是一个用于读写矢量空间数据文件的Python包,是由Sean Gillies等人写的,其底层是C++开发的OGR库(一个开源GIS库)。

安装fiona包最好是下载wheel文件(编译发布的包)进行安装,因为用源码发布的包进行安装有一个编译(build)过程,如安装机器上没有相应的编译器,则会出现安装错误。此外,fiona的依赖包GDAL也是同样情况。
GDAL和fiona的wheel文件可以从Unofficial Windows Binaries for Python Extension Packages网站上下载。安装fiona前先要安装GDAL。

在这里插入图片描述
在这里插入图片描述
要根据安装机器的操作系统是64位还是32位;以及Python版本是python 2或者 python 3选择安装对应文件;

读文件首先是要利用open()函数的读模式返回Collection对象。
Collection对象有很多属性,利用Collection对象属性可以了解文件的基本情况。

属性解释
driver文件类型,如‘ESRI Shapefile‘。
crsproj4字典形式表示的坐标参照系统。
crs_wktwkt形式表示的坐标参照系统。
bounds边界范围,最小、最大x和y组成的元组(minx, miny, maxx, maxy)。
schema数据结构,包括几何类型及字段构成。
meta包括driver、schema、crs和其它属性。
import pprint
import fiona
c = fiona.open(r"d:\data\c.shp")
pprint.pprint(c.schema)

Collection对象是可迭代对象,可以通过循环返回每条记录,也可以通过索引返回其中的一条记录。
返回的记录是一个字典,包含id、geometry、properties、type等键,类似GeoJSON中的Feature。

在一个记录中,geometry键的值包含了该记录的几何类型及坐标值。
fiona定义的几何类型和GeoJSON是一致的,包括Point、LineString、Polygon、MultiPoint、MultiLineString以及MultiPolygon。
坐标值的表示方式与几何类型是相关的.

几何类型

shapefile文件的几何(二维)类型只有4种:点(Point)、多点(MultiPoint)、线(Polyline)和多边形(Polygon),在文件层面没有区分LineString和MultiLineString、Polygon和MultiPolygon,但坐标的表示形式是不一样,可以通过每个记录的坐标进行区分。

在一个记录中,geometry键的值包含了该记录的几何类型及坐标值。
fiona定义的几何类型和GeoJSON是一致的,包括Point、LineString、Polygon、MultiPoint、MultiLineString以及MultiPolygon。
坐标值的表示方式与几何类型是相关的。

几何类型坐标值形式
Point点的x和y坐标组成的元组
LineString多个点坐标组成的列表(一层列表)
Polygon多个环线(ring)坐标组成的列表(每个环线是多个点坐标组成的列表)(两层列表)
MultiPoint多个点坐标组成的列表(一层列表)
MultiLineString多个线坐标组成的列表(两层列表)
MultiPolygon多个多边形坐标组成的列表(三层列表)

注:Polygon的坐标形式实际上是由一个外边线和多个内边线坐标组成的列表,这种形式可以表示所有多边形(包括有孔洞的多边形)。

shapefile文件的几何(二维)类型只有4种:点(Point)、多点(MultiPoint)、线(Polyline)和多边形(Polygon),在文件层面没有区分LineString和MultiLineString、Polygon和MultiPolygon,但坐标的表示形式是不一样,可以通过每个记录的坐标进行区分。

什么是shapely

shapely是一个用于几何对象构建与操作的Python包,底层是用C++写的GEOS库,很多开源的GIS软件(如QGIS)是利用GEOS库开发的。

Shapely访问地址

shapely包定义了多个类型的几何对象,其中,Point、MultiPoint、LineString、MultiLineString、Polygon以及MultiPolygon和fiona定义的几何类型是一致的,此外,还包括LinearRing、box等。
shapely包中的geometry模块提供了创建几何对象的构造函数
几何对象的构建是在构造函数中传入坐标参数。

from shapely.geometry import MultiPolygon
polygon1 = Polygon([(0, 0), (1, 1), (1, 0)])
polygon2 = Polygon([(1, 1.1), (1, 2), (2, 2)])
polygons = MultiPolygon([polygon1,polygon2])
polygons[0]

几何对象有一些共同属性,如area,尽管点和线的area属性值恒定为0;也有一些属性是和几何对象的类型相关,如Polygon有exterior和interiors属性,但其它对象没有。

coords属性只有Point、LineString和LinearRing这些几何对象才具有,返回的值是一个坐标序列对象,可以利用list()函数或切片操作([:])把坐标序列对象转换成坐标序列列表,或利用索引返回某个点的坐标。

几何对象的方法主要有几个方面:
计算对象之间的距离,返回距离值。
对几何对象进行空间操作,产生新的几何对象。
对几何对象进行空间拓扑分析,返回True或False。

方法解释
distance(other):计算与另一个几何对象的最近距离。
hausdorff_distance(other)计算与另一个几何对象的最远距离。
from shapely.geometry import Point
from shapely.geometry import LineString
point = Point([1,1])
line = LineString(((0,0), (1,1), (1,0)))
dist = point.distance(line)
h_dist = point.hausdorff_distance(line)
print(dist)
print(h_dist)
方法解释
buffer()产生缓冲多边形。
intersection(other)与另一个几何对象求交。
union(other)与另一个几何对象求和。
difference(other)去除与另一个几何对象的重叠部分。
symmetric_difference(other)与另一个几何对象求和,并去除重叠部分。
from shapely.geometry import Point
a = Point(1, 1).buffer(1.5)
b = Point(2, 1).buffer(1.5)
result = a.intersection(b)
result

在这里插入图片描述
先进行缓冲操作,然后进行求交操作.

方法解释
equals(other)和另一个几何对象是否一致。
contains(other)是否包含另一个几何对象。
within(other)是否包含在另一个几何对象中。
crosses(other)是否和另一个几何对象相交(只有线和线、线和多边形有crosses关系),不包括包含关系。
intersects(other)是否和另一个几何对象相交或有包含关系。
disjoint(other)是否和另一个几何对象没有公共点。
from shapely.geometry import LineString
from shapely.geometry import Polygon
line = LineString([(0.1,0.1), (0.9,0.9)])
polygon = Polygon([(0,0), (1,0), (1,1), (0,1)])
print(line.crosses(polygon))
print(line.intersects(polygon))

在这里插入图片描述
判断线和多边形是否crosses和intersects

代码展示

本文主要实现的是:将按矢量面范围剪裁处理的另一份矢量point数据,逐个点进行与 感兴趣区面进行intersects()叠加分析,如果是包含,那么就把矢量面包含的point数据写入到新的shp数据中,然后将矢量面数据的其中一个fileds写入到新的point数据的属性表内部;

from shapely.geometry import LineString
from shapely.geometry import Point
from shapely.geometry import Polygon
from shapely.geometry import MultiLineString
import fionakansas_c = fiona.open('D:/20200309/shujushouji/python+jupyternoterbook/zuoye1(1)/6栅格数据操作/sample_regions/sample_regions.shp', 'r')
#coordinates = kansas_c[0]['geometry']['coordinates'][0]
#polygon = Polygon(coordinates)
#kansas_c.close()high_c = fiona.open('D:/20200309/shujushouji/python+jupyternoterbook/zuoye1(1)/6栅格数据操作/sample_regions/grid_points1.shp', 'r')
source_driver = high_c.driver
source_crs = high_c.crs
source_schema = high_c.schema
target = fiona.open('D:/20200309/shujushouji/python+jupyternoterbook/zuoye1(1)/6栅格数据操作/sample_regions/shuchupoint.shp', 'w',driver=source_driver,crs=source_crs, schema=source_schema)
for record in high_c:coordinates = record['geometry']['coordinates']#print(coordinates)record_line = Point(coordinates)for i in range(0, len(kansas_c)):coordinates = kansas_c[i]['geometry']['coordinates'][0]polygon = Polygon(coordinates)if record_line.intersects(polygon):FIELDSS = kansas_c[i]['properties']['type']#print(FIELDSS)record['properties']['land_type']=FIELDSStarget.write(record)
target.close()
high_c.close()

利用了foina的那些函数

主要是:
(1)利用Collection对象的write(feature)或writerecords(feature_list)方法可以把要素或要素列表写入Collection对象中。要素写入到Collection对象后,并没并没有立即保存到文件,只有当关闭Collection对象时,Collection对象中的数据才会写到文件中。

Collection对象是可迭代对象,可以通过循环返回每条记录,也可以通过索引返回其中的一条记录。
返回的记录是一个字典,包含id、geometry、properties、type等键,类似GeoJSON中的Feature。

(2)利用sharpely空间拓扑分析
,利用空间分析的intersects()叠加分析确认是否:是否和另一个几何对象相交或有包含关系。

这篇关于python如何利用fiona判定那些point在那些polygon内的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/347628

相关文章

Python基础文件操作方法超详细讲解(详解版)

《Python基础文件操作方法超详细讲解(详解版)》文件就是操作系统为用户或应用程序提供的一个读写硬盘的虚拟单位,文件的核心操作就是读和写,:本文主要介绍Python基础文件操作方法超详细讲解的相... 目录一、文件操作1. 文件打开与关闭1.1 打开文件1.2 关闭文件2. 访问模式及说明二、文件读写1.

Python将博客内容html导出为Markdown格式

《Python将博客内容html导出为Markdown格式》Python将博客内容html导出为Markdown格式,通过博客url地址抓取文章,分析并提取出文章标题和内容,将内容构建成html,再转... 目录一、为什么要搞?二、准备如何搞?三、说搞咱就搞!抓取文章提取内容构建html转存markdown

Python获取中国节假日数据记录入JSON文件

《Python获取中国节假日数据记录入JSON文件》项目系统内置的日历应用为了提升用户体验,特别设置了在调休日期显示“休”的UI图标功能,那么问题是这些调休数据从哪里来呢?我尝试一种更为智能的方法:P... 目录节假日数据获取存入jsON文件节假日数据读取封装完整代码项目系统内置的日历应用为了提升用户体验,

Python FastAPI+Celery+RabbitMQ实现分布式图片水印处理系统

《PythonFastAPI+Celery+RabbitMQ实现分布式图片水印处理系统》这篇文章主要为大家详细介绍了PythonFastAPI如何结合Celery以及RabbitMQ实现简单的分布式... 实现思路FastAPI 服务器Celery 任务队列RabbitMQ 作为消息代理定时任务处理完整

Python Websockets库的使用指南

《PythonWebsockets库的使用指南》pythonwebsockets库是一个用于创建WebSocket服务器和客户端的Python库,它提供了一种简单的方式来实现实时通信,支持异步和同步... 目录一、WebSocket 简介二、python 的 websockets 库安装三、完整代码示例1.

揭秘Python Socket网络编程的7种硬核用法

《揭秘PythonSocket网络编程的7种硬核用法》Socket不仅能做聊天室,还能干一大堆硬核操作,这篇文章就带大家看看Python网络编程的7种超实用玩法,感兴趣的小伙伴可以跟随小编一起... 目录1.端口扫描器:探测开放端口2.简易 HTTP 服务器:10 秒搭个网页3.局域网游戏:多人联机对战4.

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

Python使用自带的base64库进行base64编码和解码

《Python使用自带的base64库进行base64编码和解码》在Python中,处理数据的编码和解码是数据传输和存储中非常普遍的需求,其中,Base64是一种常用的编码方案,本文我将详细介绍如何使... 目录引言使用python的base64库进行编码和解码编码函数解码函数Base64编码的应用场景注意

Python基于wxPython和FFmpeg开发一个视频标签工具

《Python基于wxPython和FFmpeg开发一个视频标签工具》在当今数字媒体时代,视频内容的管理和标记变得越来越重要,无论是研究人员需要对实验视频进行时间点标记,还是个人用户希望对家庭视频进行... 目录引言1. 应用概述2. 技术栈分析2.1 核心库和模块2.2 wxpython作为GUI选择的优

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.