第三章 elastix

3.1 介绍

what where
Website http://elastix.isi.uu.nl
SVN repository https://svn.bigr.nl/elastix/trunkpublic
Dashboard http://my.cdash.org/index.php?project=elastix
WIKI http://elastix.isi.uu.nl/wiki.php
FAQ http://elastix.isi.uu.nl/FAQ.php
Mailing list subscription http://lists.bigr.nl/mailman/listinfo/elastix
Mailing list elastix@bigr.nl

elastix的发展始于2003年半至晚期,旨在促进我们的配准研究。在一些初始版本之后,我们决定将elastix的单独组件放在单独的库中。 这是2004年11月的主要版本3.0。elastix 3.0也是在elastix网站上公开提供的第一个版本,大约在同一时间。 今天(2015年9月4日)的持续发展带给我们4.8版本。

what where
Website http://elastix.isi.uu.nl
SVN repository https://svn.bigr.nl/elastix/trunkpublic
Dashboard http://my.cdash.org/index.php?project=elastix
WIKI http://elastix.isi.uu.nl/wiki.php
FAQ http://elastix.isi.uu.nl/FAQ.php
Mailing list subscription http://lists.bigr.nl/mailman/listinfo/elastix
Mailing list elastix@bigr.nl

该网站还包含一个doxygen1生成的部分,提供源代码的文档。 有关所有可用类的概述,请参见http://elastix.isi.uu.nl/doxygen/classes.html.
对于每个类,给出了该类的描述,以及如何在elastix中使用它的信息。 参见
http://elastix.isi.uu.nl/doxygen/modules.html 了解所有可用组件的概述。

3.1.1 主要特征

elastix 是

  • 开源,可从http://elastix.isi.uu.nl 免费获得;
  • 基于ITK,所以代码库被彻底测试。 对原始ITK代码进行了一些修改/添加,例如使用采样器,使用组合或添加组合多个转换的转换类等。
  • 适合多种图像格式。 使用ITK意味着,所有由ITK支持的图像格式都由elastix支持。 一些常用的(医疗)图像格式有:.mhd(MetaIO),
    .hdr(Analyze),.nii(NIfTI),.gipl,.dcm(DICOM片)。 elastix不直接支持DICOM目录;
  • 多平台(至少Windows,Linux和Mac OS),多编译器(至少Visual C ++ 2008,2010,gcc 4.x,clang 3.3+),并支持32位和64位系统。 基础ITK代码建立在更多平台上,请参见www.itk.org/Wiki/ITK_Prerequisites。 所以,它是高度可移植到用户选择的平台;
  • 高度可配置:所有配准组件都有很多选择。 通过人类可读和可编辑的参数文件,选择适合您需求的配置非常简单;
  • 易于使用的大量数据,因为elastix 可以轻松地在脚本中调用;
  • 快速,由于随机子采样Klein et al.[2007],感谢多线程和代码优化Shamonin el al.[2014];
  • 相对容易扩展,即添加新的组件,因此它非常适合于研究。

3.2 开始

本节介绍如何直接从二进制文件安装elastix ,或者自己编译elastix 。

3.2.1 最简单开始

使用elastix开始的最简单的方法是使用预编译的二进制文件。

  1. 从网站下载压缩档案:http://elastix.isi.uu.nl/download.php
  2. 将存档提取到您选择的文件夹。
  3. 确保您的操作系统可以找到该程序:
    (a) Windows 7:转到控制面板,进入“系统”,进入“高级系统设置”,点击“环境变量”,将文件夹添加到变量“路径”中。
    (b) Linux:将以下行添加到.bashrc文件中:
export PATH=folder/bin:$PATH
export LD_LIBRARY_PATH=folder/lib:$LD_LIBRARY_PATH

或者使用完整路径调用elastix:fullPathToFolder \ elastix。 请注意,在Linux中,您必须设置LD_LIBRARY_PATH。

3.2.2 简单开始

也可以自己编译elastix,因为源代码是免费的。 在本节中,我们假设您使用Windows下的Microsoft Visual C ++ 2008(或更高版本)编译器,以及Linux / MacOS下的GCC编译器。

  1. 下载并安装CMake: www.cmake.org.
  2. 下载并编译ITK版本4.8:www.itk.org。 确保将以下(高级)CMake变量设置为ON: Module_ITKReview,以及可选的ITK_USE_64BITS_IDS and ITK _LEGACY_REMOVE.。 为了更快地构建ITK,您可以关闭 BUILD_EXAMPLES和BUILD_TESTING。
  3. 获取来源。 有三种可能性:
    (a) 从网站下载压缩源。 将存档提取到您选择的<your-elastix-folder>。
    (b) 使用subversion(https://subversion.tigris.org) 查看从subversion版本库的版本:
svn co --username elastixguest --password elastixguest
https://svn.bigr.nl/elastix/tagspublic/elastix_XX_X <your-elastix-folder>

其中XX_X是主版本号(前2位数字)和次要版本号(1位数字)。 所以,例如,对于版本4.8,这将是04_8。
(c) 使用subversion查看最新的开发版本。 注意:这个版本可能不稳定!

svn co --username elastixguest --password elastixguest
https://svn.bigr.nl/elastix/trunkpublic <your-elastix-folder>
  1. 运行CMake for elastix:
    (a) Windows:启动CMake。 找到具有源代码的src文件夹。 设置要创建二进制文件的文件夹。 单击“配置”,然后选择您使用的编译器。 将CMAKE_ INSTALL_PREFIX设置为要安装elastix的目录。 单击“配置”,直到所有缓存值不再为红色,然后单击“生成”。
    (b) Linux:从要创建二进制文件的文件夹运行CMake,使用提取源的文件夹作为命令行参数:ccmake <src-folder>。 将CMAKE_BUILD_TYPE设置为“Release”,将CMAKE_INSTALL_PREFIX设置为要安装elastix的目录。
    CMake将为您的编译器创建一个项目或解决方案或创建文件。
  2. 通过打开项目并在release模式下选择“编译”,或者通过运行make install在Linux上编译elastix。 您的编译器现在将创建elastix二进制文件。 在Windows上,您还可以通过“编译”INSTALL项目来执行安装步骤(将二进制文件复制到CMAKE_INSTALL_PREFIX目录)。
  3. 确保您的操作系统可以找到程序,请参见上文。

对于开发人员:运行CMake时,您可以切换“高级”选项的显示。 在此列表中,您会发现几个选项,如USE_BSplineTransform ON / OFF。 默认情况下只有最常用的组件是ON。 为了减少编译时间,您可以将某些您不打算使用的组件关闭。 但要小心,不要关闭基本组件。 release的二进制文件将编译所有ON的组件。

3.3 如何调用elastix

elastix是一个命令行程序,虽然自版本4.7(2014年2月)起初支持库接口可用,请参见第7.3.2节。 这意味着您必须打开一个命令行界面(一个DOS框,一个shell)并输入一个相应的elastix命令。 这也意味着没有图形用户界面。 使用该程序的帮助可以如下获取:elastix --help
这将给出一个强制和可选参数的列表。 运行配准的最基本的命令如下:
elastix -f fixedImage.ext -m movingImage.ext -out outputDirectory -p parameterFile.txt
'ext'是图像文件的扩展名。上述参数是强制性的。这些是运行elastix最少需要的。参数文件是一个重要的文件:它在正常文本中包含执行什么样的配准(即什么是度量,优化程序等)以及定义配准的参数。它为该过程提供了大量的灵活性和控制。有关参数文件的更多信息,请参见第3.4节。elastix的所有输出都写入输出目录,在运行elastix之前需要创建它。输出包括一个日志文件(elastix.log),与固定和运动图像相关联的变换Tμ的参数(TransformParameters.?txt),以及可能的结果的注册图像IM(Tμ(x) )(result.?mhd)。日志文件包含在配准期间打印到屏幕的所有消息。此外,parameterFile.txt也被复制到日志文件中,并且包含TransformParameters.?.txt文件的内容。因此,日志文件对故障排除特别有用。

除了强制参数外,还有一些可选参数。 可以通过将-fMask fixedMask.ext和/或-mMask movingMask.ext添加到命令行来提供掩码图像。 可以通过将-t0 TransformParameters.txt添加到命令行来为初始化转换提供有效的转换参数文件。 使用命令行选项-threads unsigned int,用户可以指定elastix将使用的最大线程数。

连续运行多个配准,每个可能是不同的类型,并且以前配准的输出作为下一个配准的输入,可以通过几种方式用elastix来完成。 第一个是在首次配准时运行elastix一次,并使用其输出(在输出目录中可以找到TransformParameter.0.txt)作为输入,使用命令行参数-t0进行新的elastix运行。 所以:

elastix -f ... -m ... -out out1 -p param1.txt
elastix -f ... -m ... -out out2 -p param2.txt -t0 out1/TransformParameters.0.txt
elastix -f ... -m ... -out out3 -p param3.txt -t0 out2/TransformParameters.0.txt

等等。 另一种可能性是将配准与一个elastix组合:
elastix ... -p param1.txt -p param2.txt -p param3.txt
使用(2.18)和(2.19)的一个等式,自动组合来自每个配准的转换。

在elastix网站上,在“关于”部分,您可以找到有关如何使用该程序的示例。 也许现在是尝试这个例子并看到一个配准的时候了。

3.4 参数文件

参数文件是定义配准组件及其参数值的文本文件。 提供参数如下:
(ParameterName value(s))
所以在括号之间提供参数,首先是名称,后跟一个或多个值。 如果值为string类型,那么需要引用这些值:
(ParameterName "value1" ... "valueN")

可以通过用'//'启动行来提供注释。 给定参数的顺序并不重要,但参数只能指定一次。 附录A给出了有效参数文件的一个最小例子。每个类的可用参数列表在http://elastix.isi.uu.nl/doxygen/parameter.html 中给出。 参数文件的示例可以在wiki下找到:http://elastix.bigr.nl/wiki/index.php/Parameter_file_database.

由于几个组件和参数值的选择定义了配准,因此明智地设置它们是非常重要的。 这些选择是什么使配准成功或灾难。 因此,单独的一章专门用于调整配准,参见第5章。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 211,743评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,296评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,285评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,485评论 1 283
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,581评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,821评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,960评论 3 408
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,719评论 0 266
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,186评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,516评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,650评论 1 340
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,329评论 4 330
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,936评论 3 313
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,757评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,991评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,370评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,527评论 2 349

推荐阅读更多精彩内容