Skip to content
GitLab
Select Git revision
  • master default
  • enhancements
Find file BlameHistoryPermalink
otbMeanResampleImageFilter.hxx 3.86 KiB
1 
/*=========================================================================
2
3 
  Copyright (c) Remi Cresson (IRSTEA). All rights reserved.
4
5
6 
     This software is distributed WITHOUT ANY WARRANTY; without even
7 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
8 
     PURPOSE.  See the above copyright notices for more information.
9
10 
=========================================================================*/
11 
#ifndef __MeanResampleImageFilter_hxx
12 
#define __MeanResampleImageFilter_hxx
13
14 
#include <otbMeanResampleImageFilter.h>
15 
#include "itkProgressReporter.h"
16
17 
namespace otb
18 
{
19 
/**
20 
 *
21 
 */
22 
template <class TImage>
23 
MeanResampleImageFilter<TImage>
24 
::MeanResampleImageFilter()
25 
 {
26 
  m_StepX = 1;
27 
  m_StepY = 1;
28 
  m_NoDataValue = 0;
29 
 }
30
31 
template <class TImage>
32 
void
33 
MeanResampleImageFilter<TImage>
34 
::GenerateOutputInformation()
35 
 {
36 
  Superclass::GenerateOutputInformation();
37
38 
  // Grab input image
39 
  ImageType * inputImage = static_cast<ImageType * >(
40 
      Superclass::ProcessObject::GetInput(0) );
41
42 
  ImageType * outputPtr = this->GetOutput();
43
44 
  // The new output image has the same origin
45 
  ImagePointType origin = inputImage->GetOrigin();
46 
  origin[0] += 0.5 * inputImage->GetSignedSpacing()[0] * (m_StepX - 1);
47 
  origin[1] += 0.5 * inputImage->GetSignedSpacing()[1] * (m_StepY - 1);
48 
  outputPtr->SetOrigin ( origin );
49
50 
  // New spacing for the output image
51 
  ImageSpacingType spacing = inputImage->GetSignedSpacing();
52 
  spacing[0] *= m_StepX;
53 
  spacing[1] *= m_StepY;
54 
  outputPtr->SetSignedSpacing (spacing);
55
56 
  // New size for the output image
57 
  ImageRegionType inRegion = inputImage->GetLargestPossibleRegion();
58 
  ImageRegionType outRegion;
59 
  outRegion.SetIndex(0, 0);
60 
  outRegion.SetIndex(1, 0);
61 
  outRegion.SetSize (0, inRegion.GetSize()[0] / m_StepX);
62 
  outRegion.SetSize (1, inRegion.GetSize()[1] / m_StepY);
63 
  outputPtr->SetLargestPossibleRegion( outRegion );
64
65 
 }
66
67 
template <class TImage>
68 
void
69 
MeanResampleImageFilter<TImage>
70 
::GenerateInputRequestedRegion()
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
{ // Output requested region const ImageRegionType outRegion = this->GetOutput()->GetRequestedRegion(); // Grab input image ImageType * inputImage = static_cast<ImageType * >( Superclass::ProcessObject::GetInput(0) ); ImageRegionType inRegion; inRegion.SetIndex(0, outRegion.GetIndex()[0] * m_StepX); inRegion.SetIndex(1, outRegion.GetIndex()[1] * m_StepY); inRegion.SetSize (0, outRegion.GetSize()[0] * m_StepX); inRegion.SetSize (1, outRegion.GetSize()[1] * m_StepY); inRegion.Crop(inputImage->GetLargestPossibleRegion()); inputImage->SetRequestedRegion(inRegion); } /** * */ template <class TImage> void MeanResampleImageFilter<TImage> ::ThreadedGenerateData(const ImageRegionType& outputRegionForThread, itk::ThreadIdType threadId) { // Support progress methods/callbacks itk::ProgressReporter progress(this, threadId, outputRegionForThread.GetNumberOfPixels() ); // Iterate through the thread region OutputImageIteratorType outputIt(this->GetOutput(), outputRegionForThread); // Grab input image ImageType * inputImage = static_cast<ImageType * >( Superclass::ProcessObject::GetInput(0) ); for ( outputIt.GoToBegin(); !outputIt.IsAtEnd(); ++outputIt) { // sum float accum = 0.0; float npix = 0.0; for (unsigned int x = 0; x < m_StepX ; x++) for (unsigned int y = 0; y < m_StepY ; y++) { ImageIndexType index = outputIt.GetIndex(); index[0] *= m_StepX; index[1] *= m_StepY; index[0] += x; index[1] += y; if (inputImage->GetLargestPossibleRegion().IsInside(index)) { float pixVal = inputImage->GetPixel(index); if (pixVal != m_NoDataValue) { accum += pixVal; npix += 1.0; } } } // normalize if (npix > 0.0) accum /= npix; outputIt.Set(accum); progress.CompletedPixel();
141142143144145146147148
} // Next pixel } } #endif