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
QmitkFixedCenterOfRotationAffineTransformView.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 "QmitkFixedCenterOfRotationAffineTransformView.h"
18 #include "mitkImageAccessByItk.h"
19 #include <QValidator>
20 #include <itkCenteredTransformInitializer.h>
21 #include <itkFixedCenterOfRotationAffineTransform.h>
22 #include <mitkImageCast.h>
23 
24 QmitkFixedCenterOfRotationAffineTransformView::QmitkFixedCenterOfRotationAffineTransformView(QWidget *parent,
25  Qt::WindowFlags f)
26  : QmitkRigidRegistrationTransformsGUIBase(parent, f), m_CenterX(0), m_CenterY(0), m_CenterZ(0)
27 {
28 }
29 
30 QmitkFixedCenterOfRotationAffineTransformView::~QmitkFixedCenterOfRotationAffineTransformView()
31 {
32 }
33 
34 mitk::TransformParameters::TransformType QmitkFixedCenterOfRotationAffineTransformView::GetTransformType()
35 {
36  return mitk::TransformParameters::FIXEDCENTEROFROTATIONAFFINETRANSFORM;
37 }
38 
39 itk::Object::Pointer QmitkFixedCenterOfRotationAffineTransformView::GetTransform()
40 {
41  if (m_FixedImage.IsNotNull())
42  {
43  AccessByItk(m_FixedImage, GetTransform2);
44  return m_TransformObject;
45  }
46  return nullptr;
47 }
48 
49 template <class TPixelType, unsigned int VImageDimension>
50 itk::Object::Pointer QmitkFixedCenterOfRotationAffineTransformView::GetTransform2(
51  itk::Image<TPixelType, VImageDimension> *itkImage1)
52 {
53  typedef typename itk::Image<TPixelType, VImageDimension> FixedImageType;
54  typedef typename itk::Image<TPixelType, VImageDimension> MovingImageType;
55 
56  // the fixedImage is the input parameter (fix for Bug #14626)
57  typename FixedImageType::Pointer fixedImage = itkImage1;
58 
59  // the movingImage type is known, use the ImageToItk filter (fix for Bug #14626)
60  typename mitk::ImageToItk<MovingImageType>::Pointer movingImageToItk = mitk::ImageToItk<MovingImageType>::New();
61  movingImageToItk->SetInput(m_MovingImage);
62  movingImageToItk->Update();
63  typename MovingImageType::Pointer movingImage = movingImageToItk->GetOutput();
64 
65  typedef typename itk::FixedCenterOfRotationAffineTransform<double, VImageDimension> CenteredAffineTransformType;
66  typename itk::FixedCenterOfRotationAffineTransform<double, VImageDimension>::Pointer transformPointer =
67  itk::FixedCenterOfRotationAffineTransform<double, VImageDimension>::New();
68  transformPointer->SetIdentity();
69  if (m_Controls.m_CenterForInitializerFixedCenterOfRotationAffine->isChecked())
70  {
71  typedef typename itk::FixedCenterOfRotationAffineTransform<double, VImageDimension>
72  FixedCenterOfRotationAffineTransformType;
73  typedef typename itk::
74  CenteredTransformInitializer<FixedCenterOfRotationAffineTransformType, FixedImageType, MovingImageType>
75  TransformInitializerType;
76  typename TransformInitializerType::Pointer transformInitializer = TransformInitializerType::New();
77  transformInitializer->SetFixedImage(fixedImage);
78  transformInitializer->SetMovingImage(movingImage);
79  transformInitializer->SetTransform(transformPointer);
80  if (m_Controls.m_MomentsFixedCenterOfRotationAffine->isChecked())
81  {
82  transformInitializer->MomentsOn();
83  }
84  else
85  {
86  transformInitializer->GeometryOn();
87  }
88  m_CenterX = transformPointer->GetCenter()[0];
89  m_CenterY = transformPointer->GetCenter()[1];
90  m_CenterZ = transformPointer->GetCenter()[2];
91  transformInitializer->InitializeTransform();
92  }
93  m_TransformObject = transformPointer.GetPointer();
94  return transformPointer.GetPointer();
95 }
96 
97 itk::Array<double> QmitkFixedCenterOfRotationAffineTransformView::GetTransformParameters()
98 {
99  itk::Array<double> transformValues;
100  transformValues.SetSize(15);
101  transformValues.fill(0);
102  transformValues[0] = m_Controls.m_UseOptimizerScalesFixedCenterOfRotationAffine->isChecked();
103  transformValues[1] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale1->text().toDouble();
104  transformValues[2] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale2->text().toDouble();
105  transformValues[3] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale3->text().toDouble();
106  transformValues[4] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale4->text().toDouble();
107  transformValues[5] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale5->text().toDouble();
108  transformValues[6] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale6->text().toDouble();
109  transformValues[7] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale7->text().toDouble();
110  transformValues[8] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale8->text().toDouble();
111  transformValues[9] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale9->text().toDouble();
112  transformValues[10] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationX->text().toDouble();
113  transformValues[11] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationY->text().toDouble();
114  transformValues[12] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationZ->text().toDouble();
115  transformValues[13] = m_Controls.m_CenterForInitializerFixedCenterOfRotationAffine->isChecked();
116  transformValues[14] = m_Controls.m_MomentsFixedCenterOfRotationAffine->isChecked();
117  return transformValues;
118 }
119 
120 void QmitkFixedCenterOfRotationAffineTransformView::SetTransformParameters(itk::Array<double> transformValues)
121 {
122  m_Controls.m_UseOptimizerScalesFixedCenterOfRotationAffine->setChecked(transformValues[0]);
123  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale1->setText(QString::number(transformValues[1]));
124  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale2->setText(QString::number(transformValues[2]));
125  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale3->setText(QString::number(transformValues[3]));
126  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale4->setText(QString::number(transformValues[4]));
127  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale5->setText(QString::number(transformValues[5]));
128  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale6->setText(QString::number(transformValues[6]));
129  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale7->setText(QString::number(transformValues[7]));
130  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale8->setText(QString::number(transformValues[8]));
131  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale9->setText(QString::number(transformValues[9]));
132  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationX->setText(
133  QString::number(transformValues[10]));
134  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationY->setText(
135  QString::number(transformValues[11]));
136  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationZ->setText(
137  QString::number(transformValues[12]));
138  m_Controls.m_CenterForInitializerFixedCenterOfRotationAffine->setChecked(transformValues[13]);
139  m_Controls.m_MomentsFixedCenterOfRotationAffine->setChecked(transformValues[14]);
140  m_Controls.m_GeometryFixedCenterOfRotationAffine->setChecked(!transformValues[14]);
141 }
142 
143 QString QmitkFixedCenterOfRotationAffineTransformView::GetName()
144 {
145  return "FixedCenterOfRotationAffine";
146 }
147 
148 void QmitkFixedCenterOfRotationAffineTransformView::SetupUI(QWidget *parent)
149 {
150  m_Controls.setupUi(parent);
151  QValidator *validatorLineEditInputFloat = new QDoubleValidator(0, 20000000, 8, this);
152  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale1->setValidator(validatorLineEditInputFloat);
153  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale2->setValidator(validatorLineEditInputFloat);
154  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale3->setValidator(validatorLineEditInputFloat);
155  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale4->setValidator(validatorLineEditInputFloat);
156  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale5->setValidator(validatorLineEditInputFloat);
157  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale6->setValidator(validatorLineEditInputFloat);
158  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale7->setValidator(validatorLineEditInputFloat);
159  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale8->setValidator(validatorLineEditInputFloat);
160  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale9->setValidator(validatorLineEditInputFloat);
161  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationX->setValidator(validatorLineEditInputFloat);
162  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationY->setValidator(validatorLineEditInputFloat);
163  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationZ->setValidator(validatorLineEditInputFloat);
164 }
165 
166 itk::Array<double> QmitkFixedCenterOfRotationAffineTransformView::GetScales()
167 {
168  itk::Array<double> scales;
169  scales.SetSize(12);
170  scales.Fill(1.0);
171  if (m_Controls.m_UseOptimizerScalesFixedCenterOfRotationAffine->isChecked())
172  {
173  scales[0] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale1->text().toDouble();
174  scales[1] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale2->text().toDouble();
175  scales[2] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale3->text().toDouble();
176  scales[3] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale4->text().toDouble();
177  scales[4] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale5->text().toDouble();
178  scales[5] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale6->text().toDouble();
179  scales[6] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale7->text().toDouble();
180  scales[7] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale8->text().toDouble();
181  scales[8] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale9->text().toDouble();
182  scales[9] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationX->text().toDouble();
183  scales[10] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationY->text().toDouble();
184  scales[11] = m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationZ->text().toDouble();
185  }
186  return scales;
187 }
188 
189 vtkTransform *QmitkFixedCenterOfRotationAffineTransformView::Transform(vtkMatrix4x4 *vtkmatrix,
190  vtkTransform *vtktransform,
191  itk::Array<double> transformParams)
192 {
193  if (m_MovingImage.IsNotNull())
194  {
195  int m = 0;
196  for (unsigned int i = 0; i < m_FixedImage->GetDimension(); i++)
197  {
198  for (unsigned int j = 0; j < m_FixedImage->GetDimension(); j++)
199  {
200  vtkmatrix->SetElement(i, j, transformParams[m]);
201  m++;
202  }
203  }
204  float center[4];
205  float translation[4];
206  center[0] = m_CenterX;
207  center[1] = m_CenterY;
208  center[2] = m_CenterZ;
209  center[3] = 1;
210  vtkmatrix->MultiplyPoint(center, translation);
211  if (m_FixedImage->GetDimension() == 2)
212  {
213  vtkmatrix->SetElement(0, 3, -translation[0] + center[0] + transformParams[4]);
214  vtkmatrix->SetElement(1, 3, -translation[1] + center[1] + transformParams[5]);
215  }
216  else if (m_FixedImage->GetDimension() == 3)
217  {
218  vtkmatrix->SetElement(0, 3, -translation[0] + center[0] + transformParams[9]);
219  vtkmatrix->SetElement(1, 3, -translation[1] + center[1] + transformParams[10]);
220  vtkmatrix->SetElement(2, 3, -translation[2] + center[2] + transformParams[11]);
221  }
222  vtktransform->SetMatrix(vtkmatrix);
223  }
224  return vtktransform;
225 }
226 
227 int QmitkFixedCenterOfRotationAffineTransformView::GetNumberOfTransformParameters()
228 {
229  if (m_FixedImage.IsNotNull())
230  {
231  if (m_FixedImage->GetDimension() == 2)
232  {
233  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale7->hide();
234  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale8->hide();
235  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale9->hide();
236  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationX->hide();
237  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationY->hide();
238  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationZ->hide();
239  m_Controls.textLabel2_7_2_2->setText("Translation Scale X:");
240  m_Controls.textLabel3_6_2_2->setText("Translation Scale Y:");
241  m_Controls.textLabel4_4_3_2->hide();
242  m_Controls.textLabel5_4_2_2->hide();
243  m_Controls.textLabel6_4_2_2->hide();
244  m_Controls.textLabel11_3_2_2->hide();
245  m_Controls.textLabel12_3_2_2->hide();
246  m_Controls.textLabel13_2_2_2->hide();
247  return 6;
248  }
249  else
250  {
251  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale7->show();
252  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale8->show();
253  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScale9->show();
254  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationX->show();
255  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationY->show();
256  m_Controls.m_ScalesFixedCenterOfRotationAffineTransformScaleTranslationZ->show();
257  m_Controls.textLabel2_7_2_2->setText("Scale 5:");
258  m_Controls.textLabel3_6_2_2->setText("Scale 6:");
259  m_Controls.textLabel4_4_3_2->show();
260  m_Controls.textLabel5_4_2_2->show();
261  m_Controls.textLabel6_4_2_2->show();
262  m_Controls.textLabel11_3_2_2->show();
263  m_Controls.textLabel12_3_2_2->show();
264  m_Controls.textLabel13_2_2_2->show();
265  return 12;
266  }
267  }
268  else
269  return 0;
270 }
mitk::Pointer
itk::SmartPointer< Self > Pointer
Definition: mitkRenderingManager.h:389
QmitkFixedCenterOfRotationAffineTransformView::GetTransformType
virtual mitk::TransformParameters::TransformType GetTransformType() override
Definition: QmitkFixedCenterOfRotationAffineTransformView.cpp:34
QmitkFixedCenterOfRotationAffineTransformView::Transform
virtual vtkTransform * Transform(vtkMatrix4x4 *vtkmatrix, vtkTransform *vtktransform, itk::Array< double > transformParams) override
Definition: QmitkFixedCenterOfRotationAffineTransformView.cpp:189
QmitkFixedCenterOfRotationAffineTransformView::m_Controls
Ui::QmitkFixedCenterOfRotationAffineTransformControls m_Controls
Definition: QmitkFixedCenterOfRotationAffineTransformView.h:62
QmitkFixedCenterOfRotationAffineTransformView::SetTransformParameters
virtual void SetTransformParameters(itk::Array< double > transformValues) override
Definition: QmitkFixedCenterOfRotationAffineTransformView.cpp:120
QmitkFixedCenterOfRotationAffineTransformView::QmitkFixedCenterOfRotationAffineTransformView
QmitkFixedCenterOfRotationAffineTransformView(QWidget *parent=nullptr, Qt::WindowFlags f=nullptr)
Definition: QmitkFixedCenterOfRotationAffineTransformView.cpp:24
mitk::TransformParameters::TransformType
TransformType
Unique integer value for every transform.
Definition: mitkTransformParameters.h:53
mitk::TransformParameters::FIXEDCENTEROFROTATIONAFFINETRANSFORM
Definition: mitkTransformParameters.h:58
AccessByItk
#define AccessByItk(mitkImage, itkImageTypeFunction)
Access a MITK image by an ITK image.
Definition: mitkImageAccessByItk.h:194
mitk::ImageToItk::New
static Pointer New()
mitk::New
static itkEventMacro(BoundingShapeInteractionEvent, itk::AnyEvent) class MITKBOUNDINGSHAPE_EXPORT BoundingShapeInteractor Pointer New()
Basic interaction methods for mitk::GeometryData.