[R-sig-Geo] Finding the highest order of neighbors

Hi rgeo list! Considering a neighbor object, how can I find the highest order of neighbors and which node corresponds to it? I think of using spdep::nblag iteratively until I find no change between the lower order and the actual one, but I find the method a bit dirty. Would you have any idea?

[R-sig-Geo] kernel2d - raster/grid

Hi, how can I convert kernel2d estimation to a raster in order to use it as a mapserver layer? I saw maptools can save kernel2d to contour to shape file with ContourLines2SLDF. Any hints? Regards -r -- *L'unica speranza di catarsi, ammesso che ne esista una, resta affidata all'istinto di

[R-sig-Geo] rgeos: gOverlaps not working properly?

Hello, I have following problem with gOverlaps function: I have two lines, A and B. line A has 57 vertexes, all identical with vertexes of line B yet, gOverlaps(A,B) gives FALSE the code follows... CooA # coordinates of line A [,1] [,2] [1,] -506157.4 -1122381 [2,] -506046.4

Re: [R-sig-Geo] rgeos: gOverlaps not working properly?

Would it be possible to provide your output in a form that we can read it? I suggest dput(). Cheers, Roman On Mon, Oct 15, 2012 at 2:20 PM, zbynek.jano...@centrum.cz zbynek.jano...@centrum.cz wrote: Hello, I have following problem with gOverlaps function: I have two lines, A and B. line

Re: [R-sig-Geo] rgeos: gOverlaps not working properly?

Absolutely, sorry for that dput(CooA) structure(c(-506157.3827, -506046.35779, -506021.42681, -505984.87158, -505953.70149, -505910.16251, -505883.4289, -505840.08242, -505800.22761, -505771.72411, -505735.3649, -505701.88969, -505651.2245,

Re: [R-sig-Geo] rgeos: gOverlaps not working properly?

‘gOverlaps’ returns TRUE when the geometries share some but not all interior points, and the intersection has the same dimension as the geometries themselves. Hence gOverlaps(A,A) == FALSE due to that 'but not all' part of the documentation. Barry

Re: [R-sig-Geo] rgeos: gOverlaps not working properly?

Thank you very much, problem solved. Zbynek -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/rgeos-gOverlaps-not-working-properly-tp7581278p7581282.html Sent from the R-sig-geo mailing list archive at Nabble.com. ___ R-sig-Geo

[R-sig-Geo] problem creating psp

Hello Everyone, I am currently trying to complete the simple task of converting a SpatialLinesDataFrame to PSP so that I can do kernel density estimation on the line pattern (density.psp). When I do this, however, I get an error that does not allow me to proceed. Here is the code I am using

Re: [R-sig-Geo] problem creating psp

It works for this example, so either something is wrong with your data or the packages are out of date. Can you provide summary(SierraLeone.roads) and/or the data file, and sessionInfo()? library(maptools) xx - readShapeLines(system.file(shapes/fylk-val.shp, package=maptools)[1],

Re: [R-sig-Geo] spml, spfeml and tolerance level

After some exploration and a hint from the package administrator, the quickest way I found to avoid this problem (without having to customize functions within the package) was setting: .Machine$double.eps - .Machine$double.eps^2 Just in case anyone faces this. cheers. AOB - Ariel

[R-sig-Geo] Subset SpatialGridDataFrame by retrieving xy coordinates from bounding box

Hello ! I am trying to subset elements from a SpatialGridDataFrame not by attribute, nor by indexing of rows/columns such as X[rows, columns], but by retrieving the xy coordinates from the bounding box. Example : bbox(SGDF) min max x 457446 457512 y 5384573 5384637 The

Re: [R-sig-Geo] Subset SpatialGridDataFrame by retrieving xy coordinates from bounding box

This can be easily done in the spatstat package if you convert the object of interest to a spatial point pattern (object of class ppp). See the vignette shapefiles for some guidance on how to do the conversion. cheers, Rolf Turner On 16/10/12 13:54, Nathalie Morin wrote: Hello !

Re: [R-sig-Geo] Subset SpatialGridDataFrame by retrieving xy coordinates from bounding box

Natalie, library(raster) b - brick(SGDF) e - extent(457446, 457506, 5384573, 5384633) bb - crop(b, e) # not sure what you want to do next, perhaps sge - as(bb, 'SpatialGridDataFrame') Robert On Mon, Oct 15, 2012 at 5:54 PM, Nathalie Morin nathalie.m.mo...@usherbrooke.ca wrote: Hello !