使用Python解決常見格式圖像讀取nii,dicom,mhd

更新時間：2022年04月18日 12:38:12 作者：派大大大星?

這篇文章主要介紹了使用Python解決常見格式圖像讀取nii,dicom,mhd,下文具體操作過程需要的小伙伴可以參考一下

1. raw，mhd 格式醫(yī)學圖像數(shù)據(jù)轉換

raw+mhd格式是常見的一種醫(yī)學圖像格式，每一個病人的數(shù)據(jù)包含一個mhd文件和一個同名的raw文件，mhd即meta header data，數(shù)據(jù)頭部信息，而raw存儲了像素信息。方法需要使用的SimpleITK庫，我們需要在自己的Python環(huán)境中安裝對應的庫

pip install SimpleITK

mhd+raw的數(shù)據(jù)往往是三維體數(shù)據(jù)，我們可以從mhd文件中讀取數(shù)據(jù)的具體信息如圖像大小、切片大小、像素大小等信息。使用sitk.ReadImage() 即可讀取圖像，使用sitk.GetArrayFromImage() 可獲取圖像矩陣。

具體代碼如下：

import os
import SimpleITK as sitk
import matplotlib.pyplot as plt
from natsort import natsorted
from tqdm import tqdm
import cv2
import numpy as np

mhd_path = './xxx.mhd'  # mhd文件需和同名raw文件放在同一個文件夾
data = sitk.ReadImage(mhd_path)  # 讀取mhd文件
# print(data)
spacing = data.GetSpacing()  # 獲得spacing大小
img_data = sitk.GetArrayFromImage(data)  # 獲得圖像矩陣
print(img_data.shape)

# 圖像數(shù)據(jù)轉換為npy保存
np.save('img_data.npy', np.array(img_data))

# 將圖像轉為png格式保存
for i in range(img_data.shape[0]):
    cv2.imwrite('./{}.png'.format(i), img_data[i, :, :])

2. dicom格式醫(yī)學圖像數(shù)據(jù)轉換

DICOM（Digital Imaging and Communications in Medicine） 即醫(yī)學數(shù)字成像和通信，是醫(yī)學圖像和相關信息的國際標準（ISO 12052）。DICOM被廣泛應用于放射醫(yī)療，心血管成像以及放射診療診斷設備（X射線，CT，核磁共振，超聲等），并且在眼科和牙科等其它醫(yī)學領域得到越來越深入廣泛的應用。所有患者的醫(yī)學圖像都以 DICOM 文件格式進行存儲。使用Python讀取dicom圖像可以使用pydicom庫和SimpleITK庫來完成。由于在醫(yī)學圖像處理領域中，對不同的圖像需要使用不同的窗寬窗位導出圖像，在下面代碼中，我增加了對應的代碼。

pip install SimpleITK
pip install pydicom

使用pydicom方法具體代碼如下：

import os
import SimpleITK as sitk
import matplotlib.pyplot as plt
from natsort import natsorted
from tqdm import tqdm
import cv2
import pydicom
import numpy as np

#調整圖像窗寬窗位
def window_transform(ct_array, window_width, window_center, normal=False):
    min_window = float(window_center) - 0.5*float(window_width)
    new_img = (ct_array - min_window) / float(window_width)
    new_img[new_img < 0] = 0
    new_img[new_img > 1] = 1
    if not normal:
        new_img = (new_img * 255).astype('uint8')
    return new_img


img_path = 'xx.dcm'
#此時讀取的是所有dicom圖像信息
image = pydicom.read_file(img_path)
#獲得圖像矩陣
image_data = image.pixel_array
#獲得dicom中的患者信息、圖像信息等
information['PatientID'] = image.PatientID
information['PatientName'] = image.PatientName
information['PatientBirthDate'] = image.PatientBirthDate
information['PatientSex'] = image.PatientSex
window_width = 1000
window_center = 30
image_data = window_transform(image_data, window_width, window_center, normal=False)
cv2.imwrite('./img.png',image_data)

使用SimpleITK方法具體代碼如下：

import os
import SimpleITK as sitk
import matplotlib.pyplot as plt
from natsort import natsorted
from tqdm import tqdm
import cv2
import pydicom
import numpy as np

img_path = 'xx.dcm'
#此時讀取的是所有dicom圖像信息
image = sitk.ReadImage(img_path)
# 轉為值為0-255的灰度圖
rescalFilt = sitk.RescaleIntensityImageFilter()
rescalFilt.SetOutputMaximum(255)
rescalFilt.SetOutputMinimum(0)
image = rescalFilt.Execute(image)
image_data = sitk.GetArrayFromImage(image_data)
image_data = np.squeeze(dicom_data)
cv2.imwrite('./img.png',image_data)

3. nii格式醫(yī)學圖像轉換

醫(yī)學影像早期使用的是DICOM標準，基本上各家廠商都會使用符合DICOM標準的產品，但是這個標準對于數(shù)據(jù)分析并不方便。在神經影像興起時就誕生了各種各樣的數(shù)據(jù)存儲標準，比如analyze。后為了便于學術交流，NIH拉著其他一些組織的專家成立了工作組，制定了新的神經影像的數(shù)據(jù)存儲標準，稱為NIFTI。讀取nii后綴的NIFTI格式圖像需要安裝nibabel庫

pip install nibabel

使用nib.load() 函數(shù)即可讀取圖像數(shù)據(jù)data, 其中dataobj就是圖像矩陣

import os 
import numpy as np 
import nibabel as nib
import matplotlib.pyplot as plt
from glob import glob

img_path = './xxx.nii'
image = nib.load(img_path)    
image_data = image.dataobj
for i in range(0, dataobj.shape[2]): 
    cv2.imwrite('./{}.png'.format(i), img_data[i, :, :])

到此這篇關于使用Python解決常見格式圖像讀取nii,dicom,mhd的文章就介紹到這了,更多相關Python圖像讀取內容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關文章希望大家以后多多支持腳本之家！

您可能感興趣的文章: