Orfeo Toolbox  4.0
Fusion/PanSharpeningExample.cxx
/*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.txt for details.
Some parts of this code are derived from ITK. See ITKCopyright.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: {QB_Toulouse_Ortho_PAN.tif}, {QB_Toulouse_Ortho_XS.tif}
// OUTPUTS: {QB_Toulouse_Ortho_PXS.tif}
// OUTPUTS: {pretty_QB_Toulouse_Ortho_PXS.png}
// OUTPUTS: {pretty_QB_Toulouse_Ortho_PAN.png}
// OUTPUTS: {pretty_QB_Toulouse_Ortho_XS.png}
// Software Guide : EndCommandLineArgs
// Software Guide : BeginLatex
//
// Here we are illustrating the use of the
// \doxygen{otb}{SimpleRcsPanSharpeningFusionImageFilter} for PAN-sharpening.
// This example takes a PAN and the corresponding XS images as input. These
// images are supposed to be registered.
//
// The fusion operation is defined as
//
// \begin{equation}
// \frac{XS}{\mathrm{Filtered}(PAN)} PAN
// \end{equation}
//
// Software Guide : EndLatex
// Software Guide : BeginLatex
//
// Figure~\ref{fig:PANSHARP_FILTER} shows the result of applying
// this PAN sharpening filter to a Quickbird image.
// \begin{figure}
// \center
// \includegraphics[width=0.44\textwidth]{pretty_QB_Toulouse_Ortho_PAN.eps}
// \includegraphics[width=0.44\textwidth]{pretty_QB_Toulouse_Ortho_XS.eps}
// \includegraphics[width=0.44\textwidth]{pretty_QB_Toulouse_Ortho_PXS.eps}
// \itkcaption[Pan sharpening]{Result of applying
// the \doxygen{otb}{SimpleRcsPanSharpeningFusionImageFilter} to
// orthorectified Quickbird
// image. From left to right : original PAN image, original XS image and the
// result of the PAN sharpening}
// \label{fig:PANSHARP_FILTER}
// \end{figure}
//
// Software Guide : EndLatex
// Software Guide : BeginLatex
//
// We start by including the required header and declaring the main function:
//
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
#include "otbImage.h"
#include "otbVectorImage.h"
// Software Guide : EndCodeSnippet
// Software Guide : BeginCodeSnippet
int main(int argc, char* argv[])
{
// Software Guide : EndCodeSnippet
if (argc < 7)
{
std::cerr << "Missing Parameters " << std::endl;
std::cerr << "Usage: " << argv[0];
std::cerr <<
" inputPanchromatiqueImage inputMultiSpectralImage outputImage outputImagePrinted panchroPrinted msPrinted"
<< std::endl;
return 1;
}
// Software Guide : BeginLatex
//
// We declare the different image type used here as well as the image reader.
// Note that, the reader for the PAN image is templated by an
// \doxygen{otb}{Image} while the XS reader uses an
// \doxygen{otb}{VectorImage}.
//
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
typedef otb::Image<double, 2> ImageType;
typedef otb::VectorImage<double, 2> VectorImageType;
typedef otb::ImageFileReader<ImageType> ReaderType;
typedef otb::ImageFileReader<VectorImageType> ReaderVectorType;
typedef otb::VectorImage<unsigned short int, 2> VectorIntImageType;
ReaderVectorType::Pointer readerXS = ReaderVectorType::New();
ReaderType::Pointer readerPAN = ReaderType::New();
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
//
// We pass the filenames to the readers
//
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
readerPAN->SetFileName(argv[1]);
readerXS->SetFileName(argv[2]);
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
//
// We declare the fusion filter an set its inputs using the readers:
//
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
<ImageType, VectorImageType, VectorIntImageType> FusionFilterType;
FusionFilterType::Pointer fusion = FusionFilterType::New();
fusion->SetPanInput(readerPAN->GetOutput());
fusion->SetXsInput(readerXS->GetOutput());
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
//
// And finally, we declare the writer and call its \code{Update()} method to
// trigger the full pipeline execution.
//
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
WriterType::Pointer writer = WriterType::New();
writer->SetFileName(argv[3]);
writer->SetInput(fusion->GetOutput());
writer->Update();
// Software Guide : EndCodeSnippet
PrintableImageType::Pointer printable = PrintableImageType::New();
typedef otb::VectorImage<unsigned char,
2>
VectorCharImageType;
PNGWriterType::Pointer pngwriter = PNGWriterType::New();
printable->SetInput(fusion->GetOutput());
printable->SetChannel(3);
printable->SetChannel(2);
printable->SetChannel(1);
pngwriter->SetFileName(argv[4]);
pngwriter->SetInput(printable->GetOutput());
pngwriter->Update();
typedef otb::PrintableImageFilter<VectorImageType> PrintableImageType2;
PrintableImageType2::Pointer printable2 = PrintableImageType2::New();
printable2->SetInput(readerXS->GetOutput());
printable2->SetChannel(3);
printable2->SetChannel(2);
printable2->SetChannel(1);
pngwriter->SetFileName(argv[6]);
pngwriter->SetInput(printable2->GetOutput());
pngwriter->Update();
VectorCastFilterType::Pointer vectorCastFilter = VectorCastFilterType::New();
PNGWriterType::Pointer pngwriterPan = PNGWriterType::New();
vectorCastFilter->SetInput(readerPAN->GetOutput());
PrintableImageType2::Pointer printable3 = PrintableImageType2::New();
printable3->SetInput(vectorCastFilter->GetOutput());
printable3->SetChannel(1);
printable3->SetChannel(1);
printable3->SetChannel(1);
pngwriterPan->SetFileName(argv[5]);
pngwriterPan->SetInput(printable3->GetOutput());
pngwriterPan->Update();
return EXIT_SUCCESS;
}

Generated at Sat Mar 8 2014 14:14:00 for Orfeo Toolbox with doxygen 1.8.3.1