• 跳至主要導覽
  • 跳至主要內容
  • 跳至主要資訊欄
Office 指南

Office 指南

辦公室工作實用教學

  • Excel
  • Word
  • PowerPoint
  • Windows
  • PowerShell
  • R

SimpleITK 使用 RescaleIntensityImageFilter 轉換影像數值範圍教學與範例

介紹如何使用 SimpleITK 的 RescaleIntensityImageFilter 對影像的像素值進行線性轉換。

相關文章:
SimpleITK 使用 IntensityWindowingImageFilter 轉換影像數值範圍教學與範例

首先從 DICOM 影像序列讀取 3D 影像資料:

import SimpleITK as sitk
import matplotlib.pyplot as plt

# DICOM 影像序列檔案所在目錄
dicom_folder = "/mnt/data/ct_image"

# 建立影像序列 Reader
reader = sitk.ImageSeriesReader()

# 取得所有 DICOM 影像的 Series IDs
series_ids = reader.GetGDCMSeriesIDs(dicom_folder)

# 取得指定 Series ID 的 DICOM 影像序列檔案名稱
dicom_filenames = reader.GetGDCMSeriesFileNames(dicom_folder, series_ids[0])

# 設定 DICOM 影像序列檔案名稱
reader.SetFileNames(dicom_filenames)

# 實際讀取影像內容
image = reader.Execute()

確認影像的大小:

# 顯示影像大小
size = image.GetSize()
print("Image size: {} x {} x {}".format(size[0], size[1], size[2]))
Image size: 512 x 512 x 361

檢查影像類型:

# 顯示影像類型
print(image.GetPixelIDTypeAsString())
16-bit signed integer

使用 StatisticsImageFilter 計算基本的影像統計資訊,取得像素值的分布範圍:

# 計算影像像素值範圍
statsFilter = sitk.StatisticsImageFilter()
statsFilter.Execute(image)
print("Minimum:", statsFilter.GetMinimum())
print("Maximum:", statsFilter.GetMaximum())
Minimum: -1000.0
Maximum: 2948.0

顯示影像切面以及像素值分布的直方圖:

# 顯示影像與像素值分布
nda = sitk.GetArrayViewFromImage(image)
fig, axs = plt.subplots(1, 2, figsize=(8, 4))
axs[0].imshow(nda[180,:,:], cmap = 'gray')
axs[0].set_title('Image')
axs[1].hist(nda.flatten(), bins=30)
axs[1].set_title('Histogram')
plt.show()
影像與像素值分布
影像與像素值分布

此處的原始影像是 16 位元的有號整數,接下來我們要使用 RescaleIntensityImageFilter將影像的像素值範圍轉換至 0 至 255,再將影像轉型為 8 位元的無號整數。

# 將影像數值範圍轉換為 0 至 255
resacleFilter = sitk.RescaleIntensityImageFilter()
resacleFilter.SetOutputMinimum(0)
resacleFilter.SetOutputMaximum(255)
imageRescale = resacleFilter.Execute(image)

# 將影像類型轉換為 8 位元的無號整數
castFilter = sitk.CastImageFilter()
castFilter.SetOutputPixelType(sitk.sitkUInt8)
imageUint8 = castFilter.Execute(imageRescale)

計算新的影像像素值範圍:

# 計算影像像素值範圍
statsFilter = sitk.StatisticsImageFilter()
statsFilter.Execute(imageUint8)
print("Minimum:", statsFilter.GetMinimum())
print("Maximum:", statsFilter.GetMaximum())
Minimum: 0.0
Maximum: 255.0

顯示新影像與像素值分布直方圖:

# 顯示影像與像素值分布
nda = sitk.GetArrayViewFromImage(imageUint8)
fig, axs = plt.subplots(1, 2, figsize=(8, 4))
axs[0].imshow(nda[180,:,:], cmap = 'gray')
axs[0].set_title('Image')
axs[1].hist(nda.flatten(), bins=30)
axs[1].set_title('Histogram')
plt.show()
影像與像素值分布
影像與像素值分布

RescaleIntensityImageFilter 會自動計算原始影像像素值範圍,將整個範圍以線性轉換的方式轉換至指定的區間,若想要以窗函數自行篩選要保留的像素值範圍,可以改用 IntensityWindowingImageFilter。

分類:Python 標籤:DICOM, ITK

讀者互動方式

發佈留言 取消回覆

發佈留言必須填寫的電子郵件地址不會公開。 必填欄位標示為 *

主要資訊欄

搜尋

近期文章

  • C++ 語言使用 Crypto++ 實作 RSA 數位簽章教學與範例
  • C++ 語言使用 Crypto++ 實作 RSA-OAEP 搭配 SHA256 加密教學與範例
  • C++ 語言使用 Crypto++ 實作 AES 加密、解密、認證加密教學與範例
  • C++ 語言使用 Crypto++ 實作 MD5、SHA1、SHA2、BLAKE2 雜湊教學與範例
  • Ubuntu Linux 安裝、使用 Crypto++ 加密函式庫教學與範例
  • C 語言使用 OpenSSL 實作橢圓曲線 ECDH 金鑰交換教學與範例
  • Python 以 eciespy 實作 ECC 非對稱式加密方法教學與範例
  • C 語言使用 OpenSSL 實作 PBKDF2 教學與範例

推薦網站

  • Udemy 線上教學課程
  • Coursera 線上教學課程

關注本站

  • 電子郵件
  • Facebook

公益

  • 家扶基金會
  • 台灣世界展望會
  • Yahoo 奇摩公益
  • igiving 公益網
  • 兒福聯盟

Copyright © 2021 · Office Guide