Step characteristics on raster maps (Mathieu Basille)
I don't know. I can't reproduce your example, and a simple example on the meuse data set: library(sp) library(fields) data(meuse) Tps(data.frame(meuse$x,meuse$y), meuse$zinc) seemed to work. You could provide a reproducable example, or you could contact the author(s) of the fields package, I believe they're not listening to this list. -- Edzer
On 10/03/2010 05:59 PM, Mathieu Basille wrote:
Dear Arnaud, This is actually on my todo list, a test with the polygon approach, although I'm not sure how to start to "polygonize" the segments... But in the end, that would allow a performance comparison with different approaches (and that would also allow the calculation of the characteristics of buffers around the segments, which might be interesting for us, too). Thanks for the suggestion, Mathieu. Le 03/10/2010 10:04, Arnaud Mosnier a ?crit :
Hi Mathieu, An easy workaround could be to convert each segment to a polygon doing a very small buffer around each segment then use the overlay method to obtain info in each of these polygon. Arnaud
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de