Discussion :

[mauriceletendu]

Bonjour,

Voici mon problème: j'ai un raster et un .shape je voudrais les supersposer.

Mais ma carte donne ceci:


Mon raster a l'air déformé ou décalé. Pour infos, c'est la région PACA (pour le raster il y a un peu plus que la région PACA, vers l'ouest et le nord, c'est normal que le raster dépasse les départements).

Mon raster a en principe une résolution de 50m, et je veux comme système de coordonnées du Lambert III

Voici mon code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
 
 
extent(q) <- c(736200,1030750,83200,329600)
projection(q) <- CRS("+proj=lcc +lat_1=49 +lat_2=44+lat_0=46.5 +lon_0=3 +x_0=700000+y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
projection(g) <- CRS("+proj=lcc +lat_1=49 +lat_2=44+lat_0=46.5 +lon_0=3 +x_0=700000+y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
plot(q)
par(new=T)
plot(g)

q est mon raster, g mon .shp

Les coordonnées extent, je l'ai ai calculées, il est possible qu'il y est une erreur.... mais je pense pas...

Voici les paramètres du raster:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
class       : RasterLayer 
dimensions  : 4929, 5891, 29036739  (nrow, ncol, ncell)
resolution  : 50, 49.98986  (x, y)
extent      : 736200, 1030750, 83250, 329650  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=lcc +lat_1=49 +lat_2=44+lat_0=46.5 +lon_0=3 +x_0=700000+y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
data source : in memory
names       : layer 
values      : 0, 4086  (min, max)

et du shp

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
 
features    : 6 
extent      : 753235.1, 1030703, 83180.27, 322270.4  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=lcc +lat_1=49 +lat_2=44+lat_0=46.5 +lon_0=3 +x_0=700000+y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
variables   : 2
names       :                   nom, code 
min values  : ALPES-DE-HTE-PROVENCE,   04 
max values  :              VAUCLUSE,   84

Ma résolution du raster n'est pas de 50,50 mais de 50, 49.98986 (x, y) , je sais pas pourquoi.... mais c'est pas normal

Quelqu'un a une idée ?

Merci