其他开发

在R(Cran)中创建和裁剪网格

我需要在gstat库中创建一个用于空间内插的遮罩网格。具体地说，我有不同的随机分布的采样点，我需要创建包围这些点的最小凸多边形。然后，我必须创建一个空间网格，它应该被计算的外壳裁剪，以便将内插限制在这个多边形的范围内。如果有人能给我解释一下详细的程序并提供一些例子，我将不胜感激。提前谢谢您。 推荐答案 我自己找到了解决方案。 library(spatstat) library(s ..
发布时间：2022-05-29 20:39:28 其他开发

在R中创建空间数据

我有一个100x200米区域内物种及其大致位置的数据集。数据框的位置部分不是我认为可用的格式。在这个100x200米的矩形中，有两个名为A到CV的10x10米正方形。在每个10x10正方形内有四个5x5米的正方形，分别命名为1、2、3和4(1在2的南边，3.4的西边，4的东边，3的北边)。我想让R知道，A是在(0，0)、(10，0)、(0，0)和(0，10)处有角的正方形，B就在A的北边，有角(0 ..
发布时间：2022-05-29 20:18:28 其他开发

在R中，如何在空间网格正方形上平均空间点数据

现在设法解决了问题 我有一个大约50,000个点的集合，它们有坐标和一个与之相关的值。我希望能够将点放入网格中，平均所有落入网格正方形的点的关联值。所以我想以一个对象结束，该对象标识每个网格正方形，并给出网格正方形内部的平均值。 我有一个空间点数据框和一个空间网格对象中的数据，如果这有帮助的话。 改进答案：我确实做了一些搜索，很抱歉问题的初始状态我只是在自己的脑海里构思了这个问题 ..
发布时间：2022-05-29 20:01:08 其他开发

用SPPLOT叠加绘制两个空间多边形数据帧

我有一堆数据，我已经在县级绘制了without borders。我想加上一些州的边界。我有一个州形状文件(多边形)，但spplot似乎没有任何方法可以添加到以前的地图上。有没有办法在不重写面板函数以获取两个SPDF的情况下做到这一点(对于其他人可能遇到的问题，这似乎是非常专门的)？ 这里有一个可重复使用的示例： library(sp) Srs1 = Polygons(list(Poly ..
发布时间：2022-05-29 19:50:20 其他开发

如何推算空间权重矩阵的缺失邻居(皇后邻接性)

我有一个很大的形状文件，大约有180.000个250m^2的多边形。 我想创建一个空间权重矩阵(皇后邻接性)。 因此，如果是邻居，则为1，否则为0。 但是，有几个没有任何相邻(岛屿)的多边形。 我如何才能为那些没有任何直接邻居的单元分配最近的邻居？ (使用SF或SP包-在R中) ................................................... ..
发布时间：2022-05-29 19:37:09 其他开发

有没有一种方法可以删除ggploy外部的空格？

我使用代码创建一个静态地图，并将其包含在一个闪亮的应用程序中。我有一个问题，那里的地块大小比实际的地块大得多，所以布局看起来不太对劲。有什么方法可以减少绘图外部的空白吗？ 如果可能，我希望输出仅为绘图本身。 # reprex library(terra) library(sf) library(magrittr) library(ggplot2) lux ..
发布时间：2022-05-29 19:17:29 其他开发

SQL几何查找半径内的所有点

我精通SQL，但不熟悉使用SQL几何特性。我有一个可能是非常基本的问题要解决，但我还没有在网上找到任何好的资源来解释如何使用几何对象。(TechNet是一种学习新事物的糟糕方式...) 我在笛卡尔平面上有一个二维点的集合，我正在尝试查找半径集合内的所有点。 我使用如下语法创建并填充了一个表： 更新[事物]set[位置]=几何：：Point(@X，@Y，0) (@X，@Y只是 ..
发布时间：2022-05-29 18:33:30 其他开发

SpatialPolygonsDataFrame对象的坐标()值？

我正在尝试获取空间多边形数据框中多边形的伪重心。今天我偶然发现了coordinates函数，它实际上为SpatialPolygonsDataFrame返回了一些东西。 不幸的是，我在coordinates的帮助下没有找到有关SpatialPolygonsDataFrame的值。谁能告诉我这些坐标是什么？ 推荐答案 它是多边形质心。找到源代码here，查找函数FindCG。计算的方 ..
发布时间：2022-05-29 18:26:51 其他开发

如何将多边形修改为更改其槽的洞(SpatialPolygons)

编辑： 正如Edzer Pebesma在评论中所建议的，向多边形添加洞的推荐方法不是修改槽，而是重建多边形，如in this related question所示。 原始问题 在SpatialPolygons-class的帮助下，我试图将一个多边形修改为另一个多边形洞，但当该“洞多边形”显示为另一个多边形的边界时，其内部的颜色与其他多边形相同。 我做错了什么？ 使用d ..
发布时间：2022-05-29 18:18:41 其他开发

计算地理质心的经度/纬度

我想对内布拉斯加州的县级作物产量数据进行一些空间统计分析，用于STAT类。为此，我需要每个县的地理质心的经度和纬度。有人知道在R里怎么做吗？我知道它可以在ArcGIS中完成，但我现在无法访问它。 推荐答案 您没有提供任何详细信息，您的shapefile是从here获得的，但我从here获得了一个，您可以使用gCentroidrgeos： library(rgdal) librar ..
发布时间：2022-05-29 18:13:02 其他开发

在R中使用mat2listw函数创建空间权重矩阵

我正在尝试使用mat2listw函数在R中创建一个Weights对象。我有一个非常大的空间权重矩阵(大约22,000 x 22,000) 它是在Excel中创建的，并读取到R中，现在我正在尝试实现： library(spdep) SW=mat2listw(matrix) 我收到以下错误： Error in if (any(x ..
发布时间：2022-05-29 18:04:49 其他开发

R消息：错误-替换有x行，数据有y行

我正在尝试使用ggmap包来计算具有地址列表的目标地址的距离。我在CSV中的数据如下： Order ID Address 1652049 435 E 70TH ST,10021 1652123 1365 YORK AVE,10021 1652053 530 E 72ND ST,10021 所以我尝试获得从我的输入地址到所有这些地址的距离，例如：哈德逊街400号，10013，我 ..
发布时间：2022-05-29 18:01:36 其他开发

查找从数据框到特定位置的最近城市

下面的DataFrame包含有关LATE、LONG、STORY和CITY的信息。我想要找到 数据框中给出的每个城市对应的三个最近城市。例如，从下面的 离阿尔伯克基最近的DataFrame、俄克拉何马城和科罗拉多州斯普林德，所以离阿尔伯克基最近的三个城市应该是 保存在名为NEAREST_AL的其他数据帧中(我不知道如何获得这个结果，所以我试图通过创建数据框来给出一个想法)。 dataframe ..
发布时间：2022-05-29 17:53:41 其他开发

联合SQL Server表中的所有几何，如Postgres中的GeomUnion

先说明一下：我说的是联合几何图形，而不是SQL关键字UNION。 我正在尝试将一些空间数据从带有PostGIS的Postgres迁移到SQL Server2008。一切都很好，直到我看到这样的声明： SELECT GeomUnion(the_geom) FROM some_table 这将合并该列中的所有几何图形，并将其作为一个结果返回(类似于COUNT的工作方式)。据我所知，SQ ..
发布时间：2022-05-29 17:45:34 其他开发

在R中创建等距空间栅格

我正在尝试使用R在给定的区域上创建一个等维度的正方形网格。我希望我的网格是1千米x 1千米的正方形。我看到了这样的示例，它们说明了相等的经度/经度网格： Creating a regular polygon grid over a spatial extent, rotated by a given angle 但这甚至不是大小。看起来我应该能够使用st_make_grid函数来创建这 ..
发布时间：2022-05-29 17:37:43 其他开发

在R中按邻近度进行空间滤波

我有一个物种的发生点，我想消除潜在的抽样偏差(其中某些地区的点密度可能比其他地区大得多)。要做到这一点，一种方法是最大化彼此不小于某一距离X的点的子集。从本质上讲，我会防止点之间过于接近。 是否有任何现有的R函数可以执行此操作？我已经搜索了各种空间包，但没有找到任何东西，也不知道自己到底如何实现。 可以下载发生点数据集示例here。 谢谢！ 推荐答案 我已经编写了此函 ..
发布时间：2022-05-29 17:29:59 其他开发