Je viens vers vous avec une question de cartographie sous R:
J'ai réalisé un krigeage pour le fun à partir de donnée imaginaire composée de latitude, longitude et d'un comptage marketing sur Paris.
Je voudrais réaliser une carte de paris avec ma prédiction:
J'ai pour cela téléchargé la carte de paris sur https://opendata.paris.fr/page/home/ (très cool d'ailleurs).
Voila mon problème, mon code marche mais je me retrouve avec une carte de prédiction + au dessus une carte de paris en blanc. Je ne vois donc pas les prédictions de paris...
Les prédictions sont cachées par la carte blanche de paris. Pourtant en changeant le fill = white, je n'ai pas trouvé de solution.
Est-ce que quelqu'un a une idée de solution? Ou une idée?
Code : Tout sélectionner
###Carte Paris
list.files(pattern="^quartier")
fr<-readOGR("quartier_paris.shp",layer="quartier_paris")
line <- fortify(fr, quartier="name")
###Transformer les sorties du krigeage
krig.output=as.data.frame(lzn.kriged)
names(krig.output)[1:3]<-c("Longitude","Latitude","var1.pred")
###la carte
plot<-ggplot(data=krig.output,aes(x=Longitude,y=Latitude))#start with the base-plot and add the Kriged data to it
layer1<-c(geom_tile(data=krig.output,aes(fill=var1.pred)))#then create a tile layer and fill with predicted
layer2<-c(geom_path(data=line,aes(long, lat, group=group),colour = "grey40", size=1))#then create an outline
plot+layer1+layer2+scale_fill_gradient(low="#FEEBE2", high="#7A0177")+coord_equal()
ggplot(krig.output) + geom_tile(aes(Latitude, Longitude, fill = var1.pred)) + geom_polygon(data = line,aes(long, lat, group = group), fill = "white",colour = "grey20", size = 0.1) +coord_quickmap(xlim = range(krig.output$Latitude),ylim = range(krig.output$Longitude)) + xlab("") + ylab("")
Nb: j'ai utilisé
library(sp)
library(rgdal)
library(raster)
library(dplyr)
library(mapview)
library(mapedit)
library(sf)
library(readr)
library(ggplot2)
library(tmap)
Merci pour l'aide.
Bonne journée,
Cordialement
Tibo