概述
pydicom是一個常用python DICOM parser。但是,沒有提供解析多幀圖的示例。本文結合相關函數和DICOM知識做一個簡單說明。
DICOM多幀數據存儲
DICOM標準中關于多幀數據存儲的最重要一部分說明是PS3.5 Annex A.4 A.4 Transfer Syntaxes For Encapsulation of Encoded Pixel Data。
無論何時,Pixel Data都存放在Pixel Data (7FE0,0010)中。有可能是直接存放的(native),也有可能是被打包存放的(encapsulated)。常見的多幀圖一般采用打包存放的方式,特別是在數據存采用壓縮格式的時候。
在打包存放時,數據被拆分成很多數據段(Data Stream Fragment),每個段保存成一個Tag為(FFFE,E000)的DICOM 項目(item)。其中第一個項目保存基礎偏移表(Basic Offset Table),保存每個幀的起始偏移。然后,后面的數據段保存真正的數據。
**理論上,一幀(Frame)可以保存在多個數據段(Fragment)中。**但是,一個數據段不會保存兩幀的數據。實際應用中,很少見到Frame跨數據段的情況,不過好的DICOM Parser應該支持這種情況。這就是導致有些DICOM庫的接口看起來有點古怪的原因。
pydicom應用
以pydicom 1.2為例,pydicom.encaps封裝了多幀有關的功能。解析最可能用到的幾個函數
多幀使用 pydicom.encaps來處理
pydicom.encaps.decode_data_sequence 返回每個(FFFE,E000)組成的item,形式是byte string的list。這樣一次拿到了所有fragment數據,可能比較耗內存。注意理論上拿到的不是幀(Frame)
pydicom.encaps.generate_pixel_data 生成一個generator。每次迭代,返回做成一個fragments tuple,每個fragment是一個bytes,每個tuple是一個Frame。如果一個Fame保存在一個Fragment中。tuple就只有一個元素。
pydicom.encaps.generate_pixel_data_frame 迭代每次返回一個bytes,這個bytes代表了幀的全部數據。這應該是最常用的函數。代碼實現就是把generate_pixel_data生成的tuple給拼接起來了。
這幾個方法的輸入都是pixel_array。相應的Pydicom提供了用fp做參數的方法。
以上這篇Python解析多幀dicom數據詳解就是小編分享給大家的全部內容了,希望能給大家一個參考,也希望大家多多支持服務器之家。
原文鏈接:https://blog.csdn.net/pkrobbie/article/details/87827145