Hi,
I made an image composite with the Composite_processing.py script. As input i have 56 images from 4 different tiles. But the result represents a image composite only on one of the 4 tiles.
Is it normal ? I need to process each tile separately for the composite ?
Thanks
In fact, I realise there are important problem with my sen2agri L3A processor in manual mode …
I read the product data structure summary and my conclusion is: I can launch a command line for the composite with 2 tiles for ex. At the end, i should have in the folder named TILES, 2 sub-folder containing the composite for each tile.
In my case, I have every time only one subfolder in TILES ending with “Tnone” --> “S2AGRI_L3A_V20170205_Tnone”
The second problem, which is linked to the first one is: when i process 4 image covering 2 different tiles, the output composite is false! When i overlap the composite with goggle earth there are no matching points
Google earth:
Composite
Projection systems are the same.
Here is the command line used (maybe there is an error, but the algo is running until the end):
python /usr/bin/composite_processing.py --syntdate 20170205 --synthalf 10 --input /mnt/archive/temp/out/test/S2A_MSIL2A_20170215T105121_N0204_R051_T31UFS_20170215T105607/S2A_OPER_SSC_L2VALD_31UFS____20170215.HDR /mnt/archive/temp/out/test/S2A_MSIL2A_20170126T105321_N0204_R051_T31UFS_20170126T105612/S2A_OPER_SSC_L2VALD_31UFS____20170126.HDR /mnt/archive/temp/out/test/S2A_MSIL2A_20170215T105121_N0204_R051_T31UFR_20170215T105607/S2A_OPER_SSC_L2VALD_31UFR____20170215.HDR /mnt/archive/temp/out/test/S2A_MSIL2A_20170126T105321_N0204_R051_T31UFR_20170126T105612/S2A_OPER_SSC_L2VALD_31UFR____20170126.HDR --res 10 --outdir /mnt/archive/temp/out/test/CompositeL3A --bandsmap /usr/share/sen2agri/bands_mapping_s2_L8.txt --scatteringcoef /usr/share/sen2agri/scattering_coeffs_10m.txt
precision: I think the same problem will happens with other processors like LA4 or L4B.
Thanks for your help, don’t know how to solve this and i really need it.
Hello,
Unfortunately the script does not supports this mode of executing several tiles in the same time and it will not also be updated in version 1.8.3. In this moment, the script can be run for only one tile time series at a time. It will be updated in the version 2.0 to manage this mode.
In this moment, you can split the creation of the product into 2 commands like below. Additionally, you need to specify in this moment the tile id:
> python /usr/bin/composite_processing.py --syntdate 20170205 --synthalf 10 --input /mnt/archive/temp/out/test/S2A_MSIL2A_20170215T105121_N0204_R051_T31UFR_20170215T105607/S2A_OPER_SSC_L2VALD_31UFR____20170215.HDR /mnt/archive/temp/out/test/S2A_MSIL2A_20170126T105321_N0204_R051_T31UFR_20170126T105612/S2A_OPER_SSC_L2VALD_31UFR____20170126.HDR --res 10 --outdir /mnt/archive/temp/out/test/CompositeL3A --bandsmap /usr/share/sen2agri/bands_mapping_s2_L8.txt --scatteringcoef /usr/share/sen2agri/scattering_coeffs_10m.txt --tileid 31UFR
> python /usr/bin/composite_processing.py --syntdate 20170205 --synthalf 10 --input /mnt/archive/temp/out/test/S2A_MSIL2A_20170215T105121_N0204_R051_T31UFS_20170215T105607/S2A_OPER_SSC_L2VALD_31UFS____20170215.HDR /mnt/archive/temp/out/test/S2A_MSIL2A_20170126T105321_N0204_R051_T31UFS_20170126T105612/S2A_OPER_SSC_L2VALD_31UFS____20170126.HDR --res 10 --outdir /mnt/archive/temp/out/test/CompositeL3A --bandsmap /usr/share/sen2agri/bands_mapping_s2_L8.txt --scatteringcoef /usr/share/sen2agri/scattering_coeffs_10m.txt --tileid 31UFS
many thanks for your response, I was not able to determine if the problem was coming from my installation of Sen2agri or if it was the right behaviour of the algo.
Thanks
