Questions tagged «intersection»

我正在尝试查找多条线与两个不同地理数据框的多边形相交的位置： from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 = LineString([(0.9, 0.9), (0.2, 0.6)]) poly_gdf = gpd.GeoDataFrame(geometry=[polygon]) line_gdf = gpd.GeoDataFrame(geometry=[line1, line2]) 这就是上面的地理数据框的样子（一个具有多边形，另一个具有两条线）。在我看来，两条线都与多边形相交： 但是，相交的输出非常令人困惑： print(line_gdf.intersects(poly_gdf)) 0真实 1错误 print(line1.intersects(polygon)) print(line2.intersects(polygon)) 真正 真正 …

11 python  intersection  shapely  geopandas 

我有一组代表公交线路的线路。一些线重叠并且走相同的路。 我能够提取节点。 但是我对仅提取这样的交叉点感兴趣： 我怎样才能做到这一点？我正在寻找QGIS或Python的方法。 我尝试了GDAL Python 的交集方法，但这基本上只返回了顶点。 如果有两条线交叉，则QGIS的“ 线交叉点”方法会返回交叉点。但是，如果两条公交线路在同一条道路上走得很远，这并不能给我他们指向合并的位置。

10 qgis  python  intersection  line  crossing 

我有两个重叠的多边形层，我需要以某种方式擦除重叠区域。-这样，我应该以两层不重叠的多边形结束。听起来很简单，但是我无法使用Qgis中的任何工具正确完成！！我应该使用“剪辑”或“多边形切割形状”，还是“差异”工具和“剪辑”？ 无论我做什么产品都是不对的。

10 qgis  intersection  overlapping-features 

我在一个区域（浅蓝色）中有一些集线器（黄色）。我需要找到指示该区域中每个点最近的枢纽的多边形。 我使用QGIS获得的voronoi多边形在图像中被涂成棕色。但是，我需要voronoi多边形来填充外部边界，因为我需要多边形来填充整个区域。 在制作voronoi多边形时添加缓冲区，然后使用相交函数将多边形裁剪到区域边界，这使我错了多边形，如下图所示： 如何生成所需的voronoi多边形？

10 qgis  intersection  voronoi-thiessen 

我有两个shapefile集合： 1-邮政编码，行政和政治区 2-各种土地利用类别（例如公园和地段）以及点（例如树） 我想按邮递区号，行政区和政治区以及以相同方式划分的点数对这些土地使用类别的面积进行统计。例如，SE15中有4.1公顷的公园，SE22中有879棵树，依此类推。 我正在使用QGIS。我已经处理了原始数据，以准备好所有shapefile，但是在此网站和其他网站上进行了无济于事的搜索之后，我看不到这样做的方法。

10 qgis  intersection  statistics 

这是一个看似简单的SQL Server几何问题，我以为可以使用一种即用型的解决方案，但是我找不到运气。 我的目的是选择一个表中具有多边形的所有记录，这些记录嵌套（包含）在另一个表的较大多边形中。我有预期的功能STWithin，STContains也有需要提供的解决方案，但是不幸的是，它们都只识别嵌套在较大多边形中的内部多边形，而不识别那些接触较大多边形边界的嵌套多边形。例如参见图片。 满足我的需求的另一个选择是STIntersection。但是，此函数的问题在于它仅返回几何列！我想获取记录ID。有没有人对如何做到这一点有何建议？ STWithin： select a.bg10 from gis.usa_10_block_group a join gis.usa_10_mkt_definition b on a.shape.STWithin(b.shape) = 1 where b.mktname = 'Loop' STContains： select a.bg10 from gis.usa_10_block_group a join gis.usa_10_mkt_definition b on b.shape.STContains(a.shape) = 1 where b.mktname = 'Loop' STIntersection： select a.shape.STIntersection(b.shape) from gis.usa_10_block_group a join gis.usa_10_mkt_definition b on a.shape.STIntersects(b.shape) = …

9 sql-server  intersection 

我尝试了解DE-9IM交集矩阵的空间谓词。我不明白谓语之间的区别封面（各自CoveredBy）和包含（相应地，在）。 我试着用画几何JTS拓扑套件，他们总是Covers = True 和 Contains = True（或两者都是假），但我不管理，使他们不同，即看到的区别Covers和Contains。是否有人举例说明了差异（最好是WKT或图片）？

9 intersection  jts-topology-suite  de-9im 

我之前曾问过一个问题，该问题用于确定具有另一个多边形边界的多边形的面积。 确定多边形边界内每个多边形的面积？ 我猜想，使用相交函数是确定1000m边界内多边形面积的最可行选择。每当我与“多边形”边界内的所有多边形shapefile（大约有40个）相交时，都会生成一个空输出。 我已将每个shapefile和数据框的坐标系设置为相同的坐标系。我正在使用具有ArcInfo许可证的ArcGIS 10.1。 产生此空输出的原因是什么？ 警告000117：警告生成空输出。

9 arcgis-desktop  arcgis-10.1  polygon  intersection  warning-000117 

我正在尝试使用PostgreSQL / PostGIS计算矢量层中每个多边形的栅格统计信息（最小值，最大值，平均值）。 这个GIS.SE答案描述了如何通过计算多边形和栅格之间的交点然后计算加权平均值来做到这一点：https : //gis.stackexchange.com/a/19858/12420 我正在使用以下查询（dem我的栅格在哪里，topo_area_su_region我的矢量在哪里，并且toid是唯一的ID： SELECT toid, Min((gv).val) As MinElevation, Max((gv).val) As MaxElevation, Sum(ST_Area((gv).geom) * (gv).val) / Sum(ST_Area((gv).geom)) as MeanElevation FROM (SELECT toid, ST_Intersection(rast, geom) AS gv FROM topo_area_su_region,dem WHERE ST_Intersects(rast, geom)) foo GROUP BY toid ORDER BY toid; 这可行，但是太慢了。我的矢量层2489k的特点，周围每个人服用90毫秒到过程-这将需要几天来处理整个层。如果仅计算最小值和最大值（这避免了对ST_Area的调用），则计算速度似乎并没有明显提高。 如果我使用Python（GDAL，NumPy和PIL）进行类似的计算，则可以显着减少处理数据所需的时间，而不是对栅格进行矢量化处理（使用ST_Intersection），可以对矢量进行栅格化。在此处查看代码：https : //gist.github.com/snorfalorpagus/7320167 我真的不需要加权平均值-“如果碰到就行了”的方法就足够了-而且我有把握地确定这是放慢速度的原因。 问题：有什么方法可以使PostGIS像这样运行？即从多边形所接触的栅格中返回所有像元的值，而不是返回精确的交集。 我是PostgreSQL / PostGIS的新手，所以也许还有其他我做不正确的事情。我正在Windows …

9 postgis  postgresql  intersection 

我正在尝试从Poly2删除到Poly1的相交。以下是2个多边形的坐标。 > coords1 [,1] [,2] [1,] 992.0161 7462.531 [2,] 950.1962 7413.532 [3,] 902.3632 7421.207 [4,] 706.1985 7472.378 [5,] 654.0139 7711.760 [6,] 657.5960 7726.950 [7,] 786.2667 7844.380 [8,] 789.5935 7842.925 [9,] 981.7046 7498.659 [10,] 983.4246 7493.271 [11,] 990.6680 7469.888 [12,] 992.0161 7462.531 > coords2 [,1] [,2] [1,] 930.3464 7607.313 [2,] 979.3528 …

9 r  polygon  intersection  symmetricdifference 

我使用ENVI EX进行了基于对象的图像分析，并将每个分类特征导出为shapefile。我将这些shapefile转移到ArcMap 10中，并尝试将它们与从GeoBase获得的相同区域的其他shapefile相交。 当我尝试运行相交工具时，收到一条错误消息，指出从ENVI EX提取的shapefile的拓扑无效。我试图打开拓扑工具来修复这些shapefile的拓扑，但是按钮显示为灰色，我似乎无法单击它们。拓扑工具不适用于我拥有的任何shapefile。 有什么办法吗？如何修复拓扑并运行相交？

9 arcgis-desktop  geometry  intersection  topology  envi 

如何在ArcGIS 10中区分这些情况？ 情况1：一条线的两个端点都碰到另一条线 情况2：两个端点都悬在其相交的线上 我正在看修剪线功能，但这不是我想要的（破坏性的）。 现实中的用例是区分连接高速公路的两条道路的服务道路和与高速公路相交的道路的其他情况。

9 arcgis-desktop  arcgis-10.0  intersection 

我需要编写PL / pgSQL存储过程来计算多层的相交。例如，对于三层A，B，C，此函数应返回ABC几何： 函数将需要相交的图层的id作为输入。谁能给我建议以构建此功能？我对此的想法是，我需要先相交，然后相交，然后相交的结果与第三层相交，等等。

9 postgis  postgresql  intersection