bonjour,
je travaille sur R sur un projet goérospatiale
je sais pas si quelqu'un peut m'aider , mais j'arrive pas
le problème que j'ai une image d'un indice MODIS EVI , et shapfile de type polygone : cercle , qui ont même projection ,
l'objectif ce que je dois visualiser le polygone sur l'image mais j'arrive pas j'ai utilisé les deux dernières lignes dans mon script mais marche plus il affiche que l'image :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21 library(raster) library(rgdal) library(tiff) r <- raster("D:/Rprojet/shapfile/MODIS/testmodis2.tif") str_name<-'testmodis.tif' imported_raster=raster(str_name) plot (r) # importer les fichiers shapefil data.shape<-readOGR(dsn="D:/Rprojet/shapfile/",layer="2001") plot(data.shape,border="red") # afficher polygone a dessus du raster plot(r) plot(data.shape,col="red",border="black",add = TRUE)
merci d'avance
Partager