1. 程式人生 > >核磁資料處理: FSL系列之flirt命令

核磁資料處理: FSL系列之flirt命令

flirt 是FSL的線性修正工具. 官方提供了非常豐富的呼叫方式:
- 命令Flirt可以啟動FLIRT的GUI + InvertXFM, ApplyXFM, ConcatXFM
- 命令flirt 是命令列工具
- misc flirt utilities – flirt附屬的命令列工具.

flirt是進行仿射變換的主要程式.引數包括了 -in 指定的輸入影象, -ref 指定的參考模板, -mat 指定了仿射變換(將輸入影象配準到ref模板上用的)4*4的仿射矩陣. -out 指定的是完成仿射變換後的輸出, 除此之外, flirt的參考模板也可以是一個某個NifTI檔案,要求NifTI檔案和輸入檔案的體素大小\影象尺寸相同.

要點:
1) 目標函式/代價函式加權(cost function weighting).
引數-refweight, -inweight可以用來指定目標函式或者代價函式的加權,使得代價函式在每個體素上有不同的加權值,當weight=0,就可以去除掉一些不感興趣的區域, 或者對感興趣的區域設定weight取值比較大. 注意weight和mask是不一樣的, 是因為mask可能會引入邊界偽跡,但是weight不會.

2) 自由度 (degrees of freedom)
引數-dof是指自由度, 可以設定為6或者7,或者9.
dof自由度指標更靈活的設定方式參考$FSLDIR/etc/flirtsch.

3) 插值演算法(interpolation methods)
最近鄰演算法, 一族基於正弦曲線的插值演算法(三種window型別:矩形視窗, 漢寧視窗或者Blackman視窗,視窗寬度可設定),還有樣條插值方法. 插值只有在最後的變換中有用(applyxfm),而不是在registration計算中.

4) 目標函式/代價函式(cost function)
within-modality函式的構造採用了附加正則項的最小二乘法, between-modality函式預設為Correction Ratio, 另外還有互資訊(Mutual Information)和正則互資訊(Normalised Mutual Information). 除此之外,還有BBR代價函式, 採用了模板影象的segmentation來定義邊界, 代價函式用輸入影象和模板邊界兩側的密度區別intensity difference作為優化目標.

下面給出了一條簡單的呼叫flirt的命令列語句:
flirt -in invol -ref refvol -out outvol -omat invol2refvol.mat -dof 6
這裡invol是輸入的T1影象, refvol應該是registration的模板,例如MNI152, invol2refvol.mat 是registration計算中用到的仿射變換矩陣, dof是指自由度(degrees of freedom), 可以設定為6或者7,或者9.
flirt -in newvol -ref refvol -out outvol -init invol2refvol.mat -applyxfm
還可以將之前的仿射變換矩陣作為輸入設定, 提供給flirt 命令. -init 選項就可以.

下面這條命令得到的結果在fsleyes中開啟, 可以看到
flirt -in brain_1203.nii.gz -ref MNI152_T1_1mm_brain.nii.gz -out brain_registrated.nii.gz -omat invol2refvol.mat

這裡寫圖片描述
這裡寫圖片描述

brain_1203 綠色
brain_registrated 紅色
MNI152 灰色

相關推薦

核磁資料處理: FSL系列flirt命令

flirt 是FSL的線性修正工具. 官方提供了非常豐富的呼叫方式: - 命令Flirt可以啟動FLIRT的GUI + InvertXFM, ApplyXFM, ConcatXFM - 命令flirt 是命令列工具 - misc flirt util

核磁資料處理: FSL系列segmentation工具FAST

fsl的segmentation工具: FAST (FMRIB’s Automated Segmentation Tool)將大腦組織的3D影響分解為不同的組織,比如灰質\白質\腦脊液等等, 同時還可以進行spatial intensity variations

核磁資料處理配準標準模板

說到registration,就必須要首先介紹一下模板.由於人腦存在個體差異, 掃描時候的影象在空間中的座標也不同, 研究時必須要首先消除個體差異,統一座標, 也就是需要將被試人員的大腦全部”訂正/配準”到標準模板上, 這樣後續才能進行統計分析. 目前國際上最

docker系列基礎命令-1

列表 http port ML 配置 狀態 如何 art commit 1.docker基礎命令 docker images 顯示鏡像列表 docker ps 顯示容器列表 docker run IMAGE_ID 指定鏡像, 運行一個容器 docker start/stop

docker系列基礎命令-2

doc tar 維護 標簽 mit font ESS star 分享圖片 一.查看本地鏡像 docker images 二.需要基礎的鏡像兩種方式 1.docker pull centos 可以直接拉起鏡像 2.直接用xshell導入就行,docker  load -i

自動化運維系列Ansible命令應用基礎(模塊的應用)【持續更新中···】

cto fff 安裝 ecb code 文件的 核心 取值 back 自動化運維系列之Ansible命令應用基礎(模塊的應用) 模塊簡介 Ansible可以使用命令行方式進行自動化管理,基本語法如下: ansible <host-pattern> [-m mo

Linux系列SAR命令使用詳解

1. CPU利用率 sar -p (檢視全天) sar -u 1 10 (1:每隔一秒,10:寫入10次) 1.1. CPU輸出項說明 輸出項 詳細說明 CPU all 表示統計資訊為

資料Hadoop學習系列Hadoop、Spark學習路線

1 Java基礎: 視訊方面:推薦畢老師《畢向東JAVA基礎視訊教程》。 學習hadoop不需要過度的深入,java學習到javase,在多執行緒和並行化多多理解實踐即可。 書籍方面:推薦李興華的《java開發實戰經典》 2 Linux基礎: 視訊方面:(1)馬哥的高薪Linux

資料入門教程系列Hive篇彙總

大資料入門教程系列之Hive篇彙總 1、Hive概述、體系架構、分割槽、桶與傾斜 2、大資料入門教程系列之Hive環境搭建--Linux(Centons7)下配置安裝Hive 3、Hive DDL DML及SQL操作 4、Hadoop2.7.4 編譯 支援5種壓縮(Hive壓縮)

資料入門教程系列Hive內建函式及自定義函式

本篇文章主要介紹Hive內建函式以及自定義UDF函式和UDFT函式,自定義UDF函式通過一個國際轉換中文的例子說明。 操作步驟: ①、準備資料和環境 ②、演示Hive內建函式 ③、自定義UDF函式編寫、演示   詳細步驟: 一、準備資料和

資料入門教程系列Hive的Java API 操作

Java 想要訪問Hive,需要通過beeline的方式連線Hive,hiveserver2提供了一個新的命令列工具beeline,hiveserver2 對 之前的hive做了升級,功能更加強大,它增加了許可權控制,要使用beeline需要先啟動hiverserver2,再使用beeline連線

資料入門教程系列Hbase安裝(偽分散式)

本篇文章主要介紹在windows下使用vmware虛擬機器中的Linux(Centons7)下配置安裝Hbase。 簡要步驟: ①、下載Hbase解壓包 ②、安裝zookeeper ③、解壓Hbase,配置環境環境 ④、修改Hbase指令碼檔案hbase-env.sh ⑤、修

資料調錯系列hadoop在開發工具控制檯上打印不出日誌的解決方法

(1)在windows環境上配置HADOOP_HOME環境變數 (2)在eclipse上執行程式 (3)注意:如果eclipse打印不出日誌,在控制檯上只顯示 1.log4j:WARN No appenders could be found 

資料調錯系列:自己總結的myeclipse連線hadoop會出現的問題

在我們學習或者工作中開始hadoop程式的時候,往往會遇到一個問題,我們寫好的程式需要打成包放在叢集中執行,這無形中在浪費我們的時間,因為程式可以需要不斷的除錯,然後把最終程式放在叢集中即可。為了解決這個問題,現在我們配置遠端連線hadoop,遠端除錯的方法。 一段程式如下:獲取更多大資料視訊資料請加QQ群

大數據系列HBase命令行操作

bin sha 面向列 運行 com regions linu 說明 delete 1、課程簡介 本文章先會介紹HBase命令行,接著會介紹java代碼對hbase中的表進行增刪改查。 本文章中所有命令均在CentOS-6.4-x86_64,hadoop-2.5.2,jd

資料入門教程系列Hadoop環境搭建、軟體準備

本篇文章主要介紹大資料入門教程環境搭建所需要的軟體、環境配置等,為後面的學習做準備。 完成目標: 1、安裝VMware及新建虛擬機器 安裝VMware及新建虛擬機器 2、VMware中安裝centons7 虛擬機器VMware中安裝linux系統CentOS

資料結構學習系列二叉搜尋樹詳解!

寫在前面 近期準備補一下資料結構,尤其是關於Tree系列的,其中,二叉樹(Binary Tree)可以算是最簡單的之一,所以打算從之入手,將各種Tree的結構和操作都進一步瞭解一遍,以來充實自己的閒餘時間! 本文主要圍繞二叉樹中最簡單的實現:二叉搜尋樹。 介紹 二叉搜尋樹(Binary Search

Linux系列 tar 命令詳解

tar 是 unix/linux下的打包器 【解壓】 輸入命令: # tar  -zxvf  filename.tar.gz 引數解釋: z :表示 tar 包是被 gzip 壓縮過的 (字尾是.tar.gz),所以解壓時需要用 gu

資料入門基礎系列初步認識大資料生態系統圈(博主推薦)

  不多說,直接上乾貨!   之前在微信公眾平臺裡寫過 大資料入門基礎系列之初步認識hadoop生態系統圈 http://mp.weixin.qq.com/s/KE09U5AbFnEdwht44FGrOA 大資料入門基礎系列之初步認識大資料生態系統圈 1.概述

資料入門基礎系列Hadoop1.X、Hadoop2.X和Hadoop3.X的多維度區別詳解(博主推薦)

  不多說,直接上乾貨! 在前面的博文裡,我已經介紹了   見下面我寫的微信公眾號博文 歡迎大家,加入我的微信公眾號:大資料躺過的坑     免費給分享 同時,大家可以關注我的個人部