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

return User::create([
        'name' => $data['name'],
        'email' => $data['email'],
        'password' => bcrypt($data['password']),
        'stuID'=> $data['stuID'],
        'user_role'=> isset($data['user_role']) ? $data['user_role'] : null
    ]);

现在找到解决方案。编辑了一些代码（删除了ord（）），一个numpy解码解决方案会很棒/更好。

code source link

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

1 个答案: