[R-sig-Geo] Signification of variance and error in "krige" and "krige.cv"

2007-10-01 Thread bertrand toupin
Hi everyone! I interpolated observed precipitation on a grid. I also kept variance (from krige) and error (from krige.cv). I have some questions about the signification of those quantities. I obtain some interesting results. For example, the variance from "krige" is around 24 and 40 (mm? s

[R-sig-Geo] Finding nearest point

2007-04-26 Thread bertrand toupin
Hi, I have 2 sets of data points. One is a grid, consisting of 87 points. The other is 300 meteorological stations. I'd like to chose 10 of the grid points and I'd like to have the nearest station for each of my 10 grid-points. Is there a function that does it? I have R 2.4.0, with gstat and

[R-sig-Geo] Surface Maps

2006-10-25 Thread bertrand toupin
Hi everyone, I'm trying to do a surface map like this one : http://img63.imageshack.us/my.php?image=terreeu4.jpg Of course, my domain is smaller :) My entry are like that (polar stereographic grid from a numerical model): latlongData 40-77.024.0 42 -77.223.0 41

Re: [R-sig-Geo] Map projection

2006-10-23 Thread bertrand toupin
ok, seems to be the path to success! Thanks again! Roger Bivand <[EMAIL PROTECTED]> wrote: On Mon, 23 Oct 2006, bertrand toupin wrote: > Hi, earlier, I used those functions to transform a lat-long database to a > Lambert Azimuthal Equal Area (LAEA) : > > Lam_proj <- p

[R-sig-Geo] Map projection

2006-10-23 Thread bertrand toupin
Hi, earlier, I used those functions to transform a lat-long database to a Lambert Azimuthal Equal Area (LAEA) : Lam_proj <- project(cbind(Grille$x, Grille$y),"+proj=laea +lat_0=latitudem_grid +lon_0=longitudem_grid +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m") Now that my calculations is don

Re: [R-sig-Geo] krige.cv, again :)

2006-10-10 Thread bertrand toupin
oh...ok! I'm indeed using it on a large area. Thanks! [EMAIL PROTECTED] wrote: > hi, me again :) > > I've read that in the help of krige.cv : > > Note: > > Leave-one-out cross validation seems to be much faster in plain > (stand-alone) gstat, apparently quite a bit of the effort is spe

[R-sig-Geo] krige.cv, again :)

2006-10-06 Thread bertrand toupin
hi, me again :) I've read that in the help of krige.cv : Note: Leave-one-out cross validation seems to be much faster in plain (stand-alone) gstat, apparently quite a bit of the effort is spent moving data around from R to gstat. Since I have 6 ensemble of 2300 stations, I was w

Re: [R-sig-Geo] Error with krige.cv

2006-10-02 Thread bertrand toupin
ovariance matrix; the most common is that you have duplicate observations (i.e. two or more observations on the same location). Use zerodist to find out. -- Edzer bertrand toupin wrote: > Hi everyone, another obstacle towards the achievement :) > > I'm trying to do a kriging of my

[R-sig-Geo] Error with krige.cv

2006-10-02 Thread bertrand toupin
Hi everyone, another obstacle towards the achievement :) I'm trying to do a kriging of my data on a regular grid. My data is organized like : Donnees (the data matrix) Period x(km) y(km) Stations z(m) 1 400 400 22.5 120 1 430

Re: [R-sig-Geo] RE : Mass conversions of Lat/Long to UTM (many zone)

2006-09-25 Thread bertrand toupin
Hi! It's me again. I've tried everything this morning and it seems to work perfectly. so a big thank you to eberyone that answered. However, I still have a last question regarding this matter. You said : "+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m with