.dcm格式文件軟件讀取及python處理的示例分析-創(chuàng)新互聯(lián)
這篇文章給大家分享的是有關(guān).dcm格式文件軟件讀取及python處理的示例分析的內(nèi)容。小編覺得挺實(shí)用的，因此分享給大家做個(gè)參考，一起跟隨小編過(guò)來(lái)看看吧。
要處理一些.DCM格式的焊接缺陷圖像，需要讀取和顯示.dcm格式的圖像。通過(guò)搜集資料收集到一些醫(yī)學(xué)影像，并通過(guò)pydicom模塊查看.dcm格式文件。
若要查看dcm格式文件，可下Echo viewer 進(jìn)行查看。
若用過(guò)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圖像卻無(wú)法進(jìn)行讀取。
感謝各位的閱讀！關(guān)于“.dcm格式文件軟件讀取及python處理的示例分析”這篇文章就分享到這里了，希望以上內(nèi)容可以對(duì)大家有一定的幫助，讓大家可以學(xué)到更多知識(shí)，如果覺得文章不錯(cuò)，可以把它分享出去讓更多的人看到吧！
另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn，海內(nèi)外云服務(wù)器15元起步，三天無(wú)理由+7*72小時(shí)售后在線，公司持有idc許可證，提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國(guó)服務(wù)器、虛擬主機(jī)、免備案服務(wù)器”等云主機(jī)租用服務(wù)以及企業(yè)上云的綜合解決方案，具有“安全穩(wěn)定、簡(jiǎn)單易用、服務(wù)可用性高、性價(jià)比高”等特點(diǎn)與優(yōu)勢(shì)，專為企業(yè)上云打造定制，能夠滿足用戶豐富、多元化的應(yīng)用場(chǎng)景需求。
本文題目：.dcm格式文件軟件讀取及python處理的示例分析-創(chuàng)新互聯(lián)
網(wǎng)站URL：http://m.jiaotiyi.com/article/cosgjj.html