otbRadiometryCorrectionParametersToAtmosphericRadiativeTerms.h

Go to the documentation of this file.

/*
 * Copyright (C) 2005-2022 Centre National d'Etudes Spatiales (CNES)
 *
 * This file is part of Orfeo Toolbox
 *
 *     https://www.orfeo-toolbox.org/
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
#ifndef otbRadiometryCorrectionParametersToAtmosphericRadiativeTerms_h
#define otbRadiometryCorrectionParametersToAtmosphericRadiativeTerms_h
 
#include "otbAtmosphericRadiativeTerms.h"
#include "otbImageMetadataCorrectionParameters.h"
#include "otbSIXSTraits.h"
#include "otbMacro.h"
 
namespace otb
{
 
class RadiometryCorrectionParametersToAtmosphericRadiativeTerms
{
public:
  static AtmosphericRadiativeTerms::Pointer Compute(AtmosphericCorrectionParameters* paramAtmo, ImageMetadataCorrectionParameters* paramAcqui)
  {
    AtmosphericRadiativeTerms::Pointer radTermsOut = AtmosphericRadiativeTerms::New();
 
    typedef AtmosphericCorrectionParameters::WavelengthSpectralBandVectorType WavelengthSpectralBandVectorType;
    WavelengthSpectralBandVectorType                                          WavelengthSpectralBandVector = paramAcqui->GetWavelengthSpectralBand();
    unsigned int                                                              NbBand                       = WavelengthSpectralBandVector->Size();
 
    radTermsOut->ValuesInitialization(NbBand);
 
    double atmosphericReflectance(0.);
    double atmosphericSphericalAlbedo(0.);
    double totalGaseousTransmission(0.);
    double downwardTransmittance(0.);
    double upwardTransmittance(0.);
    double upwardDiffuseTransmittance(0.);
    double upwardDirectTransmittance(0.);
    double upwardDiffuseTransmittanceForRayleigh(0.);
    double upwardDiffuseTransmittanceForAerosol(0.);
 
    for (unsigned int i = 0; i < NbBand; ++i)
    {
      atmosphericReflectance                = 0.;
      atmosphericSphericalAlbedo            = 0.;
      totalGaseousTransmission              = 0.;
      downwardTransmittance                 = 0.;
      upwardTransmittance                   = 0.;
      upwardDiffuseTransmittance            = 0.;
      upwardDirectTransmittance             = 0.;
      upwardDiffuseTransmittanceForRayleigh = 0.;
      upwardDiffuseTransmittanceForAerosol  = 0.;
      SIXSTraits::ComputeAtmosphericParameters(
          );
 
      radTermsOut->SetIntrinsicAtmosphericReflectance(i, atmosphericReflectance);
      radTermsOut->SetSphericalAlbedo(i, atmosphericSphericalAlbedo);
      radTermsOut->SetTotalGaseousTransmission(i, totalGaseousTransmission);
      radTermsOut->SetDownwardTransmittance(i, downwardTransmittance);
      radTermsOut->SetUpwardTransmittance(i, upwardTransmittance);
      radTermsOut->SetUpwardDiffuseTransmittance(i, upwardDiffuseTransmittance);
      radTermsOut->SetUpwardDirectTransmittance(i, upwardDirectTransmittance);
      radTermsOut->SetUpwardDiffuseTransmittanceForRayleigh(i, upwardDiffuseTransmittanceForRayleigh);
      radTermsOut->SetUpwardDiffuseTransmittanceForAerosol(i, upwardDiffuseTransmittanceForAerosol);
      radTermsOut->SetWavelengthSpectralBand(i, paramAcqui->GetWavelengthSpectralBand()->GetNthElement(i)->GetCenterSpectralValue());
    }
 
    return radTermsOut;
  }
};
} // end namespace otb
#endif