I am wondering how to decompose the data[x,y] into several elliptic-shaped functions, i.e. H/(1+((x-x1)/a1)^2+((y-y1)/b1)^2)^(3/2).
I could follow in 1D case, decomposition data[x] into H/(1+((x-x1)/a1)^2), by just following (with modification) https://mathematica.stackexchange.com/questions/26336/how-to-perform-a-multi-peak-fitting but I cannot perform in 2D case.
Does anyone have a good idea?
Thank you.