[R-sig-Geo] From ashape to spatial polygon

2013-09-03 Thread Francesco Carotenuto
Hi, is there a way to convert an ashape object (object class of package alphahull) into a spatial polygon? I need only the alpha hull object (the external contour) and not the delaunay triangles. If not, is there a way to simply create a single spatial polygon from a matrix containing xy

Re: [R-sig-Geo] From ashape to spatial polygon

2013-09-03 Thread Roman Luštrik
You can create a SpatialPolygons object from a matrix. See help file for `SpatialPolygons-class` and sp vignette (page 7): http://cran.r-project.org/web/packages/sp/vignettes/intro_sp.pdf Cheers, Roman On Tue, Sep 3, 2013 at 12:25 PM, Francesco Carotenuto carotenuto_france...@yahoo.it wrote:

Re: [R-sig-Geo] From ashape to spatial polygon

2013-09-03 Thread Michael Sumner
I found a few discussions just by searching terms: https://stat.ethz.ch/pipermail/r-sig-geo/2012-March/014409.html http://casoilresource.lawr.ucdavis.edu/drupal/book/export/html/519 On Tue, Sep 3, 2013 at 8:25 PM, Francesco Carotenuto carotenuto_france...@yahoo.it wrote: Hi, is there a way

Re: [R-sig-Geo] From ashape to spatial polygon

2013-09-03 Thread Nicole K.S. Barker
Alpha shapes to polygons: http://rpubs.com/geospacedman/alphasimple And Michael Sumner already linked to the one I'm using - alpha hull to polygon shapefile: https://stat.ethz.ch/pipermail/r-sig-geo/2012-March/014409.html Nicole K.S. Barker Ph.D Student / Étudiante au Doctorat Laval University

[R-sig-Geo] From ashape to spatial polygon

2013-09-03 Thread Francesco Carotenuto
Thanks to everyone! [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] coordinate system in ncdf4

2013-09-03 Thread Dominik Schneider
Thanks for the tip about GDAL. I was checking the outfile with gdalinfo 1.9 and couldn't figure out why it wouldn't show the coordinate system. I'll check out raster as an alternative to ncdf4. Dominik Schneider o 303.735.6296 | c 518.956.3978 On Fri, Aug 30, 2013 at 10:54 PM, Michael Sumner

[R-sig-Geo] Analyzing Large Raster Files in R

2013-09-03 Thread mehajain13
Dear R-sig-geo users, I’m using the raster package (mainly the calc function) to calculate new indices for a large MODIS raster that I am using. The code works perfectly when run over small rasters (e.g. 100 MB), but I now want to run the process on a 60 GB raster. My computer has 20 GB of RAM

[R-sig-Geo] Creation of heatmaps with geocoded data

2013-09-03 Thread Reid Wodicka
Hello all, I am working on mapping the density of calls for service for our police department in my town. Please excuse me because I'm completely self-teaching R, so I'm probably pretty inefficient in my programming. Also, this is probably a simple answer, so I apologize for that, but I hope

Re: [R-sig-Geo] Holdout Sampling Adaptive Bandwidth SPGWR

2013-09-03 Thread Roger Bivand
yOn Fri, 30 Aug 2013, Roger Bivand wrote: On Fri, 30 Aug 2013, Paul Bidanset wrote: Thank you. I'd like to subset into a specific county. Should there be further partitioning from that level? No idea. Please re-create your scenario by subsetting georgia and the coordinates to suit.

[R-sig-Geo] Fitting variogram model using NLS function in R

2013-09-03 Thread Moshood Agba Bakare
Dear All, I tried to compare the result obtained by fitting the spherical variogram model using fit.variogram and nls functions. The large difference in the results is a great concern for me knowing that the two functions use Weighted Least Squares (WLS) approach for estimating parameters. The

Re: [R-sig-Geo] Holdout Sampling Adaptive Bandwidth SPGWR

2013-09-03 Thread Paul Bidanset
Thank you very much for the example and the clarification. My hold out test is random. The vector provided by gw.adapt() allows me to see the bandwidth size for each point. Is there a way to see each regression point's bandwidth size with the correct format you just showed me? On Tue, Sep 3,

[R-sig-Geo] Merge maps and data.frames using loop

2013-09-03 Thread Arnold Salvacion
Dear Colleagues, Good day! Anyone had experience replacing map attribute/data frame (i.e. map@data) with other data frame using loop. I have maps of A,B, and C and I have a data frames of a1,b1, and c1. How can I replace/merge  the data frames them using loop? So far below is the code that I