[R-sig-Geo] Create a simple spatial grid

[Remi Genevest]

Alright !
You pointed out my goal indeed... but 100meters by 100meters is still a bit
too high for a map over Europe.
So, this is what I have done :

### create a grid raster 
ll = CRS("+init=epsg:4326") 
origin = SpatialPoints(cbind(7,40),ll) 
grid = raster() 
## make resolution ie 0.008 min par 0.008 spaced cells (equals to about
800m*800m), that is actually quite accurate 
e = extent(5,18,40,50)
extent(grid)<-e
res(grid)=c(0.008,0.008) 
grid

Now, I hope that the extract() function will work fine after adding my other
layers (that do not have the same original CRS nor resolution...!)

Many thanks !!
Rémi





--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Create-a-simple-spatial-grid-tp7586481p7586496.html
Sent from the R-sig-geo mailing list archive at Nabble.com.