python讀取dicom圖像示例(SimpleITK和dicom包實(shí)現(xiàn))
更新時間：2020年01月16日 17:15:03 作者：愿十四億神州盡舜堯
今天小編就為大家分享一篇python讀取dicom圖像示例(SimpleITK和dicom包實(shí)現(xiàn))，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧
1. 用SimpleITK讀取dicom序列：
import SimpleITK as sitk import numpy as np img_path='F:\\dataset\\pancreas\\Output\\thick\\original\\1' mask_path='F:\\dataset\\pancreas\\Output\\thick\\groundtruth\\1' reader = sitk.ImageSeriesReader() img_names = reader.GetGDCMSeriesFileNames(img_path) reader.SetFileNames(img_names) image = reader.Execute() image_array = sitk.GetArrayFromImage(image) # z, y, x reader = sitk.ImageSeriesReader() mask_names = reader.GetGDCMSeriesFileNames(mask_path) reader.SetFileNames(mask_names) mask = reader.Execute() mask_array = sitk.GetArrayFromImage(mask) # z, y, x
2. 用dicom讀取單張dicom圖像并顯示：
import dicom
import pylab
ds=dicom.read_file("F:\\dataset\\pancreas\\Output\\thick\\groundtruth\\1\\FILE0001_seg.dcm")
pixel_bytes = ds.PixelData
##CT值組成了一個矩陣
pix = ds.pixel_array
##讀取顯示圖片
pylab.imshow(ds.pixel_array, cmap=pylab.cm.bone)
pylab.show()
如果要對dicom圖像中的像素值進(jìn)行修改，繼續(xù)執(zhí)行以下代碼：
##修改圖片中的元素，不能直接使用data_array,需要轉(zhuǎn)換成PixelData
for n,val in enumerate(ds.pixel_array.flat): # example: zero anything < 300
if val < 300:
ds.pixel_array.flat[n]=0
ds.PixelData = ds.pixel_array.tostring()
ds.save_as("newfilename.dcm")
3. 此外，用pydicom也可讀取dicom圖像
以上這篇python讀取dicom圖像示例(SimpleITK和dicom包實(shí)現(xiàn))就是小編分享給大家的全部內(nèi)容了，希望能給大家一個參考，也希望大家多多支持腳本之家。
相關(guān)文章
python調(diào)用百度API實(shí)現(xiàn)人臉識別
這篇文章主要介紹了python調(diào)用百度API實(shí)現(xiàn)人臉識別，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-11-11
python實(shí)現(xiàn)爬取百度圖片的方法示例
這篇文章主要介紹了python實(shí)現(xiàn)爬取百度圖片的方法,涉及Python基于requests、urllib等模塊的百度圖片抓取相關(guān)操作技巧,需要的朋友可以參考下2019-07-07
利用Python實(shí)現(xiàn)智能合約的示例詳解
智能合約是一種由計算機(jī)程序編寫的自動化合約，它可以在沒有第三方干預(yù)的情況下執(zhí)行交易和契約條款。這篇文章主要介紹了如何利用Python實(shí)現(xiàn)智能合約，需要的可以參考一下2023-04-04
這篇文章主要介紹了一個超級簡單的python web程序,需要的朋友可以參考下2014-09-09
Python圖像處理庫PIL的ImageFont模塊使用介紹
這篇文章主要介紹了Python圖像處理庫PIL的ImageFont模塊使用介紹，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-02-02
自Transformer模型被應(yīng)用于計算機(jī)視覺領(lǐng)域后,圖像分割技術(shù)得到了進(jìn)一步的發(fā)展,但圖像分割操作復(fù)雜,特別是對張量的處理,涉及多種變換方法,其中,view/reshape用于改變數(shù)據(jù)形狀,而permute/transpose用于改變數(shù)據(jù)的維度順序2024-10-10
這篇文章主要介紹了用Python編寫分析Python程序性能的工具教程,包括了對內(nèi)存泄漏等常見問題的檢測方法,需要的朋友可以參考下2015-04-04
這篇文章主要介紹了Python合并2個字典成1個新字典的方法，本文通過實(shí)例代碼給大家分享9中方法,需要的朋友可以參考下2019-12-12
這篇文章主要介紹了五個與五子棋相關(guān)的Python文件,包括游戲邏輯、神經(jīng)網(wǎng)絡(luò)模型、訓(xùn)練代碼以及玩家對戰(zhàn)代碼,文中通過代碼介紹的非常詳細(xì),需要的朋友可以參考下2025-03-03
這篇文章主要給大家介紹了關(guān)于Matplotlib快速入門指南的相關(guān)資料,Matplotlib是一個非常強(qiáng)大的Python畫圖工具,支持跨平臺運(yùn)行,它不僅是Python常用的2D繪圖庫,同時它也提供了一部分3D繪圖接口,需要的朋友可以參考下2023-09-09