From 888eeedb05a1353e32a1a7a1721e9adff5e2a147 Mon Sep 17 00:00:00 2001
From: Manuel Grizonnet <manuel.grizonnet@orfeo-toolbox.org>
Date: Thu, 3 Mar 2011 18:43:47 +0100
Subject: [PATCH] ENH:compute toc in optical calibration application (for QB
only for now) WPI
---
Radiometry/otbOpticalCalibration.cxx | 55 +++++++++++++++++++++++++---
1 file changed, 49 insertions(+), 6 deletions(-)
diff --git a/Radiometry/otbOpticalCalibration.cxx b/Radiometry/otbOpticalCalibration.cxx
index f7dd43f916..024900315b 100644
--- a/Radiometry/otbOpticalCalibration.cxx
+++ b/Radiometry/otbOpticalCalibration.cxx
@@ -61,16 +61,16 @@ int OpticalCalibration::Execute(otb::ApplicationOptionsResult* parseResult)
typedef otb::MultiplyByScalarImageFilter<FloatImageType,ImageType> ScaleFilterType;
typedef ReflectanceToSurfaceReflectanceImageFilter<FloatImageType,
FloatImageType> ReflectanceToSurfaceReflectanceImageFilterType;
- typedef otb::PipelineMemoryPrintCalculator MemoryCalculatorType;
+ typedef ReflectanceToSurfaceReflectanceImageFilterType::FilterFunctionValuesType FilterFunctionValuesType;
+ typedef FilterFunctionValuesType::ValuesVectorType ValuesVectorType;
+ typedef AtmosphericCorrectionParameters AtmosphericCorrectionParametersType;
+ typedef AtmosphericCorrectionParametersType::AerosolModelType AerosolModelType;
+ typedef otb::PipelineMemoryPrintCalculator MemoryCalculatorType;
// calibration process
if(parseResult->IsOptionPresent("Level"))
{
- if(parseResult->GetParameterString("Level") == "toc")
- {
- std::cerr<<"Not implemented yet"<<std::endl;
- return EXIT_FAILURE;
- }
+
// Read input image information
ReaderType::Pointer reader=ReaderType::New();
@@ -88,6 +88,49 @@ int OpticalCalibration::Execute(otb::ApplicationOptionsResult* parseResult)
ScaleFilterType::Pointer scaleFilter = ScaleFilterType::New();
scaleFilter->SetInput(luminanceToReflectanceFilter->GetOutput());
scaleFilter->SetCoef(1000.);
+ if(parseResult->GetParameterString("Level") == "toc")
+ {
+ //Declare the class to store atmospheric parameters default parameters for now
+ //TODO implement accessor to those parameters
+
+ //reflectanceToSurfaceReflectanceFilter->
+ AtmosphericCorrectionParametersType::Pointer atmosphericParam = reflectanceToSurfaceReflectanceFilter->GetCorrectionParameters();
+ AerosolModelType aeroMod = AtmosphericCorrectionParametersType::NO_AEROSOL;
+
+ atmosphericParam->SetAerosolModel(static_cast<AerosolModelType>(aeroMod));
+ atmosphericParam->SetOzoneAmount(0.);
+ atmosphericParam->SetAtmosphericPressure(1030.);
+ atmosphericParam->SetAerosolOptical(0.2);
+ atmosphericParam->SetWaterVaporAmount(2.5);
+
+ itk::MetaDataDictionary dict = reader->GetOutput()->GetMetaDataDictionary();
+ OpticalImageMetadataInterface::Pointer lImageMetadataInterface = OpticalImageMetadataInterfaceFactory::CreateIMI(dict);
+
+ std::string sensorID = lImageMetadataInterface->GetSensorID();
+ if (sensorID == "QB02")
+ {
+ //Get the filter function values
+ for (unsigned int i = 0; i < reader->GetOutput()->GetNumberOfComponentsPerPixel(); i++)
+ {
+ FilterFunctionValuesType::Pointer functionValues = FilterFunctionValuesType::New();
+
+// functionValues->SetMinSpectralValue(imageMetadataInterface->GetFirstWavelengths()[i]);
+// functionValues->SetMaxSpectralValue(imageMetadataInterface->GetLastWavelengths()[i]);
+// functionValues->SetUserStep(0.0025);
+// //functionValues->
+
+ atmosphericParam->SetWavelengthSpectralBandWithIndex(i, functionValues);
+ }
+ reflectanceToSurfaceReflectanceFilter->SetFilterFunctionCoef(lImageMetadataInterface->GetSpectralSensitivity());
+ }
+
+ reflectanceToSurfaceReflectanceFilter->SetIsSetAtmosphericRadiativeTerms(false);
+ reflectanceToSurfaceReflectanceFilter->SetUseGenerateParameters(true);
+// bProgress->show();
+// Fl::check();
+ reflectanceToSurfaceReflectanceFilter->GenerateParameters();
+ }
+
//Instantiate the writer
WriterType::Pointer writer = WriterType::New();
--
GitLab