Dear list member, Although the answer to this might be straightforward I'm struggling with this issue for weeks now :
I have a shapefile (only one polygon) and a raster (.tif). I would like to know the value of each pixel on the raster which overlap with the polygon in the shapefile. I had various tries but none did work : library(rgdal) pol = readOGR(«shape.shp», layer = «shape») grid = readGDAL(«raster.tif») overlay(grid,pol) With this I obtain «NA» for all pixels Could you help me ? Thanks in advance, Laure Please find attached the shape and the raster Vous êtes étudiants ? Bénéficiez de Windows 7 à 35? ! Découvrez notre offre ! _________________________________________________________________ [[elided Hotmail spam]] [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo