astropy / astropy/specutils

Using specutils to code a 3D spectra Time-series

Open
#1,055 0 comments 0 reactions 0 assignees View on GitHub
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

Open the contributing 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

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.