如何将 10 张地图绘制为 R 中的一个图形？

编程问答 2022-05-31

问题描述

我正在使用此代码将 10 张地图绘制为 r 中的一个图形？我想要地图下方的 alabeltical 订单 (a) (b) 等，以及图上方的 2009-2019 年。 我在 R 中给出了绿色作为 ndvi 的结果颜色。它并不突出。我可以使用哪种颜色使其更显眼？

library(raster)
library(rgdal)
library(rasterVis)
library(lattice)
library(RColorBrewer)

r1 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=1)

r2 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=2)
r3 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=3)
r4 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=4)
r5 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=5)
r6 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=6)
r7 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=7)
r8 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=8)
r9 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=9)
r10 <- raster("E:/GB extent(2009-2019)/rasterstack.tif",band=10)

boudary <- readOGR(dsn= "E:/Shapefile/New folder/GB_Bondary.shp",layer =  "GB_Bondary")





raster_stack <- stack(r1,r2,r3,r4,r5,r6,r7,r8,r9,r10)

#levelplot(raster_stack,layout=c(5,2)


colr <- colorRampPalette((brewer.pal(3,'Greens')))

levelplot(raster_stack,2),margin=FALSE,#main=list("title",col="Black",fontface='bold'),colorkey=list(
            space='right',labels=list(at=c(0.5,0.6,0.7,0.8,0.9,1.0),font=5,axis.line=list(col='black'),width=0.75
          ),par.settings=list(
            strip.border=list(col='transparent'),strip.background=list(col='transparent'),axis.line=list(col='transparent')
          ),scales=list(draw=FALSE),col.regions=colr,at=seq(0.5,1.0,len=101),names.attr=(c('2009','2010','2011','2013','2014','2015','2016','2017','2018','2019')),font=list(face="bold"),cex=.4)+
  layer(sp.polygons(boudary,lwd=0.5))

###########################################################################

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

mapsplotrrr-rasterstack