Add parser for TESS DVT files

CHANGES.rst

@@ -1,6 +1,8 @@
 1.1.0 (unreleased)
 ------------------
 
+* Add support for loading TESS DVT files. [#164]
+
 
 1.0.0 (12-02-2024)
 ------------------

lcviz/helper.py

@@ -1,3 +1,4 @@
+from astropy.io.fits import getheader
 import astropy.units as u
 import ipyvue
 import os
@@ -144,7 +145,7 @@
             # already been initialized
             plugin._obj.vdocs = self.app.vdocs
 
-    def load_data(self, data, data_label=None):
+    def load_data(self, data, data_label=None, extname=None):
         """
         Load data into LCviz.
 
@@ -161,7 +162,21 @@
             automatically determined from filename or randomly generated.
             The final label shown in LCviz may have additional information
             appended for clarity.
+        extname : str or `None`
+            Used for DVT parsing if only a single TCE from a multi-TCE file should be
+            loaded. Formatted as 'TCE_1', 'TCE_2', etc.
         """
+        # Determine if we're loading a DVT file, which has a separate parser
+        if isinstance(data, str):
+            header = getheader(data)
+            if (header['TELESCOP'] == 'TESS' and 'CREATOR' in header and
+                    'DvTimeSeriesExporter' in header['CREATOR']):
+                super().load_data(data=data,
+                                  parser_reference='tess_dvt_parser',
+                                  data_label=data_label,
+                                  extname=extname)
+                return
+
         super().load_data(
             data=data,
             parser_reference='light_curve_parser',

lcviz/parsers.py

@@ -1,4 +1,6 @@
 import os
+
+from astropy.io import fits
 from glue.config import data_translator
 from jdaviz.core.registries import data_parser_registry
 import lightkurve
@@ -12,9 +14,46 @@
     'kepler': {'prefix': 'Q', 'card': 'QUARTER'},
     'k2': {'prefix': 'C', 'card': 'CAMPAIGN'},
     'tess': {'prefix': 'S', 'card': 'SECTOR'},
+    'tess dvt': {'prefix': '', 'card': 'EXTNAME'}
 }
 
 
+@data_parser_registry("tess_dvt_parser")
+def tess_dvt_parser(app, file_obj, data_label=None, show_in_viewer=True, **kwargs):
+    '''
+    Read a TESS DVT file and create a lightkurve object
+    '''
+    hdulist = fits.open(file_obj)
+    ephem_plugin = app._jdaviz_helper.plugins.get('Ephemeris', None)
+    extname = kwargs.pop('extname')
+
+    # Loop through the TCEs in the file. If we only want one (specified by
+    # `extname` keyword) then only load that one into the viewers and ephemeris.
+    for i in range(1, len(hdulist)-1):
+        data = hdulist[i].data
+        header = hdulist[i].header
+        lc = lightkurve.LightCurve(data=data,
+                                   flux=data['LC_INIT'],
+                                   flux_err=data['LC_INIT_ERR'])
+        lc.meta = hdulist[0].header
+        lc.meta['MISSION'] = 'TESS DVT'
+        lc.meta['FLUX_ORIGIN'] = "LC_INIT"
+        lc.meta['EXTNAME'] = header['EXTNAME']
+
+        if extname is not None and header['EXTNAME'] != extname:
+            show_ext_in_viewer = False
+        else:
+            show_ext_in_viewer = show_in_viewer
+
+        light_curve_parser(app, lc, data_label=data_label,
+                           show_in_viewer=show_ext_in_viewer, **kwargs)
+
+        # add ephemeris information from the DVT extension
+        if ephem_plugin is not None and show_ext_in_viewer:
+            ephem_plugin.period = header['TPERIOD']
+            ephem_plugin.t0 = header['TEPOCH']
+
+
 @data_parser_registry("light_curve_parser")
 def light_curve_parser(app, file_obj, data_label=None, show_in_viewer=True, **kwargs):
     # load a LightCurve or TargetPixelFile object:
@@ -29,7 +68,6 @@
         if data_label is None:
             data_label = os.path.splitext(os.path.basename(file_obj))[0]
 
-        # read the light curve:
         light_curve = lightkurve.read(file_obj)
 
     elif isinstance(file_obj, cls_with_translator):