(57)【要約】  陽子コンピュータ断層撮影に関連するシステムと、デバイスと、方法とが開示される。幾つかの実施態様では、陽子の検出は、各陽子についてオブジェクトの前及び後のトラック情報をもたらすことができ、それによってオブジェクト内での各陽子の、可能性が高い経路を求めることが可能になる。さらに、各陽子が受けるエネルギー損失の測定によって、所与の可能性の高い経路が所与のエネルギー損失をもたらすという判定が可能になる。こうしたデータの集合によって、オブジェクトの特徴付けが可能になる。エネルギー損失に関して、こうした特徴付けは、オブジェクトの相対阻止能の画像マップを含むことができる。限定はしないが、総変動等のメリット関数の優秀化を含むこうした画像を取得するための種々の再構成方法が開示される。幾つかの実施態様では、種々の形態の総変動優秀化方法は、計算的に効率的であり、かつ、計算時間を低減しながら、優れた結果をもたらすことができる。幾つかの実施態様では、こうした方法は、比較的低い陽子線量を使用して、高品質陽子ＣＴ画像をもたらすことができる。【選択図】図７
2013519452.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータ断層撮影を実施するための方法であって、
オブジェクトを通過する複数の陽子について測定データを取得することであって、前記測定データは、前記陽子のそれぞれについて第１のトラック及び第２のトラックに関する情報を含み、該第１のトラック及び該第２のトラックは、それぞれ、前記陽子の、前記オブジェクトの通過の前及び後の該陽子の軌跡に対応し、前記測定データは、前記陽子の、前記オブジェクトの通過によって生じる各陽子の相互作用量に関する情報を更に含む、取得することと、
各陽子について、前記第１のトラック及び前記第２のトラックに少なくとも部分的に基づいて前記オブジェクト内で前記陽子がとる経路を推定することと、
前記陽子の、前記オブジェクトの通過が、連立方程式Ａｘ＝ｂとして表されるか又は表すことが可能であるように、前記陽子の前記相互作用量及び前記推定経路を配列することであって、ここで、ｘは、前記オブジェクトに関連するパラメータの分布であり、ｂは、前記オブジェクト内の前記陽子のそれぞれの経路に沿う相互作用から生じる該陽子の前記相互作用量を表し、Ａは、ｂを得るためにｘに作用する作用素であり、該作用素Ａは、前記オブジェクト内の前記陽子の前記推定経路に関する情報を有し、前記連立方程式は、複数の解を有するように構成される、配列することと、
前記連立方程式について初期解を推定することと、
前記複数の解の中で１つ又は複数の実行可能解を探すことであって、各実行可能解は、既存の解に摂動を与えることによって得られ、前記既存の解の前記摂動なしで得られる別の解よりも、前記オブジェクトパラメータ分布の再構成に関連する量について優れた特性を有する、探すことと、
前記１つ又は複数の実行可能解のうちの選択された１つの解に基づいて前記オブジェクトパラメータ分布を計算することと、
を含む、コンピュータ断層撮影を実施するための方法。
【請求項２】
前記陽子の前記相互作用量は、前記陽子の、前記オブジェクトの通過によって生じる前記陽子のエネルギー損失を含む、請求項１に記載の方法。
【請求項３】
ｂは、前記陽子の前記推定経路に沿う前記相互作用量の積分値を表す、請求項１に記載の方法。
【請求項４】
前記連立方程式Ａｘ＝ｂは、連立一次方程式を含む、請求項１に記載の方法。
【請求項５】
前記作用素Ａは行列を含む、請求項４に記載の方法。
【請求項６】
前記選択された実行可能解は、前記複数の解の中で最適解でない実行可能解を含む、請求項１に記載の方法。
【請求項７】
前記オブジェクトパラメータ分布を計算することは、３Ｄオブジェクトパラメータ分布を計算することを含む、請求項１に記載の方法。
【請求項８】
前記計算されたオブジェクトパラメータ分布に基づいて前記オブジェクトの断層画像のアレイを形成することを更に含む、請求項１に記載の方法。
【請求項９】
前記オブジェクトパラメータ分布は、電子密度ベース量の分布を含む、請求項１に記載の方法。
【請求項１０】
前記電子密度ベース量は、実質的に一様な物質に対する相対的陽子阻止能を含む、請求項９に記載の方法。
【請求項１１】
前記実質的に一様な物質は水を含む、請求項１０に記載の方法。
【請求項１２】
前記再構成されるオブジェクトパラメータ分布は、前記電子密度ベース量の３次元分布を含む、請求項９に記載の方法。
【請求項１３】
前記経路を推定することは、前記陽子の最尤経路を推定することを含む、請求項１に記載の方法。
【請求項１４】
前記オブジェクトパラメータ分布の再構成に関連する量は、前記再構成されたオブジェクトパラメータ分布の総変動を含む、請求項１に記載の方法。
【請求項１５】
前記総変動の前記優れた特性は、前記総変動の、より低い値を含む、請求項１４に記載の方法。
【請求項１６】
前記再構成されたオブジェクトパラメータ分布の前記総変動は、ｐｋで評価される関数
【数１】
2013519452.tif 000022
の量を含み、ここで、ｐｋは、相互作用パラメータ分布ｘのベクトル表現ｘｋの２次元Ｊ×Ｊ表現である、請求項１５に記載の方法。
【請求項１７】
前記解に関連する前記総変動は、他の実行可能解に関連する総変動よりも優れている、請求項１５に記載の方法。
【請求項１８】
前記初期解を推定することは、フィルタリングされた逆投影再構成解を計算することを含む、請求項１に記載の方法。
【請求項１９】
前記１つ又は複数の実行可能解を探すことは、
前記オブジェクトパラメータ分布ｘのベクトル表現ｘｋに摂動を与えることであって、それによって、摂動付与済みベクトルｙｋを得る、摂動を与えることと、
前記摂動付与済みベクトルｙｋに関連する前記再構成されたオブジェクトパラメータ分布に関連する量を評価することと、
前記摂動付与済みベクトルｙｋに関連する量が、前記未摂動ベクトルｘｋに関連する量より優れている場合、ｙｋを投影することであって、それによって、次のベクトル表現ｘｋ＋１を得る、ｙｋを投影することと、
の反復を実施することを含む、請求項１に記載の方法。
【請求項２０】
前記ベクトルｘｋに摂動を与えることは、ｙｋ＝ｘｋ＋βｋｖｋとなるようにｙｋを計算することを含み、ここで、βｋは、摂動の大きさを表し、ｖｋは、摂動ベクトルである、請求項１９に記載の方法。
【請求項２１】
前記摂動ベクトルｖｋは、ｘｋにおける前記再構成されたオブジェクトパラメータ分布に関連する量の正規化された劣勾配を負数としたものとして計算される、請求項２０に記載の方法。
【請求項２２】
ｙｋについて評価された量が、ｘｋについて評価された量以下である場合、前記摂動付与済みベクトルｙｋに関連する量は、前記未摂動ベクトルｘｋに関連する量に対して優れている、請求項１９に記載の方法。
【請求項２３】
前記投影することは、或る関連する凸集合上へのｘｋの数学的に規定された投影としてｘｋ＋１を計算することを含む、請求項１９に記載の方法。
【請求項２４】
前記投影することは、前記測定データが複数のブロックに分割されるように、ブロック反復投影アルゴリズムを使用して、超平面、半空間、超スラブ、又は他の凸集合上に投影することを含む、請求項２３に記載の方法。
【請求項２５】
前記投影することは、前記連立方程式の個々の超平面上への直交投影とともに、対角成分ごとの緩和を可能にするように構成された対角緩和した直交投影（ＤＲＯＰ）ベースアルゴリズムを使用して投影することを含む、請求項２４に記載の方法。
【請求項２６】
ＤＲＯＰベースアルゴリズムは、対角緩和した直交行列λｋＵｔ（ｋ）を含み、ここで、λｋは、ｋ番目の反復についての緩和パラメータであり、Ｕｔ（ｋ）は、ｔ番目のブロックについて対角要素ｍｉｎ（１，１／ｈｔｊ）を有する対角行列であり、ｈｔｊは、ベクトルｘｋのｊ番目のボクセルと交差するｔ番目のブロック内の陽子履歴の数である、請求項２５に記載の方法。
【請求項２７】
前記投影することは、次の反復に進む前に、前記ブロックが全て処理されるまで、前記ブロックを通して周期的に実施される、請求項２４に記載の方法。
【請求項２８】
反復は、処理される所与のブロックについて、前記投影することが該所与のブロックについてだけ実施されるように、各ブロックについて実施される、請求項２４に記載の方法。
【請求項２９】
前記反復を実施することは、実行可能性近接性チェックを計算することであって、前記再構成されたオブジェクトパラメータ分布に関連する量を評価するという更なるタスクに関する反復が、前記測定データとの所望の一致からはずれて前記解を操向しないことを検証する、計算することを更に含む、請求項１９に記載の方法。
【請求項３０】
陽子コンピュータ断層撮影を実施する方法であって、
オブジェクトを通過する複数の陽子について測定データを取得することと、
前記測定データに対する総変動優秀化に基づいて投影ベース再構成アルゴリズムを反復して適用することであって、前記オブジェクトの相対阻止能の分布を得る、適用することと、
を含む、陽子コンピュータ断層撮影を実施する方法。
【請求項３１】
前記相対阻止能分布に基づいて前記オブジェクトの視覚画像を形成することを更に含む、請求項３０に記載の方法。
【請求項３２】
陽子コンピュータ断層撮影システムであって、
オブジェクトを通過するのに十分な、選択された平均エネルギーを有する複数の陽子を送出するように構成された陽子送出システムと、
前記陽子のそれぞれについて、前記オブジェクトの前及び後の軌跡、並びに、前記オブジェクトを通過した後のエネルギーを測定するように構成された検出器システムと、
前記検出器システムから信号を読み出し、それによって前記陽子のそれぞれの前記軌跡及び前記エネルギーを表す測定データを得るように構成されたデータ収集システムと、
前記測定データを処理して画像再構成を実施し、それによって、前記オブジェクトのコンピュータ断層画像を得るように構成されたプロセッサであって、前記画像再構成は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを含む、プロセッサと、
を備える、陽子コンピュータ断層撮影システム。
【請求項３３】
陽子治療システムであって、
第１の平均エネルギー及び第２の平均エネルギーを有する陽子ビームを送出するように構成された陽子送出システムであって、前記第１の平均エネルギーは、第１のブラッグピークが、身体の一部分の内部のターゲット領域内の場所で起こるように選択され、前記第２の平均エネルギーは、前記陽子ビームが、前記身体の前記一部分を通過するように選択される、陽子送出システムと、
前記ターゲット領域に対する前記第１のエネルギーのビームの送出を容易にするように構成された第１の検出器システムと、
前記第２のエネルギーを有し、前記身体の前記一部分を通過する前記陽子のそれぞれについて、前記身体の前記一部分の前及び後の軌跡、並びに、前記身体の前記一部分を通過した後のエネルギーを測定するように構成された第２の検出器システムと、
少なくとも前記第２の検出器システムから信号を読み出し、それによって前記第２のエネルギーの陽子のそれぞれの前記軌跡及び前記エネルギーを表す測定データを得るように構成されたデータ収集システムと、
前記測定データを処理して画像再構成を実施し、それによって前記身体の前記一部分のコンピュータ断層画像を得るように構成されたプロセッサであって、前記画像再構成は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを含む、プロセッサと、
を備える、陽子治療システム。
【請求項３４】
有形コンピュータ可読記憶媒体であって、該有形コンピュータ可読記憶媒体上に記憶されたコンピュータ実行可能命令を有し、該コンピュータ実行可能命令は、１つ又は複数のコンピューティングデバイスを備えるコンピューティングシステムによって可読であり、前記コンピュータ実行可能命令は、
複数の陽子に関するデータを取得することであって、該データは、前記陽子のそれぞれについて第１のトラック及び第２のトラックに関する情報を含み、前記データは、前記第１のトラックと前記第２のトラックとの間での各陽子のエネルギー損失に関する情報を更に含む、取得することと、
各陽子について、前記第１のトラックと前記第２のトラックとの間の経路を推定することと、
総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを使用して、前記陽子のエネルギー損失及び前記経路に基づいて相対阻止能分布の断層撮影解析を実施することと、
を含む動作を、前記コンピューティングシステムに実施させるように、該コンピューティングシステム上で実行可能である、有形コンピュータ可読記憶媒体。
【請求項３５】
粒子放射線治療システムであって、
第１の平均エネルギー及び第２の平均エネルギーを有するイオンビームを送出するように構成された粒子放射線送出システムであって、前記第１の平均エネルギーは、第１のブラッグピークが、身体の一部分の内部のターゲット領域内の場所で起こるように選択され、前記第２の平均エネルギーは、前記イオンビームが、前記身体の前記一部分を通過するように選択される、粒子放射線送出システムと、
前記ターゲット領域に対する前記第１のエネルギーのビームの送出を容易にするように構成された第１の検出器システムと、
前記第２のエネルギーを有し、前記身体の前記一部分を通過する前記イオンのそれぞれについて、前記身体の前記一部分の前及び後の軌跡、並びに、前記身体の前記一部分を通過した後のエネルギーを測定するように構成された第２の検出器システムと、
少なくとも前記第２の検出器システムから信号を読み出し、それによって前記第２のエネルギーのイオンのそれぞれの前記軌跡及び前記エネルギーを表す測定データを得るように構成されたデータ収集システムと、
前記測定データを処理して画像再構成を実施し、それによって前記身体の前記一部分のコンピュータ断層画像を得るように構成されたプロセッサであって、前記画像再構成は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを含む、プロセッサと、
を備える、粒子放射線治療システム。
【請求項３６】
前記イオンは陽子を含む、請求項３５に記載のシステム。
【請求項３７】
前記イオンは炭素イオンを含む、請求項３５に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、包括的には、医療撮像の分野に関し、より詳細には、陽子コンピュータ断層撮影のためのシステム及び方法に関する。
【０００２】
［優先権の主張］
本出願は、２０１０年２月１２日に出願された「PROTON COMPUTED TOMOGRAPHY」という名称の米国仮特許出願第６１／３０３，７８３号の米国特許法第１１９条（ｅ）による優先権の利益を主張し、この出願は、参照によりその全体が本明細書に援用され、本明細書の一部とみなされる。
【０００３】
［連邦政府による資金提供を受けた研究開発の記載］
本発明は、国立心肺血液研究所によって与えられたＮＩＨ認可番号第Ｒ０１ＨＬ０７０４７２号による政府の支援によって、一部資金提供を受けたものである。政府は、本発明において一定の権利を有する。
【背景技術】
【０００４】
コンピュータ断層撮影（ＣＴ）は、３次元オブジェクトの特定の物理的特性の再構成を可能にし、この再構成を、オブジェクトの２次元断面、すなわち「断層」画像のアレイとして配列し表示する。こうした再構成は、オブジェクトに貫入する、適切に構成されたＸ線又は粒子放射線によって容易にすることができる。こうした放射線の検出及びこうしたデータの処理は、こうした２次元画像の再構成を容易にすることができる。
【０００５】
陽子ビームは、オブジェクトを貫通するように構成することができる。こうしたオブジェクトを通過する陽子からの２次元画像の再構成は、関わる陽子の数が比較的多く、オブジェクト内で陽子が受ける相互作用のタイプが様々であるため困難である場合がある。
【発明の概要】
【０００６】
幾つかの実施態様において、本開示は、コンピュータ断層撮影を実施するための方法に関する。本方法は、オブジェクトを通過する複数の陽子について測定データを取得することを含む。測定データは、陽子のそれぞれについて第１のトラック及び第２のトラックに関する情報を含む。第１のトラック及び第２のトラックは、それぞれ、陽子の、オブジェクトの通過の前及び後の陽子の軌跡に対応する。測定データは、陽子の、オブジェクトの通過によって生じる各陽子の相互作用量に関する情報を更に含む。本方法は、第１のトラック及び第２のトラックに少なくとも部分的に基づいてオブジェクト内で陽子がとる経路を推定することを更に含む。本方法は、陽子の、オブジェクトの通過が、連立方程式Ａｘ＝ｂとして表されるか又は表すことが可能であるように、陽子の相互作用量及び推定経路を配列すること更に含み、ここで、ｘは、オブジェクトに関連するパラメータの分布であり、ｂは、オブジェクト内の陽子のそれぞれの経路に沿う相互作用から生じる陽子の相互作用量を表し、Ａは、ｂを得るためにｘに作用する作用素である。作用素Ａは、オブジェクト内の陽子の推定経路に関する情報を有する。連立方程式は、複数の解を有するように構成することができる。本方法は、連立方程式について初期解を推定することを更に含む。本方法は、複数の解の中で１つ又は複数の実行可能解を探すことであって、各実行可能解は、既存の解に摂動を与えることによって得られ、既存の解の摂動なしで得られる別の解よりも、オブジェクトパラメータ分布の再構成に関連する量について優れた特性を有する、探すことを更に含む。本方法は、１つ又は複数の実行可能解のうちの選択された１つの解に基づいてオブジェクトパラメータ分布を計算することを更に含む。
【０００７】
幾つかの実施態様において、陽子の相互作用量は、陽子の、オブジェクトの通過によって生じる陽子のエネルギー損失を含むことができる。幾つかの実施態様において、ｂは、陽子の推定経路に沿う相互作用量の積分値を表すことができる。幾つかの実施態様において、連立方程式Ａｘ＝ｂは、連立一次方程式とすることができる。こうした連立方程式において、作用素Ａは行列とすることができる。
【０００８】
幾つかの実施態様において、選択された実行可能解は、複数の解の中で最適解でない実行可能解を含むことができる。幾つかの実施態様において、オブジェクトパラメータ分布を計算することは、３Ｄオブジェクトパラメータ分布を計算することを含むことができる。幾つかの実施態様において、本方法は、計算されたオブジェクトパラメータ分布に基づいてオブジェクトの断層画像のアレイを形成することを更に含むことができる。
【０００９】
幾つかの実施態様において、オブジェクトパラメータ分布は、電子密度ベース量の分布を含むことができる。電子密度ベース量は、水等の実質的に一様な物質に対する相対的陽子阻止能を含むことができる。幾つかの実施態様において、再構成されるオブジェクトパラメータ分布は、電子密度ベース量の３次元分布を含むことができる。
【００１０】
幾つかの実施態様において、経路を推定することは、陽子の最尤経路を推定することを含むことができる。幾つかの実施態様において、オブジェクトパラメータ分布の再構成に関連する量は、再構成されたオブジェクトパラメータ分布の総変動を含むことができる。総変動の優れた特性は、総変動の、より低い値を含むことができる。
【００１１】
幾つかの実施態様において、初期解を推定することは、フィルタリングされた逆投影再構成解を計算することを含むことができる。
【００１２】
幾つかの実施態様において、１つ又は複数の実行可能解を探すことは、オブジェクトパラメータ分布ｘのベクトル表現ｘｋに摂動を与えることであって、摂動付与済みベクトルｙｋを得る、摂動を与えることと、摂動付与済みベクトルｙｋに関連する再構成されたオブジェクトパラメータ分布に関連する量を評価することと、摂動付与済みベクトルｙｋに関連する量が、未摂動ベクトルｘｋに関連する量より優れている場合、ｙｋを投影することであって、それによって次のベクトル表現ｘｋ＋１を得る、ｙｋを投影することと
の反復を実施することを含む。
【００１３】
幾つかの実施態様において、ベクトルｘｋに摂動を与えることは、ｙｋ＝ｘｋ＋βｋｖｋとなるようにｙｋを計算することを含むことができ、ここで、βｋは、摂動の大きさを表し、ｖｋは、摂動ベクトルである。幾つかの実施態様において、ｙｋについて評価された量が、ｘｋについて評価された量以下である場合、摂動付与済みベクトルｙｋに関連する量は、未摂動ベクトルｘｋに関連する量に対して優れている。幾つかの実施態様において、投影することは、或る関連する凸集合上へのｘｋの数学的に規定された投影としてｘｋ＋１を計算することを含むことができる。投影することは、測定データが複数のブロックに分割されるように、ブロック反復投影アルゴリズムを使用して、超平面、半空間、超スラブ、又は他の凸集合上に投影することを含むことができる。幾つかの実施態様において、投影することは、連立方程式の個々の超平面上への直交投影とともに、対角成分ごとの緩和を可能にするように構成された対角緩和した直交投影（ＤＲＯＰ）ベースアルゴリズムを使用して投影することを含むことができる。幾つかの実施態様において、ＤＲＯＰベースアルゴリズムは、対角緩和した直交行列λｋＵｔ（ｋ）を含むことができ、ここで、λｋは、ｋ番目の反復についての緩和パラメータであり、Ｕｔ（ｋ）は、ｔ番目のブロックについて対角要素ｍｉｎ（１，１／ｈｔｊ）を有する対角行列であり、ｈｔｊは、ベクトルｘｋのｊ番目のボクセルと交差するｔ番目のブロック内の陽子履歴の数である。
【００１４】
幾つかの実施態様において、投影することは、次の反復に進む前に、ブロックが全て処理されるまで、ブロックを通して周期的に実施することができる。幾つかの実施態様において、反復は、処置される所与のブロックについて、投影することが該所与のブロックについてだけ実施されるように、各ブロックについて実施することができる。
【００１５】
幾つかの実施態様において、本開示は、陽子コンピュータ断層撮影を実施する方法に関する。本方法は、オブジェクトを通過する複数の陽子について測定データを取得することを含む。本方法は、測定データに対する総変動優秀化に基づいて投影ベース再構成アルゴリズムを反復して適用することであって、オブジェクトの相対阻止能の分布を得る、適用することを更に含む。幾つかの実施態様において、本方法は、相対阻止能分布に基づいてオブジェクトの視覚画像を形成することを更に含むことができる。
【００１６】
幾つかの実施態様において、本開示は、陽子コンピュータ断層撮影システムに関する。本システムは、オブジェクトを通過するのに十分な、選択された平均エネルギーを有する複数の陽子を送出するように構成された陽子送出システムを備える。本システムは、陽子のそれぞれについて、オブジェクトの前及び後の軌跡、並びに、オブジェクトを通過した後のエネルギーを測定するように構成された検出器システムを更に備える。本システムは、検出器システムから信号を読み出し、それによって陽子のそれぞれの軌跡及びエネルギーを表す測定データを得るように構成されたデータ収集システムを更に備える。本システムは、測定データを処理して画像再構成を実施し、それによってオブジェクトのコンピュータ断層画像を得るように構成されたプロセッサを更に備える。画像再構成は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを含む。
【００１７】
幾つかの実施態様において、本開示は、陽子治療システムに関する。本システムは、第１の平均エネルギー及び第２の平均エネルギーを有する陽子ビームを送出するように構成された陽子送出システムを備える。第１の平均エネルギーは、第１のブラッグピークが、身体の一部分の内部のターゲット領域内の場所で起こるように選択される。第２の平均エネルギーは、陽子ビームが、身体の一部分を通過するように選択される。本システムは、ターゲット領域に対する第１のエネルギーのビームの送出を容易にするように構成された第１の検出器システムを更に備える。本システムは、第２のエネルギーを有し、身体の一部分を通過する陽子のそれぞれについて、身体の一部分の前及び後の軌跡、並びに、身体の一部分を通過した後のエネルギーを測定するように構成された第２の検出器システムを更に備える。本システムは、少なくとも第２の検出器システムから信号を読み出し、それによって第２のエネルギーの陽子のそれぞれの軌跡及びエネルギーを表す測定データを得るように構成されたデータ収集システムを更に備える。本システムは、測定データを処理して画像再構成を実施し、それによって身体の一部分のコンピュータ断層画像を得るように構成されたプロセッサを更に備える。画像再構成は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを含む。
【００１８】
幾つかの実施態様において、本開示は、有形コンピュータ可読記憶媒体であって、該有形コンピュータ可読記憶媒体上に記憶されたコンピュータ実行可能命令を有し、該コンピュータ実行可能命令は、１つ又は複数のコンピューティングデバイスを有するコンピューティングシステムによって可読であり、コンピュータ実行可能命令は、複数の陽子に関するデータを取得することを含む動作を、コンピューティングシステムに実施させるように、該コンピューティングシステム上で実行可能である、有形コンピュータ可読記憶媒体に関する。データは、陽子のそれぞれについて第１のトラック及び第２のトラックに関する情報を含む。データは、第１のトラックと第２のトラックとの間での各陽子のエネルギー損失に関する情報を更に含む。この動作は、各陽子について、第１のトラックと第２のトラックとの間の経路を推定することを更に含む。この動作は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを使用して、陽子のエネルギー損失及び経路に基づいて相対阻止能分布の断層撮影解析を実施することを更に含む。
【００１９】
幾つかの実施態様において、本開示は、粒子放射線治療システムに関する。本システムは、第１の平均エネルギー及び第２の平均エネルギーを有するイオンビームを送出するように構成された粒子放射線送出システムを備える第１の平均エネルギーは、第１のブラッグピークが、身体の一部分の内部のターゲット領域内の場所で起こるように選択される。第２の平均エネルギーは、イオンビームが、身体の一部分を通過するように選択される。本システムは、ターゲット領域に対する第１のエネルギーのビームの送出を容易にするように構成された第１の検出器システムを更に備える。本システムは、第２のエネルギーを有し、身体の一部分を通過するイオンのそれぞれについて、身体の一部分の前及び後の軌跡、並びに、身体の一部分を通過した後のエネルギーを測定するように構成された第２の検出器システムを更に備える。本システムは、少なくとも第２の検出器システムから信号を読み出し、それによって第２のエネルギーのイオンのそれぞれの軌跡及びエネルギーを表す測定データを得るように構成されたデータ収集システムを更に備える。本システムは、測定データを処理して画像再構成を実施し、それによって身体の一部分のコンピュータ断層画像を得るように構成されたプロセッサを更に備える。画像再構成は、総変動優秀化に基づく反復の投影ベース再構成アルゴリズムを含む。
【００２０】
幾つかの実施態様において、イオンは陽子を含む。幾つかの実施態様において、イオンは炭素イオンを含む。
【００２１】
本開示を要約するために、本発明の或る特定の態様、利点及び新規の特徴が本明細書で述べられている。全てのこうした利点を、本発明の任意の特定の実施形態によって達成することが必ずしもできるわけではないことが理解される。そのため、本発明は、本明細書で教示又は示唆される可能性がある他の利点を必ずしも達成することなく、本明細書で教示される１つの利点又は利点の群を達成するか又は最適化するように、具現化又は実行することができる。
【図面の簡単な説明】
【００２２】
【図１Ａ】幾つかの実施態様において、陽子コンピュータ断層撮影（ｐＣＴ）システムを撮像システムとして構成することができることを概略的に示す図である。
【図１Ｂ】幾つかの実施態様において、ｐＣＴシステムを、陽子治療システムを使用して患者の処置を容易にするように構成することができることを概略的に示す図である。
【図２Ａ】光子治療及び陽子治療についての例示的な線量プロファイルを示す図であり、陽子は、ターゲットにされる所望の領域においてより選択的な線量送出を提供するように構成することができる。
【図２Ｂ】光子治療及び陽子治療についての例示的な線量プロファイルを示す図であり、陽子は、ターゲットにされる所望の領域においてより選択的な線量送出を提供するように構成することができる。
【図３Ａ】所望のターゲットに対する陽子のブラッグピークの場所の間違った計算が、ターゲットの外の領域の望ましくない照射をもたらし得ることを示す図である。
【図３Ｂ】所望のターゲットに対する陽子のブラッグピークの場所の間違った計算が、ターゲットの外の領域の望ましくない照射をもたらし得ることを示す図である。
【図４】Ｘ線コンピュータ断層撮影（ＣＴ）を使用した所望のターゲットの識別に基づいて、陽子飛程を処置のためにどのように調整することができるかを概略的に示す図である。
【図５】関心体積を通過する陽子の特徴付けに基づいて断層撮影を実施するためのプロセスを示す図であり、こうした断層撮影から得られる情報は、所望のターゲットの識別等の多数のアプリケーションで利用することができる。
【図６】図５の断層撮影プロセスを容易にするための、陽子ビーム構成及びターゲットを通過する陽子の検出の例を概略的に示す図である。
【図７】幾つかの実施態様において、ｐＣＴシステムが、ターゲットを通過する単一陽子を追跡し、該単一陽子のエネルギー損失を確定することが可能であるように構成された検出器を含むことができることを示す図であり、こうした測定は、ターゲットの少なくとも一部分の画像の断層再構成を可能にするために、ターゲットに対して複数の配向で実施することができる。
【図８】図７の例示的なｐＣＴ検出器を搭載するための例示的なプラットフォームを示す図である。
【図９】人間等の被検体のｐＣＴ撮像を実施するように構成された例示的なシステムを示す図である。
【図１０】ターゲットの上流又は下流で単一陽子を追跡するように構成された例示的な追跡検出器の部分分解図である。
【図１１】単一陽子のエネルギーを測定することが可能であるように構成された例示的な熱量計を示す図である。
【図１２】ターゲット内の多重クーロン散乱（ＭＣＳ）等の作用によって陽子の経路がどのように偏向し得るか、及び、最尤経路（ＭＬＰ）をどのように推定することができるかの例を示す図である。
【図１３】陽子経路が真っすぐであるという仮定に基づく再構成画像と陽子のＭＬＰを組込んだ再構成画像との間の比較を示す図である。
【図１４】オブジェクトを通過する陽子の測定に基づいて、オブジェクトを特徴付けするための例示的な再構成方法を概略的に示す図である。
【図１５】図１４の再構成方法に対する例示的な離散的手法を概略的に示す図である。
【図１６】図１５の離散的手法のより詳細な例を示す図である。
【図１７】図１５の離散的相互作用手法を記述するシステム行列を推定するために、陽子経路とターゲットのボクセルとの交差をどのように特徴付けることができるかの例を示す図である。
【図１８】図１５の離散的相互作用手法を記述するシステム行列を推定するために、陽子経路とターゲットのボクセルとの交差をどのように特徴付けることができるかの例を示す図である。
【図１９Ａ】シミュレートされたデータの再構成の例示的な結果を示す図である。
【図１９Ｂ】シミュレートされたデータの再構成の例示的な結果を示す図である。
【図２０Ａ】陽子エネルギー測定が実質的に雑音なしであり、また、ＭＬＰが実質的に正確である理想的な再構成状況について複数の超平面を示す図であり、各超平面は、実質的に全ての超平面が、真のオブジェクト関数を表す点で交差するように、ターゲットの離散的ボクセルにおける陽子の相互作用を表すシステム行列による投影に対して、ターゲットのオブジェクト関数の変換を特徴付ける連立一次方程式の一方程式を表す。
【図２０Ｂ】超スラブが真のオブジェクト関数を含む凸領域で交差することができるように、陽子測定が雑音なしではなく、かつ、ＭＬＰが正確でないより現実的な再構成状況について複数の超スラブを示す図である。
【図２１Ａ】様々な投影技法についての例示的な性能結果を示す図である。
【図２１Ｂ】様々な投影技法についての例示的な性能結果を示す図である。
【図２２】様々な投影技法についての例示的な性能結果を示す図である。
【図２３】様々な投影技法についての例示的な性能結果を示す図である。
【図２４】優秀化方法を利用して陽子ＣＴ投影解を得るために実装することができるプロセスを示す図である。
【図２５】ターゲットオブジェクトを特徴付けるための、優れた陽子ＣＴ解を使用するために実装することができるプロセスを示す図である。
【図２６Ａ】優秀化方法を含む、種々の再構成スキームを調査するために陽子データをシミュレートするための例示的な陽子ＣＴ構成を示す図である。
【図２６Ｂ】優秀化方法を含む、種々の再構成スキームを調査するために陽子データをシミュレートするための例示的な陽子ＣＴ構成を示す図である。
【図２７】シミュレートされた陽子データに基づく種々の再構成スキームの例示的な結果を示す図である。
【図２８Ａ】シミュレートされた陽子データに基づく種々の再構成スキームの例示的な結果を示す図である。
【図２８Ｂ】シミュレートされた陽子データに基づく種々の再構成スキームの例示的な結果を示す図である。
【図２８Ｃ】シミュレートされた陽子データに基づく種々の再構成スキームの例示的な結果を示す図である。
【図２９】シミュレートされた陽子データに基づく種々の再構成スキームの例示的な結果を示す図である。
【図３０】実験的陽子データに適用される、選択された再構成スキームを使用して得られる例示的な陽子ＣＴ画像を示す図である。
【発明を実施するための形態】
【００２３】
本明細書に設けられる表題は、たとえあっても、単に便宜上のものであり、必ずしも特許請求される本発明の範囲にも意味にも影響を及ぼさない。
【００２４】
本明細書には、相互作用する放射線として陽子を使用して、コンピュータ断層撮影（ＣＴ）を実施するための方法及び関連するシステムが述べられる。本明細書の説明は陽子に関して行われるが、本開示の１つ又は複数の特徴は、正に帯電した他のイオンを使用したＣＴアプリケーションにも同様に実装することができることも理解されるであろう。
【００２５】
図１Ａは、幾つかの実施形態において、陽子コンピュータ断層撮影（ｐＣＴ）オペレーションを実施し、オブジェクト１１０の１つ又は複数の部分のＣＴ画像として表すことができるデータを与えるように撮像システム１００を構成することができることを示す。撮像システム１００は、オブジェクト１１０に陽子ビームを送出するように構成された陽子ビームコンポーネント１０２を含むことができる。エネルギー、方向及び強度等の陽子ビームの種々のパラメータを制御することは、複数の既知の方法で達成することができる。
【００２６】
撮像システム１００は、オブジェクト１１０に入射する陽子及びオブジェクトを通過した陽子を特徴付けるように構成された検出器コンポーネント１０４を更に含むことができる。幾つかの実施態様では、こうした検出器コンポーネント１０４は、単一陽子を特徴付けることが可能であるように構成することができる。陽子のこうした特徴付けを容易にすることができるデバイスの例は、本明細書でより詳細に述べられる。
【００２７】
撮像システム１００は、ＣＴ解析を容易にするために、検出器コンポーネント１０４から信号を読出すように構成されたデータ収集（ＤＡＱ）コンポーネント１０６を更に含むことができる。ＤＡＱコンポーネント１０６によって実施される信号処理量は、変動し得る。
【００２８】
幾つかの実施態様では、種々の検出器からの信号は、１つ又は複数のアナログ／デジタル変換器（ＡＤＣ）によってデジタル信号に変換することができ、こうしたデジタル信号は、制御コンポーネント１１２の制御下で読み出すことができる。イベントトリガー、イベント信号のタイミング及び読出し、並びに検出器のリセット等の種々の制御パラメータも、制御コンポーネント１１２が制御することができる。
【００２９】
幾つかの実施態様では、撮像システム１００は、デジタル化信号を受信し、オブジェクト１１０の上流及び下流での陽子の追跡、並びに、オブジェクト１１０を通過した下流の陽子のエネルギーの計算等の解析を実施するように構成されるプロセッサ１０８を更に含むことができる。幾つかの実施態様では、断層再構成処理もまた、プロセッサ１０８によって実施することができる。他の実施態様では、こうした断層再構成処理は、別個のプロセッサによって実施することができる。
【００３０】
幾つかの実施態様では、撮像システム１００は、システム１００の１つ又は複数のコンポーネントのオペレーションを容易にする情報及び／又は実行可能命令を記憶するように構成されたコンピュータ可読媒体１１４を更に含むことができる。幾つかの実施態様では、コンピュータ可読媒体１１４は、本明細書で述べる１つ又は複数の再構成プロセスの実施を容易にする情報及び／又は実行可能命令を含むことができる。幾つかの実施態様では、こうした情報及び／又は実行可能命令は、一時的でない方法で記憶することができる。
【００３１】
幾つかの実施態様では、本開示の１つ又は複数の特徴は、陽子線治療システム又は炭素ビーム治療システム等の放射線治療システム１２０に組込むことができる。治療システム１２０は、患者１３０に陽子線又は炭素イオンビームを送出するように構成された陽子線コンポーネント又は炭素ビームコンポーネント１２２を含むことができる。陽子又は炭素イオンのこうしたビームは、患者に対する治療効果をもたらすように構成することができる。或る特定の実施態様では、陽子ビームコンポーネント１２２は、図１Ａを参照して上述したように、断層解析を可能にするために、患者を通過することができる陽子ビームを与えるようにも構成することができる。こうしたビームを提供することができる方法の例は、本明細書で更に詳細に述べられる。
【００３２】
撮像システム１２０は、陽子ビームの処置利用を容易にするように構成された検出器コンポーネント１２４を更に含むことができる。こうした検出器コンポーネント１２４は、エネルギー、方向及び強度等の所望のパラメータを用いて患者１３０に入射する陽子を特徴付けるように構成されるデバイスを含むことができる。こうしたデバイスは、複数の既知の方法で実装することができる。
【００３３】
幾つかの実施態様では、検出器コンポーネント１２４は、図１Ａを参照して述べたｐＣＴ撮像機能のようなｐＣＴ撮像機能を容易にするように構成されるデバイスを更に含むことができる。幾つかの実施形態では、治療に関連する検出デバイスのうちの少なくとも幾つかは、ｐＣＴ撮像のためにも利用することができる。例えば、患者の上流のビーム検出器は、撮像モードのオペレーション中に患者に入射する個々の陽子を特徴付けるために利用することができる。
【００３４】
撮像システム１２０は、治療モード及び／又は撮像モードのオペレーションを容易にするように構成されたデータ収集（ＤＡＱ）コンポーネント１２６と、制御コンポーネント１３２と、プロセッサコンポーネント１２８と、コンピュータ可読媒体コンポーネント１３４とを更に含むことができる。
【００３５】
適切に構成された陽子又は炭素イオン放射線（すなわち、帯電した粒子放射線）は、或る特定のがん処置等の治療アプリケーションにおいて多くの利益を提供することができる。こうした利益のうちの１つの利益は、所与の物質内での陽子の移動の終りにおける急激なエネルギー損失による可能性がある。こうした急激なエネルギー損失は、ブラッグピークと呼ばれる比較的急峻なピークを有し、（同じビームエネルギーを有する）粒子のうちのほんの僅かの粒子が、こうした深さを超えて貫入する。こうしたブラッグピークの深さの位置は、特定のビームエネルギーに依存し得る。一般に、より深いブラッグピークは、より高いエネルギー粒子ビームによって達成することができる。治療に使用される陽子は、約７０ＭｅＶ／原子質量単位〜２５０ＭｅＶ／原子質量単位の範囲のエネルギーを有することができ、炭素イオンは、４３０ＭｅＶ／原子質量単位までの範囲のエネルギーを有することができる。
【００３６】
図２Ａは、エネルギー陽子が所与の物質（例えば、組織）内で移動するときの、深さの関数としてのエネルギー損失プロファイル１４２のブラッグピーク１４４の例を示す。比較すると、電磁放射線（例えば、Ｘ線又はガンマ線）についての相対的線量プロファイル１４０は、最大まで比較的急峻に上昇し、それに続いて、深さの関数として徐々に減少している。したがって、陽子ベース放射線は、陽子と炭素イオンの使用によって提供される同様の終了範囲制御を提供しない。
【００３７】
図２Ｂは、複数の初期（pristine）ブラッグピークが結合されて、累積線量プロファイル１５２の拡大ブラッグピーク（spread out Brag Peak）（ＳＯＢＰ）１５４をもたらすことができることを示す。こうした初期ブラッグピークは、同じ体積を、異なるエネルギーを有する陽子ビームに晒すことによって達成することができる。結果として得られる拡大ブラッグピーク１５４の位置は、ターゲット領域１５０の深さ境界にオーバーラップするように選択することができる。ビームエネルギーが適切に選択される場合、拡大ブラッグピークは、ターゲット領域の遠位境界を超えると急激に降下し得る。
【００３８】
上記に基づいて、ターゲット領域の深さ境界と、拡大ブラッグピークとを適切に一致させることは、粒子治療の重要な考慮事項であり得る。拡大ブラッグピークの遠位部分が、図３Ａの例の場合のように深すぎる場合、不必要かつ有害な放射線線量（例えば、ブラッグピーク１６４のかなりの部分）が、ターゲット領域１６０の遠位境界を越えた領域に提供される。拡大ブラッグピークの近位部分が、図３Ｂの例の場合のように浅すぎる場合、不必要で余分な放射線線量（例えば、ブラッグピーク１６２のかなりの部分）が、ターゲット領域１６０の近位境界の前の領域に提供される。他方で、深すぎる拡大ブラッグピークの近位部分及び／又は浅すぎる拡大ブラッグピークの遠位部分は、ターゲット領域の或る特定の部分が適切に照射されないことをもたらす場合がある。
【００３９】
所与の媒体内での陽子貫入深さの上記の不確実性は、媒体と異なるように相互作用する異なるプローブ放射線による媒体の特徴付けからもたらされる可能性がある。図４は、所与のエネルギーを有し、所与の媒体を横切る陽子ビームの飛程を求めるために使用することができる一般的なプロセスを示す。同様のプロセスは、炭素イオンビームについて適用することができ、その差は、陽子の相対阻止能の代わりに、炭素イオンの相対阻止能（陽子についての陽子阻止能に原子数比の２乗（ＺＣ／Ｚｐ）２＝３６を乗算することによって推定することができる）が使用されることである。オブジェクトの減衰係数測定に基づいて画像を取得するために、Ｘ線ＣＴをオブジェクト（例えば、患者）に対して実施することができる（ブロック１５０）。こうした測定は、ハウンスフィールド単位スケールに変換することができる（ブロック１５２）。媒体内での陽子又は炭素イオンの、計算された相対阻止能は、ハウンスフィールドスケール測定に基づいて取得することができる（ブロック１５４）。陽子又は炭素イオンの飛程は、陽子又は炭素イオンについて、ビームエネルギー及び物質の相対阻止能の関数として表すことができるため、計算された陽子又は炭素イオンの飛程を取得することができる（ブロック１５６）。
【００４０】
Ｘ線ＣＴ測定から陽子飛程を計算する上記のプロセスは、陽子の飛程又は炭素イオンの飛程の約３．５％の不確実性をもたらす可能性がある。こうした相対的な不確実性は、人間の身体の様々な部分についての様々な距離に等しいとすることができる。例えば、３．５％の相対的不確実性は、脳で約３ｍｍ〜５ｍｍに等しく、骨盤で約１０ｍｍ〜１２ｍｍに等しい。更なる不確実性が、未知の密度を有する物質の存在及びＸ線ＣＴ画像のストリーク状アーティファクトによって導入され得る。
【００４１】
幾つかの実施態様では、本明細書に述べるようなｐＣＴ技法に基づく陽子又は炭素イオンの飛程の確定によって、飛程の不確実性を陽子の飛程の約１％以下に低減することができる。図５は、陽子コンピュータ断層撮影を実施するために実装することができるプロセス１７０を示す。図６は、図５のプロセス１７０のために実装することができる例示的な構成を示す。
【００４２】
ブロック１７２において、陽子ビームを、関心体積を有するオブジェクトに提供することができる。幾つかの実施態様では、ビームは、ビーム内の陽子のかなりの部分がオブジェクトを通過するように構成することができる。例えば、また、図６に示すように、陽子ビームには、結果として得られるブラッグピーク１８２が関心体積の遠位境界を過ぎた所で起こるようにエネルギーを提供することができる。幾つかの実施態様では、撮像されるオブジェクトは、関心体積を画定することができる。こうした状況では、陽子は、オブジェクト全体を通過し、オブジェクトの下流で検出することができる。
【００４３】
ブロック１７４において、体積を通過する複数の陽子のそれぞれを追跡することができる。例えば、図６に示すように、こうした追跡は、関心体積１８０の上流及び下流にある陽子の追跡を含むことができる。こうした追跡の例は、本明細書でより詳細に述べられる。
【００４４】
ブロック１７６において、追跡された各陽子のエネルギー損失を測定することができる。例えば、図６に示すように、こうしたエネルギー測定は、体積１８０の上流にある陽子についての平均ビームエネルギーの測定と、熱量計等のエネルギー検出器による、下流の各陽子の測定とを含むことができる。こうしたエネルギー測定の例は、本明細書でより詳細に述べられる。
【００４５】
ブロック１７８において、コンピュータ断層撮影解析を、陽子の追跡情報及びエネルギー損失情報に基づいて実施することができる。例えば、図６に示すように、こうしたＣＴ解析は、複数の異なる向きで体積１８０に入射し（１８４）、体積１８０を通過する（１８６）陽子に関するこうした情報を測定することによって容易にすることができる。こうした向きの変化の例は本明細書でより詳細に述べられる。
【００４６】
図７は、幾つかの実施形態において、検出器の組立体を、オブジェクト（例えば、患者の一部分）２００に入射する個々の陽子及びオブジェクト２００を通過した後の個々の陽子の測定を可能にするように構成することができることを示す。例えば、単一陽子１９２が、オブジェクト２００に入射し（１９４）、オブジェクト２００を横切り（１９６）、オブジェクト２００を通過した（１９８）として示される。
【００４７】
幾つかの実施態様では、上流検出器２０２は、陽子１９２について入射トラック（１９４）情報を提供することができ、下流検出器２０４は、陽子１９２について下流トラック（１９８）情報を提供することができる。本明細書で述べるように、陽子が、オブジェクト２００を横切るときにとる経路１９６は、ＣＴ解析を容易にするために推定することができる。
【００４８】
上流追跡検出器２０２及び下流追跡検出器２０４は、複数の方法で実装することができる。上流追跡検出器２０２は、陽子がオブジェクト２００に入るときの陽子の空間ベクトル及び場所の確定を可能にするように構成することができる。同様に、下流追跡検出器２０４は、陽子がオブジェクト２００を出るときの陽子の空間ベクトル及び場所の確定を可能にするように構成することができる。こうした情報に基づいて、横断経路１９６を、本明細書で述べるように推定することができる。
【００４９】
幾つかの実施態様では、上流追跡検出器２０２及び下流追跡検出器２０４のそれぞれは２つ以上の２Ｄ位置感応性検出面を含むことができる。例えば、上流検出器２０２は、２つの２Ｄ位置感応性検出面２０２ａ、２０２ｂを含むことができ、下流検出器２０４は、２つの２Ｄ位置感応性検出面２０４ａ、２０４ｂを含むことができる。２Ｄ位置感応性検出面２０２ａ、２０２ｂ、２０４ａ、２０４ｂのそれぞれは、例えばマイクロストリップ検出器（例えば、シリコンストリップ検出器）を含む、複数の方法で実装されたＸ位置感応性面及びＹ位置感応性面を含むことができる。（例えば、高計数環境で）追跡が不明瞭である可能性がある幾つかの状況では、更なる位置感応性面（例えば、Ｕ及び／又はＶ面）を付加することができるか、又は、より高速な検出器及び／又は読出しシステムを、こうした不明瞭さを低減するために使用することができる。
【００５０】
図７を参照すると、検出器の組立体は、下流の陽子のエネルギーを検出するように構成されたエネルギー検出器２０６を更に含むことができる。こうしたエネルギー検出器は、例えば、多結晶ベース熱量計を含む、複数の方法で実装することができる。こうしたエネルギー検出器の例は、本明細書でより詳細に述べられる。
【００５１】
或る特定の実施態様では、複数の他の検出器及び／又はデバイスを、種々の機能を提供するために含めることができる。例えば、２重粒子又は入射ストレイ（stray）粒子を含むイベントを識別し拒否し得るビーム規定カウンターを含めることができる。別の例では、偽信号を拒絶しながら、良好なイベントを規定しトリガーするように構成されたカウンター及び回路も含めることができる。
【００５２】
図７を参照すると、検出器のこうした組立体及び／又はオブジェクト２００は、ＣＴ解析のために複数の向きで陽子データを取得することを容易にするために、選択された軸の回りに回転させることができる。例えば、検出器の組立体は、オブジェクト２００の軸２１２の回りに一緒に回転する（矢印２１０）ように構成することができる。
【００５３】
他の例では、検出器、オブジェクト、又はそれらの何らかの組合せが、ＣＴ機能を容易にするために、回転及び／又は並進することができることが理解されるであろう。例えば、図８は、検出器（例えば、２Ｄ位置感応性検出器２３０、２３２、２３４、２３６及びエネルギー検出器２３８）の組立体を、実質的に固定してプラットフォーム２２６に搭載し、陽子ビーム２５０を受信して下流の陽子を検出するように向けることができる構成２２０を示す。特徴付けされるオブジェクト２２２は、ＣＴ機能を提供するために、軸の回りのオブジェクト２２２の回転を可能にするように、搭載構造２２４に搭載することができる。オブジェクトが回転するこうした構成は、オブジェクト２２２が無生物オブジェクト（例えば、ファントム）であり、かつ、陽子ビームが実質的に固定方向を有する場合に、より適切であり得る。
【００５４】
別の例では、図９は、オブジェクト２６２（例えば、患者）が実質的に静止しており、かつ、検出器（例えば、上流追跡検出器２７２、下流追跡検出器２７４、及びエネルギー検出器２７６）の組立体を、オブジェクト２６２によって規定される軸２８０の回りに回転するように構成することができる構成２６０を示す。
【００５５】
例示的な検出器２７２、２７４、２７６は、軸２８０の回りの検出器組立体の回転（矢印２８２）を可能にする回転可能構造２７０に搭載することができる。静止オブジェクト２６２は、例えば、安定した構造２６６に搭載された支持構造２６４によって達成することができる。幾つかの実施形態では、支持構造２６４は、陽子２９０が上流検出器２７２と下流検出器２７４との間で移動する領域を出入りするオブジェクトの移動を可能にするように構成することができる。検出器組立体がオブジェクトの回りに回転するこうした構成は、オブジェクト２６２が、人間の患者等の生き物である場合により適切であり得る。
【００５６】
幾つかの実施態様では、図９の例示的な構成２６０の検出器の一部又は全ては、処置ガントリ設計を有する陽子治療システムのような陽子治療システムに組込むことができる。こうした設計では、患者を、ビームラインの内外に移動させることができる。ビームラインは、処置陽子の所望の送出のために患者の回りに回転することができる。
【００５７】
こうしたシステムの場合、上流検出器と、下流検出器と、エネルギー検出器とは、ビームラインとともに回転するように回転可能構造に搭載することができる。こうした検出器は、所定場所に固定することができるか、又は、治療オペレーションを妨げないように、格納可能とすることができる。幾つかの実施形態では、治療システムのビーム制御デバイスは、上流追跡検出器に関連する機能の一部又は全てを提供することができる。
【００５８】
図１０は、図７〜図９を参照して述べた上流追跡検出器又は下流追跡検出器として利用することができる例示的なトラッカモジュール３００を示す。本明細書で開示される多数のｐＣＴ結果が、こうしたデバイスを使用して得られた。図示するように、トラッカモジュール３００は、２つ以上の２Ｄ位置感応性検出器を保持するように構成されたハウジング３０２を含むことができる。こうした２Ｄ位置感応性検出器３０６のうちの１つが、部分的に除去されて示され、他の２Ｄ位置感応性検出器は、図１０では取除かれている。
【００５９】
２Ｄ位置感応性検出器３０６は、アクティブエリア３０８（例えば、シリコンストリップ検出素子を有する長方形）と、検出素子を制御し、検出素子から信号を読出すための回路３１０と、電力、制御信号、読出し信号等を供給するためのインタフェース３１２とを含むことができる。
【００６０】
ハウジング３０２は、相互作用の確率が比較的低い状態での陽子の通過を可能にするために、（設置すると）２Ｄ位置感応性検出器３０６のアクティブエリア３０８に好ましくは整列する窓３０４を含むことができる。必要であるとき、こうした窓は、マイラー等の物質の薄いシートを含むことができる。
【００６１】
ハウジング３０２は、モジュール３００が、（例えば、オブジェクトが回転可能である場合）静止状態になるか又はオブジェクト（図示せず）の回りに回転可能であることを可能にするために、モジュール３００を支持構造に搭載することを可能にするように構成することができる。
【００６２】
図１１は、図７〜図９を参照して述べた例で利用することができる例示的なエネルギー検出器３２０を示す。本明細書で開示される複数のｐＣＴ結果が、こうしたデバイスを使用して得られた。図示するように、エネルギー検出器３２０は、分割型（segmented）エネルギー検出能力を提供するように配列された複数の検出器結晶３２２（例えば、ＣｓＩ結晶）を含むことができる。図示の例では、結晶３２２は、グリッド様分割を提供するために、グリッドアレイで詰めることができる。結晶３２２の前面の寸法及び結晶３２２の深さ等の設計パラメータは、下流の陽子の予想されるエネルギー範囲及び強度に対応するように選択することができる。
【００６３】
図１１を参照すると、結晶３２２の後端が、光検出器３２６に結合されているものとして示される。こうした光検出器は、例えば、種々のタイプの光電子増倍管の使用を含む多数の方法で実装することができる。例えば、下流の陽子の飛程を検出するように構成された、固体検出器又は有機シンチレータ若しくは無機シンチレータの積重体を含む多数の他のタイプのエネルギー検出器も利用することができる。
【００６４】
図１１では、結晶３２２は、ハウジング３２４内に詰められているものとして示される。こうしたハウジングは、エネルギー検出器３２０が、（例えば、オブジェクトが回転可能である場合）静止状態になるか又はオブジェクト（図示せず）の回りに回転可能であることを可能にするために、エネルギー検出器３２０を支持構造に搭載することを可能にするように構成することができる。
【００６５】
本明細書で述べるように、特徴付けされるオブジェクト内での陽子の相互作用は通常、入射方向からの方向変化及び入射エネルギーからのエネルギー損失をそれぞれもたらす。陽子は、オブジェクトを通過するとき、ターゲット原子の核ポテンシャルにおける散乱によって、複数回の小角偏向を受ける可能性があり、こうした偏向は、元の方向からの、実質的にランダムな巨視的逸脱をもたらす可能性がある。幾つかの状況では、こうした逸脱は、数度程度の大きさである可能性があり、オブジェクトの遠位端上で（入口点及び入口方向の直線投影に対して）出口点の著しい（例えば、最大数ミリメートルの）変位をもたらす。
【００６６】
図１２は、陽子が、入射方向４００及び入口点４０２でオブジェクト４１０を横切り、出射方向４０４及び出口点４０６で出る際の、こうした小角偏向の例の２次元表現を示す。この例では、出射方向は、入る方向から逸脱しており、同様に、出口点４０６は、陽子が、直線で移動する場合に到達することになる投影点４０８からシフトしている。
【００６７】
こうした多重クーロン散乱（ＭＣＳ）作用は、陽子の軌跡Ｌに不確実性を生じる可能性がある。こうした軌跡は、陽子ＣＴ画像の空間分解能に寄与する可能性があるため、軌跡の正確な推定が望ましい。幾つかのターゲットオブジェクトにおいて、オブジェクト内での陽子の実際の経路に関する直接的な情報は存在しない。しかし、オブジェクトの前及び後の陽子の経路について知識が与えられると、オブジェクト内での最尤経路を推定するための多数の技法が存在する。例えば、準解析的計算を、複数散乱のガウス近似に関して実施することができる。別の例では、モンテカルロシミュレーションのセットを、オブジェクトのサイズ、オブジェクトの物質、ビームエネルギー、並びに検出される陽子経路及びエネルギー等の因子を含むことができる所与の構成に基づいて実施することができる。
【００６８】
図１２では、入射経路４００及び出射経路４０４に対応する最尤経路（ＭＬＰ）は、曲線４２０として示され、こうした軌跡は、ＣＴ再構成に使用することができる。
【００６９】
幾つかの実施態様では、多重クーロン散乱等の作用から生じる上記の経路逸脱は、一部の画像再構成アルゴリズムの効果を減少させる可能性がある。例えば、フィルタ補正逆投影（filtered back projection）（ＦＢＰ）は、Ｘ線ＣＴシステムで使用される一般的な再構成方法であり、特徴付けられるオブジェクトの内部で、光子が直線上に進むと仮定する。
【００７０】
こうした方法が、陽子ＣＴ画像再構成に利用されるとき、画像品質が、大幅に悪化する。図１３は、陽子ＣＴ画像を示し、ここで、上側画像は、ＦＢＰ法を使用して得られ、下側画像は、ＭＬＰを組込むことによって得られた。一般に、ＦＢＰ画像は、ＭＬＰ画像より品質が低いことがわかる。例えば、ＦＢＰ画像のかなりの部分がぼけている。
【００７１】
陽子は、その軌跡Ｌに沿って媒体を横切るにつれ、媒体との相互作用によるエネルギー損失を被る。以下の積分は、このエネルギー損失を含む。
【００７２】
【数１】
2013519452.tif 000003
【００７３】
ここで、Ｅｉｎ及びＥｏｕｔは、それぞれ、入射エネルギー及び出射エネルギーを表し、Ｓｗａｔｅｒ（Ｅ）は、水中でのエネルギー依存の陽子阻止能関数を表す。さらに、こうした積分は、
【００７４】
【数２】
2013519452.tif 000004
【００７５】
となるように、陽子経路に沿う、相対電子密度又は相対陽子阻止能の積分に実質的に等しくなるように設定することができる。幾つかの実施態様では、陽子経路は、最尤経路（ＭＬＰ）として近似することができる。
【００７６】
所与の被検出陽子の場合、方程式１の左辺は、陽子の入射エネルギー及び出射エネルギーが測定されるため、求めることができる。そのため、陽子のＭＬＰを知ることによって、オブジェクトの電子密度に密接に関連する相対阻止能ηｅをＭＬＰに沿って計算することが可能になる。十分な数のこうした陽子が与えられると、相対阻止能関数ηｅを、特徴付けされるオブジェクト内の関心領域について計算することができる。こうした関数ηｅは、高品質画像の生成を可能にするために、十分に詳細に求めることができる。
【００７７】
図１４は、幾つかの実施態様において、出力陽子５０４を得るために、オブジェクト５００を横切り、オブジェクト５００と相互作用する複数の入力陽子５０２を変換５１０  Ｍ・ｘ＝ｇとして表すことができることを示し、ここで、ｘは、オブジェクトに関連する相互作用パラメータ（例えば、電子密度又は相対阻止能（ＲＳＰ）等の電子密度ベース量）の３次元分布であり、Ｍは、ｇを得るためにｘに作用する作用素であり、ｇは、測定データを表す。そのため、オブジェクト関数ｘは、ｘ＝Ｍ−１・ｇとして求めることができ、ここで、Ｍ−１はＭの逆作用素である。本明細書で述べるように、作用素が、例えば正方行列であるとき、作用素の逆が適切であり得る。こうした正方行列が要件でないことが理解されるであろう。作用素行列が、より一般的な長方形行列である状況では、オブジェクト関数ｘは、ｘ＝Ｍ†・ｇとして求めることができ、Ｍ†（「Ｍダガー」）は、Ｍの共役転置である。同様な表記を本明細書で述べる他の作用素（例えば、図１５の作用素Ａ）に適用することができることが理解されるであろう。
【００７８】
図１５は、幾つかの実施態様において、図１４の変換の一般的な形態を、オブジェクト５２０がベクトルにデジタル化される離散的変換５２２として近似することができることを示す。こうした変換は、Ａ・ｘ＝ｂとして表すことができ、ここで、ｘは、離散的なオブジェクトベクトルであり、Ａは、ｂを得るためにｘに作用する作用素であり、ｂは、測定データのベクトルを表す。そのため、オブジェクトベクトルｘは、ｘ＝Ａ−１・ｂとして求めることができ、ここで、Ａ−１はＡの逆作用素である。
【００７９】
幾つかの実施態様では、図１５の離散的変換５２２は、オブジェクトベクトルｘが、ｎ個の要素ｘ１、ｘ２、ｘ３、．．．、ｘｎを有する列ベクトルであり、測定ベクトルｂが、ｍ個の要素ｂ１、ｂ２、ｂ３、．．．、ｂｍを有する列ベクトルであり、変換作用素Ａが、ｍ×ｎ行列であるような、連立一次方程式を表すことができる。
【００８０】
図１６は、オブジェクトベクトルｘとして表すことができるデジタル化されたオブジェクト５２０の例を示す。オブジェクト５２０は、複数のボクセル５３０を有するものとして示される。ｉ番目の陽子５３２は、その経路（例えばＭＬＰ）に沿って複数のボクセル５３０に交差するように、デジタル化されたオブジェクト５２０を横切るものとして示される。陽子５３２は、測定ベクトルｂのｉ番目の値を得るために、オブジェクト５２０を出て、検出器によって測定されるものとして更に示される。幾つかの実施態様では、変換行列Ａの要素ａｉｊは、ｉ番目の陽子とオブジェクトベクトルｘのｊ番目のボクセルとの交差長等の相互作用パラメータを表すことができる。
【００８１】
幾つかの実施態様では、変換行列Ａは、交差長（ａｉｊｓ）の推定を容易にするために、複数の方法で構成することができる。例えば、最尤経路（ＭＬＰ）は、ボクセルサイズに対するステップサイズ（例えば、半ボクセルステップ等のボクセル以下のステップサイズ）で離散化することができる。
【００８２】
幾つかの実施態様では、交差長ａｉｊは、複数の方法で推定することができる。例えば、ＭＬＰと交差する任意のボクセルに、交差長について一定値を割り当てることができる。別の例では、所与のボクセルと交差するＭＬＰのコードの長さを計算し、そのボクセルに割り当てることができる。更に別の例では、平均有効コード長係数に、ＭＬＰに沿うボクセルの全てに割り当てられた一定値を乗算することができ、こうした係数は、オブジェクトグリッド上での陽子の入射角度に基づくことができる。
【００８３】
図１７は、上記の平均有効コード長係数の例示的なプロット５４０を示し、示す入射角度は、水平入射がゼロの入射角度を有するように配向されたオブジェクトグリッド５２０（図１８）を基準とする。こうした例では、水平交差は、１の係数を与えられ、４５度交差は、最小値を与えられる。こうした割当ては、角度の付いた交差（例えば、矢印５４２によって示される区画）が、たとえ全対角長がより長くても、横寸法（例えば、正方ボクセルの水平寸法）より短い可能性が高いことを考慮する。幾つかの実施態様では、上記の単一パラメータ（入射角度）に基づく調整は、望ましいことには、処理時間を減少させることができる。
【００８４】
図１９Ａ及び図１９Ｂは、上記の交差長推定の使用から得ることができる画像品質及び阻止能特徴付けの性能の例を示す。図１９Ａでは、左上画像は、２次元陽子ビームを使用したモンテカルロシミュレーションから得られるファントムの再構成画像である。右上画像は、一定相互作用長手法を使用したファントムの再構成画像である。左下は、有効平均コード長係数手法についてのものであり、右下は、ボクセルごとに計算したコード長手法についてのものである。
【００８５】
一定相互作用長手法から得られる画像は、有効平均コード長係数手法及びボクセルごとに計算したコード長手法を使用して得られた画像に比べて、比較的雑音が多いことが示される。画像品質は別として、また、図１９Ｂに示すように、有効平均コード長係数手法及びボクセルごとに計算したコード長手法は、一定相互作用長手法の再現に比べて、ファントム内の種々の特徴部の相対阻止能のより正確な再現をもたらす。
【００８６】
幾つかの実施態様では、図１５及び図１６を参照して本明細書で述べる連立一次方程式（例えば、Ａ・ｘ＝ｂ）は、次のように表すことができる。ベクトルｂは、ｂが、ｍ個の要素ｂｉ（ｉ＝１〜ｍ）を有するように、ｍ個の陽子について測定（例えば、エネルギー損失）ベクトルを表すことができる。ベクトルｘは、特徴付けされるオブジェクトのパラメータ（例えば、相対阻止能、ＲＳＰ）の値を有する未知の画像ベクトルとすることができる。そのため、そのＭＬＰに沿ってオブジェクトを横切り、測定されるエネルギー損失をもたらすｉ番目の陽子に関して、ベクトルｂのｉ番目の要素、すなわちｂｉは、ＭＬＰに沿う積分されたＲＳＰを表す。こうした量は、例えば既知の厚さ及び相対阻止能の物質のブロックを横切った陽子に対するエネルギー検出器の応答を較正することによって測定することもできるし、
【００８７】
【数３】
2013519452.tif 000005
【００８８】
になるように、方程式１の左辺によって数値的に計算することもできる。ここで、Ｓｗａｔｅｒは、
【００８９】
【数４】
2013519452.tif 000006
【００９０】
によって与えられる水中の陽子阻止能を表す。方程式３では、ｒｅは、古典電子半径であり、ｍｅは電子質量であり、ηｗａｔｅｒ及びＩｗａｔｅｒは、それぞれ水の電子密度及び平均イオン化ポテンシャルであり、βは光の速度ｃに対する陽子の速度である。
【００９１】
幾つかの実施態様では、画像ベクトルｘは、反復再構成技法を利用することによって解くことができるか、又は推定することができる。説明のために、集合Ｉが、Ｉ＝｛１，２，．．．，ｍ｝であるような、ｍ個の指数を有すると仮定し、集合｛Ｈｉ：ｉはＩの要素｝が、ｎ次元実数座標空間Ｒｎ内の超平面の有限族であるとする。反復プロセス中にベクトルｘｋが投影される集合Ｈｉは、
【００９２】
【数５】
2013519452.tif 000007
【００９３】
であるような、ｍ×ｎ線形系Ａｘ＝ｂのｉ番目の行によって規定することができる。方程式４では、ａｉは、ＡＴ（Ａの転置）のｉ番目の列ベクトルであり、それにより、その成分は、Ａのｉ番目の行を占める。
【００９４】
幾つかの実施態様では、方程式４のｉ番目の超平面Ｈｉを満たすオブジェクト関数Ｓｉは、コスト又はメリット関数
【００９５】
【数６】
2013519452.tif 000008
【００９６】
がゼロ以下である制約を有する、ｘ（Ｒｎの要素（複数の場合もあり））の集合として表すことができる。ここで、εｉは、超スラブの半値幅を規定する閾値パラメータである。全体のオブジェクト関数Ｓは、その後、Ｓ＝∩Ｓｉ（ｉ＝１〜ｍ）であるようなｍ個のＳｉ関数の交差として求めることができる。
【００９７】
陽子エネルギー測定が実質的に雑音なしであり、かつ、ＭＬＰが実質的に正確である理想的な状況では、こうしたオブジェクト関数Ｓは、超平面の全てが交差する単一点として表すことができる。こうした理想的な状況が図２０Ａに示され、図２０Ａでは、例示的な超平面６００が、「Ｘ」として示す点で交差することが示される。こうした理想的な状況では、閾値パラメータεｉは、実質的にゼロとすることができ、＜ａｉ，ｘ＞＝ｂｉ（ｇｉ（ｘ）＝０）である実質的に正確な解を得ることができる。
【００９８】
より現実的な状況では、陽子エネルギー測定は、雑音なしの状態ではなく、かつ、ＭＬＰは正確でない。こうした状況では、超平面は、平面の代わりに、スラブとして示され、それにより、こうした超スラブの交差は、超スラブが全てオーバーラップする領域をもたらすことができる。こうした状況は図２０Ｂに示され、図２０Ｂでは、例示的な超スラブ６１０が、「Ｑ」として示す領域で交差することが示される。こうした状況では、閾値パラメータεｉは、ゼロより大きく、かつ、方程式５のメリット関数を満たす許容可能な解の外側境界を規定する量とすることができる。
【００９９】
図２０Ｂの例では、撮像される媒体を特徴付ける所望のオブジェクト関数Ｓは、オーバーラップ領域Ｑ内に存在することができる。より一般的には、撮像される媒体を特徴付ける所望のオブジェクト関数Ｓは、制約された空間ε＝｛εｉ｝内に存在することができる。
【０１００】
幾つかの実施態様では、こうした所望のオブジェクト関数は、最小化等の最適化技法によって得ることができる。最適化は、一部の撮像アプリケーションで一定の利点を有し得るが、欠点も有し得る。例えば、メリット関数によって決定される最適解は、必ずしも対象となる真のオブジェクトデータを最も良く再現する解ではない場合がある。こうした作用は、例えば、採集されたデータの不整合性及び／又はメリット関数の選択によって生じる可能性がある。さらに、メリット関数のこうした選択は、不適切であることがある引数によって、及び／又は、結果として得られる最適化問題を計算的に処理できること若しくはできないことによって影響を受ける可能性がある。幾つかの実施態様では、最適化手法は、メモリ及び処理パワー等の計算処理リソースの制限によって非実用的である可能性もある。
【０１０１】
幾つかの実施態様では、陽子ＣＴ画像再構成は、実行可能性探求法のうちの１つ又は複数を利用して実施することができる。凸集合の有限族の交差（例えば、図２０Ｂの超スラブ６１０の集合から得られる例示的な交差領域Ｑ）をもたらす撮像状況の場合、複数の解が存在し得る。凸実行可能性問題（ＣＦＰ）法では、こうした複数の解の中の解が探索されて得られ、こうした解が、交差領域内の点に対応することができる。
【０１０２】
ＣＦＰ法に基づく画像再構成（例えば、ｐＣＴ再構成）の場合、所望の性能は、効率的な実行可能性探求投影法の使用、及び、所与のメリット関数の減少した値を有する実行可能解の発見等の因子に基づくことができる。幾つかの実施態様では、優秀化（superiorization）は、対応するＣＦＰ集合の実行可能解でもある、或るメリット関数に関する優れた解を発見するプロセスを指すことができる。優れた解は、１がそれに対して優れているメリット関数の値が（the value of the merit function, with respect to which one superiorizes）、優秀化プロセスが適用されていない場合に達するであろう実行可能点で、この関数の値より小さい（必ずしも最小ではないが）、ＣＦＰの実行可能解とすることができる。
【０１０３】
幾つかの実施態様では、効率的な実行可能性探求能力を提供する上記の投影方法は、制限がある摂動弾力性（perturbation resilience）を含む１つ又は複数の因子に基づいて選択することができる。実行可能点への収束性を失うことなく、所与の投影アルゴリズムに摂動を与えることができることは、メリット関数に関して、摂動なしの状態で到達するであろう別の実行可能点より優れている実行可能点に向かってのアルゴリズムの操向を可能にすることができる。
【０１０４】
いかなる特定の理論によって拘束されることを所望することも意図することもなく、アルゴリズムＰは、以下が満たされる場合、制限がある摂動に対して弾力性があると言うことができる。シーケンス
【０１０５】
【数７】
2013519452.tif 000009
【０１０６】
（ｘから始めて、Ｐの順次的な反復適用によって得られる）が、ｎ次元実数座標空間Ｒｎ内の全てのｘについて問題Ｑの解に収束する場合、Ｒｎ内の点の任意のシーケンス
【０１０７】
【数８】
2013519452.tif 000010
【０１０８】
も、Ｑの解に収束する。ただし、全てのｋ≧０について、
【０１０９】
【数９】
2013519452.tif 000011
【０１１０】
であるときに限る。ここで、βｋｖｋは、制限がある摂動であり、βｋが、
【０１１１】
【数１０】
2013519452.tif 000012
【０１１２】
であるような非負の実数であり、ベクトルのシーケンス
【０１１３】
【数１１】
2013519452.tif 000013
【０１１４】
が制限されることを意味する。
【０１１５】
幾つかの実施態様では、優秀化方法は、次のように利用することができる。制約型最小化問題を解こうとする代わりに、優秀化手法は、幾つかの実行可能探求アルゴリズムに摂動を与え得るため、実行可能性に対するその収束性を失うことなく、減少した目的関数値（複数の場合もあり）を有する１つの点（又は複数の点）が得られることになる。そのため、幾つかの実施態様では、ｐＣＴ撮像のための１つ又は複数の実行可能探求投影アルゴリズムは、こうした摂動弾力性特性を含むことができるか、又は含むように適合することができる。こうした能力を含むように適合することができるこうした１つ又は複数の摂動弾力性投影アルゴリズムの非制限的な例は、本明細書でより詳細に述べられる。
【０１１６】
出願人は、優秀化方法のために利用することができる複数の投影アルゴリズムを解析し、こうした解析の幾つかの結果が図２１〜図２３を参照して述べられる。図２１Ａ及び図２１Ｂは、例示的な投影アルゴリズムＡＲＴ（６５０）、ＢＩＰ（６５２）、ＳＡＰ（６５４）、ＯＳＡＲＴ（６５６）、ＢＩＣＡＶ（６５８）、ＤＲＯＰ（６６０）、及びＣＡＲＰ（６６２）を使用した、シミュレートされたＧＥＡＮＴ４陽子相互作用の再構成についての性能曲線（例えば、反復サイクル数の関数としての相対誤差）を示す。図２２は、陽子相互作用のシミュレーションで使用されたファントムの画像（上の行の第１列）、及び上記の投影アルゴリズムを使用した再構成画像を示す。図２３は、例示的な投影アルゴリズム（ＤＲＯＰ）の異なる数の反復によって得られた画像の比較を示し、幾つかの状況では、或る点を超える更なる反復は、必ずしも画像品質を大幅に上げないことを示す。
【０１１７】
図２１〜図２３の例示的な性能関連プロット及び性能関連画像は、既知の最適化技法を使用して得られた。これらの例示的な性能関連図に対する更なる参照が、例示的な投影アルゴリズムの説明においてなされる。
【０１１８】
幾つかの実施態様では、摂動弾力性投影アルゴリズムは、当業者に知られている、ＡＲＴ（代数再構成技法）等の順次投影アルゴリズムとすることができる。陽子ＣＴ再構成に適用されるＡＲＴアルゴリズムは、その順次的性質により遅い可能性がある。ＡＲＴアルゴリズムに関する更なる詳細は、例えば、G.T.Herman「Fundamentals of Computerized Tomography: Image Reconstruction from Projections」（2nded., Springer, New York, NY, 2009）において見出すことができる。図２１Ａ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるＡＲＴ投影アルゴリズムは、曲線６５０及び上の行の第２列の画像として示される。
【０１１９】
幾つかの実施態様では、ＡＲＴアルゴリズムは、同時ＡＲＴ（ＳＡＲＴ）アルゴリズムとして修正することができる。幾つかの実施態様では、ＳＡＲＴアルゴリズム等のＡＲＴベースのアルゴリズムは、ブロック反復アルゴリズム（例えば、ブロック反復ＳＡＲＴ又はＯＳ−ＳＡＲＴ）として修正することができる。陽子ＣＴ再構成に適用されるＳＡＲＴアルゴリズムは、最小２乗の最小に収束するように構成することができる。しかし、アルゴリズムは、重み付け値（１／ｍ）が比較的小さいため、依然として遅い可能性がある。ここで、量ｍは、データ集合内の陽子の数を表す。図２１Ｂ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるようなＯＳ−ＡＲＴ投影アルゴリズムが、曲線６５６及び下の行の第１列の画像として示される。
【０１２０】
幾つかの実施態様では、摂動弾力性投影アルゴリズムは、ブロック反復投影（block-iterative projection）（ＢＩＰ）アルゴリズムとすることができる。こうしたアルゴリズムは、Aharoni及びCensorによって開発され、このアルゴリズムに関する更なる詳細は、例えば、R.Aharoni及びY.Censor「Block-iterative projection methods for parallel computation of solutions to convex feasibility problems」（Linear Algebra and its Applications, 120, 165-175（1989））に見出すことができる。ＢＩＰアルゴリズムの幾つかの実施態様では、同時投影は、超平面の各ブロック内で行うことができ、超平面の複数のこうしたブロックは、順次に反復的に処理することができる。幾つかの実施態様では、ＢＩＰアルゴリズムは、ブロックサイズ（複数の場合もあり）に応じて重み付けを提供するように構成することができ、それにより、１／ｎ重み付け値に関連する再構成が遅いという問題が回避される。図２１Ａ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるようなＢＩＰ投影アルゴリズムが、曲線６５２及び上の行の第３列の画像として示される。
【０１２１】
幾つかの実施態様では、摂動弾力性投影アルゴリズムは、ストリング平均化投影（string-averaging projection）（ＳＡＰ）アルゴリズムとすることができる。こうしたアルゴリズムは、Censor、Elfving及びHermanによって開発され、このアルゴリズムに関する更なる詳細は、例えばY.Censor、T.Elfving及びG.T.Herman「Averaging strings of sequential iterations for convex feasibility problems」（Inherently Parallel Algorithms in Feasibility and Optimization and Their Applications, Elsevier Science Publications, Amsterdam, The Netherlands, D. Butnariu, Y. Censor, and S. Reich （Ed）,101-114 （2001））に見出すことができる。ＳＡＰアルゴリズムの幾つかの実施態様では、それぞれが超平面のストリングを有する複数のストリングを形成することができる。順次投影は、各ストリング内で実施することができ、こうした順次投影は、複数のストリングについて並列に実施することができる。こうした投影は、ストリングの全ての凸の組合せの反復として処理することができる。図２１Ａ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるようなＳＡＰ投影アルゴリズムが、曲線６５４及び上の行の第４列の画像として示される。
【０１２２】
幾つかの実施態様では、摂動弾力性投影アルゴリズムは、コンポーネント平均化（ＣＡＶ）アルゴリズムとすることができる。こうしたアルゴリズムは、Censor、Gordon及びGordonによって開発され、このアルゴリズムに関する更なる詳細は、例えばY.Censor、D.Gordon及びR.Gordon「Component averaging: An efficient iterative parallel algorithm for large and sparse unstructured problems」（Parallel Computing, 27, 777-808 （2001））に見出すことができる。幾つかの実施態様では、ＣＡＶアルゴリズムは、実質的に完全に同時であるように構成することができ、１／ｎ重み付け係数は、対角要素がｊ番目のボクセルを交差する陽子の数に等しい対角行列の族によって置換することができる。こうした構成によって、非直交投影を得ることができる。幾つかの実施態様では、ＣＡＶアルゴリズムは、ブロック反復ＣＡＶ（ＢＩＣＡＶ）アルゴリズムとして動作するように修正することができる。幾つかの実施態様では、ＣＡＶアルゴリズムは、コンポーネント平均化行投影（component averaging row projection）（ＣＡＲＰ）アルゴリズムとして動作するように修正することができる。図２１Ｂ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるようなＢＩＣＡＶ投影アルゴリズムが、曲線６５８及び下の行の第２列の画像として示される。図２１Ｂ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるようなＣＡＲＰ投影アルゴリズムが、曲線６６２及び下の行の第４列の画像として示される。
【０１２３】
幾つかの実施態様では、摂動弾力性投影アルゴリズムは、対角緩和した直交緩和（ＤＲＯＰ）アルゴリズムとすることができる。こうしたアルゴリズムは、Censor、Herman、Elfving、及びNikazadによって開発され、このアルゴリズムに関する更なる詳細は、例えばY.Censor、T.Elfving、G.T.Herman及びT.Nikazad「On diagonally-relaxed orthogonal projection methods」（SIAM Journal of Scientific Computing, 30, 473-504（2008））に見出すことができる。幾つかの実施態様では、ＤＲＯＰアルゴリズムは、実質的に完全に同時であるように構成することができ、重み付け手法は、ＣＡＶアルゴリズムの重み付け手法と同様とすることができる。しかし、結果として得られる投影は直交とすることができる。幾つかの実施態様では、ＤＲＯＰアルゴリズムは、ブロック反復ＤＲＯＰ（ＢＩＤＲＯＰ）アルゴリズムとして動作するように修正することができる。図２１Ｂ及び図２２では、例示的なＧＥＡＮＴ４シミュレーションに適用されるようなＤＲＯＰ投影アルゴリズムが、曲線６６０及び下の行の第３列の画像として示される。
【０１２４】
幾つかの実施態様では、摂動弾力性投影アルゴリズムは、修正重み付けスキームを有する同時代数再構成技法（simultaneous algebraic reconstruction technique）（ＳＡＲＴ）アルゴリズムとすることができる。こうしたアルゴリズムは、Anderson及びKakによって開発され、このアルゴリズムに関する更なる詳細は、例えばA.H.Andersen及びA.C.Kak「Simultaneous algebraic reconstruction technique（SART）: A superior implementation of the ART algorithm」（Ultrasonic Imaging, 6, 81-94 （1984））に見出すことができる。幾つかの実施態様では、こうしたアルゴリズムは、実質的に完全に同時であるように構成することができ、重み付けは、各行列の行の非ゼロ要素の和に対して１の量とすることができる。幾つかの実施態様では、こうしたアルゴリズムは、ブロック集合反復（ＯＳＡＲＴ）アルゴリズムとして動作するように修正することができる。
【０１２５】
他の投影アルゴリズムも、本明細書で述べる例示的な総変動優秀化（total variation superiorization）方法がある状態で又はない状態で利用することができる。こうした投影アルゴリズムの例に関する更なる詳細は、（ｉ）D.Butnariu、Y.Censor及びS.Reich （Editors）「Inherently Parallel Algorithms in Feasibility and Opimization and Their Applications」（Elsevier Science Publishers, Amsterdam, The Netherlands, 2001）、（ｉｉ）Y.Censor及びS.A.Zenios「Parallel Optimization: Theory, Algorithms, and Applications」（Oxford University Press, New York, NY, USA, 1997）、及び（ｉｉｉ）Y.Censor、W.Chen、P.L.Combettes、R.Davidi及びG.T.Herman「On the effectiveness of projection methods for convex feasibility problems with linear inequality constraints」（Technical Report, December 22, 2009（http://arxiv.org/abs/0912.4367のarXivで入手可能））等の文献に見出すことができる。
【０１２６】
本明細書で述べるように、優秀化手法は、或るメリット関数に関して優れた解を見出すプロセスを含むことができ、こうした優れた解は、摂動弾力性投影から得られる集合を含む、対応する凸実行可能性問題（convex feasibility problem）（ＣＦＰ）集合の実行可能解でもある。幾つかの実施態様では、こうしたメリット関数は、優秀化を総変動優秀化（ＴＶＳ）とすることができるような総変動（ＴＶ）関数を含むことができる。
【０１２７】
総変動優秀化技法を実証するために、出願人は、本明細書で述べるように、ＤＲＯＰ投影アルゴリズム（ＢＩＤＲＯＰ）のブロック反復バージョンを使用した。こうしたＴＶＳ技法はまた、本明細書で開示したアルゴリズムを含む他の投影アルゴリズムとともに実装することができることが理解されるであろう。
【０１２８】
例示的なＢＩＤＲＯＰアルゴリズムは、超平面指数の集合Ｉを、Ｉ＝Ｉ１∪Ｉ２∪．．．∪ＩＭであるようなＭ個のブロックの和集合として分割することによって使用することができる。幾つかの実施態様では、Ｍ個のブロックは、固定ブロックとすることができる。本明細書で述べるように、ＢＩＤＲＯＰアルゴリズムは、コンポーネント依存の重み付けスキームを使用するブロック反復投影法の変形であるとみなすことができる。
【０１２９】
幾つかの実施態様では、ＢＩＤＲＯＰは、次のように構成することができる。初期化は、オブジェクトベクトルｘの第１の要素ｘ０の値を選択することを含むことができる。ｘ０のこうした値は、任意とすることもできるし、例えば、フィルタ補正逆投影（ＦＢＰ）法を使用して推定することもできる。
【０１３０】
幾つかの実施態様では、反復ステップは、所与のベクトル要素ｘｋについて、
【０１３１】
【数１２】
2013519452.tif 000014
【０１３２】
である、次の反復ｘｋ＋１の計算を含むことができる。方程式７では、対角行列Ｕｔ（ｋ）は、
【０１３３】
【数１３】
2013519452.tif 000015
【０１３４】
として表すことができる。ここで、ｈｔｊは、ｊ番目のボクセルに交差するｔ番目のブロック内の陽子履歴の数であり、
【０１３５】
【数１４】
2013519452.tif 000016
【０１３６】
は、ユーザ選択の緩和パラメータのシーケンスである。本明細書で述べるＴＶＳ技法の例示的な実証では、λは、ＢＩＤＲＯＰに関する結果に基づいて約１．９の値に維持され、データは、図２３に示すように、１２ブロックに分割された。さらに、ブロックは、ｔ（ｋ）＝ｋ ｍｏｄ Ｍ＋１であるような周期制御である制御シーケンス
【０１３７】
【数１５】
2013519452.tif 000017
【０１３８】
によるアルゴリズムによって取上げられる。例示的な陽子データ集合は、実質的に等しいサイズで、かつ、各投影角度から実質的に等しい数の陽子履歴を有する１２個のブロックに分割された。
【０１３９】
幾つかの実施態様では、メリット関数φ及び／又は近接性関数Ｐｒは、優秀化再構成を操向するように選択することができる。実行可能性近接性関数は、測定された積分相対阻止能（ＲＳＰ）値の残りと、現在の画像推定によって得られた値とを提供するように選択することができる。こうした実行可能性近接性チェックは、メリット関数φによって表される更なるタスクに関する優秀化が、測定データと一致しない解を操向しないことを確実にするために利用することができる。幾つかの実施態様では、測定データに対する現在の画像推定ｘｋの実行可能性近接性は、
【０１４０】
【数１６】
2013519452.tif 000018
【０１４１】
として計算することができる。ここで、ｍは、対象となる集合内の陽子履歴の数である。
【０１４２】
幾つかの状況では、（制約の交差の）実行可能集合は空とすることができる。こうした状況でも、方程式８の近接性関数を低減することは、制約に低い程度に「違反する（violate）」ことがより少ない点をもたらし、したがって、たとえ近接性関数がゼロの値に達しない（また、或る状況では、達することができない可能性がある）場合でも、有用であり得る。
【０１４３】
幾つかの実施態様では、再構成プロセス中に低減され、かつ、再構成された画像推定の総変動に関連するメリット関数φは、
【０１４４】
【数１７】
2013519452.tif 000019
【０１４５】
として表すことができる。ここで、ｐｋは、ｎ次元画像ベクトルｘｋの２次元Ｊ×Ｊ表現である。
【０１４６】
幾つかの実施態様では、画像推定の低減した総変動に向かって画像推定の反復シーケンスを操向するための摂動ベクトルｖを計算することができる。例えば、摂動ベクトルは、
【０１４７】
【数１８】
2013519452.tif 000020
【０１４８】
になるように、ｘｋにおける総変動の正規化された劣勾配（subgradient）を負数にしたものとして計算することができる。幾つかの実施態様では、総変動の例示的な劣勾配（複数の場合もあり）は、P.L.Combettes及びJ.Luoによって著された文献、「An adaptive level set method for nondifferentiable constrained image recovery」（IEEE Trans. Image Process, 11, 1295-1304 （2002））に述べられている方法で計算することができる。上記の例示的な摂動ベクトルに関する更なる詳細は、例えば、D. Butnariu、R. Davidi、G. T. Herman及びI. G. Kazantsevによって著された文献、「Stable convergence behavior under summable perturbations of a class of projection methods for convex feasibility and optimization problems」（IEEE J. Sel. Top. Signal Process, 1, 540-547, （2007））に見出すことができる。
【０１４９】
例１
上記の例示的なＤＲＯＰベースの優秀化技法は、複数の方法で実装することができる。第１の例では、続く反復プロシージャのための初期画像推定は、測定データからフィルタ補正逆投影（ＦＢＰ）再構成を実施することによって取得することができる。こうしたＦＢＰは、既知の超音波画像グリッドに実質的に整合するように、個々の陽子履歴をリビニングすることによって実行することができる。反復プロシージャの場合、ＤＲＯＰ投影作用素は、データ集合の全てのブロックが処理されるまで、周期的に適用することができる。幾つかの実施態様では、データ集合内の実質的に全ての履歴を含む、本明細書で述べるオプションの実行可能性近接性計算をチェックすることができる。
【０１５０】
幾つかの実施態様では、上記の第１の例示的なＤＲＯＰベースの総変動優秀化技法は、複数の方法でコード化することができる。例えば、こうしたコード化は、以下のようなアルゴリズムを含むことができる。
（１）ｋ＝０にセットする
（２）ｘｋ＝ｘＦＢＰ初期ＦＢＰ再構成、及びβｋ＝１にセットする
（３）Ｎサイクルにわたって以下を繰り返す
（４）  ｘｋにおいてｓをφの劣勾配にセットする
（５）  ‖ｓ‖＞０である場合、ｖｋ＝−ｓ／‖ｓ‖にセットする
（６）  ‖ｓ‖＞０でない場合、ｖｋ＝ｓにセットする
（７）  ｃｏｎｔｉｎｕｅ＝ｔｒｕｅにセットする
（８）  ｃｏｎｔｉｎｕｅである間、以下を行う
（９）    ｙｋ＝ｘｋ＋βｋｖｋにセットする
（１０）    方程式９によってメリット関数（総変動）を計算し、φ（ｙｋ）≦φ（ｘｋ）である場合、
（１１）      ｙｋ（方程式７）に投影作用素Ｐｔ（ｋ）を順次にＭ回適用する
（１２）      Ｍ個全てのブロックからの履歴を使用して、方程式８によって実行可能性近接性を計算し、Ｐｒ（ＰＭｙ）＜Ｐｒ（ｘｋ）である場合、
（１３）      ｘｋ＋１＝ＰＭｙにセットする
（１４）      ｃｏｎｔｉｎｕｅ＝ｆａｌｓｅにセットする
（１５）    Ｐｒ（ＰＭｙ）＜Ｐｒ（ｘｋ）でない場合、βｋ＝βｋ／２にセットする
（１６）  φ（ｙｋ）≦φ（ｘｋ）でない場合、βｋ＝βｋ／２にセットする
（１７）ｋ＝ｋ＋１にセットする
幾つかの実施態様では、サイクルの数を表す量Ｎを、例えば１０にセットすることができる。Ｎの他の値も使用することができる。
【０１５１】
例２
第２の例では、ＤＲＯＰ投影作用素を、後続のブロックからの履歴によって実施されるオプションの実行可能性近接性チェックを継続する前に、所与のブロックに適用することができる。こうした手法は、各ブロックが、各投影角度から実質的に等しい数の履歴を含み、したがって、全体としてデータを表し得るときに利用することができる。第１の例と同様に、初期画像推定は、測定データからフィルタ補正逆投影（ＦＢＰ）再構成を実施することによって取得することができる。
【０１５２】
幾つかの実施態様では、上記の第２の例示的なＤＲＯＰベースの総変動優秀化技法は、複数の方法でコード化することができる。例えば、こうしたコード化は、以下のようなアルゴリズムを含むことができる。
（１）ｋ＝０にセットする
（２）ｘｋ＝ｘＦＢＰ初期ＦＢＰ再構成、及びβｋ＝１にセットする
（３）Ｎサイクルにわたって以下を各ブロックについて繰り返す
（４）  ｘｋにおいてｓをφの劣勾配にセットする
（５）  ‖ｓ‖＞０である場合、ｖｋ＝−ｓ／‖ｓ‖にセットする
（６）  ‖ｓ‖＞０でない場合、ｖｋ＝ｓにセットする
（７）  ｃｏｎｔｉｎｕｅ＝ｔｒｕｅにセットする
（８）  ｃｏｎｔｉｎｕｅである間、以下を行う
（９）    ｙｋ＝ｘｋ＋βｋｖｋにセットする
（１０）    方程式９によってメリット関数（総変動）を計算し、φ（ｙｋ）≦φ（ｘｋ）である場合、
（１１）      ｙ（方程式７）に投影作用素Ｐｔ（ｋ）を適用する
（１２）      後続のブロックからの履歴を使用して、方程式８によって実行可能性近接性を計算し、Ｐｒ（Ｐt（k）ｙ）＜Ｐｒ（ｘｋ）である場合、
（１３）      ｘｋ＋１＝Ｐt（k）ｙにセットする
（１４）      ｃｏｎｔｉｎｕｅ＝ｆａｌｓｅにセットする
（１５）    Ｐｒ（Ｐt（k）ｙ）＜Ｐｒ（ｘｋ）でない場合、βｋ＝βｋ／２にセットする
（１６）  φ（ｙｋ）≦φ（ｘｋ）でない場合、βｋ＝βｋ／２にセットする
（１７）ｋ＝ｋ＋１にセットする
幾つかの実施態様では、サイクルの数を表す量Ｎを、例えば１０にセットすることができる。Ｎの他の値も使用することができる。
【０１５３】
例３
第３の例では、第１の例の方法は、ステップ１２の実行可能性近接性チェックが実施されないように修正することができる。こうした修正は、所与の実行可能性近接性チェックが計算的に煩雑である場合に望ましい場合がある。
【０１５４】
例４
第４の例では、第２の例の方法は、ステップ１２の実行可能性近接性チェックが実施されないように修正することができる。こうした修正は、所与の実行可能性近接性チェックが計算的に煩雑である場合に望ましい場合がある。
【０１５５】
例５
第５の例では、限定はしないが、方程式９を含む総変動（ＴＶ）関数を利用することができる。幾つかの実施態様では、ＤＲＯＰ投影技法に基づくことができる第５の例は、複数の方法でコード化することができる。例えば、こうしたコード化は、ｘ０＝選択された値、かつ、β０＝１である場合に、初期化を行うことができることを含むことができる。ｘｋ及びβｋが与えられる反復ステップの場合、ｘｋ＋１及びβｋ＋１は、以下のように生成することができる。
（１）ブロック指数ｔを通してループする
（２）摂動ベクトルｖｋを計算する
（３）ｚｋ＝ｘｋ＋βｋｖｋにセットする
（４）ＴＶ（ｚｋ）＞ＴＶ（ｘｋ）である場合、βｋ＝βｋ／２にセットして（２）へ進み、ＴＶ（ｚｋ）＞ＴＶ（ｘｋ）でない場合、（５）へ進む
（５）方程式７と同様のＤＲＯＰアルゴリズムによって、現在のブロック内で、超平面上へのｚｋの投影を計算する
（６）Ｒｅｓ（ｙｋ）＞Ｒｅｓ（ｘｋ）である場合、βｋ＝βｋ／２にセットして（２）へ進み、Ｒｅｓ（ｙｋ）＞Ｒｅｓ（ｘｋ）でない場合、（７）へ進む
（７）ｘｋ＋１＝ｙｋにセットし、停止条件が適用されるまで、（１）〜（６）を繰返す
幾つかの撮像状況では、計算処理リソース（例えば、計算パワー及び計算時間）を犠牲にしてオブジェクトの高品質画像を得ることが可能である。意図する用途にとって十分な品質を有し、リソースの許容可能な量を使用する画像を取得することが一般に望ましい。本開示の１つ又は複数の特徴が、陽子ＣＴアプリケーションに適用されると、こうした望ましい品質及び／又は性能の利点を提供することができると考えられている。
【０１５６】
図２４は、幾つかの実施態様において、プロセス７００を、所望の品質及び／又は再構成性能特徴を有する陽子ＣＴ画像の形成を可能にする解を見出すために実施することができることを示す。ブロック７０２において、オブジェクトを通過した複数の陽子が検出することができる。こうした陽子検出の例（例えば、追跡及びエネルギー損失測定）が本明細書で述べられる。ブロック７０４において、オブジェクトを通って各陽子がとる経路を、オブジェクト内での陽子の相互作用を特徴付けるために推定することができる。こうした経路推定の例が本明細書で述べられる。ブロック７０６において、陽子の相互作用履歴を、実行可能解の凸集合に変換することができる。こうした変換の例が本明細書で述べられる。ブロック７０８において、凸集合内にある解を、変換に関連するメリット関数に基づいて見出すことができ、その解は、少なくとも１つの選択された推定又は既存の解より優れている。こうした優秀化方法の例が本明細書で述べられる。
【０１５７】
図２５は、幾つかの実施態様において、例えば図２４のプロセス７００によって得られた優れた解に基づいて画像を得るためにプロセス７１０を実施することができることを示す。ブロック７１２において、凸実行可能性集合からの優れた解を得ることができ、ここで、凸実行可能性集合は、複数の陽子のそれぞれと、撮像されるオブジェクトとの相互作用を特徴付ける。ブロック７１４において、オブジェクトは、ブロック７１２で得られた優れた解に基づいて特徴付けすることができる。例えば、所望のボクセル値を有するベクトルは、優れた解から計算することができる。ブロック７１６において、オブジェクトの画像は、オブジェクトの特徴付けに基づいて形成することができる。
【０１５８】
図２６Ａ及び図２６Ｂは、ＧＥＡＮＴ４によってシミュレートされた例示的な陽子相互作用構成８００を示し、図２７〜図２９は、本明細書で述べる再構成方法から得られる例示的な結果を示す。
【０１５９】
例示的なシミュレーション構成８００では、検出器システムは、上流追跡検出器８０４及び下流追跡検出器８０６を含む。上流検出器８０４及び下流検出器８０６のそれぞれは、個々の陽子の入射及び出射の場所及び方向の確定を可能にするために２つの追跡平面を含む。各追跡平面は、約３０ｃｍ×３０ｃｍの寸法と、約０．０４ｃｍの厚さと、約１００μｍの割り当てられた空間分解能とを有する２Ｄ感応性シリコン追跡平面として構成される。
【０１６０】
シミュレーション構成８００の検出器システムは、約３２ｃｍ×３２ｃｍの横方向寸法と約１０ｃｍの厚さとを有するＣｓＩ結晶熱量計８０８を更に含む。熱量計８０８の応答をシミュレートするために、検出器の雑音源が無視され、検出器は、完全なエネルギー分解能を有するものとして構成された。
【０１６１】
例示的なシミュレーション構成８００では、２つの異なる仮想ファントム８１０が、空間分解能及び相対阻止能分解能を定量化するために使用された。その断面図が図２６Ｂに示される両方のファントムは、約１６ｃｍの直径を有する。各ファントムは、国際放射線防護委員会（ＩＣＲＰ）によって規定されたような、化学組成及び電子密度が脳８１２及び頭骸骨８１４に実質的に同等な２つの物質を含む。
【０１６２】
図２６Ｂでは、第１の例示的なファントム８１０ａは、選択された再構成ボクセル寸法にほぼ等しい、約（０．８２×０．８２）ｍｍ２の断面を有する、中心長方形プリズム構造８１６を含む。こうしたファントムは、所与の再構成方法の空間分解能の推定を可能にすることができる。第２の例示的なファントム８１０ｂは、こうした中心構造を含まない。例示的なシミュレーションのために、中心構造８１６の電子密度は、周囲の脳物質の電子密度より約２０倍大きいが、同じ化学組成を保持するように選択された。
【０１６３】
図２６Ａに示すように、約２００ＭｅＶのエネルギーを有する実質的に単一エネルギーの陽子は、２Ｄ平行ビーム幾何形状を形成した。約２度間隔を有する１８０回の投影が、各ファントムについてシミュレートされた。各投影角度について、追跡平面上の位置及び熱量計に付与されたエネルギーが、２００００個の陽子について記録された。
【０１６４】
陽子とファントムとの相互作用をシミュレートするために、ＧＥＡＮＴ４ハドロンイオン化モデルが使用された。モデルは、ファントムを横切る陽子について対象となるエネルギー範囲をカバーする、約２ＭｅＶを超える陽子エネルギーについてベーテブロッホ（Bethe-Block）関係を使用した。イオン化エネルギー損失について、ＧＥＡＮＴ４構成は、１ｋｅＶ〜１００ＴｅＶの運動エネルギーにおいて対数的に等間隔で配置された１００ステップにおける平均値の計算を含む。しかし、こうしたデフォルト構成は、陽子ＣＴアプリケーションにとって十分に正確でない場合があることを調査が示した。したがって、エネルギービニングが、２０００ステップにおいて１ｋｅＶ〜５００ＭｅＶまで計算された。ＧＥＡＮＴ４における低エネルギー弾性核衝突モデル及び低エネルギー非弾性核衝突モデルが使用可能にされた。
【０１６５】
シミュレートされた陽子データ（入口座標及び出口座標並びに熱量計に付与されたエネルギーを含む）に基づいて、各ファントムの２Ｄ画像が、種々のＤＲＯＰベースの投影アルゴリズムによって、また、優秀化方法がある状態及びない状態で再構成された。所与のブロック内の投影等の計算（例えば、他のブロックと独立であり得る方程式７の加算）の並列実行を実施するために、ＧＰＧＰＵ（汎用グラフィカル処理ユニット）が使用された。所与のＧＰＧＰＵ上でのブロック投影の完了に続いて、更なる処理のために、加算されたアレイが、プロセッサ（例えば、ＣＰＵ）に戻された。ＣＰＵはまた、ブロック反復アルゴリズムの順次部分を実施した。
【０１６６】
上記の方法で得られた画像は、以下のように更に解析された。再構成画像の空間分解能は、空間周波数の関数として撮像システムの信号伝達特性の尺度を提供することができる２Ｄ変調伝達関数（ＭＴＦ）によって定量化された。こうした尺度について、中心にある密な長方形プリズム（図２６Ｂの８１６）の画像の点拡がり関数（point spread function）（ＰＳＦ）が使用された。再構成に続いて、ＰＳＦを中心とする１６×１６の２Ｄボクセルの関心領域の２Ｄ高速フーリエ変換（ＦＴＴ）が実行された。ファントムの軸対称性を利用して、ＭＴＦは、画像の、結果として得られる空間周波数表現のｘ軸成分及びｙ軸成分の大きさを平均することによって、関心領域において得られた。
【０１６７】
低コントラスト密度分解能は、所与のサイズのオブジェクトを周囲組織から識別するのに必要とされる最小コントラストを確定するためのオブジェクト統計解析法を提供することができるコントラスト識別関数（ＣＤＦ）によって評価された。ＣＤＦは、一様なファントム（図２６Ｂの８１０ｂ）の再構成画像を、サイズが１×１〜１０×１０の２Ｄボクセルの範囲のオブジェクトのグリッドに分割することによって計算された。グリッド要素内の平均ピクセル値の分布の標準偏差は、所与の信頼レベルで検出可能な最小コントラストを確定するために使用された。９５％の信頼レベルの場合、選択されたグリッドサイズのオブジェクトと背景との間の検出可能な密度差は、平均ピクセル値分布の３．２９標準偏差として規定された。
【０１６８】
再構成された相対阻止能（ＲＳＰ）値の量的精度は、ヒストグラム解析を使用し、相対的ＲＳＰ誤差を
【０１６９】
【数１９】
2013519452.tif 000021
【０１７０】
として規定して、確定された。ここで、ｘ’ｊは、ファントムのボクセルｊ内のＲＳＰであり、ｘｎｊは、ｎサイクル後のボクセルｊ内の再構成されたＲＳＰである。
【０１７１】
図２７は、（ａ）総変動優秀化なしのＤＲＯＰと、（ｂ）ＴＶＳを用いるが、（第３の例としても述べられる）第１の例における実行可能性近接性チェックがない状態のＤＲＯＰ（ＴＶＳ１−ＤＲＯＰ＊）と、（ｃ）ＴＶＳを用いるが、（第４の例としても述べられる）第２の例における実行可能性近接性チェックがない状態のＤＲＯＰ（ＴＶＳ２−ＤＲＯＰ＊）とによって再構成された画像を示す。上の行の画像は、一様なファントム（図２６Ｂの８１０ｂ）の再構成であり、下の行の画像は、空間分解能ファントム（図２６Ａの８１０ａ）の再構成である。例示的な画像は、０．８と１．２との間にＲＳＰ値を含むように切断したＲＳＰ窓から得られる。
【０１７２】
図２８Ａ〜図２８Ｃに示す種々の例示的な結果において、実行可能性近接性チェックを用いるＤＲＯＰベースのＴＶＳ法は、「＊」シンボルなしで示される。
【０１７３】
上記の方法を用い、かつ、実行可能性近接性チェックなしで再構成された画像は、実行可能性近接性チェックを用いて再構成された画像と比較すると、より小さいか又はほぼ同等な最小相対誤差を有していた。こうした作用が図２８Ａに示される。ＤＲＯＰベースのＴＶＳ法は全て、実行可能性近接性チェックを用いても用いなくても、少なくとも１０まで延びる反復サイクル数範囲について、ＴＶＳなしのＤＲＯＰ（DROP-without-TVS）法より小さな相対誤差を有する。図２７に示す画像は、ＴＶＳなしのＤＲＯＰの場合、サイクル３であり、ＴＶＳ１−ＤＲＯＰ＊及びＴＶＳ２−ＤＲＯＰ＊の場合、サイクル１０である、最小相対誤差の各サイクルにおいて得られた画像に対応することが留意されるべきである。定性的には、ＴＶＳ２−ＤＲＯＰ＊スキームが、おそらく余分の摂動ステップのために、最も低い雑音レベルを有したことがわかる。
【０１７４】
図２８Ａは、種々のスキームについてのサイクル数の関数としての相対誤差のプロットを示す。プロットにおいて、サイクル０の相対誤差は、反復ＴＶ優秀化の初期点を生成するのに使用されたＦＢＰアルゴリズムによって生成される相対誤差に対応する。
【０１７５】
ＴＶＳ１−ＤＲＯＰ及びＴＶＳ１−ＤＲＯＰ＊スキーム（実行可能性近接性チェックを用いた状態及び用いない状態）によって再構成された画像は、一般に、定量的な相対阻止能（ＲＳＰ）精度の点で同等であり、相対誤差は、サイクル数が増加するにつれて、単調減少傾向に従う。実行可能性近接性チェックの排除は、この場合チェック条件が違反されないため、大きな差を生じない。
【０１７６】
一方、実行可能性近接性チェックの排除は、ＴＶＳ２−ＤＲＯＰスキームの場合、大きな差を生じる。図２８Ａは、実行可能性近接性チェックを含むことが、ＤＲＯＰアルゴリズムと同様の、最小に達した後の相対誤差の漸進的な増加をもたらし得ることを実証する。こうした作用は、減少したβが、摂動ステップの雑音低減効果を弱め得ることによって説明することができる。したがって、ＤＲＯＰアルゴリズムが低い相対誤差から離れるため、より厳格なＴＶＳ手法も低い相対誤差から離れる。これは、ＴＶＳ２−ＤＲＯＰに関しては起こるが、ＴＶＳ１−ＤＲＯＰに関しては起こらない。その理由は、実行可能性近接性条件の違反が、前者に関しては観測されるが、後者に関しては観測されないからである。実行可能性近接性チェックがない状態では、ＴＶＳ２−ＤＲＯＰ＊スキームの相対誤差は、１０サイクル内で単調減少傾向に従う。最初の１０サイクル内での最小相対誤差は、ＤＲＯＰに関して約２．６４％、ＴＶＳ１−ＤＲＯＰ及びＴＶＳ１−ＤＲＯＰ＊に関して約１．９６％、ＴＶＳ２−ＤＲＯＰに関して約１．６４％、及びＴＶＳ２−ＤＲＯＰ＊に関して約１．５５％である。これらの差は、異なるスキームによって再構成された画像内の様々な程度の雑音から生じ得る。
【０１７７】
上記の例示的な結果は、１２ブロックに細分されたデータによって得られた。結果の別の集合は、１８０ブロックに細分されたデータによって得られ、その結果（図示せず）は、１２ブロックの結果に非常に似ていた。したがって、幾つかの実施態様では、総変動優秀化方法によって処理されたブロック反復再構成アルゴリズムは、ブロックの数の或る変化により、実質的に影響されないことができる。
【０１７８】
図２７の上の行に提示されている画像のヒストグラムは、脳及び骨等価領域の平均再構成値を解析するために生成された。再構成雑音をモデル化するために、ガウス分布がピークに当てはめられた。全てのスキームは、ピーク当てはめ不確実性内で、脳及び骨等価領域について実質的に同じ平均ＲＳＰ値を再構成した。したがって、ＴＶＳ摂動スキームは、これらの物質の再構成値の精度に悪い影響を及ぼさなかった。
【０１７９】
ＴＶＳ１−ＤＲＯＰ＊及びＴＶＳ２−ＤＲＯＰ＊スキーム（ともに、雑音性能を改善し、再構成時間を低減した）を使用した更なる解析によって、図２８Ｂに示すＭＴＦプロットを得た。プロット内の任意の空間周波数について、ＴＶＳ１−ＤＲＯＰ＊スキームはより大きなＭＴＦ値、したがって、ＤＲＯＰスキームの空間分解能に比べて、優れた空間分解能を有する。ＴＶＳ１−ＤＲＯＰ＊及びＴＶＳ２−ＤＲＯＰ＊スキームは、空間分解能の点で同様に機能するが、ＴＶＳ１−ＤＲＯＰ＊スキームが僅かに良好である。ＴＶＳ再構成スキームに関する空間分解能の改善は、最低の相対誤差に達する前に実施されるサイクル数がより多いことに起因し得る。反復アルゴリズムとともにＭＬＰ形式を使用すると、ｐＣＴ空間分解能がサイクル数とともに改善することが過去に観測されている。これは、密度分解能を改善する一部の再構成法（図２８Ｃを参照して本明細書で述べる）が、劣った空間分解能を表示することが多いため、重要な結果であり得る。
【０１８０】
ＤＲＯＰ及び低減型（reduced）ＴＶＳ−ＤＲＯＰ＊スキームに関連するＣＤＦは、図２８Ｃにプロットされる。ＴＶＳ１−ＤＲＯＰ＊スキームは、ＤＲＯＰスキームよりほんの僅かだけ良好に機能するが、ＴＶＳ２−ＤＲＯＰ＊スキームは、他の２つのスキームよりはるかに良好に機能する。約１ｍｍ２程度の小さいオブジェクトの場合、ＴＶＳ２−ＤＲＯＰ＊スキームは、約１％と約１．５％との間のコントラスト識別を可能にする。ＴＶＳ２−ＤＲＯＰ＊スキームの優れたコントラスト識別は、画像雑音の低減と空間分解能の改善の組合せに起因し得る。
【０１８１】
図２９は、幾つかの実施態様において、総変動優秀化方法等の優秀化スキームが、相対阻止能（ＲＳＰ）の測定の改善ももたらすことができることを示す。図２９の例示的なプロットは、ファントムに関連する、既知の急峻な阻止能ピーク（実線）並びにＤＲＯＰ法（「ｘ」マーク）及びＴＶＳ−ＤＲＯＰ法（固体充填マーク）（例えば、本明細書で述べる第５の例）を使用した、シミュレートされたデータの再構成から得られた阻止能ピークを示す。ＤＲＯＰの例についてのＲＳＰ分布は、図１９Ｂを参照して述べた例（例えば、有効平均コード長技法）の最良のものと同様の性能を有する。こうしたＤＲＯＰ再構成法にＴＶＳ技法を実装することによって、結果として得られるＲＳＰの分布は、実際のＲＳＰ分布をかなり良く複製することが示される。
【０１８２】
図３０は、実験陽子データに適用されたＴＶＳ−ＤＲＯＰ再構成スキームを使用して得られた例示的な陽子ＣＴ画像を示す。画像は、空気（左上）、アクリル（右上、約１．２の相対阻止能）及び骨等価プラスチック（右下、約１．７の相対阻止能）の挿入物を有するポリスチレン（約１．０３５の相対阻止能）で作られた球ファントムの２ｍｍスライス再構成に対応する。異なる阻止能の物質をはっきりと識別することができることがわかる。
【０１８３】
本明細書で述べるように、優秀化スキームは、手近の撮像タスクに適合するようにｐＣＴ再構成に適用することができる。幾つかの実施態様では、ｐＣＴ又は他の反復画像再構成方法に適用される優秀化方法は、実行可能性探求投影方法の反復ステップの間で、計算された画像推定に摂動を与えることによって容易にすることができる。摂動方法を適切に選択することによって、再構成画像のシーケンスに対する、著しく有利な変更を達成することができる。
【０１８４】
本明細書で述べるように、総変動の低減に基づく、ＴＶＳ１−ＤＲＯＰ及びＴＶＳ２−ＤＲＯＰ等の例示的な優秀化スキームは、ＤＲＯＰ法に比べて改善された画像品質をもたらすことができる。より具体的には、ＴＶＳ２−ＤＲＯＰにおいて利用される更なる摂動ステップは、画像雑音のより大きな低減及び優れた密度分解能をもたらすことができる。
【０１８５】
幾つかの実施態様では、総変動優秀化方法の１つ又は複数の特徴に関連する計算時間は、提供される利益の観点で考慮することができる。例えば、ＴＶメリット関数（例えば、方程式９）及び摂動ベクトルｖｋの計算は、画像の寸法が比較的大きいとき、画像再構成時間を増加させ得る。別の例では、再構成時間は、実行可能性近接性関数（例えば、方程式８）の計算によって著しく増加し得る。したがって、幾つかの実施態様では、ＴＶＳスキームは、実行可能性近接性関数の省略によって減少した時間（例えば、約半分の時間）で実施することができる。ＴＶＳ２−ＤＲＯＰスキーム等の幾つかの実施態様では、画像雑音の低減を更にもたらす省略を行うことができる。ＴＶＳ１−ＤＲＯＰスキーム等の幾つかの実施態様では、省略は、他の性能パラメータに対して実質的に全く有害な影響を及ぼさない可能性がある。
【０１８６】
幾つかの実施態様において、空間分解能の改善は、両方のＴＶＳ促進性再構成スキームによって達成することができることも留意される。本明細書で述べる例示的なＴＶＳ１−ＤＲＯＰ＊スキームは、ＴＶＳ２−ＤＲＯＰ＊スキームに比べて僅かに優れた空間分解能を示す一方、ＴＶＳ２−ＤＲＯＰ＊スキームは、そのより良好な雑音低減にもかかわらず、ＤＲＯＰ再構成によって再構成された画像に比べて優れた空間分解能を依然としてもたらす。再構成画像を「平滑化する（smoothing）」ことによって密度分解能を改善しようとする以前の試みは、概して、空間分解能の低下をもたらしたことが注目されてきた。これは、例示的なＴＶＳ１−ＤＲＯＰ及びＴＶＳ２−ＤＲＯＰのスキームでは当てはまらず、ここでは、空間分解能は維持又は改善された。
【０１８７】
本明細書で述べるように、本開示の１つ又は複数の特徴は、画像等の有用な結果を得るために、陽子ＣＴデータの効果的かつ比較的高速な再構成を容易にすることができる。こうした有利な特徴は、こうした有用な結果を得るために必要とされる放射線線量を低減するためにも利用することができる。例えば、また、図５及び図６を参照すると、ＣＴに利用される陽子は、平均して、ブラッグピークが、撮像されるオブジェクトの下流の場所で起こるように構成することができる。したがって、オブジェクトは、陽子の放射線線量のかなりの部分を受けることを免れる。
【０１８８】
オブジェクトが受けることを予想され得る（ブラッグピークエネルギー損失がない）線量は、以下のように概略で推定することができる。撮像されるオブジェクトがボクセルに分割され、各ボクセルが１ｍｍ３立方体であり、スライス画像を得るために陽子投影が（図２６の例示的なシミュレーションと同様の）１８０の向きで得られると仮定する。概略推定のために、オブジェクトが、２５０ＭｅＶの陽子についてのエネルギー損失が約１００ＭｅＶであるような水によって表すことができる通常の頭部サイズのオブジェクトであるとも仮定する。このとき、こうしたスライスについて陽子当たりの線量は、約０．５ｎＧｙであると推定することができる。
【０１８９】
上記のスライスを表すベクトルの寸法は、１ｍｍ３立方体の例示的なボクセルサイズを仮定して、１８０２＝３．６×１０４として推定することができる。各ボクセルは、少なくとも１つの陽子によって横切られると仮定する。このとき、こうしたスライスをカバーするのに必要とされることが予想され得る陽子の数は、ベクトル寸法の大きさ程度のものであると推定することができる。したがって、所与のスライスが受けることが予想され得る線量は、（３．６×１０４）陽子×０．５ｎＧｙ／陽子＝１．８×１０４ｎＧｙ又は約２０μＧｙ（又は、２５０ＭｅＶ範囲内の陽子は一般に低ＬＥＴ放射線であるため、約２０μＳｖ）として推定することができる。こうした放射線量は、比較的低く、状況によっては、極度に低いとさえ考えられ得る。
【０１９０】
本開示の１つ又は複数の特徴に関する更なる詳細は、S.N.Penfold、R.W.Schulte、Y.Censor、A.B.Rosenfeldによる文献「Total variation superiorization schemes in proton computed tomography image reconstruction」（Med. Phys. 37（11）, November 2010）に見出すことができる。この文献は、参照によりその全体が本明細書に援用される。
【０１９１】
本明細書で述べる種々の例示的なプロセス及びアルゴリズムは、電子ハードウェア、コンピュータソフトウェア、又は両方の組合せとして実装することができる。種々のプロセス及びアルゴリズムを実装するために使用されるハードウェア及びデータ処理装置は、汎用のシングルチッププロセッサ若しくはマルチチッププロセッサ、デジタル信号プロセッサ（ＤＳＰ）、汎用グラフィック処理ユニット（ＧＰＧＰＵ）、特定用途向け集積回路（ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）若しくは他のプログラマブルロジックデバイス、ディスクリートゲート若しくはトランジスタロジック、ディスクリートハードウェアコンポーネント、又は、本明細書で述べた機能を実施するように設計されたそれらの任意の組合せによって実装又は実施することができる。プロセッサは、マイクロプロセッサ、又は任意の従来のプロセッサ、コントローラ、マイクロコントローラ、又はステート・マシンとすることができる。プロセッサはまた、コンピューティングデバイスの組合せ、例えば、ＤＳＰとマイクロプロセッサとの組合せ、複数のマイクロプロセッサ、ＤＳＰコアと連携する１つ若しくは複数のマイクロプロセッサ、又は、任意の他のこうした構成として実装することができる。幾つかの実施態様では、特定のステップ及び方法は、所与の機能に固有である回路要素によって実施することができる。
【０１９２】
本明細書で述べる種々の機能は、ハードウェア、デジタル電子回路要素、コンピュータソフトウェア、ファームウェアにおいて、又はそれらの任意の組合せにおいて実装することができる。本明細書で述べる再構成スキームの実施態様は、１つ又は複数のコンピュータプログラム（例えば、データ処理装置が実行するため、又は、データ処理装置の動作を制御するために、コンピュータ記憶媒体上に符号化されるコンピュータプログラム命令の１つ又は複数のモジュール）として実装することができる。
【０１９３】
ソフトウェアにおいて実装される場合、機能は、コンピュータ可読媒体上の１つ又は複数の命令又はコードとして記憶又は伝送することができる。本明細書で開示される方法又はアルゴリズムのステップは、コンピュータ可読媒体上に存在することができるプロセッサ実行可能ソフトウェアモジュールに実装することができる。コンピュータ可読媒体は、コンピュータプログラムを１つの場所から別の場所へ転送できるようにすることができる任意の媒体を含む、コンピュータ記憶媒体及び通信媒体の両方を含むことができる。記憶媒体は、コンピュータがアクセスすることができる任意の入手可能な媒体とすることができる。制限としてではなく例として、こうしたコンピュータ可読媒体は、ＲＡＭ、ＲＯＭ、ＥＥＰＲＯＭ、ＣＤ−ＲＯＭ若しくは他の光学ディスク記憶装置、磁気ディスク記憶装置若しくは他の磁気記憶デバイス、又は、命令若しくはデータ構造の形態で所望のプログラムコードを記憶するために使用されることができ、コンピュータによってアクセスすることができる任意の他の媒体を含むことができる。同様に、任意の接続を、コンピュータ可読媒体と適切に呼ぶことができる。本明細書で使用されるディスク（disk）及びディスク（disc）は、コンパクトディスク（compact disc）（ＣＤ）、レーザディスク（laser disc）、光ディスク（optical disc）、デジタル多用途ディスク（digital versatile disc）（ＤＶＤ）、フロッピーディスク（floppy disk）、及びブルーレイディスク（blu-ray disc）を含み、ディスク（disk）は通常、データを磁気的に再生する一方、ディスク（disc）は、データをレーザによって光学的に再生する。上記の組合せも、コンピュータ可読媒体の範囲内に含まれるべきである。さらに、方法又はアルゴリズムの動作は、コンピュータプログラム製品に組込むことができる、機械可読媒体及びコンピュータ可読媒体上のコード及び命令の１つ又は任意の組合せ又は集合として存在することができる。
【０１９４】
本明細書で使用される用語「約（approximately）」、「約（about）」、及び「実質的に（substantially）」は、依然として所望の機能を実施するか又は所望の結果を達成する、述べた量に近い量を表す。例えば、用語「約（approximately）」、「約（about）」、及び「実質的に（substantially）」は、述べた量の１０％未満内、５％未満内、１％未満内、０．１％未満内、及び０．０１％未満内である量を指すことができる。本明細書で使用される用語「の少なくとも一部分（at least a portion of）」は、全体の或る量を表し、この或る量は、全体を含むことができる全量を含む。例えば、用語「の一部分（a portion of）」は、全体の０．０１％より多い、０．１％より多い、１％より多い、１０％より多い、２０％より多い、３０％より多い、４０％より多い、５０％より多い、６０％より多い、７０％より多い、８０％より多い、９０％より多い、９５％より多い、９９％より多い、及び１００％の量を指すことができる。
【０１９５】
文脈上別途明確に必要としない限り、明細書及び特許請求の範囲全体を通じて、語「備える、含む（comprise）」、「備える、含む（comprising）」、及び同様なものは、排他的意味又は網羅的意味と対照的に、包含的意味で、すなわち、「含むが、それに限定されない（including, but not limited to）」という意味で解釈されるものとする。本明細書で一般に使用される語「結合した（coupled）」は、直接接続されるか又は１つ若しくは複数の中間要素によって接続することができる２つ以上の要素を指す。さらに、語「本明細書で（herein）」、「上記（above）」、「以下（below）」、及び同様の意味の語は、本出願で使用されるとき、本出願を全体として指し、本出願のいかなる特定の部分も指さないものとする。文脈が許可する場合、単数又は複数を使用する上記詳細な説明内の語はまた、複数又は単数をそれぞれ含むことができる。語「又は、若しくは（or）」は、２つ以上のアイテムのリストを参照し（is in reference to）、その語は、語の以下の解釈、すなわち、リスト内のアイテムの任意のものと、リスト内のアイテムの全てと、リスト内のアイテムの任意の組合せとの全てを包含する。
【０１９６】
本発明の実施形態の上記で詳述された説明は、網羅的であることも、上記で開示された厳密な形態に本発明を限定することも意図しない。本発明の特定の実施形態及び例が例証のために上述されるが、当業者が認識するように、本発明の範囲内で種々の等価な修正が可能である。例えば、プロセス又はブロックは所与の順序で提示されるが、代替の実施形態は、異なる順序で、ステップを有するルーチンを実施することも、ブロックを有するシステムを使用することもでき、幾つかのプロセス又はブロックは、削除し、移動し、追加し、細分し、組み合わせ、かつ／又は修正することができる。これらのプロセス又はブロックのそれぞれは、種々の異なる方法で実装することができる。同様に、プロセス又はブロックは、時として直列に実施されるものとして示されるが、これらのプロセス又はブロックは、代わりに、並列に実施することも、異なる時点で実施することもできる。
【０１９７】
本明細書で提供される本発明の教示は、必ずしも上述したシステムではなく、他のシステムに適用することができる。上述した種々の実施形態の要素及び動作は、更なる実施形態を提供するために組み合わせることができる。
【０１９８】
本発明の或る特定の実施形態が述べられたが、これらの実施形態は、単に例として提示されており、本開示の範囲を制限することを意図されない。実際には、本明細書で述べる新規な方法及びシステムは、種々の他の形態で具現化することができ、さらに、本明細書で述べる方法及びシステムの形態における種々の省略、置換及び変更を、本開示の趣旨から逸脱することなく行うことができる。添付特許請求の範囲及びその均等物は、本開示の範囲及び趣旨内に入る形態又は修正形態を包含することを意図される。
【図１Ａ】
2013519452.tif 000023
【図１Ｂ】
2013519452.tif 000024
【図２Ａ】
2013519452.tif 000025
【図２Ｂ】
2013519452.tif 000026
【図３Ａ】
2013519452.tif 000027
【図３Ｂ】
2013519452.tif 000028
【図４】
2013519452.tif 000029
【図５】
2013519452.tif 000030
【図６】
2013519452.tif 000031
【図７】
2013519452.tif 000032
【図８】
2013519452.tif 000033
【図９】
2013519452.tif 000034
【図１０】
2013519452.tif 000035
【図１１】
2013519452.tif 000036
【図１２】
2013519452.tif 000037
【図１４】
2013519452.tif 000038
【図１５】
2013519452.tif 000039
【図１６】
2013519452.tif 000040
【図１７】
2013519452.tif 000041
【図１８】
2013519452.tif 000042
【図１９Ｂ】
2013519452.tif 000043
【図２０Ａ】
2013519452.tif 000044
【図２０Ｂ】
2013519452.tif 000045
【図２１Ａ】
2013519452.tif 000046
【図２１Ｂ】
2013519452.tif 000047
【図２４】
2013519452.tif 000048
【図２５】
2013519452.tif 000049
【図２６Ａ】
2013519452.tif 000050
【図２８Ａ】
2013519452.tif 000051
【図２８Ｂ】
2013519452.tif 000052
【図２８Ｃ】
2013519452.tif 000053
【図２９】
2013519452.tif 000054
【図１３】
2013519452.tif 000055
【図１９Ａ】
2013519452.tif 000056
【図２２】
2013519452.tif 000057
【図２３】
2013519452.tif 000058
【図２６Ｂ】
2013519452.tif 000059
【図２７】
2013519452.tif 000060
【図３０】
2013519452.tif 000061
【国際調査報告】
2013519452.tif 000062
2013519452.tif 000063
2013519452.tif 000064
2013519452.tif 000065