ERS stack processing

README.md

@@ -22,7 +22,7 @@ utility in a general sense for building other types of software packages.  In
 its InSAR aspect ISCE supports data from many space-borne satellites and one
 air-borne platform.  We continue to increase the number of sensors supported.
 At this time the sensors that are supported are the following: ALOS, ALOS2,
-COSMO_SKYMED, ENVISAT, ERS, KOMPSAT5, RADARSAT1, RADARSAT2, RISAT1, Sentinel1,
+COSMO_SKYMED, ENVISAT, ERS, KOMPSAT5, LUTAN1, RADARSAT1, RADARSAT2, RISAT1, Sentinel1,
 TERRASARX, UAVSAR and SAOCOM1A.
 
 ## Contents

applications/make_raw.py

@@ -226,7 +226,7 @@ def calculateSquint(self):
         self.frame.squintAngle = math.radians(self.squint)
 
     def make_raw(self):
-        from  isceobj.Image import createRawImage, createSlcImage
+        from isceobj.Image import createRawImage, createSlcImage
         self.activateInputPorts()
 
         # Parse the image metadata and extract the image

components/isceobj/Alos2Proc/runFrameOffset.py

@@ -97,18 +97,23 @@ def frameOffset(track, image, outputfile, crossCorrelation=True, matchingMode=0)
         swath1 = track.frames[j-1].swaths[0]
         swath2 = track.frames[j].swaths[0]
 
+        #consider frame/swath azimuth sensing start differences caused by swath mosaicking
+        # tested with alos2App.py, alos2burstApp.py and alosStack
+        frame1 = track.frames[j-1]
+        frame2 = track.frames[j]
+        delta_az = -((frame2.sensingStart - frame1.sensingStart).total_seconds() - (swath2.sensingStart - swath1.sensingStart).total_seconds()) / swath1.azimuthLineInterval
 
         #offset from geometry
         offsetGeometrical = computeFrameOffset(swath1, swath2)
         rangeOffsetGeometrical.append(offsetGeometrical[0])
-        azimuthOffsetGeometrical.append(offsetGeometrical[1])
+        azimuthOffsetGeometrical.append(offsetGeometrical[1]+delta_az)
 
         #offset from cross-correlation
         if crossCorrelation:
             offsetMatching = estimateFrameOffset(swath1, swath2, image1, image2, matchingMode=matchingMode)
             if offsetMatching != None:
                 rangeOffsetMatching.append(offsetMatching[0])
-                azimuthOffsetMatching.append(offsetMatching[1])
+                azimuthOffsetMatching.append(offsetMatching[1]+delta_az)
             else:
                 print('******************************************************************')
                 print('WARNING: bad matching offset, we are forced to use')

components/isceobj/Orbit/ODR.py

@@ -184,15 +184,11 @@ def parse(self):
 
     def parseLine(self,line):        
         arc = Arc()
-        arc.number = line[0:3] # Arc number
-        arc.start = datetime.datetime.strptime(line[5:17],'%y%m%d %H:%M') # Starting time for the arc
-        arc.stop = datetime.datetime.strptime(line[20:32],'%y%m%d %H:%M') # End time for the arc
-        arc.slrResidual = line[34:38] # Satellite laser ranging residual in cm        
-        arc.crossOver = line[39:43]
-        arc.altimeter = line[45:49]
-        arc.repeat = line[51:57] # Repeat cycle in days
-        arc.version = line[58:61] # Version number
-        arc.precise = datetime.datetime.strptime(line[63:78],'%y%m%d %H:%M:%S') # Starting time of the precise segment of the arc
+        ###### Change to also work for the new arclist format in the recalcuated ODR orbit files.#######
+        arc.number = line.split()[0] # Arc number
+        arc.start = datetime.datetime.strptime(" ".join(line.split()[1:3]),'%y%m%d %H:%M') # Starting time for the arc
+        arc.stop = datetime.datetime.strptime(" ".join(line.split()[4:6]),'%y%m%d %H:%M') # End time for the arc
+        arc.precise = datetime.datetime.strptime(" ".join(line.split()[-2:]),'%y%m%d %H:%M:%S') # Starting time of the precise segment of the arc
 
         return arc
 

components/isceobj/Sensor/ALOS2.py

@@ -106,7 +106,7 @@ def __init__(self, name=''):
         self.leaderFile = None
         self.imageFile = None
 
-        #####Soecific doppler functions for ALOS2
+        #####Specific doppler functions for ALOS2
         self.doppler_coeff = None
         self.azfmrate_coeff = None
         self.lineDirection = None

components/isceobj/Sensor/CMakeLists.txt

@@ -41,6 +41,7 @@ set(installfiles
     UAVSAR_RPI.py
     UAVSAR_Stack.py
     SAOCOM_SLC.py
+    Lutan1.py
     )
 
 if(HDF5_FOUND)