Merge overlapping grids

[ckohnke]

Description of the desired feature:

Grid stitching methods to handle grid discrepancies at overlapping grid points. This may already be possible with the tools in Verde and an example may be beneficial if it is.

Using airborne magnetic data as an example:

• After processing, leveling, and gridding individual flights of a magnetic survey there is usually a difference between grids in the region where they overlap.
• Theoretically, these overlapping points are at the same time (diurnal removed), and space (equivalent source, upward continuation, etc.), and should have the same value. However, a variety of factors, such as instrument drift, can cause mismatch.
• If one is lucky, there is only a static shift between the grids in the overlapping region and can be handled by using a common datum between the datasets.
• However, if the difference in the data in the overlapping regions are more complex, stitching the grids together using the overlapping points becomes more nebulous and can easily lead to errors in the resulting final grid.

Are you willing to help implement and maintain this feature?

Potentially. I've been having trouble finding mathematical methods and processes to handle such situations.

[leouieda]

@ckohnke thanks for opening the issue! I have a few thoughts on this:

1. For gravity and magnetics, I really don't like these stitching approaches because they can heavily bias the long wavelengths. For those, I'd rather go back to the original data and use a single equivalent source model to re-grid them all.
2. You're right about there not being a good way to do this. For grids of other data types (topography, etc), maybe we could try to blend them in some way that only introduces distortions in the boundaries and doesn't propagate them to the whole grid. Or maybe warn users if there are a lot of discrepancies.
3. GMT's grdblend does something like this already. It's not yet available in PyGMT. We could try to implement something similar to what it does using xarray. It doesn't seem like it would be too difficult to do.

[leouieda]

It would also be worth looking into what xarray has for merging datasets. Perhaps that could be enough and we don't have to implement it ourselves.