Raster 01: Plot Raster Data in R

最新推荐文章于 2022-12-30 17:07:19 发布

GeoSuper

最新推荐文章于 2022-12-30 17:07:19 发布

阅读量810

点赞数

分类专栏： R

R 专栏收录该内容

27 篇文章

订阅专栏

该博客介绍了如何在R中利用raster和rgdal包加载与绘制数字表面模型（DSM），并创建山体阴影图以生成优雅的基础地图。首先，加载并展示了HARV站点的DSM数据，然后通过hist函数展示DSM的高度分布。接着，将DSM分段并用不同颜色表示，最后将DSM与hillshade图叠加，以增强地理特征的视觉效果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#Know how to plot a single band raster in R.
#Know how to layer a raster dataset on top of a hillshade to create an elegant basemap

library(pacman)
p_load(raster, rgdal)

wd <- "G:/Rdata/neon_data/NEONDSAirborneRemoteSensing/"
setwd(wd)

DSM_HARV <- raster(paste0(wd,"HARV/DSM/HARV_dsmCrop.tif"))
#plot raster object
plot(DSM_HARV,
     main = "Digital Surface Model\n NEON Harvard Forest Field Site")

#plot distribution of raster values
DSMhist <- hist(DSM_HARV,
                breaks = 3,
                main="Histogram Digital Surface Model\n NEON Harvard Forest Field Site",
                col = "wheat3",
                xlab = "Elevation(m)")

#where are breaks and how many pixels in each category
DSMhist$breaks
# [1] 300 350 400 450
DSMhist$counts
# [1] 32180 67382   438

#Plot using breaks
plot(DSM_HARV,
     breaks = c(300,350,400,450),
     col = terrain.colors(3),
     main = "Digital Surface Model(DSM)\n NEON Harvrad Forest Field Site",
     axes = F#不显示坐标轴 
     )

DSM_hill_HARV <- raster(paste0(wd,"HARV/DSM/HARV_DSMhill.tif"))
plot(DSM_hill_HARV,
     col=grey(1:100/100),
     legend=F,
     main = "DSM with Hillshade\nNEON Harvard Forest Field Site",
     axes = F)
plot(DSM_HARV,
     col = rainbow(100),
     alpha = 0.4,
     add = T,
     legend = F,
     )