问题标签 [rgdal]

我有一些地理数据点（纬度和经度），我想为它们提取 1 公里分辨率的气候数据（最大值、最小值、平均温度）。这里有一个很好的编码示例，说明如何在 R 中执行此操作：

问题是它需要“rgdal”包，不幸的是，这似乎不适用于我正在使用的 Mavericks。我尝试了在此处安装“rgdal”的解决方法，但没有成功。

有人知道在 Mavericks 上安装 rgdal 的解决方法吗？或者任何熟悉在 R 中使用 Worldclim 并知道另一个包或其他解决方法的人？

这是站点的标识符，纬度/经度：

我正在尝试从顶点位置创建一组多边形，以 X，Y 格式保存。

这是我的数据示例 - 每行代表一个多边形的顶点。多边形是正方形

我正在使用SpatialPolygons，因此我的数据需要在列表中。所以我创建了一个循环来尝试将我的数据从矩阵中转换为列表格式。

我在此站点上的其他一些问题中找到的代码之后创建了一个循环。我打破了每一步，试图理解为什么我只有一个多边形作为输出，即使我有 2 组点。

你能帮我理解我如何调整代码来写出两个多边形而不是一个多边形吗？而且，我如何将 ID 分配给每个多边形，因为我使用矩阵（正方形）作为我的起始数据集，如果我分配一个字符 ID，它会将我的所有数据转换为一个字符。

我的最终目标是SpatialPolygons对象中有两个多边形，第一个带有 ID SJER1，第二个带有SJER2存储在SpatialPolygons对象中的 ID。

然后我会把它写到一个 shapefile 中。

我正在尝试使用 de raster 包来读取具有 160 个反射率值和 160 个每个像素波长值的多层（多波段）图像（ENVI 格式 [.hdr]），但是当我使用我开发的代码时，程序仅返回 1 个波段和相关的反射值。section1=raster("./x") getValuesBlock(section1, row=1, nrows=1, col=1, ncol=1 )

我正在尝试使用开放数据（避免像谷歌这样的许可限制）计算远足路线的海拔数据。

我能够使用 readGDAL（来自 RGDAL 包）读取我国的公共 DEM（分辨率为 10 米），并且 proj4string(mygrid) 给了我：

.asc 文件的开头是：

我需要做的就是从这个网格中获取路线特定节点的高程数据，以便能够计算高程增益、负斜率、最小/最大高度......

我使用漂亮的包 OSMAR 从 OpenStreetMap 带来路线数据，所以我的路线的数据表是这样的：

但我不知道如何在 DEM 坐标参考系中转换纬度/经度坐标，然后如何带来相应的网格值（对最近点进行某种平均？）

我在谷歌上搜索到的所有文档都是为了渲染网格图，而不是从中提取值。

任何帮助将不胜感激！

干杯，MB

PS第二个问题应该是：“有几个网格图块，如果一条路线跨越两个或多个图块我该怎么办？合并它们，引用两者......”

我想使用leaflet-shiny 和我自己的geojson 文件生成一个等值线图。可以从这里访问整个数据和 r 代码。我不明白为什么地图没有显示但没有错误消息。

我正在尝试使用我在此处下载的空间数据，以便在 ggplot2 中制作地图。

使用函数 fortify 时，我不断收到此错误：

有人知道这里出了什么问题吗？最欣赏！

编辑：

如评论中所述，可以在此处找到此问题的可能重复项。据说可以在安装包gpclib中找到问题的解决方案。我不确定是不是这样，这个包已经从 CRAN 中删除了。

我正在处理可以具有多个波段/通道的某些图像（参见多光谱和高光谱遥感图像）。

让我们假设我们事先不知道给定图像中存在多少条带。所有通道都被命名为“band1”、“band2”……等等。要阅读我们使用的图像

我们将其作为data.frame第一个：

可能存在特定“像素”的所有波段的内容为零的某些情况。因此，如果它是一个四波段图像，那么我们可以使用以下方法从其中删除这些行data.frame：

但是假设我的图像有 256 个波段。如何从数据集中删除此类行？

我想将一些 Landsat 波段/tiff 文件堆叠并写入 ENVI 格式的 BIP 交错。然而，结果总是以 BSQ 的形式出现，即使我将 bandorder 更改为 BIP。

下面是我的代码：

这也不起作用

任何帮助表示赞赏。

在 RI 中有一个带有重复点（坐标和属性）的 SpatialPointsDataFrame，我想删除所有具有相同数据的点...

我在sp包中找到了该remove.duplicates()功能，但它似乎只在位置上删除......还有其他方法吗？

谢谢你

E.

提前致谢。

我有 86,000 个纬度/经度点需要转换为 UTM。此处包含 25 个样本。我找到了许多帮助完整的 SO 帖子，但继续收到有关工作示例的错误，如下所述。

样本数据在这里：

有 3 个人的 ID、时间戳（还不是 POSIXct）、纬度和经度。

我遵循了此处提出的建议并使用了以下代码：

并得到

为了按照这里的建议删除负数，我跑了

并得到

最后，正如这里使用我运行的矩阵所建议的那样

并得到

任何关于如何将这些纬度/经度点转换为 UTM 的建议将不胜感激。纬度/经度在 WGS84 中，我的 UTM 区域是 12。