This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
data:data_analysis_manual [2024/11/11 20:17] Peter Young [Line fitting] |
data:data_analysis_manual [2025/01/09 11:25] (current) eric buchlin Moved Python line fitting |
||
---|---|---|---|
Line 263: | Line 263: | ||
==== Line fitting ==== | ==== Line fitting ==== | ||
+ | For IDL Gaussian fitting software, please check: **[[data:line_fitting_idl|IDL software for emission line fitting]]**. | ||
- | === Fitting a Gaussian at each spatial pixel in a raster (IDL) === | + | Fitting in Python can be done with any fitting library, including [[https://docs.astropy.org/en/stable/modeling/|astropy.modelling]]. From astropy 7.0, there is a user-friendly way of doing the fitting in parallel ([[https://github.com/aperiosoftware/parallel-modeling-examples/blob/main/spice/Final_Tutorial.ipynb|tutorial]] by Stuart Mumford). |
- | Spectra at each spatial pixel in a raster can be fit with Gaussians using the eis\_auto\_fit IDL routine developed for the Hinode/EIS mission (you must have the EIS branch in your Solarsoft distribution). The key step is to format a SPICE wavelength window into a "windata" structure using the routine spice\_getwindata. An example call is: | ||
- | IDL> wd=spice_getwindata(file) ; choose a wavelength window | ||
- | IDL> eis_fit_template,wd,template ; create a template for the fit | ||
- | IDL> eis_auto_fit,wd,fit,template=template | ||
- | IDL> eis_fit_viewer,wd,fit | ||
- | |||
- | The routine eis\_fit\_viewer is a GUI-based routine for browsing the fit results. More details on the eis\_auto\_fit suite of routines are available in [[https://doi.org/10.5281/zenodo.6339584 | ||
- | |EIS Software Note 15]]. | ||
- | |||
- | |||
- | The EIS software enable a number of operations to be performed on windata structures, as described in [[https://sohoftp.nascom.nasa.gov/solarsoft/hinode/eis/doc/eis_notes/21_WINDATA/eis_swnote_21.pdf|EIS Software Note 21]]. For example, eis\_bin\_windata allows spatial binning of the rasters. | ||
- | |||
- | If you have problems running the EIS software on SPICE data, please contact [[https://science.gsfc.nasa.gov/sed/bio/peter.r.young|Dr. Peter Young]]. | ||
=== Line window width and spectral tilt === | === Line window width and spectral tilt === |