Kriging wiht anisotropy modeling using great circle distance in gstat [SEC=UNCLASSIFIED]

Thank you Johan for the clarification.
Given this limitation with gstat, I was wondering if there is any package in R or so that supports kriging with anisotropy modelling using great circle distance. Or is gstat going to be developed to overcome this limitation? Thank you in advance!
Kind regards,
Jin


-----Original Message-----
From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Johan Van de Wauw
Sent: Saturday, 20 August 2011 12:49 AM
To: r-sig-geo at r-project.org
Subject: Re: [R-sig-Geo] Kriging wiht anisotropy modeling using great circle distance in gstat [SEC=UNCLASSIFIED]

Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, ?:
?gstat: value not allowed for: for long/lat data, anisotropy cannot be defined

Is there anything wrong in my code or it is a bug? Any suggestions are appreciated!

As you can read from the error message, nothing is wrong with your
code or with gstat, but anisotropy is not supported when you use great
circle distances.

Sincerely Yours
Johan

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo