Как улучшить пространственную растровую карту с помощью ggplot по сравнению с spplot?

Как я могу улучшить легенду построения пространственной растровой карты с помощью ggplot по сравнению с легендой spplot ()?

Я хотел бы построить пространственные карты, используя ggplot () вместо ssplot (), однако есть несколько вещей, которые я хотел бы улучшить по сравнению со spplot:

  1. создать легенду ggplot, которая идет от маленького (внизу) к большому values ​​(вверху)
  2. Сделайте разрывы в легенде ggplot, аналогичные легенде ssplot (), чтобы я знал, каковы границы каждого цвета.

## load packages
require(raster)
require(ggplot2)
require(rgdal)
require(RColorBrewer)
set.seed(1)

r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40,
          crs="+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100
+ellps=WGS84")
r <- setValues(r,matrix(rnorm(1600, mean=0.4,sd=0.2))) 

## 1. spatial map with spplot
cuts <-seq(minValue(r),maxValue(r),length.out=8)
cuts = round(cuts,digits=2)
col.regions = brewer.pal(length(cuts)+3-1, "RdYlGn")
print( 
spplot(as(r, 'SpatialGridDataFrame'),at=cuts,
col.regions=col.regions,
colorkey=list(labels=list(at=cuts),at=cuts), pretty=TRUE,
scales=list(draw=T)
) 
)

## 2. spatial map with ggplot
p = rasterToPoints(r); df = data.frame(p)
colnames(df) = c("x", "y", "NDVI")

p  <- ggplot(data=df) + geom_tile(aes(x, y, fill=NDVI)) +
coord_equal() + labs(x=NULL, y=NULL) + 
scale_fill_gradient2(low="red", mid="yellow",high="green",
limits=c(minValue(r),maxValue(r)), midpoint = 0.4) + theme_bw() +
scale_x_continuous(expand=c(0,0)) + scale_y_continuous(expand=c(0,0))
print(p)

результат ssplot () ssplot

результат ggplot () ggplot

8
задан Janvb 29 July 2011 в 07:28
поделиться