R语言terra包,hist直方图。
# Fri Dec 25 13:58:30 2020 -
# 字符编码:UTF-8
# R 版本:R x64 4.0.3 for window 10
# cgh163email@163.com
# 个人笔记不负责任
# —— 拎了个梨🍐
.rs.restartR()
require(terra)
# hist直方图。
# 创建SpatRaster值的直方图。对于大型数据集,使用maxcell示例。
rm(list=ls());gc()
r1 <- r2 <- rast(nrows=50, ncols=50)
values(r1) <- runif(ncell(r1))
values(r2) <- runif(ncell(r1))
rs <- r1 + r2
rp <- r1 * r2
opar <- par(no.readonly =TRUE)
par(mfrow=c(2,2))
plot(rs, main='sum')
plot(rp, main='product')
hist(rs)
dev.copy(png,'hist直方图.jpg');dev.off()
a <- hist(rp)
a
x <- c(rs, rp, sqrt(rs))
hist(x)
par(opar)
dev.copy(png,'hist直方图,三个图层.jpg');dev.off()
# Fri Dec 25 14:15:16 2020 --end