特許 7631505 Ｘ線イメージングに関連する方法及びシステム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2025-02-07

(45)【発行日】2025-02-18

(54)【発明の名称】Ｘ線イメージングに関連する方法及びシステム

(51)【国際特許分類】

   A61B 6/03 20060101AFI20250210BHJP

【ＦＩ】

A61B6/03 573

【請求項の数】 32

(21)【出願番号】P 2023511887

(86)(22)【出願日】2021-07-06

(65)【公表番号】

(43)【公表日】2023-09-13

(86)【国際出願番号】 SE2021050685

(87)【国際公開番号】W WO2022039643

(87)【国際公開日】2022-02-24

【審査請求日】2023-04-17

(31)【優先権主張番号】63/067,656

(32)【優先日】2020-08-19

(33)【優先権主張国・地域又は機関】US

(73)【特許権者】

【識別番号】512120638

【氏名又は名称】プリズマティック、センサーズ、アクチボラグ

【氏名又は名称原語表記】ＰＲＩＳＭＡＴＩＣ ＳＥＮＳＯＲＳ ＡＢ

(74)【代理人】

【識別番号】100105588

【弁理士】

【氏名又は名称】小倉 博

(74)【代理人】

【識別番号】100129779

【弁理士】

【氏名又は名称】黒川 俊久

(74)【代理人】

【識別番号】100151286

【弁理士】

【氏名又は名称】澤木 亮一

(72)【発明者】

【氏名】エギザヴァル，アルマ

(72)【発明者】

【氏名】グロンバーグ，フレドリック

(72)【発明者】

【氏名】パーション，マット

【審査官】佐藤 賢斗

(56)【参考文献】

【文献】特表２０１８－５１５１６０（ＪＰ，Ａ）

【文献】特開２０２０－１０３５７１（ＪＰ，Ａ）

【文献】特開２０２０－０９９６６２（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ６／００ － ６／５８

(57)【特許請求の範囲】

【請求項1】

エネルギー分解Ｘ線データに基づく画像再構成方法であって、
エネルギー分解Ｘ線データの少なくとも１つの表現を収集するステップ（Ｓ１）と、
エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を実行して、少なくとも２つの元の基底画像表現セットを生成するステップ（Ｓ２）と、
元の基底画像表現セットの少なくとも２つから少なくとも２つの基底画像表現を取得又は選択するステップ（Ｓ３）と、
取得又は選択された少なくとも２つの基底画像表現を機械学習に基づいてデータ処理して、出力画像データの少なくとも１つの表現を生成するステップ（Ｓ４）であって、前記データ処理は、機械学習に基づいて、前記少なくとも２つの基底物質分解から得られた前記取得又は選択された少なくとも２つの基底画像表現を、前記出力画像データの少なくとも１つの表現に対応付けるように構成されている、ステップ（Ｓ４）と、
を含む、方法。

【請求項2】

機械学習に基づくデータ処理は、訓練されたニューラルネットワーク、サポートベクターマシン、または決定木に基づいて、取得又は選択された基底画像表現を出力画像データの少なくとも１つの表現にマッピングするよう構成される、請求項１に記載の方法。

【請求項3】

訓練されたニューラルネットワーク、サポートベクターマシン、または決定木は、測定された患者画像データ若しくはファントム画像データの少なくとも１つ、又は数値ファントムに基づくシミュレーション画像データで訓練される、請求項２に記載の方法。

【請求項4】

訓練されたニューラルネットワーク、サポートベクターマシン、または決定木は、第１のタイプの画像で訓練され、転移学習を使用して第２のタイプの画像で使用するように適合される、請求項２に記載の方法。

【請求項5】

訓練されたニューラルネットワーク、サポートベクターマシン、または決定木は、非エネルギー分解ＣＴ画像で訓練され、転移学習を使用してエネルギー分解光子計数ＣＴ画像で使用されるように適合される、請求項４に記載の方法。

【請求項6】

エネルギー分解Ｘ線データの少なくとも１つの表現は、Ｘ線検出器または中間記憶装置に由来し、エネルギー分解Ｘ線画像データおよび／またはエネルギー分解Ｘ線測定データを含む、請求項１に記載の方法。

【請求項7】

Ｘ線検出器が、マルチビン光子計数Ｘ線検出器とも呼ばれるエネルギー分解光子計数検出器である、請求項６に記載の方法。

【請求項8】

出力画像データの少なくとも１つの表現は、基底画像表現の出力セットと呼ばれる基底画像表現のセットを含み、出力セットにおける基底画像表現の数は、少なくとも２つの元の基底画像表現セットのいずれかにおける基底画像表現の最大数より大きい、請求項１に記載の方法。

【請求項9】

出力セットの基底画像表現は、複数の基底関数を有し、元の基底画像表現セットの各々は、それぞれ対応する数の基底関数を有し、出力セットの基底画像表現に対応する基底関数は、少なくとも２つの元の基底画像表現セットに対応する基底関数のサブセットであるか、またはそれに等しい、請求項８に記載の方法。

【請求項10】

出力セットの基底画像表現は、複数の基底関数を有し、元の基底画像表現セットの各々はそれぞれ対応する数の基底関数を有し、取得又は選択された基底画像表現の数は、出力セットの基底画像表現に対応する基底関数の数と同じかより大きい、請求項８に記載の方法。

【請求項11】

元の基底画像表現セットが２基底分解に基づき、基底画像表現の出力セットが３基底分解に基づく、請求項８に記載の方法。

【請求項12】

機械学習に基づくデータ処理がニューラル・ネットワークによって実行され、ニューラル・ネットワークが入力データとしていくつかの２基底画像表現セットと出力データとして３基底画像表現セットで訓練される、請求項１１に記載の方法。

【請求項13】

２基底画像表現セット及び３基底画像表現セットが、測定された患者画像データ又はファントム画像データから、又は数値ファントムに基づくシミュレートされた画像データから生成される、請求項１２に記載の方法。

【請求項14】

取得又は選択された基底画像表現を機械学習に基づいてデータ処理するステップ（Ｓ４）は、
取得又は選択された基底画像表現を、非基底画像データの少なくとも１つの表現と結合するステップと、
機械学習に基づくデータ処理によって、取得又は選択された基底画像表現を、非基底画像データの少なくとも１つの表現と一緒に処理するステップと、
を含む、請求項１に記載の方法。

【請求項15】

取得又は選択された基底画像表現を機械学習に基づいてデータ処理するステップ（Ｓ４）は、
取得又は選択された基底画像表現を少なくとも１つの事前情報の表現と組み合わせるステップと、
機械学習に基づくデータ処理によって、取得又は選択された基底画像表現を、少なくとも１つの事前情報の表現と共に処理するステップと、
を含む、請求項１に記載の方法。

【請求項16】

取得又は選択された基底画像表現を機械学習に基づいてデータ処理するステップ（Ｓ４）は、
取得又は選択された基底画像表現を、収集システム内の物理プロセスの少なくとも１つの計算表現と組み合わせるステップと、
機械学習に基づくデータ処理によって、物理的プロセスの少なくとも１つの計算表現を用いて、取得又は選択された基底画像表現を処理するステップと、
を含む、請求項１に記載の方法。

【請求項17】

取得又は選択された基底画像表現を機械学習に基づいてデータ処理するステップ（Ｓ４）は、
取得又は選択された基底画像表現を、測定データの統計モデルの少なくとも１つの計算表現と組み合わせ、
機械学習に基づくデータ処理によって、統計モデルの少なくとも１つの計算表現を用いて、取得又は選択された基底画像表現を処理するステップと、
を含む、請求項１に記載の方法。

【請求項18】

出力画像データの少なくとも１つの表現が、基底画像表現及び／又は物質分解再構成画像の出力セットを含む、請求項１に記載の方法。

【請求項19】

機械学習に基づくデータ処理は、学習データに基づいて構成されたディープラーニング後処理を含み、ディープラーニング後処理は、入力として偏った物質分解サイノグラムを受け取り、出力として偏りのない物質分解サイノグラム又は偏りのない物質分解再構成画像を提供するよう構成される、請求項１に記載の方法。

【請求項20】

サイノグラムのセットが、２つ以上の基底画像表現を生成するための２基底分解の入力として提供され、ディープラーニングニューラルネットワークが、得られた基底画像表現を３物質分解に変換するために使用される、請求項１９に記載の方法。

【請求項21】

エネルギー分解Ｘ線データに基づく画像再構成のためのシステム（３０、２００）であって、
エネルギー分解Ｘ線データの少なくとも１つの表現を受信するための入力（３１、２４０）と、
エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を実行して、基底画像表現の少なくとも２つのセットを生成するように構成された基底物質分解システム（３２；２１０、２２５）と、
入力－出力対の訓練データに基づいて決定された構成を有し、基底画像表現のセットの少なくとも２つから少なくとも２つの基底画像表現を入力として受信し、受信した少なくとも２つの基底画像表現を機械学習に基づいてデータ処理して、出力画像データの少なくとも１つの表現を生成するように構成された訓練された機械学習システム（３３、２１０、２２５）であって、前記機械学習システム（３３、２１０、２２５）は、前記少なくとも２つの基底物質分解から得られた前記受信した少なくとも２つの基底画像表現を、前記出力画像データの少なくとも１つの表現に対応付けるように構成された、機械学習システム（３３、２１０、２２５）と、
出力画像データの少なくとも１つの表現を出力する出力（３４、２４０）と、
を備える画像再構成システム（３０、２００）。

【請求項22】

機械学習システム（３３、２１０、２２５）は、ニューラル・ネットワーク、サポートベクターマシン、又は決定木ベースの実装に基づく、請求項２１に記載の画像再構成システム。

【請求項23】

機械学習システム（３３、２１０、２２５）は、測定された患者画像データ又はファントム画像データのうちの少なくとも１つで、または数値ファントムに基づくシミュレーション画像データで訓練される、請求項２２に記載の画像再構成システム。

【請求項24】

機械学習システム（３３、２１０、２２５）は、第１のタイプの画像で訓練され、転移学習を使用して第２のタイプの画像で使用されるように適合される、請求項２２に記載の画像再構成システム。

【請求項25】

機械学習システム（３３、２１０、２２５）は、非エネルギー分解ＣＴ画像で訓練され、転移学習を使用してエネルギー分解光子計数ＣＴ画像で使用されるように適合される、請求項２４に記載の画像再構成システム。

【請求項26】

エネルギー分解されたＸ線データの少なくとも１つの表現が、Ｘ線検出器（２０）又は中間記憶装置（２２０）から生じ、エネルギー分解Ｘ線画像データ及び／又はエネルギー分解Ｘ線測定データを含む、請求項２２に記載の画像再構成システム。

【請求項27】

Ｘ線検出器（２０）が、マルチビン光子計数Ｘ線検出器とも呼ばれるエネルギー分解光子計数検出器である、請求項２６に記載の画像再構成システム。

【請求項28】

出力画像データの少なくとも１つの表現は、基底画像表現の出力セットと呼ばれる基底画像表現のセットを含み、出力セットにおける基底画像表現の数は、少なくとも２つの元の基底画像表現セットのいずれかにおける基底画像表現の最大数より大きい、請求項２２に記載の画像再構成システム。

【請求項29】

出力セットの基底画像表現は、複数の基底関数を有し、元の基底画像表現セットの各々は、それぞれ対応する数の基底関数を有し、出力セットの基底画像表現に対応する基底関数は、少なくとも２つの元の基底画像表現セットに対応する基底関数のサブセットであるか、またはそれに等しい、請求項２８に記載の画像再構成システム。

【請求項30】

出力画像データの少なくとも１つの表現が、基底画像表現及び／又は物質分解再構成画像の出力セットを含む、請求項２２に記載の画像再構成システム。

【請求項31】

前記システムが少なくとも１つのプロセッサ及びメモリを備え、メモリが、少なくとも１つのプロセッサによって実行されると、少なくとも１つのプロセッサに画像再構成を実行させる命令を備える、請求項２２に記載の画像再構成システム。

【請求項32】

基底物質分解システムが、基底画像表現の少なくとも２つのセットを生成するように構成された処理回路を備え、訓練された機械学習システムが、出力の少なくとも１つの表現を生成するように構成された処理回路を備える、請求項２２に記載の画像再構成システム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、Ｘ線技術とＸ線イメージング、及び対応するイメージング再構成とイメージングタスクに関連する。特に、本開示は、画像再構成のための方法及びシステム、対応するＸ線イメージングシステム、ならびに対応するコンピュータプログラム及びコンピュータプログラム製品に関する。

【背景技術】

【0002】

Ｘ線イメージング（Ｘ線撮像）などの放射線イメージングは、医療用途や非破壊検査で長年使用されてきている。

【0003】

通常、Ｘ線イメージング装置は、Ｘ線源と、１つまたは多数の検出器要素（Ｘ線の強度／蛍光を測定する独立した手段）からなる複数の検出器からなるＸ線検出器アレイを含む。Ｘ線源から放射されたＸ線は、被写体を透過して検出器アレイに入射し、検出されたＸ線は検出器アレイに記録される。Ｘ線を吸収する物質と吸収されない物質があるため、被写体の像が形成される。

【0004】

Ｘ線画像検出器の課題は、検出されたＸ線から最大限の情報を抽出し、物体や被写体が密度、組成、構造などの観点から描写された画像の入力を提供することである。

【0005】

典型的な医療用Ｘ線イメージングシステムでは、Ｘ線はＸ線管によって生成される。一般的な医療用Ｘ線管のエネルギースペクトルは幅広く、０から１６０ｋｅＶまでの範囲である。したがって、検出器は通常、エネルギーが変化するＸ線を検出する。図１を参照して、全体的なＸ線イメージングシステムの概要を簡単に説明することは有用であろう。この例示的であるが非限定的な例では、Ｘ線イメージングシステム１００は、基本的に、Ｘ線源１０、Ｘ線検出器システム２０、及び関連する画像処理システム又は装置３０を備える。Ｘ線検出器システム２０は、Ｘ線源１０からの放射線を記録するように構成され、Ｘ線源１０は、任意のＸ線光学系によって任意に集束され、物体、対象、又はその一部を通過する。Ｘ線検出器システム２０は、Ｘ線検出器システム２０に少なくとも部分的に統合された適切なアナログ及び読み出し電子機器を介して画像処理システム３０に接続可能であり、画像処理システム３０によってＸ線検出器による画像処理及び／又は画像再構成を可能にする。

【0006】

例として、Ｘ線コンピュータ断層撮影（ＣＴ）システムは、Ｘ線源とＸ線検出器とを含み、Ｘ線源とＸ線検出器は、被写体又は物体の投影画像が、少なくとも１８０度。これは、最も一般的には、対象又は物体の周りを回転できる支持体に線源と検出器を取り付けることによって実現される。さまざまなビュー角度のさまざまな検出器要素に記録された投影を含む画像は、サイノグラム（sinogram）と呼ばれる。以下では、異なるビュー角度に対して異なる検出素子に記録された投影の集合は、たとえ検出器が二次元であってもサイノグラムと呼ばれ、サイノグラムは三次元画像にされる。

【0007】

Ｘ線イメージングのさらなる発展として、Ｘ線透過率を複数の異なるエネルギーレベルについて測定するエネルギー分解Ｘ線イメージング（スペクトルＸ線イメージングとも呼ばれる）がある。これは、線源を２つの異なる発光スペクトルの間で高速に切り替えるか、異なるＸ線スペクトルを発する２つ以上のＸ線源を使用するか、または入射放射線を２つ以上のエネルギーレベルで測定するエネルギー弁別検出器を使用することにより実現できる。このような検出器の例としては、マルチビン光子計数型検出器（multi-bin photon-counting detector）があり、記録された各光子が電流パルスを生成し、これを一連の閾値と比較することにより、多数のエネルギービンのそれぞれにおいて入射した光子の数が計数される。

【0008】

スペクトルＸ線投影測定により、各エネルギーレベルの投影画像が得られる。これらの投影画像の加重和は、「Tapiovaara and Wagner, "SNR and DQE analysis of broad spectrum X-ray imaging", Phys. Med. Biol. 30, 519」に記載されている。

【0009】

エネルギー分解Ｘ線イメージングによって可能になるもう１つの手法は、基底物質分解（basis material decomposition）である。この手法は、人間の組織などの原子番号の小さい元素から構成されるすべての物質が線形減衰係数を持ちμ（Ｅ）、そのエネルギー依存性を２つ（又はそれ以上）の基底関数の線形結合として近似的に表現できることを利用した手法である。
μ（Ｅ）＝ａ_１ｆ_１（Ｅ）＋ａ_２ｆ_２（Ｅ）
ここで、ｆｉは基底関数、ａｉは対応する基底係数である。イメージングに使用されるエネルギー範囲にｋ吸収端が存在するのに十分な高い原子番号を持つ１つ以上の元素がイメージングされたボリュームに存在する場合、１つの基底関数をそのような元素ごとに追加する必要がある。医用画像診断の分野では、そのようなｋ吸収端の元素は、典型的には、造影剤として使用される物質であるヨウ素またはガドリニウムであり得る。

【0010】

基底物質の分解は、Alvarez and Macovski, "Energy-selective reconstructions in X-ray computerised tomography", Phys.Med.Biol.21,733に記載されている。基底物質分解では、各基底係数の積分は、線源から検出器要素への各投影光線ｌの測定データから推測される。

【0011】

【数1】

ある実装では、これはまず、各エネルギービンにおけるカウント（計数）の予想記録数を、以下のＡｉの関数として表現することで達成される。

【0012】

【数2】

ここで、λｉはエネルギービンｉ内のカウント数の期待値，Ｅはエネルギーである。Ｓｉは応答関数で、撮像された物体に入射するスペクトル形状、検出器の量子効率、エネルギービンｉのエネルギーＥを有するＸ線に対する感度に依存する。「エネルギービン」という用語は光子計数検出器に最もよく使われるが、この式は多層検出器やｋＶｐスイッチング源など他のエネルギー分解Ｘ線システムでも使用できる。

【0013】

各ビンのカウント数がポアソン分布確率変数であるという仮定の下で、最尤法を使用してＡｉを推定できる（参照：Roessl and Proksa, K-edge imaging in x-ray computed tomography using multi-bin photon counting detectors, Phys. Med. Biol. 52 (2007), 4679-4696）。

【0014】

【数3】

ここで、ｍｉは，エネルギービンで測定されたカウント数であり、Ｍｂは、エネルギービンの数である。

【0015】

結果として得られる各投影線の推定基底係数線積分（estimated basis coefficient line integral）

【0016】

【数4】

が画像行列に配置されると、結果は各基底ｉの基底画像とも呼ばれる物質固有の投影画像になる。この基底画像は、直接表示するか（投影Ｘ線イメージングなど）、再構成アルゴリズムへの入力として取得して、オブジェクト内部の基底係数ａｉのマップを形成できる（例えばＣＴにおける）。いずれにしても、基底分解の結果は、基底係数の線積分や基底係数自体など、１つ又は複数の基底画像表現（basis image representations）と見なすことができる。

【0017】

ただし、この手法のよく知られた制限は、推定された線積分の分散（variance of the estimated line integrals）が通常、基底分解で使用される基底の数とともに増加することである。とりわけ、これは改善された組織定量化と増加した画像ノイズの間の不所望なトレードオフをもたらす。

【0018】

さらに、３つ以上の基底関数を使用した正確な基底分解は、実際には実行が難しく、アーティファクト、偏り（バイアス）、又は過剰なノイズ（artifacts, bias or excessive noise）が発生する可能性がある。このような基底分解では、正確な結果を得るために、過度なキャリブレーション測定（extensive calibration measurements）とデータの前処理も必要になる場合がある。

【0019】

また、Ｘ線画像再構成に関する改善に対する一般的な要求もある。

【先行技術文献】

【特許文献】

【0020】

【文献】米国特許第９，０３６，８８５号

【発明の概要】

【0021】

一般に、Ｘ線イメージング用途の画像再構成に関して改善を提供することが望ましい。

【0022】

Ｘ線イメージングのための基底分解に関して改善を提供することが特定の目的である。

【0023】

画像再構成の改善された方法を提供することが特定の目的である。

【0024】

画像再構成のための改良されたシステムを提供することも目的である。

【0025】

別の目的は、改良されたＸ線イメージングシステム全体を提供することである。

【0026】

対応するコンピュータ・プログラム及びコンピュータ・プログラム製品（コンピュータ・プログラム・プロダクト）を提供することがさらなる目的である。

【0027】

これら及び他の目的は、提案された技術の１つ又は複数の実施形態によって達成され得る。

【0028】

第１の態様によれば、エネルギー分解Ｘ線データに基づく画像再構成の方法が提供され、
エネルギー分解Ｘ線データの少なくとも１つの表現を収集し、
エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を実行して、少なくとも２つの元の基底画像表現セット（at least two original basis image representation sets：少なくとも２つのオリジナルの基底画像表現のセット）を生成し、
元の基底画像表現セットの少なくとも２つから少なくとも２つの基底画像表現を取得又は選択し、
取得又は選択された基底画像表現を機械学習に基づいてデータ処理して、出力画像データの少なくとも１つの表現を生成する。

【0029】

第２の態様によれば、エネルギー分解Ｘ線データに基づく画像再構成のためのシステムが提供される。システムは、エネルギー分解Ｘ線データの少なくとも１つの表現を取得するように構成されている。システムは、エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を行い、少なくとも２セットの基底画像表現を生成するように構成される。システムは、少なくとも２つの基底画像表現のセットから少なくとも２つの基底画像表現を取得又は選択するようにさらに構成される。システムはまた、取得又は選択された基底画像表現を機械学習に基づいてデータ処理して、出力画像データの少なくとも１つの表現を生成するように構成される。

【0030】

第３の態様によれば、エネルギー分解Ｘ線データに基づく画像再構成のためのシステムが提供される。システムは、
エネルギー分解Ｘ線データの少なくとも１つの表現を受信するための入力と、
エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を実行して、基底画像表現の少なくとも２つのセットを生成するように構成された基底物質分解システムと、
入力－出力対の訓練データに基づいて決定された構成（configuration determined based on input-output-paired training data）を有し、基底画像表現のセットの少なくとも２つから少なくとも２つの基底画像表現を受信し、受信した基底画像表現を機械学習に基づいてデータ処理して、出力画像データの少なくとも１つの表現を生成するように構成された訓練された機械学習システムと、
出力画像データの少なくとも１つの表現を出力する出力と、
を備える。

【0031】

第４の態様によれば、第２の態様又は第３の態様による画像再構成のためのシステムを備えるＸ線イメージングシステムが提供される。

【0032】

第５の態様によれば、プロセッサによって実行されると、プロセッサに複数のステップを行わせる命令を含むコンピュータプログラムを格納した非一時的コンピュータ可読媒体を含むコンピュータプログラム製品が提供される。前記複数のステップは、
エネルギー分解Ｘ線データの少なくとも１つの表現を収集し、
エネルギー分解Ｘ線データの前記少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を実行して、少なくとも２つの元の基底画像表現セットを生成し、
前記元の基底画像表現セットの少なくとも２つから少なくとも２つの基底画像表現を抽出し、
前記抽出された基底画像表現を機械学習に基づいてデータ処理して、出力画像データの少なくとも１つの表現を生成することを含む。

【0033】

提案された技術は、機械学習ベースのデータ処理への入力として、２つ以上の基底物質分解から得られた基底画像表現を使用することにより、大幅に改善された画像再構成を提供する。提案された技術によって提供される技術的解決策は、例として、出力画像及び対応する画像表現における偏り（バイアス）、アーティファクト、及び／又はノイズを低減することができる。

【0034】

本発明者らは、Ｘ線画像再構成を実行する際に、機械学習プロセス及び／又はシステムの効率、操作及び／又は結果を大幅に改善及び容易にすることが可能であることを認識した。生の画像データに対して機械学習を使用するという優勢な傾向に従う代わりに、本発明者らは、機械学習プロセス及び／又は機械学習システムの入力として、２つ以上の４基底物質分解から得られる基底画像表現を利用する場合に、改善された画像結果が達成できることを認識した。システム。たとえば、改善には、偏り、アーティファクト及び／又はノイズの低減、又は解像度の向上が含まれる場合がある。もう１つの利点は、提案された技術が、信頼性と品質の高い画像再構成を取得するために必要な訓練データの量を削減することである。本発明の実施形態は、そのさらなる目的及び利点とともに、添付の図面とともに以下の説明を参照することによって最もよく理解され得る。

【図面の簡単な説明】

【0035】

【図1】Ｘ線イメージングシステム全体の一例を示す概略図である。

【図2】Ｘ線イメージングシステムの別の例を示す概略図である。

【図3】Ｘ線イメージングシステムの例示的な例としてのＣＴシステムの概略ブロック図である。

【図4】Ｘ線イメージングシステムの関連部分の別の例を示す概略図である。

【図5】従来技術による光子計数回路及び／又はデバイスの概略図である。

【図6】一実施形態によるエネルギー分解Ｘ線測定データに基づく画像再構成の方法の一例を示す概略フローチャートである。

【図7】一実施形態による機械学習ベースの画像再構成の補足表現の一例を示す概略図である。

【図8】一実施形態によるＸ線画像再構成システムの一例を示す概略ブロック図である。

【図9】一実施形態による、提案されたディープラーニング偏り除去（debiase）のためのスキーム及び／又はアーキテクチャの例を示す概略図である。

【図10】残差ブロック（ＲｅｓＢｌｏｃｋｓ：residual blocks）の一例を示す模式図である。

【図11】二物質分解からの偏った結果の後の一連の残差ブロック（ＲｅｓＢｌｏｃｋｓ）を含む偏り除去ディープラーニングスキーム及び／又はアーキテクチャの例を示す概略図である。

【図12】シミュレートされたデータセットのサンプルの一例を示す概略図であり、シェップローガンファントムの異なる部分が異なる物質を模倣している。

【図13】シェップローガンによってシミュレートされたファントムに対して提案されたディープラーニングベースの偏り緩和の定性的な結果の例を示す概略図である。

【図14】機械学習システムが物質分解サイノグラムを再構成画像にマッピングする、提案された発明の例示的な実施形態を示す概略図である。

【図15】機械学習システムへの入力として使用される物質分解画像をワンステップ画像再構成システムが提供する、提案された発明の例示的な実施形態を示す概略図である。

【図16】提案された発明の例示的な実施形態を示す概略図であり、基底物質分解システムとそれに続く画像再構成システムを含む２ステップ画像再構成方法が、機械学習への入力として使用される物質分解画像を提供する。システム。

【図17】一実施形態によるコンピュータ実装の一例を示す概略図である。

【発明を実施するための形態】

【0036】

理解を深めるために、Ｘ線イメージングシステム全体の非限定的な例の紹介的な説明を続けることが有用な場合がある。

【0037】

図２は、Ｘ線を放射するＸ線源１０と、Ｘ線が物体を通過した後のＸ線を検出するＸ線検出器を備えたＸ線検出器システム２０と、Ｘ線検出器からの生の電気信号を処理してデジタル化するアナログ処理回路２５とを備えるＸ線イメージングシステム１００の一例を示す模式図である。補正の適用、一時的な保存、またはフィルタリングなど、測定データに対してさらなる処理操作を行うことができるデジタル処理回路４０と、処理されたデータを保存し、さらなる後処理および／または画像再構成を行うことができるコンピュータ５０と、を備える。本発明によれば、アナログ処理回路２５の全て又は一部は、Ｘ線検出器システム２０に実装されてもよい。

【0038】

Ｘ線検出器全体は、Ｘ線検出器システム２０、又は関連するアナログ処理回路２５と組み合わされたＸ線検出器システム２０と見なすことができる。

【0039】

デジタル処理回路４０及び／又はコンピュータ５０を含むデジタル部分は、Ｘ線検出器からの画像データに基づいて画像再構成を実行する画像処理システム３０とみなすことができる。したがって、画像処理システム３０は、コンピュータ５０、または代替的にデジタル処理回路４０とコンピュータ５０の複合システム、またはデジタル処理回路が画像処理および／または再構成のためにもさらに特化されている場合には、デジタル処理回路４０単体として見なすことができる。

【0040】

一般的に使用されるＸ線イメージングシステムの一例は、Ｘ線によるコンピュータ断層撮影（ＣＴ）システムであり、Ｘ線のファンビームまたはコーンビームを生成するＸ線管と、患者または物体を透過したＸ線の一部（fraction）を測定するＸ線検出器の対向するアレイとを含むことができる。Ｘ線管と検出器アレイは、撮影対象物の周りを回転するガントリに取り付けられている。

【0041】

図３は、Ｘ線イメージング・システムの例示的な例としてのＣＴシステムの概略ブロック図である。ＣＴシステムは、オペレータコンソール６０を介してオペレータからコマンド及び走査パラメータを受信するコンピュータ５０を含み、オペレータコンソール６０は、ディスプレイ及び何らかの形のオペレータインターフェース、例えば、キーボード及びマウスを有することができる。オペレータが供給したコマンドとパラメータは、次にコンピュータ５０によって利用され、Ｘ線コントローラ４１、ガントリコントローラ４２、及びテーブルコントローラ４３に制御信号を提供する。具体的には、Ｘ線コントローラ４１は、電力及びタイミング信号を提供する。ガントリコントローラ４２は、Ｘ線源１０及びＸ線源１０を備えるガントリ１１の回転速度及び位置を制御する。例として、Ｘ線検出器は光子計数Ｘ線検出器であってもよい。テーブルコントローラ４３は、患者テーブル１２の位置及び患者の走査範囲を制御及び決定する。検出器コントローラ４４は、検出器２０からのデータを制御及び／又は受信するように構成されている。

【0042】

一実施形態では、コンピュータ５０は、Ｘ線検出器から出力された画像データの後処理及び画像再構成も実行する。したがって、コンピュータは、図１及び図２に示される画像処理システム３０に対応する。関連付けられたディスプレイにより、オペレータはコンピュータから再構成された画像やその他のデータを観察できる。

【0043】

ガントリ１１に配置されたＸ線源１０はＸ線を放出する。例えば光子計数検出器（フォトンカウンティング検出器）の形態のＸ線検出器２０は、患者を通過した後のＸ線を検出する。Ｘ線検出器２０は、例えば、センサ又は検出器要素とも呼ばれる複数のピクセル、及び検出器モジュールに配置されたＡＳＩＣなどの関連する処理回路によって形成され得る。アナログ処理部分の一部は、ピクセルで実装されてもよいが、残りの処理部分は、例えばＡＳＩＣで実装される。一実施形態では、処理回路（ＡＳＩＣ）は、ピクセルからのアナログ信号をデジタル化する。処理回路（ＡＳＩＣ）は、デジタル処理部を備えてもよく、デジタル処理部は、測定データに対して、補正の適用、一時的な保存、及び／又はフィルタリングなどのさらなる処理操作を実行することができる。Ｘ線投影データを取得するための走査中、ガントリ及びガントリに取り付けられた構成要素はアイソセンタ１３を中心に回転する。

【0044】

現代のＸ線検出器は通常、入射したＸ線を電子に変換する必要がある。これは通常、光電効果またはコンプトン相互作用によって行われ、生成された電子はエネルギーが失われるまで通常二次可視光を生成し、この光が感光材料によって検出される。この場合、Ｘ線によって生成された電子は、適用された電場によって収集される電子－正孔ペアの観点から電荷を生成する。

【0045】

多数のＸ線から積分された信号を提供するという意味で、エネルギー積分モードで動作する検出器がある。出力信号は、検出されたＸ線によって蓄積された総エネルギーに比例する。

【0046】

光子計数機能とエネルギー分解機能を備えたＸ線検出器は、医療用Ｘ線アプリケーションで一般的になりつつある。光子計数検出器は、原理的に各Ｘ線のエネルギーを測定することができ、対象物の組成に関する追加情報をもたらすという利点を有する。この情報を使用して、画質を向上させたり、放射線量を減らしたりすることができる。

【0047】

一般に、光子計数Ｘ線検出器は、検出器材料内の光子相互作用によって生成される電気パルスの高さを一連のコンパレータ電圧と比較することによって、光子のエネルギーを決定する。これらのコンパレータ電圧は、エネルギーしきい値とも呼ばれる。一般に、比較器（コンパレータ）のアナログ電圧は、Ｄ／Ａ変換器（ＤＡＣ）によって設定される。ＤＡＣは、コントローラから送信されたデジタル設定を、光子パルスの高さを比較できるアナログ電圧に変換する。

【0048】

光子計数検出器は、測定時間中に検出器内で相互作用した光子の数をカウントする。一般に、新しい光子は、電気パルスの高さが少なくとも１つの比較器の比較器電圧を超えるという事実によって識別される。光子が識別されると、チャネルに関連付けられたデジタルカウンタをインクリメント（増分）することによってイベントが保存される。

【0049】

いくつかの異なるしきい値を使用すると、いわゆるエネルギー弁別光子計数検出器が得られ、検出された光子をさまざまなしきい値に対応するエネルギービンに分類できる。このタイプの光子計数検出器は、マルチビン検出器と呼ばれることもある。一般に、エネルギー情報により、新しい種類の画像を作成できるようになり、新しい情報が利用可能になり、従来の技術に固有の画像アーティファクトを取り除くことができる。言い換えれば、エネルギー弁別光子計数検出器の場合、パルス高は比較器のプログラム可能なしきい値（Ｔ１－ＴＮ）の数と比較され、エネルギーに比例するパルス高に従って分類される。つまり、複数の比較器からなる光子計数検出器を、ここではマルチビン型光子計数検出器と呼ぶことができる。マルチビン光子計数検出器の場合、光子計数は、典型的にはエネルギー閾値ごとに１つのカウンタの集合に格納される。例えば、カウンタは、光子パルスが超えた最も高いエネルギー閾値に対応するように割り当てられることができる。別の例では、カウンタは、光子パルスが各エネルギー閾値を横切った回数を記録している。

【0050】

一例として、エッジオン（edge-on）は、光子計数検出器用の特別な非限定的な設計であり、Ｘ線検出素子又はピクセルなどのＸ線センサが入射Ｘ線に対してエッジオンに向けられる。

【0051】

例えば、そのような光子計数検出器は、少なくとも２つの方向にピクセルを有することができ、エッジオン光子計数検出器の方向の１つは、Ｘ線の方向の成分を有する。このようなエッジオン光子計数検出器は、入射Ｘ線の方向にピクセルの２つ以上の深度セグメントを有する、深度セグメント化された光子計数検出器と呼ばれることがある。

【0052】

あるいは、ピクセルは、入射Ｘ線の方向に実質的に直交する方向にアレイ（深さ分割されていない）として配置することができ、ピクセルのそれぞれは、入射Ｘ線に対して真横に向けることができる。言い換えれば、光子計数型検出器は、入射Ｘ線に対してエッジオンに配置されたまま、非深度分割型であってもよい。

【0053】

エッジオン光子計数検出器をエッジオンに配置することにより、吸収効率を高めることができる。この場合、吸収深さは任意の長さに選択でき、エッジオン光子計数検出器は、非常に高い電圧にならなくても完全に使い果たす（fully deplete）ことができる。

【0054】

直接半導体検出器をを介してＸ線光子を検出する従来の仕組みは、基本的に次のように動作する。検出器材料におけるＸ線相互作用のエネルギーは、半導体検出器内部で電子－正孔対に変換され、電子－正孔対の数は一般に光子エネルギーに比例する。電子と正孔は、検出器の電極と裏面に向かってドリフトする（またはその逆）。このドリフトの間、電子と正孔は電極に電流を誘起し、その電流は測定されることができる。

【0055】

図４に示すように、信号（複数可）は、Ｘ線検出器の検出器要素子２１から並列処理回路（例えば、ＡＳＩＣ）２５の入力２７に送られる。特定用途向け集積回路（ＡＳＩＣ）という用語は、特定のアプリケーション用に使用及び構成された一般的な回路として広く解釈される。ＡＳＩＣは、各Ｘ線から生成された電荷を処理し、それをデジタルデータに変換し、それを用いて光子数及び／又は推定エネルギーなどの測定データを得ることができる。ＡＳＩＣは、デジタルデータ処理回路に接続するように構成されているので、デジタルデータは、さらなるデジタルデータ処理４０及び／又は１つ以上のメモリ４５に送ることができ、最終的にデータは、再構成画像を生成するための画像処理３０．５０の入力となるであろう。

【0056】

１つのＸ線イベントから発生する電子と正孔の数はＸ線光子のエネルギーに比例するため、１つの誘導電流パルスの総電荷はこのエネルギーに比例する。ＡＳＩＣでのフィルタリング工程の後、パルス振幅は電流パルスの総電荷に比例するため、Ｘ線エネルギーに比例する。次に、パルス振幅は、その値を１つ又は複数の比較器（ＣＯＭＰ：comparator）で１つ又は複数のしきい値（ＴＨＲ）と比較することによって測定できる。また、パルスがしきい値よりも大きい場合の数を記録できるカウンタが導入される。このようにして、特定の時間枠内で検出されたそれぞれの閾値（ＴＨＲ）に対応するエネルギーを超えるエネルギーを有するＸ線光子の数を計数及び／又は記録することが可能である。

【0057】

ＡＳＩＣは通常、クロックサイクルごとに１回アナログ光子パルスをサンプリングし、比較器の出力を記録する。比較器（しきい値）は、アナログ信号がコンパレータ電圧より上か下かに応じて、１又は０を出力する。各サンプルで利用可能な情報は、たとえば、比較器がトリガーされた（光子パルスがしきい値よりも高かった）かどうかを表す各比較器の１又は０である。

【0058】

光子計数検出器では、通常、新しい光子が記録（register）されたかどうかを判断し、光子をカウンタに記録する光子計数ロジックがある。マルチビン光子計数検出器の場合、通常、複数のカウンタがあり、例えば、各比較器に１つずつあり、光子計数は、光子エネルギーの推定値に従ってカウンタに記録される。このロジックは、いくつかの異なる方法で実装することができる。光子計数ロジックの最も一般的なカテゴリーは、いわゆる非パラライザブル計数モード（non-paralyzable counting modes）と、パラライザブル計数モード（paralyzable counting modes）との２つである。その他の光子計数ロジックとしては、例えば、ローカルマキシマ検出（local maxima detection）があり、電圧パルスで検出された極大値をカウントし、場合によってはパルス高も記録する極大値検出が含まれる。

【0059】

光子計数型検出器には、高い空間分解能、低い電子ノイズ、エネルギー分解能、物質分離能力（スペクトルイメージング能力）など、多くの利点があるが、これらに限定されるわけではない。しかし、エネルギー積分型検出器には、高い計数率耐性（high count-rate tolerance；計数率の許容度が高い）という利点がある。計数率の許容誤差は、光子の全エネルギーを測定するため、現在検出器に記録されている光子の量に関係なく、光子を１つ追加すれば必ず出力信号が（合理的な範囲内で）増加するという事実／認識に由来する。この利点は、エネルギー積分型検出器が今日の医療用ＣＴの標準である主な理由の１つである。

【0060】

提案された測定方法のより良い理解のために、簡単なシステムの概要及び／又は技術的問題の分析から始めることが有用であろう。この目的のために、先行技術による光子計数回路及び／又はデバイスの概略図を提供する図５を参照する。

【0061】

光子が半導体材料中で相互作用すると、電子－正孔対のクラウド（cloud：雲）が生成される。検出器材料に電界をかけると、電荷キャリアは検出器材料に取り付けられた電極に集められる。信号は、検出器から並列処理回路（ＡＳＩＣなど）の入力に送られる。特定用途向け集積回路（ＡＳＩＣ）という用語は、特定の用途に使用され構成された一般的な回路として広く解釈されることが理解されるべきである。ＡＳＩＣは、各Ｘ線から発生した電荷を処理し、それをデジタルデータに変換し、それを用いて光子数及び／又は推定エネルギーなどの測定データを得ることができる。一例として、ＡＳＩＣは、検出器材料に光子によって堆積されたエネルギー量に比例する最大高さの電圧パルスが生成されるように電荷を処理することができる。

【0062】

ＡＳＩＣは、各比較器３０２が電圧パルスの大きさを基準電圧と比較する比較器３０２のセットを含むことができる。比較器の出力は、比較された２つの電圧のどちらが大きいかに応じて、ゼロ又は１（０／１）になる。ここで、電圧パルスが基準電圧よりも高い場合、比較器出力は（１）であり、基準電圧が電圧パルスよりも高い場合（０）であると仮定する。デジタル／アナログ変換器（ＤＡＣ）３０１を使用して、ユーザー又は制御プログラムによって供給され得るデジタル設定を、比較器３０２によって使用され得る基準電圧に変換することができる。電圧パルスの高さが特定の比較器の基準電圧を超える場合、その比較器がトリガされたと呼ぶことにする。各比較器は、一般に、光子計数論理に従って比較器出力に基づいて増分されるデジタルカウンタ３０３と関連付けられる。

【0063】

一般に、基底物質分解は、人体組織のような原子番号の小さい元素から構成されるすべての物質が線形減衰係数μ（Ｅ）を持つという事実を利用する。そのエネルギー依存性は、２つ（以上）の基底関数の線形結合として近似的に表現できることを利用する。
μ（Ｅ）＝ａ_１ｆ_１（Ｅ）＋ａ_２ｆ_２（Ｅ）
ここで、ｆｉは基底関数、ａｉは対応する基底係数である。イメージングに使用されるエネルギー範囲にｋ吸収端が存在するのに十分な高い原子番号を持つ１つ以上の元素がイメージングされたボリュームに存在する場合、１つの基底関数をそのような元素ごとに追加する必要がある。医用画像診断の分野では、そのようなｋ吸収端の元素は、典型的には、造影剤として使用される物質であるヨウ素またはガドリニウムであり得る。

【0064】

前述のように、基底係数ａｉのそれぞれの線積分Ａｉは、線源から検出素子への各投影光線ｌにおける測定データから推測される。線積分Ａｉは次のように表すことができる。

【0065】

【数5】

一実施例において、基底物質の分解は、最初に、各エネルギービンにおける期待される記録されるカウント数を関数Ａｉとして表すことによって達成される。このような関数は次のような形式を取る。

【0066】

【数6】

ここでλｉはエネルギービンｉ内のカウント数の期待値、Ｅはエネルギーである。Ｓｉは応答関数で、撮像された物体に入射するスペクトル形状、検出器の量子効率、エネルギービンｉのエネルギーＥを有するＸ線に対する感度に依存する。「エネルギービン」という用語は光子計数検出器に最もよく使われるが、この式は多層検出器やｋＶｐスイッチング源など他のエネルギー分解Ｘ線システムにも適用できる。

【0067】

次に、最尤推定法（maximum likelihood method）を使用して、各ビンのカウント数がポアソン分布確率変数となるという仮定の下でＡｉを推定することができる。これは、負の対数尤度関数を最小化する（minimizing the negative log-likelihood function）ことによって達成される（参照：Roessl and Proksa, K-edge imaging in x-ray computed tomography using multi-bin photon counting detectors, Phys. Med. Biol. 52 (2007), 4679-4696）。

【0068】

【数7】

ここでｍｉはエネルギービンｉ内の測定カウント数であり、Ｍｂはエネルギービンの数である。

【0069】

線積分Ａから、基底係数を取得する断層撮影再構成ａｉが実行される場合がある。この手順工程は、別個のトモグラフィー再構成と見なすことができ、あるいは全体的な基底分解の一部とみなすこともできる。

【0070】

前述のように、各投影線の結果として得られる。

【0071】

【数8】

が画像マトリックスに配置されると、その結果、各基底ｉの基底画像とも呼ばれる物質固有の投影画像となる。この基底画像は、直接表示するか（投影Ｘ線イメージングなど）、再構成アルゴリズムへの入力として取得して、物体内部の基底係数ａｉのマップを形成できる（ＣＴなど）。いずれにしても、基底分解の結果は、基底係数の線積分や基底係数自体など、１つ又は複数の基底画像表現と見なすことができる。

【0072】

提案する技術は、２つ以上の基底物質分解から得られる基底画像表現を機械学習ベースのデータ処理への入力として使用することにより、画像再構成を大幅に改善することを目的とする。提案された技術によって提供される技術的解決策は、一例として、出力画像及び対応する画像表現における偏り、アーチファクト及び／又はノイズを低減することができる。

【0073】

図６は、実施形態によるエネルギー分解Ｘ線測定データに基づく画像再構成方法の一例を示す模式的なフロー図である。
基本的に、この方法には以下の工程が含まれる。
Ｓ１：エネルギー分解されたＸ線データの少なくとも１つの表現を収集する。
Ｓ２：エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて、少なくとも２つの基底物質分解を実行して、少なくとも２つの元の基底画像表現セットを生成する。
Ｓ３：元の基底画像表現セットの少なくとも２つから少なくとも２つの基底画像表現を取得又は選択する。
Ｓ４：機械学習に基づくデータ処理によって取得又は選択された基底画像表現を処理し、出力画像データの少なくとも１つの表現を生成する。
一例として、機械学習に基づくデータ処理は、取得又は選択された基底画像表現を、出力画像データの少なくとも１つの表現に対応付けるように構成され得る。

【0074】

特定の例では、機械学習に基づくデータ処理は、訓練されたニューラル・ネットワーク、サポートベクターマシン、又は決定木ベースの手段に基づく実装によって、取得又は選択された基底画像表現を出力画像データの少なくとも１つの表現にマッピングするように構成されてもよい。

【0075】

例えば、訓練（トレーニング）されたニューラル・ネットワーク、サポートベクターマシン、又は決定木に基づく実装は、測定された患者の画像データ又はファントム画像データの少なくとも１つに対して、又は数値ファントムに基づくシミュレートされた画像データで訓練することができる。

【0076】

任意で、ニューラル・ネットワーク、サポートベクターマシン、又は決定木ベースの実装は、第１のタイプの画像で訓練され、転移学習を使用して第２のタイプの画像で使用するように適合される場合がある。

【0077】

一例として、ニューラル・ネットワーク、サポートベクターマシン、又は決定木に基づく実装は、非エネルギー分解ＣＴ画像上にあり、転移学習を使用してエネルギー分解光子計数ＣＴ画像上で使用するように適合される場合がある。

【0078】

典型的には、エネルギー分解Ｘ線データの少なくとも１つの表現は、Ｘ線検出器又は中間記憶装置から生じさせることができ、エネルギー分解Ｘ線画像データ及び／又はエネルギー分解Ｘ線測定データを含み得る。

【0079】

例えば、Ｘ線検出器は、マルチビン光子計数Ｘ線検出器とも呼ばれるエネルギー分解光子計数検出器であってもよい。

【0080】

あるいは、エネルギー分解Ｘ線データは、マルチＸ線管収集、低速又は高速ｋＶスイッチング収集、多層検出器から、又は分割フィルタ収集によって取得することができる。

【0081】

図７は、一実施形態による機械学習ベースの画像再構成の補足表現の一例を示す概略図である。基本的に、エネルギー分解されたＸ線データは、基底物質分解への入力として提供され、少なくとも２セットの基底画像表現を生成し、その少なくともサブセットが機械学習プロセス及び／又はシステムに転送される。機械学習プロセス及び／又はシステムは、適切に訓練されると、基底画像表現を出力画像データに（通常、推論と呼ばれるプロセスによって）マッピングする。

【0082】

特定の実施例では、出力画像データの前記少なくとも１つの表現は、基底画像表現の出力セットと呼ばれる基底画像表現のセットを含み、出力セット内の基底画像表現の数は、前記少なくとも２つの元の基底画像表現セットのいずれかにおける基底画像表現の最大数より大きい。

【0083】

一例として、基底画像表現の出力セットは、多数の基底関数を有してよく、元の基底画像表現セットの各々は、それぞれの数の基底関数を有してよく、基底画像表現の出力セットに対応する基底関数は、前記少なくとも２つの元の基底画像表現セットに対応する基底関数の和の部分集合であるか、又はそれに等しくてもよい。

【0084】

例えば、基底画像表現の出力セットは多数の基底関数を有し、元の基底画像表現セットの各々はそれぞれの基底関数の数を有し、前記取得された又は選択された基底画像表現の数は、基底画像表現の前記出力セットに対応する基底関数の数と同じか又はそれより大きい。

【0085】

特定の例では、元の基底画像表現セットは２基底分解（two-basis decompositions）に基づいており、基底画像表現の出力セットは３基底分解（three-basis decompositions）に基づいている。

【0086】

このような場合、機械学習に基づくデータ処理は、例えばニューラルネットワークによって実行されてもよく、ニューラルネットワークは、入力データとしての複数の２基底画像表現セットと出力データとしての３基底画像表現セットとを用いて訓練されてもよい。

【0087】

一例として、２基底画像表現セット及び３基底画像表現セットは、測定された患者画像データ又はファントム画像データから、又は数値ファントムに基づくシミュレートされた画像データから生成され得る。

【0088】

オプションとして、取得又は選択された基底画像表現を処理するステップは、以下を含むことができる：
－取得又は選択された基底画像表現を、非基底画像データ（non-basis image data）の少なくとも１つの表現と組み合わせるステップ
－機械学習に基づくデータ処理によって、取得又は選択された基底画像表現を少なくとも１つの非基底画像データの表現と一緒に処理するステップ。

【0089】

例えば、非基底画像データは、ビン計数サイノグラムデータ、対数正規化ビン計数サイノグラムデータ、サムオーバービンサイノグラムまたは対数正規化データ、再構成ＣＴ画像データ、およびセグメンテーションマスクデータ（bin counts sinogram data, log-normalized bin counts sinogram data, sum-over-bins sinogram or log-normalized data, reconstructed CT image data, and segmentation mask data）のうちの少なくとも１つを含むことができる。

【0090】

オプションとして、取得又は選択された基本画像表現を処理するステップは、以下を含むことができる。
－前記取得又は選択された基底画像表現を事前情報の少なくとも１つの表現と組み合わせるステップ
－機械学習に基づく前記データ処理によって、事前情報の前記少なくとも１つの表現を用いて、前記取得又は選択された基底画像表現を処理するステップ。

【0091】

事前情報とは、一般に、例えば、物理法則及び原理並びに以前の測定、画像及び／又はシミュレーションに基づき、既知であると想定される利用可能な情報の別個のセットとして定義される。

【0092】

オプションとして、取得又は選択された基底画像表現を処理するステップは、以下を含むことができる：
－取得又は選択された基底画像表現を、取得システム内の物理プロセスの少なくとも１つの計算表現と組み合わせるステップ、
－取得又は選択された基底画像表現を、機械学習に基づくデータ処理によって、物理プロセスの前記少なくとも１つの計算表現で処理するステップ。

【0093】

オプションとして、取得又は選択された基底画像表現を処理するステップは、以下を含むことができる：
－取得又は選択された基底画像表現を、測定データの統計モデルの少なくとも１つの計算表現と組み合わせるステップ、
－機械学習に基づくデータ処理によって、統計モデルの前記少なくとも１つの計算表現を用いて、取得又は選択された基底画像表現を処理するステップ。

【0094】

特定の例では、出力画像データの前記少なくとも１つの表現は、基底画像表現及び／又は物質分解再構成画像の出力セットを含む。

【0095】

一例として、機械学習に基づくデータ処理は、訓練データに基づいて構成されたディープラーニング後処理を含むことができ、ディープラーニング後処理は、偏った物質分解サイノグラム（biased material-decomposed sinograms）を入力として受け取り、出力として偏りのない物質分解サイノグラム（de-biased material-decomposed sinograms）又は偏りのない物質分解再構成画像（debiased material-decomposition reconstructed imagesを提供するように構成されてもよい。

【0096】

例として、サイノグラムのセットを２基底分解の入力として提供して、２つ以上の基底画像表現を生成することができ、ディープラーニングニューラル・ネットワークを使用して、結果の基底画像表現を３物質分解に変換することができる。

【0097】

第２の態様によれば、エネルギー分解Ｘ線データに基づく画像再構成のためのシステムが提供される。システムは、エネルギー分解Ｘ線データの少なくとも１つの表現を取得するように構成されている。システムは、エネルギー分解Ｘ線データの少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を行い、少なくとも２セットの基底画像表現を生成するように構成される。システムは、少なくとも２つの基底画像表現のセットから少なくとも２つの基底画像表現を取得又は選択するようにさらに構成される。システムはまた、出力画像データの少なくとも１つの表現を生成するために、機械学習システムに基づくデータ処理によって、取得又は選択された基底画像表現を処理するように構成される。

【0098】

第３の態様によれば、エネルギー分解Ｘ線データに基づく画像再構成のためのシステムが提供される。システムは次の構成要素を備える。
－エネルギー分解Ｘ線データの少なくとも１つの表現を受信するための入力、
－エネルギー分解Ｘ線データの前記少なくとも１つの表現に基づいて少なくとも２つの基底物質分解を実行して、少なくとも２セットの基底画像表現を生成するように構成された基底物質分解システム、
－入出力ペア訓練データに基づいて決定された構成を有し、前記ベース画像表現のセットの少なくとも２つから少なくとも２つのベース画像表現を受信し、受信したベース画像表現を機械学習によって処理するように適合された、訓練済み機械学習システムと、出力画像データの少なくとも１つの表現を生成するためのベースのデータ処理。
－出力画像データの前記少なくとも１つの表現を出力するための出力。

【0099】

図８は実施形態に係るＸ線画像再構成システムの一例を示す模式的なブロック図である。図８のＸ線画像再構成システム３０は、基本的に、入力３１、基底物質分解システム３２、機械学習システム３３、及び出力３４から構成される。

【0100】

上述のように、機械学習システム３３は、例えば、ニューラル・ネットワーク、サポートベクターマシン、又は決定木ベース（a neural network, a support vector machine or a decision-tree based implementation）の実装に基づくことができる。

【0101】

例として、機械学習システム３３は、測定された患者画像データ又はファントム画像データの少なくとも１つに対して、又は数値ファントムに基づくシミュレートされた画像データから訓練され得る。

【0102】

特定の例では、機械学習システム３３は、第１のタイプの画像で訓練され、転移学習を使用して第２のタイプの画像で使用するように適合される。

【0103】

例えば、機械学習システム３３は、非エネルギー分解ＣＴ画像で訓練され、転移学習を使用してエネルギー分解光子計数ＣＴ画像での使用に適合され得る。

【0104】

特定の例では、機械学習システム３３は、取得又は選択された基底画像表現を出力画像データの前記少なくとも１つの表現にマッピングするように構成される。

【0105】

典型的には、エネルギー分解Ｘ線データの前記少なくとも１つの表現は、Ｘ線検出器又は中間記憶装置から生じ得、エネルギー分解Ｘ線画像データ及び／又はエネルギー分解Ｘ線測定データを含み得る。

【0106】

例として、Ｘ線検出器は、マルチビン光子計数Ｘ線検出器とも呼ばれるエネルギー分解光子計数検出器であってもよい。

【0107】

前述のように、出力画像データの前記少なくとも１つの表現は、基底画像表現の出力セットと呼ばれる基底画像表現のセットを含んでよく、出力セットにおける基底画像表現の数は、前記少なくとも２つの元の基底画像表現セットのいずれかにおける基底画像表現の最大数より多い。

【0108】

たとえば、基底画像表現の出力セットは、多数の基底関数を有し、各元の基底画像表現セットは、それぞれの数の基底関数を有し、基底画像表現の出力セットに対応する基底関数は、前記少なくとも２つの元の基底画像表現セットに対応する基底関数の和の部分集合であるか、又はそれに等しい。

【0109】

例えば、出力画像データの前記少なくとも１つの表現は、基底画像表現及び／又は物質分解再構成画像の出力セットを含むことができる。

【0110】

特定の例では、画像再構成システムは、少なくとも１つのプロセッサ及びメモリを備え、メモリは、少なくとも１つのプロセッサによって実行されると、少なくとも１つのプロセッサに画像再構成を実行させる命令を備える。

【0111】

例として、基底物質分解システムは、基底画像表現の少なくとも２つのセットを生成するように構成された処理回路を備えてもよく、訓練された機械学習システムは、出力画像データの少なくとも１つの表現を生成するように構成された処理回路を備えてもよい。

【0112】

より良く理解するために、次に、提案する技術を、２つ以上の基底分解から得られた基底画像表現を機械学習プロセス及び／又はシステムへの入力として用いる画像再構成に関して、非限定的な例を参照して説明する。

【0113】

Ｘ線イメージングの分野では、画像データの表現は、例えばサイノグラム、投影画像、又は再構成されたＣＴ画像を含むことがある。画像データのそのような表現は、異なるチャネルのデータが異なるエネルギー間隔で測定されたＸ線データに関連する複数のチャネルを含む場合、エネルギー分解され得る。

【0114】

エネルギー分解Ｘ線データの表現を入力とする物質分解のプロセスを通じて、基底画像表現セットを生成することができる。このようなセットは、多数の基底画像表現の集まりであり、各基底画像表現は、全Ｘ線減衰に対する特定の基底関数の寄与に関連している。このような一連の基底画像表現は、基底サイノグラムの集合、再構成された基底ＣＴ画像の集合、又は投影画像の集合であってよい。この文脈における「画像」は、例えば二次元画像、三次元画像、又は時間分解画像系列（time-resolved image series）を意味し得ることが理解されよう。

【0115】

例として、エネルギー分解Ｘ線データは、エネルギー分解Ｘ線画像データ及び／又はエネルギー分解Ｘ線測定データであってもよい。一例として、エネルギー分解されたＸ線データの少なくとも１つの表現は、Ｘ線検出器又は中間記憶装置から生じ得る。

【0116】

たとえば、エネルギー分解Ｘ線データの表現は、エネルギー・ビン・サイノグラムのコレクションを含むことができる。ここで、各エネルギーサイノグラムには、１つのエネルギービンで測定されたカウント数が含まれる。このエネルギービンサイノグラムの集まりを物質分解アルゴリズムへの入力とすることにより、基底サイノグラムのセットを生成することができる。このような基底サイノグラムは、例示的に、再構成アルゴリズムへの入力として取り込まれ、再構成された基底画像を生成することができる。

【0117】

２基底分解では、撮像された物体中の任意の物質の減衰が、２基底関数の線形結合として表現できるという近似に基づいて、２基底画像表現が生成される。３基底分解では、撮像された物体中の任意の物質の減衰が３基底画像上の線形結合として表現されるという近似に基づいて、３基底画像表現が生成される。同様に、４基底分解、５基底分解、および同様の多基底分解（four-basis decomposition, a five-basis decomposition and similar multi-basis decomposition ）を定義することができる。また、画像対象物中の全ての物質を、密度スケールファクターまで類似したエネルギー依存性を持つＸ線減衰係数を持つものとして近似し、１基底分解（one-basis decomposition）を行うことも可能である。

【0118】

２基底分解は、例えば、それぞれ水及びヨウ素の線形減弱係数によって与えられる基底関数に対応する、水サイノグラム及びヨウ素サイノグラムを含む基底サイノグラムのセットをもたらすことができる。あるいは、基底関数は、水とカルシウム、カルシウムとヨウ素、ポリ塩化ビニルとポリエチレンの減衰を表すことができる。２基底分解は、例えば、水サイノグラム、カルシウムサイノグラム、ヨウ素サイノグラムからなる基底サイノグラムのセットをもたらすことができる。あるいは、基底関数は、水、ヨウ素、ガドリニウムの減衰、またはポリ塩化ビニル、ポリエチレン、ヨウ素の減衰を表すことができる。

【0119】

発明者は、３つ以上の基底関数を用いた正確な基底分解は、実際には実行するのが難しく、アーティファクト、偏り、又は過度のノイズをもたらす可能性があることを認識した。このような基底分解では、正確な結果を得るために、大規模なキャリブレーション測定とデータの前処理も必要になる場合がある。一般に、多数の基底関数への基底分解は、少数の基底関数への分解よりも技術的に難しい場合がある。

【0120】

たとえば、２基底分解と比較して、画像の偏りやアーティファクトのレベルが低い３基底分解を行うのに十分正確なキャリブレーションを実行することは困難な場合がある。また、過度にノイズの多い基底画像を生成せずに、非常にノイズの多いデータで３基底分解を実行できる物質分解アルゴリズムを見つけるのは難しい場合がある。すなわち、クラメール・ラオの下限（Cramer-Rao lower bound）によって与えられる基底画像ノイズの理論的下限を達成することが困難で、２基底分解を行う場合にはこの下限を達成することが容易である場合がある。

【0121】

本発明者らはさらに、より多くの基底画像表現を生成するために必要な情報量は、それぞれがより少ない数の基底画像表現を有する複数の基底画像表現のセットから抽出することが可能であり得ることを理解している。例えば、水、カルシウム、及びヨウ素のサイノグラムへの３基底分解を生成するために必要な情報は、水－カルシウム分解、水－ヨウ素分解、及びカルシウム－ヨウ素分解の３つの２基底分解のセット（a set of three two-basis decompositions）から抽出することが可能であり得る。

【0122】

単一の３基底分解を正確に実行するよりも、複数の２基底分解を正確に実行する方が簡単な場合がある。開示された発明の特定の例において、この観察は、例えば、正確な３基底分解を実行するという問題を解決するために使用され得る。例として、エネルギー分解画像データを最初に使用して、水－カルシウム分解、水－ヨウ素分解、及びカルシウム－ヨウ素分解を実行することができる。次に、畳み込みニューラルネットワークなどの機械学習システムを用いて、得られた６つの（３×２）基底画像、またはそのサブセットを、水、カルシウム、およびヨウ素画像からなる3つの出力画像の集合に対応付けることができる。このような機械学習システム又はネットワークは、入力データとしてのいくつかの２基底画像表現セット及び出力データとしての３基底画像表現セットを用いて訓練することができ、前記２基底画像表現セット及び３基底画像表現セットは、測定した患者画像データ若しくはファントム画像データから、又は数値ファントムに基づくシミュレーション画像データから生成される。

【0123】

上述の方法を使用すると、エネルギー分解された画像データに対して直接実行される３基底分解と比較して、３基底画像表現セットの偏り、アーティファクト、又はノイズを大幅に削減できる。あるいは、より解像度の高い画像を得ることができる。

【0124】

元の基底画像に適用される機械学習アルゴリズムは、ニューラル・ネットワークの代わりに、サポートベクターマシン又は決定木ベースの方法などの別の機械学習方法を含むことができる。

【0125】

より一般的には、本発明者らは、Ｘ線画像再構成を行う際の機械学習プロセス及び／又はシステムの効率、動作及び／又は結果を実質的に改善及び促進することが可能であることに気付いた。生画像データに機械学習を使用する優勢な傾向に従う代わりに、本発明者らは、２つ以上の基底物質分解から得られる基底画像表現を機械学習プロセス及び／又はシステムの入力として利用する場合に、改善された画像結果を達成することができることに気がついた。例えば、改善は、偏り、アーティファクト及び／又はノイズの低減、又はより高い解像度を含むことができる。別の利点は、提案された技術が、高い信頼性と品質を有する画像再構成を得るために必要な学習データの量を減少させることである。

【0126】

元の基底画像表現を生成するために使用される基底物質分解は、例えば体積又は質量保存制約又は非負性制約などの事前情報を使用し及び／又は関与させてもよい。あるいは、そのような事前情報は、例えば、以前の検査からの画像又は全エネルギービンにおける集計カウントから再構成された画像などの事前画像表現の形態をとることができ、アルゴリズムは、この事前画像表現に対する分解された基底画像表現の逸脱にペナルティを科すことができる。別の選択肢は、例えば学習された辞書又は事前に訓練された畳み込みニューラルネットワークとして表される、学習画像のセットから学習された事前情報、又は再構成された画像が存在すると予想される学習された部分空間、即ち可能な画像のベクトル空間の部分空間を使用することである。

【0127】

物質分解は、例えば、測定された各投影光線を独立して処理することにより、投影画像データまたはサイノグラムデータに対して実施することができる。この処理は、最大尤度分解（maximum likelihood decomposition）、または撮像されたオブジェクトの物質組成に関する事前確率分布（prior probability distribution）が仮定される最大事後分解（maximum a posteriori decomposition）の形式をとることができる。また、入力カウントのセットから出力カウントのセットへの線形又はアフィン変換（affine transform）、Ａｌｖａｒｅｚによって例示的に説明されるようなＡテーブル推定器（参照：Med Phys. 2011 May; 38(5)：232－2334）、低次多項式近似（low-order polynomial approximation）例えばLee et al. (IEEE Transactions on Medical Imaging (Volume: 36, Issue: 2, Feb. 2017: 560-573)またはルックアップテーブルによって例示的に記述されたニューラルネットワーク推定器（neural network estimator）である。あるいは、物質分解法は、複数の光線を一緒に処理してもよいし、１段階または２段階の再構成アルゴリズムで構成されてもよい。

【0128】

ＷＯ２０１７／２２３５６０は、ディープラーニングアルゴリズムを使用して再構成画像を生成する断層画像再構成方法を開示している。例えば、最初に第１のアルゴリズムを適用して中間データを取得し、次に第２のアルゴリズムを中間画像に適用する。

【0129】

米国特許第９，０３６，８８５号は、例えば学習収縮関数を用いて、スパース化変換及び収縮を用いて投影データを処理する投影領域ノイズ除去処理方法を開示している。

【0130】

ＷＯ２０１５／０２８９７５は、少なくとも２つの基底物質サイノグラムを受信し、非相関フィルタを用いて基底物質サイノグラムのノイズを除去する方法を開示している。

【0131】

Optical Engineering 58(1), 013104に掲載されたChenとLiによる論文（An article by Chen and Li in Optical Engineering 58(1), 013104）では、ディープニューラルネットワークを用いてスペクトルＣＴデータの複数の物質分解を行う方法が開示されており、単にニューラルネットワークによって実際の物質分解を行うだけである。

【0132】

Scientific Reports第9巻、記事番号：17709（2019）のPoirotらによる論文（An article by Poirot et al. in Scientific Reports volume 9, Article number: 17709 (2019)）では、畳み込みニューラルネットワークを用いてデュアルエナジーＣＴ画像から非造影シングルエナジーＣＴ画像を生成する方法が開示されている。

【0133】

第１の態様によれば、エネルギー分解Ｘ線データを処理するための方法の非限定的な例が提供され、方法は以下のステップを含む。
－エネルギー分解Ｘ線データの少なくとも１つの表現を収集するステップ
－エネルギー分解Ｘ線データの前記少なくとも１つの表現に対して少なくとも２つの基底物質分解を実行し、少なくとも２つの元の基底画像表現セットをもたらすステップ
－元の基底画像表現セットの少なくとも２つから少なくとも２つの基底画像表現を取得又は選択するステップ
－前記取得又は選択された基底画像表現を、機械学習アルゴリズムに基づくデータ処理（データ処理手順など）で処理し、出力画像データの少なくとも１つの表現を生成するステップ。

【0134】

つまり、エネルギー分解されたＸ線データに基づく（Ｘ線）画像再構成法として表現することができる

【0135】

第２の態様によれば、以下の構成を有する、エネルギー分解Ｘ線データを処理するための対応するシステムの非限定的な例も提供される。
－システムは、エネルギー分解Ｘ線データの少なくとも１つの表現を収集するように構成され、
－システムは、画像データの前記表現に対して少なくとも２つの基底物質分解を実行して、少なくとも２つの元の基底画像表現セットを生成するように構成され、
－システムは、前記元の基底画像表現セットのうちの少なくとも２つから少なくとも２つの基底画像表現を取得又は選択するように構成され、
－出力画像データの少なくとも１つの表現を生成するために、機械学習に基づくデータ処理（データ処理手順など）を使用して、取得又は選択された基底画像表現を処理するように構成される。

【0136】

つまり、エネルギー分解されたＸ線データに基づいて（Ｘ線）画像を再構成するためのシステムとして表現できる。

【0137】

本発明の例示的な実施形態では、エネルギー分解Ｘ線データの少なくとも１つの表現を収集する前記ステップ又は構成は、ＣＴイメージングシステムによって行われる。

【0138】

本発明の例示的な実施形態では、エネルギー分解Ｘ線データの少なくとも１つの表現を収集する前記ステップ又は構成は、マルチビン光子計数Ｘ線検出器とも呼ばれるエネルギー分解光子計数検出器によって行われる。

【0139】

あるいは、エネルギー分解Ｘ線データの少なくとも１つの表現を収集する前記ステップまたは構成は、マルチＸ線管収集、低速または高速ｋＶスイッチング収集、多層検出器または分割フィルタ収集の方法によって行われる。

【0140】

本発明の例示的な実施形態では、前記少なくとも１つの元の基底表現又は出力基底表現は、サイノグラム、投影Ｘ線画像、又は再構成されたＣＴ画像のうちの少なくとも１つを含む。

【0141】

本発明の例示的な実施形態では、選択された基底画像を処理する前記ステップ又は構成は、以下のステップ又は構成を含む。
－前記選択された基底画像データを、非基底画像データの少なくとも１つの表現と組み合わせることであって、非基底画像データは、例えば、ビン計数サイノグラム・データ、対数正規化ビン計数サイノグラム・データ、ビンオーバーサムサイノグラム又は対数正規化データ、再構築されたＣＴ画像データ、及び／又はセグメンテーション・マスク・データ（bin counts sinogram data, log-normalized bin counts sinogram data, sum-over-bins sinogram or log-normalized data, reconstructed CT image data, and/or segmentation mask data）である。
－機械学習に基づくデータ処理方法を用いて、選択された基底画像データを非基底画像データの表現とともに処理すること。

【0142】

本発明の例示的な実施形態では、出力画像データの前記少なくとも１つの表現は、元の基底画像セットとは異なる基底関数に対応する基底画像セットの一部を形成する基底画像表現の近似値である。

【0143】

本発明の例示的な実施形態では、出力画像データの前記少なくとも１つの表現は、基底画像表現のセットからなり、前記出力セットにおける基底画像表現の数は、元の分解における基底画像のセットのいずれかにおける基底画像の最大数よりも大きい。

【0144】

本発明の例示的な実施形態では、前記出力基底画像表現セットに対応する基底関数は、前記元の基底画像表現セットに対応する基底関数の和集合のサブセットであるか、又はそれに等しい。あるいは、前記出力基底画像表現セットに対応する基底関数は、前記元の基底画像表現セットに対応する基底関数の和集合のサブセットでなくてもよい。

【0145】

本発明の例示的な実施形態では、前記選択された基底画像表現の数は、前記出力基底画像表現セットに対応する基底関数の数以上である。

【0146】

本発明の例示的な実施形態において、前記元の基底画像表現セットは、２基底分解を含む／であり、前記少なくとも１つの出力画像表現セットは、３基底分解を含む。

【0147】

本発明の例示的な実施形態では、前記基底物質分解のうちの少なくとも１つは、最尤法によって実行される。代替として、前記基底物質分解の少なくとも１つは、線形又はアフィン変換、Ａテーブル法、低次多項式近似、ニューラル・ネットワーク、又はルックアップテーブル（a linear or affine transformation, an A-table method, a low-order polynomial approximation, a neural network or a look-up table）を通じて実行される。

【0148】

本発明の例示的な実施形態では、前記基底物質分解のうちの少なくとも１つは、１段階の基底画像再構成法によって、又は２段階の基底画像再構成法によって実行される。

【0149】

本発明の例示的な実施形態では、前記機械学習は、畳み込みニューラル・ネットワークに基づくことができる機械学習アーキテクチャ及び／又はアルゴリズムを含む。あるいは、前記機械学習アーキテクチャ及び／又はアルゴリズムは、サポートベクターマシン又は決定木ベースの方法に基づくことができる。

【0150】

本発明の例示的な実施形態では、前記畳み込みニューラル・ネットワークは、残差ネットワーク（ＲｅｓＮｅｔ）、残差エンコーダ・デコーダ、Ｕ－Ｎｅｔ、ＡｌｅｘＮｅｔ、又はＬｅＮｅｔアーキテクチャ（a residual network (ResNet), residual encoder-decoder, U-Net, AlexNet or LeNet architecture）に基づく。代替として、畳み込みニューラル・ネットワークに基づく前記機械学習アルゴリズムは、勾配降下アルゴリズム、主双対アルゴリズム、又は乗数の交互方向法（ＡＤＭＭ）アルゴリズム（gradient descent algorithm, a primal-dual algorithm or an alternating direction method of multipliers (ADMM) algorithm）に基づくアンロール最適化法（an unrolled optimization method）に基づくことができる。

【0151】

本発明の例示的な実施形態では、前記畳み込みニューラル・ネットワークは、ネットワークアーキテクチャの一部として、少なくとも１つの順投影又は少なくとも１つの逆投影を含む。

【0152】

本発明の例示的実施形態では、前記畳み込みニューラル・ネットワークは、ファントム画像、患者画像、又は非エネルギー分解ＣＴで画像化された患者から生成された数値ファントムの少なくとも１つで訓練される。

【0153】

本発明の例示的な実施形態では、前記畳み込みニューラルネットワークは、第１のタイプのＣＴ画像、例えば非エネルギー分解ＣＴ画像について移転学習を使用して訓練され、第２のタイプのＣＴ画像、例えば光子計数ＣＴ画像について移転学習を使用して使用するように適応される。

【0154】

本発明の例示的な実施形態では、選択された基底画像を処理する前記ステップまたは構成は、以下を含む。
－前記選択された基底画像データを、事前画像上のボリューム制約、学習済み辞書として表される事前分布、正規化ペナルティ、学習済み正規化ニューラルネットワーク、もしくは学習済み部分空間、または統計的解剖学的形状モデル（a volume constraint on a prior image, a prior distribution represented as a learned dictionary, a regularization penalty, a learned regularization neural network, or a learned subspace, or a statistical anatomical shape model）などの、事前情報の少なくとも１つの表現と結合するステップと
－前記選択された基底画像データを、前記事前情報の表現を用いて、機械学習に基づく前記データ処理で処理するステップ。

【0155】

本発明の例示的な実施形態では、選択された基底画像を処理する前記ステップ又は構成は、
－前記選択された基底画像データを、収集システムにおいて物理的プロセスの少なくとも１つの計算表現、例えば、コンプトン散乱、蛍光、電荷共有、パイルアップ又は幾何学的ぼかしと組み合わせるステップと、
－機械学習アルゴリズムに基づく物理プロセスの計算表現を用いて、選択された基底画像データを処理するステップとを含む。

【0156】

あるいは、選択された基底画像を処理する前記ステップ又は構成は、前記選択された基底画像データを、尤度関数又は尤度関数の導関数などの測定データの統計モデルの少なくとも１つの計算表現と組み合わせることを含む。

【0157】

本発明の例示的実施形態では、前記機械学習アーキテクチャ及び／又はアルゴリズムは、元の基底画像表現よりもノイズが低い１つ又は複数の出力基底画像表現を生成するように訓練される。あるいは、前記機械学習アーキテクチャ及び／又はアルゴリズムは、元の基底画像表現よりも偏りが低い１つ又は複数の出力基底画像表現を生成するように訓練されてもよい。さらに別の代替では、前記機械学習アーキテクチャ及び／又はアルゴリズムは、元の基底画像表現よりも高い空間解像度を有する１つ又は複数の出力画像データ表現を生成するように訓練され得る。

【0158】

一般に、ディープラーニングは、人工ニューラル・ネットワーク又は表現学習による同様のアーキテクチャに基づく機械学習方法に関連している。学習は、教師あり、半教師あり、又は教師なしで行うことができる。ディープ・ニューラル・ネットワーク、ディープ・ビリーフ・ネットワーク、リカレント・ニューラル・ネットワーク、畳み込みニューラル・ネットワーク（deep neural networks, deep belief networks, recurrent neural networks and convolutional neural networks）などのディープラーニングシステムは、コンピューター・ビジョン、音声認識、自然言語処理、ソーシャルネットワーク・フィルタリング、機械翻訳、ボードゲームプログラム（computer vision, speech recognition, natural language processing, social network filtering, machine translation, and board game programs）など、さまざまな技術分野に適用されている。彼らは、人間の専門家のパフォーマンスに匹敵し、場合によってはそれを超える結果を生み出す。

【0159】

ディープラーニング（深層学習）における「ディープ：深い」という形容詞は、ネットワークに複数の層を用いることに由来する。初期の研究では、線形パーセプトロンは万能分類器になりえないこと、非多項式活性化関数を持ち、非限定幅の隠れ層を持つネットワークは万能分類器になりうることが示された。ディープラーニングは、現代的なバリエーションであり、無限の層数を持ち、最適化された実装を可能にし、かつ穏やかな条件下で理論的な普遍性を保持する。ディープラーニングでは、効率性、訓練性、理解しやすさのために、層は異種であってもよく、生物学的知見に基づくコネクショニストモデルから大きく外れてもよいことになっている。より良い理解のために、次に、提案された技術の例示的かつ非限定的な例を説明する。

【0160】

物質分解問題は、計算コストと精度の両面で解決が困難な非線形逆問題である。最も受け入れられているソリューションには、ポアソン統計を使用した最尤（ＭＬ）推定に基づく最適化が含まれる。これは、考慮されるフォワードモデルに大きく依存するモデルベースのアプローチである。一方、データ駆動型のアプローチでは、利用可能な測定データを利用する。これはディープラーニングの場合で、事前の情報を高速な方法で提供できる。

【0161】

例えば、光子計数ＣＴ（ＰＣＣＴ：Photon-Counting CT）では、通常、３つの材料（物質）が臨床解析に関係する。しかし、実際には、２物質基底のＭＬ推定は３物質基底の推定よりも収束性がよく、安定で、ノイズも少ない。しかしながら、２物質基底の推定は高いバイアス（偏り）を持つ。この偏りを補正するために、実測した学習データを用いたディープラーニングの後処理を行うことで、２物質ＭＬ推定の利点（安定性、収束性）を維持しつつ、データ駆動型の正則化を行い、偏りの補正を行うことができる。したがって、２物質基底分解の利点と偏りを制御するデータ駆動型のアプローチを組み合わせることが提案される。

【0162】

特定の非限定的な例では、アイデアは、２物質分解から３物質分解に移行することである。この例では、提案された技術は、２物質基底ＭＬ推定の安定性と収束性を利用し、偏りを補正するためにデータ駆動型の後処理を考慮する。一例として、物質分解は投影（サイノグラム）領域で行われてもよく、すなわち、ＣＴ画像の再構成は独立した第２段階であってもよい。

【0163】

スキャンされたオブジェクトが３つの物質、例えば骨、軟組織、及び造影剤（ヨウ素など）を含むと仮定して、２つの基底物質、すなわち（ｉ）骨及び軟組織、（ｉｉ）骨及びヨウ素、及び（ｉｉｉ）軟組織及びヨウ素のあらゆる可能な組み合わせを得るために一連の２物質推定を定義することにする。例えば、提案技術のディープラーニングアーキテクチャ及び／又は手順は、この６つの偏った推定値の集まりを入力として考慮し、出力として骨、軟組織、及びヨウ素の偏りのない推定値を提供することができる。

【0164】

例えば、偏った物質分解サイノグラム（２基底ＭＬ推定値から）を受け取り、３つの補正された物質分解サイノグラムを出力として提供する教師あり機械学習アーキテクチャおよび／または手順（入力－出力ペアの学習データ）を使用することが実行可能である。

【0165】

例えば、設計されたアーキテクチャは、５つの残差ブロック（ＲｅｓＢｌｏｃｋ）で構成される残差ネットワークであってもよい。各ブロックは、３つの畳み込み層、それに続くバッチ正規化および整流線形ユニット（ＲｅＬｕ：rectified linear units）活性化を含んでもよい。特定の例では、各畳み込み層の出力チャンネルは、それぞれ３２、３２、及び３である。

【0166】

図９は、提案されたディープラーニング偏り除去（proposed deep learning debiasing）のためのスキーム及び／又はアーキテクチャの一例を示す模式図である。非限定的な例では、例えば８個のサイノグラムのセットは、２つ以上の基底画像表現を生成するための２基底分解の入力として提供されてもよい。ディープラーニングニューラルネットワークは、得られた基底画像表現を、例えば投影空間における３つの物質分解に変換するために使用されてもよい。

【0167】

例として、ＲｅｓＢｌｏｃｋｓは、最適化問題の反復解の更新を模倣し、学習した畳み込みフィルタで学習データからの情報を追加することができる。各ブロックΨは、ｎ番目の更新ステップと見なすことができる。すなわち、（ｎ＋１）番目のブロックであるａｎ＋１の入力は、ｎ番目のブロックであるａｎの出力と、そのフィルタ重みθｎによってパラメタ化された最後のブロックΨｎ の演算結果との差分である。

【0168】

【数9】

【0169】

Ψは変数ａｎの関数であり、θｎによってパラメータ化されることを理解されたい。図１０は、ＲｅｓＢｌｏｃｋの一例を示す模式図である。

【0170】

図１１は、２物質分解からの偏った結果の後の一連の残差ブロック（ＲｅｓＢｌｏｃｋｓ）を含む偏り除去ディープラーニングスキーム及び／又はアーキテクチャの例を示す概略図である。

【0171】

本発明と研究開発の一環として、本技術のシミュレーションによる概念実証を行った。一例として、実現可能なセットアップを以下に紹介する。

【0172】

Ｓｈｅｐｐ－Ｌｏｇａｎヘッド（Shepp-Logan head）に基づく、シミュレートされたＰＣＣＴデータのデータセット。ＰＣＣＴシステムとポアソンノイズのフォワードモデル（forward model：順モデル）を考える。ファントムのサイズは、画像空間で６４×６４ピクセル、投影空間で１００×１００ピクセルである。考慮される物質は、ＰＶＣ（骨）、ＰＥ（軟部組織）、及びヨウ素である。データに変動性を追加するために、各サンプルのサイズ、回転、及びヨウ素領域の位置をランダムに変更した。図１２は、シミュレートされたデータセットのサンプルの例を示す概略図であり、シェップ・ローガン（Shepp-Logan）の異なる部分が異なる物質を模倣している。

【0173】

２物質分解ＭＬ関数を考慮し、推計のペアを実行する：ＰＶＣ／ＰＥ、ＰＶＣ／Ｉ、及び ＰＥ／Ｉ。これらの推定値は、３つの物質の混合物に属しているため、偏りがある。

【0174】

これらの偏った推定値は、ニューラル・ネットワークの入力となる。出力は、物質ごとに１つずつ、３つの偏りのないノイズ除去された推定値（ＰＶＣ、ＰＥ、Ｉ）で構成される。

【0175】

ニューラルネットワークの学習には、平均二乗誤差（ＭＳＥ）学習損失を用いたＡｄａｍ最適化スキーム（Adam optimization scheme）を考える。この例では、訓練に２５０サンプル、テストに１００サンプルを使用する。

【0176】

１００個のテストサンプルの評価では、ＭＳＥと結果の推定の偏りが大幅に改善されていることがわかる。以下の表１及び図１３を参照されたい。

【0177】

【表1】

【0178】

図１３（ａ）～（ｆ）は、シェップ・ローガンによるシミュレーションファントムに対して提案されたディープラーニングに基づく偏り緩和を適用した、定性的な結果の例を示す概略図である。

【0179】

図１３（ａ）は、ＰＶＣ（第１列）、ＰＥ（第２列）及びヨウ素（第３列）の物質の真のサイノグラム（true sinograms of materials）の例を示す。

【0180】

図１３（ｂ）は、偏った２物質ＭＬからの６つの解のうちの最良の３つの例を示す。

【0181】

図１３（ｃ）は、提案された偏りのないサイノグラムの例を示す。

【0182】

図１３（ｄ）は、フィルタ逆投影（ＦＢＰ）による再構成後の真の基底画像の例を示す。

【0183】

図１３（ｅ）は、ＦＢＰ後の偏った推定の例を示している。

【0184】

図１３（ｆ）は、ＦＢＰ後の提案された偏りのない例を示す。

【0185】

以下では、特にこの技術/手法を実際の臨床データに対応させるために拡張する必要がある場合に、結果を向上させる可能性のある提案技術の代替実装を紹介する。
（ｉ）ニューラル・ネットワーク
o提案する残差ネットワークのスケーリング：入力の次元が大きい場合、畳み込み層のサイズとブロックの数はそれに応じて設計することができる。すなわち、各ブロックの畳み込み層を増やす、各層の畳み込み出力チャンネルを増やす、または５つ以上のＲｅｓＢｌｏｃｋを増やすなど。
o ネットワークアーキテクチャの変更：提案した残差ネットワークは、入力のマルチスケール情報も利用するU-netに置き換えることができる。例えば、"U-Netは"O. Ronnebergerらによる「Convolutional Networks for Biomedical Image Segmentation, Medical Image Computing and Computer-Assisted Intervention－ MICCAI (2015)」を参照することができる。代替的に、提案されたネットワークは、畳み込み層と完全連結層（convolutional and fully-connected layers）を組み合わせた別の一般的なネットワークアーキテクチャ、例示的にAlexNetまたはLeNetによって置き換えられることもあり得る。
o ニューラルネットワークのアーキテクチャの強化：幅方向（投影優位）、高さ方向（検出優位）のフィルタなど、正方形ではなく長方形の畳み込みフィルタを検討する。参考文献として、”Machine Friendly Machine Learning： Interpretation of Computed Tomography Without Image Reconstruction”, by Lee, H., et al., Sci Rep 9, 15540 (2019)を参照することができる。
o最適化問題の反復更新ステップをより明示的に模倣：それゆえ、その定義にシステムのフォワードモデルに関する事前情報を含むニューラルネットワークアーキテクチャを考察する。例えば、“Learned Primal Dual Reconstruction”, by Adler, J and Oktem, O., IEEE Transactions on Medical Imaging (2018)を参照することができる。
o学習損失をカスタマイズ：基本的な平均二乗誤差（MSE：basic mean-squared error）の代わりに、偏り、コントラスト対ノイズ比（CNR：contrast-to-noise ratio）、正則化損失（regularized loss）など、学習損失において関心のある品質指標を考慮することができる。
oネットワークへの入力として、複数の形式のデータを取り入れる：例えば、エネルギービンサイノグラムと物質分解サイノグラムを組み合わせて、ネットワークへの入力として使用することができる。また、エネルギービンサイノグラム及び／又は物質分解サイノグラムと、エネルギービン画像及び／又は物質分解画像とを組み合わせて、ネットワークへの入力とすることも可能である。また、尤度関数の微分値や基底画像ノイズのCramer-Rao下限値の表現を、物質分解サイノグラムと組み合わせてネットワークに入力することも可能である。上記の入力データのすべての組み合わせは、２つ以上の基底物質サイノグラム、または少なくとも２つの異なる基底画像表現セットからの画像と共に使用することができる。
oスペクトルアテンションモジュール（spectral attention modules）でネットワークの表現能力を高める：エネルギービンデータはネットワークの入力にもなり、ネットワークアーキテクチャのスペクトルアテンションモジュールでいくつかの情報が抽出される可能性がある。例えば、"Learning Spectral-Spatial Prior Via 3DDNCNN for Hyperspectral Image Deconvolution" by X. Wang et al., ICASSP (2020)などを参照することができる。
（ｉｉ）別の分解：
o組織／診断／解剖学的分解：ディープラーニングは柔軟な学習プロセスを可能にし、そこから基礎となるモデルを定義する必要なく、入力データの興味深い変換を推論することができる。したがって、必ずしも３基底分解を得る必要はなく、むしろ（３よりもさらに大きな）組織型分解、診断領域分解、解剖学的断面分解を得ることができる場合がある。ラベル付きの入出力学習データのセットがあれば、同様の設定でこれらの目的に使用できる。
（ｉｉｉ）その過程でＣＴ再構成を検討する。
oステップの分解と再構成：提案された方法では、ＣＴ再構成（サイノグラムから画像空間へ）が独立したステップであると仮定される。ニューラルネットワークは、２つのステップを同時に処理することができる。例えば、偏った物質分解されたサイノグラムを、偏った物質分解された再構成画像に直接変換するように学習させた非転写反復最適化ネットワークが考えられる。このようにして、ニューラルネットワークは、より最適な方法（計算と精度）で、両方のステップ（分解と再構成）を解くことができる。例えば、"Learned Primal-Dual Reconstruction" by J. Adler and O. Oktem, IEEE Transactions on Medical Imaging, (2018)を参照されたい。このような展開反復最適化ネットワークは、勾配降下最適化法、プライマル－デュアル最適化法、または交互方向乗数法（a gradient descent optimization method, a primal-dual optimization method, or an alternating direction method of multipliers method）に基づくものであってよい。この種の繰り込み最適化法は、他のニューラルネットワークアーキテクチャと比較して、ネットワークの重みが少なく、少量の訓練データでより良く動作するなど、有利な特性を有する。このようなワンステップ法は、例えば、基底画像データ及びエネルギー分解された測定データの組合せを入力として取ることができる。あるいは、ワンステップ再構成法は、反復アルゴリズムとして実装されてもよい。本発明のこの実施形態の非限定的な例は、図１４に示されており、機械学習システムが、少なくとも２つの基底物質分解から物質分解サイノグラムを処理し、出力として物質分解画像のセットを生成するものである。
oあるいは、元の基底画像セットを生成するために使用される物質分解ステップは、出力として再構成された基底画像をもたらす、ワンステップ再構成方法を構成することができる。これらは、次に、機械学習アルゴリズムによって処理され、１つ以上の出力基底画像をもたらす。本発明のこの実施形態の非限定的な例は、図１５に示されており、機械学習システムは、入力エネルギービンサイノグラムを物質分解再構成画像に変換する少なくとも２つのワンステップ画像再構成から物質分解画像を処理し、出力として物質分解画像セットを生成する。さらなる代替案は、２ステップ基底分解法を用いて元の基底画像セットを生成することであり、これは、基底物質分解が各投影光線に対して別々に実施される第１のステップと、得られた基底サイノグラムが画像再構成アルゴリズムへの入力として取られて再構成された基底画像をもたらし、これが機械学習アルゴリズムによって処理されて１以上の出力基底画像をもたらす第２のステップとから構成されている。本発明のこの実施形態の非限定的な例が図１６に示されており、機械学習システムは、それぞれが基底物質分解ステップに続いて画像再構成ステップからなる少なくとも２つの２ステップ画像再構成から物質分解画像を処理し、物質分解画像のセットを出力として生成する。前記１ステップまたは２ステップ分解方法は、任意に、事前画像または体積保存制約などの事前情報を組み込んでもよいし、サイノグラムまたは画像内のノイズレベルまたは相関構造に関する情報、例示的にCramer-Rao下限に関する情報を組み込んでもよい。

【0186】

本明細書に記載の機構及び構成は、さまざまな方法で実装、組み合わせ、及び再構成できることが理解されるであろう。

【0187】

例えば、実施形態は、ハードウェアで、又は適切な処理回路による実行のために少なくとも部分的にソフトウェアで、又はそれらの組み合わせで実装され得る。

【0188】

本書に記載されたステップ、機能、手順、および／またはブロックは、汎用電子回路および特定用途向け回路の両方を含むディスクリート回路（discrete circuit）または集積回路技術などの任意の従来技術を使用してハードウェアに実装することができる。

【0189】

代替的に、または補完的に、本明細書に記載のステップ、機能、手順、および／またはブロックの少なくとも一部は、１つまたは複数のプロセッサまたは処理ユニットなどの適切な処理回路による実行のために、コンピュータプログラムなどのソフトウェアで実装することができる。

【0190】

図１７は、実施形態によるコンピュータ実装の一例を示す概略図である。この特定の例では、システム２００は、プロセッサ２１０とメモリ２２０とを備え、メモリは、プロセッサによって実行可能な命令を含み、それによって、プロセッサは、本明細書に記載のステップ及び／又はアクションを実行するように動作可能である。命令は、典型的には、コンピュータプログラム２２５、２３５として構成され、メモリ２２０に予め構成されてもよいし、外部メモリデバイス２３０からダウンロードされてもよい。任意選択で、システム２００は、入力パラメータ（複数可）及び／又は結果出力パラメータ（複数可）などの関連データの入力及び／又は出力を可能にするために、１つまたは複数のプロセッサ２１０及び／又はメモリ２２０に相互接続され得る入力／出力インタフェース２４０を具備する。

【0191】

特定の例では、メモリは、プロセッサによって実行可能な命令のセットを含み、それによって、プロセッサは、
－エネルギー分解されたＸ線データの少なくとも１つの表現を取得し、
－エネルギー分解Ｘ線データの前記少なくとも１つの表現に基づいて、少なくとも２つの基底物質分解を実行して、少なくとも２つの元の基底画像表現セットを生成すること、
－前記元の基底画像表現セットの少なくとも２つから少なくとも２つの基底画像表現を抽出し、
－出力画像データの少なくとも１つの表現を生成するために、機械学習に基づくデータ処理によって抽出された基底画像表現を処理する。

【0192】

「プロセッサ」という用語は、特定の処理、決定、又は計算タスクを実行するためにプログラムコード又はコンピュータ・プログラム命令を実行できるシステム又はデバイスとして一般的な意味で解釈されるべきである。

【0193】

したがって、１つ又は複数のプロセッサを含む処理回路は、コンピュータ・プログラムを実行するときに、本明細書で説明するような明確に定義された処理タスクを実行するように構成される。

【0194】

処理回路は、上述のステップ、機能、手順及び／又はブロックを実行するためだけに専用である必要はなく、他のタスクを実行することもできる。

【0195】

提案された手法はまた、このようなコンピュータプログラムを格納したコンピュータ読み取り可能な媒体２２０、２３０からなコンピュータプログラム製品を提供するものである。

【0196】

例として、ソフトウェアまたはコンピュータプログラム２２５、２３５は、コンピュータプログラム製品として実現されてもよく、これは通常、コンピュータ可読媒体２２０、２３０、特に不揮発性媒体上に担持または格納される。コンピュータ可読媒体は、リードオンリーメモリ（ＲＯＭ）、ランダムアクセスメモリ（ＲＡＭ）、コンパクトディスク（ＣＤ）、デジタルバーサタイルディスク（ＤＶＤ）、ブルーレイディスク、ユニバーサルシリアルバス（ＵＳＢ）メモリ、ハードディスクドライブ（ＨＤＤ）記憶装置、フラッシュメモリ、磁気テープ、または他の任意の従来のメモリデバイスを含むが、これに限定されない１以上の取り外し可能または取り外し不能メモリ装置を含んでも良い。したがって、コンピュータプログラムは、その処理回路による実行のために、コンピュータまたは同等の処理装置の動作メモリにロードされ得る。

【0197】

メソッドフローは、１つ以上のプロセッサによって実行される場合、コンピュータアクションフローと見なすことができる。対応するデバイス、システム及び／又は装置は、機能モジュールのグループとして定義されてもよく、プロセッサによって実行される各ステップは、機能モジュールに対応する。この場合、機能モジュールは、プロセッサ上で実行されるコンピュータプログラムとして実装される。したがって、デバイス、システム及び／又は装置は、代替的に、機能モジュールのグループとして定義されてもよく、機能モジュールは、少なくとも１つのプロセッサ上で実行されるコンピュータプログラムとして実装される。

【0198】

したがって、メモリに常駐するコンピュータ・プログラムは、プロセッサによって実行されると、本明細書に記載のステップ及び／又はタスクの少なくとも一部を実行するように構成された適切な機能モジュールとして編成され得る。

【0199】

また、主にハードウェアモジュールで実現することも可能であり、ハードウェアと交互に実現することも可能である。ソフトウェアとハードウェアをどの程度使い分けるかは、純粋に実装の選択である。

【0200】

上述した実施形態は、単に例として与えられたものであり、提案された技術はこれに限定されないことが理解されるべきである。添付の特許請求の範囲によって定義される本発明の範囲から逸脱することなく、実施形態に対して様々な修正、組み合わせ、変更を行うことができることは当業者には理解されよう。特に、異なる実施形態における異なる部品ソリューションは、技術的に可能であれば、他の構成で組み合わせることができる。

【符号の説明】

【0201】

１０：Ｘ線源 １１：ガントリ １２：テーブル １３：アイソセンタ ２０：Ｘ線検出器 ２１：検出器素子 ２５：アナログ処理回路 ３０：画像処理システム ３１：入力 ３２：基底物質分解システム ３３：機械学習システム ３４：出力 ４０：デジタル処理回路 ４１：Ｘ線コントローラ ４２：ガントリコントローラ ４３：テーブルコントローラ ４４：検出器コントローラ ４５：メモリ ５０：コンピュータ ６０：オペレータコンソール １００：Ｘ線イメージングシステム ２００：システム ２１０：プロセッサ ２２０：メモリ ２２５：コンピュータプログラム ２３０：外部メモリデバイス ２３５：コンピュータプログラム ２４０：入力／出力インタフェース ３０１：デジタル／アナログ変換器 ３０２：比較器 ３０３：デジタルカウンタ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】