CTスキャンの断面再構成アルゴリズム ML-EM法の実装 - Qiita

はじめに 健康診断などでCTスキャンを受ける機会がありました。 どうやって断層の画像を再構成しているのか気になったので、調べて実装してみようと思います。 X線CTの画像再構成アルゴリズムは、従来はフィルタ補正逆投影法(FBP)が用いられていましたが、X線量を減らした場合には得られる画質に限界があるそうです。 そこで、計算時間はFBPより増えるものの、明瞭な画像が得られる統計的手法の一つ、ML-EM法をここでは試していきたいと思います。 ML-EM法は、実測の投影データをもとに再構成した画像から、再び投影データを作り、それができるだけ実測値に近づくように逐次更新していくアルゴリズムです。 参考 以下を参考にしました。 篠原 広行, 断層映像法の基礎 第32回 ML-EM法とOS-EM法 田中敏幸, 解説:特集 電磁界モデリングを用いる計測技術の基礎と新展開 X線CTの原理・現状とさらなる画像