天文图像拟合程序imfit的python包装
pyimfit的Python项目详细描述
Pyimfit
这是天文图像拟合程序imfit的python包装器。
联机文档：https://pyimfit.readthedocs.io/en/latest/。
示例用法
下面假设有一个交互式python会话（例如ipython会话 或笔记本电脑）：
from astropy.io import fits
import pyimfit
imageFile = "<path-to-FITS-file-directory>/ic3478rss_256.fits"
imfitConfigFile = "<path-to-config-file-directory>/config_exponential_ic3478_256.dat"
# read in image data, convert to proper double-precisions, little-endian format
image_data = fits.getdata(imageFile)
# construct model from config file; construct new Imfit fitter based on model,;
model_desc = pyimfit.ModelDescription.load(configFile)
# create an Imfit object, using the previously loaded model configuration
imfit_fitter = pyimfit.Imfit(model_desc)
# load the image data and image characteristics and do a standard fit
# (using default chi^2 statistics and Levenberg-Marquardt solver)
imfit_fitter.fit(image_data, gain=4.725, read_noise=4.3, original_sky=130.14)
# check the fit and print the resulting best-fit parameter values
if imfit_fitter.fitConverged is True:
print("Fit converged: chi^2 = {0}, reduced chi^2 = {1}".format(imfit_fitter.fitStatistic,
imfit_fitter.reducedFitStatistic))
bestfit_params = imfit_fitter.getRawParameters()
print("Best-fit parameter values:", bestfit_params)
欢迎加入QQ群-->： 979659372
