如何旋转图像R栅格

Question

我下面的代码会将图像保存到我的电脑上。我希望将图像围绕其中心（或左下角）旋转45,90和135度，然后保存为3个不同的图像。我怎么能这样做？

library(raster)
r1 <- brick(system.file("external/rlogo.grd", package="raster"))
x <- crop(r1, extent(0,50,0,50))
plotRGB(x)
png(width=50, height=50)
par(mai=c(0,0,0,0))
image(x)
dev.off()

--------- UPDATE1 -------------------------

根据接受的答案，工作代码如下

library(raster)
r1 <- brick(system.file("external/rlogo.grd", package="raster"))
r1
x <- crop(r1, extent(0,ncol(r1),0,nrow(r1)))
plotRGB(x)

x1 <- 0:ncol(x)
y1 <- 0:nrow(x)
z <- matrix(1, nrow=length(x1), ncol=length(y1))

col.mat <- t(apply(matrix(rgb(getValues(x)/255), nrow=nrow(x), byrow=TRUE), 2, rev))

# Rotate 45 degrees
persp(x1, y1, z, zlim=c(0,2), theta = 20, phi = 90, 
      col = col.mat, scale=FALSE, border=NA, box=FALSE)
png("SaveThisPlot.png")
persp(x1, y1, z, zlim=c(0,2), theta = 20, phi = 90, 
      col = col.mat, scale=FALSE, border=NA, box=FALSE)
dev.off()

Answer 1

对于90度旋转，这是一个简单的解决方案：

image(t(flip(x, 1)))
image(t(flip(x, 2)))
plotRGB(t(flip(x, 1)))
plotRGB(t(flip(x, 2)))

对于45度和135度旋转，它会有点棘手。可能还有其他方法，但我会使用persp函数，并为theta参数赋予不同的角度。

只需正确设置对persp的调用即可。 x1，y1和z只是persp函数的输入（有关该函数参数的更多信息，请参阅?persp）。 col.mat是一个包含颜色值的矩阵。

x1 <- 0:ncol(x)
y1 <- 0:nrow(x)
z <- matrix(1, nrow=length(x1), ncol=length(y1))
col.mat <- t(apply(matrix(rgb(getValues(x)/255), nrow=nrow(x), byrow=TRUE), 2, rev))
# the transposing and reversing are just to get the colors in the same 
# spots as they're in when viewing plotRGB(x).
#
# getValues(x) is how you get the rgb colors, in the form of a 3-column matrix.
# rgb(getValues(x)/255) converts them into hex code, which is convenient enough here.

如果您发现这是您要查找的镜像，请尝试以不同方式填充颜色矩阵。例如：

col.mat <- matrix(rgb(getValues(x)/255), nrow=nrow(x))

正如您所知，正确填充颜色矩阵是使此方法适合您的关键。我将把它作为练习留给读者，以弄清楚如何对颜色矩阵进行任何其他操作。

现在，请致电persp。在这里，我设置zlim值，因此有一个范围包括1.因为我将所有z值设为1，您需要设置一个有效范围，否则persp将抛出一个关于无效限制的错误。它不喜欢从1到1的范围。

# Rotate 45 degrees
persp(x1, y1, z, zlim=c(0,2), theta = 45, phi = 90, 
    col = col.mat, scale=FALSE, border=NA, box=FALSE)

这是135度：

# Rotate 135 degrees
persp(x1, y1, z, zlim=c(0,2), theta = 135, phi = 90, 
    col = col.mat, scale=FALSE, border=NA, box=FALSE)

保存图表的方式与您在问题中显示的方式相同：

png("SaveThisPlot.png")
persp(x1, y1, z, zlim=c(0,2), theta = 135, phi = 90, 
    col = col.mat, scale=FALSE, border=NA, box=FALSE)
dev.off()

Answer 2

在不抛出错误或警告消息的情况下执行此操作（阅读帮助页面并播放示例后）：

extent(x) <- extent(0, 360, -90, 90)
rr <- rotate(x)
 png()
 plotRGB(x)
dev.off()

但你可能不喜欢它，因为它假设你可能没有特定的坐标系。如果要拉出数据元素，请使用@运算符：

str(x)  The slots are named 'data' and 'values'
image(x)
image(matrix(x@data@values[,1], 50, byrow=TRUE)) # rotate one layer