1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
5 - University of LYON http://www.universite-lyon.fr/
6 - Léon Bérard cancer center http://www.centreleonberard.fr
7 - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
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.
13 It is distributed under dual licence
15 - BSD See included LICENSE.txt file
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ===========================================================================**/
20 #include "vvToolCropImage.h"
24 #include "clitkCropImage_ggo.h"
25 #include "clitkCropImageGenericFilter.h"
26 #include "clitkExceptionObject.h"
31 #include <QApplication>
32 #include <QMessageBox>
33 #include <QSignalMapper>
36 #include <vtkVersion.h>
37 #include <vtkStreamingDemandDrivenPipeline.h>
38 #include <vtkInformation.h>
39 #include <vtkImageClip.h>
40 #include <vtkImageTranslateExtent.h>
41 #include <vtkImageData.h>
42 #include <vtkSmartPointer.h>
44 //------------------------------------------------------------------------------
45 // Create the tool and automagically (I like this word) insert it in
46 // the main window menu.
47 ADD_TOOL(vvToolCropImage);
48 //------------------------------------------------------------------------------
50 //same order of int* returned by VtkImageData::WholeExtent
51 enum sliderId {xmin, xmax, ymin, ymax, zmin, zmax, tmin, tmax};
52 //------------------------------------------------------------------------------
53 vvToolCropImage::vvToolCropImage(vvMainWindowBase * parent, Qt::WindowFlags f):
54 vvToolWidgetBase(parent, f),
55 vvToolBase<vvToolCropImage>(parent),
56 Ui::vvToolCropImage(),mSliders(8)
60 Ui_vvToolCropImage::setupUi(mToolWidget);
61 mTLabel1->setHidden(true);
62 mTLabel2->setHidden(true);
63 tminSlider->setHidden(true);
64 tmaxSlider->setHidden(true);
65 spin_tmin->setHidden(true);
66 spin_tmax->setHidden(true);
67 mLabelTimeCropping->setHidden(true);
69 // Set how many inputs are needed for this tool
70 AddInputSelector("Select one image");
72 mSliders[xmin]=xminSlider;
73 mSliders[xmax]=xmaxSlider;
74 mSliders[ymin]=yminSlider;
75 mSliders[ymax]=ymaxSlider;
76 mSliders[zmin]=zminSlider;
77 mSliders[zmax]=zmaxSlider;
78 mSliders[tmin]=tminSlider;
79 mSliders[tmax]=tmaxSlider;
81 // Record initial extend
82 mReducedExtent = new int[6];
83 mInitialExtent = new int[6];
85 //------------------------------------------------------------------------------
88 //------------------------------------------------------------------------------
89 vvToolCropImage::~vvToolCropImage()
91 delete [] mReducedExtent;
92 delete [] mInitialExtent;
94 //------------------------------------------------------------------------------
97 //------------------------------------------------------------------------------
98 void vvToolCropImage::closeEvent(QCloseEvent *event)
100 if(mCurrentSlicerManager){
102 for(int i=0; i<6; i++){
103 mReducedExtent[i] = mInitialExtent[i];
105 for(int i=0; i<mCurrentSlicerManager->GetNumberOfSlicers(); i++)
106 mCurrentSlicerManager->GetSlicer(i)->EnableReducedExtent(false);
109 mCurrentSlicerManager->GetImage()->GetTransform()[0]->SetMatrix(mConcatenedTransform);
110 for (int i=0; i<mCurrentSlicerManager->GetNumberOfSlicers(); i++) {
111 mCurrentSlicerManager->GetSlicer(i)->ResetCamera();
112 mCurrentSlicerManager->GetSlicer(i)->Render();
113 mCurrentSlicerManager->UpdateLinkedNavigation( mCurrentSlicerManager->GetSlicer(i) );
115 vvToolWidgetBase::closeEvent(event);
117 //------------------------------------------------------------------------------
120 //------------------------------------------------------------------------------
121 bool vvToolCropImage::close()
123 return vvToolWidgetBase::close();
125 //------------------------------------------------------------------------------
128 //------------------------------------------------------------------------------
129 void vvToolCropImage::reject()
131 for(int i=0; i<mExtentSize; i++)
132 mReducedExtent[i] = mInitialExtent[i];
133 for(int i=0; i<mCurrentSlicerManager->GetNumberOfSlicers(); i++)
134 mCurrentSlicerManager->GetSlicer(i)->EnableReducedExtent(false);
136 return vvToolWidgetBase::reject();
138 //------------------------------------------------------------------------------
141 //------------------------------------------------------------------------------
142 void vvToolCropImage::sliderValueChanged(int dim)
145 if(dim%2==0){//case we are minimum
146 mSliders[dim+1]->setMinimum(mSliders[dim]->value());
148 mSliders[--dimMin]->setMaximum(mSliders[dim]->value());
150 mReducedExtent[dim] = mSliders[dim]->value() + mInitialExtent[dimMin];
153 //------------------------------------------------------------------------------
156 //------------------------------------------------------------------------------
157 void vvToolCropImage::UpdateExtent()
159 for(int i=0; i<mCurrentSlicerManager->GetNumberOfSlicers(); i++) {
160 mCurrentSlicerManager->GetSlicer(i)->SetReducedExtent(mReducedExtent);
161 mCurrentSlicerManager->GetSlicer(i)->ForceUpdateDisplayExtent();
162 mCurrentSlicerManager->GetSlicer(i)->Render();
165 //------------------------------------------------------------------------------
168 //------------------------------------------------------------------------------
169 void vvToolCropImage::InputIsSelected(vvSlicerManager * slicer)
171 //Save the current transformation
172 mConcatenedTransform = vtkSmartPointer<vtkMatrix4x4>::New();
173 mConcatenedTransform->DeepCopy(slicer->GetSlicer(0)->GetConcatenatedTransform()->GetMatrix());
174 vtkSmartPointer<vtkMatrix4x4> matrix = vtkSmartPointer<vtkMatrix4x4>::New();
176 mCurrentSlicerManager->GetImage()->GetTransform()[0]->SetMatrix(matrix);
177 for (int i=0; i<mCurrentSlicerManager->GetNumberOfSlicers(); i++) {
178 mCurrentSlicerManager->GetSlicer(i)->ResetCamera();
179 mCurrentSlicerManager->GetSlicer(i)->Render();
180 mCurrentSlicerManager->UpdateLinkedNavigation( mCurrentSlicerManager->GetSlicer(i) );
183 // Change interface according to number of dimension
184 mExtentSize = 2*slicer->GetDimension();
185 if (slicer->GetDimension()<4) {
186 mTLabel1->setHidden(true);
187 mTLabel2->setHidden(true);
188 tminSlider->setHidden(true);
189 tmaxSlider->setHidden(true);
190 spin_tmin->setHidden(true);
191 spin_tmax->setHidden(true);
192 mLabelTimeCropping->setHidden(true);
194 if (slicer->GetDimension()<3) {
195 mZLabel1->setHidden(true);
196 mZLabel2->setHidden(true);
197 zminSlider->setHidden(true);
198 zmaxSlider->setHidden(true);
199 spin_zmin->setHidden(true);
200 spin_zmax->setHidden(true);
202 #if VTK_MAJOR_VERSION <= 5
203 int *a = mCurrentImage->GetFirstVTKImageData()->GetWholeExtent();
205 int *a = mCurrentImage->GetFirstVTKImageData()->GetInformation()->Get(vtkDataObject::DATA_EXTENT());
207 for(int i=0; i<6; i++){
208 mInitialExtent[i] = a[i];
209 mReducedExtent[i] = a[i];
211 for(int i=0; i<mCurrentSlicerManager->GetNumberOfSlicers(); i++) {
212 mCurrentSlicerManager->GetSlicer(i)->EnableReducedExtent(true);
213 mCurrentSlicerManager->GetSlicer(i)->SetReducedExtent(mInitialExtent);
216 // Set initial sliders values
217 int w_ext[6], imsize[3];
218 mCurrentSlicerManager->GetSlicer(0)->GetRegisterExtent(w_ext);
219 for(int dim=0; dim<slicer->GetDimension() && dim<3; ++dim){
220 imsize[dim] = w_ext[2*dim+1] - w_ext[2*dim] +1;
221 mSliders[dim*2]->setMaximum(imsize[dim]-1);
222 mSliders[dim*2+1]->setMaximum(imsize[dim]-1);
223 mSliders[dim*2+1]->setValue(imsize[dim]-1);
225 spin_xmin->setMaximum(imsize[0]-1);
226 spin_xmax->setMaximum(imsize[0]-1);
227 spin_xmax->setValue(imsize[0]-1);
229 spin_ymin->setMaximum(imsize[1]-1);
230 spin_ymax->setMaximum(imsize[1]-1);
231 spin_ymax->setValue(imsize[1]-1);
233 if (slicer->GetDimension() >2) {
234 spin_zmin->setMaximum(imsize[2]-1);
235 spin_zmax->setMaximum(imsize[2]-1);
236 spin_zmax->setValue(imsize[2]-1);
239 if (slicer->GetDimension() >3) {
240 spin_tmin->setMaximum(imsize[3]-1);
241 spin_tmax->setMaximum(imsize[3]-1);
242 spin_tmax->setValue(imsize[3]-1);
245 QSignalMapper* signalMapper = new QSignalMapper(this);
246 connect(signalMapper, SIGNAL(mapped(int)), this, SLOT(sliderValueChanged(int)));
247 for(unsigned int i=0; i<mSliders.size(); ++i){
248 signalMapper->setMapping(mSliders[i], i);
249 connect(mSliders[i], SIGNAL(valueChanged(int)), signalMapper, SLOT(map()));
253 //------------------------------------------------------------------------------
256 //------------------------------------------------------------------------------
257 void vvToolCropImage::autoCropValueChanged(double v)
261 //------------------------------------------------------------------------------
264 //------------------------------------------------------------------------------
265 void vvToolCropImage::apply()
267 if (!mCurrentSlicerManager) close();
268 QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
271 typedef args_info_clitkCropImage ArgsInfoType;
272 typedef clitk::CropImageGenericFilter CropFilterType;
275 ArgsInfoType mArgsInfo;
276 cmdline_parser_clitkCropImage_init(&mArgsInfo); // Initialisation to default
277 int n = mCurrentSlicerManager->GetDimension()*2; // 2D and 3D only
278 mArgsInfo.boundingBox_given = n;
279 mArgsInfo.boundingBox_arg = new int[n];
281 for(int dim=0; dim<mCurrentSlicerManager->GetDimension() && dim<3; ++dim){
282 mArgsInfo.boundingBox_arg[dim*2] = mSliders[dim*2]->value();
283 mArgsInfo.boundingBox_arg[dim*2+1] = mSliders[dim*2+1]->value();
286 mArgsInfo.boundingBox_arg[6] = 0;
287 mArgsInfo.boundingBox_arg[7] = mCurrentImage->GetSize()[3]-1;
289 // We MUST reset initial extend to input image before using the
290 // filter to retrieve the correct image size
291 for(int i=0; i<6; i++) {
292 mReducedExtent[i] = mInitialExtent[i];
297 CropFilterType::Pointer filter = CropFilterType::New();
298 filter->SetInputVVImage(mCurrentImage);
299 filter->SetArgsInfo(mArgsInfo);
301 // Go ! (not threaded)
305 catch(clitk::ExceptionObject & e) {
307 QApplication::restoreOverrideCursor();
308 delete [] mArgsInfo.boundingBox_arg;
311 std::ostringstream croppedImageName;
312 croppedImageName << "Cropped_" << mCurrentSlicerManager->GetSlicer(0)->GetFileName() << ".mhd";
313 // Retrieve result and display it
314 vvImage::Pointer output = filter->GetOutputVVImage();
316 output->GetTransform()[0]->SetMatrix(mConcatenedTransform);
318 AddImage(output,croppedImageName.str());
321 QApplication::restoreOverrideCursor();
322 delete [] mArgsInfo.boundingBox_arg;
325 //------------------------------------------------------------------------------