Bonjour,
J'ai fusionné vos deux messages (qui traitent du même sujet) et je passe en anglais afin de permettre à la communauté d'utilisateurs de nous comprendre.
For obtaining a smoothing estimation map by kriging you should first check your moving neighborhood parameters.
Try with a unique neighborhood just for test:
- Code: Select all
neigh = neigh.create(ndim = 2, type = 0)
If you have too much data samples, you should use an ellipsoid neighborhood (type = 2) with an ellipse size big enough and few sectors.
Could you plot your data location over your estimation map please ?
- Code: Select all
plot(data.db,add=T)
The error regarding the matrix inversion is usually due to the presence of duplicate samples (two samples at same location) in the data set.
Please, execute the function "duplicate" before calculating your experimental variogram and its model and launching kriging.
Hope this helps.
Regards