问题标签 [rgdal]

伙计们，

我需要一些专家意见。我一天中的大部分时间都在与这个问题作斗争，但我别无选择。

我正在尝试创建一张地图，描绘自上次火灾以来整个景观的时间

很明显有些地方是不会发生火灾的，这里用0表示，但是有点随意

我的代码如下所示：

数据可以在这里找到

如果你不想运行代码，这就是它在谷歌地球上的样子

这是主要问题，其中：

1）它几乎看起来像多边形边界信息被保留，因为棕褐色线是确切的边界，

2）这里并不总是报告“最大值”，尤其是在它出现的边缘附近？

3）我的目标是在物种分布建模中将此地图与其他栅格结合使用（因此使用栅格堆栈、模型和预测功能来制作地图）。当我尝试使用该图层制作地图时，由于 NA 的缘故，我最终得到了一堆洞......应该没有 NA 的地方。

如果我在没有这个时间的情况下运行模型，那么一切都很好。当我包含它时，生成的地图中的 NA 通常是您在第一张图片中看到棕褐色“多边形轮廓”的位置。

我怎样才能让这个光栅层有清晰的边缘，而不是像现在这样一团糟？？？

这就是我的最终地图物种地图的样子，而不是丢失的数据反映了火灾地图中的问题

我有一个要拼接在一起的栅格列表。除了 utm 区域外，投影是相同的。这是 UTM 区域 50 中一个 rasterStack 的坐标线：

这是另一个具有相同投影的 rasterStack，但它是 UTM 区域 51：

当我尝试在不先重新投影的情况下进行镶嵌时（就像在其他 GIS 中一样），我收到以下错误：

这让我相信我需要重新投影所有的 rasterStacks。这在计算上很耗时projectRaster，但我相信我可以这样做：

我是否正确假设我必须先重新投影？感谢您的任何指导。

我用 utm32-shapefile 制作了一个等值线图，现在，我想添加第二层，其中包含来自另一个数据帧的数据。第二个数据帧只有 Lon/Lat 信息。我想将 Lon/Lat 转换为 utm32 坐标，因此我可以将它们与 ggplot2 一起用作另一个收敛层。

数据框如下所示：

首先，我使用gsub将“，”重新编码为“.”。然后我将变量转换为数字。

它没有用。我只是得到一份我的纬度/经度变量的副本。

我没有找到有关如何执行此操作的良好文档。也许其他人可以帮助我？

编辑经过更多研究但仍然没有解决方案，我正在添加大量编辑以及指向 .shp 文件的链接。

形状文件包含在这里

我有一个包含 9 个多边形的 SpatialPolygonsDataFrame，每个多边形还包含多个嵌套多边形 - “洞”。数据摘要在这里。

数据样本structure如下。

如下例所示（九个之一），父多边形有多个孔。

这是我第一次涉足sp()、rgdal()、rgeos()和其他空间包，我发现了许多关于使用运算符提取区域等的有用帖子，但仍然存在问题。虽然这篇文章提供了接近的解决方案，但我似乎无法调整代码以满足此处描述的需求。

我想获得一个 SpatialPolygonsDataFrame，它只包含来自每组子列表（即 data@polygon）的父（最大）多边形。似乎我应该能够仅提取父多边形或“溶解”这些孔。

最终结果将是 9 个多边形，每个都是 9 个列表的父级，我可以将其导出为 ESRI shapefile。

任何建议，将不胜感激。

我正在尝试为空间多边形数据框中的多边形获取伪重心。今天我偶然发现了coordinates实际上为SpatialPolygonsDataFrame返回某些内容的函数。

不幸的是，我在SpatialPolygonsDataFramecoordinates的值的帮助中一无所获。有人能告诉我这些坐标是什么吗？

包中 SpatialPolygon 对象的 labpt 插槽sp在文档中模糊定义为显示多边形标签的点。我想知道它背后的计算以及它的属性是什么。通常我们确定它会在多边形中还是像质心一样它可以在多边形之外？

我在rgeos包中找到了一个计算类似名称对象的函数，但无法确定它是否真的是 sp 使用的函数。

作为物种分布建模项目的一部分，我正在尝试在 R 中创建背景数据点。我无法弄清楚如何将我的环境预测栅格图层（年降水量、夏季平均气温等）转换为 .grd 文件，以便我可以将它们读入 R。是否有一些基本步骤可以将栅格数据转换为 .grd文件？

我有兰伯特保角圆锥投影 x,y 信息。
我需要 WGS84 坐标。但我不知道lcc到底是什么。
我在下面提供了 lcc 信息。
有没有办法在 r 中将 lcc 转换为 WGS84？

lcc 信息：
第一个标准平行的纬度：30.0
第二个标准平行的纬度：60.0
原纬度：38.0 原经度：126.0
计算点东距：43
计算点北距：136
4 edge lon,lat 点：左上（43.3935， 123.3102), 左下(31.7944, 123.7613),
右上(43.2175, 132.7750), 右下(31.6518, 131.6423)

我尝试安装rgdal但configure.args没有运气。

安装的gdal版本：

找到 gdal 二进制工具：

找到 proj 头文件：

找到 proj lib 文件：

关于 gdal 的安装包：

编辑：我试过 install.packages('rgdal') 没有任何 configure.args。我得到了同样的错误。现在，我正在尝试使用已编译的 proj4 4.9.1 安装 rgdal，并将 LD_LIBRARY_PATH 导出到 proj4 4.9.1。不幸的是，它又失败了。这是我的新 script.sh：

一切都很好，但测试时失败：

更新：最后，我解决了这个问题。首先，通过 apt-get 删除 proj4-4.8。

其次，将我编译好的proj4-4.9 lib路径添加到'/etc/ld.so.conf.d/'路径下的一个新的conf文件（如rgdal.conf）中。

然后，使用 configure.args 安装 rgdal，引用 proj4-4.9 的编译 lib 路径。

现在，一切都很好。但是我无法使用 proj4-4.8 和默认的 configure.args 安装任何版本的 rgdal。我已经按照@Pascal 的建议清理了 proj4 的文件。所以我不确定这个问题的确切原因。

我的数据是以度、分和秒的形式进行的经纬度投影。这已经并以文本形式存储（character-class，忽略row.names）。

通过邮件列表中提供的答案，我已将数据转换为十进制形式...

之后，我将 转换data.frame为SpatialPoints对象。

然后继续通过以下方式转换投影spTransform

我已经使用此网站转换了示例中的第一个坐标（由于坐标指的是南和西，因此请舒尔提供负数），并提供以下坐标：633726.46 mE 8468757.51 mN Zone 18L，这是正确的。相比之下，sp变换后的对象是：

该spTransform手册特别提到应正确提供元数据。我的数据使用 WGS84 椭球，位于 18L 区域。更准确地说是epsg:32718。我哪里错了？