Lu Li

Just find this paper, that might give an idea of how to include different density contrast (maybe different geology) into the inversion. Haas et al 2020. https://academic.oup.com/gji/article/221/3/1896/5808816?login=true

There is an example of how this code works. https://nbviewer.org/github/LL-Geo/PFToolbox/blob/master/FFT_Filter.ipynb

I'll work on the documentation (formula) & formate later... Any suggestions are welcome 😄

Opps.. commit to wrong branch 🙈

> > Opps.. commit to wrong branch see_no_evil > > Ha! No big deal, this can happen to everyone! But since `git` is so good, you can easily solve this...

I have added all filters in the new format. @santisoler Just curious how do we do the test? Do we want to do it by comparing the FFT result with...

We might need to make the magnetic point source work first then we could test RTP and pseudo_gravity... And how could we test the Gaussian filter?..

> > > transform is wrong in Y? > > > > > > @RichardScottOZ I think there's a missing `return da_filter` in the `pseudo_gravity_kernel` function: https://github.com/LL-Geo/harmonica/blob/5416778c88a5591030c765c203aa5adc54460c7b/harmonica/filters/_filters.py#L437-L528 > > Let...

> Thanks! > > so what is the difference between https://github.com/LL-Geo/PFToolbox/blob/master/FFT_Filter.ipynb and the changes here? > > Don't think I see it at first or second look. Yup. These are...

>  [gAerogmagnetico_Reducido_al_polo.zip](https://github.com/fatiando/harmonica/files/9662238/gAerogmagnetico_Reducido_al_polo.zip) Yup, I see the problem now. I think the issue come from xrft, which will return an increased order of y coordinate when doing ifft.  ...