1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
|
/*=========================================================================
*
* Copyright NumFOCUS
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* https://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
/*=========================================================================
*
* Portions of this file are subject to the VTK Toolkit Version 3 copyright.
*
* Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
*
* For complete copyright, license and disclaimer of warranty information
* please refer to the NOTICE file at the top of the ITK source tree.
*
*=========================================================================*/
#ifndef itkRescaleIntensityImageFilter_hxx
#define itkRescaleIntensityImageFilter_hxx
#include "itkMinimumMaximumImageCalculator.h"
#include "itkMath.h"
namespace itk
{
template <typename TInputImage, typename TOutputImage>
RescaleIntensityImageFilter<TInputImage, TOutputImage>::RescaleIntensityImageFilter()
: m_Scale(1.0)
, m_Shift(0.0)
, m_InputMinimum(NumericTraits<InputPixelType>::max())
, m_InputMaximum(InputPixelType{})
, m_OutputMinimum(NumericTraits<OutputPixelType>::NonpositiveMin())
, m_OutputMaximum(NumericTraits<OutputPixelType>::max())
{}
template <typename TInputImage, typename TOutputImage>
void
RescaleIntensityImageFilter<TInputImage, TOutputImage>::BeforeThreadedGenerateData()
{
if (this->m_OutputMinimum > this->m_OutputMaximum)
{
itkExceptionMacro("Minimum output value cannot be greater than Maximum output value.");
}
using CalculatorType = MinimumMaximumImageCalculator<TInputImage>;
auto calculator = CalculatorType::New();
calculator->SetImage(this->GetInput());
calculator->Compute();
this->m_InputMinimum = calculator->GetMinimum();
this->m_InputMaximum = calculator->GetMaximum();
if (itk::Math::NotAlmostEquals(this->m_InputMinimum, this->m_InputMaximum))
{
this->m_Scale = (static_cast<RealType>(this->m_OutputMaximum) - static_cast<RealType>(this->m_OutputMinimum)) /
(static_cast<RealType>(this->m_InputMaximum) - static_cast<RealType>(this->m_InputMinimum));
}
else if (itk::Math::NotAlmostEquals(this->m_InputMaximum, typename NumericTraits<InputPixelType>::ValueType{}))
{
this->m_Scale = (static_cast<RealType>(this->m_OutputMaximum) - static_cast<RealType>(this->m_OutputMinimum)) /
static_cast<RealType>(this->m_InputMaximum);
}
else
{
this->m_Scale = 0.0;
}
this->m_Shift =
static_cast<RealType>(this->m_OutputMinimum) - static_cast<RealType>(this->m_InputMinimum) * this->m_Scale;
// Set up the functor values
this->GetFunctor().SetMinimum(this->m_OutputMinimum);
this->GetFunctor().SetMaximum(this->m_OutputMaximum);
this->GetFunctor().SetFactor(this->m_Scale);
this->GetFunctor().SetOffset(this->m_Shift);
}
template <typename TInputImage, typename TOutputImage>
void
RescaleIntensityImageFilter<TInputImage, TOutputImage>::PrintSelf(std::ostream & os, Indent indent) const
{
Superclass::PrintSelf(os, indent);
os << indent
<< "Output Minimum: " << static_cast<typename NumericTraits<OutputPixelType>::PrintType>(this->m_OutputMinimum)
<< std::endl;
os << indent
<< "Output Maximum: " << static_cast<typename NumericTraits<OutputPixelType>::PrintType>(this->m_OutputMaximum)
<< std::endl;
}
} // end namespace itk
#endif
|