TransformixFilter does not process non-float pixel type images

[thewtex]

When itk::TransformixFilter is templated with a TMovingImage that does not match the internal ElastixTemplate TMovingImage, this dynamic_cast will fail:

elastix/Core/Kernel/elxElastixTemplate.hxx

Lines 100 to 101 in a187f9a

	return dynamic_cast< MovingImageType * >(
	this->GetMovingImageContainer()->ElementAt( idx ).GetPointer() );

and no moving image will be detected here:

elastix/Core/Kernel/elxElastixTemplate.hxx

Lines 387 to 388 in a187f9a

	if( ( this->GetNumberOfMovingImageFileNames() > 0 )
	|| ( this->GetMovingImage() != nullptr ) )

and no transformation is applied here:

elastix/Core/Kernel/elxTransformixMain.cxx

Line 145 in a187f9a

errorCode = this->GetElastixBase()->ApplyTransform();

[thewtex]

A simple check for this is changing the PixelType from float to double in the test added in: #323

[thewtex]

@N-Dekker @mstaring @stefanklein is it possible to take a look at this? I think a cast of the type is intended and required somewhere, which I do not know where that should be.

[dzenanz]

I was unaware of this as I was trying to use transformix to deform an atlas (uchar image).

[N-Dekker]

Current status:

1. The pixel type of the input image should correspond with the template argument of TransformixFilter<TMovingImage>:
   

   elastix/Core/Main/itkTransformixFilter.h

   Lines 58 to 59 in d8e9e1b

   	template <typename TMovingImage>
   	class ITK_TEMPLATE_EXPORT TransformixFilter : public ImageSource<TMovingImage>

   When using Python + ITKElastix, the image type may be specified in square brackets:

ImageType = itk.Image[itk.F, 2]
transformix_filter = itk.TransformixFilter[ImageType].New()

2. The pixel type of the images may also need to be specified as parameter, to be included with the parameter map, to be added to the parameter object. In Python, something like the following (untested code):

parameter_map = 
{
  "FixedInternalImagePixelType": ("float"),
  "MovingInternalImagePixelType": ("float")
   ... [other parameters here] ...
}

parameter_object = itk.ParameterObject.New()
parameter_object.AddParameterMap(parameter_map)
transformix_filter.SetMovingImage(moving_image)
transformix_filter.SetTransformParameterObject(parameter_object)

3. The pixel type must be included during CMake configuration of elastix (ELASTIX_IMAGE_2D_PIXEL_TYPES, ELASTIX_IMAGE_3D_PIXEL_TYPES, etc)

That's it for now! Even though it may be somewhat inconvenient, I think it works well.