十二、ITK例程-医疗图像配准的HelloWorld程序

不问归期 提交于 2020-08-20 08:23:21

一、说明

  医疗图像配准是ITK的一个重要内容,而我们今天想要说的一个程序则相当于是其中的HelloWorld程序。

  程序源码位置:

InsightToolkit-5.0.1\Examples\RegistrationITKv4\\ImageRegistration1.cxx

 二、ITK的配准框架

  

三、程序的主要思路

  程序的主要思路可以用如下的框图来表示:

  

  其中我们需要输入输出的文件如下:

  

 

   但是在这里,为了便于程序阅读,这个时候我把所有的参数全部都写在了程序里面。

四、工程构建

 1-构建工程

  按照之前的方法,利用已经有的源文件构建一个工程,并且添加属性表和动态链接库文件。

 2-程序修改

  源码程序这个时候不是能用的,需要稍作修改

  2.1将每一个argv参数的位置直接修改为文件名

  参数含义如下:

参数 内容
argv[1] 参考图片
argv[2] 待配准图片
argv[3] 输出待配准图片(以参考图片格式输出)
argv[4] 两个输入图片配准之后的差值图片
argv[5] 两个输入图片未配准的时候的差值图片

   实际程序里面修改:

//参考图片位置
fixedImageReader->SetFileName("D:\\FilesSetup\\ITK\\InsightToolkit-5.0.1\\Examples\\Data\\BrainProtonDensitySliceBorder20.png"); //待配准图片位置
movingImageReader
->SetFileName("D:\\FilesSetup\\ITK\\InsightToolkit-5.0.1\\Examples\\Data\\BrainProtonDensitySliceShifted13x17y.png");
//文件写入位置,这三张图片全部都写入到本工程目录下
writer->SetFileName("D:\\Files\\ITKFiles\\ITK_6_Registration\\Data\OutputData\\Moving13x17yInputType.png");
writer2->SetFileName("D:\Files\ITKFiles\ITK_6_Registration\Data\OutputData\Moving13x17yInputType.png");
writer2->SetFileName("D:\\Files\\ITKFiles\\ITK_6_Registration\\Data\\OutputData\\DifferenceBeforeRegistration.png");




  2.2添加头文件

  因为要阅读png文件,所以需要添加png文件的IOFactory

#include "itkPNGImageIOFactory.h"

  初始化IOFactory:

itk::PNGImageIOFactory::RegisterOneFactory();

  2.3添加命令行参数

  因为程序里面有对输入参数的判断,当我们修改以后,还是要满足输入参数的个数(虽然我们已经不用这个参数了,但是这个判断条件是程序运行的触发条件)

  打开调试-属性:

  

 

   在调试-命令参数里面输入六个随意数字,以空格间隔:

 

 五、运行结果

  得到的几张图片与参考书吻合:

  

 

      控制台输出如下:

  

 

     可以看到,一共经过了36次迭代,最后计算出来的测度函数值为0.0007,然后得到的偏移量是(13,17),和图片的真实偏移量是符合的

六、参考

  InsightSoftwareGuide-Book2-5.0.1 Page189(书籍页码)

七、问题

  今天出了一个糗事,我在构建工程的时候,由于偷懒,直接复制的另外一个工程,然后程序都已经修改好了,不论怎么运行,都没有输出。

  还以为是函数复制出了问题,但是这个时候我才发现,我根本没有把这个main.cpp文件添加到工程里面。

  害的我特地对比了一下是不是函数复制的时候出了问题。

  ......

  不注重细节害死人啊。

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!