问题标签 [rastervis]

我有一个印度尼西亚水域盐度的 NetCDF 文件，有 4 个维度（经度、纬度、深度和时间）。如何在此处从我的数据下载数据创建每周复合： https://onedrive.live.com/redir?resid=6FFDD661570C7D0A% 21177 此处输出地图：https : //onedrive.live.com/redir?resid=6FFDD661570C7D0A% 21176

我想将栅格转换为矢量并使用 apply 来获得平均值，但我无法使用绘制矢量数据rasterVis

我有一个印度尼西亚水域盐度的 NetCDF 文件，有 4 个维度（经度、纬度、深度和时间），在此处下载数据：https ://onedrive.live.com/redir?resid=6FFDD661570C7D0A%21177

我想使用levelplot（rastesVis包）绘制具有相同范围（最小值 - 最大值）值的盐度数据。不幸的是，我为每个情节获得了不同的价值范围。at=seq(20, 35, 2)没有在我的情节中定义相同的范围？

我想将颜色与栅格类匹配。假设我有以下颜色和这些栅格：

我想对“s”中的栅格进行分类，以便从 0 到 0.1 的值以“深红色”着色，从 0.9 到 1 的值以“深蓝色”着色。基本上，我想将我的数据分类为“n”类，然后将每个类分配给我选择的颜色。

最后，colorkey 应该将各种栅格类作为标签。在这里可以找到一些见解1和2

我想删除 R 在格点图周围的默认边距。这意味着我想摆脱红色矩形之外的所有空白。这是示例：

我有兴趣制作两个并排显示在同一个 Xquartz 设备窗口中的 3d 地形图。使用 rgl 包在同一个窗口中显示两个 3d 图很简单——文档中有很多示例 usingmfrow3d()和其他方法。

但是，我需要使用plot3D()rasterVis 包中的函数。我使用plot3D()rasterVis 而不是 rgl的原因surface3d()是我需要使用drape参数 inplot3D()将栅格中的值显示为 3d 地形图上的颜色（并且该栅格的值与创建该栅格的值不同）图上的 z 轴）。如果有人对 rgl 函数中类似于悬垂的东西有提示，我也会对此感兴趣！

当我尝试mfrow3d()与该plot3D()功能一起使用时，它会显示一系列空白设备窗口，而不是在同一个窗口中并排显示两个图。

以下是 rasterVis 文档中使用 plot3D 制作地形图的一些代码：

在这里，我尝试使用 mfrow3d 并排绘制 2 个相同的火山图，一个蓝色，一个红色，我改编自 rgl 文档：

我正在尝试使用 rasterVis 的 plot3D 做些什么？

我有以下问题，我未能重命名层堆栈。

这是我的代码示例。

如果只使用数字

结果：名称（oo）

如图所示：

但我需要没有“X”，即：

显然不适用于数字，因为如果我尝试使用：

结果在剧情中没问题。

PD：也尝试使用参数：

欢迎所有建议。

问候！

我们可以在rasterVis::levelplot条形标签中包含下标/上标字符吗？

考虑以下RasterStack, s:

对象的默认绘图方法Raster*允许将表达式传递给main参数：

With rasterVis::levelplot，strip 名称是通过names.attr参数传递的，但似乎这些character在最终被传递给之前被强制传递给lattice::levelplotwith strip.custom(factor.levels = names.attr)。

结果是：

除了修改源代码之外，有没有办法在带状标签中使用表达式（或以其他方式实现下标/上标字符）rasterVis::levelplot？

我想：

1. 使用以下方法将栅格重新分类ras为九个类reclassify

2. 在每种颜色旁边提供一个带有值的颜色键（请参见下面的示例图）。colorkey 不应像示例图中那样拆分，而是显示为具有 9 种颜色的单一色带。

3. 使用levelplot函数显示数据（我可以这样做）

谢谢你的帮助。AEZ

请使用这些颜色：

下面是dput我的栅格图层：

我正在使用 RasterVis 和 levelplot 制作一些栅格的格状图。我目前对大多数事情都很好，但我想将每个面板的标题从文件名更改为选定的字符串（文件名很复杂而且很长，我只想使用一年，例如“2004”）。

查看levelplot 页面，这表明 levelplot 会根据参数“useRaster”寻找一些设置，它要么转到panel.levelplot 要么 panel.levelplot.raster，但我正在努力使用后面的这些功能。

非常感谢任何帮助，这里有一些示例代码；

在上面例如，我希望“layer.1.1”为“2004”，依此类推到 2007

我想使用 rasterVis 中的 levelplot 函数在栅格时间序列对象 (rts) 中显示结果。

这是来自 rts 包的简短代码：

为了从 .rts 强制转换为我使用的光栅：

r从重新投影UTM coordinates to latlong coordinates

如何读取所有 4 层myfile100.tif？ ras = raster("myfile100.tif")只读layer 1

我收到此错误：

我怎样才能解决这个问题？也许还有一种更短更合理的方法来做到这一点。感谢您的建议。