西门子CT剂量报告/患者协议PixelData

Question

我试图从西门子的剂量报告中提取dicom PixelData，但它只包含零。使用GE剂量报告，我读取数据时使用pydicom或simpleITK没有任何问题。你知道为什么西门子的报告只包含零吗？

谢谢!

已经走到这一步了，但是角色是一团糟，我不知道还能做什么。

import dicom 
import numpy as np
import matplotlib.pyplot as plt

f="patient-protocol.dcm"
ds=dicom.read(f)  
overlay_px_data = ds[0x6000, 0x3000].value  
rows = ds[0x6000, 0x0010].value  
cols=ds[0x6000, 0x0011].value

arr=np.fromstring(overlay_px_data, dtype="uint8")  
arr=np.unpackbits(arr)  
arr=np.reshape(arr, (rows, cols))

plt.imshow(arr, cmap="gray")  
plt.show() 

dose-image

Answer 1

i_overlay = 1
n_bits = 8


# On (60xx,3000) are stored ovelays.
# First is (6000,3000), second (6002,3000), third (6004,3000),
# and so on.
dicom_tag1 = 0x6000 + 0*i_overlay

overlay_raw = data[0x6000,0x3000].value

# On (60xx,0010) and (60xx,0011) is stored overlay size
rows = data[0x6000,0x0010].value # rows = 512
cols = data[0x6000,0x0011].value # cols = 512

decoded_linear = np.zeros(len(overlay_raw)*n_bits)

# Decoding data. Each bit is stored as array element
for i in range(1,len(overlay_raw)):
    for k in range (0,n_bits):
        byte_as_int = overlay_raw[i]
        decoded_linear[i*n_bits + k] = (byte_as_int >> k) & 0b1


overlay = np.reshape(decoded_linear,[rows,cols])

plt.imshow(overlay)
plt.show()

找到了目前的解决方案。编辑了一点代码(删除ord())，一个numpy解码解决方案会更好/更好。

code source link