Re: [R-sig-Geo] Newbie Question...

2010-01-04 Thread Edzer Pebesma
Michael Sumner wrote: Here's one way. I'm not sure about how to control the aspect ratio (something I've been meaning to check for a while now). somewhat of a non-newby answer to Michael: asp is set in plot.Spatial, and if the default value of NA is passed, set to: if (is.na(asp)) asp

Re: [R-sig-Geo] Newbie Question...

2010-01-04 Thread Michael Sumner
Here's one way. I'm not sure about how to control the aspect ratio (something I've been meaning to check for a while now). library(sp) data(volcano) x <- 10*(1:nrow(volcano)) y <- 10*(1:ncol(volcano)) ## image xyz list imlist <- list(x = x, y = y, z = volcano) ## SpatialGridDataFrame (image2Grid

[R-sig-Geo] Newbie Question...

2010-01-04 Thread Jason Rupert
At the suggestion of Walmes Zeviani I downloaded and am taking a look at the sp package. In an attempt to get going, can someone offer a suggestion about how to convert the below over to a version that uses the sp package? data(volcano) x <- 10*(1:nrow(volcano)) y <- 10*(1:ncol(volcano)) imag

[R-sig-Geo] PBSmapping importShapefile question

2010-01-04 Thread Julie Lee-Yaw
Hi again, Thanks to those who helped me with my earlier post. I've installed rgdal and have made some progress but I'm now stuck on an error message that I don't know how to interpret. I have imported one shapefile that is projected in the CRS that I want. I then import another shapefile that

Re: [R-sig-Geo] Creating random (but believable) geographies

2010-01-04 Thread Barry Rowlingson
Here's a link to a blog entry I've just written on my code for doing this: http://geospaced.blogspot.com/2010/01/generating-fake-geographies.html The code isn't available yet, it needs a bit of cleaning up and documenting, but if anyone wants the raw code then I can scrape it up for you! Barry

Re: [R-sig-Geo] newbie rgdal projection problem

2010-01-04 Thread Roger Bivand
On Mon, 4 Jan 2010, Peter Adler wrote: Dear R-sig-geo, I am new to RGDAL and PROJ4, and am having trouble projecting a raster dataset. The dataset is a 721 x 721 cell raster of snow cover from the NSIDC. The projection information is here: http://nsidc.org/data/atlas/epsg_3408.html My com

[R-sig-Geo] newbie rgdal projection problem

2010-01-04 Thread Peter Adler
Dear R-sig-geo, I am new to RGDAL and PROJ4, and am having trouble projecting a raster dataset. The dataset is a 721 x 721 cell raster of snow cover from the NSIDC. The projection information is here: http://nsidc.org/data/atlas/epsg_3408.html My command: D=readGDAL("annual_snow_trend_ease

Re: [R-sig-Geo] color category SpatialGridDataFrame and rworldmap

2010-01-04 Thread Andy South (Cefas)
Dear Patrick & Renaud, Categorical maps are one of the issues we have looked at in a new package 'rworldmap' that we've been developing. Coincidentally I submitted the first version to Cran for consideration yesterday. The package and source can currently be downloaded from : http://code.google