From 7433ebb05b7c0f2c17be59b7420449fe14144a58 Mon Sep 17 00:00:00 2001
From: Jonathan Guinet <jonathan.guinet@c-s.fr>
Date: Tue, 2 Aug 2011 08:45:08 +0200
Subject: [PATCH] DOC: LAI From NDVI code example.
---
Examples/Simulation/CMakeLists.txt | 21 +-
.../Simulation/LAIFromNDVIImageTransform.cxx | 199 ++++++++++++++++++
.../Simulation/otbSimulationExamplesTests.cxx | 9 +
3 files changed, 228 insertions(+), 1 deletion(-)
create mode 100644 Examples/Simulation/LAIFromNDVIImageTransform.cxx
diff --git a/Examples/Simulation/CMakeLists.txt b/Examples/Simulation/CMakeLists.txt
index c4c8a6f318..7ea96f0f7e 100644
--- a/Examples/Simulation/CMakeLists.txt
+++ b/Examples/Simulation/CMakeLists.txt
@@ -9,9 +9,15 @@ ADD_EXECUTABLE(ProsailModel ProsailModel.cxx )
TARGET_LINK_LIBRARIES(ProsailModel OTBSimulation)
+ADD_EXECUTABLE(LAIFromNDVIImageTransform LAIFromNDVIImageTransform.cxx )
+TARGET_LINK_LIBRARIES(LAIFromNDVIImageTransform OTBSimulation)
+
+
IF( NOT OTB_DISABLE_CXX_TESTING AND BUILD_TESTING )
SET(BASELINE ${OTB_DATA_ROOT}/Baseline/Examples/Simulation)
+SET(BASELINE_IMAGES_OTB ${OTB_DATA_ROOT}/Baseline/OTB/Images)
+
SET(INPUTDATA ${OTB_DATA_ROOT}/Examples)
SET(TEMP ${OTB_BINARY_DIR}/Testing/Temporary)
@@ -23,7 +29,7 @@ SET(TOL_2 0.001)
# tests#
-# ------- DEMToRainbowExampleTest ----------
+# ------- ProsailModelExampleTest ----------
ADD_TEST(siTvProsailModelExampleTest ${EXE_TESTS}
--compare-ascii ${TOL_2} ${BASELINE}/siTvProsailModelExampleTest.txt
@@ -47,6 +53,19 @@ ADD_TEST(siTvProsailModelExampleTest ${EXE_TESTS}
)
+ADD_TEST(siTvLAIFromNDVIImageTransfromExampleTest ${EXE_TESTS}
+ --compare-image ${TOL} ${BASELINE_IMAGES_OTB}/raRAndNIRVegetationIndex_LAIFromNDVIFormosat2_verySmallFSATSWcRcNir.tif
+ ${TEMP}/siTvLAIFromNDVIImageTransfromExampleTest_verySmallFSATSW.tif
+ LAIFromNDVIImageTransform
+ ${INPUTDATA}/verySmallFSATSW.tif
+ ${TEMP}/siTvLAIFromNDVIImageTransfromExampleTest_verySmallFSATSW.tif
+ ${TEMP}/verySmallFSATSW_visu.png
+ ${TEMP}/siTvLAIFromNDVIImageTransfromExampleTest_verySmallFSATSW_visu.png
+ 1 # red
+ 4 # nir
+)
+
+
INCLUDE_DIRECTORIES(${OTB_SOURCE_DIR}/Testing/Code)
ADD_EXECUTABLE(otbSimulationExamplesTests otbSimulationExamplesTests.cxx)
TARGET_LINK_LIBRARIES(otbSimulationExamplesTests OTBSimulation OTBCommon OTBTesting)
diff --git a/Examples/Simulation/LAIFromNDVIImageTransform.cxx b/Examples/Simulation/LAIFromNDVIImageTransform.cxx
new file mode 100644
index 0000000000..e7db9ef4fd
--- /dev/null
+++ b/Examples/Simulation/LAIFromNDVIImageTransform.cxx
@@ -0,0 +1,199 @@
+/*=========================================================================
+
+ 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.
+
+=========================================================================*/
+
+// Software Guide : BeginCommandLineArgs
+//
+// INPUTS: {verySmallFSATSW.tif}
+// OUTPUTS: {siTvLAIFromNDVIImageTransfromExampleTest_verySmallFSATSW.tif},{verySmallFSATSW_visu.png},{siTvLAIFromNDVIImageTransfromExampleTest_verySmallFSATSW_visu.png}
+// 1 4
+// Software Guide : EndCommandLineArgs
+//
+// Software Guide : BeginLatex
+//
+// This example presents a way to generate LAI image using formula dedicated to Formosat2.
+// LAI Image is used as an input in Image Simulation process.
+//
+// Let's look at the minimal code required to use this algorithm. First, the
+// following headers must be included.
+//
+// Software Guide : EndLatex
+
+// Software Guide : BeginCodeSnippet
+#include "otbMultiChannelRAndNIRIndexImageFilter.h"
+#include "otbVegetationIndicesFunctor.h"
+// Software Guide : EndCodeSnippet
+
+#include "otbCommandLineArgumentParser.h"
+#include "otbImage.h"
+#include "otbImageFileWriter.h"
+
+#include "otbImage.h"
+#include "otbVectorImage.h"
+#include "otbImageFileReader.h"
+#include "otbImageFileWriter.h"
+
+#include "itkRescaleIntensityImageFilter.h"
+
+int main(int argc, char *argv[])
+{
+ if (argc != 7)
+ {
+ std::cerr << "Wrong Parameters " << std::endl;
+ return EXIT_FAILURE;
+ }
+
+ const unsigned int Dimension = 2;
+ typedef otb::VectorImage<double, Dimension> InputImageType;
+ typedef otb::Image<double, Dimension> OutputImageType;
+ typedef otb::Image<unsigned char, Dimension> ImageVisuType;
+ typedef otb::ImageFileReader<InputImageType> ReaderType;
+ typedef otb::ImageFileWriter<OutputImageType> WriterType;
+ typedef otb::ImageFileWriter<ImageVisuType> VisuWriterType;
+ typedef otb::ImageFileWriter<InputImageType> InWriterType;
+ // Software Guide : BeginLatex
+ //
+ // Filter type is a generic \doxygen{otb}{MultiChannelRAndNIRIndexImageFilter} using Formosat2 specific LAI
+ // \doxygen{otb}{LAIFromNDVIFormosat2Functor}.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ typedef otb::Functor::LAIFromNDVIFormosat2Functor
+ <InputImageType::InternalPixelType,
+ InputImageType::InternalPixelType, OutputImageType::PixelType> FunctorType;
+ typedef otb::MultiChannelRAndNIRIndexImageFilter
+ <InputImageType, OutputImageType, FunctorType>
+ MultiChannelRAndNIRIndexImageFilterType;
+ // Software Guide : EndCodeSnippet
+
+ // Instantiating object
+
+ // Software Guide : BeginLatex
+ //
+ // Next the filter is created by invoking the \code{New()}~method and
+ // assigning the result to a \doxygen{itk}{SmartPointer}.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ MultiChannelRAndNIRIndexImageFilterType::Pointer filter
+ = MultiChannelRAndNIRIndexImageFilterType::New();
+ // Software Guide : EndCodeSnippet
+
+ ReaderType::Pointer reader = ReaderType::New();
+ WriterType::Pointer writer = WriterType::New();
+ InWriterType::Pointer inWriter = InWriterType::New();
+ VisuWriterType::Pointer visuWriter = VisuWriterType::New();
+
+ char * InputName = argv[1];
+ char * OutputName1 = argv[2];
+ char * OutputName2 = argv[3];
+ char * OutputName3 = argv[4];
+ reader->SetFileName(InputName);
+
+
+ // Software Guide : BeginLatex
+ //
+ // filter input is set with input image
+ //
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ filter->SetInput(reader->GetOutput());
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // then red and nir channels index are set using \code{SetRedIndex()} and \code{SetNIRIndex()}
+ //
+ //
+ // Software Guide : EndLatex
+ // Software Guide : BeginCodeSnippet
+ unsigned int redChannel = static_cast<unsigned int> (atoi(argv[5]));
+ unsigned int nirChannel = static_cast<unsigned int> (atoi(argv[6]));
+ filter->SetRedIndex(redChannel);
+ filter->SetNIRIndex(nirChannel);
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ // The invocation of the \code{Update()} method triggers the
+ // execution of the pipeline.
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ filter->Update();
+ // Software Guide : EndCodeSnippet
+
+
+
+ writer->SetFileName(OutputName1);
+ writer->SetInput(filter->GetOutput());
+ writer->Update();
+
+ //rescale data
+
+ inWriter->SetFileName(OutputName2);
+ //typedef itk::RescaleIntensityImageFilter<InputImageType,
+ // OutputImageType> RescalerType;
+
+ //RescalerType::Pointer rescaler = RescalerType::New();
+ //rescaler->SetInput(reader->GetOutput());
+ //rescaler->SetOutputMinimum(0);
+ //rescaler->SetOutputMaximum(255);
+
+ inWriter->SetInput(reader->GetOutput());
+ inWriter->Update();
+
+
+ visuWriter->SetFileName(OutputName3);
+ typedef itk::RescaleIntensityImageFilter<OutputImageType,
+ ImageVisuType> RescalerTypeOut;
+
+ RescalerTypeOut::Pointer rescaler = RescalerTypeOut::New();
+ rescaler->SetInput(filter->GetOutput());
+ rescaler->SetOutputMinimum(0);
+ rescaler->SetOutputMaximum(255);
+
+ visuWriter->SetInput(rescaler->GetOutput());
+ visuWriter->Update();
+
+ // Software Guide : BeginLatex
+ // \begin{figure}
+ // \center
+ // \includegraphics[width=0.44\textwidth]{verySmallFSATSW_visu.eps}
+ // \includegraphics[width=0.44\textwidth]{siTvLAIFromNDVIImageTransfromExampleTest_verySmallFSATSW_visu.eps}
+ // \itkcaption[LAIFromNDVIImageTransform Filter]{LAI generation \emph{(right)} from NDVI applied on Formosat 2 Image \emph{(left)} .}
+ // \label{fig:LAIFromNDVIImageTransform}
+ // \end{figure}
+ //
+ // Figure \ref{fig:LAIFromNDVIImageTransform} illustrates the LAI generation using Formosat 2 data.
+ //
+ //
+ // Software Guide : EndLatex
+
+
+
+
+ return EXIT_SUCCESS;
+}
+
+
+
+
+
diff --git a/Examples/Simulation/otbSimulationExamplesTests.cxx b/Examples/Simulation/otbSimulationExamplesTests.cxx
index b27b524277..b77609a631 100644
--- a/Examples/Simulation/otbSimulationExamplesTests.cxx
+++ b/Examples/Simulation/otbSimulationExamplesTests.cxx
@@ -24,9 +24,18 @@
void RegisterTests()
{
REGISTER_TEST(ProsailModel);
+ REGISTER_TEST(LAIFromNDVIImageTransform);
}
#undef main
#define main ProsailModel
#include "ProsailModel.cxx"
+
+#undef main
+#define main LAIFromNDVIImageTransform
+#include "LAIFromNDVIImageTransform.cxx"
+
+
+
+
--
GitLab