(57)【要約】  本開示は、心臓壁（２１４、２１８）の運動をスコアリングするためのシステム（１００）に関する。システム（１００）は、最初に心臓壁（２１４、２１８）の第１の画像（２３０）を、次に心臓壁の第２の画像（２４０）を取得可能な撮像システム（１０２）を含む。プロセッサ（１０８）は、第１の画像（２３０）の輪郭データ第１のセットと、第２の画像（２４０）の輪郭データ第２のセットを識別するために提供される。プロセッサは、輪郭データの第１および第２のセットに依存して、心臓壁（２１４、２１８）のセクションにおける心臓周期的変化を表す少なくとも１つの要素（Ｅｉ）を定義する。各要素（Ｅｉ）を分析し、参照データモデルと比較される少なくとも１つのメトリックを生成し、心臓壁の対応するセクションの動きをスコアリングする。本開示はまた、心臓壁の運動記録方法および非一時的なコンピュータ可読媒体に関する。
2021507433.tif 000002
【特許請求の範囲】
【請求項１】
心臓壁の運動をスコアリングするシステムであって、
１回目に前記心臓壁の第１の画像を取得し、２回目に前記心臓壁の第２の画像を取得するように動作可能な撮像システムと、
プロセッサであって：
前記第１の画像の輪郭データの第１のセットを識別し、
前記第２の画像の輪郭データの第２のセットを識別し、
輪郭データの前記第１および前記第２のセットに依存して、前記心臓壁のセクションにおける心臓の周期的変化を表す少なくとも１つの要素を定義し、
各要素を分析して、少なくとも1つのメトリックを生成し、
前記少なくとも1つのメトリックを参照データモデルと比較して、前記心臓壁の対応する前記セクションの前記運動をスコアリングする
ように構成された前記プロセッサと
を含むシステム。
【請求項２】
前記プロセッサは、輪郭点の複数の対を識別するように構成され、前記輪郭点の各対は、輪郭データの前記第１および前記第２のセット内の対応する特徴に関連する、請求項１に記載のシステム。
【請求項３】
各要素が輪郭点の少なくとも第１および第２の対を含み、各対の前記輪郭点が前記要素の頂点を定義する、請求項２に記載のシステム。
【請求項４】
各要素の前記分析が、各対の前記輪郭点間の距離を決定することと、各要素の前記輪郭点間の平均距離を計算することとを含む、請求項３に記載のシステム。
【請求項５】
各要素が、前記輪郭点の運動の軌跡に対応する対向する第１および第２の側面を含む、請求項２、３または４のいずれか一項に記載のシステム。
【請求項６】
各要素の前記分析が、前記要素の面積を決定することを含む、請求項１乃至５のいずれか一項に記載のシステム。
【請求項７】
各要素の前記分析が、各要素の矩形度を決定することを含む、請求項１乃至６のいずれか一項に記載のシステム。
【請求項８】
各要素の前記分析が、各要素の硬度を決定することを含む、請求項１乃至７のいずれか一項に記載のシステム。
【請求項９】
前記プロセッサは、連続的なスケールでスコアリングを実行するように構成される、請求項１乃至８のいずれか一項に記載のシステム。
【請求項１０】
前記撮像システムが二次元画像を生成し、前記第１及び第２の画像は前記心臓壁の二次元画像を含むように構成されている、請求項１乃至９のいずれか一項に記載のシステム。
【請求項１１】
各要素が多角形を含む、請求項１乃至１０のいずれか一項に記載のシステム。
【請求項１２】
前記プロセッサが、安静状態中およびストレス状態中の前記心臓壁の前記運動をスコアリングするように構成されている、請求項１乃至１１のいずれか一項に記載のシステム。
【請求項１３】
前記第１の画像が収縮末期画像を含み、前記第２の画像が拡張末期画像を含み、輪郭データの前記第１セットが収縮末期輪郭データを含み、輪郭データの前記第２セットが拡張末期データを含む、請求項１乃至１２のいずれか一項に記載のシステム。
【請求項１４】
心臓壁の運動をスコアリングする方法であって、
１回目の前記心臓壁の第１の画像および２回目の前記心臓壁の第２の画像を取得するステップと、
前記第１の画像の輪郭データの第１のセットを識別するステップと、
前記第２の画像の輪郭データの第２のセットを識別するステップと、
輪郭データの前記第１および前記第２のセットに依存して、前記心臓壁のセクションにおける心臓の周期的変化を表す少なくとも１つの要素を定義するステップと、
各要素を分析して、少なくとも1つのメトリックを生成するステップと、
前記少なくとも1つのメトリックを参照データモデルと比較して、前記心臓壁の対応する前記セクションの前記運動をスコアリングするステップと
を含む、方法。
【請求項１５】
輪郭点の複数の対を識別することを含み、輪郭点の各対は、輪郭データの前記第１および前記第２のセットにおける対応する画像要素に関連する、請求項１４に記載の方法。
【請求項１６】
各要素は、輪郭点の少なくとも第１および第２の対を含み、各対の前記輪郭点が前記要素の頂点を定義する、請求項１５に記載の方法。
【請求項１７】
各要素の前記分析が、各対の前記輪郭点間の距離を決定することと、各要素の前記輪郭点間の平均距離を計算することとを含む、請求項１６に記載の方法。
【請求項１８】
各要素が、前記輪郭点の運動の軌跡に対応する対向する第１および第２の側面を含む、請求項１５、１６または１７のいずれか一項に記載の方法。
【請求項１９】
各要素の前記分析が、前記要素の面積を決定することを含む、請求項１４乃至１８のいずれか一項に記載の方法。
【請求項２０】
各要素の前記分析が、各要素の矩形度を決定することを含む、請求項１４乃至１９のいずれか一項に記載の方法。
【請求項２１】
各要素の前記分析が、各要素の硬度を決定することを含む、請求項１４乃至２０のいずれか一項に記載の方法。
【請求項２２】
スコアリングが連続的なスケールで実行される、請求項１４乃至２１のいずれか一項に記載の方法。
【請求項２３】
前記第１および第２の画像は、前記心臓壁の二次元画像を含む、請求項１４乃至２２のいずれか一項に記載の方法。
【請求項２４】
各要素が多角形を含む、請求項１４乃至２３のいずれか一項に記載の方法。
【請求項２５】
安静状態中およびストレス状態中の前記心臓壁の前記運動をスコアリングすることを含む、請求項１４乃至２４のいずれか一項に記載の方法。
【請求項２６】
前記第１の画像が収縮末期画像を含み、前記第２の画像が拡張末期画像を含み、輪郭データの前記第１セットが収縮末期輪郭データを含み、輪郭データの前記第２セットが拡張末期輪郭データを含む、請求項１４乃至２５のいずれか一項に記載の方法。
【請求項２７】
実行時にプロセッサに請求項１４乃至２６のいずれか一項に記載の方法を実行させる一連の命令が格納された非一時的なコンピュータ可読媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、画像を分析するための方法および装置に関する。より具体的には、限定ではないが、本開示は、心臓壁の運動をスコアリングするためのシステムおよび方法に関する。
【背景技術】
【０００２】
２次元（２Ｄ）心臓超音波検査は、心臓の運動を様々な条件下、例えば、安静時やストレス時において評価できる画像処理技術である。分析は、例えば、中間ストレス段階及び／又は回復段階を含む他の条件下で実行されてもよい。これは、機能低下または機能不全の心臓の領域を強調することができ、従って、医療介入が必要となる可能性のある患者を識別することができる。左心室の典型的なモデルは、心臓の異なる２Ｄ画像を使用して表示できる１６個のセグメントを含む。左心室の他のモデルは、例えば、１７個のセグメントを含んでもよい。心尖部下壁セグメント、中部下壁セグメント、基部下壁セグメント、心尖部前壁セグメント、中部前壁セグメント、及び基部前壁セグメントは、心尖部二腔画像で見ることができる。心尖部中隔セグメント、中部中隔セグメント、基部中隔セグメント、心尖部側壁セグメント、中部側壁セグメント、及び基部側壁セグメントは、心尖部四腔画像で見ることができる。前壁中隔セグメント、下壁中隔セグメント、中部下壁セグメント、中部前壁セグメント、前側側壁セグメント、及び下側壁セグメントは、傍胸骨短軸像で見ることができる。心尖部側壁セグメント、中部下側壁セグメント、基部下側壁セグメント、心尖部側壁セグメント、中部中隔セグメント、及び基部中隔セグメントは、心尖部三腔画像（または傍胸骨長軸像）で見ることができる。各セグメントの動作は、左心室のさまざまなセクションで見ることができる。（安静時やストレス時など）さまざまな条件下での心筋の各セグメントの運動は、現在、熟練の心臓専門医による２Ｄ心臓超音波検査データの判読によって決定される。これは、カテゴリ化するという方法に基づいて実行される。例えば、心筋壁の各セクションは、以下の基準壁運動スコアのうちの１つを有するものとして分類され得る：正常（「１」）、運動低下（「２」）、無動（「３」）、運動異常（「４」）、及びスコアリング不能（「Ｘ」）。例えば、各セグメントに５つまたは７つの離散スコアを定義するといった他の分類を使用してもよい。既知の手法は主観的な分類に依存しており、多大な時間を要する作業となり得る。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
画像定量化ツールは、以下を可能にする必要があることが認識されている：（ｉ）基礎疾患の病態生理学の経時的な変化、（ｉｉ）地理的位置および患者集団の性質の変化に伴う疾患病理の変動、（ｉｉｉ）何が病気と定義されているか、何が病気を引き起こしているのかについての理解の変化。
【０００４】
少なくとも特定の実施形態では、本発明は、画像を分析するための改善された方法および装置を提供しようと試みるものである。
【課題を解決するための手段】
【０００５】
本発明の態様は、添付の特許請求の範囲に記載の、心臓壁の運動をスコアリングするためのシステム、心臓壁の運動をスコアリングする方法、および非一時的なコンピュータ可読媒体に関する。
【０００６】
本発明の一態様によれば、心臓壁の運動をスコアリングするためのシステムが提供され、このシステムは、１回目に心臓壁の第１の画像を取得し、２回目に心臓壁の第２の画像を取得するように動作可能な撮像システムと、プロセッサであって、第１の画像の輪郭データの第１のセットを識別し、第２の画像の輪郭データの第２のセットを識別し、輪郭データの第１および第２のセットに依存して心臓壁のセクションにおける心臓の周期的変化を表す少なくとも１つの要素を定義し、各要素を分析して、少なくとも1つのメトリックを生成し、少なくとも1つのメトリックを参照データモデルと比較して、心臓壁の対応するセクションの運動をスコアリングする、プロセッサとを含む。少なくとも特定の実施形態では、本明細書に記載のシステムは、第１および第２の画像に応じて、壁運動のスコアリングを部分的または完全に自動化することができる。さらに、本明細書に記載されているシステムは、カテゴリーモデルではなく、連続的なスコアリングシステムを実装することができる。少なくとも特定の実施形態では、心臓壁の運動を、実質的に連続的なスケールでスコアリングすることができる。スコアリングは、１つのメトリックに依存して、または複数のメトリックに依存して実行され得る。
【０００７】
各要素は、心臓壁のセクションの心臓周期パラメータの変化を表す。要素はそれぞれ少なくとも２つの次元で定義される。変形例では、要素はそれぞれ３次元で定義できる。
【０００８】
少なくとも１つの要素は、心臓壁のセクションの心臓の周期的変化を表す。心臓壁は、複数のセグメントによってモデル化されてもよい。各要素は、セグメントの１つにおける心臓の周期的変化を表してもよい。あるいは、各要素は、セグメントのうちの１つの一部（すなわち、サブセグメント）における心臓の周期的変化を表してもよい。要素は、従来のモデルのセグメントと無関係である場合があることが理解されよう。例えば、要素は心臓壁を８個の領域に細分してもよい。
【０００９】
システムは、心臓の心室の運動をスコアリングしてもよい。例えば、システムは、心臓の左心室の運動をスコアリングしてもよい。本明細書で説明するように、心臓の左心室は、１６個のセグメントによってモデル化されてもよい。各要素は、左心室の内壁のセグメントの１つにおける心臓の周期的変化を表してもよい。あるいは、各要素は、左心室の内壁のセグメントのうちの１つの部分（すなわち、サブセグメント）の心臓の周期的変化を表してもよい。
【００１０】
参照データモデルは事前定義されてもよい。参照データモデルは、例えば、複数の個人に対して取得された複数の画像の分析から導出された、履歴データの分析に依存して生成されてもよい。各要素に対して生成された少なくとも1つのメトリックは、心臓壁の対応するセクションの参照データモデルと比較されてもよい。参照データモデルは、心臓壁のモデルの各セグメントについて生成されてもよい。少なくとも１つのメトリックは、心臓壁の対応するセグメントの参照データモデルと比較されてもよい。
【００１１】
輪郭データの第１のセットは、複数の第１の輪郭点を含んでもよい。第１の輪郭点は、第１の輪郭を形成するために互いに接続されてもよい。第１の輪郭点は、直線及び／又は曲線によって互いに接続されてもよい。第１の輪郭は、第１の連続的な曲線を含んでもよい。
【００１２】
輪郭データの第２のセットは、複数の第２の輪郭点を含んでもよい。第２の輪郭点は、第２の輪郭を形成するために互いに接続されてもよい。第２の輪郭点は、直線及び／又は曲線によって互いに接続されてもよい。第２の輪郭は、第２の連続的な曲線を含んでもよい。
【００１３】
輪郭データの第１のセットは、ユーザ入力、例えば第１および第２の画像内の輪郭点を指定するためのユーザ入力に応じて識別されてもよい。あるいは、プロセッサは、第１および第２の画像を分析して、輪郭データの第１および第２のセットを識別するように構成されてもよい。プロセッサは、例えば、画像処理技術を利用して、輪郭データの第１および第２のセットを識別してもよい。画像処理技術は、例えば、第１および第２の輪郭点を識別するための動作セグメンテーションおよび／または第１および第２の輪郭点および／または第１および第２の輪郭を識別するための境界検出技術を含んでもよい。
【００１４】
プロセッサは、輪郭点の複数の対を識別するように構成されてもよい。輪郭点の各対は、輪郭データの第１および第１のセット内の対応する特徴に関連付けられてもよい。各対は、対応する第１および第２の輪郭点で構成されてもよい。第１および第２の輪郭点は、第１および第２の輪郭データセットに存在する単一の特徴に関連してもよい。
【００１５】
各要素は、輪郭点の少なくとも第１および第２の対を含んでもよい。各対を構成する輪郭点は、要素の頂点を定義してもよい。
【００１６】
各要素の分析には、各対の輪郭点間の距離を決定することと、各要素の輪郭点間の平均距離を計算することが含まれる。各要素の平均距離は、輪郭データの第１又は第２のセットの全周に正規化されてもよい。
【００１７】
各要素は、輪郭点の運動の軌跡に対応する対向する第１および第２の側面を含んでもよい。
【００１８】
各要素の分析は、要素の面積を決定することを含んでもよい。面積は、輪郭データの第１及び第２のセットによって囲まれた面積に対応してもよい。この面積は、例えば、拡張末期の輪郭と収縮末期の輪郭との間の面積に対応してもよい。決定された面積は正規化されてもよい。例えば、各要素の面積は、拡張末期データ点の総面積に正規化されてもよい。
【００１９】
各要素の分析は、各要素の矩形度を決定することを含んでもよい。各要素の矩形度は、各要素の面積とその最小外接矩形の面積との比率として計算されてもよい。
【００２０】
各要素の分析は、各要素の硬度を決定することを含んでもよい。各要素の硬度は、要素の面積と凸包の面積の比率として計算されてもよい。
【００２１】
プロセッサは、連続的なスケールでスコアリングを実行するように構成されてもよい。プロセッサは、スコアを出力して心臓壁運動を分類または等級付けするように構成されてもよい。あるいは、又はさらに、プロセッサは、記録を色として表すためにＲＧＢ値を出力するように構成されてもよい。
【００２２】
撮像システムは、３次元画像を生成するように構成されてもよい。第１および第２の画像は、３次元画像であってもよい。少なくとも１つの要素は、３次元で定義されてもよい。撮像システムは、２次元画像を生成するように構成されてもよい。第１および第２の画像は、心臓壁の２次元画像を含んでもよい。少なくとも１つの要素は、２次元で定義されてもよい。
【００２３】
各要素は、例えば１つ以上の湾曲した側面を含む、不規則な形状を含んでもよい。あるいは、各要素は多角形を含んでもよい。
【００２４】
プロセッサは、少なくとも１つの状態の間に心臓壁の運動をスコアリングするように構成されてもよい。プロセッサは、複数の状態の間に心臓壁の運動をスコアリングするように構成されてもよい。スコアリングは、次のセットの１つ以上の間に実行されてもよい：安静状態、ストレス状態、中程度のストレス状態、回復状態など。撮像システムは、ストレス条件下で心臓壁の第１および第２の画像を取得するように動作可能であってよい。撮像システムは、安静条件下で心臓壁の第１および第２の画像を取得するように動作可能であってよい。スコアリングは、様々な条件下、例えばストレスおよび休息条件下で取得された画像に関して実行されてもよい。
【００２５】
第１の画像は収縮末期画像を含んでもよく、第２の画像は拡張末期画像を含んでもよい。輪郭データの第１のセットは収縮末期輪郭データを含んでもよく、輪郭データの第２のセットは拡張末期データを含んでもよい。
【００２６】
本発明のさらなる態様によれば、心臓壁の運動をスコアリングする方法が提供され、この方法は、１回目に心臓壁の第１の画像を取得し、２回目に心臓壁の第２の画像を取得するステップと、第１の画像の輪郭データの第１のセットを識別するステップと、第２の画像の輪郭データの第２のセットを識別するステップと、輪郭データの第１と第２のセットに依存して心臓壁のセクションにおける心臓の周期的変化を表す少なくとも１つの要素を定義するステップと、各要素を分析して、少なくとも1つのメトリックを生成するステップと、少なくとも１つのメトリックを参照データモデルと比較して、心臓壁の対応するセクションの運動をスコアリングするステップとを含む。
【００２７】
この方法は、輪郭点の複数の対を識別することを含んでもよい。輪郭点の各対は、輪郭データの第１および第２のセットの対応する特徴に関連してもよい。
【００２８】
各要素は、輪郭点の少なくとも第１及び第２の対を含んでもよい。各対の輪郭点は、要素の頂点を定義してもよい。
【００２９】
各要素の分析は、各対の輪郭点間の距離を決定することと、各要素の輪郭点間の平均距離を計算することを含んでもよい。
【００３０】
この方法は、各輪郭点の運動の軌跡を決定することを含んでもよい。運動の軌跡は、例えば、各対の輪郭点の位置を比較することにより、第１および第２の画像を参照して決定してもよい。あるいは、またはさらに、運動の軌跡は、第１と第２の時間の間に取得された１つまたは複数の中間画像に応じて生成されてもよい。各要素は、画像要素の運動の軌跡に対応する対向する第１および第２の側面を含んでもよい。
【００３１】
各要素の分析は、要素の面積を決定することを含んでもよい。面積は、輪郭データの第１及び第２のセットによって囲まれた面積に対応してもよい。この面積は、例えば、拡張末期の輪郭と収縮末期の輪郭との間の面積に対応してもよい。決定された面積は正規化されてもよい。例えば、各要素の面積は、拡張末期データ点の総面積に正規化されてもよい。
【００３２】
各要素の分析は、各要素の矩形度を決定することを含んでもよい。各要素の矩形度は、各要素の面積とその最小外接矩形の面積との比率として計算されてもよい。
【００３３】
各要素の分析は、各要素の硬度を決定することを含んでもよい。各要素の硬度は、要素の面積と凸包の面積の比率として計算されてもよい。
【００３４】
記録は、連続的なスケールで実行され得る。
【００３５】
第１および第２の画像は、心臓壁の３次元画像を含み得る。第１および第２の画像は、心臓壁の２次元画像を含み得る。
【００３６】
各要素は多角形を含んでもよい。多角形は２次元または３次元で定義してもよい。
【００３７】
この方法は、安静状態中およびストレス状態中の心臓壁の運動をスコアリングすることを含んでもよい。
【００３８】
第１の画像は収縮末期画像を含み得、第２の画像は拡張末期画像を含んでもよい。輪郭データの第１のセットは収縮末期輪郭データを含んでもよく、輪郭データの第２のセットは拡張末期データを含んでもよい。
【００３９】
本発明のさらなる態様によれば、参照データモデルを生成するためのシステムが提供され、このシステムは、以下のように構成されたプロセッサを備える：複数セットの心臓超音波検査画像を分析し、各セットは少なくとも第１及び第２の心臓超音波検査画像を含み、分析は各セット内の収縮末期画像および拡張末期画像を識別することを含み、各セットの収縮末期の画像と拡張末期の画像を比較して、心臓の周期的変化を特定し、各セットに関して識別された心臓の周期的変化に応じて少なくとも1つのメトリックを計算し、心臓超音波検査画像のセットの少なくともいくつかに関連する結果データを取得し、計算されたメトリックと関連する結果データに応じて参照データモデルをコンパイルする。
【００４０】
本発明のさらなる態様によれば、参照データモデルを生成する方法が提供され、この方法は、複数セットの心臓超音波検査画像を分析するステップであって、各セットは少なくとも第１及び第２の心臓超音波検査画像を含み、分析は各セット内の収縮末期画像および拡張末期画像を識別することを含む、ステップと、各セットの収縮末期の画像と拡張末期の画像を比較して、心臓の周期的変化を特定するステップと、各セットに関して識別された心臓の周期的変化に応じて少なくとも1つのメトリックを計算するステップと、心臓超音波検査画像のセットの少なくともいくつかに関連する結果データを取得するステップと、計算されたメトリックと関連する結果データに応じて参照データモデルをコンパイルするステップとを含む。
【００４１】
本発明のさらなる態様によれば、実行時にプロセッサに本明細書に記載された方法を実行させる一連の命令が格納された非一時的なコンピュータ可読媒体が提供される。
【００４２】
本明細書で説明される任意の制御ユニットまたはコントローラは、１つまたは複数の電子プロセッサを有する計算デバイスを適切に含んでもよい。システムは、単一の制御ユニットまたは電子コントローラを備えることができ、あるいは、コントローラの異なる機能を、異なる制御ユニットまたはコントローラに具体化またはホストすることができる。本明細書で使用する場合、「コントローラ」または「制御ユニット」という用語は、単一の制御ユニットまたはコントローラと、述べられた制御機能を提供するために集合的に動作する複数の制御ユニットまたはコントローラの両方を含むと理解される。コントローラまたは制御ユニットを構成するために、実行されると、制御ユニットまたは計算デバイスに本明細書で指定された制御技術を実施させる適切な命令のセットが提供されてもよい。命令のセットは１つまたは複数の電子プロセッサに適切に埋め込まれてもよい。あるいは、命令のセットは、計算装置で実行される前記コントローラに関連する１つまたは複数のメモリに保存されたソフトウェアとして提供されてもよい。制御ユニットまたはコントローラは、１つまたは複数のプロセッサ上で実行されるソフトウェアで実装され得る。１つまたは複数の他の制御ユニットまたはコントローラは、１つまたは複数のプロセッサ、任意選択で第１のコントローラと同じ１つまたは複数のプロセッサ上で実行されるソフトウェアで実装することができる。他の適切な配置も使用できる
【００４３】
本出願の範囲内で、前の段落、特許請求の範囲および／または以下の説明および図面、ならびに特にそれらの個々の特徴に記載されたさまざまな態様、実施形態、例および代替案が、単独で、または任意の組み合わせで撮影され得る。すなわち、すべての実施形態および／または任意の実施形態の特徴は、そのような特徴に互換性がない限り、任意の方法および／または組み合わせで組み合わせることができる。申請者は、最初に提出された特許請求の範囲を変更するか、または新しい特許請求の範囲を提出する権利を留保する。
【図面の簡単な説明】
【００４４】
次に、本発明の１つまたは複数の実施形態を、例としてのみ、添付の図面を参照して説明する。
【図１】図１は、本発明の実施形態による心臓超音波検査システムの概略図である。
【図２】図２は、心臓の四腔断層像を概略的に示す図である。
【図３】図３は、図１のシステムによって実行される診断方法の主なステップを示すフロー図である。
【図４】図４Ａは、図１に示されている心臓超音波検査システムによって取得された収縮末期の画像を示し、図４Ｂは、図１に示す心臓超音波検査システムによって取得された拡張末期の画像を示す。
【図５】図５Ａは、収縮末期の輪郭点と拡張末期の輪郭点で構成される輪郭データセットを示し、図５Ｂは、図５Ａに示される等高線データセットから生成された複数の要素を示す。
【図６】図６Ａは、図５Ｂに示されている各要素の面積メトリックの生成を示し、図６Ｂは、図５Ｂに示されている各要素の平均距離メトリックの生成を示し、図６Ｃは、図５Ｂに示されている各要素の矩形度メトリックの生成を示し、図６Ｄは、図５Ｂに示されている各要素の硬度メトリックの生成を示しています。
【図７】図７Ａは、安静状態の面積メトリックの分析を示し、図７Ｂは、安静状態の平均距離メトリックの分析を示し、図７Ｃは、安静状態の矩形度メトリックの分析を示し、図７Ｄは、安静状態の硬度メトリックの分析を示す。
【図８】図８Ａは、ストレス状態の面積メトリックの分析を示し、図８Ｂは、ストレス状態の距離メトリックの分析を示し、図８Ｃは、ストレス状態の矩形度メトリックの分析を示し、図８Ｄは、ストレス条件の硬度メトリックの分析を示す。
【図９】図９Ａは、正規化された面積および平均距離メトリックの二変量解析に基づく参照データモデルを示し、図９Ｂは、図９Ａに示されている参照データモデルに適用された正規分布データセットを示す。
【図１０】図１０Ａは、二腔心尖部画像における第１の安静状態の各メトリックの多変量解析を示し、図１０Ｂは、四腔心尖部画像における第２の安静状態の各メトリックの多変量解析を示す。
【図１１】図１１Ａは、二腔心尖部画像における第１のストレス状態の各メトリックの多変量解析を示し、図１１Ｂは、四腔心尖部画像における第２のストレス状態の各メトリックの多変量解析を示す。
【図１２】図１２Ａは、安静状態の要素に適用された記録を示し、図１２Ｂは、ストレス状態の要素に適用された記録を示す。
【図１３】図１３は、しきい値を定義する一連の決定点で構成される決定木を示す。
【発明を実施するための形態】
【００４５】
次に、本発明の実施形態による心臓超音波検査システム１００を、添付の図面を参照して説明する。心臓超音波検査システム１００は、心臓２００の画像を分析し、心臓の周期運動をスコアリングするように動作可能である。
【００４６】
図１に示されるように、心臓超音波検査システム１００は、患者１０４の身体の近くに、典型的には可能な限り心臓の近くとなるように配置されたトランスデューサアレイ１０２、デジタル電子プロセッサであり得るプロセッサ１０８とハードディスクなどのメモリ１１０とを含む処理ユニット１０６、およびフラットスクリーンモニタまたはＬＥＤディスプレイなどのディスプレイ１１２を備える。システムはさらに、ユーザ入力装置、例えば、ディスプレイ１１２に統合されたタッチスクリーン１１４を含んでもよく、これは、ユーザが心臓超音波検査システム１００に入力を提供することを可能にするユーザ入力を提供する。マウス、タッチパッドまたはキーボードなどの他のユーザ入力は、言うまでもなく使用可能である。処理ユニット１０６は、トランスデューサアレイ１０２に接続され、トランスデューサアレイがフェーズドアレイとして一連のパルスが患者を透過してスキャンする超音波ビームを放射し、各パルスによって心臓からの反射超音波を検出するようにトランスデューサアレイを制御するように構成される。心臓の１回のスキャンで単一の画像が作成され、スキャンは通常１秒あたり２５〜５０枚の画像で繰り返され、心臓周期中の運動を示す心臓のリアルタイムビデオ画像が作成される。各画像は、例えば、画像を構成する各ピクセルの輝度値を含み得る画像データセットとしてメモリ１１０に記憶されてもよい。本明細書ではシステムを一般的な用語で説明しているが、適切な心臓超音波検査システムには、例えば、フィリップスエピックｉＥ３３、ＧＥビビッドｅ９、またはフィリップスＣＸ５０などのポータブルシステム、またはハンドヘルドシステムが含まれる。
【００４７】
心臓超音波検査のプロセスはよく知られており、本明細書では詳細に説明しない。いくつかの異なる撮像方法があるが、本実施形態による心臓超音波検査システム１００は、２次元撮像を使用する。心臓から複数の異なる平面に画像が提供されることが知られており、この画像は、心臓の４つの主な心室である、左心室（ＬＶ）、右心室（ＲＶ）、左心房（ＬＡ）および右心房（ＲＡ）のさまざまな側面を示す。そのような画像には、例えば、心尖部四腔画像、心尖部二腔画像、心尖部三腔画像、および胸骨傍の長軸および短軸画像が含まれる。いずれの場合も、単一の静止画像を取得できるが、通常、心臓の周期にわたって一連の画像が取得されるため、その動作を記録して分析することができる。心臓超音波検査システム１００は、前述の画像の１つまたは複数を利用して、心臓２００の心臓周期運動を記録することができる。
【００４８】
例として、心臓２００の四腔心尖部画像が図２に示されている。画像は、心臓２００の２Ｄ平面を含む。画像は、左心室（ＬＶ）２０２、右心室（ＲＶ）２０４、左心房２０６、右心房２０８、および中隔２１０を示す。心尖部２１２、側方左心室２０２の側壁２１４、基部２１６、および内壁２１８もまた、四腔心尖部画像で見ることができる。左心室２０２の長手方向軸２２０は、心尖部２１２を通って延びる。左心室２０２は、長手方向軸２２０の対向する側に配置された第１および第２の側面２２２、２２４を有する。
【００４９】
処理ユニット１０６は、本明細書に記載されているスコアリング技術を実施するために四腔心尖部画像を分析する。あるいは、またはさらに、処理ユニット１０６は、以下の１つ以上を利用してもよい：二腔心尖部画像、胸骨傍短軸画像および三腔心尖部画像。本明細書に記載されたスコアリング技術を実施するために、他の心臓超音波検査画像が処理ユニット１０６によって使用されてもよい。処理ユニット１０６は、心臓超音波検査画像の様々な組み合わせを使用して、左心室の１６のセグメントに対する記録を提供することができる。処理ユニット１０６は、複数の画像を分析し、複数の画像の測定基準に応じて心臓周期運動にスコア付けすることができる。処理ユニット１０６は、有効な画像を定性的に評価し、心臓壁の特定のセクションの心臓周期運動のより明確な表現を提供すると決定された画像に優先順位を付けることができる。直交座標系は、左心室２０２の心尖部２１２を通って伸び、その縦軸に沿って伸びる垂直軸（ここではｙ軸と呼ばれる）と、心尖部２１２と基部２１６の中間の左心室２０２の中点を通る水平軸（ここではｘ軸と呼ばれる）を含むと定義される。
【００５０】
心臓超音波検査システム１００の動作を表すブロック図を図３に示す。心臓超音波検査システム１００は、一連の２Ｄ画像を取得し、それらをメモリ１１０（ブロック３００）に記憶するように構成される。画像は、単一の心臓周期にわたって取得されてもよく、例えば、１つの周期をカバーする１０から５０の画像を含んでもよい。心臓超音波検査システム１００は、単一のスキャンまたは複数のスキャンを実行することができる。例えば、心臓超音波検査システム１００は、第１および第２のスキャンを実行することができる。第１のスキャンは、患者が安静状態にあるときに実行されてもよく、第２のスキャンは、患者がストレス状態にあるときに実行されてもよい。心臓超音波検査システム１００は、例えば、心拍数がストレスを受けた後、心拍数が正常に戻るときの回復段階の間に、休止状態とストレス状態との間で１つ以上の中間スキャンを任意に実行してもよい。画像の取得は、従来の心臓超音波検査システムで実行できる。画像のその後の分析は、図１に示すように心臓超音波検査システムの一部を形成する同じ処理ユニット１０６を使用して実行することができる。ここで、画像は、プロセッサ、メモリ、ユーザ入力、およびディスプレイを有し、この目的のために処理ユニット１０６と同様に動作するラップトップまたはＰＣなどのコンピュータにダウンロードすることができ、画像のさらなる分析は、専用ソフトウェアの制御下でそのコンピュータ上で実行することができる。画像はＰＡＣＳ（画像アーカイブおよび通信システム）から検索されてもよいことが理解されよう。あるいは、またはさらに、画像は、処理のために外部サーバーに送信されてもよい。画像は送信前に匿名化されてもよい。
【００５１】
左心室２０２について、収縮末期に最も近い画像、すなわち心臓周期中の最大収縮、および拡張末期に最も近い画像、すなわち心臓周期中の最大容積が識別される（ブロック３０２）。これは、ユーザがディスプレイ１１２上の画像を見て、収縮末期に最も近いものとして第１の画像２３０（ここでは収縮末期画像２３０と呼ぶ）、および拡張末期に最も近いものとして第２の画像２４０（ここでは拡張末期画像２４０と呼ぶ）を選択することによって実施される。収縮末期画像２３０および拡張末期画像２４０は、それぞれ心臓周期の最初および２回目に取得される。例示的な収縮末期画像２３０が図４Ａに示され、例示的な拡張末期画像２４０が図４Ｂに示されている。収縮末期画像２３０および拡張末期画像２４０の選択は、目視により判断した各画像における左心室２０２の容積の評価および比較に基づいて、または僧帽弁の開閉点に注意することにより、又はＥＣＧプロット上のＱＲＳ複合体の使用により、又はこれらの任意の組み合わせによりユーザが行うことができる。あるいは、プロセッサ１０８は、収縮末期画像２３０および拡張末期画像２４０を識別するために画像処理技術を使用するように構成されてもよい。画像処理技術は、例えば、画像のそれぞれにおける左心室２０２の容積を決定してもよい。プロセッサは、左心室容積が最小の画像を収縮末期画像２３０として、及び左心室容積が最大の画像を拡張末期画像２４０として識別してもよい。あるいは、画像処理技術は、複数の画像にわたって持続する画像要素の動きを識別および追跡して、収縮末期画像２３０および拡張末期画像２４０を識別してもよい。例えば、持続性画像要素の移動方向の反転を使用して収縮末期および拡張末期を識別することができる。収縮末期画像２３０および拡張末期画像２４０は、例えば、適切なフラグでマークされてメモリ１１０内で識別され、それにより、ユーザによって選択および閲覧されてもよい。
【００５２】
左心室２０２の内壁２１８は、収縮末期画像２３０では収縮末期で、拡張末期画像２４０では拡張末期で識別される（ブロック３０４）。左心室２０２は、収縮末期画像２３０の拡張末期および拡張末期画像２４０の収縮末期で輪郭が描かれる（またはマッピングされる）（ブロック３０６）。左心室２０２の輪郭を描くことは、収縮末期画像２３０の内壁２１８の周りの複数の収縮末期輪郭点２３２−ｎ、及び拡張末期画像２４０の内壁２１８の周囲の複数の拡張末期輪郭点２４２−ｎを識別することを含む。第１の連続曲線が収縮末期輪郭点２３２−ｎの点の間に描かれて収縮末期輪郭線２３３が形成され、第２の連続曲線が拡張末期輪郭点２４２−ｎの点の間に描かれて拡張末期輪郭線２４３が形成される。収縮末期輪郭線２３３および拡張末期輪郭線２４３は、直線および／または曲線を含んでもよい。収縮末期輪郭線２３３および拡張末期輪郭線２４３は、例えば、収縮末期画像２３０および拡張末期画像２４０でそれぞれ識別された境界と一致するように描かれてもよい。
【００５３】
収縮末期輪郭点２３２−ｎおよび収縮末期輪郭線２３３は、収縮末期輪郭データセット２３４を形成し、拡張末期輪郭点２４２−ｎおよび拡張末期輪郭線２４３は、拡張末期輪郭データセット２４４を形成する。収縮末期輪郭データセット２３４における各収縮末期輪郭点２３２−ｎは、拡張末期輪郭データセット２４４における拡張末期輪郭点２４２−ｎのうちの対応する点と対になっている。結果として生じる収縮末期および拡張末期輪郭点２３２−ｎ、２４２−ｎの対は、心臓周期中の心臓２００の壁の運動の変化を表す。収縮末期および拡張末期の輪郭点２３２−ｎ、２４２−ｎの対は、心臓周期中の壁運動により収縮末期画像２３０および拡張末期画像２４０の異なる位置ではあるが、左心室２０２の同じ特徴に対応してもよい。本実施形態では、１３個の収縮末期および拡張末期の輪郭点が、収縮末期画像２３０および拡張末期画像２４０において識別される。拡張末期輪郭点および収縮末期輪郭点は、心臓内膜に沿った位置（すなわち、ｎ＝１，２，３...１３）。に応じて１〜１３とラベル付けされる。図５Ａに示すように、収縮末期および拡張末期の等高線データセット２３４、２４４が組み合わされる。
【００５４】
図５Ｂに示されるように、収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎは、内壁の心臓周期運動を表す複数の要素Ｅｉ（ｉは整数）を形成する（ブロック３０８）。本明細書で説明されるように、要素Ｅｉは、左心室２０２の内壁２１８の対応するセクションの心臓周期運動を等級付けするために分析およびスコアリングされる。図示の構成では、収縮末期輪郭点２３２−ｎと拡張末期輪郭点２４２−ｎとの３個の対（すなわち、収縮末期および拡張末期輪郭データセット２３４、２４４のそれぞれからの３つの収縮末期輪郭点２３２−ｎおよび３つの拡張末期輪郭点２４２−ｎ）がそれぞれの要素Ｅｉを定義する。図示の例では、収縮末期画像２３０および拡張末期画像２４０は、二腔心尖部画像である。内壁の心臓周期運動の変化は、６個の要素Ｅ１−６で表される。要素Ｅ１−６のそれぞれは、平面（２次元）多角形の形状をなす。要素Ｅ１−６は、左心室のモデルのそれぞれのセグメントに対応する。具体的には、最初の要素Ｅ１は、基底の下位セグメントに対応し、第２の要素Ｅ２は、中間下位セグメントに対応し、第３の要素Ｅ３は、頂端下部セグメントに対応し、第４の要素Ｅ４は、頂端前部に対応し、第５の要素Ｅ５は、中間前部に対応し、第６の要素Ｅ６は、前部基底部に対応する。心周期変化は、例えば６つ未満の要素または６つを超える要素など、異なる数の要素Ｅｉによって表すことができることが理解されよう。
【００５５】
要素Ｅｉを分析して、心臓２００の対応するセクションの心臓周期運動をスコアリング（すなわち、分類または格付け）するための少なくとも１つの壁運動メトリックを生成する（ブロック３１０）。要素Ｅｉの分析は、本明細書でより詳細に説明される。生成されたメトリックは、事前定義された参照データモデルと比較され、壁の運動をスコアリングする（ブロック３１２）。次に、スコアリングの結果が、例えば画面またはディスプレイに出力される（ブロック３１４）。スコアリングは臨床医によって見直されてもよい。
【００５６】
次に、左心室２０２の輪郭をより詳細に説明する。輪郭付けは心臓超音波検査技師によって実行されてもよく、または適切な画像処理技術を使用して実行されてもよい。造影剤で取得された左心室２０２の超音波画像が図４Ａおよび図４Ｂに示されている。収縮末期画像２３０が図４Ａに示されており、拡張末期画像２４０が図４Ｂに示されている。左心室２０２の心尖部２１２は、左心室２０２の上部として配置することができ、各側面２２２、２２４の基部２１６は、内壁２１８の形状から配置することができる。縦（Ｙ）軸は、心尖部２１２と２つの側面２２２、２２４の基部の間の中点を通る参照線として定義される。ｘ軸は、心尖部と基部２１６の２つの辺の間の中点の間のｙ軸に垂直な線として定義できる。各側面２２２、２２４の中点は、ｘ軸がその側面２２２、２２４の側壁と交差する点として識別できる。中間の収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎは、心尖部２１２と各側面２２２、２２４の中間点との間の領域を細分することによって識別されてもよく、各側面２２２、２２４の中間点と基部との間の領域を細分することによって識別されてもよい。
【００５７】
上述したように、これらの収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎのそれぞれは、ユーザによって識別されてもよい。あるいは、画像処理を使用して、収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎを識別してもよい。画像処理が使用される場合、左心室２０２の輪郭は、最初に左心室２０２内の明るい領域とその周囲の壁を形成する心筋の暗い領域との間の境界として識別される（または、造影剤を使用せずに取得した画像においては逆となる）。そのような境界を識別するための適切なアルゴリズムはよく知られている。境界が識別されると、アルゴリズムは、次に、境界の最高点（最大ｙ値）を心尖部２１２として識別し、境界が下端で方向を変える点を基部２１６として識別するように構成されてもよい。曲率の半径と方向、およびそれが境界の周りでどのように変化するかを分析するためのアルゴリズムを使用して、これらの点、および心尖部２１２の下端の点を識別できる。収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎのそれぞれの座標は、座標系を参照して決定される。心臓超音波検査システム１００によって取得された画像のスケールは既知である。したがって、収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎのそれぞれの座標は、対応する画像の平面内の点の位置を定義する。各対の輪郭点間の距離は、収縮末期と拡張末期の間の心臓２００の対応するセクションによって移動された距離を示す。
【００５８】
次に、図６Ａから図６Ｄを参照して、壁運動メトリックを生成するための要素Ｅｉの分析について説明する。図６Ａに示すように、各要素Ｅｉは、ｎ個の辺を有する多角形の形状をなす。図６Ａに示すように、各要素Ｅｉの要素面積Ａは、座標法を用いて、輪郭データ２３４、２４４の第１および第２のセットから計算される。
【００５９】
【数１】
2021507433.tif 000003
（式１）
ここで、ｎは多角形の辺の数であり、（ｘｉ，ｙｉ）は多角形の頂点である（ｉ＝１，２，．．．，６）。
【００６０】
次に、各要素Ｅｉの計算された面積Ａは、拡張末期の輪郭点の総面積によって表される総面積の割合として正規化される。
【００６１】
図６Ｂに示すように、拡張末期および収縮末期の収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎの各対間のユークリッド距離（ｄ）は、次の式を使用して計算される：
【００６２】
【数２】
2021507433.tif 000004
（式２）
ここで、ｎは次元の数であり、ｐは拡張末期輪郭点の座標であり、ｑは収縮末期輪郭点の座標である。
【００６３】
次に、各要素Ｅｉの平均距離（ｄ）（すなわち、最初の要素Ｅ１では
【００６４】
【数３】
2021507433.tif 000005
２番目の要素Ｅ２では
【００６５】
【数４】
2021507433.tif 000006
など）が計算される。その後、平均距離は、拡張末期輪郭点の全周囲距離の割合として正規化される。
【００６６】
図６Ｃに示すように、各要素Ｅｉの矩形度は、各要素の面積（Ａｉ）とその最小外接矩形の面積（Ｒｉ）の間の比率として計算され：
【００６７】
【数５】
2021507433.tif 000007
（式３）
ここで、Ａｉは各要素Ｅｉの面積であり、Ｒｉは最小外接矩形の面積である。
【００６８】
図６Ｄに示すように、各要素Ｅｉの硬度Ｓｉは、要素の面積（Ａｉ）とその凸包の面積（Ｈｉ）の間の比率として計算され：
【００６９】
【数６】
2021507433.tif 000008
（式４）
ここで、Ａｉは各要素Ｅｉの面積であり、Ｈｉは対応する凸包の面積である。
【００７０】
計算されたメトリクスと参照壁運動スコアの間の相関を評価するために、生の（すなわち、未処理の）二次元心臓超音波検査データを含む参照データセットが分析された。参照データセットは、患者グループの拡張末期画像と収縮末期画像を含む履歴データで構成された。拡張末期画像および収縮末期画像は、本明細書に記載の技法に従って分析され、収縮末期輪郭点２３２−ｎおよび拡張末期輪郭点２４２−ｎが識別された。これにより、左心室２０２の標準モデルの各セグメントに対応する要素Ｅｉが特定された。要素Ｅｉは、以下に示すメトリックを計算するために、本明細書に記載されている手法を使用して分析された：正規化された面積Ａ、正規化された平均距離ｄ、矩形度、および硬度Ｓｉである。メトリックは、各要素Ｅｉの休息およびストレス状態について生成された。要素Ｅｉは、参照壁運動スコアからなる標準スコアリングシステムを使用して、２人の心臓専門医によって互いに干渉することなくスコアリングされた：正常（「１」）、運動低下（「２」）、無動（「３」）、運動異常（「４」）、及びスコアリング不能（「Ｘ」）。２つの参照データセット間で不一致だったスコアは再調査され、一致するに至った。壁運動スコアが「Ｘ」の要素は、参照データセット（ｎ＝２）から削除された。壁運動スコアが「４」（ｎ＝２）である有効な参照データセット内の要素の数が少ないため、これらも分析から削除された。したがって、本実施形態では、参照データセットからの各要素Ｅｉは、正常（「１」）、運動低下（「２」）、無動（「３」）としてスコアリングされた。各患者の安静状態とストレス状態について分析が繰り返された。本明細書で説明されるように、生の参照データセットの分析を通じて計算されたメトリックは、計算されたメトリックが比較され得る参照データモデルを生成するために使用される。参照データモデルは、要素Ｅｉごとに生成される。参照データモデルは、単変量モデルまたは多変量モデルである。参照データモデルは、心臓超音波検査システム１００のメモリ１１０に格納され得る。あるいは、参照データセットは、メモリ１１０に格納され得、参照データモデルは処理ユニット１０６によって生成される。このアプローチは、例えば、特定のメトリックまたはメトリックの組み合わせに関して、動的比較を可能にし得る。
【００７１】
各要素Ｅｉの計算されたメトリックと、安静状態に対して心臓専門医によって割り当てられた壁運動スコア（つまり、正常（「１」）、運動低下（「２」）、無動（「３」））のプロットを図７Ａから７Ｄに示す。メトリックの第２四分位数は、各参照壁運動スコアについて、そのメトリックのセットの中央線を含む箱ひげ図によって表される。各参照壁運動スコアの正規化された要素面積Ａを図７Ａに示す。各参照壁運動スコアの正規化された平均距離ｄを図７Ｂに示す。各参照壁運動スコアの計算された硬度を図７Ｃに示す。各参照壁運動スコアの計算された矩形度を図７Ｄに示す。Ｐ＜０．０５のＷｉｌｋｓ−Ｌａｍｂｄａノンパラメトリック多変量テスト統計が、各壁運動スコアグループで決定される。計算されたメトリックと割り当てられた壁運動スコアとの間に統計的に有意な相関関係が特定される。
【００７２】
各要素Ｅｉの計算されたメトリックと、ストレス状態に対して心臓専門医によって割り当てられた壁運動スコア（つまり、正常（「１」）、運動低下（「２」）、無動（「３」））のプロットを図８Ａから８Ｄに示す。メトリックの第２四分位数は、各参照壁運動スコアについて、そのメトリックのセットの中央線を含む箱ひげ図によって表される。各参照壁運動スコアの正規化された要素面積Ａを図８Ａに示す。各参照壁運動スコアの正規化された平均距離ｄを図８Ｂに示す。各参照壁運動スコアの計算された硬度を図８Ｃに示す。各参照壁運動スコアの計算された矩形度を図８Ｄに示す。Ｐ＜０．０５のＷｉｌｋｓ−Ｌａｍｂｄａノンパラメトリック多変量テスト統計が、各壁運動スコアグループで決定される。計算されたメトリックと割り当てられた壁運動スコアとの間に統計的に有意な相関関係が特定される。
【００７３】
各要素Ｅｉのスコアリングは、計算されたメトリックの1つに基づいて計算できる。処理ユニット１０６は、一変量分布、例えば、一変量正規分布を定義するように構成されてもよい。例として、スコアリングは、計算されたメトリックのうちの１つに対するＺスコア（標準スコア）に対応できる。Ｚスコアは、計算されたメトリックの母平均からの標準偏差の数を標準偏差の単位で示す。処理ユニット１０６は、決定されたＺスコアに対応する各要素Ｅｉにスコアを割り当てるように構成されてもよい。ただし、各要素Ｅｉに対して計算されたスコアの精度は、計算されたメトリックの２つ以上を参照することで改善される場合がある。処理ユニット１０６は、多変量分布、例えば、多変量正規分布を定義するように構成されてもよい。処理ユニット１０６は、二変量分布またはより高次元の分布を定義するように構成されてもよい。処理ユニット１０６は、参照母集団から計算されたメトリックの「距離」を計算するように構成されてもよい。この手法により、より高い次元で分析を実行できる。１つのアプローチは、データＰＣ１の最初の主成分からの距離を使用することである。あるいは、またはさらに、Ｍａｈａｌａｎｏｂｉｓ距離は、処理ユニット１０６によって計算されてもよい。他の統計分析技法も適切である。
【００７４】
平均距離と各要素Ｅｉの面積の間の相関を決定することにより、各要素Ｅｉに対応する壁セクションのスコアを決定できる。本実施形態では、スコアリングは、離散値に依存する従来のスコアリングシステムではなく、連続スケールを含む。図９Ａに、Ｚでスケーリングされた要素の面積と平均距離の間の相関の模式的な表現を示す。左下象限Ｑ３は、潜在的に異常な壁運動があると特定された要素Ｅｉを表す。３つの楕円Ｅ１からＥ３は、データポイントの分布の信頼区間（ＣＩ）を表す：内側の楕円Ｅ１はＣＩ＝６８％を表し、中央の楕円Ｅ２はＣＩ＝９５％を表し、外側の楕円Ｅ３はＣＩ＝９９％を表す。データの重心０が表示され、線ＰＣ１は、データ（すなわち、データＰＣ１の第１主成分）を通る直交回帰線を表す。マーカーＭ１〜Ｍ３は、スコアリングされているデータポイントを表す。ランダムに生成された、０．５６（ｎ＝１，０００）の共分散を有する正規分布データを図９Ｂに示す。各データポイントは、ここで説明する分析手法に従って決定された連続スコアリングに従って色分けされる。
【００７５】
本実施形態による処理ユニット１０６は、各要素Ｅｉをスコアリングするための連続的アプローチを実施する。本明細書で説明するように、主成分モデルは、各要素ＥｉのＺでスケーリングされたメトリックに応じて構築される。これは、収縮末期画像２３０および拡張末期画像２４０から導出される各要素Ｅｉに対して実行される。ここでの説明は、二腔心尖部画像の標準モデルで見えるセグメントに対応する６つの要素Ｅｉに焦点を合わせる。左心室２０２の他のセグメントに対応する追加の要素Ｅｉに関して、例えば三腔心尖部画像および／または四腔心尖部画像の分析によって、同じ技法を実施できることが理解されよう。分析は、安静状態およびストレス状態について取得された収縮末期および拡張末期の画像２３０、２４０に関して独立して行われる。処理ユニット１０６は、安静状態およびストレス状態に関して分析の結果を比較することができる。
【００７６】
２つの計算されたメトリックに基づく二変量分析に基づく要素Ｅｉのスコアリングは、図９Ａおよび９Ｂに視覚化されている。この場合のメトリックは、各要素Ｅｉの正規化された面積と平均距離である。参照データモデルは、参照データセットの分析を通じて特定された要素Ｅｉの正規化された領域と平均距離に応じて生成される。この場合、割り当てられた壁運動スコアが「１」である要素のみが、参照データモデルの生成に含まれていた。所与の要素Ｅｉのスコアを生成するために、処理ユニット１０６は、その要素Ｅｉの対応するメトリックを計算する。処理ユニット１０６は、患者の収縮末期画像２３０および拡張末期画像２４０の分析を通じて識別された要素Ｅｉの正規化された面積および平均距離を計算する。図９Ａおよび９Ｂを参照して本明細書で説明されている実装は、各要素Ｅｉの正規化された面積および平均距離に基づく二変量分析を利用する。メトリックの他の組み合わせを使用して、各要素Ｅｉをスコアリングすることができることが理解されよう。たとえば、二変量分析では、平均距離と硬度メトリック、または、正規化された面積と矩形度のメトリックを組み合わせることができる。
【００７７】
処理ユニット１０６は、多変量分析を実行するように構成されてもよい。処理ユニット１０６は、本明細書に記載される各メトリック、すなわち、正規化された面積Ａ、正規化された平均距離ｄ、矩形度、および硬度Ｓｉを組み合わせるように構成されてもよい。各要素Ｅｉのスコアは、４つの計算されたメトリックの多変量解析に基づいて計算されてもよい。第１主成分ＰＣ１および割り当てられた壁運動スコア（すなわち、正常な「１」、運動低下「２」、無動「３」）に応じて計算されたスコアのプロットを図１０、１０Ｂ、１１Ａおよび１１Ｂに示す。二腔心尖部画像における安静状態の多変量解析の第１主成分ＰＣ１に応じて計算されたスコアのプロットを図１０Ａに示す。四腔心尖部画像における安静状態の多変量解析の第１主成分ＰＣ１に応じて計算されたスコアのプロットを図１０Ｂに示す。二腔心尖部画像におけるストレス状態の多変量解析の第１主成分ＰＣ１に応じて計算されたスコアのプロットを図１１Ａに示す。四腔心尖部画像におけるストレス状態の多変量解析の第１主成分ＰＣ１に応じて計算されたスコアのプロットを図１１Ｂに示す。
【００７８】
各要素Ｅｉの計算されたメトリックは、対応する要素Ｅｉの参照データモデルと比較される。各要素Ｅｉのスコアは、この比較に基づいて計算される。スコアは、新しいデータの最初の主成分の値、すなわち、計算されたメトリックがデータの重心からどれだけ離れているか、したがって、それらが参照データとどれほど異なっているかを表す。スコアは連続的なスケールで計算できる。この例は、図１２Ａと１２Ｂに示されており、各要素には、ここで説明する連続スコアリングスケールに従って陰影が付けられている。図１２Ａに示すスコア付けされた画像は、安静状態の二腔データを表しており、図１２Ｂに示すスコア付けされた画像は、ストレス状態の二腔データを表している。
【００７９】
処理ユニット１０６は、各要素Ｅｉをスコアリングするための異なるメトリックを計算するように構成されてもよい。これらのメトリックは、単変量および多変量解析のために本明細書に記載されているメトリックに加えて、またはその代わりに使用されてもよい。処理ユニット１０６は、例えば、以下のメトリックのうちの１つ以上を計算してもよい：剪断；ひずみ;要素内の距離の変動係数;アスペクト比。処理ユニット１０６はまた、本明細書で説明される点の対間の距離以外の距離メトリックを計算してもよい。例えば、処理ユニット１０６は、同じ要素内の正反対の点の間の距離、または、異なる要素の対応する点間の距離（すなわち、ＥｉとＥｉ＋ｊの間）を計算してもよい。
【００８０】
処理ユニット１０６は、メトリックを計算するための１つの画像における要素Ｅｉの分析に特に重点を置いて説明されてきた。処理ユニット１０６は複数の画像を分析し得ることが理解されるであろう。異なる画像に同じ要素Ｅｉが含まれてもよい。処理ユニット１０６は、異なる画像の分析に応じて、特定の要素Ｅｉに対して生成されたスコアを比較するように構成されてもよい。スコア間に不一致が検出された場合、これは画質などの潜在的な問題として報告されることが出来る。これにより、さまざまな画像の品質をチェックできる。同様に、特にきめの細かい要素Ｅｉでは、互いに近接して配置された要素Ｅｉについて計算されたスコアは通常、関連している。予想される関係が特定された場合、これは潜在的な画質の問題としてフラグを立てることができる。これらの技法は、より正確である可能性が高い要素Ｅｉの１つまたは複数の画像の識別を可能にし得る。
【００８１】
処理ユニット１０６は、１つまたは複数の画像内に見ることが出来る要素Ｅｉを分析して、見ることが出来ない要素Ｅｉの挙動を推測することができる。例えば、スコアは、見えない要素Ｅｉに隣接してまたは近位に配置された少なくとも１つの要素Ｅｉについての計算されたスコアに応じて、見えない要素Ｅｉについて推定されてもよい。有効な表示に基づいて左心室２０２の完全なモデルを構築するために、さらなる分析が実行されてもよい。結果として得られるモデルにより、１つまたは複数の近くの要素Ｅｉからスコアリングを推測することができる。
【００８２】
ここで説明する統計分析は、ｇｇｐｌｏｔ２、ｄｐｌｙｒ、ｇｇｂｅｅｓｗａｒｍ、Ｍｏｍｏｃｓ、ｐａｔｈｍａｐｐｉｎｇ、およびｇｇｐｕｂｒパッケージを使用して、Ｒ統計的計算環境（ｖ３．４．１）内で実行された。各比較におけるグループ間の観測数の不均衡のため、多変量ノンパラメトリック仮説検定を採用し、ｎｐｍｖパッケージを使用してグループ平均を比較した。タイプＩのエラー率（α）０．０５がすべての比較に使用された。
【００８３】
処理ユニット１０６は、左心室１０２の標準モデルのセグメントに対応する要素Ｅｉのメトリックを計算するものとして本明細書で説明されている。本明細書で説明される技法は、要素Ｅｉがセグメントに対応することを必要としないことが理解されよう。たとえば、要素Ｅｉは標準モデルのセグメントよりも小さい場合がある。要素Ｅｉは、標準モデルのサブセグメントに対応し得る。要素Ｅｉのサイズを小さくすることにより、スコアリングは、心臓の周期運動の異常の位置をより正確に示すことができる。例えば、スコアリングは、例えば、セグメントの境界の近くの位置または中央の場所を強調して、標準モデルのセグメントの１つ内の異常な機能の場所を示すことが想定される。
【００８４】
本明細書で説明されるように、参照データモデルは、生の２次元心臓超音波検査データを含む参照データセットを分析することによって生成される。上記の実施形態では、参照データセットは、患者のグループの拡張末期画像および収縮末期画像を含む履歴データを含む。さらなる発展において、参照データモデルを生成するために使用される参照データセットは、繰り返し更新されてもよい。たとえば、新しい心臓超音波検査の分析を参照データセットに組み込んで、有効なデータ母集団を増やすことができる。したがって、追加のデータが有効になると、参照データモデルは引き続き改良される。参照データモデルの反復的な発展が、病理学的変化と患者の進展を可能にし得る。
【００８５】
ここで説明する分析は、要素Ｅｉごとに独立して実行される。ここで、分析は、複数の要素Ｅｉ間の関係を考慮するように修正されてもよいことが理解されるであろう。例えば、分析は、例えば、左心室２０２の反対側に、互いに隣接して、または互いに反対に配置された第１および第２の要素Ｅｉの運動を同時にスコア付けすることができる。
【００８６】
本明細書で説明される参照データモデルは、参照データセットの一部またはすべてに関して有効な結果データに応じて変更することもできる。本明細書では、「結果データ」という用語は、診断情報を指すために使用される。結果データは、参照データセット内の対応する記録またはデータセットに関連付けられる。診断情報は、例えば、患者の血管造影データおよび／または心臓事象に関連してもよい。結果データは、経過時間中に、患者が心臓状態、例えば冠状動脈疾患の有無の陽性または陰性の診断を下したかどうかを示してもよい。結果データは、例えば、心臓超音波検査データの取得の１年後、２年後、または３年後に生成されてもよい。参照データセット内のデータの重み付けは、結果データに応じて調整されてもよい。例えば、結果データが有効である参照データセット内のデータに適用される重み付けは、その統計的有意性を変更するために増加または減少されてもよい。重み付けは、心臓超音波画像の取得とその後の診断イベントとの間の経過時間に応じて調整することができる。変形例では、参照データモデルは、結果データが有効なデータにのみ依存して生成されてもよい。参照データモデルは、結果データが冠動脈疾患などの特定の状態の有無を示すデータにのみ依存して生成されてもよい。結果データを使用して、参照データセットをフィルター処理し、異なる参照データモデルを生成してもよい。
【００８７】
処理ユニット１０６はまた、診断出力を生成する診断機能を提供するように構成されてもよい。診断システムは、出願人の国際特許出願ＰＣＴ／ＧＢ２０１７／０５１７２０に開示されており、その内容は参照により全体として本明細書に組み込まれる。診断機能は、本明細書に記載されている結果データを利用できることが認識されている。診断機能は、壁運動スコアに対して生成された1つ以上のメトリックに依存する場合もある。例として、各要素Ｅｉの矩形度は、診断モデルにおける特徴として使用されてもよい。診断モデルを生成する際の結果データの使用は、さまざまな疾患の割合と特性を経時的に、および/またはさまざまな医療現場で考慮するのに役立つことがある。たとえば、異なる現場では、ポジティブ（「疾患」）とネガティブ（「正常」）の結果の比率が異なる場合がある。診断モデルを生成する際に結果データを利用することにより、これらのタイプのバリエーションを考慮に入れることができる。ストレスエコーテストの結果（テスト中/直後に心臓専門医が決定）は、常に正確であるとは限らない。出願人が行った１年の結果の正確さの分析は、複数のデータセットにわたるストレスエコー結果の平均７．２％の不正確さを示している。一定期間にわたる結果データを参照することにより、診断モデルの精度が向上し、個人が冠動脈疾患などの疾患を発症するかどうかをモードで正確に予測できるようになる。結果データの使用は、単独で特許を受けることができると考えられる。この強化された診断機能は、前の実施形態の発展として説明される。同様の参照番号は、同様のコンポーネントに使用される。
【００８８】
本明細書で説明するように、各収縮末期輪郭点２３２−ｎは、拡張末期輪郭データセット２４４の拡張末期輪郭点２４２−ｎの対応する1つと対になっている。収縮末期と拡張末期輪郭点２３２−ｎと２４２−ｎの結果の対は、心臓周期中の心臓２００の壁の運動の変化を表す。収縮末期および拡張末期の輪郭点２３２−ｎ、２４２−ｎが識別されると、デカルト座標系におけるそれらのｘおよびｙ座標は、例えば、収縮末期画像上の点の座標を含む収縮末期座標セットと拡張末期画像上の点の座標を含む拡張末期座標セットとしてメモリ１１０に格納されてもよい。プロセッサは、２つの座標セットから、収縮末期と拡張末期との間の左心室２０２の形状の変形を計算するように構成されてもよい。
【００８９】
処理ユニット１０６は、収縮末期と拡張末期の間の左心室２０２の動きを定量化する様々なパラメータの値を計算するように構成されている。計算には、ｘ軸とｙ軸の両方に沿った位置の変化（拡張末期―収縮末期）を計算することにより、各点がｘ方向とｙ方向のそれぞれにどれだけ移動したかを計算することが含まれてもよい。これにより、収縮末期と拡張末期の輪郭点２３２−ｎ、２４２−ｎの対応する対のそれぞれについて、ｘ軸の移動Δｘのセットとｙ軸の移動Δｙのセットが得られる。これらの各値は、方向を示さない単純な距離であってよい。次に、ｘ軸（ΔＸ）とｙ軸（ΔＹ）の両方のすべての点の変化量の平均値を別々に計算して、左心室２０２全体に対して、平均のΔｘ値またはｘ方向の移動ΔＸ、および平均のΔｙ値またはｙ方向の移動ΔＹを提供してもよい。個々の移動値のそれぞれが純粋な距離であり、正または負のｘまたはｙ方向にあるかどうかの表示がない場合、これらの平均値は総移動量を示すが、ＬＶ壁の異なる部分が同じ方向に移動しているか、反対の方向に移動しているかは示されていない。
【００９０】
収縮末期と拡張末期の輪郭点２３２−ｎ、２４２−ｎの各対に対して計算できる別のパラメータは、ｘおよびｙ方向の移動ΔｘおよびΔｙの平均であって、各点の平均値Δｘｙ＝（Δｘ＋Δｙ）／２である。次に、すべての点についてのΔｘｙのすべての値の平均を計算して、心室全体の値ΔＸＹを計算できる。この計算は、せん断ひずみの計算に似ているため、ここではせん断変形と呼ぶ。所与の移動距離の場合、このパラメータは、ｘ軸とｙ軸の両方に対して４５度の移動で最大になり、軸の１つに沿った移動で最小になることが理解されよう。
【００９１】
計算できる別のパラメータは、ｘとｙのひずみ成分から計算できる主変形に似ているため、ここでは主変形と呼び、
主変形＝Ｃ１（ΔＸ＋ΔＹ−√（ΔＸ＋ΔＹ）＾２＋Ｃ２ΔＸＹ＾２）
で与えられる。ここで、Ｃ１とＣ２は定数である。定数Ｃ１は、例えば、１／２であってよく、定数Ｃ２は、４であってよい。これらの値は、以下に記載される例において使用される。
【００９２】
この変換は、せん断変形と密接に関連しているため、そのパラメータと同様に変化する傾向があるが、心臓の収縮を示す負の値を有する。ただし、以下のテスト結果で示されているように、主な変形値は、一部のケース、特に冠動脈疾患（ＣＡＤ）で、より信頼できる診断を提供できる。
【００９３】
これらのパラメータのそれぞれは、単一の冠状動脈周期における収縮末期と拡張末期の間の変化に関連していることが理解されよう。ただし、ストレス心臓超音波検査（または他のイメージング方法で実行される対応するテスト）では、安静時の心臓に対し各パラメータに１つの値とストレス時の心臓に対し１つの値がある。それらの値を比較すること、たとえば、それらの値の違いを判断することで、診断に使用できる心臓機能についてのさらなる情報を得ることができる。
【００９４】
xとyの移動、せん断と主変形の値が計算されると、プロセッサは、これらをメモリ１１０に格納されている参照値と比較して、1つまたは複数の特定の心臓の状態の診断を行い、診断の出力を生成するように構成される。出力は、陽性または陰性の診断を示す単純なバイナリ出力であってもよい。処理ユニット１０６は、ディスプレイ１１２上に出力を表示するように構成されてもよい。あるいは、またはそれに加えて、例えば出力データを追加することにより、出力が基づいた画像に関連するデータとして、診断を示す出力を画像が保存されているファイルへ格納するように構成されてもよい。
【００９５】
基準値は、例えば、処理ユニット１０６上で実行することができ、従来の方法によって決定された関連診断を伴うストレスエコー画像のデータベースを使用する学習アルゴリズムによって決定することができ、これは、メモリ１１０に格納されてもよい。具体的には、データベースは、多数の画像のセットを含んでもよく、各セットは、安静状態およびストレス状態の両方についての収縮末期画像および拡張末期画像を、各画像セットについて、冠動脈疾患（ＣＡＤ）の陽性または陰性の診断などといった関連する診断とともに含む。学習アルゴリズムは、画像を分析して上記の様々なパラメータの値を計算し、次いで診断と様々なパラメータのそれぞれの値との間の相関を決定するように構成されてもよい。
【００９６】
分析は７０人の被験者からのサンプル画像で行われた。生成されたすべての結果は、心尖部四腔画像によるものである。最初に、ＤＳＥの結果から決定された正と負の結果について値を比較した。次に、ＤＳＥ結果で確認された偽陽性を修正するためにＤＳＥ結果を修正して、比較を繰り返した。
【００９７】
表１に、心尖部四腔画像でのＤＳＥの結果（１＝Ｐｏｓ、２＝Ｎｅｇ）の安静状態およびストレス状態の主変形（ｍｍ）、せん断変形値（ｍｍ）、平均ΔＸ（ｍｍ）の値を示す。
【００９８】
【表１】
2021507433.tif 000009
【００９９】
表２に、調整後のＤＳＥの結果（１＝Ｐｏｓ、２＝Ｎｅｇ）の安静状態およびストレス状態の主変形値（ｍｍ）、せん断変形（ｍｍ）およびＸ変形（ｍｍ）の平均を示す。
【０１００】
【表２】
2021507433.tif 000010
【０１０１】
表３に、変数対調整後のＤＳＥの独立標本Ｔ検定を示す。
【０１０２】
【表３】
2021507433.tif 000011
【０１０３】
サンプルデータから取得したさまざまなパラメータの値から、機械学習を使用して、調整されたドブタミン負荷心エコー（ＤＳＥ）の結果の指標として各パラメータの精度を決定してもよい。上記のデータを用い、１０分割交差検証法を備えたＪ４８剪定決定木を使用してデータが分類された。診断結果の指標としての各パラメータの精度は、以下の表に要約されており、ここでは、次の略語が使用されている：ＴＰ＝真陽性、ＦＰ＝偽陽性、ＦＮ＝偽陰性、ＴＮ＝真陰性、ＰＰＶ＝陽性的中率、およびＮＰＶ＝陰性的中率。
【０１０４】
表４、コンサルタントの解釈の精度。
【０１０５】
【表４】
2021507433.tif 000012
【０１０６】
表５、調整されたＤＳＥの結果に対するストレス状態の主変形の精度。
【０１０７】
【表５】
2021507433.tif 000013
【０１０８】
表６、調整されたＤＳＥの結果に対する安静状態の主変形の精度。
【０１０９】
【表６】
2021507433.tif 000014
【０１１０】
表７、調整されたＤＳＥの結果に対するストレス状態のせん断変形の精度。
【０１１１】
【表７】
2021507433.tif 000015
【０１１２】
表８、調整されたＤＳＥの結果に対する安静状態のせん断変形の精度。
【０１１３】
【表８】
2021507433.tif 000016
【０１１４】
次に、すべての変数から、機械学習を使用して、データから正確な診断を提供するために、図１３に示す決定木が導き出された。決定木は一連の決定点を定義し、各決定点はパラメータの参照値またはしきい値を定義する。決定木は、次のように動作する単純なアルゴリズムの概要を示す。最初に、上記のような左心室２０２の主変形が心臓のストレス状態について決定される。変形が−５．９５ｍｍ未満（つまり、大きさが５．９５ｍｍより大きい負の値）の場合、診断は陰性である。値が−５．９５ｍｍより大きい場合（つまり、大きさが５．９５ｍｍより大きい負の値）、安静状態とストレス状態の間の主変形の差は１２．２７８０５３ｍｍより大きい場合、診断は負になるが、その距離より短い場合、診断は陽性となる。決定木の構造、および決定木内の各決定点における参照またはしきい値は、行われる診断に依存することが理解されよう。
【０１１５】
上述の処理ユニット１０６は、冠状動脈疾患を診断するための固定（静的）診断モデルを実装する。図１３に示すように、決定木は一連の決定点を定義し、各決定点はパラメータの参照値またはしきい値を定義する。処理ユニット１０６は、動的診断モデルを実装してもよい。決定木の各決定点における参照値または閾値は、例えば、新しい心臓超音波検査データおよび／または結果データを反映するように動的に修正されてもよい。少なくとも特定の実施形態では、これは、改善された診断機能を提供することができる。
【０１１６】
新しい心臓超音波検査データは、既存の参照データのセットに段階的に追加されてもよい。新しいデータはデータ母集団を拡大するために使用され、診断モデルを徐々に変更する可能性がある。決定木で使用される参照値またはしきい値は、有効な心臓超音波検査データを反映するように更新されてもよい。参照データの反復的な発展により、診断モデルを時間の経過とともに変化させることができる。本明細書で説明される決定木は、監視された機械学習モデルなどの他の分析ツールで置き換えられてもよいことが理解されよう。
【０１１７】
結果データは、例えば血管造影データおよび／または心臓事象に関する、各患者の診断情報を含む。本実施形態における結果データは、参照データセットで使用される収縮末期画像および拡張末期画像の取得後の経過時間中に冠動脈疾患の有無が検出されたかどうかを示す。結果データは、例えば、心臓超音波検査データの取得の１年、２年、３年、またはそれより長期間後に生成されてもよい。本実施形態における結果データは、心エコー検査データを取得してから１年後に生成される。結果データは、経過時間中に発生した血管造影データおよび心臓事象を考慮することによって編集される。結果データは時間とともに進化し続けることが理解されよう。したがって、結果データは継続的に、たとえば年次ベースで、または分類が変更されたときに更新され得る。結果データを更新することにより、参照データに基づいて生成された診断ツールと診断モデルを動的に調整して、病理学的変化と患者の進展を表すことができる。
【０１１８】
動的診断モデルを実装するために、教師付き機械学習アルゴリズムを使用して分類モデルが構築される。結果データは、機械学習アルゴリズムによってアクセスされる参照データにラベルを付けるために使用される。機械学習アルゴリズムは、ラベルを使用して異なる分類を区別する。本実施形態では、分類は、冠動脈疾患の有無に対応する。あるいは、またはさらに、分類は、例えば、特定された狭窄レベルまたは割合に依存して、特定の状態を等級付けし得る。分類は他の条件を区別し得ることが理解されよう。診断モデルを生成するとき、機械学習アルゴリズムは、結果データから導出されたラベルに応じて参照データの相対的重み付けを調整することができる。少なくとも特定の実施形態では、結果データに応じて参照データを更新することにより、ストレスエコー図に基づいて診断精度を向上させることができる。
【０１１９】
診断モデルを構築するために、輪郭データから一連の特徴が計算される。特徴は、セグメントごとに（たとえば、本明細書で説明する要素Ｅｉの１つまたは複数を分析することによって）、任意選択で左心室２０２全体に関して計算される。有効な特徴セットを分析して、最も関連性のある特徴を識別する。これにより、最も適切な特徴を識別して、診断モデルを構築できる。（複数の決定木で構成される）ランダムフォレストの場合、識別された特徴が決定ノードを形成する。最も関連性の高い特徴は、地域によって異なる場合や、疾患の進展に伴って変化する場合があり、モデルで使用するために識別された特徴は変更される場合がある。特徴が同じであっても、しきい値と重み付けは変更される場合がある。本明細書の表１および２に示されるように、ピークストレスでの駆出率として、最重要の特徴は変化しない。ただし、次に重要な特徴は、さまざまな条件によって異なる。最初のデータセット（表１）では、収縮末期と拡張末期の間の体積変化が次に最も関連している。ただし、結合されたデータセットでは、二腔画像で安静時の特定のセグメントの領域が次に最も関連している。モデルをトレーニングするには、参照データにラベルを付ける必要がある。潜在的な不正確さを考慮して、ラベルとして（心臓専門医によって決定された）ストレスエコーの結果を使用しても、必ずしも正確なモデルが導出されるとは限らない。参照データの取得後一定期間（たとえば、心臓超音波画像の取得後1年）に収集された結果データを使用すると、これらの欠陥の少なくとも一部を克服または改善できる。
【０１２０】
結果データは、さまざまな期間にわたって収集できる。結果データは、少なくとも特定の実施形態では、冠状動脈疾患の影響をどれだけ前もって特定できるかを示す指標を提供することができる。さらに、複数のクラスのラベルを使用して、さまざまな疾患の重症度を予測できる。疾患が進展したり人口が変化する可能性があるため、より多くの結果データが蓄積されると、診断モデルが更新され、分類が可能な限り正確に保たれるようになる。これは、新しい結果データが受信されるたびにモデル全体を再トレーニングすることで実行できる。実際には、これには時間がかかる場合がある。別の方法として、診断モデルを継続的に更新するために、プロセッサによって増分機械学習技術を実装することができる。
【０１２１】
次に、分類モデルの実装について、第１の参照データセットおよび第２の参照データセットを参照して説明する。最初のデータセットは、１２４のストレス心エコー図の最初のセット（２０１１年５月から２０１３年８月の間にオックスフォードで収集）で構成される。２番目のデータセットは、個別の研究（２０１５年３月から２０１６年８月の間にテムズバレー全体の６つの異なる病院で収集）から得られた３３９のストレス心エコー図のセットで構成されている。結果データは、ストレス心エコー図を取得してから１年後に編集される。結果データは、バイナリの結果値を生成する。特に、経過した１年の間に次のイベントのいずれかが識別された場合、結果は肯定的であると見なされる：
（ｉ）心臓事象（例、心筋梗塞）；
（ｉｉ）７０％を超える狭窄を示した血管造影。
経過した1年の間に前述の事象（ｉ）または（ｉｉ）のいずれも発生しなかった場合、結果は否定的と見なされる。最初のデータセットでは、１０の肯定的な結果が識別され、２番目のデータセットでは、１３の肯定的な結果が識別された。
【０１２２】
結果を予測するための最も関連性の高い特徴を評価するＲ統計的計算環境のＢｏｒｕｔａパッケージ。Ｂｏｒｕｔａパッケージは、属性の重要性を可能なものとランダムに比較することにより、特徴の選択を実行する。５００本の木を持つランダムフォレストを含む標準実装が実装された。表９は、最初のデータセットの最も重要な特徴とその平均重要度スコアを示す。２番目のデータセットが最初のデータセットに追加された。表１０は、１番目と２番目のデータセットを組み合わせた場合の最も重要な特徴とその平均重要度スコアを示す。処理に使用できるデータが増えると、最も関連する特徴が変化する。これは、参照データが追加されると分類モデルが変わる可能性があることを示している。これらの変更は、追加の参照データが後日取得された場合、および/またはより広範囲の地理的領域にわたって取得された場合に、より顕著になると考えられる。本明細書ではランダムフォレストモデルの使用について説明したが、別のモデル、または実際にはモデルの組み合わせを使用できることは理解されよう。
【０１２３】
表９、最初のデータセットを使用する最も関連性の高い特徴。
【０１２４】
【表９】
2021507433.tif 000017
【０１２５】
表１０、１番目と２番目のデータセットの結合データを使用した最も関連性の高い特徴。
【０１２６】
【表１０】
2021507433.tif 000018
【０１２７】
新しい参照データを組み込むことができる継続的な学習戦略の実装により、より堅牢で正確な診断モデルを実現できる。新しい参照データを段階的に組み込むことにより、新しいデータが有効になるたびに、（特に参照データセットのサイズが大きくなるにつれて時間のかかるプロセスであることを示す）モデル全体を再トレーニングする必要性を減らすか回避することができる。さらに、診断モデルは、時間とともに変化する疾患の特徴に適応できる。人々の環境やライフスタイルの変化により、最も関連性の高いバイオマーカーが時間とともに変化する可能性があり、モデルはこれらに対応する必要があるため、これは特に重要である。動的診断モデルは、心血管疾患の様相と特性の変化に適応できるため、堅牢で正確な予測モデルを提供できる。
【０１２８】
ここで説明する動的診断モデルは、1年間にわたって取得された結果データを利用する。結果データは、異なる期間にわたって蓄積されてもよいことが理解されよう。より長い期間にわたる結果データを組み合わせることにより、より長い期間にわたる動的診断モデルの予測力を評価することができる。
【０１２９】
本出願は、心血管疾患を参照して説明されてきた。しかしながら、本明細書に記載される方法および装置は他の用途を有し得ることが理解されるであろう。例えば、環境が変化して腫瘍のサイズ、外観、または石灰化が変化した場合に、変化する腫瘍のイメージングバイオマーカーに適応する診断ツールを開発できる。さらに、この技法は心臓超音波検査以外の画像システムにも適用し得る。
【０１３０】
添付の特許請求の範囲から逸脱することなく、本明細書に記載された実施形態に様々な変更を加えることができることが理解されよう。
【図１】
2021507433.tif 000019
【図２】
2021507433.tif 000020
【図３】
2021507433.tif 000021
【図４Ａ】
2021507433.tif 000022
【図４Ｂ】
2021507433.tif 000023
【図５Ａ】
2021507433.tif 000024
【図５Ｂ】
2021507433.tif 000025
【図６Ａ】
2021507433.tif 000026
【図６Ｂ】
2021507433.tif 000027
【図６Ｃ】
2021507433.tif 000028
【図６Ｄ】
2021507433.tif 000029
【図７Ａ】
2021507433.tif 000030
【図７Ｂ】
2021507433.tif 000031
【図７Ｃ】
2021507433.tif 000032
【図７Ｄ】
2021507433.tif 000033
【図８Ａ】
2021507433.tif 000034
【図８Ｂ】
2021507433.tif 000035
【図８Ｃ】
2021507433.tif 000036
【図８Ｄ】
2021507433.tif 000037
【図９Ａ】
2021507433.tif 000038
【図９Ｂ】
2021507433.tif 000039
【図１０Ａ】
2021507433.tif 000040
【図１０Ｂ】
2021507433.tif 000041
【図１１Ａ】
2021507433.tif 000042
【図１１Ｂ】
2021507433.tif 000043
【図１２Ａ】
2021507433.tif 000044
【図１２Ｂ】
2021507433.tif 000045
【図１３】
2021507433.tif 000046
【国際調査報告】
2021507433.tif 000047
2021507433.tif 000048
2021507433.tif 000049