Using specutils to code a 3D spectra Time-series
- Dominant language
- Python
- Stars
- 204
- Forks
- 135
- Avg merge
- 12h 55m
- Merged PRs (30d)
- 1
Description
I have several FITS files that I have downloaded from DACE CCF spectra
(https://dace.unige.ch/sun/?)
I would like to create a 3d-spectra-timeseries-with-python
This is an example of one of the filenames on the DACE database
r.HARPN.2018-01-18T14:20:10.645_CCF_A.fits
XTENSION= 'IMAGE ' / IMAGE extension
BITPIX = -32 / number of bits per data pixel
NAXIS = 2 / number of data axes
NAXIS1 = 49 / length of data axis 1
NAXIS2 = 70 / length of data axis 2
PCOUNT = 0 / required keyword; must = 0
GCOUNT = 1 / required keyword; must = 1
EXTNAME = 'SCIDATA '
CHECKSUM= 'T9gHV8fFT8fFT8fF' / HDU checksum updated 2020-08-25T17:28:02
DATASUM = '2349888738' / data unit checksum updated 2020-08-25T17:28:02
I have attempted to follow the coding on the following site
(https://stellartrip.net/3d-spectra-timeseries-with-python/?lang=en)
However I am getting an errors, is there a standard code that I can use to add several FITS files to generate a MESHplot.
KeyError: "Keyword 'CDELT1' not found."
More info on the file
SIMPLE = T / file does conform to FITS standard
BITPIX = 8 / number of bits per data pixel
NAXIS = 0 / number of data axes
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
DATE = '2020-08-25T17:28:02' / file creation date (YYYY-MM-DDThh:mm:ss UT)
VERSION = 120203 / Application parameter value
DWELL = 300000 / Application parameter value
RD_TIME = 0 / Application parameter value
RS_TIME = 383 / Application parameter value
FRM_TYPE= 1 / Application parameter value
NUM_EXPS= 1 / Application parameter value
EL_GAIN = 1 / Application parameter value
SPEED = 1 / Application parameter value
SER_REG = 0 / Application parameter value
CLK_MON1= 0 / Application parameter value
CLK_MON2= 1 / Application parameter value
TELESCOP= 'TNG '
TIMESYS = 'UTC ' / Explicit time scale specification
DATE-OBS= '2018-01-18T14:20:10.645' / Date and time of start of obs. in UTC
PROCESS = 'Single raw read' / Data processing mode
FILENAME= 'HARPN.2018-01-18T14-20-10.645.fits' / FITS file generatedName
TIME = '14:25:31' / Time FITS file was written
Contributor guide
Research direction
Start by reproducing the CDELT1 KeyError with the DACE HARPN FITS file and inspect the primary header and SCIDATA extension described in the issue. Compare the file structure with the linked 3D spectra time-series example and the relevant specutils entry points. Done means establishing whether the file is supported and documenting or implementing a clear path to load multiple files and create the requested mesh plot.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- data
- Issue type
- Bug
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 20/100