import numpy as np import textkernel import tabulation import cspyce ROOT = '/Users/mark/GitHub/pds-migration/COUVIS/' LSK = '/Volumes/Data-HD/Resources/SPICE/leapseconds.ker' EOM = '2017-09-16/00:00:00' cspyce.furnsh(LSK) tk = textkernel.from_file(ROOT + 'cas00172.tsc') coeffts = np.array(tk['SCLK'][1]['COEFFICIENTS_82']).reshape(-1, 3) sclk = (coeffts[:, 0] + tk['SCLK_PARTITION_START_82']) / 256. tdt = coeffts[:, 1] tdb = [cspyce.unitim(t, 'TDT', 'TDB') for t in tdt] tdb_end = cspyce.utc2et(EOM) tdt_end = cspyce.unitim(tdb_end, 'TDB', 'TDT') rate_end = coeffts[-1, 2] sclk_end = sclk[-1] + (tdt_end - tdt[-1]) / rate_end sclk = np.array(list(sclk) + [sclk_end]) tdb = np.array(list(tdb) + [tdb_end]) cas00172 = tabulation.Tabulation(tdb, sclk) tk = textkernel.from_file(ROOT + 'cas00171.tsc') coeffts = np.array(tk['SCLK'][1]['COEFFICIENTS_82']).reshape(-1, 3) sclk = (coeffts[:, 0] + tk['SCLK_PARTITION_START_82']) / 256. tdt = coeffts[:, 1] tdb = [cspyce.unitim(t, 'TDT', 'TDB') for t in tdt]
# Syntax: # python coiss_production.py COISS/COISS_2xxx/index/index.lbl ################################################################################ import oops import oops.inst.cassini.iss as iss import numpy as np import os, sys, traceback import cspyce import metadata as meta iss.initialize() PREFIX = '/Users/mark/Desktop/CK-cruise/' cspyce.furnsh(PREFIX + '97288_98002rc.bc') cspyce.furnsh(PREFIX + '98002_98091rc.bc') cspyce.furnsh(PREFIX + '98091_98185rc.bc') cspyce.furnsh(PREFIX + '98185_98279rc.bc') cspyce.furnsh(PREFIX + '98279_99001rc.bc') cspyce.furnsh(PREFIX + '99001_99093rc.bc') cspyce.furnsh(PREFIX + '99093_99182rc.bc') cspyce.furnsh(PREFIX + '99182_99274rc.bc') cspyce.furnsh(PREFIX + '99274_00001rc.bc') cspyce.furnsh(PREFIX + '00001_00092rc.bc') cspyce.furnsh(PREFIX + '00092_00183rc.bc') cspyce.furnsh(PREFIX + '00183_00275rc.bc') cspyce.furnsh(PREFIX + '00275_01001rc.bc') cspyce.furnsh(PREFIX + '01001_01091rc.bc') cspyce.furnsh(PREFIX + '01091_01182rc.bc') cspyce.furnsh(PREFIX + '01182_01274rc.bc')
# Syntax:
# NHxxLO_pluto_metadata.py path/to/NHPELO_2xxx_supplemental_index.lbl
################################################################################
import oops
import oops.inst.nh.lorri as lorri
import numpy as np
import os, sys, traceback
import metadata as meta
PLANET = 'PLUTO'
BIGMOON = 'CHARON'
import cspyce
cspyce.furnsh('/Resources/SPICE/General/PCK/pck00010_edit_v01.tpc')
# This reload last overrides some bug that produces this error message:
# SPICE(INSUFFICIENTANGLES) -- TISBOD: Insufficient number of nutation/precession angles for body 599 at time 2.8472754946623E+07.
############################################
# Key parameters of run
############################################
execfile('COLUMNS_PLUTO.py') # define the columns, tiles, etc.
SAMPLING = 8 # pixel sampling density
SELECTION = 'S' # summary files only
############################################
# Construct the meshgrid for each LORRI FOV
############################################