Add regularised visibility calculation
Opened this issue · 0 comments
samaloney commented
- Visibilities with an observer coordinate frame (d_sun)
- Reorder vis so for each (u,v) as a function of energy
- 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)
- Reject low energy bins if error too large (magic number)
- Need the SRM for vis energies (attenuator)
- get bremsstrahlung cross section
- estimate electron flux simple multiplication -> 1st estimate of electron vis
- SVD on matrix SRM * cross section * energy (cm^2 keV-1) [ for photon could do on SRM * energy (cm^2 keV-1)]
- regaulisation using SVD, with both discrepancy and 3-sigman -> reg electron vis
- confidence strip - perturb the count spectra repeat above to obtain an estimate of confidence