Re: [R-sig-Geo] Using multicore processors with gstat

2010-07-20 Thread Edzer Pebesma
On 07/21/2010 05:27 AM, Pierre Roudier wrote: > Dear geostatisticans, > > I'm going to deal with pretty big datasets in one of my projects (1G > pixels +), and I'd like to keep using gstat for the geostatistics. > > Do someone know if the gstat package supports parralelisation of > functions li

[R-sig-Geo] Using multicore processors with gstat

2010-07-20 Thread Pierre Roudier
Dear geostatisticans, I'm going to deal with pretty big datasets in one of my projects (1G pixels +), and I'd like to keep using gstat for the geostatistics. Do someone know if the gstat package supports parralelisation of functions like cariogram() and fit.variogram() on multicore processors? C

Re: [R-sig-Geo] display projected NetCDF files

2010-07-20 Thread Robert J. Hijmans
Walt, You should be able to read and plot this file like this: library(raster) b = brick("air.mon.mean.nc", varname= 'air') plot(b, 1) #or r = raster("air.mon.mean.nc", varname= 'air', band=1) plot(r) #or sp = as(b, 'SpatialGridDataFrame') #etc The projection info is stored in r...@prj (for

Re: [R-sig-Geo] NA color composite

2010-07-20 Thread Robert J. Hijmans
Carsten, Thanks for reporting this. I fixed the NA handling in plotRGB in raster 1.3-1, which should be available from R-Forge later today, and from CRAN later this week. In stead of first reading the data with rgdal, a more direct way to get at what you want with 'raster' is: d1 <- stack('Hasen

Re: [R-sig-Geo] Bug in rasterToPoints() ?

2010-07-20 Thread Robert J. Hijmans
Thanks Agus, I fixed this in raster version 1.3-0. Robert On Mon, Jul 19, 2010 at 11:46 PM, Agustin Lobo wrote: > > While > test <-rasterToPoints(Br025,fun=function(x){x>0},asSpatialPoints=T) > does not work, the following works: >> test <-rasterToPoints(Br025,fun=function(x){x>0},asSpatialPoints

Re: [R-sig-Geo] error in errorsarlm

2010-07-20 Thread Roger Bivand
On Tue, 20 Jul 2010, elaine kuo wrote: Hello, Thanks for the always concerns. answering the previous questions first. 1. distance threshold But 9 what, metres, degrees? Should your code be: cbind(datam$lat,datam$lon), as the normal ordering is eastings then northings, not the opposite?

[R-sig-Geo] NA color composite

2010-07-20 Thread carstennh
Dear list members, I'm having trouble creating a rgb color composite based one a multiple grid object. Both the SGDF2PCT function from rgdal package and the plotRGB function from raster package can not handle NA's. For my work it is necessary to include No Data values for masking the study ar

Re: [R-sig-Geo] rasterImage support in R > 2.11.0 (problem for Windows in SDI mode)

2010-07-20 Thread Renaud Lancelot
Thank you for this information. I was able to reproduce the problem even after the first dev.off(). > sessionInfo() R version 2.11.1 Patched (2010-06-08 r52227) i386-pc-mingw32 locale: [1] LC_COLLATE=French_France.1252 LC_CTYPE=French_France.1252 LC_MONETARY=French_France.1252 LC_NUMERIC=C [5] L

Re: [R-sig-Geo] categoric and dummy variables

2010-07-20 Thread Roger Bivand
On Tue, 20 Jul 2010, Pinar Aslantas Bostan wrote: Dear all, I am working with gwr function to predict precipitation distribution. I have two datasets; first one is consisted of observation points (225 meteorological stations), and the second one is DEM points used to obtain precipitation p

[R-sig-Geo] categoric and dummy variables

2010-07-20 Thread Pinar Aslantas Bostan
Dear all, I am working with gwr function to predict precipitation distribution. I have two datasets; first one is consisted of observation points (225 meteorological stations), and the second one is DEM points used to obtain precipitation predictions over it. I am using these functions: bw1=

Re: [R-sig-Geo] error in errorsarlm

2010-07-20 Thread elaine kuo
Hello, Thanks for the always concerns. answering the previous questions first. >> 1. distance threshold >> > > But 9 what, metres, degrees? Should your code be: > cbind(datam$lat,datam$lon), as the normal ordering is eastings then > northings, not the opposite? => unit is degree, and the la

[R-sig-Geo] rasterImage support in R > 2.11.0 (problem for Windows in SDI mode)

2010-07-20 Thread Michael Sumner
Hello, this is just a note for new raster plotting support - and only affects users of R for Windows in SDI mode (window application mode- which is not the default install option). It affects both 32- and 64-bit R on Windows. The behaviour is a bit weird and subtle, so this should only be of inte