问题描述
我是python的新手。我正在尝试导入dicom文件。在许多教程中,以下代码用于使用pydicom导入dicom文件(例如https://www.kaggle.com/gzuidhof/full-preprocessing-tutorial):
INPUT_FOLDER = '../input/sample_images/'
patients = os.listdir(INPUT_FOLDER)
patients.sort()
将扫描结果加载到给定的文件夹路径中
def load_scan(path):
slices = [dicom.read_file(path + '/' + s) for s in os.listdir(path)]
slices.sort(key = lambda x: float(x.ImagePositionPatient[2]))
try:
slice_thickness = np.abs(slices[0].ImagePositionPatient[2] - slices[1].ImagePositionPatient[2])
except:
slice_thickness = np.abs(slices[0].SliceLocation - slices[1].SliceLocation)
for s in slices:
s.SliceThickness = slice_thickness
return slices
first_patient = load_scan(INPUT_FOLDER + patients[0])
每次运行此命令时,都会出现以下错误:
AttributeError Traceback (most recent call last)
<ipython-input-14-883826d71e95> in <module>()
----> 1 first_patient = load_scan(INPUT_FOLDER + patients[0])
2 frames
/usr/local/lib/python3.6/dist-packages/dicom/dataset.py in __getattr__(self,name)
255 if tag not in self:
256 raise AttributeError("Dataset does not have attribute "
--> 257 "'{0:s}'.".format(name))
258 else: # do have that dicom data_element
259 return self[tag].value
AttributeError: Dataset does not have attribute 'ImagePositionPatient'.
对每个人来说都很明显,但是这段代码导入了单独的dicom文件,这些文件一起构成了CT扫描,并将它们组合成单个数组“切片”。当我将“切片”定义为单个dicom文件而不是尝试将所有切片合并时,我可以编写slices.ImagePositionPatient并返回类似[x,y]的内容,因此原始文件肯定具有此属性。似乎当我将它们全部导入单个数组时,它将以某种方式丢失该属性。据我所知,我是唯一遇到此问题的人。谁能提供任何见识?
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)