cost 64 ms
当经度为 c(-180, 180) 但 griddap 数据集经度为 0-360° 时,如何使用 library(griddap) 中的 rerddap function 获取数据? - How to get data using the rerddap function in library(griddap) when longitude is c(-180, 180) but the griddap dataset longitude is 0-360°?

R 中的库 (rerddap) 中的 griddap function 提供了一个优秀的 function,用于以简单直接的方式从 ERDDAP 服务器访问卫星数据。 下面是一个可重现的示例,它使用ncdcOisst21Agg_LonPM180数据集提取大西洋的最新海面温度 (NOAA OISS ...

使用 geom_sf() 绘制地理数据 - Plotting geographic data with geom_sf()

我正在尝试学习(或更具体地说,记住并学习实施)R 中的空间数据分析。我在处理地理数据方面遇到了一些基本困难,这妨碍了任何方法的学习。 切入正题,当我试图产生 map 的财富时,我得到了一些明显奇怪的东西:财富和填充之间没有美学联系,许多国家缺失。 我希望由于部分完整的数据集而丢失一些,但不会像这个那 ...

在 R 中找到运动矢量的平均轨迹 - find average trajectory of movement vectors in R

我有一堆带有方向(罗盘角度从 0-360)和速度(例如公里/小时)的运动矢量,我想将它们汇总为一个平均矢量,但我似乎找不到好办法。 基本上我想能够说“这个组的平均运动是 ___ 度 ___ 公里/小时。” 我找到了一些零碎的方法来做到这一点,但必须有一种更简单的方法。 我正在尝试扩大它,它最终将应用 ...

MCAR(多变量条件自回归)model in R - MCAR (Multivariate Conditional Autoregressive) model in R

我正在使用 R 根据学生的原籍城市对学生成绩进行空间分析。我有几个协变量,例如贫困、教育和社会文化指数。 到目前为止,我已经拟合了单变量模型,例如:线性回归、加权线性回归和 CAR(条件自回归)。 现在,我正在阅读 Banerjee、Carlin 和 Gelfand 的“ Hierarchical ...

raster::gridDistance() 在大中型光栅文件上失败 - raster::gridDistance() fails on medium to large raster files

我有一些中等大小的光栅文件(最大大小 ~190 MB),我想计算使用raster::gridDistance()的网格距离我发现我的最大文件的操作很慢和/或 R 崩溃。 请注意:我不是在寻求 memory 管理建议(例如最大化memory.limit() ,分解成更小的栅格或追求并行处理方法),因为 ...

根据 R 中的值更改色标 - Change color scale based on value in R

我制作了一个空间 plot,色标从 -100 到 100。但条件是我绘制的值的范围是 -100 到 1500。结果是白色区域(例如红色箭头)描绘的超过 100 的值. 我的问题是如何将值大于 100 的区域概括为与色标 (90-100) 深绿色相同的颜色。 它的脚本: 如果有人可以帮助我,我可以提 ...

使用 rast(..., type=xyz) 将 data.frame 转换为 terra 中的 SpatRaster 时指定图层 package (R) - Specifying layers when using rast(..., type=xyz) to convert data.frame to SpatRaster in the terra package (R)

我有一个简单的 csv 文件,其中包含 lon/lat/time/values 列: 我想使用签名 'data.frame' rast(x, type="xyz", crs="", digits=6, extent=NULL)的 S4 方法将其转换为 SpatRaster,其中每个“日期”将是一个 ...

如何计算地图上半径范围内的度量平均值 - How to Calculate Avg of Metric within Radius on Map

我正在处理一个使用类似于以下代码片段的数据的问题: 我希望能够选择 X 公里的半径,假设本例为 15,并使用该半径首先计算落在我数据每一行范围内的观测值的数量,其中每一行都作为半径的中心。 然后,我想对半径内的所有观测值取一个度量的平均值。 理想情况下,最终数据集看起来像这样: 这些数字是编造的,但 ...

使用空间索引mysql获取附近的用户 - get nearby users using spatial index mysql

我在用户表中有超过 20K 个用户,我实现了一次显示 10 个附近用户的查询。 我试过这个 这个查询很慢,如果我做 Explain 那么 我的表结构是这样的 现在为了让它更快,我做了一些研究并了解了spatial index 。 所以我做了这个之后,我添加了名为polygon的新列,并添加了数据类型 ...

如何从堆叠栅格图层(asc 格式)中提取值以出现(csv)数据? - How to extract values from stacked raster layers(asc format) to occurence(csv) data?

现在,当我尝试将值提取到 occ 坐标时,出现以下错误。 经过大量研究后,我没有得到我需要的信息,大多数文章建议转换成我不需要的 csv 格式。 [我试过的链接没有给我任何想法(https://gis.stackexchange.com/questions/403555/raster-subse ...

如果 df 列表的名称匹配,则在 Function 上使用 Map - Using Map on Function if name of df lists match

我给了一个名称向量列: #a 数据框列表,以 col 的元素作为列,并在特定研究区域之后命名。 (见下文)列表结果包含 19 个不同的文件(不同年份) 具有 6 个元素(对应国家)的 shapefile 列表: 我想从 prevR 库中预制 function: 但实际上对于名称适合的列表中的每 ...

空间加入SQL服务器 - Spatial Join in SQL Server

我来自采矿和勘探背景。 我有一个点列表(带有 ID 以及 X 和 Y 坐标)。 例子我有另一个包含多边形(索赔编号)的表。 我试图找出哪些点属于哪些多边形。 我的坐标在 UTM 投影中。 先感谢您! 我正在尝试以下代码,但是我的查询没有返回结果。 我期望从多边形表中获得一个包含 claim_numb ...

WinBUGS 错误 - 我的“未定义变量”是什么? - WinBUGS Error - What is my "undefined variable"?

这是我第一次运行 WinBUGS。 我的 model 相当简单,只有一个截距和两个随机效应。 但是,当我运行 model 时,我得到了一个完整的错误列表(见下文)。 第一个错误是“未定义的变量”。 我无法确定我的哪些变量在语法中未定义。 任何帮助是极大的赞赏。 这是我的 R 代码。 我正在使用 R ...


 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM