ITK 配准框架示例

图像配准的基本过程如下：
1.指定用于评估配准效果的相似度或误差测度；
2.指定一个变换模型，如刚体变换、仿射变换、弹性变换(elastic)、流体变换或B-样条等；
3.指定插值策略，如最邻近插值(nearst neighbour)、三线性插值(trilinear)、sinc插值等；
4.寻找变换参数以最大化相似性测度。
如下图所示：

 

配准框架的基本流程如下：
1．输入待配准的两幅图像，参考图 Fixed Image,浮动图 Moving Image.
2．对参考图的指定区域进行几何坐标变换(Transform) 得到新的区域 ;
3．通过插值方法(Image Interpolator)得到浮动图在上一步新区域的坐标；
4．相似性测度模块计算参考图和插值图之间的相似度，是一个关于几何变换参数的函数；
5．相似度函数输入到优化模块中进行最优化计算得到最终变换参数，这个过程一般通过
   迭代来实现，即重复2~4步直到取得最大值；
6．整个配准算法模块输出浮动图在最优变换下的插值图像。 
   配准过程是一个优化问题，配准过程每进行一次迭代，得到一测度值，将该测度值与
我们所设定的值进行比较，如果达到预期的效果则停止迭代，得到最终配准结果。当然，
迭代可能无限制进行，所以我们还需要设置一迭代上限。

下面以一个简单例子对 ITK 中的配准框架进行说明：
ITK 的配准框架由如下几个组件组成：变换组件、相似性测度组件、插值组件和优化组件；各个组件通过一个称为“配准方法”的组件连接到一起，形成一个一个管道结构。

//代码：

 

 #if defined(_MSC_VER) #pragma warning ( disable : 4786 ) #endif //本例演示了使用刚体变换进行 2D 图像配准, 使用的变换为 itk::CenteredRigid2DTransform #include "itkImageRegistrationMethod.h" #include "itkMeanSquaresImageToImageMetric.h" #include "itkLinearInterpolateImageFunction.h" #include "itkRegularStepGradientDescentOptimizer.h" #include "itkImage.h" #include "itkCenteredRigid2DTransform.h" #include "itkImageFileReader.h" #include "itkImageFileWriter.h" #include "itkResampleImageFilter.h" #include "itkSubtractImageFilter.h" #include "itkRescaleIntensityImageFilter.h" #include "itkCommand.h" /******************************************************************** *类 CommandIterationUpdate：创建 Command/Observer 设计模式的简单方式 *本例中它的主要功能为监视配准过程的执行，每迭代一次，调用一次Execute() 方法，输出相应参数 *此处涉及三个类: itk::Object, itk::Command，itk::EventObject。 *Object 是大多数类的基类，该类维护一个着一个指针链表,指针指向事件观察者(event observer). *Observer 的角色由 Command 类扮演, Observers 通过一个 Object 注册自己,声明当一特定事件发 *生时. 配准过程通过一 itk::Optimizer 执行一个迭代过程来控制. 大多数 optimizer 在每次迭代 *后都会调用一 itk::IterationEvent. 当一个事件发生时, 该对象遍历 registred observers(Command) *链表,并检查哪一个对该事件感兴趣, 当找到一个对该事件感兴趣的 observer 时,则调用与其关联的 *Execute() 方法. 在这方面, Execute() 应该看做是一个 callback,所以它应该 遵循一般 callback *的原则,如不应该执行计算机量大的任务,应该只完成一些快速且简单的任务. ********************************************************************/ class CommandIterationUpdate : public itk::Command { public: typedef CommandIterationUpdate Self; typedef itk::Command Superclass; typedef itk::SmartPointer<Self> Pointer; //itkNewMacro: 该宏包装了 New() 方法的所有代码 itkNewMacro( Self ); protected: CommandIterationUpdate() {}; public: typedef itk::RegularStepGradientDescentOptimizer OptimizerType; typedef const OptimizerType * OptimizerPointer; void Execute(itk::Object *caller, const itk::EventObject & event) { Execute( (const itk::Object *)caller, event); } //object: 指向激活事件的对象；event: 需要被激活的事件 void Execute(const itk::Object * object, const itk::EventObject & event) { OptimizerPointer optimizer = dynamic_cast< OptimizerPointer >( object ); //CheckEvent() 检测是否是观察的事件 if( ! itk::IterationEvent().CheckEvent( &event ) ) { return; } std::cout << optimizer->GetCurrentIteration() << " "; std::cout << optimizer->GetValue() << " "; std::cout << optimizer->GetCurrentPosition() << std::endl; } }; int main( int argc, char *argv[] ) { if( argc < 4 ) { std::cerr << "Missing Parameters " << std::endl; std::cerr << "Usage: " << argv[0]; std::cer