CT画像の再構成（dicom画像)
アキシャル断面の連続スライス画像(dicom形式)を用いて、再構成を行い、サジタル断面とコロナル断面のスライスを取得したい。
そのためのPythonコードや手法などを教えていただきたいです。
0
1Answer
Qiitaで"dicom"で検索するとたくさん記事が見つかりました。
次の記事が代表的で網羅的な気がします。
https://qiita.com/fukuit/items/ed163f9b566baf3a6c3f
Windows Copilotに問い合わせ内容そのままプロンプトに投入すると次のような回答でした。
回答の正誤は確認していませんが、そこまで複雑ではなさそうなので、検索されてみてはいかがでしょうか
アキシャル断面のDICOM画像をサジタル断面やコロナル断面に再構成するためのPythonコードを提供します。以下の手順を参考にしてください。
DICOM画像の読み込み: pydicomライブラリを使用してDICOM画像を読み込みます。
画像の再構成: numpyライブラリを使用して画像データを再構成します。
スライスの取得: 再構成された画像からサジタル断面やコロナル断面のスライスを取得します。
以下は、サンプルコードです：
import pydicom
import numpy as np
# DICOM画像の読み込み
ds = pydicom.dcmread('path_to_your_dicom_file.dcm')
# 画像データの取得
image_data = ds.pixel_array
# 画像の再構成（例：サジタル断面）
sagittal_slice = image_data[:, :, int(image_data.shape[2] / 2)]
# 画像の再構成（例：コロナル断面）
coronal_slice = image_data[:, int(image_data.shape[1] / 2), :]
# スライスの表示（例）
import matplotlib.pyplot as plt
plt.imshow(sagittal_slice, cmap='gray')
plt.title('サジタル断面')
plt.show()
plt.imshow(coronal_slice, cmap='gray')
plt.title('コロナル断面')
plt.show()
1Like