参考:https://stackoverflow.com/questions/33227182/how-to-set-use-ggplot2-to-map-a-raster
levelplot方法:https://www.rdocumentation.org/packages/rasterVis/versions/0.46/topics/levelplot-methods
利用rasterVis::levelplot:绘制地图raster
# 1. load packages
library(rgdal)
library(RColorBrewer)
library(rasterVis)
# 2. read data
oregon = readOGR("D:/数据/行政区划/中国底图数据/省级行政区.shp")
r = raster("F:/大论文/PM25反演/应用数据/追算tif/年/1999MEAN追算.tif")
修改颜色、查看属性信息(数据类型、空间范围、地理坐标系统)
# define a color rmap palette
colr = colorRampPalette(brewer.pal(7,"RdYlBu"))
colr = colorRampPalette(c("blue","red"))
display.brewer.pal(7,"RdYlBu")
display.brewer.all(type = "seq")
display.brewer.all(type = "div")
display.brewer.all(ty