張量分析在腦醫學圖像的應用范文

本站小編為你精心準備了張量分析在腦醫學圖像的應用參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

【摘要】醫學影像中的CT、MRI圖像診斷是目前醫生確診疾病的重要依據。另一方面，大腦作為人類最復雜也是最重要的器官之一，對于腦部圖像的特征提取和分類具有重要意義。傳統上對圖像特征提取習慣從向量的角度出發，這樣忽略了圖像結構特點。為了解決這個問題，本文結合高維空間數據結構，以數據張量化為重點將傳統的PCA算法升級為Tensor-PCA，并選擇當下最流行的SimulatedBrainphantomDatabase數據集作為本文的仿真對象，經過實驗證明，數據張量化的方法在提取圖像的特征上具有良好的適用性。

【關鍵詞】張量模型；腦部醫學圖像；主成分分析；特征

提取隨著腦成像技術和計算機的快速發展，使得人們能夠更加深入的探索大腦奧秘。MRI成像技術已成為當代最佳的輔助工具之一，幫助了醫生對病情診斷的同時也提供給學者們作為最佳的實驗選擇對象，從而獲得了相關領域的關注。但由于MRI成像技術復雜，存在著相當大的數據運算量，傳統算法很難解決這一難題，本文以數據張量化的方法來打開突破口[1]。張量是目前能夠解決此類問題的主要方法之一，其最大的特點是能表達高維空間的數據特性，因此，如何將醫學圖像進行數據張量化，并將張量化方法進行有效的推廣，是本文的主旨所在[2]。

1張量代數

1.1張量的概念

張量（tensor）是一個多線性函數，其可以被直觀的理解為一個多維數組，它是由標量和向量擴展而來，標量、向量、矩陣是張量的特殊形式。標量可以理解成0階張量，向量是1階張量，而矩陣為2階張量[3]。

1.2高維圖像的張量表示

由張量概念可知，二維圖像可以用二階張量進行描述，而高維圖像需要用三階張量來表示。第一階表示了圖像的寬度，第二階表示圖像的長度，而第三階則表示圖像的矢量維數[4]。1.3張量模型首先要定義張量模型數學框架，為了表達簡潔，以二階張量的代數為例，它很容易拓展到高階張量。定義1.張量加法：給定尺寸均為的張量At和Bt，兩者之和Ct=At+Bt是相同尺寸的張量，其中表示二階數組Ct的（i;j）的元素，R則代表實數域。

2張量的主成份分析（Tensor-PCA）

在傳統運算（PCA）的基礎上，推廣到張量的變體TPCA[6]。

3實驗

3.1實驗來源

為了評估TPCA的性能，本文所使用的數據集是關于腦部醫學的開放數據庫：SimulatedBrainphantomDatabase,brainWeb網站提供了模擬數據操控臺可以對腦圖像數據進行仿真，此數據庫是用于對腦部MRI圖像做定量腦部分析所使用的數據集是一個3D立體腦部圖像卷。對于本實驗來講，數據集包涵有“groundtruth”，并且數據集將腦部結構分為10類。除了實驗背景外，選用剩下的9類進行分類實驗，同時選取數據集中的腦橫斷位上的第90幀圖像作為本實驗的仿真對象。3.2實驗數據處理在圖像的預處理上，首先利用Matlab軟件提取出數據集（BPD）的腦橫斷面第90幀，圖像尺寸大小為217181，像元大小為1mm1mm1mm，無噪聲，且識別率極高一幀圖像。

3.2.1PCA本文選取圖像數據集上的第90幀圖像，其每個位置的體素是一個樣本，每個樣本屬于9類結構中的一種，具體見表1所示。除去圖像背景外，一共有28043個像素點，現將圖像的像素點擴展成為33鄰域，此時圖像大小為3×3×28043，隨機抽取20%作為訓練集（5607個樣本），其余的作為測試樣本（22436個）。將所分好的數據帶入PCA進行特征提取處理，最后選用KNN（最近鄰）分類器對其分類。

3.2.2TPCA同理，在PCA的基礎上，我們將數據張量化，在每個標量上再找一個3×3的鄰域，經過實驗驗證后3×3的鄰域為最佳，大的鄰域效果并不明顯，反而鄰域過大會出現過飽和現象。擴充出來的位置用0來代替，此時數據集的大小為3×3×9×28043，這樣就完成了數據張量化，能夠帶入TPCA中進行特征提取。由于TPCA的運算是循環卷積，由于計算量過大，這里將數據帶入快速傅里葉變換（2DFFT）中進行，這樣就大大的降低了運算時間。

3.3實驗結果

首先由于NN分類器的局限性，還處理不了張量數據，這里我們將TPCA所得到結果進行切片操作，這樣就可以用NN分類器進行分類了。表2、表3分別給出了PCA和TPCA在NN分類器得到的分類混淆矩陣的對比。從主對角線上黑白程度上對混淆矩陣進行分析可以看出，在采用NN分類器進行分類時，PCA這種算法對于圖像上相似的類別仍有一定的錯分現象[7]。反觀TPCA算法，其主對角線黑色程度較深，說明準確率高，而且在其他分錯的類別上，百分比較低。

4結論

本文的主要工作是將磁共振大腦結構數據通過張量化的方法將其張量化，并利用基于循環卷積的張量模型，將大腦結構圖像進行有監督的分類對比，其主要工作如下：（1）首先BrainWeb在網上提供了一套現實的模擬腦MR的圖像卷，允許受控。（2）本文結合張量模型的優點將主成分分析(PCA)算法進行了全面的升級，傳統的PCA算法在對數據進行壓縮時，必須首先要將數據轉為向量，這對于圖像來說就破壞了數據的內在結構，因此從PCA的基本角度出發，我們將經典的PCA（主成分分析）算法升級為張量變體，稱為TPCA（張量主成分分析）。通過大量的仿真實驗后發現，基于循環卷積的張量代數框架下的TPCA算法具有非常優秀表現，實現了對張量數據的主成份分析，從而可以擴展到高維圖像數據。（3）結合TPCA算法獲得具有張量結構的數據特征，為了兼容已有的分類器，在TPCA的分類過程中還進行了切片操作，使得張量化后的數據能夠被一般分類器接受，這既體現了張量算法的優勢，也能夠擴展其應用范圍。

作者:廖亮；葉海昌；王新強 單位：中原工學院電子信息學院