PythonでDICOMヘッダファイルを読み込む例-Python チュートリアル-php.cn

PythonでDICOMヘッダファイルを読み込む例

不言

リリース： 2018-05-07 16:28:19

オリジナル

3671 人が閲覧しました

这篇文章主要介绍了关于python 读取DICOM头文件的实例，有着一定的参考价值，现在分享给大家，有需要的朋友可以参考一下

用dicompyler软件打开dicom图像，头文件如图所示：

当然也可以直接读取：

ds = dicom.read_file(&#39;H:\Data\data\\21662\\2.16.840.1.113662.2.0.105002416.1489146183.701\CT\\CT#0#21662#E7AB693D.dcm&#39;)
print ds
>>
(0008, 0008) Image Type       CS: [&#39;ORIGINAL&#39;, &#39;SECONDARY&#39;, &#39;AXIAL&#39;]
(0008, 0016) SOP Class UID      UI: CT Image Storage
(0008, 0018) SOP Instance UID     UI: 2.16.840.1.113662.2.0.105002416.1489146183.501.0
(0008, 0020) Study Date       DA: &#39;20170310&#39;
(0008, 0021) Series Date       DA: &#39;20170310&#39;
(0008, 0023) Content Date      DA: &#39;20060505&#39;
(0008, 0030) Study Time       TM: &#39;1942&#39;
(0008, 0031) Series Time       TM: &#39;1942&#39;
(0008, 0033) Content Time      TM: &#39;&#39;
(0008, 0050) Accession Number     SH: &#39;63071&#39;
(0008, 0060) Modality       CS: &#39;CT&#39;
(0008, 0070) Manufacturer      LO: &#39;NOMOS&#39;
(0008, 0090) Referring Physician&#39;s Name   PN: &#39;&#39;
(0008, 1010) Station Name      SH: &#39;&#39;
(0008, 1090) Manufacturer&#39;s Model Name   LO: &#39;CORVUS 6.4&#39;
(0010, 0010) Patient&#39;s Name      PN: &#39;*M32-2^CHENJUN^^^&#39;
(0010, 0020) Patient ID       LO: &#39;21662&#39;
(0010, 0030) Patient&#39;s Birth Date    DA: &#39;&#39;
(0010, 0040) Patient&#39;s Sex      CS: &#39;&#39;
(0018, 0050) Slice Thickness      DS: &#39;4.75&#39;
(0018, 0060) KVP         DS: &#39;0&#39;
(0018, 1020) Software Version(s)     LO: &#39;CORVUS 6.4&#39;
(0018, 5100) Patient Position     CS: &#39;HFS&#39;
(0020, 000d) Study Instance UID     UI: 2.16.840.1.113662.2.0.105002416.1489146183.701
(0020, 000e) Series Instance UID     UI: 2.16.840.1.113662.2.0.105002416.1489146183.501
(0020, 0010) Study ID       SH: &#39;63071&#39;
(0020, 0011) Series Number      IS: &#39;0&#39;
(0020, 0012) Acquisition Number     IS: &#39;0&#39;
(0020, 0013) Instance Number      IS: &#39;0&#39;
(0020, 0020) Patient Orientation     CS: [&#39;L&#39;, &#39;P&#39;]
(0020, 0032) Image Position (Patient)   DS: [&#39;0.73437356948853&#39;, &#39;0&#39;, &#39;3.25&#39;]
(0020, 0037) Image Orientation (Patient)   DS: [&#39;1&#39;, &#39;0&#39;, &#39;0&#39;, &#39;0&#39;, &#39;1&#39;, &#39;0&#39;]
(0020, 0052) Frame of Reference UID    UI: 2.16.840.1.113662.2.0.105002416.1489146183.601
(0020, 0060) Laterality       CS: &#39;&#39;
(0020, 1040) Position Reference Indicator  LO: &#39;&#39;
(0020, 1041) Slice Location      DS: &#39;3.25&#39;
(0028, 0002) Samples per Pixel     US: 1
(0028, 0004) Photometric Interpretation   CS: &#39;MONOCHROME2&#39;
(0028, 0010) Rows        US: 330
(0028, 0011) Columns        US: 339
(0028, 0030) Pixel Spacing      DS: [&#39;0.734375&#39;, &#39;0.734375&#39;]
(0028, 0100) Bits Allocated      US: 16
(0028, 0101) Bits Stored       US: 12
(0028, 0102) High Bit       US: 11
(0028, 0103) Pixel Representation    US: 0
(0028, 1052) Rescale Intercept     DS: &#39;-1024&#39;
(0028, 1053) Rescale Slope      DS: &#39;1&#39;
(7fe0, 0010) Pixel Data

ログイン後にコピー

发现一套图的第一张和最后一张的Slice Thickness往往和中间层的值是不同的。

Path = &#39;H:\Data\data\\21662\\2.16.840.1.113662.2.0.105002416.1489146183.701\CT\\a.dcm&#39;
slices = dicom.read_file(path)
spacing = slices.PixelSpacing

ログイン後にコピー

读取头文件信息方法：

Slices.上图头文件中的name列。

这时要注意的是，name大小写不变，去掉空格，去掉符号，比如括号。

举例：

origin = slices.SoftwareVersions
print origin
>>CORVUS 6.4

ログイン後にコピー

spacing = slices[1].PixelSpacing
print spacing
>>[&#39;0.734375&#39;, &#39;0.734375&#39;]

ログイン後にコピー