Higher Level Level 2 ImproPhe

[kribeetz]

Hello everyone,

I'm trying to use the force higher level Level 2 ImproPhe especially for resampling my Landsat data to 10m resolution and to spectral harmonize both. I am working with Landsat 8,9 and Sentinel 2A/B data with about 16 images from mid July to mid August of 2022. There are already in Level 2 ARD.

Therefore I created a copy of some data (6 days, from 11/07 to 28/07 for both satellites) into a lower directory folder with a .prj file and a higher directory where the parameter file is laying. Everytime I try to use the command force-higher-level l2imp.prm the processing is done within a second and doesn't process any image. It only creates this CITEME.txt and another .prj file in the higher level directory. You can see the screenshot of it here:

The .prj file is also in the X0033_Y0025 folder. When I put it into the lower directory (one folder above) it will give me the error message: Reading datacube definition failed. Copying datacube definition failed.

And a screenshot of the command line output here:

I specified all the following option in the parameter file:

• Input and Output Directory
• No Masking
• Output Options: GTiff
• Adjusted the parallel processing
• Processing Extent: Tiles are only covering X0033_Y0025
• Set Sensor Allow_list to Sen2A, Sen2B, LND08, LND09
• Selected some QAI for screening
• Adjusted the timeframe to DOY = 91 273 (so this should contain all the days I'm covering) and DATE_RANGE = 2022-01-01 2022-12-31
• Sensor Low Resolution: LND08 LND09

With the ImproPhe Parameters I'm not to sure what to put exactly. Maybe you can take a look into the parameter file and give me a suggestion on what to change. (I exported the .prm file to .txt so it could be uploaded)

Thank you!
Kristina
l2imp.txt

[davidfrantz]

Hi Kristina,

I have to admit that I didn't look into ImproPhe for a while, but please try following things:

DIR_LOWER needs to be the directory of the whole datacube, not one tile:

DIR_LOWER = /home/user/X0033_Y0025 --> DIR_LOWER = /home/user/lower_dir

lower_dir should contain datacube-definition.prj and X0033_Y0025. In X0033_Y0025, there shouldn't be a datacube-definition.prj file.

DIR_HIGHER should be the same folder as DIR_LOWER

SENSORS should only contain the high-res Sentinel-2 satellites, i.e. sen2a sen2b (small letters), and please disable SPECTRAL_ADJUST

The SEASONAL_WINDOW does not make too much sense, as you are only working with a seasonal subset anyway. The idea behind ImproPhe is to use all base images within a year, and generate seasonal slices. You could try to use July only, but in this case, only one window will suffice. For the other ImproPhe parameters, I guess the default values are mostly OK, but you might need to adapt KERNEL_SIZE to be a multiple of 30m, e.g. 3000

I hope this helps,
David