mitkHistogramGenerator.cpp
Go to the documentation of this file.
00001 /*========================================================================= 00002 00003 Program: Medical Imaging & Interaction Toolkit 00004 Language: C++ 00005 Date: $Date$ 00006 Version: $Revision$ 00007 00008 Copyright (c) German Cancer Research Center, Division of Medical and 00009 Biological Informatics. All rights reserved. 00010 See MITKCopyright.txt or https://www.mitk.org/copyright.html for details. 00011 00012 This software is distributed WITHOUT ANY WARRANTY; without even 00013 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 00014 PURPOSE. See the above copyright notices for more information. 00015 00016 =========================================================================*/ 00017 00018 #if(_MSC_VER==1200) 00019 #include <itkFixedCenterOfRotationAffineTransform.h> 00020 #endif 00021 00022 #include "mitkHistogramGenerator.h" 00023 #include "mitkImageTimeSelector.h" 00024 #include "mitkImageAccessByItk.h" 00025 00026 00027 // 00028 // The new ITK Statistics framework has 00029 // a class with the same functionality as 00030 // MITKScalarImageToHistogramGenerator.h, but 00031 // no longer has the classis the MITK class depends on. 00032 #if !defined(ITK_USE_REVIEW_STATISTICS) 00033 #include "itkMITKScalarImageToHistogramGenerator.h" 00034 #else 00035 #include "itkScalarImageToHistogramGenerator.h" 00036 #endif 00037 00038 mitk::HistogramGenerator::HistogramGenerator() : m_Image(NULL), m_Size(256), m_Histogram(NULL) 00039 { 00040 00041 } 00042 00043 mitk::HistogramGenerator::~HistogramGenerator() 00044 { 00045 } 00046 00047 template < typename TPixel, unsigned int VImageDimension > 00048 void InternalCompute(itk::Image< TPixel, VImageDimension >* itkImage, const mitk::HistogramGenerator* mitkHistoGenerator, mitk::HistogramGenerator::HistogramType::ConstPointer& histogram) 00049 { 00050 00051 #if !defined(ITK_USE_REVIEW_STATISTICS) 00052 typedef itk::Statistics::MITKScalarImageToHistogramGenerator< 00053 itk::Image< TPixel, VImageDimension >, 00054 double > HistogramGeneratorType; 00055 #else 00056 typedef itk::Statistics::ScalarImageToHistogramGenerator< itk::Image< TPixel, VImageDimension > > 00057 HistogramGeneratorType; 00058 00059 #endif 00060 typename HistogramGeneratorType::Pointer histogramGenerator = HistogramGeneratorType::New(); 00061 00062 histogramGenerator->SetInput( itkImage ); 00063 00064 histogramGenerator->SetNumberOfBins( mitkHistoGenerator->GetSize() ); 00065 // histogramGenerator->SetMarginalScale( 10.0 ); 00066 histogramGenerator->Compute(); 00067 00068 histogram = histogramGenerator->GetOutput(); 00069 } 00070 00071 void mitk::HistogramGenerator::ComputeHistogram() 00072 { 00073 if((m_Histogram.IsNull()) || (m_Histogram->GetMTime() < m_Image->GetMTime())) 00074 { 00075 const_cast<mitk::Image*>(m_Image.GetPointer())->SetRequestedRegionToLargestPossibleRegion(); //@todo without this, Image::GetScalarMin does not work for dim==3 (including sliceselector!) 00076 const_cast<mitk::Image*>(m_Image.GetPointer())->Update(); 00077 mitk::ImageTimeSelector::Pointer timeSelector=mitk::ImageTimeSelector::New(); 00078 timeSelector->SetInput(m_Image); 00079 timeSelector->SetTimeNr( 0 ); 00080 timeSelector->UpdateLargestPossibleRegion(); 00081 AccessByItk_2( timeSelector->GetOutput() , InternalCompute, this, m_Histogram); 00082 } 00083 00084 // debug code 00085 /* 00086 MITK_INFO << "Histogram modfied 1" << m_Histogram->GetMTime() << std::endl; 00087 m_Histogram->Modified(); 00088 MITK_INFO << "Histogram modfied 2" << m_Histogram->GetMTime() << std::endl; 00089 MITK_INFO << "Image modfied" << m_Image->GetMTime() << std::endl; 00090 const unsigned int histogramSize = m_Histogram->Size(); 00091 00092 MITK_INFO << "Histogram size " << histogramSize << std::endl; 00093 00094 HistogramType::ConstIterator itr = GetHistogram()->Begin(); 00095 HistogramType::ConstIterator end = GetHistogram()->End(); 00096 00097 int bin = 0; 00098 while( itr != end ) 00099 { 00100 MITK_INFO << "bin = " << GetHistogram()->GetBinMin(0,bin) << "--" << GetHistogram()->GetBinMax(0,bin) << " frequency = "; 00101 MITK_INFO << itr.GetFrequency() << std::endl; 00102 ++itr; 00103 ++bin; 00104 } 00105 */ 00106 } 00107 00108 float mitk::HistogramGenerator::GetMaximumFrequency() const { 00109 return CalculateMaximumFrequency(this->m_Histogram); 00110 }; 00111 00112 float mitk::HistogramGenerator::CalculateMaximumFrequency(const HistogramType* histogram) 00113 { 00114 HistogramType::ConstIterator itr = histogram->Begin(); 00115 HistogramType::ConstIterator end = histogram->End(); 00116 00117 float maxFreq = 0; 00118 while( itr != end ) 00119 { 00120 maxFreq = vnl_math_max(maxFreq, 00121 // get rid of ambiguity with type signature 00122 // for vnl_math_max 00123 static_cast<float>(itr.GetFrequency())); 00124 ++itr; 00125 } 00126 return maxFreq; 00127 };
