cross_correlation.py

#!/usr/bin/env python
"""Script to compute cross correlation between DSHIP data and LIMRAD94 data to check for possible time shift
input: DSHIP data
output: print to console of time shifts detected on given days and with used version
author: Johannes Roettenbacher
"""

import sys
sys.path.append('/projekt1/remsens/work/jroettenbacher/Base/larda')
sys.path.append('.')
import datetime as dt
import logging
import functions_jr as jr

log = logging.getLogger('jr')
log.setLevel(logging.INFO)
log.addHandler(logging.StreamHandler())

versions = [1, 2]
date = dt.datetime(2020, 2, 15, 0, 0, 0)
seapath = jr.read_seapath(date)
seapath = jr.calc_heave_rate(seapath)
for version in versions:
    t_shift, shift, seapath = jr.calc_time_shift_limrad_seapath(seapath, version, plot_xcorr=True)
    print(f"Time shift calculated from {date:%Y-%m-%d}: {t_shift:.4f} with version {version}.")
    # 1.9362, 1.9363

date = dt.datetime(2020, 2, 16, 0, 0, 0)
seapath = jr.read_seapath(date)
seapath = jr.calc_heave_rate(seapath)
for version in versions:
    t_shift, shift, seapath = jr.calc_time_shift_limrad_seapath(seapath, version)
    print(f"Time shift calculated from {date:%Y-%m-%d}: {t_shift:.4f} with version {version}.")
    # 1.9353, 1.9353

begin_dt = dt.datetime(2020, 1, 25, 6, 0, 5)
end_dt = dt.datetime(2020, 1, 25, 18, 0, 0)
seapath = jr.read_seapath(begin_dt)
seapath = jr.calc_heave_rate(seapath)
for version in versions:
    t_shift, shift, seapath = jr.calc_time_shift_limrad_seapath(seapath, version, plot_xcorr=True,
                                                                begin_dt=begin_dt, end_dt=end_dt)
    print(f"Time shift calculated from {begin_dt:%Y-%m-%d}: {t_shift:.4f} with version {version}.")
    # 0, 1.671, depending on the time chosen