From 1b352372bc24b927f999c2094dd0bd45a35e324e Mon Sep 17 00:00:00 2001 From: Bryan Marfito Date: Sun, 10 Sep 2023 11:29:46 +0800 Subject: [PATCH 1/4] Use default values instead of using DOPIQ for unpackFrame_ROIPAC_raw.py --- .../stripmapStack/unpackFrame_ROIPAC_raw.py | 147 ++++++++---------- 1 file changed, 64 insertions(+), 83 deletions(-) diff --git a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py index 10988dc4..439752a8 100755 --- a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py +++ b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py @@ -1,83 +1,64 @@ -#!/usr/bin/env python3 - -import isce -from isceobj.Sensor import createSensor -import shelve -import argparse -import glob -from isceobj.Util import Poly1D -from isceobj.Planet.AstronomicalHandbook import Const -import os -from mroipac.dopiq.DopIQ import DopIQ -import copy - -def cmdLineParse(): - ''' - Command line parser. - ''' - - parser = argparse.ArgumentParser(description='Unpack raw data and store metadata in pickle file.') - parser.add_argument('-i','--input', dest='rawfile', type=str, - required=True, help='Input ROI_PAC file') - parser.add_argument('-r','--hdr', dest='hdrfile', type=str, - required=True, help='Input hdr (orbit) file') - parser.add_argument('-o', '--output', dest='slcdir', type=str, - required=True, help='Output data directory') - - return parser.parse_args() - - -def unpack(rawname, hdrname, slcname): - ''' - Unpack raw to binary file. - ''' - - if not os.path.isdir(slcname): - os.mkdir(slcname) - - date = os.path.basename(slcname) - obj = createSensor('ROI_PAC') - obj.configure() - obj._rawFile = rawname - obj._hdrFile = hdrname - obj.output = os.path.join(slcname, date+'.raw') - - print(obj._rawFile) - print(obj._hdrFile) - print(obj.output) - obj.extractImage() - obj.frame.getImage().renderHdr() - - - #####Estimate doppler - dop = DopIQ() - dop.configure() - - img = copy.deepcopy(obj.frame.getImage()) - img.setAccessMode('READ') - - dop.wireInputPort('frame', object=obj.frame) - dop.wireInputPort('instrument', object=obj.frame.instrument) - dop.wireInputPort('image', object=img) - dop.calculateDoppler() - dop.fitDoppler() - fit = dop.quadratic - coef = [fit['a'], fit['b'], fit['c']] - - print(coef) - obj.frame._dopplerVsPixel = [x*obj.frame.PRF for x in coef] - - pickName = os.path.join(slcname, 'raw') - with shelve.open(pickName) as db: - db['frame'] = obj.frame - -if __name__ == '__main__': - ''' - Main driver. - ''' - - inps = cmdLineParse() - if inps.slcdir.endswith('/'): - inps.slcdir = inps.slcdir[:-1] - - unpack(inps.rawfile, inps.hdrfile, inps.slcdir) +#!/usr/bin/env python3 + +import isce +from isceobj.Sensor import createSensor +import shelve +import argparse +import glob +from isceobj.Util import Poly1D +from isceobj.Planet.AstronomicalHandbook import Const +import os +import copy + +def cmdLineParse(): + ''' + Command line parser. + ''' + + parser = argparse.ArgumentParser(description='Unpack raw data and store metadata in pickle file.') + parser.add_argument('-i','--input', dest='rawfile', type=str, + required=True, help='Input ROI_PAC file') + parser.add_argument('-r','--hdr', dest='hdrfile', type=str, + required=True, help='Input hdr (orbit) file') + parser.add_argument('-o', '--output', dest='slcdir', type=str, + required=True, help='Output data directory') + + return parser.parse_args() + + +def unpack(rawname, hdrname, slcname): + ''' + Unpack raw to binary file. + ''' + + if not os.path.isdir(slcname): + os.mkdir(slcname) + + date = os.path.basename(slcname) + obj = createSensor('ROI_PAC') + obj.configure() + obj._rawFile = os.path.abspath(rawname) + obj._hdrFile = hdrname + obj.output = os.path.join(slcname, date+'.raw') + + print(obj._rawFile) + print(obj._hdrFile) + print(obj.output) + obj.extractImage() + obj.frame.getImage().renderHdr() + obj.extractDoppler() + + pickName = os.path.join(slcname, 'raw') + with shelve.open(pickName) as db: + db['frame'] = obj.frame + +if __name__ == '__main__': + ''' + Main driver. + ''' + + inps = cmdLineParse() + if inps.slcdir.endswith('/'): + inps.slcdir = inps.slcdir[:-1] + + unpack(inps.rawfile, inps.hdrfile, inps.slcdir) From 49269d7f4f5620a767dc3ea930611fd83ca0d239 Mon Sep 17 00:00:00 2001 From: Bryan Marfito Date: Sun, 10 Sep 2023 11:37:20 +0800 Subject: [PATCH 2/4] Removed unused libraries --- contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py | 4 ---- 1 file changed, 4 deletions(-) diff --git a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py index 439752a8..c4a0dc13 100755 --- a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py +++ b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py @@ -4,11 +4,7 @@ from isceobj.Sensor import createSensor import shelve import argparse -import glob -from isceobj.Util import Poly1D -from isceobj.Planet.AstronomicalHandbook import Const import os -import copy def cmdLineParse(): ''' From 6ee519e0a6cf4a1dfa03b2105f90505f3df3bd64 Mon Sep 17 00:00:00 2001 From: Bryan Marfito Date: Tue, 9 Apr 2024 17:58:16 +0800 Subject: [PATCH 3/4] Update unpackFrame_ROIPAC_raw.py --- .../stripmapStack/unpackFrame_ROIPAC_raw.py | 29 ++++++++++++++++--- 1 file changed, 25 insertions(+), 4 deletions(-) diff --git a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py index c4a0dc13..9635144f 100755 --- a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py +++ b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py @@ -19,17 +19,19 @@ def cmdLineParse(): parser.add_argument('-o', '--output', dest='slcdir', type=str, required=True, help='Output data directory') + parser.add_argument('-d','--default', dest='defaultDoppler', type=str, action='store_true' + required=False, help='use DEFAULT to estimate doppler') + return parser.parse_args() -def unpack(rawname, hdrname, slcname): +def unpack(rawname, hdrname, slcname, default): ''' Unpack raw to binary file. ''' if not os.path.isdir(slcname): os.mkdir(slcname) - date = os.path.basename(slcname) obj = createSensor('ROI_PAC') obj.configure() @@ -42,7 +44,26 @@ def unpack(rawname, hdrname, slcname): print(obj.output) obj.extractImage() obj.frame.getImage().renderHdr() - obj.extractDoppler() + + if defaultDoppler == True: + obj.extractDoppler() + else: + dop = DopIQ() + dop.configure() + + img = copy.deepcopy(obj.frame.getImage()) + img.setAccessMode('READ') + + dop.wireInputPort('frame', object=obj.frame) + dop.wireInputPort('instrument', object=obj.frame.instrument) + dop.wireInputPort('image', object=img) + dop.calculateDoppler() + dop.fitDoppler() + fit = dop.quadratic + coef = [fit['a'], fit['b'], fit['c']] + + print(coef) + obj.frame._dopplerVsPixel = [x*obj.frame.PRF for x in coef] pickName = os.path.join(slcname, 'raw') with shelve.open(pickName) as db: @@ -57,4 +78,4 @@ def unpack(rawname, hdrname, slcname): if inps.slcdir.endswith('/'): inps.slcdir = inps.slcdir[:-1] - unpack(inps.rawfile, inps.hdrfile, inps.slcdir) + unpack(inps.rawfile, inps.hdrfile, inps.slcdir, inps.defaultDoppler) From 358a6f5fcb5b2b39cfbec7c85bfb600cf041c2fe Mon Sep 17 00:00:00 2001 From: Bryan Marfito Date: Tue, 9 Apr 2024 17:59:21 +0800 Subject: [PATCH 4/4] Update unpackFrame_ROIPAC_raw.py --- contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py index 9635144f..8584d87c 100755 --- a/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py +++ b/contrib/stack/stripmapStack/unpackFrame_ROIPAC_raw.py @@ -5,6 +5,8 @@ import shelve import argparse import os +from mroipac.dopiq.DopIQ import DopIQ +import copy def cmdLineParse(): '''