[Insight-users] ITK is N-D

[Luis Ibanez]

Hi Martin,

ITK is designed and implemented to suport N-Dimensional
images and N-Dimensional Meshes.

Some particular filters have been implemented for a
specific dimension and in those cases the dimension makes
part of the filter name.

About your questions:

1) Yes, you can import from a volume (a 3D image) into
    a VTK data set.

2) Yes, segmentation algorithms can be applied in 3D.
    Watersheds, Regions growing, Classifiers... all are
    implemented for N-D images.


Support for 3D images (volumes) has been from the
beggining one of the design requirements of ITK.

We went one step further and make it N-D   :-)



Please let us know if you encounter any problem
applying these methods to 3D images.


Thanks

Luis

===========================================================

martin.garcia@epfl.ch wrote:
> (Sorry for the first message it was truncated!!!!!)
> 
> HI everyone, 
> I'm trying to develop an application in order to produce medical images 
> segemntation and volume rendering. For visualisation I use vtk and I'd like to 
> use itk for segmentation. 
> 
> To do this I've tried to use vtk2itk example , but I only could do it with 
> single 2D images. 
> 
> 2 questions: 
> 
> 1.Can I do import export from a volume file? 
> (I've got a vtkImageData object called m_pImageData with all the volumic data 
> (spacing, voxel size, image size,slice thickness...)that I' like to import into 
> itk and then export back to vtk. 
> 
> sample code that uses itkvtktoitk example 
>          
> typedef itk::Image<float, 3> ImageType; 
> typedef itk::VTKImageImport<ImageType> ImageImportType; 
> typedef itk::VTKImageExport<ImageType> ImageExportType; 
> 
> vtkImageExport* vtkExporter = vtkImageExport::New(); 
> vtkImageCast * imagecast = vtkImageCast::New(); 
> imagecast->SetOutputScalarTypeToFloat(); 
> imagecast->SetInput(m_pImageData->GetOutput()); 
> vtkExporter->SetInput(imagecast->GetOutput()); 
> 
> ImageImportType::Pointer itkImporter = ImageImportType::New();
> ConnectPipelines(vtkExporter, itkImporter);
> //back to vtk
> ImageExportType::Pointer itkExporter = ImageExportType::New();
>   itkExporter->SetInput(itkImporter->GetOutput());
>   
>   // Create the vtkImageImport and connect it to the
>   // itk::VTKImageExport instance.
>   vtkImageImport* vtkImporter = vtkImageImport::New();  
>   ConnectPipelines(itkExporter, vtkImporter);
>   vtkImageCast * imagecastbk = vtkImageCast::New();
>   imagecastbk->SetOutputScalarTypeToUnsignedChar();
>   imagecastbk->SetInput(vtkImporter->GetOutput());
>   m_pImageData->DeepCopy(imagecastbk->GetOutput());
> 
> Execution errors appear when Itry to visualize back the information (concerning 
> the extent update and the whole extent values)
> 
> 
> 2. Can I apply itk segmentation filters directly on 3 dimensional data
>    Image<float,3>(watershed, region growing,...?
> 
> Thanks, and sorry again for the truncated message
> 
> 
> 
> 
> 2.
> 
> 
> 
> 
> 
> -------------------------------------------------
> This mail sent through IMP: http://horde.org/imp/
> _______________________________________________
> Insight-users mailing list
> Insight-users@public.kitware.com
> http://public.kitware.com/mailman/listinfo/insight-users
> 
>