Re: [R-sig-Geo] gridded time series analysis

2010-11-30 Thread Robert J. Hijmans
Advait, The is an error in the file, I think, see below: > library(ncdf) > filename = "c:/downloads/regcm.cumul.hdd.10k.96_05.new.nc" > nc <- open.ncdf(filename) > > zvar <- 'hdd' > nc$var[[zvar]]$dim[[1]]$len [1] 507 > nc$var[[zvar]]$dim[[2]]$len [1] 356 > > xx <- nc$var[[zvar]]$dim[[1]]$vals >

Re: [R-sig-Geo] geostatistic memory issue

2010-11-30 Thread Robert J. Hijmans
Oops, that has been fixed on r-forge, but not CRAN. Please try to install from here (you need the latest version of R for that) like this: install.packages("raster", repos="http://R-Forge.R-project.org";) I will update to CRAN this week. Robert On Tue, Nov 30, 2010 at 12:32 PM, Fernando Santo

[R-sig-Geo] geostatistic memory issue

2010-11-30 Thread Fernando Santo
Dear Giuseppe, I am using 2 processor of 2 Gb. I did not know about this problem. Is there a way to allocate more memory in linux? I believe that the answer is in R help ?Memory. But I did not understand how to do it. Starting R from the shell with R --min-vsize=100M --max-vsize=3G --min-

Re: [R-sig-Geo] How to draw a kernel surface map in R?

2010-11-30 Thread Michael Sumner
See ?image, ?persp in the graphics package, ?levelplot and ?wireframe in the lattice package, and ?surface3d in the rgdal package. The sp package has methods for the graphics and lattice functions if you use Spatial objects. If you want functions for kernel density, see ?kde2d in MASS, and the sp

[R-sig-Geo] geostatistic memory issue

2010-11-30 Thread Fernando Santo
Dear Robert, Thanks a lot for your attention. I downloaded the raster package and manual. I am still digesting all this information... I just tried to do a quick adaptation in my code and it did not worked. I am doing something wrong... Fernando # read the data dat <- read.csv("/home/fern

Re: [R-sig-Geo] gridded time series analysis

2010-11-30 Thread Advait Godbole
Dear Robert, here is the link to the file: http://dl.dropbox.com/u/4030944/regcm.cumul.hdd.10k.96_05.new.nc. It is a large file ~400 MB. Let me know if dropbox is not convenient and you would rather do ftp - but our machines are down for maintenance and I may not be immediately able to get access

[R-sig-Geo] back-transforming lognormal kriging

2010-11-30 Thread Burton Shank
Hello List (and Jorge), I’m looking for advice on back-transforming kriged values. I know this has previously been discussed on the list but haven’t found a satisfactory solution for my application. Apologies if I’ve missed it somewhere. I’m trying to estimate the density of two species at

Re: [R-sig-Geo] geostatistic memory issue

2010-11-30 Thread Robert J. Hijmans
On the multi-core topic: Matteo Matiuzzi and I have started building multi-core processing into 'raster' functions. Only testing it with a few functions right now. (see ?beginCluster in the R-Forge version of raster). 'interpolate' or 'predict' do not have this capacity yet, but they are very high

[R-sig-Geo] How to draw a kernel surface map in R?

2010-11-30 Thread Everton Emanuel
Dear all, I am writing a paper for publication and I have problems to draw kernel surface maps in R. In the descriptive part of the paper I would like to place some thematic kernel surface maps, but I don't no how to do it. Is there a way to do with the traditional libraries like: spdep or maps?

Re: [R-sig-Geo] warning at creating SpatialPolygonsDataFrame and error at saving through rgdal

2010-11-30 Thread Agustin Lobo
The error in wirteOGR() was caused by the long names of the variables in the table of the SPDF object that were listed in the warnings. Fixed by renaming some columns in the data slot to shorter names. The error was not related to the warnings issued at creating the SPDF object by my function, tho

Re: [R-sig-Geo] geostatistic memory issue

2010-11-30 Thread Giuseppe Amatulli
Yes, in almost all "prediction" command you can apply to a tiles grids, as is suggested by Robert. Anyway be sure that you are using all your 4 G. If you have 2 processor of 2 G, and you run R in a shell you will use only the 2 G of one processor. In this case the multicore/parallel packages can

Re: [R-sig-Geo] raster::crop

2010-11-30 Thread steven mosher
Ya, the codes done, I was just checking out the changes to the source before doing an update and saw that he had changed extent. On Mon, Nov 29, 2010 at 11:48 PM, Agustin Lobo wrote: > Thanks, > > and summarizing messages exchanged with Robert, > > subs = bN3[1000:1100,1000:1100, drop=FALSE] > >