第四章 transformix

4.1 介绍

到目前为止,您至少可以通过正确调用elasticix来进行配准。 将配准发现的转换应用于另一个图像通常也是有用的。 也许您想将转换应用于原始(较大)图像以获得分辨率。 或者您可能需要转换以将其应用于标签图像(细分)。 为了这些目的,可以使用一个名为transformix的程序。 它与elastix同时开发。

4.2 如何调用transformix

像elastix一样,transformix是一个命令行驱动的程序。 您可以通过以下方式获得有关如何调用的基本帮助:transformix --help 这将给出一个强制和可选参数的列表。

最基本的命令如下:
transformix -in inputImage.ext -out outputDirectory -tp TransformParameters.txt
此调用将转换输入图像并将其与日志文件transformix.log一起写入输出目录。 您要应用的转换在transform参数文件中定义。 转换参数文件可能是以前运行elastix的结果(见第3.3节),但也可以由您自己编写。 第4.3节说明了转换参数文件应具有的结构和内容。

除了使用transformix进行图像变形外,还可以使用transformix来评估某些点x∈ΩF的变换Tμ(x)。 这意味着输入点在固定的图像区域(!)中被指定,因为变换方向是从固定到运动图像,如第2.6节所述。 如果要变形一组用户指定的点,则相应的调用是:
transformix -def inputPoints.txt -out outputDirectory -tp TransformParameters.txt
这将创建一个文件outputpoints.txt,其中包含输入点x和变换点Tμ(x)(给定为固定图像的体素索引,另外作为物理坐标),位移矢量Tμ(x)-x(物理坐标) ,并且如果还指定了-in inputImage.ext,则转换后的输出将作为输入图像的索引。 inputPoints.txt文件应具有以下结构:

<index, point>
<number of points>
point1 x point1 y [point1 z]
point2 x point2 y [point2 z]
. . .
```
第一行表示是否将点作为固定图像的“索引”,或作为“点”(物理坐标)给出。 第二行存储要指定的点数。 之后给出点数据。
而不是输入点的自定义.txt格式,transformix还支持.vtk文件:
`transformix -def inputPoints.vtk -out outputDirectory -tp TransformParameters.txt`
输出然后保存为outputpoints.vtk。 对.vtk文件的支持仍然有限。 目前,仅支持ASCII文件,带有三角形网格。 任何元数据在输出文件中丢失。

如果您想知道固定图像的所有体素的变形,只需使用-def all:
`transformix -def all -out outputDirectory -tp TransformParameters.txt`

变形字段存储为矢量图像deformField.mhd。 每个体素在物理坐标中包含位移向量Tμ(x) - x。 向量的元素存储为浮点值。

除了计算变形场外,transformix还具有计算变换空间雅可比的能力。 空间Jacobian的决定因素识别局部压缩或扩张的量,并且可能是非常有用的,例如在肺通气研究中。 可以在整个图像上计算空间雅可比的行列式:
`transformix -jac all -out outputDirectory -tp TransformParameters.txt`
也可以计算完整的空间雅可比矩阵:
`transformix -jacmat all -out outputDirectory -tp TransformParameters.txt`
其中每个体素填充有d×d矩阵,其中d是图像维度,而不是简单的标量值。

使用命令行选项-threads unsigned int,用户可以指定transformix将使用的最大线程数。

# 4.3 转换参数文件
配准的结果是与固定和运动图像相关的变换Tμ。 此变换的参数存储在TransformParameters.?.txt-file中。 2D刚体变换的结构示例在附录B中给出。文本文件包含将输入图像(运动图像)重新取样到文件中指定的区域(默认为固定图像区域)所需的所有信息。

转换参数文件可以手动编辑或创建,方便用户使用。 通过迭代提供另一个转换参数文件与InitialTransformParametersFileName标签组成多个转换。 最后的转换将是初始转换参数文件名设置为“NoInitialTransform”的转换。

转换参数文件中的一个重要参数是FinalBSplineInterpolationOrder。 通常设置为3,因为在配准后产生最佳质量的结果图像,请参见第5.3.4节。 但是,如果您使用transformix来变形运动图像的分割(因此,二进制图像),则需要手动将FinalBSplineInterpolationOrder更改为0。这将确保变形分割仍是二进制标签图像。 如果使用三阶插值,则变形的分割图像将包含垃圾。 这与高阶B样条插值的“过冲属性”有关。

# 4.4 一些细节
### 4.4.1 运行
transformix的运行时间由以下部分组成:
  1. 计算输入图像的B样条分解(如果您选择了FinalBSplineInterpolator);
  2. 计算每个体素的变换;
  3. 插入每个体素的输入图像。

我们从未进行过测试来测量每个步骤的计算复杂度,但是我们认为第1步是最耗时的任务。 显然可以通过使用最近邻或线性内插器来避免这一步骤。 步骤2取决于变换的选择,其中线性变换(例如刚性和仿射变换)比非线性变换(例如B样条变换)快得多。 步骤3取决于具体的内插器。 按照复杂度顺序,最近邻,线性,一阶B样条,二阶B样条等。

### 4.4.2 内存消耗
有关内存消耗的更多信息,请参见第5.5.3节,并且:`http://elastix.bigr.nl/wiki/index.php/Memory_consumption_transformix`










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

推荐阅读更多精彩内容

  • 本章介绍了基于elastix的基本配准概念。 更高级的配准主题将在第6章中讨论。图像配准是医学影像领域的重要工具。...
    peterpan_hai阅读 9,134评论 1 10
  • 6.1 度量 6.1.1 使用多个指标和/或图像进行图像配准 到目前为止,我们将图像配准视为找到一个固定图像和一个...
    peterpan_hai阅读 1,239评论 1 1
  • 《ilua》速成开发手册3.0 官方用户交流:iApp开发交流(1) 239547050iApp开发交流(2) 1...
    叶染柒丶阅读 10,356评论 0 11
  • 《ijs》速成开发手册3.0 官方用户交流:iApp开发交流(1) 239547050iApp开发交流(2) 10...
    叶染柒丶阅读 5,016评论 0 7
  • 点滴银线敲打我的窗 像 丝缕回忆轻叩我的心 碧波之莲亭亭在银锭河畔 像 记忆之花盛开在遗忘彼岸 我的心波澜不惊 未...
    不上相姑娘阅读 258评论 0 4