Comments (3)
I am happy to say that this longstanding issue is now solved. Thank you for reporting it again, otherwise I would have delayed the investigation further.
There was a very subtle bug in the theoretical variogram expressions where the discontinuity jump at the origin was not being applied. I fixed all the variogram models accordingly as shown here: juliohm/Variography.jl@84b540c
I have tested with your data set again, and the solution looks correct now. As a general rule, when you find ill-conditioned systems in the future, the solution is to add a nugget effect as you did above. This is equivalent to adding a positive term to the diagonal of the covariance which improves the condition number.
The new release of Variography.jl should come out soon and you can get it by updating with the Julia package manager.
Marking as solved. 👍
from geostats.jl.
Thank you @tianzeliu for reporting the issue. I encountered similar problems in the past (#12) and delayed the resolution for too long. The problem with the Gaussian variogram is known in the literature. This variogram has a behavior at the origin (it is very flat) that leads to quite ill-conditioned covariance matrices. Other software such as SGeMS and GSLIB add a small positive term to the diagonal of the Kriging system to improve the condition number, and I think I will do the same here in GeoStats.jl
I was trying to avoid the solution described above because this type of preconditioning inhibits the interpolation property of Kriging. But because the implementation in GeoStats.jl never touches locations with data, this shouldn't be an issue. I will use your data set to debug the problem further, thank you.
from geostats.jl.
@tianzeliu could you please provide a minimum example so that I can reproduce the exact same plot?
Thank you for reporting the issue again.
from geostats.jl.
Related Issues (20)
- Fix `show` of very large `GeoTable` HOT 1
- Observable transform pipelines HOT 1
- Add support for `geotable[bools, cols]` HOT 1
- Allow other Kriging variants in sequential Gaussian simulation
- Online generation of large realizations
- Gaussian Processes with Variogram vs Covariance HOT 4
- pairwise! with Covariance types
- Add `Aggregate` transform HOT 1
- missing "(" , ")" in 8.4 Applying expressions HOT 3
- Add `join` for non-geospatial join HOT 1
- Review colorbar creation with `viz` and `viewer` HOT 1
- Rename `colorscheme` --> `colormap` in `viz` and `viewer` HOT 1
- Add new `RegularRefinement` method HOT 1
- Add better error message in `viz` and `viewer` HOT 1
- Add `InterpolateMissing` transform
- Refactor `showfacets` and `facetcolor` in `viz` HOT 1
- Visualization of 3D `Ellipsoid` geometry HOT 1
- Add `DropNaN` and `InterpolateNaN` HOT 1
- Allow `georef(x, ::Symbol)` HOT 6
- Idea / Question: Kriging under constraints HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from geostats.jl.