@@ -38,23 +38,23 @@ TestDisplacementJacobianDeterminantValue()
38
38
using VectorImageType = FieldType;
39
39
40
40
41
- std::cout << " Create the dispacementfield image pattern." << std::endl;
41
+ std::cout << " Create the displacementField image pattern." << std::endl;
42
42
VectorImageType::RegionType region;
43
43
// NOTE: Making the image size much larger than necessary in order to get
44
44
// some meaningful time measurements. Simulate a 256x256x256 image.
45
45
constexpr VectorImageType::SizeType size = { { 4096 , 4096 } };
46
46
region.SetSize (size);
47
47
48
- auto dispacementfield = VectorImageType::New ();
49
- dispacementfield ->SetLargestPossibleRegion (region);
50
- dispacementfield ->SetBufferedRegion (region);
51
- dispacementfield ->Allocate ();
48
+ auto displacementField = VectorImageType::New ();
49
+ displacementField ->SetLargestPossibleRegion (region);
50
+ displacementField ->SetBufferedRegion (region);
51
+ displacementField ->Allocate ();
52
52
53
53
VectorType values;
54
54
values[0 ] = 0 ;
55
55
values[1 ] = 0 ;
56
56
using Iterator = itk::ImageRegionIteratorWithIndex<VectorImageType>;
57
- for (Iterator inIter (dispacementfield , region); !inIter.IsAtEnd (); ++inIter)
57
+ for (Iterator inIter (displacementField , region); !inIter.IsAtEnd (); ++inIter)
58
58
{
59
59
const unsigned int i = inIter.GetIndex ()[0 ];
60
60
const unsigned int j = inIter.GetIndex ()[1 ];
@@ -96,7 +96,7 @@ TestDisplacementJacobianDeterminantValue()
96
96
#endif
97
97
ITK_TEST_SET_GET_BOOLEAN (filter, UseImageSpacing, useImageSpacing);
98
98
99
- filter->SetInput (dispacementfield );
99
+ filter->SetInput (displacementField );
100
100
filter->Update ();
101
101
102
102
@@ -127,7 +127,7 @@ TestDisplacementJacobianDeterminantValue()
127
127
128
128
// Define physical point to test
129
129
itk::Point<double , 2 > physPt;
130
- dispacementfield ->TransformIndexToPhysicalPoint (index, physPt);
130
+ displacementField ->TransformIndexToPhysicalPoint (index, physPt);
131
131
132
132
// Use this function to get the determinant at a specified physical point (it will be a different index between tests)
133
133
auto GetDeterminantAtPoint = [](FieldType::Pointer image, const itk::Point<double , 2 > & pt) -> float {
@@ -141,7 +141,7 @@ TestDisplacementJacobianDeterminantValue()
141
141
return filter->GetOutput ()->GetPixel (mappedIdx);
142
142
};
143
143
144
- const float detOriginal = GetDeterminantAtPoint (dispacementfield , physPt);
144
+ const float detOriginal = GetDeterminantAtPoint (displacementField , physPt);
145
145
146
146
// Check that this is the same as above, just to be sure the function above works
147
147
if (itk::Math::abs (detOriginal - expectedJacobianDeterminant) > epsilon)
@@ -162,24 +162,24 @@ TestDisplacementJacobianDeterminantValue()
162
162
order[0 ] = 1 ;
163
163
order[1 ] = 0 ; // swap X <-> Y
164
164
permute->SetOrder (order);
165
- permute->SetInput (dispacementfield );
165
+ permute->SetInput (displacementField );
166
166
permute->Update ();
167
167
168
- auto dispacementfieldPermuted = permute->GetOutput ();
168
+ auto displacementFieldPermuted = permute->GetOutput ();
169
169
170
170
// Adjust direction to match permutation
171
- itk::Matrix<double , 2 , 2 > origDir = dispacementfield ->GetDirection ();
171
+ itk::Matrix<double , 2 , 2 > origDir = displacementField ->GetDirection ();
172
172
itk::Matrix<double , 2 , 2 > newDirPermute;
173
173
newDirPermute[0 ][0 ] = origDir[1 ][0 ];
174
174
newDirPermute[0 ][1 ] = origDir[1 ][1 ];
175
175
newDirPermute[1 ][0 ] = origDir[0 ][0 ];
176
176
newDirPermute[1 ][1 ] = origDir[0 ][1 ];
177
177
178
- dispacementfieldPermuted ->SetDirection (newDirPermute);
179
- dispacementfieldPermuted ->SetOrigin (dispacementfield ->GetOrigin ());
180
- dispacementfieldPermuted ->SetSpacing (dispacementfield ->GetSpacing ());
178
+ displacementFieldPermuted ->SetDirection (newDirPermute);
179
+ displacementFieldPermuted ->SetOrigin (displacementField ->GetOrigin ());
180
+ displacementFieldPermuted ->SetSpacing (displacementField ->GetSpacing ());
181
181
182
- const float detPermuted = GetDeterminantAtPoint (dispacementfieldPermuted , physPt);
182
+ const float detPermuted = GetDeterminantAtPoint (displacementFieldPermuted , physPt);
183
183
184
184
using FlipFilterType = itk::FlipImageFilter<FieldType>;
185
185
auto flip = FlipFilterType::New ();
@@ -188,22 +188,22 @@ TestDisplacementJacobianDeterminantValue()
188
188
flipAxes[0 ] = true ; // Flip X
189
189
flipAxes[1 ] = false ; // Keep Y
190
190
flip->SetFlipAxes (flipAxes);
191
- flip->SetInput (dispacementfield );
191
+ flip->SetInput (displacementField );
192
192
flip->FlipAboutOriginOff ();
193
193
flip->Update ();
194
194
195
- auto dispacementfieldFlip = flip->GetOutput ();
195
+ auto displacementFieldFlip = flip->GetOutput ();
196
196
197
197
// Adjust direction to compensate flip
198
198
itk::Matrix<double , 2 , 2 > flipMat;
199
199
flipMat.SetIdentity ();
200
200
flipMat[0 ][0 ] = -1.0 ;
201
201
202
- itk::Matrix<double , 2 , 2 > newDirFlip = flipMat * dispacementfield ->GetDirection ();
202
+ itk::Matrix<double , 2 , 2 > newDirFlip = flipMat * displacementField ->GetDirection ();
203
203
204
- dispacementfieldFlip ->SetDirection (newDirFlip);
204
+ displacementFieldFlip ->SetDirection (newDirFlip);
205
205
206
- const float detFlipped = GetDeterminantAtPoint (dispacementfieldFlip , physPt);
206
+ const float detFlipped = GetDeterminantAtPoint (displacementFieldFlip , physPt);
207
207
208
208
std::cout << " Determinant at point " << physPt << " :" << std::endl;
209
209
std::cout << " Original: " << detOriginal << std::endl;
0 commit comments