mitkShapeBasedInterpolationAlgorithm.cpp
Go to the documentation of this file.
00001 /*========================================================================= 00002 00003 Program: Medical Imaging & Interaction Toolkit 00004 Language: C++ 00005 Date: $Date$ 00006 Version: $Revision: 1.12 $ 00007 00008 Copyright (c) German Cancer Research Center, Division of Medical and 00009 Biological Informatics. All rights reserved. 00010 See MITKCopyright.txt or https://www.mitk.org/copyright.html for details. 00011 00012 This software is distributed WITHOUT ANY WARRANTY; without even 00013 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 00014 PURPOSE. See the above copyright notices for more information. 00015 00016 =========================================================================*/ 00017 00018 #include "mitkShapeBasedInterpolationAlgorithm.h" 00019 #include "mitkImageCast.h" 00020 00021 #include "ipSegmentation.h" 00022 00023 mitk::Image::Pointer 00024 mitk::ShapeBasedInterpolationAlgorithm::Interpolate( 00025 Image::ConstPointer lowerSlice, unsigned int lowerSliceIndex, 00026 Image::ConstPointer upperSlice, unsigned int upperSliceIndex, 00027 unsigned int requestedIndex, 00028 unsigned int /*sliceDimension*/, // commented variables are not used 00029 Image::Pointer resultImage, 00030 unsigned int /*timeStep*/, 00031 Image::ConstPointer /*referenceImage*/) 00032 { 00033 // convert these slices to the ipSegmentation data type (into an ITK image) 00034 itk::Image< ipMITKSegmentationTYPE, 2 >::Pointer correctPixelTypeLowerITKSlice; 00035 CastToItkImage( lowerSlice, correctPixelTypeLowerITKSlice ); 00036 assert ( correctPixelTypeLowerITKSlice.IsNotNull() ); 00037 00038 itk::Image< ipMITKSegmentationTYPE, 2 >::Pointer correctPixelTypeUpperITKSlice; 00039 CastToItkImage( upperSlice, correctPixelTypeUpperITKSlice ); 00040 assert ( correctPixelTypeUpperITKSlice.IsNotNull() ); 00041 00042 // correct direction info 00043 itk::Image< ipMITKSegmentationTYPE, 2 >::DirectionType imageDirection; 00044 imageDirection.SetIdentity(); 00045 correctPixelTypeLowerITKSlice->SetDirection(imageDirection); 00046 correctPixelTypeUpperITKSlice->SetDirection(imageDirection); 00047 00048 // back-convert to MITK images to access a mitkIpPicDescriptor 00049 Image::Pointer correctPixelTypeLowerMITKSlice = Image::New(); 00050 CastToMitkImage( correctPixelTypeLowerITKSlice, correctPixelTypeLowerMITKSlice ); 00051 mitkIpPicDescriptor* lowerPICSlice = correctPixelTypeLowerMITKSlice->GetSliceData()->GetPicDescriptor(); 00052 00053 Image::Pointer correctPixelTypeUpperMITKSlice = Image::New(); 00054 CastToMitkImage( correctPixelTypeUpperITKSlice, correctPixelTypeUpperMITKSlice ); 00055 mitkIpPicDescriptor* upperPICSlice = correctPixelTypeUpperMITKSlice->GetSliceData()->GetPicDescriptor(); 00056 00057 // calculate where the current slice is in comparison to the lower and upper neighboring slices 00058 float ratio = (float)(requestedIndex - lowerSliceIndex) / (float)(upperSliceIndex - lowerSliceIndex); 00059 00060 mitkIpPicDescriptor* ipPicResult = ipMITKSegmentationInterpolate( lowerPICSlice, upperPICSlice, ratio ); // magic 00061 if (!ipPicResult) return NULL; 00062 00063 Geometry3D::Pointer originalGeometry = resultImage->GetGeometry(); 00064 resultImage->Initialize( ipPicResult ); 00065 resultImage->SetPicSlice( ipPicResult ); 00066 resultImage->SetGeometry( originalGeometry ); 00067 00068 mitkIpPicFree( ipPicResult ); 00069 00070 return resultImage; 00071 } 00072
