def runSimamp(outdir, hname='z.rdr'): from iscesys.StdOEL.StdOELPy import create_writer #####Run simamp stdWriter = create_writer("log", "", True, filename='sim.log') objShade = isceobj.createSimamplitude() objShade.setStdWriter(stdWriter) hgtImage = isceobj.createImage() hgtImage.load(os.path.join(outdir, hname) + '.xml') hgtImage.setAccessMode('read') hgtImage.createImage() simImage = isceobj.createImage() simImage.setFilename(os.path.join(outdir, 'simamp.rdr')) simImage.dataType = 'FLOAT' simImage.setAccessMode('write') simImage.setWidth(hgtImage.getWidth()) simImage.createImage() objShade.simamplitude(hgtImage, simImage, shade=3.0) simImage.renderHdr() hgtImage.finalizeImage() simImage.finalizeImage() return
def runShadecpx2rg(self): imageAmp = self._insar.getResampAmpImage() widthAmp = imageAmp.getWidth() endian = self._insar.getMachineEndianness() filenameSimAmp = self._insar.getSimAmpImageName() objSimAmp = isceobj.createImage() widthSimAmp = widthAmp objSimAmp.initImage(filenameSimAmp, 'read', widthSimAmp, 'FLOAT') imageSimAmp = isceobj.createImage() IU.copyAttributes(objSimAmp, imageSimAmp) self._insar.setSimAmpImage(imageSimAmp) objSimAmp.setAccessMode('write') objSimAmp.createImage() filenameHt = self._insar.getHeightFilename() widthHgtImage = widthAmp # they have same width by construction objHgtImage = isceobj.createImage() objHgtImage.initImage(filenameHt, 'read', widthHgtImage, 'FLOAT') imageHgt = isceobj.createImage() IU.copyAttributes(objHgtImage, imageHgt) self._insar.setHeightTopoImage(imageHgt) objHgtImage.createImage() logger.info("Running Shadecpx2rg") objShade = isceobj.createSimamplitude() #set the tag used in the outfile. each message is precided by this tag #is the writer is not of "file" type the call has no effect self._stdWriter.setFileTag("simamplitude", "log") self._stdWriter.setFileTag("simamplitude", "err") self._stdWriter.setFileTag("simamplitude", "out") objShade.setStdWriter(self._stdWriter) shade = self._insar.getShadeFactor() objShade.simamplitude(objHgtImage, objSimAmp, shade=shade) # Record the inputs and outputs from isceobj.Catalog import recordInputsAndOutputs recordInputsAndOutputs(self._insar.procDoc, objShade, "runSimamplitude", \ logger, "runSimamplitude") objHgtImage.finalizeImage() objSimAmp.finalizeImage() objSimAmp.renderHdr()
def run(widthAmp, infos, stdWriter, catalog=None, sceneid='NO_ID'): logger.info("Running shadecpx2rg: %s" % sceneid) endian = infos['machineEndianness'] filenameSimAmp = infos['outputPath'] + '.' + infos['simAmpImageName'] filenameHt = infos['outputPath'] + '.' + infos['heightFilename'] shade = infos['shadeFactor'] objSimAmp = isceobj.createImage() widthSimAmp = widthAmp objSimAmp.initImage(filenameSimAmp, 'read', widthSimAmp, 'FLOAT') imageSimAmp = isceobj.createImage() IU.copyAttributes(objSimAmp, imageSimAmp) objSimAmp.setAccessMode('write') objSimAmp.createImage() widthHgtImage = widthAmp # they have same width by construction objHgtImage = isceobj.createImage() objHgtImage.initImage(filenameHt, 'read', widthHgtImage, 'FLOAT') imageHgt = isceobj.createImage() IU.copyAttributes(objHgtImage, imageHgt) objHgtImage.createImage() objShade = isceobj.createSimamplitude() #set the tag used in the outfile. each message is precided by this tag #if the writer is not of "file" type the call has no effect objShade.stdWriter = stdWriter.set_file_tags("simamplitude", "log", "err", "out") objShade.simamplitude(objHgtImage, objSimAmp, shade=shade) if catalog is not None: # Record the inputs and outputs isceobj.Catalog.recordInputsAndOutputs(catalog, objShade, "runSimamplitude.%s" % sceneid, logger, "runSimamplitude.%s" % sceneid) objHgtImage.finalizeImage() objSimAmp.finalizeImage() objSimAmp.renderHdr() return imageSimAmp, imageHgt
def runShadecpx2rg(self): width = self.insar.formSLC1.slcImage.getWidth() filenameSimAmp = self._insar.getSimAmpImageName() objSimAmp = isceobj.createImage() widthSimAmp = width objSimAmp.initImage(filenameSimAmp,'read',widthSimAmp,'FLOAT') imageSimAmp = isceobj.createImage() IU.copyAttributes(objSimAmp, imageSimAmp) self._insar.setSimAmpImage(imageSimAmp) objSimAmp.setAccessMode('write') objSimAmp.createImage() filenameHt = self._insar.getHeightFilename() widthHgtImage = width # they have same width by construction objHgtImage = isceobj.createImage() objHgtImage.initImage(filenameHt,'read',widthHgtImage,'DOUBLE') imageHgt = isceobj.createImage() IU.copyAttributes(objHgtImage, imageHgt) self._insar.setHeightTopoImage(imageHgt) objHgtImage.setCaster('read','FLOAT') objHgtImage.createImage() logger.info("Running Shadecpx2rg") objShade = isceobj.createSimamplitude() #set the tag used in the outfile. each message is precided by this tag #is the writer is not of "file" type the call has no effect self._stdWriter.setFileTag("simamplitude", "log") self._stdWriter.setFileTag("simamplitude", "err") self._stdWriter.setFileTag("simamplitude", "out") objShade.setStdWriter(self._stdWriter) shade = self._insar.getShadeFactor() objShade.simamplitude(objHgtImage, objSimAmp, shade=shade) objHgtImage.finalizeImage() objSimAmp.finalizeImage() objSimAmp.renderHdr()