Medical Imaging Interaction Toolkit  2016.11.0
Medical Imaging Interaction Toolkit
mitkAbstractTransformGeometryTest.cpp
Go to the documentation of this file.
1 /*===================================================================
2 
3 The Medical Imaging Interaction Toolkit (MITK)
4 
5 Copyright (c) German Cancer Research Center,
6 Division of Medical and Biological Informatics.
7 All rights reserved.
8 
9 This software is distributed WITHOUT ANY WARRANTY; without
10 even the implied warranty of MERCHANTABILITY or FITNESS FOR
11 A PARTICULAR PURPOSE.
12 
13 See LICENSE.txt or http://www.mitk.org for details.
14 
15 ===================================================================*/
16 
17 #include "mitkExternAbstractTransformGeometry.h"
18 #include "mitkImage.h"
19 #include "mitkPlaneGeometry.h"
20 #include "mitkSlicedGeometry3D.h"
21 
22 #include <vtkSphericalTransform.h>
23 
24 #include <fstream>
25 
26 int mitkAbstractTransformGeometryTest(int /*argc*/, char * /*argv*/ [])
27 {
28  mitk::Point3D origin;
29  mitk::Vector3D right, bottom;
30  mitk::ScalarType width, height;
31  mitk::ScalarType widthInMM, heightInMM;
32 
33  std::cout << "Initializing an x-/y-plane (xyPlane) as parameter plane by InitializeStandardPlane(rightVector, "
34  "downVector, spacing = NULL): "
35  << std::endl;
36  mitk::PlaneGeometry::Pointer xyPlane = mitk::PlaneGeometry::New();
37  width = 100;
38  widthInMM = width;
39  height = 200;
40  heightInMM = height;
41  mitk::ScalarType bounds[6] = {0, width, 0, height, 0, 1};
42  mitk::FillVector3D(origin, 4.5, 7.3, 11.2);
43  mitk::FillVector3D(right, widthInMM, 0, 0);
44  mitk::FillVector3D(bottom, 0, heightInMM, 0);
45  xyPlane->InitializeStandardPlane(right, bottom);
46  xyPlane->SetOrigin(origin);
47  xyPlane->SetSizeInUnits(width, height);
48 
49  std::cout << "Creating AbstractTransformGeometry: " << std::endl;
50  mitk::ExternAbstractTransformGeometry::Pointer abstractgeometry = mitk::ExternAbstractTransformGeometry::New();
51 
52  std::cout << "Setting xyPlane as parameter plane of AbstractTransformGeometry: " << std::endl;
53  abstractgeometry->SetPlane(xyPlane);
54 
55  std::cout
56  << "Testing whether the bounds of xyPlane and the parametric bounds of AbstractTransformGeometry are equal: ";
57  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
58  const_cast<mitk::BoundingBox *>(xyPlane->GetBoundingBox())->GetMinimum()) == false) ||
59  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
60  const_cast<mitk::BoundingBox *>(xyPlane->GetBoundingBox())->GetMaximum()) == false))
61  {
62  std::cout << "[FAILED]" << std::endl;
63  return EXIT_FAILURE;
64  }
65  std::cout << "[PASSED]" << std::endl;
66 
67  std::cout << "Testing whether the parametic bounds of AbstractTransformGeometry and the bounds of the plane accessed "
68  "from there are equal: ";
69  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
70  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
71  false) ||
72  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
73  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
74  false))
75  {
76  std::cout << "[FAILED]" << std::endl;
77  return EXIT_FAILURE;
78  }
79  std::cout << "[PASSED]" << std::endl;
80 
81  std::cout << "Change parametic bounds of AbstractTransformGeometry and test whether they are equal to the bounds of "
82  "the plane accessed from there: "
83  << std::endl;
84  height = 300;
85  bounds[3] = height;
86  abstractgeometry->SetParametricBounds(bounds);
87  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
88  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
89  false) ||
90  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
91  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
92  false))
93  {
94  std::cout << "[FAILED]" << std::endl;
95  return EXIT_FAILURE;
96  }
97  std::cout << "[PASSED]" << std::endl;
98 
99  std::cout << "Initializing an phi-/theta-plane (sphereParameterPlane) as parameter plane by "
100  "InitializeStandardPlane(rightVector, downVector, spacing = NULL): "
101  << std::endl;
102  mitk::PlaneGeometry::Pointer sphereParameterPlane = mitk::PlaneGeometry::New();
103  width = 100;
104  widthInMM = 2 * vnl_math::pi;
105  height = 200;
106  heightInMM = vnl_math::pi;
107  mitk::ScalarType radiusInMM = 2.5;
108 
109  mitk::FillVector3D(origin, radiusInMM, 0, widthInMM);
110  mitk::FillVector3D(right, 0, 0, -widthInMM);
111  mitk::FillVector3D(bottom, 0, heightInMM, 0);
112  sphereParameterPlane->InitializeStandardPlane(right, bottom);
113  sphereParameterPlane->SetOrigin(origin);
114  sphereParameterPlane->SetSizeInUnits(width, height);
115 
116  std::cout << "Creating an vtkSphericalTransform (sphericalTransform) to use with sphereParameterPlane: " << std::endl;
117  vtkSphericalTransform *sphericalTransform = vtkSphericalTransform::New();
118 
119  std::cout << "Setting sphereParameterPlane as parameter plane and sphericalTransform as transform of "
120  "AbstractTransformGeometry: "
121  << std::endl;
122  abstractgeometry->SetPlane(sphereParameterPlane);
123  abstractgeometry->SetVtkAbstractTransform(sphericalTransform);
124 
125  std::cout << "Testing whether the bounds of sphereParameterPlane and the parametric bounds of "
126  "AbstractTransformGeometry are equal: ";
127  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
128  const_cast<mitk::BoundingBox *>(sphereParameterPlane->GetBoundingBox())->GetMinimum()) == false) ||
129  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
130  const_cast<mitk::BoundingBox *>(sphereParameterPlane->GetBoundingBox())->GetMaximum()) == false))
131  {
132  std::cout << "[FAILED]" << std::endl;
133  return EXIT_FAILURE;
134  }
135  std::cout << "[PASSED]" << std::endl;
136 
137  std::cout << "Testing whether the parametic bounds of AbstractTransformGeometry and the bounds of the plane accessed "
138  "from there are equal: ";
139  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
140  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
141  false) ||
142  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
143  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
144  false))
145  {
146  std::cout << "[FAILED]" << std::endl;
147  return EXIT_FAILURE;
148  }
149  std::cout << "[PASSED]" << std::endl;
150 
151  std::cout << "Testing mapping Map(pt2d_mm(phi=Pi,theta=Pi/2.0), pt3d_mm) and compare with expected (-radius, 0, 0): ";
152  mitk::Point2D pt2d_mm;
153  mitk::Point3D pt3d_mm, expected_pt3d_mm;
154  pt2d_mm[0] = vnl_math::pi;
155  pt2d_mm[1] = vnl_math::pi_over_2;
156  mitk::FillVector3D(expected_pt3d_mm, -radiusInMM, 0, 0);
157  abstractgeometry->Map(pt2d_mm, pt3d_mm);
158  if (mitk::Equal(pt3d_mm, expected_pt3d_mm) == false)
159  {
160  std::cout << "[FAILED]" << std::endl;
161  return EXIT_FAILURE;
162  }
163  std::cout << "[PASSED]" << std::endl;
164 
165  std::cout << "Testing mapping Map(pt3d_mm, pt2d_mm) and compare with expected: ";
166  mitk::Point2D testpt2d_mm;
167  abstractgeometry->Map(pt3d_mm, testpt2d_mm);
168  if (mitk::Equal(pt2d_mm, testpt2d_mm) == false)
169  {
170  std::cout << "[FAILED]" << std::endl;
171  return EXIT_FAILURE;
172  }
173  std::cout << "[PASSED]" << std::endl;
174 
175  std::cout << "Testing IndexToWorld(pt2d_units, pt2d_mm) and compare with expected: ";
176  mitk::Point2D pt2d_units;
177  pt2d_units[0] = width / 2.0;
178  pt2d_units[1] = height / 2.0;
179  pt2d_mm[0] = widthInMM / 2.0;
180  pt2d_mm[1] = heightInMM / 2.0;
181  abstractgeometry->IndexToWorld(pt2d_units, testpt2d_mm);
182  if (mitk::Equal(pt2d_mm, testpt2d_mm) == false)
183  {
184  std::cout << "[FAILED]" << std::endl;
185  return EXIT_FAILURE;
186  }
187  std::cout << "[PASSED]" << std::endl;
188 
189  std::cout << "Change parametic bounds of AbstractTransformGeometry and test whether they are equal to the bounds of "
190  "the plane accessed from there: "
191  << std::endl;
192  height = 300;
193  bounds[3] = height;
194  abstractgeometry->SetParametricBounds(bounds);
195  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
196  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
197  false) ||
198  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
199  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
200  false))
201  {
202  std::cout << "[FAILED]" << std::endl;
203  return EXIT_FAILURE;
204  }
205  std::cout << "[PASSED]" << std::endl;
206 
207  std::cout << "Testing IndexToWorld(pt2d_units, pt2d_mm) and compare with expected: ";
208  pt2d_units[0] = width / 2.0;
209  pt2d_units[1] = height / 2.0;
210  pt2d_mm[0] = widthInMM / 2.0;
211  pt2d_mm[1] = heightInMM / 2.0;
212  abstractgeometry->IndexToWorld(pt2d_units, testpt2d_mm);
213  if (mitk::Equal(pt2d_mm, testpt2d_mm) == false)
214  {
215  std::cout << "[FAILED]" << std::endl;
216  return EXIT_FAILURE;
217  }
218  std::cout << "[PASSED]" << std::endl;
219 
220  // std::cout << "Testing availability and type (PlaneGeometry) of first geometry in the SlicedGeometry3D: ";
221  // mitk::PlaneGeometry* accessedplanegeometry3 =
222  // dynamic_cast<mitk::PlaneGeometry*>(slicedWorldGeometry->GetGeometry2D(0));
223  // if(accessedplanegeometry3==NULL)
224  //{
225  // std::cout<<"[FAILED]"<<std::endl;
226  // return EXIT_FAILURE;
227  //}
228  // std::cout<<"[PASSED]"<<std::endl;
229 
230  // std::cout << "Testing whether the first geometry in the SlicedGeometry3D is identical to planegeometry3 by axis
231  // comparison: "<<std::endl;
232  // if((mitk::Equal(accessedplanegeometry3->GetAxisVector(0), planegeometry3->GetAxisVector(0))==false) ||
233  // (mitk::Equal(accessedplanegeometry3->GetAxisVector(1), planegeometry3->GetAxisVector(1))==false) ||
234  // (mitk::Equal(accessedplanegeometry3->GetAxisVector(2), planegeometry3->GetAxisVector(2))==false))
235  //{
236  // std::cout<<"[FAILED]"<<std::endl;
237  // return EXIT_FAILURE;
238  //}
239  // std::cout<<"[PASSED]"<<std::endl;
240 
241  sphericalTransform->Delete();
242 
243  std::cout << "[TEST DONE]" << std::endl;
244  return EXIT_SUCCESS;
245 }
mitk::BoundingBox
itk::BoundingBox< unsigned long, 3, ScalarType > BoundingBox
Standard 3D-BoundingBox typedef.
Definition: mitkBaseGeometry.h:45
mitkAbstractTransformGeometryTest
int mitkAbstractTransformGeometryTest(int, char *[])
Definition: mitkAbstractTransformGeometryTest.cpp:26
mitk::Pointer
itk::SmartPointer< Self > Pointer
Definition: mitkRenderingManager.h:389
mitk::ScalarType
double ScalarType
Definition: mitkNumericConstants.h:24
mitk::FillVector3D
void FillVector3D(Tout &out, mitk::ScalarType x, mitk::ScalarType y, mitk::ScalarType z)
Definition: mitkArray.h:110
mitk::Equal
MITKNEWMODULE_EXPORT bool Equal(mitk::ExampleDataStructure *leftHandSide, mitk::ExampleDataStructure *rightHandSide, mitk::ScalarType eps, bool verbose)
Returns true if the example data structures are considered equal.
Definition: mitkExampleDataStructure.cpp:60
mitk::PlaneGeometry::New
static Pointer New()
mitk::New
static itkEventMacro(BoundingShapeInteractionEvent, itk::AnyEvent) class MITKBOUNDINGSHAPE_EXPORT BoundingShapeInteractor Pointer New()
Basic interaction methods for mitk::GeometryData.