From f905f2cd3fad62ca32c93b811f1de4d54f179e46 Mon Sep 17 00:00:00 2001
From: arigdon-usgs <arigdon@usgs.gov>
Date: Fri, 20 Jul 2018 11:29:15 -0600
Subject: [PATCH] just putting some things into place for testing
---
geomagio/TimeseriesUtility.py | 19 +++++++++++++++++++
1 file changed, 19 insertions(+)
diff --git a/geomagio/TimeseriesUtility.py b/geomagio/TimeseriesUtility.py
index da7086d9..aea308c5 100644
--- a/geomagio/TimeseriesUtility.py
+++ b/geomagio/TimeseriesUtility.py
@@ -195,17 +195,36 @@ def merge_streams(*streams):
stream with contiguous traces merged, and gaps filled with numpy.nan
"""
merged = obspy.core.Stream()
+ # masked_trace = None
+
# add unmasked, split traces to be merged
for stream in streams:
merged += mask_stream(stream)
+
+ # if trace is completely masked separate out to be added back in later
+ # for trace in merged:
+ # if trace.data.mask.all():
+ # if not masked_trace:
+ # masked_trace = trace
+ # else:
+ # masked_trace += trace
+
# split traces that contain gaps
merged = merged.split()
+
# merge data
merged.merge(
# 1 = do not interpolate
interpolation_samples=1,
# 1 = when there is overlap, use data from trace with last endtime
method=1)
+
+ # trim masked trace to the same size as other traces and add back to merged stream
+ # if masked_trace:
+ # masked_trace.trim(merged[0].stats.starttime,merged[0].stats.endtime)
+ # merged += masked_trace
+
+
# convert back to NaN filled array
merged = unmask_stream(merged)
return merged
--
GitLab