From 6e36a73fb0a0e34cc2e30bb06ca5995bd14deb3a Mon Sep 17 00:00:00 2001
From: Zhang Yunjun <yunjunz@users.noreply.github.com>
Date: Mon, 18 Mar 2019 19:00:18 -0400
Subject: [PATCH] stripmapStack/unwrap: enable corrlooks for SNAPHU

Stack:
+ add alks and rlks while writing config file for unwrap

unwrap:
+ calculate the az/rg spacing/resolution and the equivalent number of independent looks in extractInfoFromPickle()
+ uncomment corrlooks setting in runUnwrap()
---
 contrib/stack/stripmapStack/Stack.py  |  2 ++
 contrib/stack/stripmapStack/unwrap.py | 16 +++++++++++++---
 2 files changed, 15 insertions(+), 3 deletions(-)

diff --git a/contrib/stack/stripmapStack/Stack.py b/contrib/stack/stripmapStack/Stack.py
index 74d0df4..b0a116e 100755
--- a/contrib/stack/stripmapStack/Stack.py
+++ b/contrib/stack/stripmapStack/Stack.py
@@ -197,6 +197,8 @@ class config(object):
         self.f.write('nomcf : ' + self.noMCF + '\n')
         self.f.write('master : ' + self.master + '\n')
         self.f.write('defomax : ' + self.defoMax + '\n')
+        self.f.write('alks : ' + self.alks + '\n')
+        self.f.write('rlks : ' + self.rlks + '\n')
         self.f.write('method : ' + self.unwMethod + '\n')
         self.f.write('##########################'+'\n')
 
diff --git a/contrib/stack/stripmapStack/unwrap.py b/contrib/stack/stripmapStack/unwrap.py
index bae4893..f1de7f9 100755
--- a/contrib/stack/stripmapStack/unwrap.py
+++ b/contrib/stack/stripmapStack/unwrap.py
@@ -113,9 +113,19 @@ def extractInfoFromPickle(pckfile, inps):
     data['earthRadius'] = elp.local_radius_of_curvature(llh.lat, hdg)
     
     #azspacing  = burst.azimuthTimeInterval * sv.getScalarVelocity()
-    azres = 20.0 
+    #azres = 20.0 
+    azspacing = sv.getScalarVelocity() / burst.PRF
+    azres = burst.platform.antennaLength / 2.0
+    azfact = azres / azspacing
+
+    burst.getInstrument()
+    rgBandwidth = burst.instrument.pulseLength * burst.instrument.chirpSlope
+    rgres = abs(SPEED_OF_LIGHT / (2.0 * rgBandwidth))
+    rgspacing = burst.instrument.rangePixelSize
+    rgfact = rgres / rgspacing
 
     #data['corrlooks'] = inps.rglooks * inps.azlooks * azspacing / azres
+    data['corrlooks'] = inps.rglooks * inps.azlooks / (azfact * rgfact)
     data['rglooks'] = inps.rglooks
     data['azlooks'] = inps.azlooks
 
@@ -149,7 +159,7 @@ def runUnwrap(infile, outfile, corfile, config, costMode = None,initMethod = Non
     altitude   = config['altitude']
     rangeLooks = config['rglooks']
     azimuthLooks = config['azlooks']
-    #corrLooks = config['corrlooks']
+    corrLooks = config['corrlooks']
     maxComponents = 20
 
     snp = Snaphu()
@@ -163,7 +173,7 @@ def runUnwrap(infile, outfile, corfile, config, costMode = None,initMethod = Non
     snp.setAltitude(altitude)
     snp.setCorrfile(corfile)
     snp.setInitMethod(initMethod)
-   # snp.setCorrLooks(corrLooks)
+    snp.setCorrLooks(corrLooks)
     snp.setMaxComponents(maxComponents)
     snp.setDefoMaxCycles(defomax)
     snp.setRangeLooks(rangeLooks)