134 lines
4.1 KiB
Python
134 lines
4.1 KiB
Python
|
#
|
||
|
|
#
|
||
|
|
|
||
|
|
import isce
|
||
|
|
import isceobj
|
||
|
|
import stdproc
|
||
|
|
from isceobj.Util.Poly2D import Poly2D
|
||
|
|
import logging
|
||
|
|
from isceobj.Util.decorators import use_api
|
||
|
|
|
||
|
|
import os
|
||
|
|
import numpy as np
|
||
|
|
import shelve
|
||
|
|
|
||
|
|
logger = logging.getLogger('isce.insar.runResampleSlc')
|
||
|
|
|
||
|
|
def runResampleSlc(self, kind='coarse'):
|
||
|
|
'''
|
||
|
|
Kind can either be coarse, refined or fine.
|
||
|
|
'''
|
||
|
|
|
||
|
|
if kind not in ['coarse', 'refined', 'fine']:
|
||
|
|
raise Exception('Unknown operation type {0} in runResampleSlc'.format(kind))
|
||
|
|
|
||
|
|
if kind == 'fine':
|
||
|
|
if not self.doRubbersheeting:
|
||
|
|
print('Rubber sheeting not requested, skipping resampling ....')
|
||
|
|
return
|
||
|
|
|
||
|
|
logger.info("Resampling slave SLC")
|
||
|
|
|
||
|
|
slaveFrame = self._insar.loadProduct( self._insar.slaveSlcCropProduct)
|
||
|
|
masterFrame = self._insar.loadProduct( self._insar.masterSlcCropProduct)
|
||
|
|
|
||
|
|
inimg = isceobj.createSlcImage()
|
||
|
|
inimg.load(slaveFrame.getImage().filename + '.xml')
|
||
|
|
inimg.setAccessMode('READ')
|
||
|
|
|
||
|
|
prf = slaveFrame.PRF
|
||
|
|
|
||
|
|
doppler = slaveFrame._dopplerVsPixel
|
||
|
|
coeffs = [2*np.pi*val/prf for val in doppler]
|
||
|
|
|
||
|
|
dpoly = Poly2D()
|
||
|
|
dpoly.initPoly(rangeOrder=len(coeffs)-1, azimuthOrder=0, coeffs=[coeffs])
|
||
|
|
|
||
|
|
rObj = stdproc.createResamp_slc()
|
||
|
|
rObj.slantRangePixelSpacing = slaveFrame.getInstrument().getRangePixelSize()
|
||
|
|
rObj.radarWavelength = slaveFrame.getInstrument().getRadarWavelength()
|
||
|
|
rObj.dopplerPoly = dpoly
|
||
|
|
|
||
|
|
# for now let's start with None polynomial. Later this should change to
|
||
|
|
# the misregistration polynomial
|
||
|
|
|
||
|
|
misregFile = os.path.join(self.insar.misregDirname, self.insar.misregFilename)
|
||
|
|
if ((kind in ['refined','fine']) and os.path.exists(misregFile+'_az.xml')):
|
||
|
|
azpoly = self._insar.loadProduct(misregFile + '_az.xml')
|
||
|
|
rgpoly = self._insar.loadProduct(misregFile + '_rg.xml')
|
||
|
|
else:
|
||
|
|
print(misregFile , " does not exist.")
|
||
|
|
azpoly = None
|
||
|
|
rgpoly = None
|
||
|
|
|
||
|
|
rObj.azimuthOffsetsPoly = azpoly
|
||
|
|
rObj.rangeOffsetsPoly = rgpoly
|
||
|
|
rObj.imageIn = inimg
|
||
|
|
|
||
|
|
#Since the app is based on geometry module we expect pixel-by-pixel offset
|
||
|
|
#field
|
||
|
|
offsetsDir = self.insar.offsetsDirname
|
||
|
|
rgname = os.path.join(offsetsDir, self.insar.rangeOffsetFilename)
|
||
|
|
if kind in ['coarse', 'refined']:
|
||
|
|
azname = os.path.join(offsetsDir, self.insar.azimuthOffsetFilename)
|
||
|
|
else:
|
||
|
|
azname = os.path.join(offsetsDir, self.insar.azimuthRubbersheetFilename)
|
||
|
|
|
||
|
|
rngImg = isceobj.createImage()
|
||
|
|
rngImg.load(rgname + '.xml')
|
||
|
|
rngImg.setAccessMode('READ')
|
||
|
|
|
||
|
|
aziImg = isceobj.createImage()
|
||
|
|
aziImg.load(azname + '.xml')
|
||
|
|
aziImg.setAccessMode('READ')
|
||
|
|
|
||
|
|
width = rngImg.getWidth()
|
||
|
|
length = rngImg.getLength()
|
||
|
|
|
||
|
|
|
||
|
|
flatten = True
|
||
|
|
rObj.flatten = flatten
|
||
|
|
rObj.outputWidth = width
|
||
|
|
rObj.outputLines = length
|
||
|
|
rObj.residualRangeImage = rngImg
|
||
|
|
rObj.residualAzimuthImage = aziImg
|
||
|
|
|
||
|
|
if masterFrame is not None:
|
||
|
|
rObj.startingRange = slaveFrame.startingRange
|
||
|
|
rObj.referenceStartingRange = masterFrame.startingRange
|
||
|
|
rObj.referenceSlantRangePixelSpacing = masterFrame.getInstrument().getRangePixelSize()
|
||
|
|
rObj.referenceWavelength = masterFrame.getInstrument().getRadarWavelength()
|
||
|
|
|
||
|
|
|
||
|
|
# preparing the output directory for coregistered slave slc
|
||
|
|
coregDir = self.insar.coregDirname
|
||
|
|
|
||
|
|
if os.path.isdir(coregDir):
|
||
|
|
logger.info('Geometry directory {0} already exists.'.format(coregDir))
|
||
|
|
else:
|
||
|
|
os.makedirs(coregDir)
|
||
|
|
|
||
|
|
# output file name of the coregistered slave slc
|
||
|
|
img = slaveFrame.getImage()
|
||
|
|
|
||
|
|
if kind == 'coarse':
|
||
|
|
coregFilename = os.path.join(coregDir , self._insar.coarseCoregFilename)
|
||
|
|
elif kind == 'refined':
|
||
|
|
coregFilename = os.path.join(coregDir, self._insar.refinedCoregFilename)
|
||
|
|
elif kind == 'fine':
|
||
|
|
coregFilename = os.path.join(coregDir, self._insar.fineCoregFilename)
|
||
|
|
else:
|
||
|
|
print('Exception: Should not have gotten to this stage')
|
||
|
|
|
||
|
|
imgOut = isceobj.createSlcImage()
|
||
|
|
imgOut.setWidth(width)
|
||
|
|
imgOut.filename = coregFilename
|
||
|
|
imgOut.setAccessMode('write')
|
||
|
|
|
||
|
|
rObj.resamp_slc(imageOut=imgOut)
|
||
|
|
|
||
|
|
imgOut.renderHdr()
|
||
|
|
|
||
|
|
return
|
||
|
|
|