python 2d-gaussian fitting

#目的：对一个二维的data做二维高斯拟合

#需加载的包，另外还需要numpy等

from astropy.io import fits

from astropy.modeling import Fittable2DModel, Parameter, models, fitting

from astropy.wcs import WCS

import astropy.units as u

#读入XXXX.fits

fitsname = 'XXXX.fits'

FITS_1 = fits.open(fitsname)

data = FITS_1[0].data

hdr = FITS_1[0].header

wcs = WCS(m0_fits[0].header)

#二维高斯模型&输入参数初值和限制范围：

g_init = models.Gaussian2D(amplitude=amp_start, x_mean=x_mu_start, y_mean=y_mu_start, x_stddev=x_sig_start, y_stddev=y_sig_start)

g_init.amplitude.fixed=True #固定amp值

g_init.x_mean.bounds=[x_mu_start-5,x_mu_start+5] #限制x_mu变化范围，下同

g_init.y_mean.bounds=[y_mu_start-5,y_mu_start+5]

g_init.x_stddev.bounds=[0.1,5]

g_init.y_stddev.bounds=[0.1,5]

<