I am looking for an IDL function which enables a 2-dimensional fit to gridded data, i.e. similar to the GAUSS2DFIT or the SFIT function. In the particular case the data is a grid array, let's say a=dblarr(x,y), where x, y and z=a(x,y) are the coordinates in the 3d space. The function I want to fit to my grid array could be a hemisphere, cone, or even any other mathematical area function, which can be described as z=f(x,y). CURVEFIT is something which is very to close to that what I need, but CURVEFIT is limited to vectors f(x), and cannot fit arrays f(x,y). How can I solve the problem?
|