Medical Imaging Interaction Toolkit  2016.11.0
Medical Imaging Interaction Toolkit
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
mitkAbstractTransformPlaneGeometryTest.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 "mitkImage.h"
18 #include "mitkPlaneGeometry.h"
19 #include "mitkSlicedGeometry3D.h"
20 #include "mitkVtkAbstractTransformPlaneGeometry.h"
21 
22 #include <vtkSphericalTransform.h>
23 
24 #include <vnl/vnl_quaternion.h>
25 #include <vnl/vnl_quaternion.txx>
26 
27 #include <fstream>
28 
29 int mitkVtkAbstractTransformPlaneGeometryTest(int argc, char *argv[])
30 {
31  mitk::Point3D origin;
32  mitk::Vector3D right, bottom;
33  mitk::ScalarType width, height;
34  mitk::ScalarType widthInMM, heightInMM;
35 
36  std::cout << "Initializing an x-/y-plane (xyPlane) as parameter plane by InitializeStandardPlane(rightVector, "
37  "downVector, spacing = NULL): "
38  << std::endl;
39  mitk::PlaneGeometry::Pointer xyPlane = mitk::PlaneGeometry::New();
40  width = 100;
41  widthInMM = width;
42  height = 200;
43  heightInMM = height;
44  mitk::ScalarType bounds[6] = {0, width, 0, height, 0, 1};
45  mitk::FillVector3D(origin, 4.5, 7.3, 11.2);
46  mitk::FillVector3D(right, widthInMM, 0, 0);
47  mitk::FillVector3D(bottom, 0, heightInMM, 0);
48  xyPlane->InitializeStandardPlane(right, bottom);
49  xyPlane->SetOrigin(origin);
50  xyPlane->SetSizeInUnits(width, height);
51 
52  std::cout << "Creating VtkAbstractTransformPlaneGeometry: " << std::endl;
53  mitk::VtkAbstractTransformPlaneGeometry::Pointer abstractgeometry = mitk::VtkAbstractTransformPlaneGeometry::New();
54 
55  std::cout << "Setting xyPlane as parameter plane of VtkAbstractTransformPlaneGeometry: " << std::endl;
56  abstractgeometry->SetPlane(xyPlane);
57 
58  std::cout << "Testing whether the bounds of xyPlane and the parametric bounds of VtkAbstractTransformPlaneGeometry "
59  "are equal: ";
60  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
61  const_cast<mitk::BoundingBox *>(xyPlane->GetBoundingBox())->GetMinimum()) == false) ||
62  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
63  const_cast<mitk::BoundingBox *>(xyPlane->GetBoundingBox())->GetMaximum()) == false))
64  {
65  std::cout << "[FAILED]" << std::endl;
66  return EXIT_FAILURE;
67  }
68  std::cout << "[PASSED]" << std::endl;
69 
70  std::cout << "Testing whether the parametic bounds of VtkAbstractTransformPlaneGeometry and the bounds of the plane "
71  "accessed from there are equal: ";
72  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
73  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
74  false) ||
75  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
76  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
77  false))
78  {
79  std::cout << "[FAILED]" << std::endl;
80  return EXIT_FAILURE;
81  }
82  std::cout << "[PASSED]" << std::endl;
83 
84  std::cout << "Change parametic bounds of VtkAbstractTransformPlaneGeometry and test whether they are equal to the "
85  "bounds of the plane accessed from there: "
86  << std::endl;
87  height = 300;
88  bounds[3] = height;
89  abstractgeometry->SetParametricBounds(bounds);
90  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
91  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
92  false) ||
93  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
94  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
95  false))
96  {
97  std::cout << "[FAILED]" << std::endl;
98  return EXIT_FAILURE;
99  }
100  std::cout << "[PASSED]" << std::endl;
101 
102  std::cout << "Initializing an phi-/theta-plane (sphereParameterPlane) as parameter plane by "
103  "InitializeStandardPlane(rightVector, downVector, spacing = NULL): "
104  << std::endl;
105  mitk::PlaneGeometry::Pointer sphereParameterPlane = mitk::PlaneGeometry::New();
106  width = 100;
107  widthInMM = 2 * vnl_math::pi;
108  height = 200;
109  heightInMM = vnl_math::pi;
110  mitk::ScalarType radiusInMM = 2.5;
111 
112  mitk::FillVector3D(origin, radiusInMM, 0, widthInMM);
113  mitk::FillVector3D(right, 0, 0, -widthInMM);
114  mitk::FillVector3D(bottom, 0, heightInMM, 0);
115  sphereParameterPlane->InitializeStandardPlane(right, bottom);
116  sphereParameterPlane->SetOrigin(origin);
117  sphereParameterPlane->SetSizeInUnits(width, height);
118 
119  std::cout << "Creating an vtkSphericalTransform (sphericalTransform) to use with sphereParameterPlane: " << std::endl;
120  vtkSphericalTransform *sphericalTransform = vtkSphericalTransform::New();
121 
122  std::cout << "Setting sphereParameterPlane as parameter plane and sphericalTransform as transform of "
123  "VtkAbstractTransformPlaneGeometry: "
124  << std::endl;
125  abstractgeometry->SetPlane(sphereParameterPlane);
126  abstractgeometry->SetVtkAbstractTransform(sphericalTransform);
127 
128  std::cout << "Testing whether the bounds of sphereParameterPlane and the parametric bounds of "
129  "VtkAbstractTransformPlaneGeometry are equal: ";
130  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
131  const_cast<mitk::BoundingBox *>(sphereParameterPlane->GetBoundingBox())->GetMinimum()) == false) ||
132  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
133  const_cast<mitk::BoundingBox *>(sphereParameterPlane->GetBoundingBox())->GetMaximum()) == false))
134  {
135  std::cout << "[FAILED]" << std::endl;
136  return EXIT_FAILURE;
137  }
138  std::cout << "[PASSED]" << std::endl;
139 
140  std::cout << "Testing whether the parametic bounds of VtkAbstractTransformPlaneGeometry and the bounds of the plane "
141  "accessed from there are equal: ";
142  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
143  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
144  false) ||
145  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
146  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
147  false))
148  {
149  std::cout << "[FAILED]" << std::endl;
150  return EXIT_FAILURE;
151  }
152  std::cout << "[PASSED]" << std::endl;
153 
154  std::cout << "Testing mapping Map(pt2d_mm(phi=Pi,theta=Pi/2.0), pt3d_mm) and compare with expected (-radius, 0, 0): ";
155  mitk::Point2D pt2d_mm;
156  mitk::Point3D pt3d_mm, expected_pt3d_mm;
157  pt2d_mm[0] = vnl_math::pi;
158  pt2d_mm[1] = vnl_math::pi_over_2;
159  mitk::FillVector3D(expected_pt3d_mm, -radiusInMM, 0, 0);
160  abstractgeometry->Map(pt2d_mm, pt3d_mm);
161  if (mitk::Equal(pt3d_mm, expected_pt3d_mm) == false)
162  {
163  std::cout << "[FAILED]" << std::endl;
164  return EXIT_FAILURE;
165  }
166  std::cout << "[PASSED]" << std::endl;
167 
168  std::cout << "Testing mapping Map(pt3d_mm, pt2d_mm) and compare with expected: ";
169  mitk::Point2D testpt2d_mm;
170  abstractgeometry->Map(pt3d_mm, testpt2d_mm);
171  if (mitk::Equal(pt2d_mm, testpt2d_mm) == false)
172  {
173  std::cout << "[FAILED]" << std::endl;
174  return EXIT_FAILURE;
175  }
176  std::cout << "[PASSED]" << std::endl;
177 
178  std::cout << "Testing IndexToWorld(pt2d_units, pt2d_mm) and compare with expected: ";
179  mitk::Point2D pt2d_units;
180  pt2d_units[0] = width / 2.0;
181  pt2d_units[1] = height / 2.0;
182  pt2d_mm[0] = widthInMM / 2.0;
183  pt2d_mm[1] = heightInMM / 2.0;
184  abstractgeometry->IndexToWorld(pt2d_units, testpt2d_mm);
185  if (mitk::Equal(pt2d_mm, testpt2d_mm) == false)
186  {
187  std::cout << "[FAILED]" << std::endl;
188  return EXIT_FAILURE;
189  }
190  std::cout << "[PASSED]" << std::endl;
191 
192  std::cout << "Change parametic bounds of VtkAbstractTransformPlaneGeometry and test whether they are equal to the "
193  "bounds of the plane accessed from there: "
194  << std::endl;
195  height = 300;
196  bounds[3] = height;
197  abstractgeometry->SetParametricBounds(bounds);
198  if ((mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMinimum(),
199  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMinimum()) ==
200  false) ||
201  (mitk::Equal(const_cast<mitk::BoundingBox *>(abstractgeometry->GetParametricBoundingBox())->GetMaximum(),
202  const_cast<mitk::BoundingBox *>(abstractgeometry->GetPlane()->GetBoundingBox())->GetMaximum()) ==
203  false))
204  {
205  std::cout << "[FAILED]" << std::endl;
206  return EXIT_FAILURE;
207  }
208  std::cout << "[PASSED]" << std::endl;
209 
210  std::cout << "Testing IndexToWorld(pt2d_units, pt2d_mm) and compare with expected: ";
211  pt2d_units[0] = width / 2.0;
212  pt2d_units[1] = height / 2.0;
213  pt2d_mm[0] = widthInMM / 2.0;
214  pt2d_mm[1] = heightInMM / 2.0;
215  abstractgeometry->IndexToWorld(pt2d_units, testpt2d_mm);
216  if (mitk::Equal(pt2d_mm, testpt2d_mm) == false)
217  {
218  std::cout << "[FAILED]" << std::endl;
219  return EXIT_FAILURE;
220  }
221  std::cout << "[PASSED]" << std::endl;
222 
223  // std::cout << "Testing availability and type (PlaneGeometry) of first geometry in the SlicedGeometry3D: ";
224  // mitk::PlaneGeometry* accessedplanegeometry3 =
225  // dynamic_cast<mitk::PlaneGeometry*>(slicedWorldGeometry->GetGeometry2D(0));
226  // if(accessedplanegeometry3==NULL)
227  //{
228  // std::cout<<"[FAILED]"<<std::endl;
229  // return EXIT_FAILURE;
230  //}
231  // std::cout<<"[PASSED]"<<std::endl;
232 
233  // std::cout << "Testing whether the first geometry in the SlicedGeometry3D is identical to planegeometry3 by axis
234  // comparison: "<<std::endl;
235  // if((mitk::Equal(accessedplanegeometry3->GetAxisVector(0), planegeometry3->GetAxisVector(0))==false) ||
236  // (mitk::Equal(accessedplanegeometry3->GetAxisVector(1), planegeometry3->GetAxisVector(1))==false) ||
237  // (mitk::Equal(accessedplanegeometry3->GetAxisVector(2), planegeometry3->GetAxisVector(2))==false))
238  //{
239  // std::cout<<"[FAILED]"<<std::endl;
240  // return EXIT_FAILURE;
241  //}
242  // std::cout<<"[PASSED]"<<std::endl;
243 
244  std::cout << "[TEST DONE]" << std::endl;
245  return EXIT_SUCCESS;
246 }
mitk::BoundingBox
itk::BoundingBox< unsigned long, 3, ScalarType > BoundingBox
Standard 3D-BoundingBox typedef.
Definition: mitkBaseGeometry.h:45
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
mitkVtkAbstractTransformPlaneGeometryTest
int mitkVtkAbstractTransformPlaneGeometryTest(int argc, char *argv[])
Definition: mitkAbstractTransformPlaneGeometryTest.cpp:29
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.