00001 /* 00002 // 00003 // Copyright 2004-2010 SRI International 00004 // 00005 // Copyright 1997-2009 Torsten Rohlfing 00006 // 00007 // This file is part of the Computational Morphometry Toolkit. 00008 // 00009 // http://www.nitrc.org/projects/cmtk/ 00010 // 00011 // The Computational Morphometry Toolkit is free software: you can 00012 // redistribute it and/or modify it under the terms of the GNU General Public 00013 // License as published by the Free Software Foundation, either version 3 of 00014 // the License, or (at your option) any later version. 00015 // 00016 // The Computational Morphometry Toolkit is distributed in the hope that it 00017 // will be useful, but WITHOUT ANY WARRANTY; without even the implied 00018 // warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00019 // GNU General Public License for more details. 00020 // 00021 // You should have received a copy of the GNU General Public License along 00022 // with the Computational Morphometry Toolkit. If not, see 00023 // <http://www.gnu.org/licenses/>. 00024 // 00025 // $Revision: 2022 $ 00026 // 00027 // $LastChangedDate: 2010-07-21 15:26:03 -0700 (Wed, 21 Jul 2010) $ 00028 // 00029 // $LastChangedBy: torstenrohlfing $ 00030 // 00031 */ 00032 00033 #include "cmtkTypedArrayFunctionHistogramEqualization.h" 00034 00035 #include <vector> 00036 #include <algorithm> 00037 00038 cmtk::TypedArrayFunctionHistogramEqualization 00039 ::TypedArrayFunctionHistogramEqualization 00040 ( const TypedArray& variableArray, const size_t numberOfHistogramBins ) 00041 { 00042 this->m_Histogram = Histogram<unsigned int>::SmartPtr( variableArray.GetHistogram( numberOfHistogramBins ) ); 00043 (*this->m_Histogram)[0] = 0; // this effectively stretches the distribution 00044 this->m_Histogram->ConvertToCumulative(); 00045 00046 const Types::DataItemRange range = variableArray.GetRange(); 00047 this->m_MinValue = range.m_LowerBound; 00048 this->m_ScaleFactor = 1.0 * range.Width() / (*this->m_Histogram)[numberOfHistogramBins-1]; 00049 } 00050 00051 cmtk::Types::DataItem 00052 cmtk::TypedArrayFunctionHistogramEqualization 00053 ::operator()( const cmtk::Types::DataItem valueIn ) const 00054 { 00055 return this->m_MinValue + this->m_ScaleFactor * (*this->m_Histogram)[ this->m_Histogram->ValueToBin( valueIn ) ]; 00056 }