问题标签 [r-raster]

我正在拼命尝试，writeRaster()但是由于栅格很大，我需要大量临时内存。我的C:/驱动器上的空间有限，因此我想将临时目录更改为D:/TEMP/. 我尝试了在互联网上找到的不同方法，例如：

更改临时目录

http://r.789695.n4.nabble.com/How-do-I-set-the-Windows-temporary-directory-in-R-td876483.html

等等

rasterOptions(tmpdir = "D:/RTEMP/")没有解决问题。即使它在计算栅格时起作用，也不会影响writeRaster()功能。

如果有人可以帮助我，我将非常感激。

这是我的sessionInfo()：

我在等面积 Behrmann 投影中有一个栅格，我想将它投影到 Mollweide 投影和绘图。

但是，当我使用以下代码执行此操作时，绘图似乎不正确，因为地图延伸到两侧，并且有各种陆地的轮廓，我不希望它们出现。此外，地图超出了情节窗户。

谁能帮我把这个画得很好？

谢谢！

使用的数据文件可以从此链接下载。

这是我到目前为止的代码：

在将栅格的值提取到点时，我发现我有几个NA's，而不是使用函数的buffer和fun参数extract，而是我想将最近的非NAPixel 提取到重叠的点NA。

我正在使用基本的提取功能：

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

下面是我的代码：

这也不起作用

任何帮助表示赞赏。

我想知道你们中是否有人可以帮助我完成以下处理 R raster 包中的 focus() 函数的任务。

默认情况下，focus() 函数将遍历给定栅格（以下称为“背景”栅格）的每个像元，并将给定函数应用于用户定义的移动窗口定义的相邻值。为了优化和加快我在大型栅格/研究区域上的计算，我想仅在“背景”栅格在“移动”覆盖的范围内具有某些值（例如大于零）时应用此功能（过滤器）窗口”并跳过所有其他焦点单元格。这样，过滤器就不会花时间计算任何不需要的焦点值。

下面是一个可重现的小示例和内嵌注释：

我应该如何更改此功能以获得预期的结果？

我有五个 netcdf 文件，其中每个文件都包含一个时间段的数据。我想单独计算每个单元格的整个时间跨度的第 98 个百分位数。netcdf 文件的累积文件大小约为 250 MB。

我的做法是：

但是，我收到此错误：

我的代码的堆叠部分有效，崩溃发生在 calc 函数期间。你知道这可能来自哪里吗？这可能是数据存储位置（内存/磁盘）的问题吗？

我正在使用一些环境数据（climond）对 35 个栅格进行 PCA。一切正常，我在 R 中使用此命令在光栅堆栈上执行 PCA：

所有的栅格看起来都很好，但是每个分量都解释了完全相同的方差比例（0.029），所以最后它们总和为 1。这对我来说有点奇怪，因为我已经习惯了例如前三个 pca 的结果轴解释了例如 50% 的方差，而其余的组件解释的越来越少。那么我的结果是正确的还是应该对princcomp进行一些修改？这就是它的样子：

然而，特征值（或此处的标准差）正在减少：

我希望这是一个相关的问题！提前感谢您的回答。

我正在使用如下行所示的栅格函数。我的最后一行产生了一些输出。该输出有一行写着dimensions : 240, 320, 76800 (nrow, ncol, ncell)。我想重新打印该图像，但只说前 200 行和前 300 列。我怎样才能做到这一点？下面的倒数第二行绘制了整个图像

===============================update1

我确实png(filename = '~/x.png');par(mar=rep(0, 4), xpd = TRUE, oma=rep(0, 4),bty='n') ; plot(r,xlim=c(0,200),ylim=c(0,200),legend=FALSE,axes=FALSE); dev.off()保存了裁剪后的图像。我能够摆脱图例和轴以及黑匣子。但问题是保存的图像包含的不仅仅是裁剪部分 - 例如图像周围的白色部分。我只想保存原始的裁剪部分（保持图像大小 200*200 像素）。请让我知道该怎么做？

此外，如何将与上述裁剪部分相对应的红色正方形添加到原始图像中？我的意思是我想在原始图像的顶部获得一个红色正方形（只有边缘），然后将这个（原始图像+正方形）保存为新图像。

我怎么能这样做？

更新2++++++++++++++++++++++++++++++++++++++++++++++++++++

添加可重复的示例以显示我所说的白色背景

下面的最后一行绘制裁剪图像。我希望该图像为 100*100，因为我的 xlim 和 ylim 为 100。但我看到一个白色背景，如下例所示。（你看不到背景。但是如果你在你的机器上运行代码并打开图像，你会看到它）

我有一系列二进制栅格图层（ascii 文件），显示欧洲和非洲某个物种的存在/不存在。该文件基于未投影的纬度/经度 (WGS84) 数据。我的目标是使用 R 计算存在区域（我无权访问 ArcGIS）。

我知道 raster 包有计算面积的功能，但我担心这对于未投影的数据不准确。我也看过 raster 包中的 cellStats 函数，可以用它来“求和”占用的单元格数，但我觉得这有同样的问题。

有谁知道准确计算存在区域的方法，考虑地球曲率？

谢谢！

自 R 3.xx 发布以来，我在绘图时遇到了问题。到目前为止，我总是能够回到旧版本来避免这个问题，但现在我使用 R 2.15.x 不可用的包。

小例子：

我在用着：

我已经能够找到有关此错误的可能原因的一些信息（R 3.0.0 已发布！（有什么新功能，以及如何升级）），但我从来没有找到任何提示如何解决这个问题。

提前感谢您的帮助！