问题标签 [geopandas]

我正在测试 geopandas 以使事情变得非常简单：使用差异方法删除一个圆内的 GeoDataFrame 的一些点。

这是我的脚本的开头：

这是 points_df 的第一行：

然后，我添加一个点，该点将包含第一个 GeoDF 的几个点：

这是 iPython notebook 显示的结果：

现在，目标是删除绿色圆圈内的红点。为此，我认为差异方法就足够了。但是当我写：

结果是...... points_df 没有任何改变：

我猜diff（）方法仅适用于多边形GeoDataFrames，并且点和多边形之间的混合是不可能的。但也许我错过了什么！

在这种情况下，测试圆中是否存在点的函数会比差分方法更好吗？

我正在尝试在 shapefile 的多边形内找到一个点。

我需要编写一个循环来遍历多边形并返回该点所在的多边形的索引。

我将如何编写一个循环来找出该点所在的多边形？

这是我到目前为止所写的：

输出可能是这样的：'3001100'（正确多边形的 BCTCB2010）

我正在尝试 geopandas 来操纵一些点数据。我的最终 GeoDataFrame 在那里表示：

为了使用另一个 Python 模块，它使用 OSM 数据计算两点之间的最短路径，我必须像 tour 一样对我的点进行排序。

如果不是，则下一个 Python 模块计算最短道路，但不一定在最近点之间。主要问题是旅游的约束。

如果我的点只是在一条线上，那么每个点的纬度和经度的基本排序功能就足够了，例如：

如果我们在此排序之后从上一张图片的“上”点开始，DataFrame 的第二个点将是它的最近点，等等......直到第五个点，它在图片的右侧（所以不是最近的不再）...

所以我的问题是：有人知道如何实现这种特殊的排序，还是我必须手动更改索引？

我试过这个，但我的前景很混乱

{ "type": "FeatureCollection", "features":[{"type":"Feature","geometry":{"type":"Point","coordinates": [-73.87057,40.773757]},"properties ":{"tpep_dropoff_datetime":"2013-04- 02T16:00:00.000","trip_distance":"11.27999999999999","dropoff_longitude":"-73.870570000000001","pickup_latitude":"40.73289700000000001","toll ","tip_amount":"0","payment_type":"CRD","fare_amount":"37","pickup_longitude":"-73.991167000000004","passenger_count":"6","store_and_fwd_flag":null,"额外":"0","供应商":"VTS","pickup_location":{"type":"Point","coordinates":[-73.991167,40.732897]},"total_amount":"37.5","tpep_pickup_datetime":"2013-04-02T15:22:00.000 ","dropoff_latitude":"40.773757000000003","ratecodeid":"1","mta_tax":"0.5"}},{"type":"Feature","geometry":{"type":"Point", "坐标":[-74.00296,40.722112]},"属性":{"tpep_dropoff_datetime":"2013-07-19T07:52:00.000","trip_distance":"5.5","dropoff_longitude":"-74.002960000000002", "pickup_latitude":"40.766105000000003","tolls_amount":"0","tip_amount":"3.8999999999999999","payment_type":"CRD","fare_amount":"19.5","pickup_longitude":"-73.954407000000003","passenger_count":"1","store_and_fwd_flag":null,"extra":"0","vendorid" :"VTS","pickup_location":{"type":"Point","coordinates":[-73.954407,40.766105]},"total_amount":"23.899999999999999","tpep_pickup_datetime":"2013-07-19T07:33 :00.000","dropoff_latitude":"40.722112000000003","ratecodeid":"1","mta_tax":"0.5"}},{"type":"Feature","geometry":954407000000003","passenger_count":"1","store_and_fwd_flag":null,"extra":"0","vendorid":"VTS","pickup_location":{"type":"Point","coordinates": [-73.954407,40.766105]},"total_amount":"23.899999999999999","tpep_pickup_datetime":"2013-07-19T07:33:00.000","dropoff_latitude":"40.722112000000003","ratecodeid":"1","mta_tax ":"0.5"}},{"type":"Feature","geometry":954407000000003","passenger_count":"1","store_and_fwd_flag":null,"extra":"0","vendorid":"VTS","pickup_location":{"type":"Point","coordinates": [-73.954407,40.766105]},"total_amount":"23.899999999999999","tpep_pickup_datetime":"2013-07-19T07:33:00.000","dropoff_latitude":"40.722112000000003","ratecodeid":"1","mta_tax ":"0.5"}},{"type":"Feature","geometry":766105]},"total_amount":"23.899999999999999","tpep_pickup_datetime":"2013-07-19T07:33:00.000","dropoff_latitude":"40.722112000000003","ratecodeid":"1","mta_tax":"0.5 "}},{"type":"Feature","geometry":766105]},"total_amount":"23.899999999999999","tpep_pickup_datetime":"2013-07-19T07:33:00.000","dropoff_latitude":"40.722112000000003","ratecodeid":"1","mta_tax":"0.5 "}},{"type":"Feature","geometry":

我有一个读入 geopandas 数据框的 .shp 文件。我将坐标参考系统更改为 2163，因为我正在制作一些矩形地图并希望它们看起来有些正常。

数据框如下所示：

现在我正在计算一个成对的距离矩阵，但由于这些坐标以米为单位，所以数字非常大，我遇到了溢出问题。所以我试试这个：

但是我的多边形的坐标并没有变小（只是移动了一点，这很好）。打印据称重新缩放的数据框给出了这个：

我希望得到这样的东西。

也许我只是不明白shapely 的缩放功能应该做什么，尽管它看起来很简单。

我正在尝试用某种图案而不是颜色填充多边形（因为我必须打印项目 B&W）。到目前为止，我只设法用灰度填充它们，但如果有任何方法，我想采用模式方法。

我的一段代码如下：

这个 if 语句在一个遍历 info 和 shape 的 for 循环中。如果某些信息事实与字符串匹配，我只需为多边形着色并将其作为补丁添加到 shapefile 上。

我还想以某种方式将此模式添加到图例中。

提前致谢。

我试图根据点地面测量创建河流横截面剖面。当尝试LineString从具有公共 ID 的一系列点创建 Shapely 时，我意识到给定点的顺序确实很重要，因为它LineString只会“按索引”连接给定点（以列表给定的顺序连接点）。下面的代码说明了默认行为：

这将产生图像：

问题： Shapely 中是否有任何内置方法可以通过给定的随机 2D 点列表自动创建最合乎逻辑的（又名：最短、最简单、最不交叉……）线？

与默认（红色）相比，您可以在下面找到所需的线（绿色）。

我对python有点陌生，所以我希望我的问题的答案相对简单。

我正在尝试使用 geopandas 制作等值线图。但是，由于我正在制作多个需要相互比较的地图，因此使用自定义数据分类方案（而不是分位数或詹克斯）是必不可少的。因此，我一直在尝试使用 User_Defined 方案，并且我能够创建 bin，但我不知道如何将它们应用于地图本身。

这是我为创建分类方案所做的：

（其中“值”是我在地图中绘制的列）

然后当我尝试绘制等值线图时，我不知道该方案是什么意思

如果有人可以提供帮助，我将不胜感激！

谢谢 x

让我们说，有两个df，df1和df2。两者都有一列（geometry_1和geometry_2分别）表示线串类型的几何图形。

两个 df 都有更多的行，但现在我想专注于以下问题。有什么方法可以评估两条线是否相似。类似的意思是，如果线的各个点之间的距离不高于有效值（例如100m），则两条线被认为是相同的。

我是python和qgis的初学者。

我尝试使用交集来使两个 shapefile 相交。（多边形）

首先，我尝试了与 QGIS 的交集，它可以工作。(EPSG : 4326) QGIS 图片

但问题是，当我试图在 python 处相交时，只返回一个值。怎么了？？

这是我的代码

有 27000 个值。