Re: [R-sig-Geo] External Envelope of a set of points in space - not convex hull

2012-06-27 Thread Mathieu Rajerison
Hi, alphahull is nice but having a SpatialPolygonsDataFrame from the result is quite tricky If anyone has the solution.. 2012/6/22 Michael Sumner mdsum...@gmail.com Try the alphahull package, which is on CRAN. I don't believe it will smooth the result for you, but you could take it from

[R-sig-Geo] breaking lines at intersection but without GRASS

2012-06-27 Thread Mathieu Rajerison
Hi, My final goal is to generate a graph from a road network. I didn't manage to install shp2graph correctly, so I want to do it from scratch. Is there a way to break lines at intersection but without rgrass? Also, apart from shp2graph which seems a very good package, does one exist that would

Re: [R-sig-Geo] interpreting SAR and CAR outputs

2012-06-27 Thread Corey Sparks
Hi, I would examine the Lagrange multiplier specification tests from a lm() fit. like: fit1-lm(yield~Temp+precept, data=data) lm.LMtest(fit1, listw=ND_queenW, test=all) and see which of the tests has the largest significant test statistic. That will let you know if your data follow a lag or error

Re: [R-sig-Geo] ordinary kriging using 'krige.conv' function

2012-06-27 Thread Giuseppe Calamita
I never worked with geoR package for interpolation pourposes. Have you already checked the gstat package ? Further, do you have any idea of which model would be the best for fitting the data you're working on (exponential, spherical, gaussian, etc ?). Give me more information so I can help you

Re: [R-sig-Geo] External Envelope of a set of points in space - not convex hull

2012-06-27 Thread Michael Sumner
Can you share what you've tried? A triangulation is really at odds with the representation available with polygons in their usual GIS form, but there are ways to do it. The basic options (to me) are either a set of the triangles as individual polygons, or even a set of shapes in a multipolygon, or

[R-sig-Geo] Problem with crop() on raster stacks

2012-06-27 Thread Klaus Jacobi
Dear Robert, I have a script where I am cutting large raster stacks into slices by using the crop function. Until recently, the script worked without a problem. In the mean time, I upgraded from R version 15.0 to 15.1 and from raster version 1.9-x to 2.0-x. Today I noticed, that I still can

[R-sig-Geo] costDistance with gdistance error

2012-06-27 Thread Ian Craig
Hello geo-savvy R people, I am new to this list so please be forgiving if I don't include everything that is needed. I am having trouble with some code I wrote using the gdistance package. #I loaded an ESRI ascii grid file cost - read.asciigrid(filepath, as.image = F, plot.image=F,, proj4string

Re: [R-sig-Geo] costDistance with gdistance error

2012-06-27 Thread Jacob van Etten
The raster is probably too big? -- difficult too know without access to your file.  A limitation of gdistance is that everything should be done in memory. You could try to aggregate the grid. You should remove all NA values from the raster. Running trim() is not enough.  You can use standard