]> Creatis software - clitk.git/blob - vv/vvSegmentationDialog.cxx
remove debug messages
[clitk.git] / vv / vvSegmentationDialog.cxx
1 /*=========================================================================
2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
3
4   Authors belong to:
5   - University of LYON              http://www.universite-lyon.fr/
6   - Léon Bérard cancer center       http://oncora1.lyon.fnclcc.fr
7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
8
9   This software is distributed WITHOUT ANY WARRANTY; without even
10   the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
11   PURPOSE.  See the copyright notices for more information.
12
13   It is distributed under dual licence
14
15   - BSD        See included LICENSE.txt file
16   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ======================================================================-====*/
18 #ifndef _vvSegmentationDialog_CXX
19 #define _vvSegmentationDialog_CXX
20 #include <QtGui>
21 #include <Qt>
22
23 #include "vvSegmentationDialog.h"
24 #include "vvProgressDialog.h"
25 #include "vvImageWriter.h"
26 #include "vvLandmarks.h"
27 #include "vvInteractorStyleNavigator.h"
28 #include "vvSlicer.h"
29
30 #include "vtkMarchingCubes.h"
31 #include "vtkMarchingSquares.h"
32 #include "vtkImageClip.h"
33 #include "vtkCamera.h"
34 #include "vtkRenderer.h"
35 #include "vtkProperty.h"
36 #include "vtkLookupTable.h"
37 #include "vtkClipPolyData.h"
38 #include "vtkImageToPolyDataFilter.h"
39 #include "vtkLookupTable.h"
40 #include "vtkDoubleArray.h"
41 #include "vtkPointData.h"
42 #include "vtkCellData.h"
43 #include "vtkImageMapToWindowLevelColors.h"
44 #include "vtkImageContinuousErode3D.h"
45 #include "vtkImageContinuousDilate3D.h"
46 #include "vtkImageLogic.h"
47 #include "vtkInteractorStyleTrackballCamera.h"
48 #include "vtkImageSeedConnectivity.h"
49 #include "vtkConnectivityFilter.h"
50 #include "vtkPolyData.h"
51 #include <vtkPolyDataMapper.h>
52 #include <vtkImageData.h>
53 #include "vtkInformation.h"
54 #include "vtkInformationVector.h"
55 #include "vtkStreamingDemandDrivenPipeline.h"
56 #include <vtkPolyDataWriter.h>
57
58 #include <QMessageBox>
59
60 //====================================================================
61 vvSegmentationDialog::vvSegmentationDialog(QWidget * parent, Qt::WindowFlags f)
62   :QDialog(parent,f), Ui::vvSegmentationDialog()
63 {
64
65   // initialization
66   setupUi(this);
67   mManager = new vvSlicerManager(1);
68
69   mClipper = vtkImageClip::New();
70   mSquares1 = vtkMarchingSquares::New();
71   mSquaresMapper1 = vtkPolyDataMapper::New();
72   mSquaresActor1 = vtkActor::New();
73
74   mSquares2 = vtkMarchingSquares::New();
75   mSquaresMapper2 = vtkPolyDataMapper::New();
76   mSquaresActor2 = vtkActor::New();
77
78   //m3DMapper = vtkPolyDataMapper::New();
79   //m3DActor = vtkActor::New();
80   m3DExtractor = vtkMarchingCubes::New();
81   m3DExtractor->ComputeScalarsOff();
82   m3DMappers.clear();
83   m3DActors.clear();
84
85   mBinaireImages.clear();
86   mKernelValue = 2;
87
88   connect(clipping1Slider,SIGNAL(valueChanged(int)),this,SLOT(clippingvaluechanged(int)));
89   connect(clipping2Slider,SIGNAL(valueChanged(int)),this,SLOT(clippingvaluechanged(int)));
90   connect(binaryButton,SIGNAL(clicked()),this,SLOT(BinariseSurface()));
91   connect(saveButton,SIGNAL(clicked()),this,SLOT(Save()));
92   connect(erodeButton,SIGNAL(clicked()),this,SLOT(Erode()));
93   connect(dilateButton,SIGNAL(clicked()),this,SLOT(Dilate()));
94   connect(dimButton,SIGNAL(clicked()),this,SLOT(ChangeDimRendering()));
95   connect(kernelSpinBox,SIGNAL(valueChanged(int)),this,SLOT(KernelValueChanged(int)));
96
97   binaryButton->setEnabled(0);
98   erodeButton->setEnabled(0);
99   dilateButton->setEnabled(0);
100   infoLabel->setText("Select Up and Down threshold before clicking binarise !");
101 }
102
103 vvSegmentationDialog::~vvSegmentationDialog()
104 {
105   mClipper->Delete();
106
107   mSquaresActor1->Delete();
108   mSquaresMapper1->Delete();
109   mSquares1->Delete();
110
111   mSquaresActor2->Delete();
112   mSquaresMapper2->Delete();
113   mSquares2->Delete();
114
115   //m3DMapper->Delete();
116   //m3DActor->Delete();
117   m3DExtractor->Delete();
118
119   for (unsigned int i = 0; i < mBinaireImages.size(); i++)
120     mBinaireImages[i]->Delete();
121
122   for (unsigned int i = 0; i < m3DActors.size(); i++)
123     m3DActors[i]->Delete();
124
125   for (unsigned int i = 0; i < m3DMappers.size(); i++)
126     m3DMappers[i]->Delete();
127
128   delete mManager;
129 }
130
131 //----------------------------------------------------------------------------
132 // This templated function executes the filter for any type of data.
133 // Handles the one input operations
134 template <class T>
135 void vvImageBinarize(vtkImageData *in1Data, T *in1Ptr,
136                      int outExt[6],int clampMin, int clampMax)
137 {
138   int idxR, idxY, idxZ;
139   int maxY, maxZ;
140   vtkIdType inIncX, inIncY, inIncZ;
141   int rowLength;
142
143   // find the region to loop over
144   rowLength =
145     (outExt[1] - outExt[0]+1)*in1Data->GetNumberOfScalarComponents();
146   // What a pain. Maybe I should just make another filter.
147
148   maxY = outExt[3] - outExt[2];
149   maxZ = outExt[5] - outExt[4];
150
151   // Get increments to march through data
152   in1Data->GetContinuousIncrements(outExt, inIncX, inIncY, inIncZ);
153
154   for (idxZ = 0; idxZ <= maxZ; idxZ++) {
155     for (idxY = 0; idxY <= maxY; idxY++) {
156       for (idxR = 0; idxR < rowLength; idxR++) {
157         if (static_cast<double>(*in1Ptr) > clampMin && static_cast<double>(*in1Ptr) <= clampMax)
158           *in1Ptr = static_cast<T>(1);
159         else
160           *in1Ptr = static_cast<T>(0);
161         in1Ptr++;
162       }
163       in1Ptr += inIncY;
164     }
165     in1Ptr += inIncZ;
166   }
167 }
168
169 void vvSegmentationDialog::SetImage(vvImage::Pointer image)
170 {
171
172   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
173   mManager->SetImage(image);
174   mManager->SetSlicerWindow(0,viewWidget->GetRenderWindow());
175   vvInteractorStyleNavigator* style = vvInteractorStyleNavigator::New();
176   mManager->SetInteractorStyleNavigator(0,style);
177   style->Delete();
178
179   double range[2];
180   mManager->GetImage()->GetScalarRange(range);
181   mManager->GetSlicer(0)->SetColorWindow(range[1]-range[0]);
182   mManager->GetSlicer(0)->SetColorLevel((range[1]+range[0])/2);
183
184   clipping1Slider->setMinimum(range[0]);
185   clipping1Slider->setMaximum(range[1]);
186   clipping2Slider->setMinimum(range[0]);
187   clipping2Slider->setMaximum(range[1]);
188   clipping1Slider->setValue(range[0]);
189   clipping2Slider->setValue(range[1]);
190
191   mClipper->SetInput(mManager->GetSlicer(0)->GetInput());
192   mSquares1->SetValue(0,clipping1Slider->value());
193   mSquares2->SetValue(0,clipping2Slider->value());
194   mSquares1->SetInput(mClipper->GetOutput());
195   mSquares2->SetInput(mClipper->GetOutput());
196
197   mSquaresMapper1->SetInput(mSquares1->GetOutput());
198   mSquaresMapper2->SetInput(mSquares2->GetOutput());
199   mSquaresMapper1->ScalarVisibilityOff();
200   mSquaresMapper2->ScalarVisibilityOff();
201
202   mSquaresActor1->SetMapper(mSquaresMapper1);
203   mSquaresActor2->SetMapper(mSquaresMapper2);
204   mSquaresActor1->GetProperty()->SetColor(1.0,0,0);
205   mSquaresActor2->GetProperty()->SetColor(0,0,1.0);
206   mSquaresActor1->SetPickable(0);
207   mSquaresActor2->SetPickable(0);
208
209   mManager->GetSlicer(0)->GetRenderer()->AddActor(mSquaresActor1);
210   mManager->GetSlicer(0)->GetRenderer()->AddActor(mSquaresActor2);
211
212   mSquares1->Update();
213   mSquares2->Update();
214
215   UpdateSlice(0,mManager->GetSlicer(0)->GetSlice());
216
217   connect(mManager,SIGNAL(UpdateTSlice(int,int)),this,SLOT(UpdateSlice(int, int)));
218   connect(mManager,SIGNAL(UpdateSlice(int,int)),this,SLOT(UpdateSlice(int, int)));
219   connect(mManager,SIGNAL(UpdateSliceRange(int,int,int,int,int)),this,SLOT(UpdateSlice(int, int)));
220   connect(mManager,SIGNAL(LandmarkAdded()),this,SLOT(InsertSeed()));
221   QApplication::restoreOverrideCursor();
222 }
223
224 void vvSegmentationDialog::UpdateSlice(int slicer,int slices)
225 {
226   int slice = mManager->GetSlicer(0)->GetSlice();
227   int tslice = mManager->GetSlicer(0)->GetTSlice();
228   mClipper->SetInput(mManager->GetSlicer(0)->GetInput());
229   int* extent = mManager->GetSlicer(0)->GetImageActor()->GetDisplayExtent();
230   mClipper->SetOutputWholeExtent(extent[0],extent[1],extent[2],extent[3],extent[4],extent[5]);
231   int i;
232   for (i = 0; i < 6; i = i+2) {
233     if (extent[i] == extent[i+1]) {
234       break;
235     }
236   }
237
238   switch (i) {
239   case 0:
240     if (mManager->GetSlicer(0)->GetRenderer()->GetActiveCamera()->GetPosition()[0] > slice) {
241       mSquaresActor1->SetPosition(1,0,0);
242       mSquaresActor2->SetPosition(1,0,0);
243     } else {
244       mSquaresActor1->SetPosition(-1,0,0);
245       mSquaresActor2->SetPosition(-1,0,0);
246     }
247     break;
248   case 2:
249     if (mManager->GetSlicer(0)->GetRenderer()->GetActiveCamera()->GetPosition()[1] > slice) {
250       mSquaresActor1->SetPosition(0,1,0);
251       mSquaresActor2->SetPosition(0,1,0);
252     } else {
253       mSquaresActor1->SetPosition(0,-1,0);
254       mSquaresActor2->SetPosition(0,-1,0);
255     }
256     break;
257   case 4:
258     if (mManager->GetSlicer(0)->GetRenderer()->GetActiveCamera()->GetPosition()[2] > slice) {
259       mSquaresActor1->SetPosition(0,0,1);
260       mSquaresActor2->SetPosition(0,0,1);
261     } else {
262       mSquaresActor1->SetPosition(0,0,-1);
263       mSquaresActor2->SetPosition(0,0,-1);
264     }
265     break;
266   }
267   mSquares1->Update();
268   mSquares2->Update();
269
270   if (m3DActors.size()) {
271     for (unsigned int i =0; i < m3DActors.size(); i++) {
272       if (m3DActors[i]->GetVisibility()) {
273         m3DActors[i]->VisibilityOff();
274       }
275     }
276     std::cout << "display " << tslice << " on " << m3DActors.size()  << std::endl;
277     m3DActors[tslice]->VisibilityOn();
278   }
279
280   mManager->Render();
281 }
282
283
284 void vvSegmentationDialog::clippingvaluechanged(int value)
285 {
286   binaryButton->setEnabled(1);
287   int min = (clipping1Slider->value() < clipping2Slider->value() ) ?
288             clipping1Slider->value():clipping2Slider->value();
289   int max = (clipping1Slider->value() > clipping2Slider->value() ) ?
290             clipping1Slider->value():clipping2Slider->value();
291   mSquares1->SetValue(0,min);
292   mSquares2->SetValue(0,max);
293
294   QString textMin = "<b> Min : </b>";
295   textMin += QString::number(min);
296   QString textMax = "\n <b> Max : </b>";
297   textMax += QString::number(max);
298   minLabel->setText(textMin);
299   maxLabel->setText(textMax);
300
301   if (mSquares1->GetInput()) {
302     mSquares1->Update();
303     mSquares2->Update();
304     mManager->Render();
305   }
306 }
307
308
309 void vvSegmentationDialog::BinariseSurface()
310 {
311   infoLabel->setText("Click erode then space on desired organ !");
312
313   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
314   int clampMin = (clipping1Slider->value() < clipping2Slider->value() ) ?
315                  clipping1Slider->value():clipping2Slider->value();
316   int clampMax = (clipping1Slider->value() > clipping2Slider->value() ) ?
317                  clipping1Slider->value():clipping2Slider->value();
318   vtkImageData* outputImage = vtkImageData::New();
319
320   for (unsigned int numImage = 0; numImage < mManager->GetSlicer(0)->GetImage()->GetVTKImages().size(); numImage++) {
321     vtkImageData* image = mManager->GetSlicer(0)->GetImage()->GetVTKImages()[numImage];
322     int ext[6];
323     image->GetWholeExtent(ext);
324     void *in1Ptr;
325     in1Ptr = image->GetScalarPointerForExtent(ext);
326
327     switch (image->GetScalarType()) {
328       vtkTemplateMacro(
329         vvImageBinarize(image, static_cast<VTK_TT *>(in1Ptr),
330                         ext,clampMin,clampMax));
331     default:
332       std::cerr << "Error, unknown pixel format : " << image->GetScalarTypeAsString() << std::endl;
333       return;
334     }
335
336     outputImage->Initialize();
337     outputImage->SetExtent(ext);
338     outputImage->SetOrigin(image->GetOrigin());
339     outputImage->SetSpacing(image->GetSpacing());
340     outputImage->SetScalarTypeToUnsignedChar();
341     outputImage->CopyAndCastFrom(image,ext);
342     outputImage->Update();
343
344     image->DeepCopy(outputImage);
345     image->UpdateInformation();
346     image->PropagateUpdateExtent();
347
348     vtkImageData* imageBin = vtkImageData::New();
349     imageBin->DeepCopy(image);
350     imageBin->Update();
351     mBinaireImages.push_back(imageBin);
352   }
353
354   outputImage->Delete();
355   erodeButton->setEnabled(1);
356   QApplication::restoreOverrideCursor();
357   mManager->SetColorWindow(2);
358   mManager->SetColorLevel(0.5);
359   mManager->Render();
360 }
361
362 void vvSegmentationDialog::Erode()
363 {
364   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
365   vtkImageContinuousErode3D* erode = vtkImageContinuousErode3D::New();
366   erode->SetKernelSize(mKernelValue,mKernelValue,mKernelValue);
367   for (unsigned int numImage = 0; numImage < mManager->GetSlicer(0)->GetImage()->GetVTKImages().size(); numImage++) {
368     vtkImageData* image = mManager->GetSlicer(0)->GetImage()->GetVTKImages()[numImage];
369     erode->SetInput(image);
370     erode->Update();
371     image->DeepCopy(erode->GetOutput());
372     image->Update();
373   }
374   erode->Delete();
375   dilateButton->setEnabled(1);
376   mManager->Render();
377   QApplication::restoreOverrideCursor();
378 }
379
380 void vvSegmentationDialog::Dilate()
381 {
382   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
383   vtkImageContinuousDilate3D* dilate = vtkImageContinuousDilate3D::New();
384   vtkImageLogic* And = vtkImageLogic::New();
385   And->SetOperationToAnd();
386   dilate->SetKernelSize(mKernelValue,mKernelValue,mKernelValue);
387   for (unsigned int numImage = 0; numImage < mManager->GetSlicer(0)->GetImage()->GetVTKImages().size(); numImage++) {
388     vtkImageData* image = mManager->GetSlicer(0)->GetImage()->GetVTKImages()[numImage];
389     dilate->SetInput(image);
390     vtkImageData* mask = mBinaireImages[numImage];
391     And->SetInput1(dilate->GetOutput());
392     And->SetInput2(mask);
393     And->Update();
394     image->DeepCopy(And->GetOutput());
395     image->Update();
396   }
397   And->Delete();
398   dilate->Delete();
399   mManager->Render();
400   QApplication::restoreOverrideCursor();
401 }
402
403 void vvSegmentationDialog::InsertSeed()
404 {
405   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
406   int point4D[4];
407   point4D[0] = mManager->GetLandmarks()->GetCoordinates(
408                  mManager->GetLandmarks()->GetNumberOfPoints()-1)[0];
409   point4D[1] = mManager->GetLandmarks()->GetCoordinates(
410                  mManager->GetLandmarks()->GetNumberOfPoints()-1)[1];
411   point4D[2] = mManager->GetLandmarks()->GetCoordinates(
412                  mManager->GetLandmarks()->GetNumberOfPoints()-1)[2];
413   point4D[3] = mManager->GetLandmarks()->GetCoordinates(
414                  mManager->GetLandmarks()->GetNumberOfPoints()-1)[3];
415
416   point4D[0] = (point4D[0]-mManager->GetSlicer(0)->GetImage()->GetVTKImages()[0]->GetOrigin()[0])/mManager->GetSlicer(0)->GetImage()->GetVTKImages()[0]->GetSpacing()[0];
417   point4D[1] = (point4D[1]-mManager->GetSlicer(0)->GetImage()->GetVTKImages()[0]->GetOrigin()[1])/mManager->GetSlicer(0)->GetImage()->GetVTKImages()[0]->GetSpacing()[1];
418   point4D[2] = (point4D[2]-mManager->GetSlicer(0)->GetImage()->GetVTKImages()[0]->GetOrigin()[2])/mManager->GetSlicer(0)->GetImage()->GetVTKImages()[0]->GetSpacing()[2];
419
420   vtkImageSeedConnectivity* seed = vtkImageSeedConnectivity::New();
421   seed->SetInputConnectValue(1);
422   seed->SetOutputConnectedValue(1);
423   seed->SetOutputUnconnectedValue(0);
424   seed->AddSeed(point4D[0],point4D[1],point4D[2]);
425
426   for (unsigned int numImage = 0; numImage < mManager->GetSlicer(0)->GetImage()->GetVTKImages().size(); numImage++) {
427     vtkImageData* image = mManager->GetSlicer(0)->GetImage()->GetVTKImages()[numImage];
428     seed->SetInput(image);
429     seed->Update();
430     image->DeepCopy(seed->GetOutput());
431     image->Update();
432   }
433
434   seed->Delete();
435   QApplication::restoreOverrideCursor();
436 }
437
438 void vvSegmentationDialog::ChangeDimRendering()
439 {
440   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
441
442   if (dimButton->text() == "3D") {
443     if (m3DActors.size() == 0) {
444       m3DExtractor->SetValue(0,0.5);
445       for (unsigned int numImage = 0; numImage < mManager->GetSlicer(0)->GetImage()->GetVTKImages().size(); numImage++) {
446         vtkActor* actor = vtkActor::New();
447         m3DExtractor->SetInput(mManager->GetSlicer(0)->GetImage()->GetVTKImages()[numImage]);
448         m3DExtractor->Update();
449
450         vtkPolyDataMapper* mapper = vtkPolyDataMapper::New();
451         mapper->SetInput(m3DExtractor->GetOutput());
452         m3DMappers.push_back(mapper);
453
454         actor->SetMapper(mapper);
455         actor->GetProperty()->SetColor(1.0,0.7,0.2);
456         actor->VisibilityOff();
457
458         mManager->GetSlicer(0)->GetRenderer()->AddActor(actor);
459         m3DActors.push_back(actor);
460       }
461     }
462
463     mManager->GetSlicer(0)->GetRenderer()->SetBackground(0.5,0.6,0.9);
464     m3DActors[0]->VisibilityOn();
465
466     vtkInteractorStyleTrackballCamera* style = vtkInteractorStyleTrackballCamera::New();
467     mManager->SetInteractorStyleNavigator(0,style);
468     style->Delete();
469
470     mManager->GetSlicer(0)->GetImageActor()->VisibilityOff();
471     mSquaresActor1->VisibilityOff();
472     mSquaresActor2->VisibilityOff();
473     mManager->Render();
474     dimButton->setText("2D");
475   } else {
476     mManager->GetSlicer(0)->GetRenderer()->SetBackground(0.0,0.0,0.0);
477     vvInteractorStyleNavigator* style = vvInteractorStyleNavigator::New();
478     mManager->SetInteractorStyleNavigator(0,style);
479     style->Delete();
480
481     mManager->GetSlicer(0)->SetSliceOrientation(2);
482     m3DActors[mManager->GetSlicer(0)->GetTSlice()]->VisibilityOff();
483
484     mManager->GetSlicer(0)->GetImageActor()->VisibilityOn();
485     mSquaresActor1->VisibilityOn();
486     mSquaresActor2->VisibilityOn();
487     dimButton->setText("3D");
488   }
489   QApplication::restoreOverrideCursor();
490 }
491
492 void vvSegmentationDialog::KernelValueChanged(int kernel)
493 {
494   mKernelValue = kernel;
495 }
496
497 void vvSegmentationDialog::Save()
498 {
499   if (dimButton->text() == "2D") { //If in *3D* mode, save the mesh
500     QString fileName = QFileDialog::getSaveFileName(this,
501                        tr("Save Mesh As"),
502                        QDir::home().dirName(),
503                        "Mesh Files (*.vtk *.vtp)");
504     if (!fileName.isEmpty()) {
505       vtkSmartPointer<vtkPolyDataWriter> w = vtkSmartPointer<vtkPolyDataWriter>::New();
506       w->SetInput(m3DExtractor->GetOutput());
507       w->SetFileName(fileName.toStdString().c_str());
508       w->Write();
509     }
510   } else {
511     QStringList OutputListeFormat;
512     OutputListeFormat.clear();
513     int dimension = mManager->GetDimension();
514     if (dimension == 1) {
515       OutputListeFormat.push_back(".mhd");
516     }
517     if (dimension == 2) {
518       OutputListeFormat.push_back(".bmp");
519       OutputListeFormat.push_back(".png");
520       OutputListeFormat.push_back(".jpeg");
521       OutputListeFormat.push_back(".tif");
522       OutputListeFormat.push_back(".mhd");
523       OutputListeFormat.push_back(".hdr");
524       OutputListeFormat.push_back(".vox");
525     } else if (dimension == 3) {
526       OutputListeFormat.push_back(".mhd");
527       OutputListeFormat.push_back(".hdr");
528       OutputListeFormat.push_back(".vox");
529     } else if (dimension == 4) {
530       OutputListeFormat.push_back(".mhd");
531     }
532     QString Extensions = "AllFiles(*.*)";
533     for (int i = 0; i < OutputListeFormat.count(); i++) {
534       Extensions += ";;Images ( *";
535       Extensions += OutputListeFormat[i];
536       Extensions += ")";
537     }
538     QString fileName = QFileDialog::getSaveFileName(this,
539                        tr("Save As"),
540                        QDir::home().dirName(),
541                        Extensions);
542     if (!fileName.isEmpty()) {
543       std::string fileformat = vtksys::SystemTools::GetFilenameLastExtension(fileName.toStdString());
544       if (OutputListeFormat.contains(
545             fileformat.c_str())) {
546         QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
547         vvProgressDialog progress("Saving "+fileName.toStdString());
548         qApp->processEvents();
549         vvImageWriter *writer = new vvImageWriter;
550         writer->SetOutputFileName(fileName.toStdString());
551         writer->SetInput(mManager->GetSlicer(0)->GetImage());
552         writer->Update(dimension,"unsigned_char");
553         QApplication::restoreOverrideCursor();
554         if (writer->GetLastError().size()) {
555           QString error = "Saving did not succeed\n";
556           error += writer->GetLastError().c_str();
557           QMessageBox::information(this,tr("Saving Problem"),error);
558           Save();
559         }
560       } else {
561         QString error = fileformat.c_str();
562         if (error.isEmpty())
563           error += "no file format specified !";
564         else
565           error += " format unknown !!!\n";
566         QMessageBox::information(this,tr("Saving Problem"),error);
567         Save();
568       }
569     }
570   }
571 }
572
573 #endif /* end #define _vvSegmentationDialog_CXX */
574