Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018 #include "QmitkVolumetryView.h"
00019
00020 #include "mitkNodePredicateDataType.h"
00021
00022 #include "QmitkDataStorageComboBox.h"
00023 #include "QmitkStdMultiWidget.h"
00024
00025 #include <QMessageBox>
00026
00027
00028
00029
00030 #include "itkImage.h"
00031 #include <limits>
00032 #include "mitkDataNode.h"
00033 #include "mitkRenderingManager.h"
00034 #include "mitkVolumeCalculator.h"
00035
00036 #include <qdir.h>
00037 #include <qfiledialog.h>
00038 #include <QTextStream>
00039
00040 const std::string QmitkVolumetryView::VIEW_ID = "org.mitk.views.volumetry";
00041
00042 QmitkVolumetryView::QmitkVolumetryView()
00043 : QmitkFunctionality(),
00044 m_Controls(NULL),
00045 m_MultiWidget(NULL)
00046 {
00047
00048 }
00049
00050 QmitkVolumetryView::~QmitkVolumetryView()
00051 {
00052 }
00053
00054 void QmitkVolumetryView::CreateQtPartControl(QWidget *parent)
00055 {
00056 if (!m_Controls)
00057 {
00058
00059 m_Controls = new Ui::QmitkVolumetryViewControls;
00060 m_Controls->setupUi(parent);
00061 this->CreateConnections();
00062
00063
00064 m_Controls->m_ImageSelector->SetDataStorage( this->GetDefaultDataStorage() );
00065 m_Controls->m_ImageSelector->SetPredicate( mitk::NodePredicateDataType::New("Image") );
00066 }
00067 }
00068
00069 void QmitkVolumetryView::StdMultiWidgetAvailable (QmitkStdMultiWidget &stdMultiWidget)
00070 {
00071 m_MultiWidget = &stdMultiWidget;
00072 }
00073
00074 void QmitkVolumetryView::StdMultiWidgetNotAvailable()
00075 {
00076 m_MultiWidget = NULL;
00077 }
00078
00079 void QmitkVolumetryView::CreateConnections()
00080 {
00081 if ( m_Controls )
00082 {
00083
00084 connect( m_Controls->m_ImageSelector, SIGNAL(OnSelectionChanged(const mitk::DataNode*)), this, SLOT(OnImageSelected(const mitk::DataNode*)) );
00085
00086
00087 connect( (QObject*)(m_Controls->m_ThresholdSlider), SIGNAL(valueChanged(int)), this, SLOT(OnThresholdSliderChanged(int)) );
00088
00089
00090 connect( (QObject*)(m_Controls->m_CalcButton), SIGNAL(clicked()), this, SLOT(OnCalculateVolume()) );
00091 connect( (QObject*)(m_Controls->m_TimeSeriesButton), SIGNAL(clicked()), this, SLOT(OnTimeSeriesButtonClicked()) );
00092 connect( (QObject*)(m_Controls->m_SaveCsvButton), SIGNAL(clicked()), this, SLOT(OnSaveCsvButtonClicked()) );
00093 }
00094 }
00095
00096 void QmitkVolumetryView::Activated()
00097 {
00098 QmitkFunctionality::Activated();
00099 }
00100
00101 void QmitkVolumetryView::Deactivated()
00102 {
00103 QmitkFunctionality::Deactivated();
00104 }
00105
00106
00107 void QmitkVolumetryView::OnImageSelected(const mitk::DataNode* item)
00108 {
00109
00110 if( item == NULL || item->GetData() == NULL )
00111 return;
00112
00113 if( item->GetData()->GetSource() == NULL )
00114 {
00115 m_SelectedDataNode = 0;
00116 return;
00117 }
00118
00119 m_SelectedDataNode = const_cast<mitk::DataNode*>(item);
00120
00121 std::cout << "Selected mitk::Image at address " << m_SelectedDataNode->GetData() << std::endl;
00122
00123
00124 if(item == m_OverlayNode.GetPointer())
00125 return;
00126
00127 if (m_OverlayNode)
00128 {
00129 this->GetDefaultDataStorage()->Remove(m_OverlayNode);
00130
00131
00132
00133 m_OverlayNode = NULL;
00134 }
00135
00136 this->CreateOverlayChild();
00137
00138
00139
00140
00141
00142
00143
00144 m_Controls->m_CalcButton->setEnabled(false);
00145 m_Controls->m_TimeSeriesButton->setEnabled(false);
00146 m_Controls->m_SaveCsvButton->setEnabled(false);
00147 m_Controls->m_TextEdit->clear();
00148
00149 if (m_SelectedDataNode.IsNotNull() )
00150 {
00151 mitk::Image* image = dynamic_cast<mitk::Image*>(m_SelectedDataNode->GetData());
00152 image->Update();
00153 if (image && image->IsInitialized())
00154 {
00155 if (image->GetDimension() == 4)
00156 {
00157 m_Controls->m_TimeSeriesButton->setEnabled(true);
00158 }
00159 else
00160 {
00161 m_Controls->m_CalcButton->setEnabled(true);
00162 }
00163 int minVal = (int)image->GetScalarValue2ndMin();
00164 int maxVal = (int)image->GetScalarValueMaxNoRecompute();
00165 if (minVal == maxVal)
00166 --minVal;
00167 m_Controls->m_ThresholdSlider->setMinimum(minVal);
00168 m_Controls->m_ThresholdSlider->setMaximum(maxVal);
00169 m_Controls->m_ThresholdSlider->setEnabled(true);
00170 this->UpdateSlider();
00171 mitk::RenderingManager::GetInstance()->RequestUpdateAll();
00172 }
00173 }
00174 }
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00238
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248
00249
00250
00251
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264 void QmitkVolumetryView::OnCalculateVolume()
00265 {
00266 if (m_SelectedDataNode.IsNotNull() )
00267 {
00268 mitk::Image* image = dynamic_cast<mitk::Image*>(m_SelectedDataNode->GetData());
00269 if (image)
00270 {
00271 std::cout << "Dimension:" << image->GetDimension() << std::endl;
00272 std::cout << "Dimension[3]:" << image->GetDimension(3) << std::endl;
00273 mitk::VolumeCalculator::Pointer volCalc = mitk::VolumeCalculator::New();
00274 volCalc->SetImage(image);
00275 volCalc->SetThreshold(m_Controls->m_ThresholdSlider->value());
00276 volCalc->ComputeVolume();
00277 std::stringstream vs;
00278 vs << volCalc->GetVolume() << " ml";
00279 m_Controls->m_Result->setText(vs.str().c_str() );
00280 }
00281 }
00282 }
00283
00284 void QmitkVolumetryView::OnTimeSeriesButtonClicked()
00285 {
00286 if (m_SelectedDataNode.IsNotNull() )
00287 {
00288 mitk::Image* image = dynamic_cast<mitk::Image*>(m_SelectedDataNode->GetData());
00289 if (image)
00290 {
00291 mitk::VolumeCalculator::Pointer volCalc = mitk::VolumeCalculator::New();
00292 volCalc->SetImage(image);
00293 volCalc->SetThreshold(m_Controls->m_ThresholdSlider->value());
00294 volCalc->ComputeVolume();
00295 std::vector<float> volumes = volCalc->GetVolumes();
00296 std::stringstream vs;
00297 int timeStep = 0;
00298 for (std::vector<float>::iterator it = volumes.begin(); it != volumes.end(); it++)
00299 {
00300 vs << timeStep++ << "\t" << *it << std::endl;
00301 }
00302 m_Controls->m_TextEdit->setText(vs.str().c_str());
00303 m_Controls->m_TextEdit->setTabStopWidth(20);
00304 m_Controls->m_SaveCsvButton->setEnabled(true);
00305 }
00306 }
00307 }
00308
00309 const mitk::DataNode* QmitkVolumetryView::GetImageNode()
00310 {
00311 return m_SelectedDataNode;
00312 }
00313
00314 void QmitkVolumetryView::UpdateSlider()
00315 {
00316 if (m_SelectedDataNode.IsNotNull() && dynamic_cast<mitk::Image*>(m_SelectedDataNode->GetData()))
00317 {
00318 int intSliderValue = (int)m_Controls->m_ThresholdSlider->value();
00319 QString stringSliderValue;
00320 stringSliderValue.setNum(intSliderValue);
00321 m_Controls->m_ThresholdLineEdit->setText(stringSliderValue);
00322 }
00323 }
00324
00325
00326 void QmitkVolumetryView::UpdateSliderLabel()
00327 {
00328 int sliderValue = atoi(m_Controls->m_ThresholdLineEdit->text().toLatin1());
00329 m_Controls->m_ThresholdSlider->setValue(sliderValue);
00330 this->UpdateSlider();
00331
00332 }
00333
00334 void QmitkVolumetryView::OnThresholdSliderChanged( int value)
00335 {
00336 if (m_OverlayNode)
00337 {
00338 m_OverlayNode->SetLevelWindow(mitk::LevelWindow(value,1));
00339 mitk::RenderingManager::GetInstance()->RequestUpdateAll();
00340
00341 this->UpdateSlider();
00342 }
00343 }
00344
00345 void QmitkVolumetryView::CreateOverlayChild()
00346 {
00347 if (m_SelectedDataNode.IsNotNull())
00348 {
00349 m_OverlayNode = mitk::DataNode::New();
00350 mitk::StringProperty::Pointer nameProp = mitk::StringProperty::New("volume threshold overlay image" );
00351 m_OverlayNode->SetProperty( "reslice interpolation", m_SelectedDataNode->GetProperty("reslice interpolation") );
00352 m_OverlayNode->SetProperty( "name", nameProp );
00353 m_OverlayNode->SetData(m_SelectedDataNode->GetData());
00354 m_OverlayNode->SetColor(0.0,1.0,0.0);
00355 m_OverlayNode->SetOpacity(.25);
00356 int layer = 0;
00357 m_SelectedDataNode->GetIntProperty("layer", layer);
00358 m_OverlayNode->SetIntProperty("layer", layer+1);
00359 m_OverlayNode->SetLevelWindow(mitk::LevelWindow(m_Controls->m_ThresholdSlider->value(),1));
00360
00361
00362 this->GetDefaultDataStorage()->Add(m_OverlayNode);
00363
00365
00366
00367
00368
00369
00370
00371
00372
00373
00374
00375
00376 }
00377 }
00378
00379
00380 mitk::DataNode* QmitkVolumetryView::GetOverlayNode()
00381 {
00382 return m_OverlayNode;
00383 }
00384
00385
00386
00387
00388 mitk::Image* QmitkVolumetryView::GetImage()
00389 {
00390 return dynamic_cast<mitk::Image*>(m_SelectedDataNode->GetData());
00391 }
00392 void QmitkVolumetryView::OnSaveCsvButtonClicked()
00393 {
00394 static QString lastSavePath = QDir::homePath();
00395
00396
00397 QString s = QFileDialog::getSaveFileName( this->m_Parent, "Save as..", lastSavePath, "CSV Files (*.csv)");
00398
00399
00400
00401
00402 if (! s.isEmpty())
00403 {
00404 lastSavePath = s;
00405 QFile saveFile(s);
00406
00407 if ( saveFile.open(QIODevice::WriteOnly))
00408 {
00409 QTextStream stream( &saveFile );
00410 stream << m_Controls->m_TextEdit->toPlainText().replace('\t',';');
00411 saveFile.close();
00412 }
00413 else
00414 {
00415
00416
00417 }
00418 }
00419 }
00420
00421
00422
00423
00424