Dessiner des carrés de 0.5°*0.5° sur carte

[Arnaud PRUVOST]

Bonjour,

Je souhaiterai dessiner des carrés de 0.5°*0.5° sur une carte (un peu comme des pixels).

Mes données sont les suivantes

Code : Tout sélectionner

> head(F1_column_R1)
   lon  lat     x
1 44.5 38.5  3600
2 45.0 38.5 10800
3 45.5 39.0  7200
4 46.0 39.0 10800
5 53.0 39.0 10800
6 53.5 39.0 18000


lon et lat représentant respectivement les longitudes et les latitudes et x représentant un temps de résidence en seconde.

En fait, chaque coordonnée représente une aire :
par exemple pour la première ligne les coordonnées sont lon=44.5 et lat=38.5 mais cela correspond à la surface suivante :
longitude de 44.5 à 45(exclu)
latitude de 38.5 à 39(exclu).

Je ne pense pas que ce soit possible sous geom_point. La carte utilisée jusqu'à présent provient de openmap. Peut etre faut-il que je télécharge une carte shp (si oui comment?) pour plus de facilité ?

Avez-vous une idée ?

[Mickael Canouil]

Bonjour,

les fonctions geom_tile() et geom_raster() ne conviennent pas?

[Arnaud PRUVOST]

En supposant que je ne me trompe pas dans le code,

Code : Tout sélectionner

p<-autoplot(map2,zoom=10) +
  geom_raster(data=F1_column_R1,aes(x=lon,y=lat,col=x),cex=1.2)+
scale_color_gradient(low=c("white","violet","darkblue","skyblue","yellow","orange"),high="red",guide=guide_colorbar(barwidth=1.5,barheight=9),space="Lab",name="t (s)", limits = c(0,1200000))+
xlab("Longitude") + ylab("Latitude")


J'obtiens ça (donc impossible à voir si le résultat est celui attendu).


J'ai pas modifier hjust et vjust car ils sont par défaut à 0.5
Pourquoi l'échelle de couleur ne fonctionne plus dans ce cas là ?

EDIT : en faisant cela, je me suis aperçu qu'il fallait faire du kriging. Je vais créer un nouveau topic pour celà.

[Renaud Lancelot]

Il y a bcp de manières différentes d'obtenir le résultat selon le type (classe) d'objet manipulé. Par exemple:

Code : Tout sélectionner

## Utilisation d'objets de classe map
library(maps)
library(mapdata)
library(mapproj)
## exact et fill nécessaires pour la conversion ultérieure (2e graphe)
fr <- map(region = "France", data = "worldHires", exact = TRUE, fill = TRUE, col = NA)
class(fr)
## Etendue de la grille
rg <- c(-5, 8.5, 42, 51.5)
## Trace la grille
map.grid(lim = rg,
         nx = diff(rg[3:4]) / .5,
         ny = diff(rg[1:2]) / .5,
         labels = FALSE, lty = 1, col = "grey")

## Utilisation d'objets de classe SpatialPolygons
library(sp)
library(raster)
library(maptools)
## Convertit map en SpatialPolygons
FR <- map2SpatialPolygons(fr, IDs = "France", proj4string = CRS("+proj=longlat"))
class(FR)

## crée une grille à partir de l'étendue
e <- extent(rg)
r <- raster(e,
            ncols = diff(rg[3:4]) / .5,
            nrows = diff(rg[1:2]) / .5,
            crs = CRS("+proj=longlat"))
R <- as(r, "SpatialPolygons")

## Trace
plot(FR)
plot(R, add = TRUE, border = "grey")