visualiser un shapfile polygone au dessus d une image ( R geospacial)
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:
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
