介绍了如何使用R树提高空间交叉点查询的效率。 提供了使用R树分别计算点到点和线的最近距离的方法。然而,我还没有找到一个解决方案来提高使用R-Tree从点到多边形的距离查询。我遵循了上面提到的例子,但仍然没有想出如何在从点到面的最近距离查询中实现R-Tree。有人能帮我解决这个问题吗?非常感谢。在下面的示例中,我使用了南美洲和GeoPandas的一个城市子集作为示例:
# Import libraries
import geopandas as gpd
import matplotlib.pyplot as plt
# Load data from GeoPandas
world = gpd.
目前,我使用Pandas和ArcPy结合使用Python进行地理信息系统工作。我最近听说了Geopandas,我很有兴趣学习使用它作为基本地质处理操作(多边形的空间连接点、相交多边形等)的替代ArcPy。
我已经安装了Geopandas及其依赖项(gdal、shapely、fiona和pyproj)。
我找到了一个关于空间连接()的Geopandas教程,我一直在“两个GeoDataFrames之间的空间连接”下学习这个教程。
前4行导入必要的模块:
import os
from shapely.geometry import Point
from geopandas import GeoD
我有一个数据与多边形,并希望循环,以创建一个新的列,检查一个特定的点是否在一个或多个多边形(因为他们不是排他性的)。我试着遵循一个地质公园教程,但这似乎不起作用,因为所有的列返回一个“假”。有人能指出我做错了什么吗?非常感谢。以防万一,文件"quartier_paris.geojson“是一个经典的geojson文件,它被读取没有问题(形状出现,我可以去所有的quartiers.explore(),并看到地图上的所有多边形,我百分之百肯定点在这些多边形内)。但是文件的大小太大,不能加载,这是它的链接,以防)
import pandas as pd
import geopandas
im
我需要使用python合并shapefile的多边形。
多边形内部有一些自交,我的python代码总是导致错误。
import geopandas as gpd
from shapely.geometry import Polygon
from shapely.validation import make_valid
from shapely.ops import cascaded_union
from shapely.validation import explain_validity
pz32 = gpd.read_file("B://_Shp_robocze//test//t
我有一个由id和由2D点填充的几何列组成的geopandas dataframe。我想连接每个唯一id的点来创建一个多边形,这样我的新数据框就会有多边形作为它的几何体。我的代码目前看起来像这样:
polygons = geopandas.GeoDataFrame()
for i in id:
group = df[df['id']== i]
polygon = {'type': 'Polygon', 'coordinates': group['geometry']}
polygon[
我有一个形状文件,我想把这个形状分成"n“垂直条。是否有一个纯python方法来完成这个操作,并将输出保存为每个垂直条的形状文件?我有美国的形状文件,想把多边形分成几个垂直条。
import geopandas as gpd
gdf = gpd.read_file('cb_2017_us_nation_20m/cb_2017_us_nation_20m.shp')
AFFGEOID GEOID NAME geometry
0100000US US United States (POLYGON ((-136.563223 58.035052,
我有一个Feature Collection of Polygons和MultiPolygons,我必须先把它写进一个临时文件中,然后用Feature Collection加载它,我正在寻找一种没有临时文件的方法。我尝试过使用geopandas.GeoDataFrame.from_feature(),它对于简单多边形的特性集合很好,但是我不能让它在Feature Collection of Polygons和MultiPolygons上工作,我想做下面这样的事情,但它还没有起作用。
features_collection = []
for feature in json_data['
我正在创建一个地图可视化。我正在映射的内容之一是大量地区的多边形;邮政编码、计数、州……这些实体中的一些将具有数百个经度/经度组合来组成多边形。有没有什么替代方法可以用“new google.maps.latlng”来呈现它们,有没有速记方法?把它写出来几万次似乎是多余的。
我想我可以将其存储为字符串,然后使用val + lat + ','+long + ']‘。
这是一个问题,因为我将这些预先生成的多边形存储在数据库中,它占用了大量空间
像这样的东西,但它不起作用。
var stuff = new google.maps.LatLng;
var Poly14591;
我正在尝试将GeoPandas用于一个(仅稍微)复杂的项目,但目前我无法以投影的方式写出一个简单的shapefile,其中只有一个点。 下面的代码生成了一个通常看起来不错的shapefile --但是.prj是空的: import pandas as pd
from geopandas import GeoDataFrame
from shapely.geometry import Point
df=pd.read_csv("richmond.csv")
geometry = [Point(xy) for xy in zip(df.x, df.y)]
crs = {'