[R-sig-Geo] Changing of colour scale

Good morning to everybody. I used pointsToRaster in raster package to do a geographic plot. Green represented the maximum count in the colour bar. But I wish to reverse this and choose a different colour scale where red would present the maximum counts. I would be glad if anyone could guide me on

Re: [R-sig-Geo] Producing ENFA derived suitability maps with adehabitat

Dear Mathieu, Thanks for the clarification. I have underestimated the complexity, but I admit - habitat analysis is not really 'my cup of tee' (maybe this question was more suited for the Environmentrics R-sig). As least I was correct about the whole concept - it is doable ;) I would be

Re: [R-sig-Geo] Changing of colour scale

if you do ?plot you can see that the default plot for raster is plot(x, col=rev(terrain.colors(255)) ) so what you want is not the reverse (rev), but: plot(x, col=terrain.colors(255) ) you can have a look at ?terrain.colors and the links therein to learn about other color palettes Robert On

Re: [R-sig-Geo] DEMs

Hi Abe, Did you try using the ETOPO1? This (0.1 degree) DEM can be well loaded to R: URL - http://spatial-analyst.net/worldmaps/; # list of maps: map.list - c(globedem, chlo08, dcoast) # download the zipped maps one by one: for(i in 1:length(map.list)) { download.file(paste(URL,

Re: [R-sig-Geo] DEMs

I consider GTOPO30 to be obsolete because of SRTM30, at least for tropical mountain areas where GTOPO30 is pretty bad in some areas (although not perhaps not that much different when aggregated to a 1 degree resolution)... There are a few different SRTM30 products available. On the worldclim

[R-sig-Geo] Spatial clustering of points within a given line

Julius Tesoro jutes...@yahoo.com writes: I have a dilemma, Considering I have a number of points and a line (shapefile). How do I compute for the shortest distance between a line and a point? In package 'spatstat' you would use nncross(X,Y) or project2segment(X, Y)

Re: [R-sig-Geo] Producing ENFA derived suitability maps with adehabitat

Dear Tom (and other list members), First of all, sorry not replying earlier to you Tomislav but I didn't receive your message (certainly because of a local serveur problem). For the moment, the method described by Mathieu applied to our dataset seems coherent with what we know about changes

Re: [R-sig-Geo] Spatial clustering of points within a given line

Thanks for that. Is the distance of the point to this 'line segment' perpendicular? Also, is there a package where you can do automated spatial clustering? --- On Mon, 9/14/09, adrian.badde...@csiro.au adrian.badde...@csiro.au wrote: From: adrian.badde...@csiro.au adrian.badde...@csiro.au

Re: [R-sig-Geo] Classes for spatio-temporal data?

Dear Christian and Edzder, I'm glad to hear about this, because I'm also thinking about extending in some way the SpatialPointsDataFrames class of the 'sp' package, for making easier the spatio-temporal interpolations within R. So far I'm only working with regular time point data, with

Re: [R-sig-Geo] DEMs

Hi Tomislav, Thanks for the detailed response. I tried to use the globedem.zip from your website, but then when I call overlay() with my data, there are many NAs produced in the output: it seems the DEM covers X= -180 to 180 Y= -65 to 65, my dataset covers X= -180 to 180 Y= -90 to 90 So when

Re: [R-sig-Geo] DEMs

Hi Robert, Thank you for the detailed information. I downloaded the 10 min SRTM30 DEM but was unable to run the aggregate command.. The error message was : Error in .local(x, ...) : object 'startrow' not found I tried to use readGDAL() but then ran into the same issue with NAs after the call