概述

pydicom是一個常用python DICOM parser。但是，沒有提供解析多幀圖的示例。本文結(jié)合相關(guān)函數(shù)和DICOM知識做一個簡單說明。

DICOM多幀數(shù)據(jù)存儲

DICOM標(biāo)準(zhǔn)中關(guān)于多幀數(shù)據(jù)存儲的最重要一部分說明是PS3.5 Annex A.4 A.4 Transfer Syntaxes For Encapsulation of Encoded Pixel Data。

無論何時，Pixel Data都存放在Pixel Data (7FE0,0010)中。有可能是直接存放的(native)，也有可能是被打包存放的(encapsulated）。常見的多幀圖一般采用打包存放的方式，特別是在數(shù)據(jù)存采用壓縮格式的時候。

在打包存放時，數(shù)據(jù)被拆分成很多數(shù)據(jù)段(Data Stream Fragment)，每個段保存成一個Tag為(FFFE,E000)的DICOM 項目(item)。其中第一個項目保存基礎(chǔ)偏移表(Basic Offset Table)，保存每個幀的起始偏移。然后，后面的數(shù)據(jù)段保存真正的數(shù)據(jù)。

**理論上，一幀（Frame）可以保存在多個數(shù)據(jù)段(Fragment)中。**但是，一個數(shù)據(jù)段不會保存兩幀的數(shù)據(jù)。實際應(yīng)用中，很少見到Frame跨數(shù)據(jù)段的情況，不過好的DICOM Parser應(yīng)該支持這種情況。這就是導(dǎo)致有些DICOM庫的接口看起來有點古怪的原因。

pydicom應(yīng)用

以pydicom 1.2為例，pydicom.encaps封裝了多幀有關(guān)的功能。解析最可能用到的幾個函數(shù)

多幀使用 pydicom.encaps來處理

pydicom.encaps.decode_data_sequence 返回每個(FFFE,E000)組成的item，形式是byte string的list。這樣一次拿到了所有fragment數(shù)據(jù)，可能比較耗內(nèi)存。注意理論上拿到的不是幀(Frame)

pydicom.encaps.generate_pixel_data 生成一個generator。每次迭代，返回做成一個fragments tuple，每個fragment是一個bytes，每個tuple是一個Frame。如果一個Fame保存在一個Fragment中。tuple就只有一個元素。

pydicom.encaps.generate_pixel_data_frame 迭代每次返回一個bytes，這個bytes代表了幀的全部數(shù)據(jù)。這應(yīng)該是最常用的函數(shù)。代碼實現(xiàn)就是把generate_pixel_data生成的tuple給拼接起來了。

這幾個方法的輸入都是pixel_array。相應(yīng)的Pydicom提供了用fp做參數(shù)的方法。

以上這篇Python解析多幀dicom數(shù)據(jù)詳解就是小編分享給大家的全部內(nèi)容了，希望能給大家一個參考，也希望大家多多支持腳本之家。

最新評論