MRI座標系と DICOMデータ2012/02/04松田圭司

MRI計測座標系• ＰＳ３．３－２００１翻訳　医療におけるデジタル画像と通信（ＤＩＣＯＭ）巻３：情報オブジェクト定義PDF

• ｘ軸は患者の左手側へ増加する。ｙ軸は患者の背面側に増加する。ｚ軸は患者の頭に向って増加する（ 169p)• 画像位置 (0020,0032),画像方向 (0020,0037)

Pixel-MRI座標系対応左上 (0020,0032)(RL) -60.0[mm](AP) -80.0[mm](IS) -35.0[mm]

Ｘ方向方向余弦（ 0020,0037) 前半3 個(RL) 1.00(AP) 0.00(IS) 0.00

Y 方向方向余弦（ 0020,0037) 後半3 個(RL) 0.00(AP) 1.00(IS) 0.00

Z方向　方向余弦Z 方向方向余弦画像 1 枚目左上（ 0020,0037)画像 n(n!=1) 枚目左上（ 0020,0037)の差を正規化する。(RL) 0.00(AP) 0.00(IS) 1.00

正規化は各値を、スライス間隔 (0018,0088)*(n-1)で割る。

方向余弦と鏡像変換• 方向余弦から 3x3の行列を作成し、その行列式（ determinant）を求める。

• 行列式が負となる場合は、鏡像変換を行っていることになる。• その場合は、画像の重ね順を反転する。

鏡像変換と変換行列• 画像の重ね順を変更した場合は、以下の行列を使用する。

座標変換まとめPixel座標系 LPS座標系

表示座標系 (サルは LIP)

あらかじめわかる

前出の変換行列

合成して逆行列を求めるPixel<->LPS座標系の変換が鏡像変換を含む場合は、スタック方向を変更（例 :S->Iを I->Sの順番で読み込む）。変換行列の行列式の符号で判断する。

表示座標系• ヒト (RPI)• サル (LIP) X:(R->L)

Y:(S->I)

Pixelと LPS座標変換Pixel->LPS変換行列

vsは、 voxelの幅（ =高、 =奥行）Lefttop_(L,P,S)は、 (0020,0037)の値。Delterminant>0の場合は、 1枚目。Determinant<0の場合は、最後。

上記の 4x4の逆行列を求めて、 Xlps,Ylps,Zlpsに (0,0,0)を代入すると、ピクセル空間における原点を求めることができる。