Skip to content

Updated process_lislipG.py to handle files with missing data records #215

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
sflynn-itsc merged 1 commit into from
Mar 12, 2026
Merged

Updated process_lislipG.py to handle files with missing data records #215

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
39 changes: 28 additions & 11 deletions mdx/granule_metadata_extractor/processing/process_lislipG.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -36,9 +36,15 @@ def get_variables_min_max(self):
def get_nc_metadata(self):
datafile = Dataset(self.file_path)
ref_time = datetime(1993,1,1) #TAI93 time
utc = datafile['bg_info_TAI93_time'][:]
minTime, maxTime = [ref_time+timedelta(seconds=utc.min()),
ref_time+timedelta(seconds=utc.max())]
try:
utc = datafile['bg_info_TAI93_time'][:]
minTime, maxTime = [ref_time+timedelta(seconds=utc.min()),
ref_time+timedelta(seconds=utc.max())]
except: #if bg_info_TAI93_time unavailable, use bg_orbit_summary_TAI93_xxxxx
timevar0 = datafile.variables['bg_orbit_summary_TAI93_start'][:]
timevar1 = datafile.variables['bg_orbit_summary_TAI93_end'][:]
minTime, maxTime = [ref_time + timedelta(seconds=timevar0.item()),
ref_time + timedelta(seconds=timevar1.item())]

#pre-define bounding box same as on-prem bounding box hard-coded by Nathan
maxlat, minlat, maxlon, minlon = [35., -35., 180., -180.]
Expand All @@ -61,16 +67,27 @@ def get_nc_metadata(self):
def get_hdf_metadata(self):
f = HDF(self.file_path)
vs = f.vstart()
data = vs.attach('bg_info')
bg_info = data[:]#['TAI93_time', 'alert_summary',......,'bg_value', 'noise_index', 'event_count']

bg_info_TAI93_time = [bg_info[i][0] for i in range(len(bg_info))]
utc = np.array(bg_info_TAI93_time)

ref_time = datetime(1993,1,1) #TAI93 time
minTime, maxTime = [ref_time+timedelta(seconds=utc.min()),
ref_time+timedelta(seconds=utc.max())]
data.detach()
data = vs.attach('bg_info')
if data._nrecs > 0:
bg_info = data[:]#['TAI93_time', 'alert_summary',......,'bg_value', 'noise_index', 'event_count']
bg_info_TAI93_time = [bg_info[i][0] for i in range(len(bg_info))]
utc = np.array(bg_info_TAI93_time)
minTime, maxTime = [ref_time+timedelta(seconds=utc.min()),
ref_time+timedelta(seconds=utc.max())]
data.detach() #detach 'bg_info'
else: #if 'bg_info' vdata has no records
data.detach() #detach 'bg_info'
vd = vs.attach('bg_orbit_summary')
#extract time info from 'bg_orbit_summary'
#'bg_orbit_summary' has one record
#data fields: ['id_number', 'TAI93_start', 'UTC_start', 'GPS_start', 'TAI93_end',...]
recs = vd[0]
minTime, maxTime = [ref_time + timedelta(seconds=recs[1]),
ref_time + timedelta(seconds=recs[4])]
vd.detach() #detach 'bg_orbit_summary'


#pre-define bounding box same as on-prem bounding box hard-coded by Nathan
maxlat, minlat, maxlon, minlon = [35., -35., 180., -180.]
Expand Down
Loading