快捷導(dǎo)航

.dcm格式文件軟件讀取及python處理詳解

更新時(shí)間：2020年01月16日 17:03:43 作者：西涯俠（Python）

今天小編就為大家分享一篇.dcm格式文件軟件讀取及python處理詳解，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧

要處理一些.DCM格式的焊接缺陷圖像，需要讀取和顯示.dcm格式的圖像。通過搜集資料收集到一些醫(yī)學(xué)影像，并通過pydicom模塊查看.dcm格式文件。

若要查看dcm格式文件，可下Echo viewer 進(jìn)行查看。

若用過pycharm進(jìn)行處理，可選用如下的代碼：

# -*-coding:utf-8-*-
import cv2
import numpy
import dicom
from matplotlib import pyplot as plt

dcm = dicom.read_file("dcm")
dcm.image = dcm.pixel_array * dcm.RescaleSlope + dcm.RescaleIntercept

slices = []
slices.append(dcm)
img = slices[int(len(slices) / 2)].image.copy()
ret, img = cv2.threshold(img, 90, 3071, cv2.THRESH_BINARY)
img = numpy.uint8(img)

im2, contours, _ = cv2.findContours(img, cv2.RETR_LIST, cv2.CHAIN_APPROX_SIMPLE)
mask = numpy.zeros(img.shape, numpy.uint8)
for contour in contours:
  cv2.fillPoly(mask, [contour], 255)
img[(mask > 0)] = 255

kernel = cv2.getStructuringElement(cv2.MORPH_ELLIPSE, (2, 2))
img = cv2.morphologyEx(img, cv2.MORPH_OPEN, kernel)

img2 = slices[int(len(slices) / 2)].image.copy()
img2[(img == 0)] = -2000

plt.figure(figsize=(12, 12))
plt.subplot(131)
plt.imshow(slices[int(len(slices) / 2)].image, 'gray')
plt.title('Original')
plt.subplot(132)
plt.imshow(img, 'gray')
plt.title('Mask')
plt.subplot(133)
plt.imshow(img2, 'gray')
plt.title('Result')
plt.show()

也可用如下代碼：

import pydicom
import os
import numpy
from matplotlib import pyplot, cm
# 用lstFilesDCM作為存放DICOM files的列表
PathDicom = "dicom/2" #與python文件同一個(gè)目錄下的文件夾
lstFilesDCM = []
for dirName,subdirList,fileList in os.walk(PathDicom):
  for filename in fileList:
    if ".dcm" in filename.lower(): #判斷文件是否為dicom文件
      print(filename)
      lstFilesDCM.append(os.path.join(dirName,filename)) # 加入到列表中
## 將第一張圖片作為參考圖
RefDs = pydicom.read_file(lstFilesDCM[0])  #讀取第一張dicom圖片
# 建立三維數(shù)組
ConstPixelDims = (int(RefDs.Rows),int(RefDs.Columns),len(lstFilesDCM)) # 得到spacing值 (mm為單位)
ConstPixelSpacing = (float(RefDs.PixelSpacing[0]), float(RefDs.PixelSpacing[1]), float(RefDs.SliceThickness))
# 三維數(shù)據(jù)
x = numpy.arange(0.0, (ConstPixelDims[0]+1)*ConstPixelSpacing[0], ConstPixelSpacing[0]) # 0到（第一個(gè)維數(shù)加一*像素間的間隔），步長(zhǎng)為constpixelSpacing
y = numpy.arange(0.0, (ConstPixelDims[1]+1)*ConstPixelSpacing[1], ConstPixelSpacing[1]) #
z = numpy.arange(0.0, (ConstPixelDims[2]+1)*ConstPixelSpacing[2], ConstPixelSpacing[2]) #
ArrayDicom = numpy.zeros(ConstPixelDims, dtype=RefDs.pixel_array.dtype)
for filenameDCM in lstFilesDCM:
  ds = pydicom.read_file(filenameDCM)
  ArrayDicom[:, :, lstFilesDCM.index(filenameDCM)] = ds.pixel_array # 軸狀面顯示
  pyplot.figure(dpi=300)
  pyplot.axes().set_aspect('equal', 'datalim')
  pyplot.set_cmap(pyplot.gray())
  pyplot.pcolormesh(x, y, numpy.flipud(ArrayDicom[:, :, 2])) # 第三個(gè)維度表示現(xiàn)在展示的是第幾層
  pyplot.show()

這兩個(gè)代碼都是可以進(jìn)行讀取的。但是不知道為什么在焊接檢測(cè)中的dcm圖像卻無法進(jìn)行讀取。

以上這篇.dcm格式文件軟件讀取及python處理詳解就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。

您可能感興趣的文章: