From edea69d4b6216f4ac729eba78f12547807a2751a Mon Sep 17 00:00:00 2001 From: Cheryl Date: Wed, 20 May 2020 13:33:50 +0800 Subject: [PATCH] added buffer margin to orbit times in get_orbit --- contrib/stack/topsStack/Stack.py | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/contrib/stack/topsStack/Stack.py b/contrib/stack/topsStack/Stack.py index aa4105c..8290ab8 100644 --- a/contrib/stack/topsStack/Stack.py +++ b/contrib/stack/topsStack/Stack.py @@ -912,7 +912,8 @@ class sentinelSLC(object): self.SNWE=[min(lats),max(lats),min(lons),max(lons)] - def get_orbit(self, orbitDir, workDir): + def get_orbit(self, orbitDir, workDir, margin=60.0): + margin = datetime.timedelta(seconds=margin) datefmt = "%Y%m%dT%H%M%S" orbit_files = glob.glob(os.path.join(orbitDir, self.platform + '*.EOF')) if len(orbit_files) == 0: @@ -922,10 +923,10 @@ class sentinelSLC(object): for orbit in orbit_files: orbit = os.path.basename(orbit) fields = orbit.split('_') - orbit_start_date_time = datetime.datetime.strptime(fields[6].replace('V',''), datefmt) - orbit_stop_date_time = datetime.datetime.strptime(fields[7].replace('.EOF',''), datefmt) + orbit_start_date_time = datetime.datetime.strptime(fields[6].replace('V',''), datefmt) + margin + orbit_stop_date_time = datetime.datetime.strptime(fields[7].replace('.EOF',''), datefmt) - margin - if self.start_date_time > orbit_start_date_time and self.stop_date_time < orbit_stop_date_time: + if self.start_date_time >= orbit_start_date_time and self.stop_date_time < orbit_stop_date_time: self.orbit = os.path.join(orbitDir,orbit) self.orbitType = 'precise' match = True @@ -944,7 +945,7 @@ class sentinelSLC(object): fields = os.path.basename(orbitFile).split('_') orbit_start_date_time = datetime.datetime.strptime(fields[6].replace('V',''), datefmt) orbit_stop_date_time = datetime.datetime.strptime(fields[7].replace('.EOF',''), datefmt) - if self.start_date_time > orbit_start_date_time and self.stop_date_time < orbit_stop_date_time: + if self.start_date_time >= orbit_start_date_time and self.stop_date_time < orbit_stop_date_time: print ("restituted orbit already exists.") self.orbit = orbitFile self.orbitType = 'restituted'