geospatial 相关问题

与地球表面,地上或地下的自然或构造特征和边界的地理位置和特征有关; ESP。指的是地理和空间性质的数据。

如何画椭圆弧?

我正在尝试在图像内的新月周围绘制一个弧形,并遮盖其周围的所有内容。我使用了plotrix包中的函数draw.arc: 情节(1:10,asp = 1,main =“测试绘制...

回答 1 投票 0

R 中按距离对空间点进行分组:如何对空间点进行分组/聚类,使各组之间相隔设定的距离

我有一组 52 个空间点。我试图根据 11 公里的距离将这些点分成几组。即取出所有点并将它们分成组,其中组内的点为 n...

回答 2 投票 0

使用 IndexedPointInAreaLocator 和 STRTree(JTS 拓扑套件)无法实现多边形中的点性能

私人列表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 数据库,或更简单的 在内存存储中建立索引。

回答 1 投票 0

仅预测单个类别的概率(二元分类)

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

回答 1 投票 0

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

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

回答 1 投票 0

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

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

回答 1 投票 0

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

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

回答 1 投票 0

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

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

回答 1 投票 0

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

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

回答 1 投票 0

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

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

回答 3 投票 0

sf 绘图时关闭点位置

我在 csv 中有一系列来自墨西哥的纬度和经度点,我正在将其转换为 sf 对象。我能够使用 http://projfinder.com/ 识别 crs。 图书馆(SF) 库(spData)

回答 1 投票 0

R 传单的局限性?数据太多?

我找到了一种解决方法,实际上对于我正在从事的项目来说效果更好,但我仍然很好奇。 我正在使用 leaflet 包在 R 中创建一些地图。代码运行良好,但我...

回答 2 投票 0

Python PyCharm - 无法上下文安装

我正在使用 Python 和 PyCharm,我正在尝试运行我在这里找到的代码: 从 libpysal 导入权重,示例 从上下文导入 add_basemap 将 matplotlib.pyplot 导入为 plt 导入

回答 1 投票 0

使用 Matplotlib 的抛物线连接图

我正在尝试使用Python创建一个连接图。我希望英国各地的一系列点会聚在一个目标上,并带有突出的抛物线弧。 我有一个代码尝试正在绘制 stra...

回答 1 投票 0

如何将 vue3openlayer 坐标保存到 Django SRID4326

我正在尝试在我的 OpenStreetMap 中制作标记。目前,当我使用 Django 坐标时,它显示了曼谷的一栋建筑。 “SRID=4326;多点(100.53007692708017 13.721325314629256)”...

回答 1 投票 0

为多个组创建动态 Alpha 船体时出现问题

你好吗?我需要你的帮助,在 R 中为不同的组创建 alpha 多边形。我希望这些多边形被我的“物种”列分隔开。 我明白...

回答 1 投票 0

如何使用 terra/sf R 包创建到最近要素的距离栅格?

我有一个由线(道路)组成的简单要素对象,我想从中创建一个栅格,其中每个单元格代表到最近线的距离。就我而言,带有

回答 1 投票 0

alphashape Python 库中用于凹壳计算的 Alpha 参数

我目前正在使用 Python 进行地理空间分析项目,其中我需要计算一组地理点的凹壳(alpha 形状)。我正在为此使用 alphashape 库...

回答 1 投票 0

MySQL 8 上 GeoDjango 的 SRID 问题

GeoDjango 不尊重 MySQL 8 的 SRID。 我创建了一个新专栏 geo_code = PointField(srid=4326) 并运行 makemigrations。在生成的迁移上运行 sqlmigrate 给出 开始; -- -- ...

回答 1 投票 0

如何从我的 Google 地图自动提取数据(KMZ?JSON?)

我正在寻找一种关于如何自动从“我的谷歌地图”检索数据的通用方法。虽然我可以轻松地将任何给定图层导出到 KML/KMZ,但我正在寻找一种方法来在...

回答 1 投票 0

最新问题
© www.soinside.com 2019 - 2024. All rights reserved.