Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-12 Thread v v
Hi, Paul, Thank you for the reply. I have attempted to write some code to test out what you've suggested but ran into some errors. library(fields) k - Krig(ozone$x, ozone$y, theta=20) # define objective function kfunc - function(para,k) { xcoord=para[1] ycoord=para[2] minimum -

Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-12 Thread Paul Hiemstra
On 06/12/2010 09:11 AM, v v wrote: Hi, Paul, Thank you for the reply. I have attempted to write some code to test out what you've suggested but ran into some errors. library(fields) k - Krig(ozone$x, ozone$y, theta=20) # define objective function kfunc - function(para,k) {

Re: [R-sig-Geo] How to extract coordinates values from a shapefile?

2010-06-12 Thread Paul Hiemstra
On 06/11/2010 11:43 AM, Paul Hiemstra wrote: On 06/09/2010 11:01 PM, Dylan Beaudette wrote: On Wednesday 09 June 2010, Hadley Wickham wrote: On Wed, Jun 9, 2010 at 3:24 PM, Edzer Pebesma edzer.pebe...@uni-muenster.de wrote: The example provided by Matt assumes that each polygon consists of

[R-sig-Geo] Combining multiple polylines into a single polyline

2010-06-12 Thread Don MacQueen
I have an input shapefile representing roads in a city, read into R where it becomes a SpatialLInesDataFrame. Looking at just one road at a time, it looks fine when plotted. But in the underlying structure, the single real-world road has been broken up into many shorter polylines, and they