1 pièce(s) jointe(s)
Décalage raster avec shape
Bonjour,
Voici mon problème: j'ai un raster et un .shape je voudrais les supersposer.
Mais ma carte donne ceci:
Pièce jointe 163395
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). :aie:
Mon raster a en principe une résolution de 50m, et je veux comme système de coordonnées du Lambert III :lol:
Voici mon code:
Code:
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:
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:
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.... :calim2: mais c'est pas normal :aie:
Quelqu'un a une idée ? :oops:
Merci
