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
|
/*=========================================================================
*
* 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.
*
*=========================================================================*/
#include "itkScalarImageToHistogramGenerator.h"
#include "itkMinimumMaximumImageFilter.h"
#include "itkImageFileReader.h"
#include "itkImageFileWriter.h"
#include "itkHistogramToEntropyImageFilter.h"
#include "itkTestingMacros.h"
int
itkHistogramToEntropyImageFilterTest1(int argc, char * argv[])
{
if (argc < 3)
{
std::cerr << "Missing parameters." << std::endl;
std::cerr << "Usage: " << itkNameOfTestExecutableMacro(argv);
std::cerr << " inputScalarImageFileName outputImage" << std::endl;
return EXIT_FAILURE;
}
constexpr unsigned int Dimension = 2;
using PixelComponentType = unsigned char;
using ScalarImageType = itk::Image<PixelComponentType, Dimension>;
using ReaderType = itk::ImageFileReader<ScalarImageType>;
auto reader = ReaderType::New();
reader->SetFileName(argv[1]);
ITK_TRY_EXPECT_NO_EXCEPTION(reader->Update());
itk::MinimumMaximumImageFilter<ScalarImageType>::Pointer minmaxFilter =
itk::MinimumMaximumImageFilter<ScalarImageType>::New();
minmaxFilter->SetInput(reader->GetOutput());
minmaxFilter->Update();
const ScalarImageType::PixelType imageMin = minmaxFilter->GetMinimum();
const ScalarImageType::PixelType imageMax = minmaxFilter->GetMaximum();
using HistogramGeneratorType = itk::Statistics::ScalarImageToHistogramGenerator<ScalarImageType>;
auto histogramGenerator = HistogramGeneratorType::New();
histogramGenerator->SetInput(reader->GetOutput());
const int NumberOfBins = static_cast<unsigned int>(imageMax - imageMin + 1);
histogramGenerator->SetNumberOfBins(NumberOfBins);
histogramGenerator->SetMarginalScale(1.0);
histogramGenerator->SetHistogramMin(imageMin);
histogramGenerator->SetHistogramMax(imageMax);
histogramGenerator->Compute();
using HistogramType = HistogramGeneratorType::HistogramType;
const HistogramType * histogram = histogramGenerator->GetOutput();
using HistogramToImageFilterType = itk::HistogramToEntropyImageFilter<HistogramType>;
auto histogramToImageFilter = HistogramToImageFilterType::New();
ITK_EXERCISE_BASIC_OBJECT_METHODS(histogramToImageFilter, HistogramToEntropyImageFilter, HistogramToImageFilter);
histogramToImageFilter->SetInput(histogram);
using OutputImageType = HistogramToImageFilterType::OutputImageType;
using WriterType = itk::ImageFileWriter<OutputImageType>;
auto writer = WriterType::New();
writer->SetFileName(argv[2]);
writer->SetInput(histogramToImageFilter->GetOutput());
ITK_TRY_EXPECT_NO_EXCEPTION(writer->Update());
return EXIT_SUCCESS;
}
|