ITK配准--Transform 使用的注意事项

1,最简单的方法,直接使用registration->GetTransform()

如果配准之前，没有对Transform进行初始化的，则配准之后的registration->GetTransform()可以直接使用，例子如下：

    using RegistrationType = itk::ImageRegistrationMethodv4<FixedImageType,
                                                         MovingImageType,
                                                         TransformType>;
	registration->SetMetric(metric);
	registration->SetOptimizer(optimizer);
    registration->SetFixedImage(fixed_image);
    registration->SetMovingImage(moving_image);
    registration->Update();
    TransformType::ConstPointer transform = registration->GetTransform();

2,调用SetMovingInitialTransform对Moving图像进行初始化情况，需要通过CompositeTransform拼接Transform

由于需要配准的目标在两个图像中的位置相差较远，为了提高配准的匹配度和速度，往往会在配准之前对Moving图像进行一定的平移或者旋转，而这个操作则可以通过SetMovingInitialTransform函数传入配准过程。但是这个InitialTransform是不会在配准之后同步到配准的RegistrationTransform中，需要人为的将二者合并在一起，例子如下：

    using RegistrationType = itk::ImageRegistrationMethodv4<FixedImageType,
                                                         MovingImageType,
                                                         TransformType>;
    // 初始化MovingTransform
    auto movingInitialTransform = TransformType::New();
    TransformType::ParametersType initialParameters(movingInitialTransform
    							->GetNumberOfParameters());
    initialParameters[0] = -5; // offset in mm along x
    initialParameters[1] = -5; // offset in mm along y
    movingInitialTransform->SetParameters(initialPatameters);
    registration->SetMovingInitialTransform(movingInitialTransform);
    
	registration->SetMetric(metric);
	registration->SetOptimizer(optimizer);
    registration->SetFixedImage(fixed_image);
    registration->SetMovingImage(moving_image);
    registration->Update();
    
    // 与上个例子的关键的差异处
    using CompositeTransformType = itk::CompositeTransform<double, Dimension>;
    CompositeTransformType::Pointer outputCompositeTransform =
    CompositeTransformType::New();
    outputCompositeTransform->AddTransform( movingInitialTransform );
    outputCompositeTransform->AddTransform( registration->
                                            GetModifiableTransform() );

调用SetInitialTransform,直接使用registration->GetTransform()

配准框架的SetInitialTransform和例子2的SetMovingInitialTransform看似结构一样，但是实现的过程却是完全不一样的。如例子2所示，采用SetMovingInitialTransform对Transform进行初始化的时候，配准直接输出的Transform是不能直接使用的，需要结合MovingInitialTransform；而函数SetInitialTransform是直接修改了配准输出的Transform，所以配准之后，Transform已经结合了InitialTransform，例子如下：

    using RegistrationType = itk::ImageRegistrationMethodv4<FixedImageType,
                                                         MovingImageType,
                                                         TransformType>;
 	// 与例子1的关键差异
 	registration->SetInitialTransform(initialTransform);
 	
	registration->SetMetric(metric);
	registration->SetOptimizer(optimizer);
    registration->SetFixedImage(fixed_image);
    registration->SetMovingImage(moving_image);
    registration->Update();
    TransformType::ConstPointer transform = registration->GetTransform();