(MATLAB) Volume Dataの補間

ボリュームデータのサブピクセル処理

interp3

[imaVOL, scaninfo, dcminfo] = loaddcm(loaddcmdir);

で読んだimaVOLについて補間を行なう。

interp3に渡すボリュームはdoubleでなければエラーが出る。 ついでに512x512x176のデータをまるまる使うとメモリが足りないので100x100x176の部分で切り出す。

smallVOL=dimaVOL(101:200,101:200,1:176);

受け皿を用意する。z方向だけ4倍の精度にしてみる。

[xi,yi,zi]=meshgrid(1:100,1:100,1:0.25:176);

interp3による補間

vi=interp3(smallVOL,xi,yi,zi);

座標に直交する断面の表示

slice(xi,yi,zi,vi,[20 60],20,[20 100])