X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvFromITK.h;h=d44cf728e84a91df68605f0028e3318bc76732a1;hb=1fea92ac237c16c1f027ffabdb42066f6f4f6114;hp=c7178fe3e36c417cb85d37182c6edbea2840a1a0;hpb=f181e8fc01dec00bc6b0c82e0e273e1bb0d292c5;p=clitk.git diff --git a/common/vvFromITK.h b/common/vvFromITK.h index c7178fe..d44cf72 100644 --- a/common/vvFromITK.h +++ b/common/vvFromITK.h @@ -32,7 +32,6 @@ template vvImage::Pointer vvImageFromITK(type { assert(Dim < 5 && Dim > 0); // We don't handle anything higher than 4-dimensional (for the moment :-p) vvImage::Pointer vv_image=vvImage::New(); - vv_image->Init(); //Delete any existing images typedef itk::Image< PixelType, Dim > InputImageType; if (Dim == 4 || time_sequence) //The time sequence case: create a series of VTK images @@ -41,6 +40,7 @@ template vvImage::Pointer vvImageFromITK(type typedef itk::ExtractImageFilter FilterType; //extract the 3D slices and put them in a std::vector + input->UpdateOutputInformation(); typename InputImageType::RegionType inputRegion = input->GetLargestPossibleRegion(); typename InputImageType::SizeType inputSize = inputRegion.GetSize(); typename InputImageType::IndexType start = inputRegion.GetIndex(); @@ -57,15 +57,6 @@ template vvImage::Pointer vvImageFromITK(type filter->SetExtractionRegion(extractedRegion); filter->SetInput(input); filter->ReleaseDataFlagOn(); - - try { - filter->Update(); - } - catch ( itk::ExceptionObject & err ) { - std::cerr << "Error while setting vvImage from ITK (Dim==4) [Extract phase]" - << " " << err << std::endl; - return vv_image; - } vv_image->AddItkImage(filter->GetOutput()); } vv_image->SetTimeSpacing(input->GetSpacing()[Dim-1]);