Re: [R-sig-Geo] Wiki?

Dear Jonathan, please try this link: http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data Best, Friderike -Original Message- From: r-sig-geo-boun...@stat.math.ethz.ch [mailto:r-sig-geo-boun...@stat.math.ethz.ch] On Behalf Of Jonathan Greenberg Sent: 20 February 2009 05:39 To: r

Re: [R-sig-Geo] world map data

For data of a world map of administrative boundaries, I can also recommend the "Global Administrative Unit Layers" from the FAO Geonetwork catalog http://www.fao.org/geonetwork/ Just enter "GAUL" in the search window. The data exist from 1990 to 2008. Cheers, Friderike -Original Message-

Re: [R-sig-Geo] select polygons in polygons

Try ?overlay Regards, Friderike -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Murray Richardson Sent: 25 September 2008 06:01 To: r-sig-geo@stat.math.ethz.ch Subject: [R-sig-Geo] select polygons in polygons Hello users, There is probably a very sim

[R-sig-Geo] R-SAGA-GRASS-ArcGIS strengths/weaknesses

If a non-experienced person wanted to decide for one program in which to invest time and if one was to highlight the strengths and weaknesses of R, SAGA, GRASS and ArcGIS in comparison to each other in a simplistic way ... R SAGAGRASS ArcGIS geostatistical anal

Re: [R-sig-Geo] How to remove polygons according to their attributes?

Hi Lin, I think you can simply do your selection as with any dataframe, e.g. require(rgdal) Ctries=readOGR(dsn='C:/data/datasets/HornOfAfrica', layer='HornOfAfrica', p4s=NULL) str([EMAIL PROTECTED]) try<-Ctries[Ctries$AREA>100,] spplot(try, zcol='ADM0_CODE') # or: Kenya<-Ctries[Ctries$ADM0_NAME==

Re: [R-sig-Geo] Defining a grid for interpolations ?

Dear Mauricio, with some help of Dr. Rossiter, I once created a grid covering my study area by importing a vector file delimiting the borders of the area and overlaying it over the expanded grid: # import the vector shape-file area=readShapePoly('/.../my_area') # in case that the vector file inc

Re: [R-sig-Geo] readOGR and proj4 string

Not sure whether this is helpful, but I found that the same problem has been reported elsewhere: http://osdir.com/ml/lang.r.geo/2006-09/msg00040.html Importing with readShapePoly does not recognize existant .prj files neither, but I guess that in this case it is not supposed to do so? Friderike (

Re: [R-sig-Geo] newbie: spplot/ colour ramps

Thanks Edzer, I shall try again: 1) I would like to use spplot to map my factor "TYP" as dots of different colours, however the resulting plot uses the same colours for the first and last value (10,40). I guess that my use of the "cuts" argument is wrong, but I can't find any better: LAT <-c(-6.

[R-sig-Geo] newbie: spplot/ colour ramps

Dear list members, I am still struggeling with spplot. Sorry for being that unprofessional, but honestly I find the spplot help page quite hard to understand. (Many thanks Edzer for the very helpful hint a few days ago.) 1) My first plot actually works all fine, but the colours seem to loop, i.e.