def prepareRawLandsat(self): """ Get path to scenes from metadata and convert all to radiance """ tempPath = os.path.dirname(os.path.abspath(__file__)) + '\\temp\\' metadataPath = self.ui.satTabMTLPathLine.text() mtl_dict = l8_lst_swa_common_lib.readBasicMetadata(metadataPath) metadataBasePath = os.path.dirname(metadataPath) #print 'step1' ### B4, B5, B10, B11 to radiance B4Radiance = l8_lst_swa_common_lib.Landsat8_DN_to_radiance(1,metadataBasePath + '/' + mtl_dict['mtl_band4'],4,metadataPath) B5Radiance = l8_lst_swa_common_lib.Landsat8_DN_to_radiance(1,metadataBasePath + '/' + mtl_dict['mtl_band5'],5,metadataPath) #print 'step2' B10BrightnessTemperature = l8_lst_swa_common_lib.Landsat8_simple_temperature(1,metadataBasePath + '/' + mtl_dict['mtl_band10'],10,metadataPath) B11BrightnessTemperature = l8_lst_swa_common_lib.Landsat8_simple_temperature(1,metadataBasePath + '/' + mtl_dict['mtl_band11'],11,metadataPath) #print 'step3' ### NDVI Array NDVIArray = (B5Radiance - B4Radiance) / (B5Radiance + B4Radiance) del B4Radiance del B5Radiance #print 'done' getLSEByNDVIVectorized = np.vectorize(l8_lst_swa_core.getLSEByNDVI) LSE10Array = getLSEByNDVIVectorized(NDVIArray, 10) LSE11Array = getLSEByNDVIVectorized(NDVIArray, 11) print 'done' return mtl_dict, B10BrightnessTemperature, B11BrightnessTemperature, LSE10Array, LSE11Array
def mtlCheck(self): """ Check if metadata file is correct and all scenes are available """ if not self.ui.satTabMTLPathLine.text(): QtGui.QMessageBox.critical(None, "Error", 'Input MTL file not specified') return metadataPath = self.ui.satTabMTLPathLine.text() metadataBasePath = os.path.dirname(metadataPath) # Check if all files available try: mtl_dict = l8_lst_swa_common_lib.readBasicMetadata(metadataPath) except: QtGui.QMessageBox.critical(None, "Error", 'Invalid metadata') self.readyInterface() return # Check if extra data available #if not (os.path.isfile(metadataBasePath + '/' + mtl_dict['mtl_band4']) and os.path.isfile( # metadataBasePath + '/' + mtl_dict['mtl_band5']) and os.path.isfile( # metadataBasePath + '/' + mtl_dict['mtl_band10']) and os.path.isfile(metadataBasePath + '/' + mtl_dict['mtl_band11'])): # QtGui.QMessageBox.critical(None, "Error", 'Some of needed landsat files not exists') # return self.ui.statusSatelliteCheckBox.setChecked(True)