A python-3 implementation of cone model Xie et al.,2004 for Solar Coronal Mass Ejection.
The file Xie's_cone_v1.py provides a python programme to get necessary CME parameters as discussed in the model. Folder Data consist of the datafiles used in the programme, while folder img consist of the output iamges.
- Dowload processed images of SOHO/LASCO observation from Helioviewer.org or SunPy.
- Read those data files using (data,header) variables with sequential number (In order to get running difference image; see line number 39 to 54).
- Select particular data file to infer CME parameters at particular time-stamp. (See line number 60 to 64)
- Compile the code.
- After compilation a popup window will appear in which difference image of CME is shown. Select 30 to 35 points on CME leading edge and press Esc key. (To fit ellipse we require only 5 points; but we select as many as points to minimize variation in parameters obtained from sequential frames)
- Programme will return CME paramters with fitted halo on top of SOHO/LASCO image.
A more detailed description of cone model, this python implementation and its validation is given in my M.Sc. thesis.
Currently this python program is in early stage. It only works with the SOHO/LASCO coronagraph data. So if you find any bug or have some suggestions for improvement, please notify me with Github Issues or Pull requests.
if you use this program for calculation and scientific publication, please site it using the DOI generated by zenodo. I would also appreciate it if you drop me an email on ronish.m22@gmail.com and tell me what you used it for.