geospatial 相关问题

私人列表findZonesForShipPosition（点点）{ List zoneDefinitionCacheHits = new ArrayList<>(); strTree.query(point.getEnvelopeIn... private List<ZoneDefinitionCacheMemory> findZonesForShipPosition(Point point) { List<Object> zoneDefinitionCacheHits = new ArrayList<>(); strTree.query(point.getEnvelopeInternal(), zoneDefinitionCacheHits::add); return zoneDefinitionCacheHits.stream() .filter(ZoneDefinitionCacheMemory.class::isInstance) .map(ZoneDefinitionCacheMemory.class::cast) .filter(zoneDefinition -> isPointInsideGeometry(point, zoneDefinition.getGeometry())) .toList(); } private boolean isPointInsideGeometry(Point point, Geometry geometry) { PointOnGeometryLocator pointOnGeometryLocator = new IndexedPointInAreaLocator(geometry); int location = pointOnGeometryLocator.locate(point.getCoordinate()); boolean isInside = location == Location.INTERIOR; return isInside; } 我使用上面的代码来搜索多边形中的点，在查询之前将所有多边形和多边形插入到 strTree 中，并且我还使用 IndexedPointInAreaLocator 在 STR Tree 的顶部添加一个附加过滤器来验证是否给定点是否确实位于几何中。此外，我的多边形非常大，这意味着某些多边形中有近 160 万个点。因此，通过此实现，我面临着性能问题，这种方法至少需要 100 毫秒才能完成，这对于我的用例来说是不期望的。这里可能出现什么问题？或任何以某种方式加快速度的建议。我的用例：我将获得很多点来检测它们是否位于给定几何形状（多边形和多多边形）的内部或外部我尝试使用 Polygon.contains(point) 谓词，但它比 IndexedPointInAreaLocator 方法花费更多时间。如果你的多边形很大（从某种意义上说它们覆盖了很大的区域）或者形状很奇怪。然后，细分它们可能会提高索引性能，请记住索引首先在多边形的边界框上工作，以节省在多边形测试中进行昂贵的点。如果一个点落入许多边界框但只有一个多边形，则索引不会有太大帮助。 Paul Ramsey 的这篇博客文章在 PostGIS 上下文中讨论了这一点，但这里的底层实现是相同的。这个answer将帮助您划分多边形（因为我认为JTS没有细分方法）。您还可以从 JTS 代码周围使用 GeoTools 中受益，因为这将允许您切换到使用具有内置索引的 PostGIS 或 GeoPackage 数据库，或更简单的在内存存储中建立索引。

geospatial point-in-polygon jts

回答 1 投票 0

仅预测单个类别的概率（二元分类）

我正在为数百个物种生成物种分布模型，其中包含 45 个特征的 vrt 堆栈，总计约 110GiB。这些模型以 0-1 的范围预测合适栖息地的概率，

r geospatial predict terra

回答 1 投票 0

将 JTS 几何扩展到北极/南极

我正在开发一个项目，需要将任意几何图形扩展任意米数。我的解决方案基于 Geotools:bounding box for a buffer in wgs84 的答案。我...

java geospatial geotools jts

回答 1 投票 0

项目光栅错误，包 terra 1.7.55 xy 不是数字，文档中的示例代码产生相同的错误

我不知道如何解决这个问题，有人有什么建议吗？我在投影光栅时遇到问题。我收到消息警告：PROJ 支持由 sf 和 terra 包等提供...

geospatial spatial projection terra

回答 1 投票 0

多边形三角剖分不适用于使用 Delaunay 三角剖分方法的某些多边形

我原来的多边形WKT是这样的：多边形 ((98.51786360976723 78.60265688535823, 98.51786360976723 78.67546052678281, 98.78608990953424 78.67546052678281, 98.7860 8990953424 78.60265688535823, 98.

geospatial triangulation jts

回答 1 投票 0

PostGIS for PostgreSQL 中 ST_Distance 函数使用的算法名称？

我正在使用 PostgreSQL 的 PostGIS 扩展，并使用 ST_Distance 函数来计算两点之间的距离。我很好奇 ST_Distance 我们的底层算法...

postgresql gis geospatial postgis

回答 1 投票 0

在雪花中进行地理空间搜索的最快方法是什么？

我是地理空间领域和数据的新手，最近需要地理空间搜索，我需要将其作为模型中的功能。这个想法是对于给定的点 A，如何找到 20 英里内最近的 10 个点。基于

sql optimization snowflake-cloud-data-platform geospatial

回答 1 投票 0

如何计算R中的土地覆盖面积

基本上，我以 ASCII 的形式计算了一个全局分布概率模型，比如：国内生产总值。 gdpm的值都在0到1之间。然后我从形状文件导入了本地地图：形状 <- file.