Pléiade OpticalCalibration Issue

Dear All,
I tried to do OpticalCalibration with TOC correction on two Pleiade datasets.
The first dataset (located in Canada) didn’t worked, the resulting images are black with red and yellow contour lines.
The second dataset (located in patagonia) worked (no more blue shadow), but when I am trying to orthorectify it using MicMac (the French IGN photogrammetry opensource software), I just have a black image.

Do you know those problems, especially concerning the first dataset ? Is there a way to solve it ?
Thanks for your help !

Anaïs

Dear @AnaisE,
Could you provide us with more details on your setup?

  • Which operating system are you using?
  • Which version of the OTB? Did you build it? Or are you using the provided binary files?
  • Are you using QGIS plugins? Or Monteverdi? Or the command line?
    Sincerely.
    Julien.

Thanks for your answer, I forgot a lot of parameters ^^’
I am operating under Ubuntu, and OTB version 7.1.0, and I guess I am using the binary files, because I am launching Monteverdi using mapla.sh…
I installed OTB using OTB-7.1.0-Linux64.run
Does it help you ?

Anaïs

Dear @AnaisE,
If I understand correctly, you are using Monteverdi to process your datasets? Could you share with us a screenshot of the parameters you set for the processing? And also the logs you obtain?
Sincerely.
Julien.

Hello !
I joined pictures of the parameter with the maritime aerosol model, and a picture of the result with no aerosol model choosed (which is the better choice for now).
OTB_image OTB_param.pdf (244.5 KB)

Here is the log :

2020-09-10 11:01:22 (INFO) OpticalCalibration:
File: /data/Labrador/Traitements/Dossiers_images/zip/2018-006-Sci-P004-Univ-Rouen-C.Colange_SO18010794-6-01_DS_PHR1A_201909131530380_FR1_PX_W062N56_0918_01258/TPP1601254289/IMG_PHR1A_MS_004/IMG_PHR1A_MS_201909131530380_SEN_4527881101-004_R1C1.TIF
Sensor detected: PHR 1A
Parameters extract from input image:
Acquisition Day: 13
Acquisition Month: 9
Acquisition Year: 2019
Acquisition Sun Elevation Angle: 36.7725
Acquisition Sun Azimuth Angle: 170.305
Acquisition Viewing Elevation Angle: 75.3594
Acquisition Viewing Azimuth Angle: 356.039
Acquisition gain (per band): 10.46 9.34 9.38 15.69
Acquisition bias (per band): 2.4056 3.57564 3.35568 1.1145
Solar Irradiance (per band): 1594 1830 1915 1060
2020-09-10 11:03:14 (INFO) OpticalCalibration: Default RAM limit for OTB is 256 MB
2020-09-10 11:03:14 (INFO) OpticalCalibration: GDAL maximum cache size is 6446 MB
2020-09-10 11:03:14 (INFO) OpticalCalibration: OTB will use at most 36 threads
2020-09-10 11:03:14 (INFO) OpticalCalibration: Compute Top of Atmosphere reflectance
2020-09-10 11:03:14 (INFO) OpticalCalibration: Clamp values between [0, 100]
2020-09-10 11:04:41 (INFO) OpticalCalibration: Execution took 86.552 sec

Thanks !

Dear @AnaisE,
I don’t see what could cause your issue. Could you use the command line to run your processing?

./otbcli_OpticalCalibration -in /data/Labrador/Traitements/Dossiers_images/zip/2018-006-Sci-P004-Univ-Rouen-C.Colange_SO18010794-6-01_DS_PHR1A_201909131530380_FR1_PX_W062N56_0918_01258/TPP1601254289/IMG_PHR1A_MS_004 -out /data/Labrador/Traitements/Dossiers_images/OpticalCalibration_Marit.tif -level toa -acqui.minute 30 -acqui.hour 15 -acqui.day 13 -acqui.month 9 -acqui.year 2019 -acqui.sun.elev 36.7725 -acqui.sun.azim 170.305 -acqui.view.elev 75.3594 -acqui.view.azim 356.039 -atmo.aerosol maritime -atmo.pixsize 0.008255

You should have more log to share.

Also, I noticed you set the parameter level to toa, but in your original question you were talking about TOC correction. (but I don’t think it is the source of your problem)

Hi,

I personally don’t see the error, I suppose that the values of dates and angles calculated by OTB are not the good ones, but I don’t know how to check that.

Here is the log from the command line, the result is the same :

2020-09-10 17:26:01 (INFO) OpticalCalibration:
File: /data/Labrador/Traitements/Dossiers_images/zip/2018-006-Sci-P004-Univ-Rouen-C.Colange_SO18010794-6-01_DS_PHR1A_201909131530380_FR1_PX_W062N56_0918_01258/TPP1601254289/IMG_PHR1A_MS_004/IMG_PHR1A_MS_201909131530380_SEN_4527881101-004_R1C1.TIF
Sensor detected: PHR 1A
Parameters extract from input image:
Acquisition Day: 13
Acquisition Month: 9
Acquisition Year: 2019
Acquisition Sun Elevation Angle: 36.7725
Acquisition Sun Azimuth Angle: 170.305
Acquisition Viewing Elevation Angle: 75.3594
Acquisition Viewing Azimuth Angle: 356.039
Acquisition gain (per band): 10.46 9.34 9.38 15.69
Acquisition bias (per band): 2.4056 3.57564 3.35568 1.1145
Solar Irradiance (per band): 1594 1830 1915 1060

2020-09-10 17:26:01 (INFO) OpticalCalibration: Default RAM limit for OTB is 256 MB
2020-09-10 17:26:01 (INFO) OpticalCalibration: GDAL maximum cache size is 6446 MB
2020-09-10 17:26:01 (INFO) OpticalCalibration: OTB will use at most 36 threads
2020-09-10 17:26:01 (INFO) OpticalCalibration: Compute Top of Canopy reflectance
2020-09-10 17:26:02 (INFO) OpticalCalibration: Atmospheric correction parameters compute by 6S :
AtmosphericCorrectionParameters (0x1922480)
Atmospheric pressure : 1030
Water vapor amount : 2.5
Ozone amount : 0
Aerosol model : 0
Aerosol optical : 0.2

AtmosphericRadiativeTerms (0x1923260)
4 channels:
Channel 0 :
Intrinsic Atmospheric Reflectance : 0.0202705
Shperical Albedo of the Atmosphere : 0.0470062
Total Gaseous Transmission : 0.974899
Atmosphere Downward Transmittance : 0.954343
Atmosphere Upward Transmittance : 0.970974
Upward Diffuse Transmittance : 0.0248838
Upward Direct Transmittance : 0.94609
Rayleigh Upward Diffuse Transmittance : 0.0269163
Aerosols Upward Diffuse Transmittance : 0
Wavelength : 0.6475
Channel 1 :
Intrinsic Atmospheric Reflectance : 0.0366999
Shperical Albedo of the Atmosphere : 0.0813503
Total Gaseous Transmission : 0.98512
Atmosphere Downward Transmittance : 0.925331
Atmosphere Upward Transmittance : 0.952411
Upward Diffuse Transmittance : 0.0470214
Upward Direct Transmittance : 0.905389
Rayleigh Upward Diffuse Transmittance : 0.047068
Aerosols Upward Diffuse Transmittance : 0
Wavelength : 0.555
Channel 2 :
Intrinsic Atmospheric Reflectance : 0.0599408
Shperical Albedo of the Atmosphere : 0.113361
Total Gaseous Transmission : 0.997648
Atmosphere Downward Transmittance : 0.852682
Atmosphere Upward Transmittance : 0.902662
Upward Diffuse Transmittance : 0.0532171
Upward Direct Transmittance : 0.849445
Rayleigh Upward Diffuse Transmittance : 0.0750858
Aerosols Upward Diffuse Transmittance : 0
Wavelength : 0.495
Channel 3 :
Intrinsic Atmospheric Reflectance : 0.00705798
Shperical Albedo of the Atmosphere : 0.0181052
Total Gaseous Transmission : 0.887334
Atmosphere Downward Transmittance : 0.983855
Atmosphere Upward Transmittance : 0.989931
Upward Diffuse Transmittance : 0.00937074
Upward Direct Transmittance : 0.98056
Rayleigh Upward Diffuse Transmittance : 0.0094555
Aerosols Upward Diffuse Transmittance : 0
Wavelength : 0.8325

2020-09-10 17:26:02 (INFO) OpticalCalibration: Clamp values between [0, 100]
2020-09-10 17:26:02 (INFO): Estimated memory for full processing: 13675.1MB (avail.: 256 MB), optimal image partitioning: 54 blocks
2020-09-10 17:26:02 (INFO): File /data/Labrador/Traitements/Dossiers_images/OpticalCalibration_Marit.tif will be written in 28 blocks of 9938x256 pixels
Writing /data/Labrador/Traitements/Dossiers_images/OpticalCalibration_Marit.tif

I don’t see why your processing results in a black image. Did you check the dynamic of the output image?

Just to be sure the parameters are correct, you can try to run the program with the command line and without the parameters, and see what the program finds :
./otbcli_OpticalCalibration -in /data/Labrador/Traitements/Dossiers_images/zip/2018-006-Sci-P004-Univ-Rouen-C.Colange_SO18010794-6-01_DS_PHR1A_201909131530380_FR1_PX_W062N56_0918_01258/TPP1601254289/IMG_PHR1A_MS_004 -out /data/Labrador/Traitements/Dossiers_images/OpticalCalibration_Marit.tif -level toa