2

我正在尝试复制我在raster使用 new(ish) terraR 包时所做的情节,但alphain 中的论点terra::plot似乎与在raster.

我正在尝试使用 alpha 层来指示基础(彩色)层中的不确定性。下面是一些表明我的问题的简化代码。使用该raster包,alpha 层会根据每个单独像素的值更改透明度。使用时terra似乎不起作用。

我想我正在使用最新版本的terra(1.3.4),并且我在 Mac 上(版本 15.5.11)。

library(terra)
library(raster)

par(mfrow=c(1,2))

plot_cols<-terrain.colors(n=5)

## raster package
raster.template<-raster(xmn=1, xmx=5, ymn=1, ymx=5, nrows=5, ncols=5)

legend.ras.comm<-raster.template
legend.ras.comm[]<-rep(c(5,4,3,2,1), 5)

legend.ras.alpha<-raster.template
legend.ras.alpha[]<-rep(c(0.9,0.7,0.45,0.2,0), 5)

plot(t(legend.ras.comm), col=plot_cols, legend=FALSE, axes=FALSE, box=FALSE, main="Raster package")
plot(t(legend.ras.comm), col="#000000", alpha=legend.ras.alpha, add=TRUE, legend=FALSE)


## terra package
raster.template<-rast(xmin=1, xmax=5, ymin=1, ymax=5, nrows=5, ncols=5)

legend.ras.comm<-raster.template
legend.ras.comm[]<-rep(c(5,4,3,2,1), 5)

legend.ras.alpha<-raster.template
legend.ras.alpha[]<-rep(c(0.9,0.7,0.45,0.2,0), 5)

plot(t(legend.ras.comm), col=plot_cols, legend=FALSE, axes=FALSE, main="Terra package")
plot(t(legend.ras.comm), col="#000000", alpha=legend.ras.alpha, add=TRUE, legend=FALSE, axes=FALSE)

绘图输出

对此的任何帮助将不胜感激。

4

1 回答 1

3

您现在可以使用terra 1.3-22

我得到,使用你的代码:

## raster package
library(raster)
legend.ras.comm <- raster(xmn=1, xmx=5, ymn=1, ymx=5, nrows=5, ncols=5, vals=rep(c(5,4,3,2,1), 5))
legend.ras.alpha <- setValues(raster.template, rep(c(0.9,0.7,0.45,0.2,0), 5))
legend.ras.comm <- t(legend.ras.comm)
## terra package
library(terra)
legend.ter.comm <- rast(legend.ras.comm)
legend.ter.alpha <- rast(legend.ras.alpha)

par(mfrow=c(1,2))
plot_cols <- terrain.colors(n=5)

plot(legend.ras.comm, col=plot_cols, legend=FALSE, axes=FALSE, box=FALSE, main="Raster package")
plot(legend.ras.comm, col="#000000", alpha=legend.ras.alpha, add=TRUE, legend=FALSE)

plot(legend.ter.comm, col=plot_cols, legend=FALSE, axes=FALSE, main="Terra package")
plot(legend.ter.comm, col="#000000", alpha=legend.ter.alpha, add=TRUE, legend=FALSE, axes=FALSE)

在此处输入图像描述

您还可以使用terra影响所有颜色的单个值

plot(legend.ter.comm, col=plot_cols, alpha=.75)

或者通过颜色改变透明度。像这样的东西:

plot_cols_alpha <- terrain.colors(n=5, alpha=seq(.75, 1, .05))
plot(legend.ter.comm, col=plot_cols_alpha)
于 2021-08-19T14:50:57.127 回答