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

Office 指南

辦公室工作實用教學

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

Python SimpleITK 影像處理教學:Paste 貼上影像

介紹如何在 Python 中使用 SimpleITK 將一張影像貼在另外一張影像中。


SimpleITK 的 Paste 函數可以將一張影像貼在另外一張影像中指定的位置,將原本的影像蓋過去。

在影像對準的應用中,我們可以利用這樣的方式擴展基準影像的邊界,方便進行後續的影像對準工作。

原始影像

這裡我們以 Allen Mouse Brain Atlas 老鼠標準腦的影像為範例,先讀取原始影像,並檢查其影像切面。

import SimpleITK as sitk
import matplotlib.pyplot as plt

# 讀取 Allen Mouse Brain Atlas 影像
image = sitk.ReadImage("atlasVolume/atlasVolume.mhd")

# 轉為 NumPy 陣列
nda = sitk.GetArrayViewFromImage(image)

# 顯示影像切面
z_index = 228
slice = nda[z_index,:,:]
plt.imshow(slice)
plt.show()
原始影像
原始影像

從切面上可以看到影像內容非常貼近邊界,接下來我們要加寬影像周圍的邊界。

貼上影像

加寬影像周圍邊界的方法就是另外新增一張比較大的空白影像,然後將原始影像貼在新影像的正中央。

# 取得原始影像大小
imageSize = image.GetSize()

# 增添空白框邊
padding = 50

# 建立空白影像
expand = sitk.Image(imageSize[0] + padding * 2,
                    imageSize[1] + padding * 2,
                    imageSize[2] + padding * 2,
                    sitk.sitkUInt8)

# 貼上原始影像
result = sitk.Paste(expand, image, imageSize, (0, 0, 0), (padding, padding, padding))

# 轉為 NumPy 陣列
nda = sitk.GetArrayViewFromImage(result)

# 顯示影像切面
slice = nda[z_index + padding,:,:]
plt.imshow(slice)
plt.show()
貼上後的影像
貼上後的影像

這樣就完成加寬影像邊界的動作了,接著設定好影像的 voxel size 與位置,即可儲存影像。

# 設定 Voxel Size
result.SetSpacing(image.GetSpacing())

# 儲存影像
sitk.WriteImage(result, "atlasVolume/atlasVolume_expand50.mhd")

分類:Python 標籤: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