How to make "meuse" data generates "meuse.grid" data?
Edzer Pebesma wrote:
Jessica, I don't know the exact details, either I wasn't there or it is too long ago. If you have a polygon that encloses your study region, look into using method ovelay to find out which grid points of the full (rectangular) grid fall inside your polygon.
Hi all, or alternatively, sample directly from the polygon using spsample, see also the e-mail I replied ealier to Jessica's question. quote: Hi, Loading the meuse.grid dataset is done using: data(meuse.grid) to generate meuse.grid from meuse use something like: data(meuse) coordinates(meuse) = ~x+y meuse_grid = spsample(meuse, type = "regular", cellsize = c(40,40)) gridded(meuse_grid) = TRUE unquote
for the sig/wiki: is this a FAQ as well? -- Edzer jessica243 wrote:
Dear Edzer, I am sorry to writing again,last time I ask you how to make"meuse" data generates "meuse.grid" data? you told me use "SpatialPixels( Spatialpoints(makegrid(meuse)))".I found using this method, the results do not the same as meuse.grid.
Could you define "not the same". cheers, Paul
could you tell me the details?
Sincerely
Jessica
Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul