我想用ggplot2实现光栅二进制地图的绘制,也想用ggmap添加一个石蕊地图作为背景。根据不同的帖子,我实现了这些代码行:'
library(sp)
library(raster)
library(rasterVis)
library(ggplot2)
library(ggmap)
library(tmaptools)
library(rgdal)
test <- raster('C:/France_accuracy/img_fig/image_france_2009.tif')
test_df <- rasterToPoints(test)
test_df <- data.frame(test_df)
colnames(test_df) <- c("X","Y","Values")
head(test_df)
#Simple plot using ggplot
p1 <- ggplot() + geom_raster(data=france_df, mapping=aes(X, Y, fill= factor(Values)))
#Background and final map
background <- get_stamenmap(bbox = c(-0.7744, 44.2001,-0.5286, 44.3017),
maptype='toner-background', zoom = 5)
finalmap <- ggmap(background) +
geom_raster(data=france_df, mapping=aes(X, Y, fill= factor(Values)))
我得到了这个错误信息
"Error: geom_raster only works with Cartesian coordinates Run
rlang::last_error()
以查看错误发生的位置。此外。 警告信息。删除了169728行,其中包含缺失的值(geom_raster)""
这里是小样图的链接 https:/drive.google.comdriveu0folders1945TBCzW9lmKjOaN_4CwoD_XfT9PXoO5。
我对这个软件包不是很熟悉,一直用的都是经典的绘图功能。谁能帮帮我吗?
尝试添加 coord_cartesian()
到 finalmap
:
rngX <- range(test_df$X)
rngY <- range(test_df$Y)
background <- get_stamenmap(bbox = c(rngX[1], rngY[1], rngX[2], rngY[2]),
maptype='toner-background', zoom=13)
finalmap <- ggmap(background) +
geom_raster(data=test_df, mapping=aes(X, Y, fill= factor(Values)), alpha=0.3) +
coord_cartesian()
print(finalmap)