Comment tracer une carte de chaleur sur une carte spatiale

Je suis nouveau à l'analyse spatiale des données dans R et voudrait faire quelque chose de facile, je suis encore avoir des difficultés...
J'ai une grande table avec latitudes et longitudes

sample = structure(list(Longitude = c(-0.19117, -0.211708, -0.206458, 
-0.173862, -0.156618), Latitude = c(51.489096, 51.520075, 51.525301, 
51.482442, 51.495752), Location_Easting_OSGR = c(525680L, 524170L, 
524520L, 526900L, 528060L), Location_Northing_OSGR = c(178240L, 
181650L, 182240L, 177530L, 179040L)), .Names = c("Longitude", 
"Latitude", "Location_Easting_OSGR", "Location_Northing_OSGR"
), row.names = c(NA, -5L), class = c("data.table", "data.frame"
))

J'ai reçu une carte du royaume-UNI de GADM (niveau 2 du royaume-UNI sur la carte).

Comment tracer une carte de chaleur sur une carte spatiale

Je voudrais être en mesure de

  1. parcelle de points définis par des coordonnées latitude/longitude sur la carte
  2. construire une carte de chaleur qui montre où les points sont plus concentrés...

Est-il facile ? Si non, avez-vous des pointeurs (seuls royaume-UNI s'il vous plaît)
Cheers

source d'informationauteur statquant

  1. 19

    Est-ce que vous aviez en tête?

    Votre sample était trop petite pour démontrer une carte de chaleur, j'ai donc créé un échantillon de plus grande taille avec des clusters au (long,lat) = (-1,52), (-2,54) et (-4.5,56). De l'OMI, la carte serait plus instructif sans les points.

    Aussi, j'ai téléchargé le fichier de formes, pas les .Rdata, et importés. La raison en est que vous êtes beaucoup plus susceptible de trouver des fichiers de formes dans d'autres projets, et il est facile de les importer dans R.
    Comment tracer une carte de chaleur sur une carte spatiale

    setwd("< directory with all your files>")
library(rgdal)         # for readOGR(...)
library(ggplot2)
library(RColorBrewer)  # for brewer.pal(...)

sample <- data.frame(Longitude=c(-1+rnorm(50,0,.5),-2+rnorm(50,0,0.5),-4.5+rnorm(50,0,.5)),
                     Latitude =c(52+rnorm(50,0,.5),54+rnorm(50,0,0.5),56+rnorm(50,0,.5)))
UKmap  <- readOGR(dsn=".",layer="GBR_adm2")
map.df <- fortify(UKmap)

ggplot(sample, aes(x=Longitude, y=Latitude)) + 
  stat_density2d(aes(fill = ..level..), alpha=0.5, geom="polygon")+
  geom_point(colour="red")+
  geom_path(data=map.df,aes(x=long, y=lat,group=group), colour="grey50")+
  scale_fill_gradientn(colours=rev(brewer.pal(7,"Spectral")))+
  xlim(-10,+2.5) +
  coord_fixed()

    Explication:

    Cette approche utilise la ggplot paquet, qui vous permet de créer des couches et ensuite rendre la carte. Les appels procédez de la manière suivante:

    ggplot -         establish `sample` as the default dataset and define (Longitude,Latitude) as (x,y)
stat_density2d - heat map layer; polygons with fill color based on relative frequency of points
geom_point -     the points
geom_path -      the map (boundaries of the admin regions)
scale_fill_gradientn - defines which colors to use for the fill
xlim -           x-axis limits
coord_fixed -    force aspect ratio = 1, so map is not distorted