# Fri Sep 04 09:45:09 2020 -
# 字符编码：UTF-8
# R 版本：R x64 4.0.2
# cgh163email@163.com
#个人笔记不负责任
# —— 拎了个梨🍐
rm(list=ls());gc()
.rs.restartR()
library(raster)
# Fri Sep 04 10:14:30 2020 --addLayer添加或放置图层----------------------------
file <- system.file("external/test.grd", package="raster")
s <- stack(file, file, file) # 添加2个图层
plot(s)
r <- raster(file) # 读取栅格数据
plot(r)
s <- addLayer(s, r/2, r*2) # 添加两个图层。同时叠加计算高程
plot(s)
nlayers(s) # 图层数量
s <- dropLayer(s, c(3, 5))
nlayers(s) # 返回图层数量
plot(s)
# Fri Sep 04 20:27:46 2020 --DIY
file.name <- '惠东水底山环线徒步20米精度高程.grd' # 新科图新等导出grd高程
map.gaocheng <- stack(file.name,file.name)
plot(map.gaocheng)
map.gaocheng1 <- raster(file.name)
plot(map.gaocheng1)
# 对高程做简单的加减乘除计算：
map.gaocheng2 <- addLayer(map.gaocheng,map.gaocheng1/2,map.gaocheng1*2)
plot(map.gaocheng2) # 并列显示
stack函数
对高程做简单的加减乘除计算addLayer函数展示
