[R-sig-Geo] spatiotemporal co-kriging?

2016-11-08 Thread Santiago Beguería Portugués
Dear R-sig-geo list members, I am working on a spatio-temporal kriging model (using gstat::krigeST) of two closely related variables. I’m getting good results on one of them, and not so good on the other. Since both dependent variables are inter-related, I expect that co-kringing would be a

Re: [R-sig-Geo] Kriging with uncertain data

2016-10-11 Thread Santiago Beguería Portugués
Hi, Thank you for your suggestions. I gave the gstat with weights argument a try, but I’m not sure about the results / my implementation. I have worked out an example with the meuse dataset so I can share it with you: library(gstat) library(sp) data(meuse) coordinates(meuse) <- ~x+y

[R-sig-Geo] rgdal - problem with reading grib files

2012-01-24 Thread Santiago Beguería Portugués
Hi, I am stuck when trying to read data from certain grib files with GDAL. When I try to read the data I get the following error: tmp - readGDAL('./data/ghcn_cams/ghcn_cams_1948_cur.grb') ... Error in `proj4string-`(`*tmp*`, value = S4 object of class CRS) : Geographical CRS given to