我正在尝试使用栅格注释功能在ggplot地图上叠加图像。我已按照此帖子How do I overlay an image on to a ggplot?中的说明操作,但我使用“ggplot”代替“qplot”来绘制我的数字。
然而,图像要么没有出现在图上(当我运行绘图函数的前三行时:ggplot,geom_polygon和annotation_raster),或者如果我运行整个序列,则会产生以下错误:
Error in if (nrow(layer_data) == 0) return() : argument is of length zero
我已经尝试了各种各样的事情,所以如果有人有任何想法如何解决问题,我将非常感激!
### Reproducible example
library(png)
library(ggplot2)
#Create data
scores <-c(5,10,-12,20)
area <-c(2630, 3970, 1550, 7485)
lat <- c(-16, -23, -30, -27)
long <- c(132, 143, 120, 140)
data <- as.data.frame(scores)
data[,2] <- as.data.frame(area)
data[,3] <- as.data.frame(long)
data[,4] <- as.data.frame(lat)
colnames(data)<-c("scores","area","lon","lat")
#Load PGN image for raster annotation
mypngfile <- download.file('http://upload.wikimedia.org/wikipedia/commons/thumb/c/c1/Rlogo.png/200px-Rlogo.png', destfile = 'mypng.png', mode = 'wb')
mypng <- readPNG('mypng.png')
#Breaks and colours for mapping
brks <- c(-20,-10,0,10,20)
colors <- rev(c(rgb(0,0,0.6),rgb(0,0.6,1),rgb(0,0.8,1),rgb(0.8,1,0.2)))
plot_data <- as.numeric(cut(data$scores,breaks=brks))
col_data <- as.factor(plot_data)
#Map of Australia
Australia<-map_data("world",region = c("Australia", "Australia:Tasmania"))
Oz<-subset(Australia,long>110 & lat>-50 & long<155)
#Create plot
ggplot() +
geom_polygon( data=Oz, aes(x=long, y=lat, group=group), colour="black", fill="grey68") +
annotation_raster(mypng,xmin=120, xmax=130, ymin=-45, ymax=-49) +
geom_point(data=data,aes(x=long,y=lat,color=col_data, size=area),show_guide=F) +
scale_size_area(name="area",max_size=20)
答案 0 :(得分:1)
只需要进行小的修正(我也调整了R徽标的坐标以使其可见):
ggplot(data=Oz, aes(x=long, y=lat)) +
geom_polygon(aes(group=group), colour="black", fill="grey68") +
annotation_raster(mypng, xmin=120, xmax=130, ymin=-37, ymax=-44) +
geom_point(data=data, aes(colour=col_data, size=area), show_guide=F) +
scale_size_area(name="area", max_size=20)