File: elxResampleInterpolatorGTest.cxx

package info (click to toggle)
elastix 5.2.0-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 42,480 kB
  • sloc: cpp: 68,403; lisp: 4,118; python: 1,013; xml: 182; sh: 177; makefile: 33
file content (128 lines) | stat: -rw-r--r-- 5,194 bytes parent folder | download
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
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
/*=========================================================================
 *
 *  Copyright UMC Utrecht and contributors
 *
 *  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
 *
 *        http://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.
 *
 *=========================================================================*/

// First include the header files to be tested:
#include "BSplineResampleInterpolator/elxBSplineResampleInterpolator.h"
#include "BSplineResampleInterpolatorFloat/elxBSplineResampleInterpolatorFloat.h"
#include "LinearResampleInterpolator/elxLinearResampleInterpolator.h"
#include "NearestNeighborResampleInterpolator/elxNearestNeighborResampleInterpolator.h"
#include "RayCastResampleInterpolator/elxRayCastResampleInterpolator.h"
#include "RDBSplineResampleInterpolator/elxRDBSplineResampleInterpolator.h" // For ReducedDimensionBSplineResampleInterpolator.

#include "elxElastixTemplate.h"
#include "elxGTestUtilities.h"
#include "../Core/Main/GTesting/elxCoreMainGTestUtilities.h"

// ITK header file:
#include <itkImage.h>

// GoogleTest header file:
#include <gtest/gtest.h>

// Standard C++ header file:
#include <typeinfo>


using ParameterValuesType = itk::ParameterFileParser::ParameterValuesType;
using ParameterMapType = itk::ParameterFileParser::ParameterMapType;


// Using-declarations:
using elx::CoreMainGTestUtilities::CheckNew;
using elx::GTestUtilities::ExpectAllKeysUnique;
using elx::GTestUtilities::MakeMergedMap;


namespace
{

template <unsigned NDimension>
using ElastixType = elx::ElastixTemplate<itk::Image<float, NDimension>, itk::Image<float, NDimension>>;


// All tests specific to a dimension.
template <unsigned NDimension>
struct WithDimension
{
  template <template <typename> class TInterpolatorTemplate>
  struct WithInterpolator
  {
    using InterpolatorType = TInterpolatorTemplate<ElastixType<NDimension>>;

    static void
    Test_CreateTransformParameterMap_for_default_interpolator(const ParameterMapType & expectedDerivedParameterMap)
    {
      SCOPED_TRACE(std::string("Function = ")
                     .append(__func__)
                     .append("\n  InterpolatorType = ")
                     .append(typeid(InterpolatorType).name()));

      const auto                                                     newInterpolator = CheckNew<InterpolatorType>();
      const elx::ResampleInterpolatorBase<ElastixType<NDimension>> & interpolator = *newInterpolator;

      ParameterMapType actualParameterMap;
      interpolator.CreateTransformParameterMap(actualParameterMap);

      const ParameterMapType expectedBaseParameterMap = { { "ResampleInterpolator",
                                                            { interpolator.elxGetClassName() } } };

      ExpectAllKeysUnique(expectedDerivedParameterMap, expectedBaseParameterMap);
      EXPECT_EQ(actualParameterMap, MakeMergedMap(expectedDerivedParameterMap, expectedBaseParameterMap));
    }
  };


  static void
  Test_CreateTransformParameterMap_for_default_interpolator()
  {
    using namespace elx;

    const std::string expectedFinalBSplineInterpolationOrderKey = "FinalBSplineInterpolationOrder";
    const std::string expectedZero = "0";

    WithInterpolator<BSplineResampleInterpolator>::Test_CreateTransformParameterMap_for_default_interpolator(
      { { expectedFinalBSplineInterpolationOrderKey, { "3" } } });
    WithInterpolator<BSplineResampleInterpolatorFloat>::Test_CreateTransformParameterMap_for_default_interpolator(
      { { expectedFinalBSplineInterpolationOrderKey, { "3" } } });
    WithInterpolator<LinearResampleInterpolator>::Test_CreateTransformParameterMap_for_default_interpolator({});
    WithInterpolator<NearestNeighborResampleInterpolator>::Test_CreateTransformParameterMap_for_default_interpolator(
      {});

    const auto skippedTest = [expectedZero] {
      // Note: The following crashes when trying to retrieve "PreParameters" by `m_PreTransform->GetParameters()`.
      WithInterpolator<RayCastResampleInterpolator>::Test_CreateTransformParameterMap_for_default_interpolator(
        { { "FocalPoint", ParameterValuesType(NDimension, expectedZero) },
          { "PreParameters", { expectedZero } },
          { "Threshold", { expectedZero } } });
    };
    (void)skippedTest;

    WithInterpolator<ReducedDimensionBSplineResampleInterpolator>::
      Test_CreateTransformParameterMap_for_default_interpolator(
        { { expectedFinalBSplineInterpolationOrderKey, { "1" } } });
  }
};

} // namespace


GTEST_TEST(ResampleInterpolator, CreateTransformParameterMapForDefaultInterpolator)
{
  WithDimension<2>::Test_CreateTransformParameterMap_for_default_interpolator();
  WithDimension<3>::Test_CreateTransformParameterMap_for_default_interpolator();
}