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