XCTスライス画像を​再構成して別方向から​のスライス画像を作り​たい．

あるものをx軸方向に0.5mm間隔でXCTをかけて，その様子を保存したdcmファイルの複数枚のスライス画像があります．

このスライス画像をmatlab内で読み込み再構成し，3D化を行いそれぞれy軸，z軸方向のスライス画像を作りたいのですが，どのようなスクリプトを書けばよいのでしょうか？

ご教示いただけますと幸いです．

Answer 1

[Shohei]

こんにちわ。

あまりdicomを扱ったことはないのですが、例えば以下のような感じでどうでしょうか。

% 読みたいファイルパスを設定
filePath = '*'

% スライスポジションを任意に設定
x = 1;
y = 2;
z = 3;

% Dicomの読み込み。3次元データと仮定してますが4次元以上かもしれません。
vol = dicomread(filePath);

% スライス画像の取り出し
slice_xy = vol(:,:,z);
slice_yz = squeeze(vol(x,:,:));
slice_xz = squeeze(vol(:,y,:));

単一スライス画像が格納された複数のdicomファイルならこちらかも

numSlices = 10;
vol = [];
for ii = 1 : numSlices
    vol(:,:,ii) = dicomread(filePath); % ファイルパスはループ中で適宜変えてください
end