问题标签 [geosphere]

我有研究站点，我收集了数据和附近的气象站，其中包含有关温度和降水的信息。我想将我在研究地点的日常数据与最近气象站的天气信息配对。我认为，要做到这一点，我需要一个两步过程，首先选择离研究地点最近的气象站，然后使用天气数据创建一个新变量。

这是我的数据的快照：

我正在尝试使用这条线来确定最近的车站，但我只得到一排距离。

计算完所有距离后，我对接下来的步骤有点不确定，但我想我需要一些东西来选择最近的车站和匹配日期。这是我想出的最好的：

最终结果是我的研究站点数据框，其中包含最近气象站的温度附加列。

我有两个数据集 A 和 B，它们给出了英国不同点的位置：

A 有 400 行，B 有 1800 行。
对于 A 中的所有行，我想找到 A 中的一个点与 B 中三个最近点中的每一个之间的最短公里距离，以及 B 中这些点的经纬度的参考和坐标。

我尝试使用这篇文章

R - 查找给定半径内最近的相邻点和相邻点的数量，坐标经纬度

但是，即使我按照所有说明进行操作，主要是使用distmpackage中的命令geosphere，距离仍以不可能为公里的单位出现。我看不出代码中有什么要改变的，特别是因为我对这些geo包一点也不熟悉。

我有一个带有两个邮政编码和相应的纬度和经度的大型数据集（2.6M 行），我正在尝试计算它们之间的距离。我主要使用该包geosphere来计算邮政编码之间的文森蒂椭圆体距离，但我的数据集需要大量时间。什么是实现这一点的快速方法？

我试过的

这两种方法都非常缓慢。我知道 210 万行永远不会是一个“快速”的计算，但我认为它可以做得更快。我在较小的测试数据上尝试了以下方法，但没有任何运气，

任何人都可以帮助我以正确的方式使用doParallelwithgeosphere或更好的方式来处理这个问题吗？

编辑：来自（一些）回复的基准

一个简单的all.equal测试表明，对于我的数据集，haversine 方法等于 vincenty 方法，但与geodist包中的“便宜”方法具有“平均相对差异：0.01002573”。

我试图用来geosphere计算两个大圆之间的交点，这些圆以这样的数据框格式给出：

我尝试使用将矩阵作为输入的gcintersect函数。geosphere为了在我使用的数据框中使用它cbind，但似乎 mutate 不适用于此：

错误可能是由于输出比预期的要长（数据帧的行数）。所以我试着把它放在一个列表中：

再次在这里，我看到输出是所有组合，而不是获取每行 2 个交叉点的 4 个坐标。

预期的输出将是新列中单元格中的列表，其中包含两个点的坐标。

是否有不使用循环（或purrr）的解决方案，因为这将明显慢于mutate.

第一行中 int_points 的值应与以下输出相同：

我想确定一些地址的地理距离，并确定这些地址的平均值（平均距离）。

如果数据框只有一行，我找到了一个解决方案：

不幸的是，我没有实现一个为整个数据集执行代码的函数。我目前的问题是，该函数不会按行计算距离平均值，而是计算数据集所有行的平均值。

你知道我犯了什么错误吗？

澄清我的问题：我正在尝试创建一个像这样的数据框（V5）：

例如，每行距离的平均值。

我正在尝试从起点和向南和向东行驶的距离计算纬度 lng 点。我将使用什么地圈函数来计算？

我试过从互联网上手动计算，但似乎没有什么对我有用。我确保使用弧度。

52.28588 -101.49886

认为它应该在附近

57.760 -99.764

（根据谷歌地图推测东 300 公里 250 公里。）

我想通过替换两个最近车站的温度平均值来估算车站温度的缺失值。

这用所有站点的非缺失温度的平均值替换缺失值。

使用此代码查找两个站点之间的距离

无法弄清楚如何使用 tidyverse 计算距离？

.pointsToMatrix(x) 中的错误：向量的长度错误，应为 2

我正在做一个地图项目，并在传单地图（闪亮的应用程序）中绘制了两个标记数据集。一个数据集比另一个大，我正在绘制连接多到少的大圆弧。一切正常。我将 R 从 3.5 更新到 3.6，现在我得到“如果错误：缺少 TRUE/FALSE 需要的值”，我已将其隔离到运行大圆弧的 for 语句中。

硬编码到 gcIntermediate() 中的值可以工作，但我只能有 1 个固定弧。使用值创建一个新的数据框是可行的，但是从我导入的数据构建的 data.frame 不再有效。

上面的代码有效，作为 NetData 提供的数据是导入的 NetData 的 head(10)。num 和 chr 方面与单独留在该庄园中时的正确文件相同。当我不硬输入它并将 100 行数据输入 for 循环时，我收到警告：如果错误：需要 TRUE/FALSE 的缺失值。更具体地说，我得到：

第 136 行是在 server.R 中调用 arcbuilder() 的地方，第 193 行是在 global.R 中的 arcbuilder() 中使用 gcIntermediate() 的地方

最后，当将整个 ​​NetData data.frame 运行到带有地图设置以在闪亮之外工作的 arcbuilder 函数中时，它运行没有错误。我检查了来自 NA 或其他异常的数据，我有理由相信没有任何数据。

我错过了什么？

I am trying to calculate the area of a polygon defined in terms of latitude and longitude:

using the function areaPolygon from the geosphere package:

However, I have noticed that if I reorder the points in the polygon I get a different result:

My question is what is the reason for this disparity in the results (the function documentation does not mention that ordering matters)? Is there a correct way to order the points in a polygon, and what is it?

我有一个包含欧洲国际投资信息的数据集，并协调了有关 NUTS3 的信息。对于每项投资，我都有城市和坐标（lat1，long1）。我想计算每个城市到我拥有的每个 NUTS 3 的距离——> EG Paris 到 Paris、Paris_Lyone、Paris_Orly、Paris_Maidenhead 等。我想为我拥有的所有城市循环这个机制，所以最后我对每个城市都有一个矩阵，其中包括到每个 NUTS 的距离。我尝试使用 geosphere，但它只给了我行之间的距离。

这不起作用，因为它逐行返回我的距离，但我实际上想要从每个城市到我拥有的所有 NUTS3 坐标的距离

有人可以帮我吗？

我不确定如何发布我的 dt，我认为这可能有助于获得更多建议。