DEM Regularised Inversion Calculation in Python (Hannah & Kontar 2012)
Project description
DEMREG-PY
This is a python implementation of Hannah & Kontar (2012)'s regularised inversion method. The code is tightly based on the IDL mapping version of the DEM reg-inv code found at https://github.com/ianan/demreg in addition, the code enforces a positivity constraint on the DEM (hence pos)
The philosophy was to produce as similar a piece of software as the original version and as such, this python version has been shown to recover the same DEM as the IDL version (to within approximately 4 significant figures). It is likely this philosophy has lead to performance hits and I plan to go back and address the more hacky parts of the code at a later date.
To calculate a DEM you first need:
Data and associated error for a range of channels: e.g. in dn/s/px or counts/s
Temperature dependent channel response: How sensitive are your channels to plasma of each temperature?
To use: simply call dn2dem_pos with either a single pixel, 1d slice or 2d map of DN values as a function of filter (and associated error on DN), an array of temperatures over which to perform the DEM analysis and a temperature response for those filters.
For large datasets, beyond 200 pixels, the code switches to a parallel execution providing significant speedups.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.