freq高程频率统计.jpg
# Thu Dec 24 14:00:16 2020 -
# 字符编码：UTF-8
# R 版本：R x64 4.0.3 for window 10
# cgh163email@163.com
# 个人笔记不负责任
# —— 拎了个梨🍐
.rs.restartR()
require(terra)
# freq频率表,高程数量统计。
rm(list=ls());gc()
# SpatRaster值的频率表。 不计算资产净值。
r <- rast(nrow=10, ncol=10)
set.seed(2)
values(r) <- sample(5, ncell(r), replace=TRUE)
plot(r)
freq(r, FALSE)
freq(r)
x <- c(r, r/3)
freq(x, FALSE)
freq(x)
freq(x, digits=1)
freq(x, digits=-1)
freq(x, value=5)
# Thu Dec 24 14:03:04 2020 --
#
rs.img <- rast('https://s1.ax1x.com/2020/10/27/BQCc2q.jpg')
tmp <- terra::freq(rs.img$BQC_1)
plot(rs.img$BQC_1,main='before')
plot(tmp[,3],type='l',main='freq高程频率统计')
dev.copy(png,'freq高程频率统计.jpg');dev.off()
# Thu Dec 24 14:10:42 2020 --end
