コンピューター断層撮影(CTまたは猫スキャン)は臓器,骨,血管の3 d画像を生成するため,単なるレントゲン写真よりもはるかに高い診断性能を有します。しかしそれと同時に,人体に有害である放射線の被ばく量が増えるという代償が伴います。CTスキャンでは2 dのX線画像をコンピューターソフトウェアで重ねて組み合わせることで3 d画像を生成します。したがって1回の胸部CTスキャンの平均実効線量は7ミリシーベルト(mSv)となり,胸部X線の場合の0.02 mSvと比較すると350倍も高くなります。[1]放射線被ばくは,がんのリスクにつながります。ガイドラインでは,小児へのCTスキャンの線量は1.5 mSvに制限されています。
医工学研究者は,医師に必要な画像の鮮明さを提供しながら,放射線被ばく量を抑えることを目指しています。有望視されている1つのアプローチは,超低線量CTを使用する方法です。この方法では1回の胸部画像の平均実効線量が約0.13 mSvとなります。[1]超低線量CTスキャンの主なデメリットは,解像度が比較的低く,ノイズが多くなるため,臓器,脂肪,間質組織が見えにくくなることです(図1)。
図 1.超低線量CT(左)と従来のCT(右)の画質の比較
筆者は,超低線量CTスキャンを入力として使用しながら通常線量のCTスキャンと同程度の画質の画像を生成する,畳み込みニューラルネットワーク(CNN)回帰をベースとしたMATLAB®ソフトウェアシステムを開発しました。このシステムでは同程度の診断情報を提供しながら,患者の放射線被ばく量を最大95%低減します。
超解像処理とCNNの比較
低線量CTの画質を向上させる方法として,初めに超解像技術を適用しました。このアプローチではMATLABを使用し,CT画像を小さな局所領域の画像に分割して,低線量の領域と通常線量の領域をペアにすることで画像ディクショナリを作成しました。このシステムで新しい低線量画像を分析する際は,このディクショナリの中から類似した低線量領域を検索し,それに対応する通常線量の画像をユーザーに表示します。
この方法の有効性は,比較対象を取得するための大きなディクショナリがあるかどうかに左右されます。しかし,ディクショナリのサイズが大きくなると,システムリソースへの要求が増え,さらに重要なことには,類似した低線量領域を見つけるための検索時間が長くなります。CNNは学習に時間がかかりますが,新しい画像を提示されたときに結果を生成するスピードが,筆者の開発した超解像処理アプローチよりもはるかに速くなります。一例として,学習済みCNNでは1人の患者の結果を約20分で得られるのに対し,超解像処理を使用した場合は同等の結果を得るために約2時間かかりました。
超解像処理の欠点に対処するためにCNN回帰の検討を始めたところ,超解像処理の方が非常に優れているケースもありました。たとえば,診断する画像のパターンがディクショナリ内のある画像に非常に似ている場合,超解像処理によるアプローチは非常に正確な結果を生成します。そのため,CNN回帰と超解像処理を組み合わせたハイブリッドシステムの構築を計画しました。
画像の入手とCNNの構築
超低線量胸部CTスキャンの鮮明度を向上させるため,2つのCNNを使用するアプローチを採用しました。1つはCT画像の肺の領域を対象としたもので,もう1つは肺以外の領域を対象としたものです(図2)。CNNの学習に使用した画像データセットは,三重大学の研究チームから提供を受けました。このデータセットは12の画像ペアで構成されており,各画像ペアには同じ患者の通常線量のスキャンと超低線量のスキャンが含まれています(追加で超低線量の画像を撮影することで患者の被ばく量が増えるため,今回の研究では被験者を比較的少なく制限しました)。各画像は512 x 512ピクセルで,各スキャンには250枚の画像(断層)が含まれています。
図 2.超低線量CTの肺の領域と肺以外の領域で学習したCNN
CNNの初期構造は,超解像処理での結果を踏まえたものとしました。筆者が行った超解像処理の研究では,7 x 7サイズの局所領域で最適な結果が得られることがわかっています。そこで,ディープラーニングモデルではそのサイズの局所領域を使用することから始めました。その後,局所領域を5 x 128 x 128との間で変えながら実験し,それぞれが生成する結果の鮮明度を確認したうえで,肺の領域では32 x 32,肺以外の領域では64 x 64を設定しました。また,MATLABで入力サイズとフィルター,畳み込み層の数を変えた128種類の異なるCNNを評価しました。
CNNの学習と検証
交差検証を用いて11人の患者の画像を使用してモデルの学習を行い,残った1人の患者の画像を使用して評価しました。学習セットと評価画像を入れ替え,この手順を12回繰り返しました。過程を迅速に進めるため,并行计算工具箱™を使用して複数のNVIDIA®GeForceシリーズのGPUで並列に学習を行いました。学習の進行状況をモニタリングするため,深度学习工具箱™のモニタリング可視化機能で精度と損失をプロットしました(図3)。
図 3.深度学习工具箱で生成した学習の進行状況のプロット例
それぞれの超低線量テスト画像の結果は,対応する通常線量画像に照らして,二乗平均平方根(RMS)レベルと画質のメトリクスを測定するための構造類似性指数(SSIM)を用いて評価しました。
次のステップ
今回のCNNベースのシステムを実際の臨床環境で使用する計画を立てています。また,医療用画像管理システム(图片实现和通信;pac)サーバーにこのシステムを導入して,医療画像を簡単に保存したりアクセスしたりできる方法も模索しています。MATLABで医療画像処理ソフトウェアを開発するメリットは数多くあります。たとえば,以前に開発した他のMATLABベースのシステムでも行ったように,基礎となるアルゴリズムへのインターフェイスを作成して全体のパッケージを医師に配布するというプロセスが簡単になる環境です。
立命館大学は全校規模でMATLABと仿真软万博1manbetx件を利用している世界中で1000校を超える大学の1つです。总学术员工(发)ライセンスを導入しており,研究者,教職員,学生が最新リリースレベルの製品の共通のコンフィギュレーションを,講義室や自宅,研究室,現場などあらゆるところで使用しています。
