X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=vv%2FvvMidPosition.cxx;h=d4192d601ed5877478671d0669e2607c2ca63cc0;hb=a9d632132fa5c0da47d7b3079ccca9d673381a1e;hp=3e4ca3d30f906046794bc67463eb5ad5a1f6db1d;hpb=c1370c27d5fc65b499d6be1c5dc49a96da0bb8cc;p=clitk.git
diff --git a/vv/vvMidPosition.cxx b/vv/vvMidPosition.cxx
index 3e4ca3d..d4192d6 100644
--- a/vv/vvMidPosition.cxx
+++ b/vv/vvMidPosition.cxx
@@ -1,26 +1,20 @@
/*=========================================================================
+ Program: vv http://www.creatis.insa-lyon.fr/rio/vv
- Program: vv
- Language: C++
- Author : Joel Schaerer (joel.schaerer@insa-lyon.fr)
+ Authors belong to:
+ - University of LYON http://www.universite-lyon.fr/
+ - Léon Bérard cancer center http://www.centreleonberard.fr
+ - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
-Copyright (C) 2008
-Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
-CREATIS-LRMN http://www.creatis.insa-lyon.fr
+ This software is distributed WITHOUT ANY WARRANTY; without even
+ the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+ PURPOSE. See the copyright notices for more information.
-This program is free software: you can redistribute it and/or modify
-it under the terms of the GNU General Public License as published by
-the Free Software Foundation, version 3 of the License.
+ It is distributed under dual licence
-This program is distributed in the hope that it will be useful,
-but WITHOUT ANY WARRANTY; without even the implied warranty of
-MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-GNU General Public License for more details.
-
-You should have received a copy of the GNU General Public License
-along with this program. If not, see .
-
-=========================================================================*/
+ - BSD See included LICENSE.txt file
+ - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
+===========================================================================**/
#include
#include
@@ -33,22 +27,21 @@ along with this program. If not, see .
#include "clitkInvertVFFilter.h"
vvMidPosition::vvMidPosition() :
- slicer_manager(0), error(false),
- reference_image_index(0),
- p_bar("Computing mid-position...","Cancel", 0,6),
- progress(0)
+ slicer_manager(0), error(false),
+ reference_image_index(0),
+ p_bar("Computing mid-position...","Cancel", 0,6),
+ progress(0)
{
}
void vvMidPosition::Update()
{
- this->start();
- while (this->isRunning())
- {
- this->wait(50);
- this->update_progress();
- qApp->processEvents();
- }
+ this->start();
+ while (this->isRunning()) {
+ this->wait(50);
+ this->update_progress();
+ qApp->processEvents();
+ }
}
//Common typedefs
@@ -66,141 +59,135 @@ template vvImage::Pointer WarpRefImage(OutputVFType::Point
void vvMidPosition::run()
{
- error=true;
- if(slicer_manager->GetImage()->GetNumberOfDimensions() != 4)
- error_message="Computation of midposition is only supported for 4D image sequences.";
- else if (!slicer_manager->GetVF())
- error_message="A VF is required for midposition computation";
- else if (slicer_manager->GetVF()->GetScalarTypeAsString() != "float")
- error_message="Please use a vector field of type float.";
- else
- {
- VFType::ConstPointer vf = vvImageToITK(slicer_manager->GetVF());
- OutputVFType::Pointer avg_vf=AverageField(vf,this->progress);
- clitk::InvertVFFilter::Pointer inv_filter=
- clitk::InvertVFFilter::New();
- inv_filter->SetInput(avg_vf);
- inv_filter->Update();
- progress++;
- if (slicer_manager->GetImage()->GetScalarTypeAsString() == "short")
- this->output=WarpRefImage(inv_filter->GetOutput(),slicer_manager->GetImage(),reference_image_index);
- else
- {
- error_message="Unsupported image pixel type.";
- return;
- }
- progress++;
- error=false;
+ error=true;
+ if(slicer_manager->GetImage()->GetNumberOfDimensions() != 4)
+ error_message="Computation of midposition is only supported for 4D image sequences.";
+ else if (!slicer_manager->GetVF())
+ error_message="A VF is required for midposition computation";
+ else if (slicer_manager->GetVF()->GetScalarTypeAsITKString() != "float")
+ error_message="Please use a vector field of type float.";
+ else {
+ VFType::ConstPointer vf = vvImageToITK(slicer_manager->GetVF());
+ OutputVFType::Pointer avg_vf=AverageField(vf,this->progress);
+ clitk::InvertVFFilter::Pointer inv_filter=
+ clitk::InvertVFFilter::New();
+ inv_filter->SetInput(avg_vf);
+ inv_filter->Update();
+ progress++;
+ if (slicer_manager->GetImage()->GetScalarTypeAsITKString() == "short")
+ this->output=WarpRefImage(inv_filter->GetOutput(),slicer_manager->GetImage(),reference_image_index);
+ else {
+ error_message="Unsupported image pixel type.";
+ return;
}
+ progress++;
+ error=false;
+ }
}
template
vvImage::Pointer WarpRefImage(OutputVFType::Pointer vf,vvImage::Pointer image,int reference_image_index)
{
- typedef itk::Image ImageType;
- typedef itk::WarpImageFilter > FilterType;
-
- typename ImageType::ConstPointer input = vvSingleFrameToITK<3,ImagePixelType>(image,reference_image_index);
-
- //We resample the VF because itk's warp filter doesn't like it when the vf and the image have
- //different spacings
- typename itk::VectorResampleImageFilter::Pointer
- resampler =itk::VectorResampleImageFilter::New();
- resampler->SetInput(vf);
- resampler->SetOutputSpacing(input->GetSpacing());
- resampler->SetOutputOrigin(vf->GetOrigin());
- //Calculate the new size so that it contains the vf
- typename ImageType::SizeType newSize;
- for (unsigned int i=0 ; i <3; i++)
- newSize[i]=(unsigned int) (vf->GetLargestPossibleRegion().GetSize()[i]*vf->GetSpacing()[i]/input->GetSpacing()[i]);
- resampler->SetSize( newSize);
-
- typename FilterType::Pointer warp_filter = FilterType::New();
- warp_filter->SetInput(input);
- warp_filter->SetDeformationField(resampler->GetOutput());
- warp_filter->SetOutputSpacing(input->GetSpacing());
- warp_filter->SetOutputOrigin(input->GetOrigin());
- warp_filter->SetOutputSize(input->GetLargestPossibleRegion().GetSize());
- warp_filter->SetEdgePaddingValue(-1000);
- warp_filter->Update();
- return vvImageFromITK<3,ImagePixelType>(warp_filter->GetOutput());
+ typedef itk::Image ImageType;
+ typedef itk::WarpImageFilter > FilterType;
+
+ typename ImageType::ConstPointer input = vvSingleFrameToITK<3,ImagePixelType>(image,reference_image_index);
+
+ //We resample the VF because itk's warp filter doesn't like it when the vf and the image have
+ //different spacings
+ typename itk::VectorResampleImageFilter::Pointer
+ resampler =itk::VectorResampleImageFilter::New();
+ resampler->SetInput(vf);
+ resampler->SetOutputSpacing(input->GetSpacing());
+ resampler->SetOutputOrigin(vf->GetOrigin());
+ //Calculate the new size so that it contains the vf
+ typename ImageType::SizeType newSize;
+ for (unsigned int i=0 ; i <3; i++)
+ newSize[i]=(unsigned int) (vf->GetLargestPossibleRegion().GetSize()[i]*vf->GetSpacing()[i]/input->GetSpacing()[i]);
+ resampler->SetSize( newSize);
+
+ typename FilterType::Pointer warp_filter = FilterType::New();
+ warp_filter->SetInput(input);
+ warp_filter->SetDisplacementField(resampler->GetOutput());
+ warp_filter->SetOutputSpacing(input->GetSpacing());
+ warp_filter->SetOutputOrigin(input->GetOrigin());
+ warp_filter->SetOutputSize(input->GetLargestPossibleRegion().GetSize());
+ warp_filter->SetEdgePaddingValue(-1000);
+ warp_filter->Update();
+ return vvImageFromITK<3,ImagePixelType>(warp_filter->GetOutput());
}
itk::Image,3>::Pointer AverageField(itk::Image,4>::ConstPointer vf, int& progress)
{
- progress++;
-
- VFType::RegionType region4D=vf->GetLargestPossibleRegion();
- VFType::RegionType::SizeType size4D=region4D.GetSize();
- VFType::IndexType index4D=region4D.GetIndex();
- VFType::SpacingType spacing4D=vf->GetSpacing();
- VFType::PointType origin4D=vf->GetOrigin();
-
- OutputVFType::RegionType region;
- OutputVFType::RegionType::SizeType size;
- OutputVFType::IndexType index;
- OutputVFType::SpacingType spacing;
- OutputVFType::PointType origin;
- for (unsigned int i=0; i< 3; i++)
- {
- size[i]=size4D[i];
- index[i]=index4D[i];
- spacing[i]=spacing4D[i];
- origin[i]=origin4D[i];
+ progress++;
+
+ VFType::RegionType region4D=vf->GetLargestPossibleRegion();
+ VFType::RegionType::SizeType size4D=region4D.GetSize();
+ VFType::IndexType index4D=region4D.GetIndex();
+ VFType::SpacingType spacing4D=vf->GetSpacing();
+ VFType::PointType origin4D=vf->GetOrigin();
+
+ OutputVFType::RegionType region;
+ OutputVFType::RegionType::SizeType size;
+ OutputVFType::IndexType index;
+ OutputVFType::SpacingType spacing;
+ OutputVFType::PointType origin;
+ for (unsigned int i=0; i< 3; i++) {
+ size[i]=size4D[i];
+ index[i]=index4D[i];
+ spacing[i]=spacing4D[i];
+ origin[i]=origin4D[i];
+ }
+ region.SetSize(size);
+ region.SetIndex(index);
+ OutputVFType::Pointer output= OutputVFType::New();
+ output->SetRegions(region);
+ output->SetSpacing(spacing);
+ output->SetOrigin(origin);
+ output->Allocate();
+ progress++;
+
+
+ // Region iterators
+ typedef itk::ImageRegionConstIterator IteratorType;
+ std::vector iterators(size4D[3]);
+ for (unsigned int i=0; i< size4D[3]; i++) {
+ VFType::RegionType regionIt=region4D;
+ VFType::RegionType::SizeType sizeIt=regionIt.GetSize();
+ sizeIt[3]=1;
+ regionIt.SetSize(sizeIt);
+ VFType::IndexType indexIt=regionIt.GetIndex();
+ indexIt[3]=i;
+ regionIt.SetIndex(indexIt);
+ iterators[i]=IteratorType(vf, regionIt);
+ }
+ progress++;
+
+ typedef itk::ImageRegionIterator OutputIteratorType;
+ OutputIteratorType avIt(output, output->GetLargestPossibleRegion());
+
+ // Average
+ VFPixelType vector;
+ VFPixelType zeroVector;//=itk::NumericTraits::Zero;
+ for(unsigned int i=0;i SetRegions(region);
- output->SetSpacing(spacing);
- output->SetOrigin(origin);
- output->Allocate();
- progress++;
-
-
- // Region iterators
- typedef itk::ImageRegionConstIterator IteratorType;
- std::vector iterators(size4D[3]);
- for (unsigned int i=0; i< size4D[3]; i++)
- {
- VFType::RegionType regionIt=region4D;
- VFType::RegionType::SizeType sizeIt=regionIt.GetSize();
- sizeIt[3]=1;
- regionIt.SetSize(sizeIt);
- VFType::IndexType indexIt=regionIt.GetIndex();
- indexIt[3]=i;
- regionIt.SetIndex(indexIt);
- iterators[i]=IteratorType(vf, regionIt);
- }
- progress++;
-
- typedef itk::ImageRegionIterator OutputIteratorType;
- OutputIteratorType avIt(output, output->GetLargestPossibleRegion());
-
- // Average
- VFPixelType vector;
- VFPixelType zeroVector=itk::NumericTraits::Zero;
-
- while (!(iterators[0]).IsAtEnd())
- {
- vector=zeroVector;
- for (unsigned int i=0; i