More Checks Against Data-loss in Correlator

* Fixed a bug that was causing slight data-loss during unstacked cross-correlation, in the way blocks of data were being fetched * Added more tests to capture heretofore untested scenarios
GeoscienceAustralia · Apr 6, 2023 · f47b069 · f47b069
1 parent 5bd2986
commit f47b069
Show file tree

Hide file tree

Showing 5 changed files with 19,209 additions and 8,427 deletions.
diff --git a/seismic/xcorqc/Readme.md b/seismic/xcorqc/Readme.md
@@ -49,7 +49,7 @@ The input parameters above dictate the following:
 ## Default Mode
 In the default mode (without interval-stacking),  
 
-DRS = WS * (1 - OLAP) * RAW + WS * WBL * 2 + OLAP * WS
+DRS = WS * (1 - OLAP) * RAW + WS * WBL * 2 + OLAP * WS * 2
 
 A table showing example parameterizations and their implications is as follows:
 

diff --git a/seismic/xcorqc/correlator.py b/seismic/xcorqc/correlator.py
@@ -577,7 +577,7 @@ def main(data_source1, data_source2, output_path, window_seconds, window_overlap
 
         read_ahead_window_seconds = window_seconds * (1 - window_overlap) * read_ahead_windows + \
                                     window_seconds * window_buffer_length * 2 + \
-                                    window_overlap * window_seconds
+                                    window_overlap * window_seconds * 2
         interval_seconds = read_ahead_window_seconds
         apply_stacking = False
         #print(read_ahead_window_seconds)