2016年9月28日 更新

【医用画像の解析処理 ①】X線胸部画像

医用画像の解析処理についてご紹介します。参考にするのは、新潟大学の医用画像処理工学演習の授業ページ( http://www.clg.niigata-u.ac.jp/~medimg/practice_medical_imaging/imgproc_scion/1sampling/index.htm)です。

5,874 view お気に入り 0
医用画像には,CTやMRIのようなディジタル撮影装置を利用して撮影されるディジタル画像と,X線フィルムに代表されるアナログ画像があります。アナログ画像をコンピュータで解析するためにはディジタル化(=A/D変換:Analog-to-Digital conversion)を行わなければなりません。近年では、装置の精度管理を除くと、アナログ画像を撮影することはほとんどなくなりました。

ディジタル化の基本原理は,”標本化”と”量子化”で成りなっており、最初に標本化を行い,その後に量子化を行います。この順番を入れ替えることはできません。まず、標本化では、一定間隔ごとに画像の信号を切り分けます。次に、量子化では、切り分けた信号がどの程度の大きさであるかをしきい値によって分類します。これにより、X線画像がモニタ上に表示され、コンピュータ上で処理出来る形になります。

また、その画像がどの程度のFOV(Field of View:視野)で見られるかによって、画像のピクセルサイズを推測することができます。

胸部X線撮影の場合、大角(14×14inch)か、半切(14×17inch)で撮影されるため、512×512のマトリクスサイズであれば、大角の場合のピクセルサイズは、
(14(inch)×25.4(mm/inch))÷512=0.694...(mm)
となります、つまり、画面上で見ることの出来る1pixelの大きさは1mm以下であるということになります。

現在、病院で撮影されているX線撮影は、イメージングプレートを用いたCR(Computed Radiography)装置や、半導体素子を用いたFPD(Flat Panel Detector)装置が用いられ、撮影された画像は直接コンピュータ画面上で表示されるようになっています。

ここで、試しに胸部X線画像(新潟大学のページから拝借)の画像処理を行ってみます。
 (1085)

Lineを引き、Analyze->Plot Profileで、画像の輝度を調べると、以下のようになります。
(以下は信号が反転しており、白い部分が低信号として表示されてしまっています)
 (1087)

 (1088)

Lineの幅を広げると、プロファイルが滑らかになります。

次に、心胸郭比(CTR)を調べます。まず、心臓の右第2弓〜左第4弓までを直線で引き(1)、次に、胸部肋骨弓両下縁を直線で引きました。(2)
 (1090)

この時、以下の値がCTRとなります。


(1)の長さ/(2)の長さ=CTR(心胸郭比)


CTRの値が、0.5を越えると、心肥大であると考えられます。

tボタンをクリックし、ROI Managerで(1)と(2)で記録した直線の長さを調べます。ROI ManagerのMeasureボタンを押すと、以下の結果が表示されます。
 (1092)

Lengthを見ると、(1)が152、(2)が336であります。計算すると、

(1)の長さ/(2)の長さ=152/336=0.452...

となり、0.5を下回っているので、正常であると考えられます。
12 件

関連する記事 こんな記事も人気です♪

医用画像位置合わせの基礎⑤~Fiji pluginを用いた画像位置合わせ~

医用画像位置合わせの基礎⑤~Fiji pluginを用いた画像位置合わせ~

前回紹介した、LPixel社のImageJ registration pluginに引き続き、今回は、FijiというImageJの拡張版ソフトを使った画像位置合わせをご紹介します。
木田智士 | 22,496 view
YoutubeでImageJ  - インストール・セットアップ・プラグイン

YoutubeでImageJ - インストール・セットアップ・プラグイン

Fiji(ImageJにあらかじめ様々なプラグインがインストールされているもの)のインストールの仕方。新しいプラグインのインストールの方法を動画として紹介します。
北村 旭 | 3,428 view
ImageJ Plugin で数値計算をしてみる #5

ImageJ Plugin で数値計算をしてみる #5

非情報科学研究者 (特に生物系研究者) が ImageJ plugin を作るために超えるべき壁やTipsをまとめます。今回は、面積の平均値や中央値、分散を出力するためのコードをご紹介します。
湖城 恵 | 10,223 view
ImageJインストール(Mac)#ImageJの初心者から中級者まで#

ImageJインストール(Mac)#ImageJの初心者から中級者まで#

ImageJのインストールの仕方についてまとめました! 「ImageJってなに?」から「ImageJって面白い!!」へ!
ImageJ まとめ

ImageJ まとめ

「このページを見れば,ImageJ のこと全てが解決する」をモットーに,日々更新を行っています.

この記事のキーワード

この記事のキュレーター

エルピクセル編集部 エルピクセル編集部