【问题标题】：Is there a raster function to reverse the x and y axis?是否有反转 x 和 y 轴的栅格函数？
【发布时间】：2020-01-20 22:26:14
【问题描述】：

我有一个从 iwd 创建的栅格数据集。我绘制了一个填充轮廓图，但我想反转 x 轴和 y 轴，以便数字减少，消除空白并垂直夸大 y 轴。像在 ggplot 或 plot 中那样设置 xlim 和 ylim 不起作用。

如果无法反转栅格数据集的 x 轴和 y 轴，我如何在转换为 s3 后保持 s4 类数据集的分辨率？例如，如果我使用filled.contour 而不是filledContour。

这是我的代码和情节，这是非常基本的，因为我尝试过的没有产生任何结果：

idw.out <- gstat::idw(Z ~ 1, core2, grd, idp = 1.5)
r <- raster(idw.out[1])
r.contour <- filledContour(r)
r.contour

我正在寻找的比例示例如下：

干杯

【问题讨论】：

core2 是什么？
数据是从洪泛区沉积岩心收集的。不过，我确实解决了我的问题并在下面发布了答案。

标签： r ggplot2 raster

【解决方案1】：

我用下面列出的代码解决了我的问题。因为我有 S4 类数据，所以我必须手动设置轴范围以与 xlim 和 ylim 重合。 'filled.contour' 绘图不需要此步骤，但适用于 'filledContour' 绘图。

来自 core2 的插值粒度数据

idw.out <- gstat::idw(Z ~ 1, core2, grd, idp = 1.5)

转换为栅格

r=raster(idw.out[1], layer = 1, values=TRUE)

添加配色方案和绘图

b = c(0,5,10,15,20,30,40,50)
col = rev(bpy.colors(length(b)-1))
r.contour = filledContour(r, zlim=c(0,50), xlim=c(9,-4), ylim=c(5.3,0),
                          asp = NA, xaxs = "i", yaxs = "i", las = 1,
                          col=col, levels=b,
                          xlab="grain-size (phi)", ylab="core depth (m)", main="Core 2")

【讨论】：