问题标签 [sp]

我正在尝试使用 Package graphics中的通用函数绘制一个 SpatialPoints 对象。图的范围应该正好是最外层数据点的大小。这是一个简单的示例，其中 delta x 与 delta y（最外层值）的比率为 1：

据我了解，上述脚本生成的图应该生成一个图形，其中只有 1 个点可见（坐标 1.5,1.5），但所有 4 个“斜线”都应该是不可见的。为什么不是这样？

我有几个数据集（大约 25 个具有所有相同的列），总共有大约 400,000 个观察值。我正在运行一个循环来提取每个数据调查点的道路类别。当我在一个小数据集（~ 4,000 obs）上运行我的脚本时，一切正常。但是，当我使用相同的脚本运行循环时，在前 47,000 个 obs 之后出现错误（停止循环）。我用谷歌搜索了那个错误，没有发现任何结果。有什么建议么？

谢谢马丁

.Call("rgeos_intersects_prepared", .RGEOS_HANDLE, spgeom1, spgeom2, 中的错误：不允许负长度向量

这是我用于分析的部分代码

我想创建扇区，然后为它们着色，下面的代码确实可以，但我还想添加一个图例来表示从提供的色调中选择的颜色光谱。它使用 colour.heat 中的颜色和为颜色提供的颜色始终介于 0 和 1 之间，作为其比率。

因此，当我运行以下代码时，它会绘制组成一个圆圈的扇区。

如何创建一个图例来表示颜色并将其添加到图表中？

我已经尝试过，但不幸的是无法添加图例。

我有一个在不同日期收集的空间多边形数据框列表。我想在列表中的每一天添加时间，这样我就可以对“Z”（DDD.LFRP）变量进行子集化和基于时间/日期的方法。以下是空间多边形列表的子集：

请问如何添加添加时间/日期？时间戳是从文件名中的一组字符串子集创建的：

然后我尝试了@RobertH 编写的代码

我在将此空间多边形数据框转换为栅格时遇到问题。当我进行转换时，栅格将 NA 作为其值。如下所示：

这就是我栅格化它的方式：

我做错了什么？我需要一种方法来确保数据框中的值反映在栅格中，请。

我有一个 0.25 * 0.25 度网格分辨率的数据集。我有另一个是 1 * 1 度分辨率，并希望使它们具有可比性（都更改为 1 * 1 分辨率）。该区域是 28.5 到 36.5 经度和 -4.5 到 4.5 纬度，2003 年到 2012 年的月度数据。我附上了前几行，例如，如果以下内容没有帮助，我附上了完整数据的 excel 文件！https://onedrive.live.com/redir?resid=9E74848E574367C6!2625&authkey=!AOcpftTyHVqU8CA&ithint=file%2ccsv我一直在研究“sp”（聚合）和“光栅”包，但无法弄清楚。一如既往地感谢任何帮助！

我正在使用 R 处理空间数据集。

资料说明

我的主数据集是SpatialPointsDataFrame格式的，并且每个月都有地表温度数据（列名 - “ruralLSTday”、“ruralLSTnight”）。数据片段如下所示：

主数据 - （SpatialPointsDataFrame 格式）

坐标：

我想在上述数据中添加有关降雨量和气温的列 - 该数据存在于SpatialGridDataFrame每个月的“secondary_data”表中。“secondary_data”的片段如下所示：

辅助数据 - （以 SpatialGridDataFrame 格式）

坐标：

问题

如何通过匹配纬度经度和月份将辅助数据中的列添加到我的主数据中？目前上述两个表格中的纬度/经度信息不会完全匹配，因为主数据是一组点，辅助数据是网格。

有没有办法在我的主数据的纬度/经度落入的“辅助数据”上找到网格的正方形并进行插值？

我有一个EuropeanCities纬度和经度的数据框

目前这些值属于数字类，但是我需要将它们设置为空间数据的空间坐标

如何使用该sp函数转换这些值的类？

下面是一个非常有用的例子。我已经拼凑了在 ggplot 中创建的带有插图的地图的代码。

这是我遇到问题的地方。使用 grid.newpage 和 viewport 我的电脑在打印插图时会停止。我究竟做错了什么？您会注意到，我什至努力降低插入插图的国家形状文件的分辨率。

提前致谢。

我有一个数据框，其中包含韩国首尔（首都）空气中 PM10 的浓度。请看一下。我想从这个数据集中绘制半变异函数。由于这里的 LAT/LON 数据是度数，所以我对这些数据进行了投影。我以这种方式预测数据：

投影后我得到了seoul311，如下所示

Q1：我发现投影后，LON/LAT 的值与之前的数据帧显示的值几乎相同。我的问题是这个proj4string(seoul311) = "+proj=longlat +datum=WGS84"命令的实际功能是什么。在这里，LON/LAT(degree) 转换为 km/m 或类似的东西？

我尝试使用 rgdal 包编写另一个代码，如下所示：

投影后我得到了seoul3112，如下所示

Q2。在这里，我可以看到 LON/LAT 值转换为某个较大的值！这些值的含义是什么？米/公里或类似的东西？在上面的代码中，北意味着什么？北半球？

Q3。正如我之前提到的，我想在韩国首尔（utm zone 52）绘制半变异函数。那么，我应该使用哪种投影规则？我应该考虑utm区吗？我什么时候应该考虑utm区？

我对投影数据有很多困惑。你能详细回答我的三个问题吗？