small fix in EnFireMap ingest_data script

enmapboxprocessing/algorithm/translaterasteralgorithm.py

@@ -234,8 +234,12 @@ def isBadBand(bandNo: int) -> bool:
                     spectralBandList = [i + 1 for i in range(spectralRaster.bandCount())]
 
                 wavelength = np.array([reader.wavelength(bandNo) for bandNo in bandList])
-                bandList = [int(np.argmin(np.abs(wavelength - spectralReader.wavelength(bandNo))) + 1)
-                            for bandNo in spectralBandList]
+
+                bandList = list()
+                for targetBandNo in spectralBandList:
+                    targetWavelength = spectralReader.wavelength(targetBandNo)
+                    bandNo = int(np.argmin(np.abs(wavelength - targetWavelength)) + 1)
+                    bandList.append(bandNo)
 
             if bandList is None:
                 nBands = raster.bandCount()

scripts/EnFireMAP/bad_wavelengths.txt

@@ -3,3 +3,5 @@
 1366.69
 1378.6
 1390.48
+1461.1
+1759.51

scripts/EnFireMAP/ingest_data.py

@@ -15,7 +15,6 @@
 from enmapboxprocessing.rasterreader import RasterReader
 from enmapboxprocessing.rasterwriter import RasterWriter
 from enmapboxprocessing.utils import Utils
-
 from qgis.core import QgsPointXY, QgsCoordinateReferenceSystem, QgsCoordinateTransform, \
     QgsProject, QgsVectorLayer, QgsGeometry, Qgis, QgsRectangle
 
@@ -60,7 +59,7 @@ def prepareSpectralImages():
                 alg = ImportEnmapL2AAlgorithm()
                 parameters = {
                     alg.P_FILE: xmlFilename,
-                    alg.P_SET_BAD_BANDS: True,
+                    alg.P_SET_BAD_BANDS: False,
                     alg.P_EXCLUDE_BAD_BANDS: False,
                     alg.P_DETECTOR_OVERLAP: alg.OrderByWavelengthOverlapOption,
                     alg.P_OUTPUT_RASTER: xmlFilename.replace('METADATA.XML', 'SPECTRAL_IMAGE_FULL.vrt'),
@@ -220,4 +219,3 @@ def auxDeriveSceneBoundingPolygon(xmlFilename) -> QgsGeometry:
 prepareSpectralImages()
 copyMetadataXml()
 ingestData()
-print('done')

tests/enmap-box/enmapboxprocessing/algorithm/test_TranslateRasterAlgorithm.py

@@ -3,11 +3,11 @@
 import numpy as np
 from osgeo import gdal
 
-from enmapboxtestdata import enmap, hires
 from enmapboxprocessing.algorithm.testcase import TestCase
 from enmapboxprocessing.algorithm.translaterasteralgorithm import TranslateRasterAlgorithm
 from enmapboxprocessing.rasterreader import RasterReader
 from enmapboxprocessing.utils import Utils
+from enmapboxtestdata import enmap, hires
 from enmapboxtestdata import water_mask_30m, enmap_grid_300m
 from qgis.core import QgsRectangle, QgsCoordinateReferenceSystem, QgsRasterLayer, QgsRasterRenderer, Qgis