[R-sig-Geo] problem with readOGR

2009-03-30 Thread gianni lavaredo
thanks for help. I have a polygon shapefile but I have a problem to read with readOGR sample <- readOGR(".", "sample") OGR data source with driver: ESRI Shapefile Source: ".", layer: "sample" with 766 rows and 4 columns Feature type: wkbPolygon with 2 dimensions Errore in validityMethod(objec

Re: [R-sig-Geo] Issues importing BIG shapefile - Limitations of R or something else?

2009-03-30 Thread Roger Bivand
On Mon, 30 Mar 2009, Tyler Dean Rudolph wrote: Hello Roger, This certainly does shed a lot of light on the mechanism behind the error and your code is greatly appreciated. Now for the question: did your final line of code work??? When I tried it on my machine it took over 4 hours without comp

Re: [R-sig-Geo] Repost: Find boundaries between two SpatialPolygonsDataFrame objects?

2009-03-30 Thread Jonathan Boyd Thayn
This may be useful. It is part of a function that I put together a while ago to measure the length of common borders. It's been a while since I used it. common.border <- function(shp){ if(length(slot(shp,"polygons"))!=2) stop('This function only handles two polygons at a tim

[R-sig-Geo] Repost: Find boundaries between two SpatialPolygonsDataFrame objects?

2009-03-30 Thread Peter Jepsen
Dear R-helpers I have two SpatialPolygonsDataFrame objects, but I want to plot only the boundary between them. The boundary object is to be overlaid I plot that I have made, so I would strongly prefer a method that does not rely on grid graphics. Thank you in advance for any pointers. K

Re: [R-sig-Geo] spatio-temporal cluster: R package

2009-03-30 Thread Matt Oliver
did you try require(class) ?knn there are many cluster packages available On Mon, Mar 30, 2009 at 2:04 PM, Donal Bisanzio wrote: > I'm a Italian PhD student. My research is focused on using GIS system to > study vector-borne diseases. I use open-source software, GRASS and R, to > make my analysi

[R-sig-Geo] spatio-temporal cluster: R package

2009-03-30 Thread Donal Bisanzio
I'm a Italian PhD student. My research is focused on using GIS system to study vector-borne diseases. I use open-source software, GRASS and R, to make my analysis.Now I'm trying to evaluate spatio-temporal cluster in my data. I know that in proprietary software, as Cluster-Seer, there are different

Re: [R-sig-Geo] Convergence problems with spgwr

2009-03-30 Thread Roger Bivand
On Mon, 30 Mar 2009, Adalberto Pineda wrote: Thank you for your replies but it's my understanding that (for an adaptive model) the percentage should be entered at the "adapt" argument (based on the spgwr package manual and Roger's reply to this posting: https://stat.ethz.ch/pipermail/r-sig-ge

Re: [R-sig-Geo] Convergence problems with spgwr

2009-03-30 Thread Adalberto Pineda
Thank you for your replies but it's my understanding that (for an adaptive model) the percentage should be entered at the "adapt" argument (based on the spgwr package manual and Roger's reply to this posting: https://stat.ethz.ch/pipermail/r-sig-geo/2005-August/000535.html). Also, the gwr

[R-sig-Geo] question about spatial bootstrapping

2009-03-30 Thread Satterthwaite, Paul
Hello R sig Geo list, Can anyone suggest methods for performing spatial bootstrap to obtain standard error of the semivariogram range parameter for geostatistical data? I've estimated the range parameter in semivariograms of herb cover by generating 2.5 and 97.5 percentile bounds from 1

Re: [R-sig-Geo] import ASCII file via RGDAL to SpatialPixelsDataFrame'

2009-03-30 Thread Edzer Pebesma
Try: Zone_2644149 = readGDAL("Zone_2644149.asc") This returns a SpatialGridDataFrame.Then: fullgrid(Zone_2644149) = FALSE converts into a SpatialPixelsDataFrame, or, in a single command: Zone_2644149 = as(readGDAL("Zone_2644149.asc"), "SpatialPixelsDataFrame") -- Edzer gianni lavaredo wrot

Re: [R-sig-Geo] import ASCII file via RGDAL to SpatialPixelsDataFrame'

2009-03-30 Thread Roger Bivand
On Mon, 30 Mar 2009, gianni lavaredo wrote: Dear R user, I have an ascii file (It's a raster in Arcmap) and I wish to import in R to do some analisys. I am using like help Zone_2644149 = readGDAL("Zone_2644149.asc")$band1 Zone_2644149.asc has GDAL driver AAIGrid and has 246 rows and 381 column

[R-sig-Geo] import ASCII file via RGDAL to SpatialPixelsDataFrame'

2009-03-30 Thread gianni lavaredo
Dear R user, I have an ascii file (It's a raster in Arcmap) and I wish to import in R to do some analisys. I am using like help Zone_2644149 = readGDAL("Zone_2644149.asc")$band1 Zone_2644149.asc has GDAL driver AAIGrid and has 246 rows and 381 columns str(Zone_2644149) num [1:93726] 35.1 34.7 3

Re: [R-sig-Geo] Issues importing BIG shapefile - Limitations of R or something else?

2009-03-30 Thread Roger Bivand
On Sun, 29 Mar 2009, Tyler Dean Rudolph wrote: As you've indicated Roger there does appear to be an error in readShapeLines() under maptools. However I can't help but wonder why I do not have any problems using the same command sequence for a subset of the same data (even after updating my vers