使用ggplot和gmap绘制栅格图像的问题

问题描述 投票:1回答:1

我想用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。

我对这个软件包不是很熟悉,一直用的都是经典的绘图功能。谁能帮帮我吗?

r ggplot2 ggmap r-raster
1个回答
0
投票

尝试添加 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)

enter image description here

© www.soinside.com 2019 - 2024. All rights reserved.