第四章 transformix

4.1 介紹

到目前為止瘪板,您至少可以通過正確調(diào)用elasticix來進(jìn)行配準(zhǔn)甸陌。 將配準(zhǔn)發(fā)現(xiàn)的轉(zhuǎn)換應(yīng)用于另一個圖像通常也是有用的疗韵。 也許您想將轉(zhuǎn)換應(yīng)用于原始(較大)圖像以獲得分辨率权薯。 或者您可能需要轉(zhuǎn)換以將其應(yīng)用于標(biāo)簽圖像(細(xì)分)嗓化。 為了這些目的,可以使用一個名為transformix的程序搅方。 它與elastix同時開發(fā)比吭。

4.2 如何調(diào)用transformix

像elastix一樣,transformix是一個命令行驅(qū)動的程序姨涡。 您可以通過以下方式獲得有關(guān)如何調(diào)用的基本幫助:transformix --help 這將給出一個強(qiáng)制和可選參數(shù)的列表衩藤。

最基本的命令如下:
transformix -in inputImage.ext -out outputDirectory -tp TransformParameters.txt
此調(diào)用將轉(zhuǎn)換輸入圖像并將其與日志文件transformix.log一起寫入輸出目錄。 您要應(yīng)用的轉(zhuǎn)換在transform參數(shù)文件中定義绣溜。 轉(zhuǎn)換參數(shù)文件可能是以前運(yùn)行elastix的結(jié)果(見第3.3節(jié))慷彤,但也可以由您自己編寫。 第4.3節(jié)說明了轉(zhuǎn)換參數(shù)文件應(yīng)具有的結(jié)構(gòu)和內(nèi)容怖喻。

除了使用transformix進(jìn)行圖像變形外底哗,還可以使用transformix來評估某些點x∈ΩF的變換Tμ(x)。 這意味著輸入點在固定的圖像區(qū)域(C小)中被指定跋选,因為變換方向是從固定到運(yùn)動圖像,如第2.6節(jié)所述哗蜈。 如果要變形一組用戶指定的點前标,則相應(yīng)的調(diào)用是:
transformix -def inputPoints.txt -out outputDirectory -tp TransformParameters.txt
這將創(chuàng)建一個文件outputpoints.txt坠韩,其中包含輸入點x和變換點Tμ(x)(給定為固定圖像的體素索引,另外作為物理坐標(biāo))炼列,位移矢量Tμ(x)-x(物理坐標(biāo)) 只搁,并且如果還指定了-in inputImage.ext,則轉(zhuǎn)換后的輸出將作為輸入圖像的索引俭尖。 inputPoints.txt文件應(yīng)具有以下結(jié)構(gòu):

<index, point>
<number of points>
point1 x point1 y [point1 z]
point2 x point2 y [point2 z]
. . .
```
第一行表示是否將點作為固定圖像的“索引”氢惋,或作為“點”(物理坐標(biāo))給出。 第二行存儲要指定的點數(shù)稽犁。 之后給出點數(shù)據(jù)焰望。
而不是輸入點的自定義.txt格式,transformix還支持.vtk文件:
`transformix -def inputPoints.vtk -out outputDirectory -tp TransformParameters.txt`
輸出然后保存為outputpoints.vtk已亥。 對.vtk文件的支持仍然有限熊赖。 目前,僅支持ASCII文件虑椎,帶有三角形網(wǎng)格震鹉。 任何元數(shù)據(jù)在輸出文件中丟失。

如果您想知道固定圖像的所有體素的變形绣檬,只需使用-def all:
`transformix -def all -out outputDirectory -tp TransformParameters.txt`

變形字段存儲為矢量圖像deformField.mhd足陨。 每個體素在物理坐標(biāo)中包含位移向量Tμ(x) - x嫂粟。 向量的元素存儲為浮點值娇未。

除了計算變形場外,transformix還具有計算變換空間雅可比的能力星虹。 空間Jacobian的決定因素識別局部壓縮或擴(kuò)張的量零抬,并且可能是非常有用的,例如在肺通氣研究中宽涌。 可以在整個圖像上計算空間雅可比的行列式:
`transformix -jac all -out outputDirectory -tp TransformParameters.txt`
也可以計算完整的空間雅可比矩陣:
`transformix -jacmat all -out outputDirectory -tp TransformParameters.txt`
其中每個體素填充有d×d矩陣平夜,其中d是圖像維度,而不是簡單的標(biāo)量值卸亮。

使用命令行選項-threads unsigned int忽妒,用戶可以指定transformix將使用的最大線程數(shù)。

# 4.3 轉(zhuǎn)換參數(shù)文件
配準(zhǔn)的結(jié)果是與固定和運(yùn)動圖像相關(guān)的變換Tμ兼贸。 此變換的參數(shù)存儲在TransformParameters.?.txt-file中段直。 2D剛體變換的結(jié)構(gòu)示例在附錄B中給出。文本文件包含將輸入圖像(運(yùn)動圖像)重新取樣到文件中指定的區(qū)域(默認(rèn)為固定圖像區(qū)域)所需的所有信息溶诞。

轉(zhuǎn)換參數(shù)文件可以手動編輯或創(chuàng)建鸯檬,方便用戶使用。 通過迭代提供另一個轉(zhuǎn)換參數(shù)文件與InitialTransformParametersFileName標(biāo)簽組成多個轉(zhuǎn)換螺垢。 最后的轉(zhuǎn)換將是初始轉(zhuǎn)換參數(shù)文件名設(shè)置為“NoInitialTransform”的轉(zhuǎn)換喧务。

轉(zhuǎn)換參數(shù)文件中的一個重要參數(shù)是FinalBSplineInterpolationOrder赖歌。 通常設(shè)置為3,因為在配準(zhǔn)后產(chǎn)生最佳質(zhì)量的結(jié)果圖像功茴,請參見第5.3.4節(jié)庐冯。 但是,如果您使用transformix來變形運(yùn)動圖像的分割(因此坎穿,二進(jìn)制圖像)肄扎,則需要手動將FinalBSplineInterpolationOrder更改為0。這將確保變形分割仍是二進(jìn)制標(biāo)簽圖像赁酝。 如果使用三階插值犯祠,則變形的分割圖像將包含垃圾。 這與高階B樣條插值的“過沖屬性”有關(guān)酌呆。

# 4.4 一些細(xì)節(jié)
### 4.4.1 運(yùn)行
transformix的運(yùn)行時間由以下部分組成:
  1. 計算輸入圖像的B樣條分解(如果您選擇了FinalBSplineInterpolator);
  2. 計算每個體素的變換;
  3. 插入每個體素的輸入圖像衡载。

我們從未進(jìn)行過測試來測量每個步驟的計算復(fù)雜度,但是我們認(rèn)為第1步是最耗時的任務(wù)隙袁。 顯然可以通過使用最近鄰或線性內(nèi)插器來避免這一步驟痰娱。 步驟2取決于變換的選擇,其中線性變換(例如剛性和仿射變換)比非線性變換(例如B樣條變換)快得多菩收。 步驟3取決于具體的內(nèi)插器梨睁。 按照復(fù)雜度順序,最近鄰娜饵,線性坡贺,一階B樣條,二階B樣條等箱舞。

### 4.4.2 內(nèi)存消耗
有關(guān)內(nèi)存消耗的更多信息遍坟,請參見第5.5.3節(jié),并且:`http://elastix.bigr.nl/wiki/index.php/Memory_consumption_transformix`










最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末晴股,一起剝皮案震驚了整個濱河市愿伴,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌电湘,老刑警劉巖隔节,帶你破解...
    沈念sama閱讀 222,590評論 6 517
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異寂呛,居然都是意外死亡怎诫,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 95,157評論 3 399
  • 文/潘曉璐 我一進(jìn)店門昧谊,熙熙樓的掌柜王于貴愁眉苦臉地迎上來刽虹,“玉大人,你說我怎么就攤上這事呢诬∮空埽” “怎么了胖缤?”我有些...
    開封第一講書人閱讀 169,301評論 0 362
  • 文/不壞的土叔 我叫張陵,是天一觀的道長阀圾。 經(jīng)常有香客問我哪廓,道長,這世上最難降的妖魔是什么初烘? 我笑而不...
    開封第一講書人閱讀 60,078評論 1 300
  • 正文 為了忘掉前任涡真,我火速辦了婚禮,結(jié)果婚禮上肾筐,老公的妹妹穿的比我還像新娘哆料。我一直安慰自己,他們只是感情好吗铐,可當(dāng)我...
    茶點故事閱讀 69,082評論 6 398
  • 文/花漫 我一把揭開白布东亦。 她就那樣靜靜地躺著,像睡著了一般唬渗。 火紅的嫁衣襯著肌膚如雪典阵。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 52,682評論 1 312
  • 那天镊逝,我揣著相機(jī)與錄音壮啊,去河邊找鬼。 笑死撑蒜,一個胖子當(dāng)著我的面吹牛歹啼,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播减江,決...
    沈念sama閱讀 41,155評論 3 422
  • 文/蒼蘭香墨 我猛地睜開眼染突,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了辈灼?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 40,098評論 0 277
  • 序言:老撾萬榮一對情侶失蹤也榄,失蹤者是張志新(化名)和其女友劉穎巡莹,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體甜紫,經(jīng)...
    沈念sama閱讀 46,638評論 1 319
  • 正文 獨居荒郊野嶺守林人離奇死亡降宅,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 38,701評論 3 342
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了囚霸。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片腰根。...
    茶點故事閱讀 40,852評論 1 353
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖拓型,靈堂內(nèi)的尸體忽然破棺而出额嘿,到底是詐尸還是另有隱情瘸恼,我是刑警寧澤,帶...
    沈念sama閱讀 36,520評論 5 351
  • 正文 年R本政府宣布册养,位于F島的核電站东帅,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏球拦。R本人自食惡果不足惜靠闭,卻給世界環(huán)境...
    茶點故事閱讀 42,181評論 3 335
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望坎炼。 院中可真熱鬧愧膀,春花似錦、人聲如沸谣光。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,674評論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽抢肛。三九已至狼钮,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間捡絮,已是汗流浹背熬芜。 一陣腳步聲響...
    開封第一講書人閱讀 33,788評論 1 274
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點兒被人妖公主榨干…… 1. 我叫王不留福稳,地道東北人涎拉。 一個月前我還...
    沈念sama閱讀 49,279評論 3 379
  • 正文 我出身青樓,卻偏偏與公主長得像的圆,于是被迫代替她去往敵國和親鼓拧。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 45,851評論 2 361

推薦閱讀更多精彩內(nèi)容

  • 本章介紹了基于elastix的基本配準(zhǔn)概念越妈。 更高級的配準(zhǔn)主題將在第6章中討論季俩。圖像配準(zhǔn)是醫(yī)學(xué)影像領(lǐng)域的重要工具。...
    peterpan_hai閱讀 9,846評論 1 10
  • 6.1 度量 6.1.1 使用多個指標(biāo)和/或圖像進(jìn)行圖像配準(zhǔn) 到目前為止梅掠,我們將圖像配準(zhǔn)視為找到一個固定圖像和一個...
    peterpan_hai閱讀 1,327評論 1 1
  • 《ilua》速成開發(fā)手冊3.0 官方用戶交流:iApp開發(fā)交流(1) 239547050iApp開發(fā)交流(2) 1...
    葉染柒丶閱讀 10,808評論 0 11
  • 《ijs》速成開發(fā)手冊3.0 官方用戶交流:iApp開發(fā)交流(1) 239547050iApp開發(fā)交流(2) 10...
    葉染柒丶閱讀 5,194評論 0 7
  • 點滴銀線敲打我的窗 像 絲縷回憶輕叩我的心 碧波之蓮?fù)ねぴ阢y錠河畔 像 記憶之花盛開在遺忘彼岸 我的心波瀾不驚 未...
    不上相姑娘閱讀 266評論 0 4