[R-sig-Geo] Local singularity

[Edzer Pebesma]

José,

you're sending this message to the wrong list, as it has nothing to do
with R -- it seems you are using gstat stand-alone, a program that
predates the gstat R package. The more appropriate list would be

http://list.52north.org/mailman/listinfo/geostatistics

With best regards,

On 12/07/2011 09:29 PM, José Matos wrote:
> Dear All,
> 
> I'm testing different interpolation settings and methods for a rainfall
> field. From the complete field I take some points, perform the
> interpolation and then compute the errors. The variogram was computed
> for the original data.
> While trying local neighborhoods, I got some unexpected results and am
> fairly confident that I'm doing something wrong, although I not sure
> about what it is.
> 
> It is as if outside the predetermined radius, gstat is computing the
> estimate gathering points up to the predefined minimum but without the
> closest point.
> 
> One of the command sets I used (made to enhance this effect) was the
> following:
> data(TRMM): 'dataKrig', x=1, y=2, v=3, radius=2, min=4, max=4, force;
> predictions(TRMM): 'predKrig';
> mask: 'maskKrig.asc';
> variogram(TRMM):   0.00000401 Nug(0) + 1232.26185112 Pow(0.54167,0,1);
> 
> A more elaborate description and all the files used to run gstat can be
> found at:
> https://documents.epfl.ch/users/g/ga/gamito/www/gstat/gstatHelp.zip
> 
> All help on this issue is much appreciated.
> 
> Thank you in advance.
> 
> Kind regards,
> José Matos
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de