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
|
/*=========================================================================
*
* 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 "itkRandomImageSource.h"
#include "itkMeanImageFilter.h"
#include "itkTextOutput.h"
#include "itkTestingMacros.h"
#include "itkVectorImage.h"
int
itkMeanImageFilterTest(int, char *[])
{
using VectorImageType = itk::VectorImage<float, 3>;
itk::MeanImageFilter<VectorImageType, VectorImageType>::New();
// Comment the following if you want to use the itk text output window
itk::OutputWindow::SetInstance(itk::TextOutput::New());
constexpr unsigned int Dimension = 2;
using PixelType = float;
using FloatImage2DType = itk::Image<PixelType, Dimension>;
itk::RandomImageSource<FloatImage2DType>::Pointer random;
random = itk::RandomImageSource<FloatImage2DType>::New();
random->SetMin(0.0);
random->SetMax(1000.0);
FloatImage2DType::SizeValueType randomSize[2];
randomSize[0] = randomSize[1] = 8;
random->SetSize(randomSize);
FloatImage2DType::SpacingValueType spacing[2] = { 0.7, 2.1 };
random->SetSpacing(spacing);
FloatImage2DType::PointValueType origin[2] = { 15, 400 };
random->SetOrigin(origin);
// Create a mean image
itk::MeanImageFilter<FloatImage2DType, FloatImage2DType>::Pointer mean =
itk::MeanImageFilter<FloatImage2DType, FloatImage2DType>::New();
ITK_EXERCISE_BASIC_OBJECT_METHODS(mean, MeanImageFilter, BoxImageFilter);
ITK_TRY_EXPECT_NO_EXCEPTION(random->Update());
mean->SetInput(random->GetOutput());
// Define the neighborhood size used for the mean filter (5x5)
FloatImage2DType::SizeType neighRadius;
neighRadius[0] = 1;
neighRadius[1] = 1;
mean->SetRadius(neighRadius);
ITK_TEST_SET_GET_VALUE(neighRadius, mean->GetRadius());
ITK_TRY_EXPECT_NO_EXCEPTION(mean->Update());
itk::ImageRegionIterator<FloatImage2DType> it;
it = itk::ImageRegionIterator<FloatImage2DType>(random->GetOutput(), random->GetOutput()->GetBufferedRegion());
std::cout << "Input image" << std::endl;
unsigned int i;
for (i = 1; !it.IsAtEnd(); ++i, ++it)
{
std::cout << '\t' << it.Get();
if ((i % 8) == 0)
{
std::cout << std::endl;
}
}
std::cout << "Output image" << std::endl;
it = itk::ImageRegionIterator<FloatImage2DType>(mean->GetOutput(), mean->GetOutput()->GetBufferedRegion());
for (i = 1; !it.IsAtEnd(); ++i, ++it)
{
std::cout << '\t' << it.Get();
if ((i % 8) == 0)
{
std::cout << std::endl;
}
}
return EXIT_SUCCESS;
}
|