(57)【要約】  非一時的コンピュータ可読媒体は、画像解釈方法１００を実行するための命令であって、少なくとも１つの電子プロセッサ２０を有するワークステーション１８によって読み取り可能かつ実行可能な命令を記憶する。方法は、患者のターゲット部分の第１及び第２の画像を共通の画像空間に空間的に位置合わせするステップであって、第１及び第２の画像は、それぞれ異なる画像セッションから得られたものであり、標準取り込み値（ＳＵＶ）単位のピクセル値を有する、ステップと、空間的に位置合わせされた第１及び第２の画像の対応するピクセルのＳＵＶペアを決定するステップと、決定されたＳＵＶペアの２次元散布図を表示するように表示装置２４を制御するステップであって、２次元散布図は、第１の画像に関する第１のＳＵＶ軸及び第２の画像に関する第２のＳＵＶ軸を有する、ステップと、を有する。
2021513054.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも１つの電子プロセッサを有するワークステーションによって読み取り可能かつ実行可能な、画像解釈方法を実行するための命令を記憶した非一時的コンピュータ可読媒体であって、前記画像解釈方法が、
患者のターゲット部分の第１及び第２の画像を共通の画像空間に空間的に位置合わせするステップであって、前記第１の画像及び第２の画像は、それぞれ異なる撮像セッションから取得され、標準取り込み値（ＳＵＶ）単位でのピクセル値を有する、ステップと、
前記空間的に位置合わせされた第１及び第２の画像の対応するピクセルのＳＵＶペアを決定するステップと、
前記決定されたＳＵＶペアの２次元散布図を表示するよう前記表示装置を制御するステップであって、前記２次元散布図は、前記第１の画像に関する第１のＳＵＶ軸と、前記第２の画像に関する第２のＳＵＶ軸とを有する、ステップと、
を有する、非一時的コンピュータ可読媒体。
【請求項２】
前記方法が更に、前記第１の画像と前記第２の画像との間のＳＵＶスケーリングシフトを決定するステップを有する、請求項１に記載の非一時的コンピュータ可読媒体。
【請求項３】
前記方法が更に、前記ＳＵＶスケーリングシフトを表示するステップを有する、請求項２の記載の非一時的コンピュータ可読媒体。
【請求項４】
前記方法が更に、前記第１の画像と前記第２の画像との間の前記ＳＵＶスケーリングシフトを補正するために、前記決定されたＳＵＶスケーリングシフトを用いて前記２次元散布図を調整するステップを有する、請求項２又は３に記載の非一時的コンピュータ可読媒体。
【請求項５】
前記ＳＵＶスケーリングシフトを決定する前記ステップが、前記ＳＵＶスケーリングシフトを決定するために、前記２次元散布図に対して線形回帰解析を実行することを含む、請求項４に記載の非一時的コンピュータ読み取り可能媒体。
【請求項６】
前記線形回帰解析が、前記空間的に位置合わせされた第１及び第２の画像のピクセルｉにわたって合計され又は平均されたＳＵＶペア間の二乗距離を最小化するようにｍを調整することを含み、ここでｍはＳＵＶスケーリングシフトを表す、請求項５に記載の非一時的コンピュータ可読媒体。
【請求項７】
線形回帰解析が、
2021513054.tif 000010
をｍについて解くことによって実行され、ここでｘｉ及びｙｉは、ピクセルｉについてＳＵＶペアのＳＵＶ値を表し、ｍはＳＵＶのスケーリングシフトを表す、請求項５に記載の非一時的コンピュータ可読媒体。
【請求項８】
前記線形回帰解析が、前記２Ｄ散布図における各ＳＵＶペアから、前記空間的に位置合わせされた第１及び第２の画像のピクセルｉにわたって合計され又は平均された傾きｍを有するラインまでの距離を最小化するようにｍを調整し、ここで、ｍはＳＵＶスケーリングシフトを表す、請求項５に記載の非一時的コンピュータ可読媒体。
【請求項９】
前記線形回帰解析が、
2021513054.tif 000011
をｍについて解くことによって実行され、ここで、ｘｉ及びｙｉは、ピクセルｉについてＳＵＶペアのＳＵＶ値を表し、ｍはＳＵＶのスケーリングシフトを表す、請求項５に記載の非一時的コンピュータ可読媒体。
【請求項１０】
前記方法が更に、ユーザ入力装置を介して前記２次元散布図の選択を受信するステップと前記２次元散布図の選択された部分のＳＵＶペアの診断プロットを表示するステップを有する、請求項１乃至９のいずれか１項に記載の非一時的コンピュータ可読媒体。
【請求項１１】
前記２次元散布図の選択部分を受信する前記ステップは、（ｉ）表示された２次元散布図の領域の境界線を受信するステップと、（ｉｉ）選択基準を定義するクエリを受信するステップと、のうち１つを含む、請求項１０に記載の非一時的コンピュータ可読媒体。
【請求項１２】
前記方法が更に、前記空間的に位置合わせされた第１及び第２の画像の軸方向スライスの関数として、前記２次元散布図の選択された部分のＳＵＶペアのヒストグラムを生成するステップを有する、請求項１０又は１１に記載の非一時的コンピュータ可読媒体。
【請求項１３】
前記第１の画像及び前記第２の画像が、ＳＵＶ単位でのポジトロンエミッショントモグラフィ画像である、請求項１乃至１２のいずれか１項に記載の非一時的コンピュータ可読媒体。
【請求項１４】
患者のターゲット部分の第１の画像と第２の画像との間のＳＵＶスケーリングシフトを決定する方法であって、前記第１及び第２の画像は、それぞれ異なる画像セッションから得られたものであり、標準取り込み値（ＳＵＶ）単位でのピクセル値を有し、前記方法が、
前記第１の画像及び前記第２の画像を共通の画像空間に空間的に位置合わせするステップと、
前記空間的に位置合わせされた第１及び第２の画像の対応するピクセルのＳＵＶペアを決定するステップと、
前記第１の画像に関する第１のＳＵＶ軸及び前記第２の画像に関する第２のＳＵＶ軸を有する２次元空間において、前記決定されたＳＵＶペアに対して線形回帰解析を実行することにより、前記第１の画像と前記第２の画像との間のＳＵＶスケーリングシフトを決定するステップと、
ｉ）前記ＳＵＶスケーリングシフトを表示装置に表示するステップ、又は（ｉｉ）前記ＳＵＶスケーリングシフトに応じて前記第１の画像又は前記第２の画像のＳＵＶ値をスケーリングすることによって、前記ＳＵＶスケーリングシフトを補正するステップ、の少なくとも一方を含むステップと、
を有する方法。
【請求項１５】
前記線形回帰解析が、前記空間的に位置合わせされた第１の画像及び第２の画像のピクセルｉについて合計され又は平均されたＳＵＶペア間の二乗距離を最小化するようにｍを調整し、ｍはＳＵＶスケーリングシフトを表す、請求項１４に記載の方法。
【請求項１６】
前記線形回帰解析が、
2021513054.tif 000012
をｍについて解くことによって実行され、ここで、ｘｉ及びｙｉは、ピクセルｉに対するＳＵＶペアのＳＵＶ値を表し、ｍは、ＳＵＶスケーリングシフトを表す、請求項１５に記載の方法。
【請求項１７】
前記線形回帰解析が、前記２次元散布図における各ＳＵＶペアから、前記空間的に位置合わせされた第１及び第２の画像のピクセルｉにわたって合計され又は平均された傾きｍを有する直線までの距離を、最小化するようにｍを調整し、ｍがＳＵＶスケーリングシフトを表す、請求項１４に記載の方法。
【請求項１８】
前記線形回帰解析が、
2021513054.tif 000013
をｍについて解くことによって実行され、ここで、ｘｉ及びｙｉは、ピクセルｉについてＳＵＶペアのＳＵＶ値を表し、ｍはＳＵＶのスケーリングシフトを表す、請求項１７に記載の方法。
【請求項１９】
前記方法が更に、前記決定されたＳＵＶペアの２次元散布図を表示するように表示装置を制御するステップであって、前記２次元散布図は、前記第１の画像に関する第１のＳＵＶ軸と、前記第２の画像に関する第２のＳＵＶ軸とを有する、ステップを有する、請求項１４乃至１８のいずれか１項に記載の方法。
【請求項２０】
表示装置と、
少なくとも１つのユーザ入力デバイスと、
少なくとも１つの電子プロセッサと、
を有するシステムであって、前記少なくとも１つの電子プロセッサは、
患者のターゲット部分の第１及び第２の画像を共通の画像空間に空間的に位置合わせするステップであって、前記第１及び第２の画像は、それぞれ異なる画像セッションから取得されたものであり、標準取り込み値（ＳＵＶ）単位のピクセル値を有する、ステップと、
前記空間的に位置合わせされた第１及び第２の画像の対応するピクセルのＳＵＶペアを決定するステップと、
前記第１の画像に関する第１のＳＵＶ軸と前記第２の画像に関する第２のＳＵＶ軸とを有する２次元空間において、前記決定されたＳＵＶペアに対して線形回帰解析を実行することにより、前記第１の画像と前記第２の画像の間のＳＵＶスケーリングシフトを決定するステップと、
前記ＳＵＶスケーリングシフトに従って前記第１の画像又は前記第２の画像のＳＵＶ値をスケーリングすることにより、前記ＳＵＶスケーリングシフトを補正するステップと、
ｉ）前記決定されたＳＵＶペアの２次元散布図であって、前記第１の画像に関する第１のＳＵＶ軸と、前記第２の画像に関する第２のＳＵＶ軸とを有する、２次元散布図と、（ｉｉ）前記ＳＵＶのスケーリングシフトと、を表示するように表示装置を制御するステップと、
を実行するようにプログラムされている、システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般的には、医用イメージング技術、ポジトロンエミッショントモグラフィ（ＰＥＴ）イメージング技術、医用画像解釈技術、画像再構成技術、及び関連技術に関する。
【背景技術】
【０００２】
今日、癌のポジトロンエミッショントモグラフィ（ＰＥＴ）／コンピュータトモグラフィ（ＣＴ）イメージングは、腫瘍の診断及びステージング（病期分類）の標準的なコンポーネントとなっている。更に、治療応答の定量的なモニタリング又は新薬開発の評価ツールとしても重要性が高まっている。
【０００３】
癌治療に対する患者の応答を評価するために、臨床医は、少なくとも２セットの画像（以前の画像と現在の画像）を読み、所見を相関させる。標準取り込み値（standardized uptake values、ＳＵＶ）の使用は、フルデオキシグルコース（ＦＤＧ）ＰＥＴ／ＣＴ腫瘍臨床イメージングでは通例であり、ＳＵＶは、治療に対する患者応答を評価する上で特定の役割を担っている。ＳＵＶは、
2021513054.tif 000003
のように計算されることができる。ここで、ｉはＰＥＴ画像のボクセルのインデックスであり、ｖｉは、画像内のボクセルｉの値（ボクセルｉの組織内の放射性物質活性濃度として、例えばＭＢｑ／ｍＬ又は同等の単位で表され、放射性物質ファントムの較正とピクセルボリュームに基づいて生のピクセル値から計算される）であり、ＳＵＶ値に変換され、Ｄは放射性医薬品の投与量、Ｍは、患者の身体質量（又は体重）であり、ｔは放射性医薬品の投与とＰＥＴ画像データ取得の間の待ち時間、ｔ１／２は、放射性医薬品の半減期である。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
臨床医は画像をスクロールして、対応する前の画像と比較して、現在の画像で病変のＳＵＶが改善されているか、あるいは悪化しているかを判断する。しかし、生物学的効果、患者の準備、トレーサの投与などにより、ＳＵＶには変動性があることが知られている。ＳＵＶの変動性を低減し、制御するための対策がとられてきた。専門学会は標準とガイドラインを作成した（例：R.L. Wahl, H. Jacene, Y. Kasamon, and M.A. Lodge, From RECIST to PERCIST: evolving considerations for PET response criteria in solid tumors, Journal of Nuclear Medicine, vol.50, no.5, 122S-150S, May 2009を参照されたい）。研究者はベストプラクティスを共有している（例えば、P.E. Kinahan and J.W. Fletcher, PET/CT standardized uptake values in clinical practice and assessing response to therapy, Seminars in Ultrasound, CT and MRI, vol.31, no.6, pp.496-505, Dec. 2010を参照されたい）。スキャナベンダーもまた、製品をリリースしている（例えば、Q.Check）。これらの努力にもかかわらず、ＳＵＶの変動性は、実際には依然として懸念されている（例えば、M.A. Lodge, Repeatability of SUV in oncologic 18F-FDG PET, Journal of Nuclear Medicine, vol.58, no.4, pp.523-532, April 2017を参照されたい）。
【０００５】
臨床医は、多くの場合、ＳＵＶの変動性に対処するために基準組織を使用する。大動脈弓部血流活性又は健康な肝臓が基準として広く使用されており、シリアル検査では腫瘍対バックグラウンド比が比較される。この基準組織アプローチは、正常組織における取り込みの安定性を前提としており、比率は、基準組織における変動を明示的に補正する。しかしながら、基準組織の変動性が最近報告されている（例えば、R.R. Boktor, G. Walker, R. Stacey, S. Gledhill, and A.G. Pitman, Reference range for intra-patient variability on blood-pool and liver SUV for 18F-FDG PET, Journal of Nuclear Medicine, vol.54, no.5, pp.677-682, May 2013）、これは、血液プールのケースでは、撮像時間の広がりが大きく、クリアランスに個人差があることに起因することがある（例えば、J.A. Thie, Understanding SUV variability in reference tissue for 18F-FDG PET with a simple measurement model, Journal of Nuclear Medicine, vol.55, no,2, pp.352-352, Feb. 2014を参照されたい）。
【０００６】
ＰＥＴ撮像検査の解釈は、一般的に放射線科医によって実施される。臨床環境において、放射線科医は、現在の放射線検査をレビューし、前回の放射線検査と比較し、前回の放射線検査に関する放射線学的レポートをレビューし、前回の放射線検査との比較を含む現在の放射線検査の臨床所見を提示する放射線レポートを作成してファイルするために、わずか数分又は数十分しか割り当てられないことがある。この作業環境において、放射線検査の臨床学的な品質と効率的なスループットの両方を維持することが重要な課題となっている。
【０００７】
以下は、これらの問題を克服するための新規の及び改善されたシステム及び方法を開示する。
【課題を解決するための手段】
【０００８】
開示された１つの側面において、非一時的コンピュータ可読媒体は、画像解釈方法を実行するための命令であって、少なくとも１つの電子プロセッサを有するワークステーションによって読み取り可能かつ実行可能な命令を記憶する。この方法は、患者のターゲット部分の第１及び第２の画像を共通の画像空間に空間的に位置合わせするステップであって、第１及び第２の画像は、異なる画像セッションから得られ、標準取り込み値（ＳＵＶ）単位のピクセル値を有する、ステップと、空間的に位置合わせされた第１及び第２の画像の対応するピクセルに対するＳＵＶペアを決定するステップと、決定されたＳＵＶペアの２次元（２Ｄ）散布図を表示すように表示装置を制御するステップであって、２次元散布図は、第１の画像に関する第１のＳＵＶ軸及び第２の画像に関する第２のＳＵＶ軸を有する、ステップと、を有する。
【０００９】
別の開示された側面において、異なる画像セッションから得られ、標準取り込み値（ＳＵＶ）単位でのピクセル値を有する患者のターゲット部分の第１及び第２の画像間のＳＵＶスケーリングシフトを決定する方法は、第１の画像及び第２の画像を共通の画像空間に空間的に位置合わせするステップと、空間的に位置合わせされた第１の画像及び第２の画像の対応するピクセルのＳＵＶペアを決定するステップと、第１の画像に関する第１のＳＵＶ軸及び第２の画像に関する第２のＳＵＶ軸を有する２次元（２Ｄ）空間において、決定されたＳＵＶペアについて線形回帰解析を実行することにより、第１の画像と第２の画像との間のＳＵＶスケーリングシフトを決定するステップと、（ｉ）ＳＵＶスケーリングシフトを表示装置上に表示するステップ、又は（ｉｉ）ＳＵＶスケーリングシフトに従って第１の画像又は第２の画像のＳＵＶ値をスケーリングすることによりＳＵＶスケーリングシフトを補正するステップ、の少なくとも一方を含むステップと、を有する。
【００１０】
開示された別の側面において、システムは、表示装置と、少なくとも１つのユーザ入力装置とを有する。少なくとも１つの電子プロセッサは、患者のターゲット部分の第１及び第２の画像を共通の画像空間に空間的に位置合わせするステップであって、第１及び第２の画像は、異なる画像セッションから得られ、標準取り込み値（ＳＵＶ）単位でのピクセル値を有する、ステップと、空間的に位置合わせされた第１及び第２の画像の対応するピクセルに対するＳＵＶペアを決定するステップと、第１の画像に関する第１のＳＵＶ軸及び第２の画像に関する第２のＳＵＶ軸を有する２次元（２Ｄ）空間において、決定されたＳＵＶペアに対して線形回帰解析を実行することにより、第１の画像と第２の画像との間のＳＵＶスケーリングシフトを決定するステップと、前記ＳＵＶスケーリングシフトに従って前記第１の画像又は前記第２の画像のＳＵＶ値をスケーリングすることによって前記ＳＵＶスケーリングシフトを補正するステップと、（ｉ）前記決定されたＳＵＶペアの２次元（２Ｄ）散布図であって、前記第１の画像に関する第１のＳＵＶ軸及び前記第２の画像に関する第２のＳＵＶ軸を有する２次元散布図と、（ｉｉ）前記ＳＵＶスケーリングシフトと、を表示するように表示装置を制御するステップと、を実行するようにプログラムされる。
【００１１】
１つの利点は、現在と前回のＰＥＴイメージング検査のＳＵＶ値を比較するＳＵＶプロットを表示する可視化装置を提供することにあり、それにより、現在の検査と前回の検査の間のＳＵＶ変化を読み取って解析する際に臨床医を支援する。
【００１２】
別の利点は、ＳＵＶスケーリング差補正を計算して適用することにあり、それにより、手動で識別された基準組織を使用してそのようなスケーリングを実行する必要がなくなる。
【００１３】
別の利点は、単一の基準組織における変動性の影響を受けにくく且つ位置合わせ誤差の影響を受けにくい計算されたＳＵＶスケーリングを生成することにある。
【００１４】
別の利点は、２つのイメージングセッション間で補正されるＳＵＶ値を生成することにある。
【００１５】
別の利点は、提案された方法は、個々の異なる装置及びアルゴリズムに起因する系統的なバイアスを補正することができるので、患者のフォローアップ検査が変動性を制御するために同じスキャナで実行されるという制約又は好ましさを低減し又は除去することにある。
【００１６】
別の利点は、従来の線形回帰技術よりもロバストな線形回帰アプローチを提供することにある。
【００１７】
所与の実施形態は、上述の利点のどれも提供せず、又は上述の利点の１つ、２つ、より多く、又はすべてを提供することができ、及び／又は、本開示を読み理解することにより、当業者には明らかになるであろう他の利点を提供することができる。
【図面の簡単な説明】
【００１８】
【図１】１つの側面に従う画像解釈システムを示す図。
【図２】図１のシステムの例示的なフローチャート動作を示す図。
【図３Ａ】図１のシステムによって生成されたデータの例示的なプロットを示す図。
【図３Ｂ】図１のシステムによって生成されたデータの例示的なプロットを示す図。
【図４Ａ】図１のシステムによって生成されたデータのヒストグラムの例を示す図。
【図４Ｂ】図１のシステムによって生成されたデータのヒストグラムの例を示す図。
【図５Ａ】図１のシステムによって生成されたデータの例示的なヒストグラムを示す図。
【図５Ｂ】図１のシステムによって生成されたデータの例示的なヒストグラムを示す図。
【図６】本明細書に開示されている線形回帰検定の結果のプロットを示す図。
【図７】本明細書に開示されている線形回帰検定の結果のプロットを示す図。
【図８】本明細書に開示されている線形回帰検定の結果のプロットを示す図。
【発明を実施するための形態】
【００１９】
本開示は、様々なコンポーネント及びコンポーネントの取り合わせ、並びに様々なステップ及びステップの取り合わせの形をとることができる。図面は、好ましい実施形態を例示する目的のためだけにあり、本開示を限定するものとして解釈されるべきでない。
【００２０】
臨床ＰＥＴにおいて、状態（例えば、腫瘍、転移）が増加しているか又は減少しているかを観察することを主な目的として、複数回にわたって画像を取得することが一般的である。撮影セッション間の比較可能性を提供するために、患者質量、放射性医薬品の投与量、待ち時間、及びおそらく他のファクタのカウントを正規化する標準取り込み値（ＳＵＶ値）を使用することが知られている。実際には、そのような正規化は不完全であり（例えば、想定される放射性医薬品量が患者に投与された実際の投与量と一致しない場合がある、放射性医薬品の活性レベルがその公称値と異なる場合がある、体重測定が誤っている場合がある、待ち時間が公称値と異なる場合がある、及び／又はその他）、更に基準領域（一般には肝臓が視野（ＦＯＶ）内にあるときは肝臓が採用される）おけるＳＵＶ値が参照される。しかしながら、この基準組織の正規化が実行される場合でも、セッション間のＳＵＶ変動性があり得る。更に、撮像セッション間のＳＵＶの変化を評価する場合、通常は、２つのセッションから一致する画像を表示し、視覚的に比較することであるが、これは、表示された強度の臨床医の視覚的知覚に依存するだけでなく、ＳＵＶが大きく変化した各領域を検出する臨床医に依存するので、主観的になる可能性がある。
【００２１】
本明細書に開示された実施形態において、一致する画像は、空間的に位置合わせされ、各ピクセルごとに、「前」及び「後」のＳＵＶペア（ＳＵＶ１、ＳＵＶ２）が集計される。一実施形態において、これらの値は、ｘ座標及びｙ座標としてプロットされ、２次元のＳＵＶ−ＳＵＶ散布図をもたらす。撮像セッション間にＳＵＶの変化がなく、ＳＵＶの誤った較正もない理想的な状況では、２次元ＳＵＶ−ＳＵＶプロットは、傾き＝１の直線であるべきである。一方、ＳＵＶ２＞ＳＵＶ１となる領域がある場合、それらはプロット上で視覚的に観察可能な集まりとして現れるはずである。いくつかのＳＵＶの誤った較正がある場合、これは、１とは異なる「変更されていない」ＳＵＶ値ペアの傾きとして現れるはずである。
【００２２】
一実施形態において、２次元ＳＵＶ−ＳＵＶデータペアが、行列データ構造として生成され、回帰解析が、ＳＵＶシフト補正を決定するために適用される。線形回帰勾配ｍは、シフト補正（シフトがない場合はｍ＝１）である。しかしながら、従来の線形回帰は、空間レジストレーション誤差に過度に敏感であり、望ましくないことに回帰方向に依存することが本明細書において認識される。これを考慮して、誤った位置合わせに対する感度を実質的に低下させ、回帰方向に関して対称的である代替的な線形回帰アプローチが本明細書に開示されている。これらの線形回帰アプローチは、ＳＵＶ解析への例示的なアプリケーションに関して本明細書に開示されているが、ここに開示される線形回帰アプローチは、実験データにラインをフィットさせるために線形回帰が実行される任意の状況において、より一般的に適用可能であることに留意されたい。結果として得られる傾きｍは、シフトを実証するために、２次元ＳＵＶ−ＳＵＶプロット上にプロットされることができ、又は代替的に、１つのデータセットが、シフトに関して補正されることができ、例えばＳＵＶ２←（１／ｍ）＊ＳＵＶ１である。シフト補正ｍは更に、放射線レポートに報告されてもよく、臨床医が利用可能なすべての情報を評価することができるように、シフト補正の有り／無しで定量的な結果がレポートされることができる。
【００２３】
本明細書に開示されている他の実施形態は、ユーザインタフェースに関連する。
【００２４】
この側面において、２次元ＳＵＶ−ＳＵＶプロットが表示される。ユーザは、例えばマウスポインタを使用して集まりを囲むことによって、プロットの領域を選択してもよく、選択されたデータに対して様々な解析情報が生成されることができる。１つのアプローチは、各スライスビンの値が、当該スライスに属する選択された領域内のデータのカウントである複数スライスのヒストグラムをプロットすることである。これにより、選択されたデータに寄与する軸方向領域のスライスピークを有するプロットが生成される。過去及び現在のＰＥＴイメージングセッションからの個々のスライスは、視覚による精査を可能にするために横に並べて表示されることができる。別の提示アプローチは、表示されたＰＥＴ画像内の選択されたデータに属するボクセルを強調表示することである。選択されたデータを含む領域を描出するために、クラスタリング（すなわち接続性）解析が実行されてもよい。クラスタリングされたデータの中心を通る３つの十字切断面（横断方向、矢状方向、及び冠状方向）が表示されてもよい。他の解析も企図される。
【００２５】
本明細書においてＰＥＴイメージングシステムについて記載されているが、開示されるアプローチは、放射性医薬品が患者に投与される他の放出イメージングモダリティ、例えば、シングルフォトンエミッションコンピュータトモグラフィ（ＳＰＥＣＴ）イメージングシステム、ハイブリッドＰＥＴ／ＣＴ又はＳＰＥＣＴ／ＣＴイメージングシステム等においても開示されうる。
【００２６】
図１を参照して、例示的な医用イメージングシステム又は装置１０が示されている。図１に示すように、システム１０は、画像取得装置１２を有する。一例において、画像取得装置１２は、ＰＥＴ／ＣＴイメージングシステムのＰＥＴガントリを含むことができ，ＰＥＴ／ＣＴイメージングシステムは、コンピュータトモグラフィ（ＣＴ）ガントリ１３も含む。他の例において、画像取得装置１２は、ＣＴコンポーネントを含まないスタンドアロンＰＥＴスキャナであってもよい。患者テーブル１４は、患者をＰＥＴガントリ１２又はＣＴガントリ１３の検査領域１６に装填するように配置されている。ＰＥＴガントリ１２は、放射線検出器１７のアレイを有する。
【００２７】
システム１０は更に、少なくとも１つの電子プロセッサ２０、少なくとも１つのユーザ入力装置（例えば、マウス、キーボード、トラックボール、放射線レポートを口述するための口述マイク、及び／又は同様のもの）２２、及び表示装置２４などの典型的なコンポーネントを有するコンピュータ又はワークステーション又は他の電子データ処理装置１８を有する。いくつかの実施形態において、表示装置２４は、コンピュータ１８とは別個のコンポーネントであり得る。一般的な臨床インプリメンテーションにおいて、少なくとも１つの電子データ処理装置１８は、撮像装置コントローラ（例えば、ＰＥＴスキャナコントローラ）として機能する第１の電子データ処理装置１８１と、放射線ワークステーションとして機能する第２の電子データ処理装置１８２とを有する。典型的なワークフローにおいて、放射線技師又は他の医療専門家は、ＰＥＴコントローラ１８１を使用してＰＥＴスキャナ１２を操作してＰＥＴ画像を取得し、ＳＵＶ値をもつ放射線画像又はＰＥＴ画像をＳＵＶ値に変換することを可能にする情報が、ピクチャ・アーカイビング・アンド・コミュニケーション・システム（ＰＡＣＳ）２６に記憶される。なお、ＰＡＣＳは、ＲＩＳ（Radiology Information Information System）等の別の命名法を用いてもよい。
【００２８】
その後、放射線科医は、放射線科ワークステーション１８２を操作することにより、現在のＰＥＴ検査及び前回のＰＥＴ検査からのＰＥＴ画像を（ＰＡＣＳ２６から）取り出し、比較することを含む、ＰＥＴ画像の読影を実行する。例えば、前回のＰＥＴ検査は、化学治療、放射線治療、又は他の腫瘍治療の開始前に実施されたものでありえ、現在のＰＥＴ検査は、そのような治療の後に実施されたものでありうる。別の例として、分割化学治療又は放射線治療中に、前回のＰＥＴ検査と現在のＰＥＴ検査が、進行中の分割治療中の異なる時間に実施されてもよい。図１に示すように、ＰＥＴコントローラ１８１及び放射線学ワークステーション１８２の各々は、１つ又は複数の表示装置２４を有し、例示の放射線ワークステーション１８２は、例示的な２つのディスプレイ２４、例えば、画像を表示するための一方のディスプレイと、ドラフト又は他のテキスト情報の下で放射線レポートを表示するための他方のディスプレイとを有し、表示タスクは、様々なディスプレイ２４の間で別のやり方で分散されることができる。
【００２９】
少なくとも１つの電子プロセッサ２０は、画像解釈方法又はプロセス１００を実行することを含む開示された処理を実行するための命令であって、少なくとも１つの電子プロセッサ２０によって読み取り可能かつ実行可能な命令を記憶する１つ又は複数の非一時的記憶媒体（図示せず；例えば、磁気ディスク、ＲＡＩＤ、又は他の磁気記憶媒体；ソリッドステートドライブ、フラッシュドライブ、電子的に消去可能な読み取り専用メモリ（ＥＥＲＯＭ）又は他の電子メモリ；光ディスク又は他の光記憶媒体；それらの様々な組み合わせ；など）と動作可能に接続されている。いくつかの例において、画像解釈方法又はプロセス１００は、放射線科ワークステーション１８２を操作する放射線科医によって実行され、少なくとも部分的にクラウド処理によって実行されることができる。
【００３０】
図２を参照して、画像解釈方法１００の例示的な実施形態がフローチャートとして図示されている。図２に描かれたプロセスの開始に先立って、画像取得装置１２（例えば、ＰＥＴイメージング装置）は、少なくとも１つの電子プロセッサ２０（具体的には、図１の例示的な例ではＰＥＴコントローラ１８１）によって、ＰＥＴ撮像データを取得する、ＰＥＴ撮像データをＰＥＴ画像に再構成する、例えば身体の質量又は体重（Ｍ）、放射性医薬品の投与量（Ｄ）、及び放射性医薬品の投与とＰＥＴ撮像データの取得との間の待ち時間（ｔ）を一般に含む正規化情報を考慮に入れる上記式（１）を使用して、ボクセル値をＳＵＶ値に変換する、ように構成され又は制御される。これは、現在のＰＥＴイメージング検査のために行われ、前回のＰＥＴイメージング検査のために行われ、前回のＰＥＴイメージング検査及び現在のＰＥＴイメージング検査は、ＰＡＣＳ２６に記憶される。１０２において、少なくとも１つの電子プロセッサ２０（より具体的に、説明的な例では放射線ワークステーション１８２）は、（ＰＡＣＳ２６から）患者のターゲット部分の第１及び第２の画像（例えば、第１及び第２のＰＥＴ画像）を取り出し、共通の画像空間に空間的に位置合わせするようにプログラムされている。ちょうど議論されたように、第１及び第２の画像は、典型的には、異なるＰＥＴ画像セッションから得られたものであり、ＳＵＶ単位のピクセル値を有する。画像の空間位置合わせは、任意の適切な剛体又は（好ましくは）非剛体の空間位置合わせ技術を採用してもよい。例えば、１つのアプローチにおいて、ユーザは、第１及び第２の画像内の対応するランドマークを手動でラベル付けし、空間変形フィールドが、一方の画像に適用されることにより、一方の画像は、他方の画像と空間的に位置合わせされる。更に、ユーザは、２つの画像内の１つ又は複数の臓器、腫瘍、又は関心のある他の特徴の周りに輪郭を定義することができ、これらは空間的に位置合わせされる。完全に自動化されたアプローチにおいて、ランドマーク及び／又は輪郭が、エッジ及び／又は点検出アルゴリズムを使用して自動的に識別される。また、画像は、明示的な特徴検出を行わずに、画像内容に基づいて自動的に位置合わせされることもできる。これらは単に説明的な例であり、より一般的には、任意の空間位置合わせアルゴリズム又はアルゴリズムの組み合わせが、第１及び第２の画像を空間的に位置合わせするために採用されることができる。
【００３１】
１０４において、少なくとも１つの電子プロセッサ２０は、空間的に位置合わせされた第１の画像及び第２の画像の対応するピクセルのＳＵＶペアを決定するようにプログラムされている。２つの画像が空間的に位置合わせされた状態で、対応するピクセル（又はボクセル）のペアを識別することは、それらが空間的にアラインされているので、容易である。しかしながら、いかなる空間位置合わせアルゴリズムも不完全であり、前回の撮像セッションと現在の撮像セッションとの間の臓器又は腫瘍のサイズ又は形状の変化（例えば、腫瘍の収縮又は成長、膀胱の拡張又は収縮など）、臓器／腫瘍／その他の回転などの交絡因により、第１の画像と第２の画像との間の完全な位置合わせを提供することに失敗することがあることに留意されたい。
【００３２】
１０６において、少なくとも１つの電子プロセッサ２０は、決定されたＳＵＶペアの２次元（２Ｄ）散布図を表示するように表示装置２４を制御するようプログラムされる。この２次元散布図は、第１の画像に関する第１のＳＵＶ軸と、第２の画像に関する第２のＳＵＶ軸とを有する。図３Ａ及び図３Ｂは、そのような２次元ＳＵＶ−ＳＵＶプロットの一例を示す。このように表示することにより、放射線技師は、前回のイメージング検査と今回のイメージング検査との間でＳＵＶ値が大きく変化したかどうかを容易に把握することができる。例えば、変化がない場合は、散布図は、傾きｍ＝１のラインとして並ぶはずである。一方、ＳＵＶ値が概ね増加しているか、又は概ね低下している場合には、点のかなりの割合がこの傾きｍ＝１の線の上又は下に位置し、その方向（上又は下）は、前回と今回の撮像検査の間にＳＵＶが上昇したか、又は低下したかに応じて決定される。ほとんどの場合、ＳＵＶが全体的に増大又は低下したとしても、これらのＳＵＶの変化は、通常、腫瘍又は他の悪性組織におけるものであり、一方、正常組織は、前回の撮像検査と今回の撮像検査の間でＳＵＶの変化をほとんど又は全く示さない可能性が高いことに注意すべきである。その結果、腫瘍又は他の悪性組織のＳＵＶ値が実質的に変化している場合でも、ＳＵＶが変化していないこれらの領域に対応する強く定義されたｍ＝１の線が、通常はまだ存在するであろう。したがって、２次元ＳＵＶ−ＳＵＶプロットは、通常、変化していないＳＵＶ値の「ベースライン」を画定する「基準」ｍ＝１ラインをなお示す。更に、ＳＵＶスケーリングが、前回と今回のイメージング検査の間で変化している場合、この「基準」ラインは、ＳＵＶスケーリングの変化を定量化するための差であるｍ＝１とは異なる勾配をもつことになる。このすべてのことは、放射線技師が２次元ＳＵＶ−ＳＵＶ散布図を見ることで、一目で理解することができる。
【００３３】
図２に戻り、１０８において、少なくとも１つの電子プロセッサ２０は、第１の画像と第２の画像との間のＳＵＶスケーリングのシフトを決定するようにプログラムされる。いくつかの実施形態において、決定されたＳＵＶスケーリングシフトは、表示装置２４に表示される。
【００３４】
幾つかの実施形態において、少なくとも１つの電子プロセッサ２０は、ＳＵＶスケーリングシフトを決定するために、２Ｄ散布図上で線形回帰解析を実行するようにプログラムされる。一例において、線形回帰解析は、空間的に位置合わせされた第１及び第２の画像シフトのピクセル「ｉ」に亘って回帰、合計、又は平均されるべきラインと、ペアにされたＳＵＶ座標の間の二乗距離を最小化するように、値「ｍ」を調整する。これは、ＳＵＶスケーリングシフト（「ｍ」で表される）について、式（１）を解くことによって実行することができる。
2021513054.tif 000004
ここで、ｘｉ及びｙｉは、ピクセルｉに関するＳＵＶペアのＳＵＶ値を表す。
【００３５】
別の例において、線形回帰解析は、各ＳＵＶペアから、空間的に位置合わせされた第１及び第２の画像のピクセルｉにわたって合計又は平均された勾配ｍを有する直線までの、２Ｄ散布図における結合された残差距離を最小化するようｍを調整する。これは、ＳＵＶスケーリングシフト（「ｍ」で表される）について、式（２）を解くことによって実行されることができる。
2021513054.tif 000005
ここで、ｘｉ及びｙｉは、ピクセルｉについてＳＵＶペアのＳＵＶ値を表す。本明細書に開示されるように、式（１）及び（２）で示される線形回帰アプローチは、従来の線形回帰アプローチと比較して、第１及び第２の画像の位置合わせにおける誤差に対してよりロバストである。
【００３６】
１１０において、少なくとも１つの電子プロセッサ２０は、第１の画像と第２の画像との間のＳＵＶスケーリングシフトを補正するために、決定されたＳＵＶスケーリングシフトにより２Ｄ散布図を調整し又は補正するように、プログラムされる。これは、例えば、第２のＳＵＶ値（ｙｉ）のＳＵＶスケーリングに一致するように、第１のＳＵＶ値（ｘｉ）を係数ｍでスケーリングすることによって実行されることができる。代替として、これは、第１のＳＵＶ値（ｘｉ）のＳＵＶスケーリングに一致するように、第２のＳＵＶ値（ｙｉ）を係数（１／ｍ）でスケーリングすることによって実行されることができる。
【００３７】
１１２において、少なくとも１つの電子プロセッサ２０は、表示された２Ｄ散布図から情報を決定するようにプログラムされる。そのために、少なくとも１つの電子プロセッサ２０は、ユーザ入力装置２２を介して２Ｄ散布図の一部の選択を受信するようにプログラムされる。選択は、ユーザ入力装置２２を介して、表示された２次元散布図の領域の境界線を受信すること、又はユーザ入力装置を介して、選択基準を規定するクエリを受信すること、を含むことができる。例えば、クエリは、ＳＵＶ２がＳＵＶ１よりも少なくとも２０％高いすべてのペアの選択を要求することができる。少なくとも１つの電子プロセッサ２０は、２次元散布図の選択された部分のＳＵＶペアの診断プロットを表示するように表示装置２４を制御するようにプログラムされる。いくつかの例において、少なくとも１つの電子プロセッサ２０は、空間的に位置合わせされた第１及び第２の画像の軸方向スライスの関数として、２Ｄ散布図の選択された部分のＳＵＶペアのヒストグラムを生成するようにプログラムされる。表示された診断プロットは、ヒストグラムを含む。
【００３８】
例
以下、処理１０２−１１２の幾つかの例について、より詳細に説明する。２つのＰＥＴ画像１０２は、同じ空間座標系に位置合わせされる。位置合わせは、剛体であっても非剛体であってもよい。ＰＥＴ画像は、直接的に位置合わせされることもできるし、関連する２つのＣＴ画像を最初に位置合わせすることによって、間接的に位置合わせされることもできる（同じ検査のＰＥＴ及びＣＴは、同じ座標空間にある）。位置合わせには、ボリューム全体又はユーザ定義のサブボリューム（例えば、関心のあるボリューム）を使用することができる。
【００３９】
画像が位置合わせされた後、画像の差分又は比率を計算して、変化を強調することができる。しかし、ここでは、変化は、処理１０４及び１０６において、２次元散布図又はグラフで可視化される。２Ｄグラフは、視覚化が容易であり、差又は比は、２Ｄグラフ上でなお評価されることができ、シリアル検査（すなわち、前の画像と現在の画像を比較すること）におけるＳＵＶスケーリング差が評価されることができる。
【００４０】
図３Ａ及び図３Ｂは、位置合わせ後の同じ空間位置にある２つのＰＥＴ画像からのＳＵＶ値を示している。図３Ａ及び図３Ｂに示すプロットは、２つのＰＥＴボリュームにわたるＳＵＶ及びそれらの関係を単一の２Ｄグラフにまとめたものである。それらのプロットを作成する際、プロットを生成するためのデータ量は、粗い画像の使用又はボクセルグリッドのサブサンプリングによって、任意に低減される。図３Ａ及び図３Ｂは、同じ２次元ＳＵＶ−ＳＵＶ散布図を示しており、以下に本明細書で説明するように、重畳されるラインに関して異なるだけである。
【００４１】
図３Ａ及び３Ｂにおいて、「１」とラベル付けされたラインは、ＳＵＶに変化がない場所を表している。この例において、「１」とラベル付けされたラインは、勾配ｍ＝１を有するが、より一般的には、前の画像と現在の画像のＳＵＶのスケーリングの違いに基づいて、異なる勾配を有してもよい（ただし、これは、勾配ｍ＝１を回復するために、本明細書で開示されるように補正されることができる）。図３Ａにおいて、「２」とラベル付けされた線の上のすべてのドットは、ＳＵＶ２≧ＳＵＶ１＋αである場所を示し、ここでαはユーザが較正可能なパラメータであり、ここでは０．５に設定される；「３」とラベル付けされた線の下のすべてのドットは、ＳＵＶ２≦ＳＵＶ１−αである場所を示している。第２及び第３のライン並びにαは、従来の差分画像と同様の目的を果たし、それらは、ＳＵＢが悪くなる又は良くなる領域を示す。第２のラインの上のドットは、そのロケーションではＳＵＶが悪くなることを示し、第３のラインの下のドットは、そのロケーションではＳＵＶが改善することを示している。（注、「より悪い」及び「改善する」という語は、ＳＵＶ２が現在のＰＥＴ画像のＳＵＶ値であり、ＳＵＶ１が前のＰＥＴ画像のＳＵＶ値であるという慣例を想定する）。
【００４２】
図３Ｂにおいて、同様に、第２のラインより上のすべてのドットは、ＳＵＶ２≧（１＋β）×ＳＵＶ１を示し、ここで、βはユーザが設定可能なパラメータであり、ここでは０．１に設定される。第３のラインより下のすべてのドットは、ＳＵＶ２≦（１−β）×ＳＵＶ１である場所を示している。第２及び第３のライン並びにβは、ＳＵＶ値が悪くなる領域又は良くなる領域を示している。
【００４３】
ユーザは、さらなる解析のために、２次元のＳＵＶ散布図に描かれた特定のデータ部分を選択することができる。一例において、ユーザは、２Ｄグラフからデータの一部を直接的に選択することができ、システムは、いくつかのデータ解析を実行する。別の例として、ユーザは或る数値選択ステートメント（例えば、「ＳＵＶ２＞ＳＵＶ１＋０．５＆ＳＵＶ２＞２．５」）を記述することができ、電子プロセッサ２０は、基準を満たすデータを抽出し、それらに対していくつかの解析を実行する。
【００４４】
これらの解析を実行するために、データの選択又はクエリが必要である。一例において、データ選択は、２次元ＳＵＶ−ＳＵＶプロット上でピッキングし又は描画することによって、直接行うことができる。別の例において、データ選択は、単純な選択ステートメント（例えば、「ＳＵＶ２＞ＳＵＶ１＋α＆ＳＵＶ２＞μ」は、ＳＵＶが悪くなる場所を示すことができ、「ＳＵＶ２＜ＳＵＶ１−α＆ＳＵＶ１＞μ」は、ＳＵＶが改善される場所を示すことができ、ここでμは閾値であり、例えば２．５に設定される）により実行されることができる。
【００４５】
いくつかの例において、データクエリによって指定された通りにデータポイントが抽出され、ヒストグラム解析のような何らかの解析が実行される。
【００４６】
図４Ａ及び図４Ｂは、ヒストグラムの例を示している。図４Ａは、ＳＵＶが悪くなるヒストグラムを示している。これらのデータポイントは、その画像スライスインデックスによって集められる。「４」とラベル付けされたピークは、膀胱がある位置に対応する。「５」とラベル付けされたピークは心臓に対応する。ユーザ選択（例えばヒストグラムのピークをクリックするなどして）されると、システムは、両方のＰＥＴ検査のスライスを表示し、それにより、医師がレビューを行い、臨床的な判断を下すことができる。
【００４７】
図４Ｂは、ＳＵＶが改善されたヒストグラムを示している。ここでも、ユーザがヒストグラムのピークをクリックすると、システムは両方のＰＥＴ検査のスライスを表示すことができ、それにより、医師は、レビューを行い、臨床的な推論を行うことができる。
【００４８】
図３Ａ及び図３Ｂに示されている２次元ＳＵＶ−ＳＵＶプロットから、この例示的な例において、２つの集まりがあることが明らかである。ＳＵＶ１は低い（約０．５）が、ＳＵＶ２は高い；ＳＵＶ１は約２．５であり、ＳＵＶ２は約３である。ユーザは、２つの集まりの中のデータを選択することができ、システムはいくつかの解析を実行する。図５Ａ及び図５Ｂは、この解析の可能な結果を示す。図５Ａは、どのスライスからボクセルが第１の集まり（心臓に関連する）を形成しているかを示している。図５Ｂは、どのスライスからのボクセルが第２の集まり（膀胱に関連する）を形成しているかを示している。ヒストグラムのそれらのピークをクリックすると、システムは、例えば、マルチプレーンリフォーマット（ＭＰＲ）画像を含む関連するスライスを表示することができ、臨床医は適切な判断を下すことができる。
【００４９】
いくつかの実施形態において、ヒストグラム解析は、それらがどこにあるかを大まかに示しているが、ＳＵＶが悪化するデータポイントは、それらのロケーションをピンポイントで特定するために、更にクラスタ化されることができる。例えば、ＳＵＶが悪化するボクセルは、より大きなクラスタを形成するように接続される。小さなクラスタ、例えば、１つのボクセルのみを持つクラスタは、任意には無視することができる。ＳＵＶが悪くなる位置は、２値ボリュームを形成する。分水嶺などのセグメンテーションツールを使用して、それらを、それぞれ異なる関心ボリュームにクラスタ化することができる。それらの関心ボクセルの重心が計算される。それらの重心を横切るＭＰＲ平面が、表示のために生成されて、臨床医は、ＳＵＶの変化を評価することができる。
【００５０】
ＳＵＶ１の値とＳＵＶ２の値が異なるＳＵＶスケーリングを有する場合、データの傾きはｍ＝１から逸脱することが予想される。ＳＵＶスケーリングの差（もしあれば）を決定するために、ＳＵＶ−ＳＵＶ関係について回帰解析が実行されることができる。任意には、ＳＵＶが悪くなっているか又は良くなっている外れ値を除外した後に、回帰解析が実行されることができる。例えば、心臓領域及び膀胱領域のデータは解析から除外されることができる。臨床医は、ＳＵＶ−ＳＵＶプロットに基づく回帰解析から追加の領域を除外することができる。
【００５１】
ＳＵＶ−ＳＵＶ関係は、切片のない線形回帰（すなわち、ＳＵＶ２＝ｍＳＵＶ１）を用いてフィットされ、ここで、ｍはスケーリング補正係数である。しかし、従来の線形回帰は、このアプリケーションではいくつかの困難に苦しむことが本明細書において認識されている。
【００５２】
従来の線形回帰は位置合わせ誤差に敏感であり、更に従来の線形回帰の結果は、どのＳＵＶが独立変数として選択されるかに依存する。これらの問題を解決するために、本明細書に開示されているよりロバストな線形回帰アプローチにおいて、組み合わされた（すなわち、平均された）二乗残差が、ｘ方向及びｙ方向の両方で最小化される。ペアのＳＵＶ座標から回帰ラインまでの二乗距離を最小化することにより、式１が得られる：
2021513054.tif 000006
ｘ方向とｙ方向の両方の結合した二乗残差を最小化すると、式２が得られる。
2021513054.tif 000007
画像の空間位置合わせ誤差に対する様々な線形回帰法の感度を調べるために、同じ画像取得からのものであるが、それぞれ異なるイベント数をもつ異なる２つの画像（フルドーズ（full-dose）画像とロードーズ（low-dose）画像と呼ばれる）が再構成された。ロードーズ画像は、フルドーズ画像の１／１０倍のイベントを用いて再構成された。位置合わせの影響を検討するために、１つの画像を−４０から４０ｍｍの範囲内で水平方向に２ｍｍのステップでシフトさせた。フィットされた従来の回帰直線（独立変数としてＳＵＶ１を使用）は、ｍ＝０．６４７０（ゼロ切片を持つ）の傾きを得た。このフィッティングは、特定の条件、例えば、誤った位置合わせ誤差が２０ｍｍのような条件で行われる。対照的に、ＳＵＶ２をＳＵＶ１にフィッティングした場合（すなわち、ＳＵＶ２を独立変数として使用）、得られた傾きはｍ＝０．６１７６であった。いずれのフィッティングにおいても、Ｒ２＝０．３９９６であった。このように、独立変数の選択に対する依存性が見られる。更に、得られた傾きは、課される画像シフトを除いて得られるであろうｍ＝１よりもはるかに小さく、位置合わせ誤差が従来の回帰ラインに対しかなり影響していることを示している。
【００５３】
位置合わせ誤差の影響は、−４０ｍｍから４０ｍｍまでの水平方向に位置合わせ誤差をスイープすることにより更に検討され、結果が図６に示されるが、従来の線形回帰が、位置合わせ誤差の影響を強く受けるという結論を支持している。
【００５４】
これらの問題（独立変数の選択に対する依存性、及び空間位置合わせ誤差に対する感度）を修正するために、よりロバストな線形回帰アプローチがここに開示される（式（２）及び（３））。式（２）の線形回帰アプローチは、ｘ方向及びｙ方向の両方において、結合された（すなわち平均）二乗残差を最小化する。式（３）の線形回帰アプローチは、点からフィット線までの距離を最小化する。フィット線までの二乗距離を最小化することは、２次方程式を解くことになる。最小化目的関数は次式である：
2021513054.tif 000008
式（４）をｍについて解くと、２次式（２）になる。位置合わせ誤差の関数としてのフィット線の傾きが、図７に示されている。従来の線形回帰（図６）と比較して、位置合わせ誤差に対する改善は明らかであり、位置合わせ誤差の全範囲にわたって、傾きは０．９４９８から０．９７６５の範囲にあるが、グラウンドトルースは１である（中が白い丸）。傾きが１．１に意図的に変更される場合、フィットされた傾きが同じ図に示される（中が黒い丸を持つ上の曲線）。ＳＵＶの役割が入れ替わった場合、すなわちどちらが独立変数であるかを変えた場合、２つの傾きの積は必ず１になる。
【００５５】
ｘとｙの両方向において結合された二乗残差を最小化することは、解析的な解を持つ４次関数を解くことになる。最小化目的関数は次式である：
2021513054.tif 000009
式（５）をｍについて解くと、４次式（３）になる。位置合わせ誤差の関数としてのフィット線の傾きが、図８に示されている。従来の線形回帰と比較して、位置合わせ誤差に対する改善は顕著である−位置合わせ誤差の全範囲にわたって、勾配は０．９７７６から０．９８５１の範囲にあり、グラウンドトルースは１である（中が白い丸）。傾きが意図的に１．１に変更された場合、フィットした傾きが同じ図に示される（中が黒い丸をもつ上の曲線）。ＳＵＶの役割が入れ替えられると、つまりどちらを独立変数とするかを変えると、２つの傾きの積は必ず１になる。
【００５６】
臨床現場において、プロトコルは変動性の管理と密接に関連していることに注意が必要である。したがって、臨床現場において、誤った位置合わせによる差は、上記の例においてシミュレーションしたものよりもはるかに小さいことが予想される。更に、先に述べたように、線形回帰を実行する前に、ラインに沿った外れ値が除去（すなわち、「剪定」）されることができる。線形回帰の出力は、例示された「ＳＵＶ変化なし」のラインとして、２次元ＳＵＶ散布図上に対角線としてプロットされることができる。
【００５７】
本開示は、好ましい実施形態を参照して説明されてきたが、本発明はこれに限定されるものではない。当業者であれば、先行する詳細な説明を読み理解することにより、変形及び変更が思いつくであろう。本発明は、それら変形及び変更が添付の特許請求の範囲又はその等価なものの範囲内にある限り、そのようなすべての変形及び変更を含むものと解釈されることが意図される。
【図１】
2021513054.tif 000014
【図２】
2021513054.tif 000015
【図３Ａ】
2021513054.tif 000016
【図３Ｂ】
2021513054.tif 000017
【図４Ａ】
2021513054.tif 000018
【図４Ｂ】
2021513054.tif 000019
【図５Ａ】
2021513054.tif 000020
【図５Ｂ】
2021513054.tif 000021
【図６】
2021513054.tif 000022
【図７】
2021513054.tif 000023
【図８】
2021513054.tif 000024
【国際調査報告】
2021513054.tif 000025