From 21e12102a703c453b6128a4a58ed06d1b2eb78ae Mon Sep 17 00:00:00 2001
From: Manuel Grizonnet <manuel.grizonnet@orfeo-toolbox.org>
Date: Mon, 9 Jul 2012 16:33:06 +0200
Subject: [PATCH] DOC: add output DEM estimation figure in stereo
reconstruction example
---
Examples/DisparityMap/CMakeLists.txt | 1 +
.../StereoReconstructionExample.cxx | 44 ++++++++++++++++---
2 files changed, 40 insertions(+), 5 deletions(-)
diff --git a/Examples/DisparityMap/CMakeLists.txt b/Examples/DisparityMap/CMakeLists.txt
index 0a4b325bd2..b566db59d1 100644
--- a/Examples/DisparityMap/CMakeLists.txt
+++ b/Examples/DisparityMap/CMakeLists.txt
@@ -52,6 +52,7 @@ ADD_TEST(dmTeStereoReconstructionExampleTest ${EXE_TESTS}
${INPUTDATA}/sensor_stereo_left.tif
${INPUTDATA}/sensor_stereo_right.tif
${TEMP}/elevationOutput.tif
+ ${TEMP}/elevationOutputPrintable.png
140
)
diff --git a/Examples/DisparityMap/StereoReconstructionExample.cxx b/Examples/DisparityMap/StereoReconstructionExample.cxx
index 01602eafac..2c7b514f9c 100644
--- a/Examples/DisparityMap/StereoReconstructionExample.cxx
+++ b/Examples/DisparityMap/StereoReconstructionExample.cxx
@@ -19,7 +19,7 @@
// Software Guide : BeginCommandLineArgs
// INPUTS: {sensor_stereo_left.tif}, {sensor_stereo_right.tif}
-// OUTPUTS: {elevationOutput.tif}
+// OUTPUTS: {elevationOutput.tif},{elevationOutputPrintable.png}
// 140
// Software Guide : EndCommandLineArgs
@@ -58,15 +58,16 @@
#include "itkVectorCastImageFilter.h"
#include "otbImageList.h"
#include "otbImageListToVectorImageFilter.h"
+#include "itkRescaleIntensityImageFilter.h"
// Software Guide : EndCodeSnippet
int main(int argc, char* argv[])
{
- if (argc != 5)
+ if (argc != 6)
{
std::cerr << "Usage: " << argv[0];
- std::cerr << " sensorImage1 sensorImage2 outputDEM ";
+ std::cerr << " sensorImage1 sensorImage2 outputDEM outputDEMPNG";
std::cerr << "averageElevation " << std::endl;
return EXIT_FAILURE;
}
@@ -86,6 +87,14 @@ int main(int argc, char* argv[])
typedef otb::StreamingImageFileWriter
<FloatImageType> WriterType;
+ typedef unsigned char OutputPixelType;
+ typedef otb::Image<OutputPixelType, 2> OutputImageType;
+
+ typedef itk::RescaleIntensityImageFilter<FloatImageType,
+ OutputImageType> RescalerType;
+
+ typedef otb::StreamingImageFileWriter
+ <OutputImageType> OutputWriterType;
// Software Guide : BeginLatex
// This example demonstrates the use of the following filters :
// \begin{itemize}
@@ -150,7 +159,7 @@ int main(int argc, char* argv[])
FloatImageType> DisparityToElevationFilterType;
// Software Guide : EndCodeSnippet
- double avgElevation = atof(argv[4]);
+ double avgElevation = atof(argv[5]);
ImageReaderType::Pointer leftReader = ImageReaderType::New();
ImageReaderType::Pointer rightReader = ImageReaderType::New();
@@ -382,7 +391,32 @@ int main(int argc, char* argv[])
m_DEMWriter->SetInput(m_DispToElev->GetOutput());
m_DEMWriter->SetFileName(argv[3]);
m_DEMWriter->Update();
-// Software Guide : EndCodeSnippet
+
+ RescalerType::Pointer fieldRescaler = RescalerType::New();
+ fieldRescaler->SetInput(m_DispToElev->GetOutput());
+ fieldRescaler->SetOutputMaximum(255);
+ fieldRescaler->SetOutputMinimum(0);
+
+ OutputWriterType::Pointer fieldWriter = OutputWriterType::New();
+ fieldWriter->SetInput(fieldRescaler->GetOutput());
+ fieldWriter->SetFileName(argv[4]);
+ fieldWriter->Update();
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // Figure~\ref{fig:STEREORECONSTRUCTIONOUTPUT} shows the result of applying
+ // terrain reconstruction based using pixel-wise block matching, sub-pixel
+ // interpolation and DEM estimation using a pair of Pleiades images over the
+ // \textit{Stadium} in Toulouse, France.
+ // \begin{figure}
+ // \center
+ // \includegraphics[width=0.4\textwidth]{elevationOutputPrintable.eps}
+ // \itkcaption[From stereo pair to elevation]{DEM image estimated from the disparity.}
+ // \label{fig:STEREORECONSTRUCTIONOUTPUT}
+ // \end{figure}
+ //
+ // Software Guide : EndLatex
return EXIT_SUCCESS;
}
--
GitLab