Re: [R-sig-Geo] Points density in R

Kernel density can be derived in spatstat package, but also in adehabitat and many other packages. You need to loop your operations (output as a list). Note that I use the same file names to save space - you only need the output of course! *** library(maptools) library(rgdal)

[R-sig-Geo] info about package divagis

Hello List, I was hoping that someone on this list might have some more information about the 'divagis' a package. http://cran.r-project.org/web/packages/divagis/ The site says that it 'Provides tools for quality checks of georeferenced plant species accessions'. I am very interest in this

[R-sig-Geo] Geonames lookup (was lat/lon to timezone)

My R interface to GeoNames is now on R-forge: http://geonames.r-forge.r-project.org/ Assuming I haven't mucked anything up the package binaries should get built later, otherwise you can check out the SVN and play with it: svn checkout svn://svn.r-forge.r-project.org/svnroot/geonames Barry

[R-sig-Geo] map2SpatialPolygons' 'proj4string' in spmaps

Hi, I don't understand what the 'proj4string' argument in map2SpatialPolygons() does. map() in the maps package has a 'projection' argument: ,-[ *help[R](map)* ] | projection: character string that names a map projection to use. See | 'mapproject' (in the 'mapproj' library). The

[R-sig-Geo] Uniras routin in R

Hi, Uniras has a nice interpolation routin that is build on Akimas ideas (a set of irregularty distrubated data points to a regular grid of points). The akima package, is quite close but I missing some variables to adjust my data (or the look of my output grid). In uniras the gridding is done

Re: [R-sig-Geo] map2SpatialPolygons' 'proj4string' in spmaps

On Wed, 15 Oct 2008, Sebastian P. Luque wrote: Hi, I don't understand what the 'proj4string' argument in map2SpatialPolygons() does. map() in the maps package has a 'projection' argument: ,-[ *help[R](map)* ] | projection: character string that names a map projection to use. See |

Re: [R-sig-Geo] map2SpatialPolygons' 'proj4string' in spmaps

On Wed, 15 Oct 2008 19:24:18 +0200 (CEST), Roger Bivand [EMAIL PROTECTED] wrote: [...] No, not at all. The coordinates returned by map() if projection= is used are only for plotting, are in arbitrary units, and are only documented in code. map2SpatialPolygons() assumes that the data are in

Re: [R-sig-Geo] map2SpatialPolygons' 'proj4string' in spmaps

On Wed, 15 Oct 2008, Sebastian P. Luque wrote: On Wed, 15 Oct 2008 19:24:18 +0200 (CEST), Roger Bivand [EMAIL PROTECTED] wrote: [...] No, not at all. The coordinates returned by map() if projection= is used are only for plotting, are in arbitrary units, and are only documented in code.

Re: [R-sig-Geo] map2SpatialPolygons' 'proj4string' in spmaps

On Wed, 15 Oct 2008 20:20:16 +0200 (CEST), Roger Bivand [EMAIL PROTECTED] wrote: On Wed, 15 Oct 2008, Sebastian P. Luque wrote: On Wed, 15 Oct 2008 19:24:18 +0200 (CEST), Roger Bivand [EMAIL PROTECTED] wrote: [...] If the fill=TRUE argument to map() is omitted, the rings are not built:

[R-sig-Geo] R+SAGA formula to subtract two raster

Hi All, I need help to write a code in R to subtract two raster in SAGA. I have DEM_1_f.sgrd and DSM_1_f.sgrd in SAGA format and I wish to subtract (DSM_1_f.sgrd - DEM_1_f.sgrd) to obtain a new layer A part of the code is this: #A# DEM - 1) Shapes to Grid

Re: [R-sig-Geo] points on polygons

On Wed, Oct 15, 2008 at 6:06 PM, Thomas Lumley [EMAIL PROTECTED] wrote: I am trying to draw a US map with color-coded cities superimposed on color-coded states, showing survey results from the Behavioral Risk Factor Surveillance System. You can see the sort of thing I'm trying to do at

Re: [R-sig-Geo] points on polygons

On Wed, 15 Oct 2008, hadley wickham wrote: It's relatively easy to do this with ggplot2 (the tricky part is converting the polygons to a data frame) - I can provide some code if you're interested. That would be great. -thomas Thomas Lumley Assoc. Professor,