def recompose_image(layer, working_directory): bands = layer.bands image_in = layer.get_source() num_band_pir = bands['pir'] num_band_red = bands['red'] num_band_green = bands['green'] band_pir = gdal_translate_get_one_band(image_in, num_band_pir, working_directory) band_red = gdal_translate_get_one_band(image_in, num_band_red, working_directory) band_green = gdal_translate_get_one_band(image_in, num_band_green, working_directory) logger.debug("pir" + band_pir) logger.debug("red" + band_red) logger.debug("green" + band_green) output_filename = os.path.join( working_directory, os.path.splitext(os.path.basename(image_in))[0] + "pir_red_green" + os.path.splitext(os.path.basename(image_in))[1]) logger.debug("recomposed image" + output_filename) if not os.path.isfile(output_filename): OTBApplications.concatenateImages_cli([band_pir, band_red, band_green], output_filename, "uint16") return output_filename
def threshold(layer, working_directory, forms): logger.debug("threshold {}".format(forms)) image_in = layer.get_source() temp = [] i = 1 if len(forms) == 1: output_filename = os.path.join( working_directory, os.path.basename(os.path.splitext(image_in)[0]) + "_threshold" + os.path.splitext(image_in)[1]) OTBApplications.bandmath_cli([image_in], forms[0], output_filename) else: for formula in forms: output_filename = os.path.join( working_directory, os.path.basename(os.path.splitext(image_in)[0]) + "_threshold" + str(i) + os.path.splitext(image_in)[1]) OTBApplications.bandmath_cli([image_in], formula, output_filename) i += 1 temp.append(output_filename) output_filename = os.path.join( working_directory, os.path.basename(os.path.splitext(image_in)[0]) + "_threshold" + os.path.splitext(image_in)[1]) OTBApplications.concatenateImages_cli(temp, output_filename) return output_filename
def recompose_image(layer, working_directory): bands = layer.bands image_in = layer.get_source() num_band_pir = bands['pir'] num_band_red = bands['red'] num_band_green = bands['green'] band_pir = gdal_translate_get_one_band(image_in, num_band_pir, working_directory) band_red = gdal_translate_get_one_band(image_in, num_band_red, working_directory) band_green = gdal_translate_get_one_band(image_in, num_band_green, working_directory) logger.debug("pir" + band_pir) logger.debug("red" + band_red) logger.debug("green" + band_green) output_filename = os.path.join(working_directory, os.path.splitext(os.path.basename(image_in))[0] + "pir_red_green" + os.path.splitext(os.path.basename(image_in))[1]) logger.debug("recomposed image" + output_filename) if not os.path.isfile(output_filename): OTBApplications.concatenateImages_cli([band_pir, band_red, band_green], output_filename, "uint16") return output_filename
def threshold(layer, working_directory, forms): logger.debug("threshold" + str(forms)) image_in = layer.get_source() temp = [] i = 1 if len(forms) == 1: output_filename = os.path.join(working_directory, os.path.basename(os.path.splitext(image_in)[0]) + "_threshold" + os.path.splitext(image_in)[1]) OTBApplications.bandmath_cli([image_in], forms[0], output_filename) else: for formula in forms: output_filename = os.path.join(working_directory, os.path.basename(os.path.splitext(image_in)[0]) + "_threshold" + str(i) + os.path.splitext(image_in)[1]) OTBApplications.bandmath_cli([image_in], formula, output_filename) i += 1 temp.append(output_filename) output_filename = os.path.join(working_directory, os.path.basename(os.path.splitext(image_in)[0]) + "_threshold" + os.path.splitext(image_in)[1]) OTBApplications.concatenateImages_cli(temp, output_filename) return output_filename