otbAtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms.cxx

Go to the documentation of this file.

/*=========================================================================

  Program:   ORFEO Toolbox
  Language:  C++
  Date:      $Date$
  Version:   $Revision$


  Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
  See OTBCopyright.txt for details.


     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#include "otbAtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms.h"
#include "otbSIXSTraits.h"

namespace otb
{
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms()
{
  this->ProcessObject::SetNumberOfRequiredInputs(1);
  this->ProcessObject::SetNumberOfRequiredOutputs(1);
  // Create the output. We use static_cast<> here because we know the default
  // output must be of type TOutputPointSet
  AtmosphericRadiativeTermsPointer output
    = static_cast<AtmosphericRadiativeTermsType*>(this->MakeOutput(0).GetPointer());

  this->ProcessObject::SetNthOutput(0, output.GetPointer());

}
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms::DataObjectPointer
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::MakeOutput(unsigned int)
{
  return static_cast<itk::DataObject*>(AtmosphericRadiativeTermsType::New().GetPointer());
}

void
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::GraftOutput(itk::DataObject *graft)
{
  this->GraftNthOutput(0, graft);
}

void
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::GraftNthOutput(unsigned int idx, itk::DataObject *graft)
{
  if (idx >= this->GetNumberOfOutputs())
    {
    itkExceptionMacro(<< "Requested to graft output " << idx <<
                      " but this filter only has " << this->GetNumberOfOutputs() << " Outputs.");
    }

  if (!graft)
    {
    itkExceptionMacro(<< "Requested to graft output that is a NULL pointer");
    }

  itk::DataObject * output = this->GetOutput(idx);

  // Call Graft on the PointSet in order to copy meta-information, and containers.
  output->Graft(graft);
}

AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms::AtmosphericRadiativeTermsType *
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::GetOutput(void)
{
  if (this->GetNumberOfOutputs() < 1)
    {
    return 0;
    }
  return static_cast<AtmosphericRadiativeTermsType *> (this->ProcessObject::GetOutput(0));
}

AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms::AtmosphericRadiativeTermsType *
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::GetOutput(unsigned int idx)
{
  return static_cast<AtmosphericRadiativeTermsType*>
         (this->itk::ProcessObject::GetOutput(idx));
}

void
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::SetInput(const AtmosphericCorrectionParametersType *object)
{
  // A single input image
  this->itk::ProcessObject::SetNthInput(0, const_cast<AtmosphericCorrectionParametersType*>(object));
}

AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms::AtmosphericCorrectionParametersType *
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::GetInput(void)
{
  // If there is no input
  if (this->GetNumberOfInputs() != 1)
    {
    // exit
    return 0;
    }
  // else return the first input
  return static_cast<AtmosphericCorrectionParametersType *>
         (this->itk::ProcessObject::GetInput(0));
}

void
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::GenerateData()
{

  AtmosphericCorrectionParametersPointer input = this->GetInput();
  AtmosphericRadiativeTermsPointer       output = this->GetOutput();

  output->GetValues().clear();
  typedef AtmosphericCorrectionParameters::WavelengthSpectralBandVectorType WavelengthSpectralBandVectorType;
  WavelengthSpectralBandVectorType WavelengthSpectralBandVector = input->GetWavelengthSpectralBand();
  unsigned int                     NbBand = WavelengthSpectralBandVector->Size();

  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(
      input->GetSolarZenithalAngle(),                  
      input->GetSolarAzimutalAngle(),                  
      input->GetViewingZenithalAngle(),                
      input->GetViewingAzimutalAngle(),                
      input->GetMonth(),                               
      input->GetDay(),                                 
      input->GetAtmosphericPressure(),                 
      input->GetWaterVaporAmount(),                    
      input->GetOzoneAmount(),                         
      input->GetAerosolModel(),                        
      input->GetAerosolOptical(),                      
      input->GetWavelengthSpectralBand()->GetNthElement(i), 
      atmosphericReflectance,                          
      atmosphericSphericalAlbedo,                      
      totalGaseousTransmission,                        
      downwardTransmittance,                           
      upwardTransmittance,                             
      upwardDiffuseTransmittance,                      
      upwardDirectTransmittance,                       
      upwardDiffuseTransmittanceForRayleigh,           
      upwardDiffuseTransmittanceForAerosol             
      );

    output->SetIntrinsicAtmosphericReflectance(i, atmosphericReflectance);
    output->SetSphericalAlbedo(i, atmosphericSphericalAlbedo);
    output->SetTotalGaseousTransmission(i, totalGaseousTransmission);
    output->SetDownwardTransmittance(i, downwardTransmittance);
    output->SetUpwardTransmittance(i, upwardTransmittance);
    output->SetUpwardDiffuseTransmittance(i, upwardDiffuseTransmittance);
    output->SetUpwardDirectTransmittance(i, upwardDirectTransmittance);
    output->SetUpwardDiffuseTransmittanceForRayleigh(i, upwardDiffuseTransmittanceForRayleigh);
    output->SetUpwardDiffuseTransmittanceForAerosol(i, upwardDiffuseTransmittanceForAerosol);
    output->SetWavelengthSpectralBand(i, input->GetWavelengthSpectralBand()->GetNthElement(i)->GetCenterSpectralValue());
    }
}

void
AtmosphericCorrectionParametersTo6SAtmosphericRadiativeTerms
::PrintSelf(std::ostream& os, itk::Indent indent) const
{
  Superclass::PrintSelf(os, indent);
}

} // end namespace otb