search

TCDSolar / stixpy

STIX data analysis in python
https://stixpy.readthedocs.io/
BSD 3-Clause "New" or "Revised" License
17 stars 20 forks source link

Add regularised visibility calculation #103

Open samaloney opened 5 months ago

samaloney commented 5 months ago
  1. Visibilities with an observer coordinate frame (d_sun)
  2. Reorder vis so for each (u,v) as a function of energy
  3. Rescale spectra by fitting a single, broken or sum power laws pick by lowest chi-square (divide by fit to have similar magnitude number to make inversion possible/faster)
  4. Reject low energy bins if error too large (magic number)
  5. Need the SRM for vis energies (attenuator)
  6. get bremsstrahlung cross section
  7. estimate electron flux simple multiplication -> 1st estimate of electron vis
  8. SVD on matrix SRM cross section energy (cm^2 keV-1) [ for photon could do on SRM * energy (cm^2 keV-1)]
  9. regaulisation using SVD, with both discrepancy and 3-sigman -> reg electron vis
  10. confidence strip - perturb the count spectra repeat above to obtain an estimate of confidence