我的多边形具有以下结构: 我正在尝试使用spatstat.geom package 将t1转换为 class owin object: t1中的坐标是纬度/经度: 我如何将t1投影到 UTM 才能获得我的owin object?! ...
我的多边形具有以下结构: 我正在尝试使用spatstat.geom package 将t1转换为 class owin object: t1中的坐标是纬度/经度: 我如何将t1投影到 UTM 才能获得我的owin object?! ...
多边形之间的重叠可以使用 SQL Server spatial .STContains function ( https://learn.microsoft.com/en-us/sql/t-sql/spatial-geometry/stcontains-geometry-data-type?view ...
我的问题遵循之前的问题,该问题已得到部分回答: Reading a grib2 dataset with 4 dimensions and 2 variables with R 我正在尝试读取带有 R 的 GRIB2 文件。此文件是概率气象预报:10 个变量、1 个提前期、17 个经度、23 个纬度 ...
R 中的库 (rerddap) 中的 griddap function 提供了一个优秀的 function,用于以简单直接的方式从 ERDDAP 服务器访问卫星数据。 下面是一个可重现的示例,它使用ncdcOisst21Agg_LonPM180数据集提取大西洋的最新海面温度 (NOAA OISS ...
我正在尝试学习(或更具体地说,记住并学习实施)R 中的空间数据分析。我在处理地理数据方面遇到了一些基本困难,这妨碍了任何方法的学习。 切入正题,当我试图产生 map 的财富时,我得到了一些明显奇怪的东西:财富和填充之间没有美学联系,许多国家缺失。 我希望由于部分完整的数据集而丢失一些,但不会像这个那 ...
我有一堆带有方向(罗盘角度从 0-360)和速度(例如公里/小时)的运动矢量,我想将它们汇总为一个平均矢量,但我似乎找不到好办法。 基本上我想能够说“这个组的平均运动是 ___ 度 ___ 公里/小时。” 我找到了一些零碎的方法来做到这一点,但必须有一种更简单的方法。 我正在尝试扩大它,它最终将应用 ...
我正在使用 R 根据学生的原籍城市对学生成绩进行空间分析。我有几个协变量,例如贫困、教育和社会文化指数。 到目前为止,我已经拟合了单变量模型,例如:线性回归、加权线性回归和 CAR(条件自回归)。 现在,我正在阅读 Banerjee、Carlin 和 Gelfand 的“ Hierarchical ...
我有一些中等大小的光栅文件(最大大小 ~190 MB),我想计算使用raster::gridDistance()的网格距离我发现我的最大文件的操作很慢和/或 R 崩溃。 请注意:我不是在寻求 memory 管理建议(例如最大化memory.limit() ,分解成更小的栅格或追求并行处理方法),因为 ...
非常抱歉,但我觉得这个问题没有很好的示例数据集,因此我无法提供可重现的示例。 如果这让这个问题站不住脚那么我一定会删除。 我正在尝试使用 geom_raster() 将 plot 中的 RasterLayer object 添加到 ggplot2,但出现错误。 我已经看到您似乎必须将 RasterL ...
我制作了一个空间 plot,色标从 -100 到 100。但条件是我绘制的值的范围是 -100 到 1500。结果是白色区域(例如红色箭头)描绘的超过 100 的值. 我的问题是如何将值大于 100 的区域概括为与色标 (90-100) 深绿色相同的颜色。 它的脚本: 如果有人可以帮助我,我可以提 ...
我有一个简单的 csv 文件,其中包含 lon/lat/time/values 列: 我想使用签名 'data.frame' rast(x, type="xyz", crs="", digits=6, extent=NULL)的 S4 方法将其转换为 SpatRaster,其中每个“日期”将是一个 ...
我在填充颜色条刻度时遇到问题。 由于我的值在 -85 到 15 之间不是对称的,间隔为 5,所以我希望正值使用红色,负值使用蓝色。 这是我使用的代码: 这是我目前的情节: 我希望你能用 colorRampPalette 为我解释,因为我需要使用这个功能。 谢谢! ...
我正在处理一个使用类似于以下代码片段的数据的问题: 我希望能够选择 X 公里的半径,假设本例为 15,并使用该半径首先计算落在我数据每一行范围内的观测值的数量,其中每一行都作为半径的中心。 然后,我想对半径内的所有观测值取一个度量的平均值。 理想情况下,最终数据集看起来像这样: 这些数字是编造的,但 ...
我在用户表中有超过 20K 个用户,我实现了一次显示 10 个附近用户的查询。 我试过这个 这个查询很慢,如果我做 Explain 那么 我的表结构是这样的 现在为了让它更快,我做了一些研究并了解了spatial index 。 所以我做了这个之后,我添加了名为polygon的新列,并添加了数据类型 ...
现在,当我尝试将值提取到 occ 坐标时,出现以下错误。 经过大量研究后,我没有得到我需要的信息,大多数文章建议转换成我不需要的 csv 格式。 [我试过的链接没有给我任何想法(https://gis.stackexchange.com/questions/403555/raster-subse ...
我想获得每个缓冲区的邻居列表。 但是返回的nb列表是空的。 当我在us$geometry上运行poly2nb()时,一切都很好: ...
我给了一个名称向量列: #a 数据框列表,以 col 的元素作为列,并在特定研究区域之后命名。 (见下文)列表结果包含 19 个不同的文件(不同年份) 具有 6 个元素(对应国家)的 shapefile 列表: 我想从 prevR 库中预制 function: 但实际上对于名称适合的列表中的每 ...
我来自采矿和勘探背景。 我有一个点列表(带有 ID 以及 X 和 Y 坐标)。 例子我有另一个包含多边形(索赔编号)的表。 我试图找出哪些点属于哪些多边形。 我的坐标在 UTM 投影中。 先感谢您! 我正在尝试以下代码,但是我的查询没有返回结果。 我期望从多边形表中获得一个包含 claim_numb ...
这是我第一次运行 WinBUGS。 我的 model 相当简单,只有一个截距和两个随机效应。 但是,当我运行 model 时,我得到了一个完整的错误列表(见下文)。 第一个错误是“未定义的变量”。 我无法确定我的哪些变量在语法中未定义。 任何帮助是极大的赞赏。 这是我的 R 代码。 我正在使用 R ...
我想根据坐标(即纬度和经度)提取每个示例站点的生物地理信息 通过不同的资源搜索,我发现了一个可能有效的代码,不幸的是,我遇到了错误。 我从以下位置获得了生物地理区域的形状文件: https://www.eea.europa.eu/data-and-maps/data/biogeographical- ...