一、说明
医疗图像配准是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文件添加到工程里面。
害的我特地对比了一下是不是函数复制的时候出了问题。
......
不注重细节害死人啊。
来源:oschina
链接:https://my.oschina.net/u/4347242/blog/4437647