В настоящее время я работаю с астрономическими данными, среди которых у меня есть изображения комет. Я хотел бы удалить фоновый градиент неба в этих изображениях из-за времени захвата (сумерки). Первая программа, которую я разработал для этого, взяла выбранные пользователем точки из Matplotlib "ginput" (x, y), вытащила данные для каждой координаты (z), а затем зациклировала данные в новом массиве с помощью ScidPy "griddata".
Поскольку предполагается, что фон изменяется незначительно, я хотел бы поместить трехмерный полином нижнего порядка в этот набор точек (x, y, z). Однако "griddata" не учитывает порядок ввода:
griddata(points,values, (dimension_x,dimension_y), method='nearest/linear/cubic')
Любые идеи о другой функции, которая может быть использована, или метод разработки набора квадратов leas-квадратов, который позволит мне контролировать порядок?