(57)【要約】  がん悪性度マッピングを実行する超音波システムは、超音波撮像データを取得する超音波撮像装置１０を含む。電子データ処理装置３０は、前記超音波撮像データから超音波画像３４を生成し、がん悪性度マップ４２を生成するようにプログラムされ、がん悪性度マップ４２は、（ｉ）前記超音波撮像データから、前記がん悪性度マップのマップ画素を表すローカルフィーチャのセットを抽出し、（ｉｉ）前記がん悪性度マップの前記マップ画素に対するがん悪性度を生成するのにがん格付け分類器４６を使用して前記ローカルフィーチャのセットを分類することにより生成される。表示コンポーネント２０は、例えば色符号化がん悪性度マップオーバレイとして前記超音波画像に重ねられた、前記がん悪性度マップを表示する。前記がん格付け分類器は、生検場所における超音波撮像データから抽出され、組織病理学的がん悪性度でラベル付けされたローカルフィーチャのセットを有するトレーニングデータセット６４から学習される。
2018516135.tif 000002
【特許請求の範囲】
【請求項１】
超音波撮像データを取得する超音波撮像装置と、
がん悪性度マップを生成するようにプログラムされた電子データ処理装置であって、前記がん悪性度マップが、（ｉ）前記超音波撮像データから、前記がん悪性度マップのマップ画素を表すローカルフィーチャのセットを抽出し、（ｉｉ）前記がん悪性度マップの前記マップ画素に対するがん悪性度を生成するのにがん格付け分類器を使用して前記ローカルフィーチャのセットを分類することにより生成される、当該電子データ処理装置と、
前記がん悪性度マップを表示する表示コンポーネントと、
を有する超音波システム。
【請求項２】
前記電子データ処理装置が、ＲＦ時系列超音波撮像データからマップ画素を表す前記ローカルフィーチャのセットを抽出するようにプログラムされる、請求項１に記載の超音波システム。
【請求項３】
前記超音波撮像装置は、組織振動を誘起するように前記超音波装置により低周波数の超音波パルスが印加されるエラストグラフィ撮像データを含む超音波撮像データを取得し、
前記電子データ処理装置が、前記エラストグラフィ撮像データからマップ画素を表す前記ローカルフィーチャのセットを抽出するようにプログラムされる、
請求項１に記載の超音波システム。
【請求項４】
前記電子データ処理装置が、前記超音波撮像データから超音波画像を生成し、前記がん悪性度マップから、前記超音波画像とアラインされるがん悪性度マップオーバレイを生成するようにプログラムされ、
前記表示コンポーネントが、前記超音波画像及び前記がん悪性度マップを結合する融合画像を表示する、
請求項１乃至３のいずれか一項に記載の超音波システム。
【請求項５】
前記電子データ処理装置が、ＲＦ時系列超音波撮像データを有する超音波撮像データから輝度（ｂモード）画像として前記超音波画像を生成するようにプログラムされる、請求項４に記載の超音波システム。
【請求項６】
前記電子データ処理装置は、前記がん悪性度マップのがん悪性度が色符号化により表される色符号化がん悪性度マップオーバレイと重ねられた前記超音波画像として前記融合画像を生成するようにプログラムされる、請求項４乃至５のいずれか一項に記載の超音波システム。
【請求項７】
前記超音波システムが、超音波撮像データを連続的に取得し、前記連続的に取得された超音波撮像データを使用して前記超音波画像、前記がん悪性度マップ、及び前記融合画像をリアルタイムで更新する、請求項４乃至６のいずれか一項に記載の超音波システム。
【請求項８】
前記がん悪性度マップの各マップ画素が、前記取得された超音波撮像データから生成された超音波画像の画素の隣接するｎ×ｎのアレイからなり、ｎ≧１である、請求項１乃至７のいずれか一項に記載の超音波システム。
【請求項９】
前記電子データ処理装置が、（１）テクスチャフィーチャ、（２）ウェーブレットベースのフィーチャ、及び（３）スペクトルフィーチャの１以上を含む前記がん悪性度マップのマップ画素を表す前記ローカルフィーチャのセットを抽出するようにプログラムされる、請求項１乃至８のいずれか一項に記載の超音波システム。
【請求項１０】
前記電子データ処理装置が、前記超音波撮像データを取得するように前記超音波撮像装置を制御するようにプログラムされた前記超音波撮像装置のマイクロプロセッサ又はマイクロコントローラを有する、請求項１乃至９のいずれか一項に記載の超音波システム。
【請求項１１】
前記超音波撮像装置と接続された直腸超音波プローブであって、前記超音波撮像装置が、前記直腸超音波プローブを使用して前立腺器官の超音波撮像データを取得し、前記電子データ処理装置が、前立腺がん悪性度マップを生成するようにプログラムされ、前記前立腺がん悪性度マップが、（ｉ）前記超音波撮像データから、前立腺がん悪性度マップのマップ画素を表すローカルフィーチャのセットを抽出し、（ｉｉ）前記前立腺がん悪性度マップの前記マップ画素に対する前立腺がん悪性度を生成するのに前立腺がん格付け分類器を使用して前記ローカルフィーチャのセットを分類することにより生成され、前記表示コンポーネントが、前記前立腺がん悪性度マップを表示する、請求項１乃至１０のいずれか一項に記載の超音波システム。
【請求項１２】
前記超音波システムが、前記直腸超音波プローブと接続され、前立腺組織生検標本を収集する直腸生検ツールを有し、
前記電子データ処理装置が、前記超音波撮像データから前立腺超音波画像を生成するようにプログラムされ、前記表示コンポーネントが、前記前立腺超音波画像及び前記前立腺がん悪性度マップを結合する融合画像を表示する、
請求項１１に記載の超音波システム。
【請求項１３】
生検場所における超音波撮像データから抽出され、組織病理学的がん悪性度でラベル付けされたローカルフィーチャのトレーニングセットを有するラベル付けされたトレーニングデータセットに対する機械学習により前記がん格付け分類器を生成するようにプログラムされた電子データ処理装置、
を有する、請求項１乃至１２のいずれか一項に記載の超音波システム。
【請求項１４】
超音波撮像データを取得するステップと、
前記超音波撮像データから超音波画像を生成するステップと、
前記超音波撮像データから抽出されたローカルフィーチャのセットにがん格付け分類器を使用することにより前記超音波撮像データからがん悪性度マップを生成するステップと、
（ｉ）前記がん悪性度マップ及び（ｉｉ）前記超音波画像及び前記がん悪性度マップを結合する融合画像の少なくとも１つを表示するステップと、
を有する超音波方法。
【請求項１５】
前記超音波撮像データが、ＲＦ時系列超音波撮像データを含み、
前記超音波画像が、前記ＲＦ時系列超音波撮像データから生成された輝度モード（ｂモード）画像を有し、
前記がん悪性度マップが、前記ＲＦ時系列超音波撮像データから生成される、
請求項１４に記載の超音波方法。
【請求項１６】
前記表示するステップが、前記がん悪性度マップの色符号化オーバレイ表現と重ねられた前記超音波画像を有する融合画像を表示するステップを有する、請求項１４乃至１５のいずれか一項に記載の超音波方法。
【請求項１７】
前記表示されるがん悪性度マップ又は融合画像をリアルタイムで更新するように、前記取得するステップ、前記超音波画像を生成するステップ、前記がん悪性度マップを生成するステップ、及び前記表示するステップを反復的に繰り返すステップを有する、請求項１４乃至１６のいずれか一項に記載の超音波方法。
【請求項１８】
前記がん悪性度マップにより示される高い悪性度のがん組織の領域に生検ツールをガイドするように前記表示されるがん悪性度マップ又は融合画像を使用して疑われるがん亜組織の超音波ガイド生検を実行するステップ、
を有する、請求項１４乃至１６のいずれか一項に記載の超音波方法。
【請求項１９】
前記がん悪性度マップにより示される高い悪性度のがん組織の領域を標的にしてがん治療を実行するステップ、
を有する、請求項１４乃至１６のいずれか一項に記載の超音波方法。
【請求項２０】
生検場所における超音波撮像データから抽出され、組織病理学的がん悪性度でラベル付けされたローカルフィーチャのトレーニングセットを有するラベル付けされたトレーニングデータセットに対して前記がん格付け分類器をトレーニングするステップ、
を有する、請求項１４乃至１９のいずれか一項に記載の超音波方法。
【請求項２１】
超音波撮像データからがん悪性度マップのマップ画素を表すローカルフィーチャのセットを抽出するステップと、
前記がん悪性度マップの対応するマップ画素に対するがん悪性度を生成するようにがん格付け分類器を使用してローカルフィーチャの各セットを分類するステップと、
を有するがん悪性度マッピング方法を実行するように電子データ処理装置により読取可能かつ実行可能な命令を記憶する不揮発性記憶媒体において、
前記がん悪性度マップが、それぞれのマップ画素に対して生成されたがん悪性度に等しいマップ画素値を持つ前記マップ画素を有する、
不揮発性記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
以下は、広くは、腫瘍診断及び治療分野、生検及び組織標本収集分野、画像ガイド医療処置分野、及び関連分野に関する。これは、特に前立腺がん診断及び治療を参照して記載されるが、肝臓がん、又は乳がん等のような他のタイプのがんの診断及び治療において応用を見つける。
【背景技術】
【０００２】
２０１４年現在、前立腺がんは、米国において、男性における最も一般的なタイプのがんであり、二番目の主ながん関連の死因である。毎年、２３００００を超える米国人男性が、前立腺がんと診断され、３００００近くの人が、この疾病で亡くなっている。血液中の前立腺特異抗原（ＰＳＡ）の増大したレベル、触診可能な結節、前立腺がんの家族歴又は前立腺の超音波画像において見られる低エコー領域が存在する場合に、前立腺がんが、疑われる。しかしながら、血液ＰＳＡ試験結果は、高い偽陽性率を生じ、これは、関連した起こりうる合併症を伴う不要な治療処置をもたらす可能性がある。
【０００３】
より確実な前立腺がん診断は、従来、経直腸超音波撮像によってガイドされる直腸ツールを用いて取得された生検標本の組織病理学的解析の方法によるものである。不幸なことに、前立腺がんは、散在している悪性領域として形成する傾向にあり、この試験に対する偽陰性率は、貧弱な標的化のために高い。この意味での「偽陰性」には、完全なミス（癌がないと誤って示す）、又は前立腺に実際に存在する最も悪性のがんより低いがん悪性度を含む。より具体的には、経直腸超音波ガイド生検は、典型的には、低い感度を持ち、４０％乃至６０％の範囲の陽性の予測値は、効果的な治療計画及び標的化を妨げる。生検は高価かつ侵襲的であり、合併症の可能性がある。したがって、反復生検は、ワークフローの観点から非効率的であるだけでなく、望ましくない。
【０００４】
前立腺癌の診断が行われた後、適切な治療法が展開される。高密度焦点式超音波（ＨＩＦＵ）、凍結療法、無線周波数アブレーション（ＲＦＡ）、又は光線力学療法（ＰＤＴ）のような焦点療法は、一般に、前立腺がんの散在領域を標的にしながら前立腺器官に最小限の影響を及ぼすように設計された最小侵襲技術である 。しかしながら、典型的な前立腺がんの散在性質は、焦点療法による高い悪性度のがん領域の効果的な標的化を困難な課題とする。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
以下は、上記の課題等に対処する改良されたシステム及び方法を開示する。
【課題を解決するための手段】
【０００６】
開示される一態様において、超音波システムは、超音波撮像データを取得するように構成された超音波撮像装置と、（ｉ）前記超音波撮像データから、がん悪性度（cancer grade）マップのマップ画素を表すローカルフィーチャのセットを抽出し、（ｉｉ）前記がん悪性度マップのマップ画素に対するがん悪性度を生成するようにがん格付け（grading）分類器を使用して前記ローカルフィーチャのセットを分類することにより、がん悪性度マップを生成するようにプログラムされた電子データ処理装置と、前記がん悪性度マップを表示するように構成された表示コンポーネントとを有する。
【０００７】
他の開示される態様において、超音波方法は、超音波撮像データを取得するステップと、前記超音波撮像データから超音波画像を生成するステップと、前記超音波撮像データから抽出されたローカルフィーチャのセットにがん格付け分類器を適用することにより前記超音波撮像データからがん悪性度マップを生成するステップと、（ｉ）前記がん悪性度マップ並びに（ｉｉ）前記超音波画像及び前記がん悪性度マップを結合する融合画像の少なくとも１つを表示するステップとを有する。
【０００８】
他の開示される態様において、不揮発性記憶媒体は、がん悪性度マッピング方法を実行するように電子データ処理装置により読取可能かつ実行可能な命令を記憶し、前記方法は、超音波撮像データから、がん悪性度マップのマップ画素を表すローカルフィーチャのセットを抽出するステップと、前記がん悪性度マップの対応するマップ画素に対するがん悪性度を生成するようにがん格付け分類器を使用してローカルフィーチャの各セットを分類するステップとを有する。前記がん悪性度マップは、それぞれのマップ画素に対して生成されたがん悪性度に等しいマップ画素値を持つマップ画素を有する。
【０００９】
１つの利点は、超音波により取得されたがん悪性度マップを提供することにある。
【００１０】
他の利点は、このようながん悪性度マップをリアルタイムで提供することにある。
【００１１】
他の利点は、このようながん悪性度マップを使用する改良された生検標本収集を提供することにある。
【００１２】
他の利点は、このようながん悪性度マップを使用する改良されたがん治療標的化を提供することにある。
【００１３】
所定の実施例は、前述の利点の１つ、２つ、それ以上、又は全てを提供してもよく、１つも提供しなくてもよく、及び／又は本開示を読み、理解すると当業者に明らかになるように他の利点を提供してもよい。
【００１４】
本発明は、様々なコンポーネント及びコンポーネントの組み合わせ、様々なステップ及びステップの組み合わせの形を取りうる。図面は、好適な実施例を説明する目的のみであり、本発明を限定すると解釈されるべきではない。
【図面の簡単な説明】
【００１５】
【図１】ここに開示されるがん悪性度マップを提供する経直腸超音波システムを図示する。
【図２】ｂモード超音波画像上に重ねられたがん悪性度マップを表示することを含む、図１のシステムを使用して適切に実行される超音波撮像方法を図示する。
【図３】図１のシステムにおいて採用されるがん格付け分類器を生成するように、コンピュータ又は他の電子データ処理装置により適切に実行される、オフライン処理を図示する。
【発明を実施するための形態】
【００１６】
前立腺がんの格付けは、典型的には、経直腸超音波ガイド生検により取得された標本を使用する組織病理学による。しかしながら、超音波は、典型的には、（良くても）前立腺の疑わしい領域の場所を示すが、これらの領域のがん悪性度を（又はこれらがそもそもがん性であるか否かさえ）決定することができない。したがって、生検が前立腺に存在する最も悪性のがんの標本を取得する保証は、存在しない。更に、この処置の経直腸性質は、実際に収集されることができる標本の数を制限する傾向にある。繰り返される経直腸生検処置も、望ましくない。
【００１７】
ここに開示されるのは、経直腸超音波撮像により取得された既存の生の「ＲＦ」時系列データを利用して、経直腸超音波画像からがん悪性度マップを生成することである。（用語「ＲＦ」は、従来、「無線周波数」を示す。超音波との関連で、撮像超音波パルスは、典型的には無線周波数と比較可能であるメガヘルツ範囲である音波周波数であり、したがって、超音波との関連で用語「ＲＦ」時系列も同様である。）典型的な超音波撮像において、超音波パルスは、毎秒３０乃至５０回のオーダで印加され、したがって、毎秒３０乃至５０の輝度画像（いわゆる２Ｄ超音波撮像における「ｂモード」画像）を生成する。これらの画像が、前記ｂモード画像の各画素に対して、前記ＲＦ時系列から対応する時間変化する信号を生成することができるように、組織加熱又は音響機械効果のような様々な可能なメカニズムによって経時的に変化しうることが既知である。これらの時間変化する信号は、一部の例において組織タイプと関連するように示されている。
【００１８】
ここに開示されるように、画素レベルのＲＦ時系列情報は、２Ｄ画像（例えばｂモード画像）又は（３Ｄ超音波システムに対する）３Ｄ画像上に重ねられることができるがん悪性度マップを生成するのに使用される。ＲＦ時系列データにおける組織コントラストをもたらす（現在）貧弱に理解される身体的メカニズムの観点から、機械学習アプローチが、開示される実施例において採用される。このために、テクスチャ又はウェーブレットのようなローカルフィーチャが、各マップ画素に対して抽出される。これらのマップ画素は、前記超音波画像の画素解像度であってもよく、又はより粗いマッピング解像度であってもよい。（更に、個々で使用される用語「画素」は、「絵素」を示し、前記ＲＦ時系列データが２Ｄ超音波システムを使用して取得されるのか又は３Ｄ超音波システムを使用して取得されるのかに依存して、２Ｄ画素又は３Ｄ画素のいずれであってもよい。）前記ローカルフィーチャは、各マップ画素を表すフィーチャベクトルを形成し、前記マップ画素に対してがん悪性度を割り当てるようにがん格付け分類器に入力される。前記がん格付け分類器（又は複数の分類器）は、組織病理学的悪性度が割り当てられている実際の生検場所の超音波画像を有するラベル付けされたトレーニングデータに対する機械学習を使用してトレーニングされる。前記がん悪性度マップは、ｂモード画像上にカラーオーバレイとして重ねられてもよく、又は他の形で前記超音波画像と融合されてもよい。
【００１９】
前記がん悪性度マップ生成は、高速である。前記トレーニングされた分類器は、計算的に効率が良く、前記トレーニングは、オフラインで実行されることができる。前記超音波がん悪性度マッピングは、また、従来の（例えばｂモード）超音波撮像中に既に生成された「生の」ＲＦ時系列データも使用する。したがって、開示されるがん悪性度マッピングは、リアルタイム超音波撮像中に容易に使用される。前記がん悪性度マップは、これにより、直腸プローブ再配置、不用意な患者の動き、又は超音波撮像設定（例えば解像度、焦点）の変化等を説明するようにリアルタイムで更新される。経直腸超音波ガイド生検処置中に使用されることに加えて、このアプローチは、小線源療法シード埋め込み、又は反転変調放射線療法（ＩＭＲＴ）に対する計画画像の取得等の間の使用に対して期待される。
【００２０】
ＲＦ時系列データが、前記がん悪性度マッピングデータを生成する超音波撮像メカニズムとして開示されるが、より一般的に、（低周波数の超音波パルスが組織振動を誘起するように印加される）エラストグラフィのような他のコントラストメカニズムにより生成されたマッピングデータが、使用されてもよい。更に、図示の実施例が、前立腺がん診断及び治療に対して経直腸超音波撮像を使用するが、このアプローチは、肝臓又は乳がんのような他のタイプのがんのリアルタイム格付けに対して容易に使用される。
【００２１】
図１を参照すると、経直腸超音波システムは、患者１４の直腸に挿入され、ケーブルにより超音波撮像システム１０と接続される直腸超音波プローブ１２を持つ超音波撮像システム１０（例えば、オランダ、アイントホーフェンのKoninklijke Philips N.V.から入手可能な図示されたＥＰＩＱ（登録商標）超音波撮像システム、又は他の業務用の又は特注の超音波撮像システム）を含む。（図１が、図式的表現であり、超音波プローブ１２が、実際には、患者の直腸に挿入される場合には視界から隠されると理解される）。例示的な超音波プローブは、生検標本を収集する一体化された生検針１６を含み、代替的には、別個の生検ツールが使用されてもよく、又は前記経直腸超音波システムは、生検ツールを使用しない他の処置に対して、例えばＩＭＲＴ計画画像取得中に、使用されてもよい。前記経直腸超音波撮像処置に対して、患者１４は、（図１に図示されるように）適切な枕又は他のサポート（図示されない）を持つ図示された患者寝台又は支持台上の彼の側に位置する。例示的な超音波撮像システム１０は、超音波画像を表示する表示コンポーネント２０、及びユーザインタフェースディスプレイ２２及びユーザ入力制御部２４（例えばボタン、トラックボール等）のような１以上のユーザインタフェースコンポーネントを含む。
【００２２】
超音波撮像システム１０は、更に、ＲＦ時系列超音波撮像データを収集し、前記ＲＦ時系列超音波撮像データの各フレームから従来の輝度（ｂモード）画像３４を生成するようにプログラムされたＲＦ時系列撮像データ取得コントローラ３２を実現する、図１に図示されたマイクロプロセッサ、マイクロコントローラ、又は他の電子データ処理コンポーネント３０を含む。典型的なセットアップにおいて、コントローラ３２は、前記超音波プローブに選択された周波数（典型的にはメガヘルツ乃至数十メガヘルツ範囲であるが、この範囲の外の周波数及び／又はマルチ周波数パルスも考えられる）における音波パルス（又はパルスパケット）を投入させ、各このようなパルス又はパルスパケットに応答して（「フレーム」として知られる）撮像データを取得させる。このようにして、典型的には毎秒３０乃至５０フレームを含むフレームのＲＦ時系列が、取得される（他のフレームレートも考えられる）。各フレームのデータは、二次元画像、例えばｂモード画像を形成するように処理されることができ、又は３Ｄ超音波プローブの場合に、３Ｄ輝度画像を形成するように処理されることができる。前記ｂモード画像は、（深度と関連する）エコー遅延及び（例えば超音波プローブ１２のフェーズドアレイ又はビーム形成設定に基づいて、又は前記プローブとともに含まれる物理的レンズを使用して決定される）方向に基づいて生成される。前記ｂモード画像は、例えば、表示コンポーネント２０上に表示され、前記ｂモード画像がリアルタイム画像であるようにフレームごとに若しくはフレームのセットごとに（例えば選択された数の連続したフレームを平均して）更新されてもよい。
【００２３】
前記ＲＦ時系列超音波撮像データは、がん悪性度マップを生成するように、超音波撮像システム１０の電子データ処理コンポーネント３０の適切なプログラミングにより実現される、がん悪性度マッパコンポーネント４０によっても処理される。がん悪性度マップ４２は、（ｂモード画像３４と同じ解像度、又はより粗い解像度であってもよく、例えば各マップ画素が、ｂモード画像画素の隣接するｎ×ｎのアレイ、例えばｂモード画像画素の３×３のアレイ又はｂモード画素の１６×１６のアレイ等に対応してもよい）マップ画素のアレイに分割される。各マップ画素に対して、がん悪性度マッパ４０のフィーチャ抽出器４４は、前記マップ画素を表すフィーチャベクトルを生成し、このフィーチャベクトルは、前記マップ画素に対するがん悪性度を生成するようにがん格付け分類器（又はがん格付け分類器のセット）４６に入力される。前記がん悪性度は、好ましくは、前立腺がんの組織病理学的格付けに対して一般に使用されるグリソンスコアのような標準的ながん格付けスキームと一致する。グリソンスコアリングシステムは、グレード１（正常な前立腺細胞、すなわち良性）から、細胞の増加する割合が不規則であるグレード２乃至４を通って、細胞が概ね異常であり、ランダムに秩序付けられる最高のグレード５までの範囲を取る。変形アプローチにおいて、２つの最も一般的な細胞パターンが格付けされ、２つのスコアが、２乃至１０のグリソンスコアを生成するように結合される。超音波撮像システム１０は、細胞レベルで撮像することができないが、しかしながら、がん格付け分類器４６は、これらの生検標本に対する組織病理学的結果と対にされた生検標本場所の超音波画像領域を有するトレーニングデータを使用して以前にトレーニングされており、これにより、分類器４６の出力は、前記マップ画素の場所から取られた標本の組織病理学的分析により割り当てられるがん悪性度と高い相関を持つ。一部の実施例において、前記分類器は、単純化された又は減少された格付けスケールを採用してもよく、例えば、がん格付け分類器４６は、値３が、グリソンスコアのグレード２乃至４に及ぶように１、３又は５の値を出力してもよい。
【００２４】
この超音波ベースのがん格付けのアプローチは、がん悪性度が増大するにつれて増大する細胞異常性及び細胞秩序の増大されたランダム性が、超音波に誘起される組織加熱の変化及び組織の音響機械的応答の変化を生じやすいという認識を前提としている。このような現象は、前記ＲＦ時系列における時間変化を生じると理解されるので、前記ＲＦ時系列超音波データは、異なるがん悪性度の悪性組織に対するコントラストを示すことを合理的に期待される。同様に、超音波エラストグラフィにおいて、異なるがん悪性度の悪性組織が、前記がん悪性度が増大するにつれて増大された細胞障害及び細胞レベルの変化により異なる弾性を示すことが期待され、したがって、超音波エラストグラフィは、異なるがん悪性度の悪性組織に対するコントラストを示すことを合理的に期待される。開示された超音波がん格付け技術は、おおよそ前記マップ画素解像度の解像度でがん格付けを提供するがん悪性度マップ４２を生じるようにこのようながん悪性度コントラストを利用する。
【００２５】
超音波撮像システム１０の電子データ処理コンポーネント３０は、更に、超音波撮像システム１０の表示コンポーネント２０に適切に表示される融合画像を生成するためにｂモード画像３４及びがん悪性度マップ４２を（必要であれば）空間的に位置合わせする空間的位置合わせ及び／又は画像融合コンポーネント４８を実現するようにプログラムされる。空間的位置合わせは、ｂモード画像３４が前記ＲＦ時系列データから生成される態様に依存して、必要とされるかもしれないし、されないかもしれず、これが、サイズ変更又は再サンプリングを伴う場合、空間的位置合わせが必要とされうる。前記画像融合は、２つの画像３４、４２を結合するいかなる適切なアプローチを採用することもできる。１つのアプローチにおいて、前記がん悪性度（例えば、グリソンスケールのグレード１乃至５）は、グレード１＝透明、グレード２＝黄色、グレード３＝黄色がかった橙色、グレード４＝橙色及びグレード５＝赤のような色コードを割り当てられる（これらは、単に説明用の色符号化である）。色符号化されたがん格付けマップは、例えば、（前記がん格付けマップオーバレイの透明度を制御するアルファ値が、随意に、ユーザ選択可能パラメータでありうる）アルファチャンネルを使用して半透明オーバレイとしてｂモード画像３４と適切に融合される。
【００２６】
２つの画像３４、４２を融合する他の考えられる画像処理技術は、以下の通りである。
【００２７】
画像融合は、例示的に図１に記載されるが、表示コンポーネント２０上に並べてｂモード画像３４及びがん格付けマップ４２を表示するような、他のディスプレイ提示形式が、使用されてもよい。この表示は、随意に、他のフィーチャを含んでもよく、例えば、生検針１６が、前記超音波画像に表れることを可能にする追跡フィーチャを含む場合に、その場所が、前記融合画像上に示されてもよい。このような場合に、可聴インジケータは、随意に、がん格付けマップ４２により示されるように、いつ追跡された生検針先端が悪性度の高いがんの領域に入るかを示すように提供されることができる（例えば前記可聴インジケータは、周波数及び／又は音の大きさが、前記針により貫通されるがん悪性度の増大とともに増大するビープ音であることができ、点滅するインジケータ光は、同様に作動されることができる）。更に、２Ｄ超音波撮像が記載されているが、３Ｄ撮像に対する拡張が容易であり、この場合、表示される画像は、三次元レンダリング、投影画像、又は他の画像表現でありうる。
【００２８】
図１を参照し続け、更に図２を参照すると、図１のシステムにより適切に実行される処理が、記載されている。動作Ｓ１において、取得コントローラ３２は、ＲＦ時系列超音波データを取得するように超音波撮像システム１０及びプローブ１２を作動する。これらのデータは、ｂモード画像３４を生成するように動作Ｓ２において処理される。（代替的には、他のタイプの画像表現が生成されてもよい。）動作Ｓ３において、フィーチャ抽出器４４は、各マップ画素に対するフィーチャのセット（すなわちベクトル）を抽出するように使用される。この処理は、（１）前記時系列データから前記画像の各画素に対する値の時系列を生成すること、（２）前記マップ画素を形成するように画像画素の隣接するｎ×ｎのグループを連結すること（ｎ＝１、すなわち前記マップ画素が前記画像画素と同じサイズである場合を除く）、及び（３）各マップ画素（すなわち、画像画素の各ｎ×ｎのグループ）に対して、前記フィーチャのセットを抽出することを伴う。前記マップ画素フィーチャは、ローカルフィーチャであるべきであり、画像画素のｎ×ｎのグループに関連付けられたローカルフィーチャの各セットは、マップ画素を形成する。いくつかの適切なローカルフィーチャは、実例として、（Haralick et al., "Textural Features for Image Classification", IEEE Transactions on Systems, Man, and Cybernetics vol. SMC 3, No. 6 pp. 610 621, 1973の標準的なテクスチャフィーチャ、又はその変形のような）テクスチャフィーチャ、ウェーブレットベースのフィーチャ、及び／又はスペクトルフィーチャを含む。動作Ｓ３の出力は、各マップ画素を表す（すなわち関連付けられた）フィーチャセット（すなわちフィーチャベクトル）ｘである。動作Ｓ４において、トレーニングされたがん悪性度分類器４６は、前記マップ画素に対するがん悪性度を生成するように各マップ画素のフィーチャベクトルｘに使用され、これらのマップ画素がん悪性度は、この場合、がん悪性度マップ４２を集合的に規定する。動作Ｓ５において、空間的位置合わせ／画像融合コンポーネント４８が、ｂモード画像３４及びがん悪性度マップ４２を（必要に応じて）空間的に位置合わせし、２つの画像３４、４２を融合するように使用され、動作Ｓ６において表示コンポーネント２０上に表示される融合画像を形成する。前記空間的位置合わせは、必要である場合、剛体又は弾性位置合わせを使用して画像３４、４２をアラインすることを適切に伴う。ｂモード及びＲＦモダリティに対して、ＲＦからｂモードへの既知の処理及びスキャン変換ステップは、前記位置合わせに対して使用されることができる。前記空間的位置合わせは、ｂモード画像３４とアラインするようにがん格付けマップ４２を調節する、又はその逆を行うことができる。フィーチャ抽出及び分類動作Ｓ３、Ｓ４を実行する前にｂモード画像３４又は前記取得されたＲＦ時系列データを調節するように前記空間的位置合わせを実行することも考えられる（すなわち、前記ＲＦ時系列データからがん格付けマップ４２を生成する前に前記ＲＦ時系列データ及び前記ｂモード画像を空間的に位置合わせすることが考えられる）。
【００２９】
図２に示されるループする矢印Ｓ７により示されるように、この処理は、ｂモード画像３４、がん格付けマップ４２、及びこれらの融合をリアルタイムで更新するように反復的に繰り返されてもよい。前記ＲＦ時系列は、迅速に、例えば毎秒３０乃至５０フレームで取得され、このようなリアルタイム更新を容易に実行可能にする。例示的な図２は、ループＳ７の各反復において同時に更新されるｂモード画像３４及びがん悪性度マップ４２の両方を示すが、これは、必要ではない。例えば、ｂモード画像３４は、がん格付けマップ４２より頻繁に更新されることができ、例えば前記ｂモード画像は、１０フレームごとに更新されることができるのに対し、がん格付けマップ４２は、１００フレームごとに更新されることができる。変形例の重複技術は、同じレートで前記ｂモード及びがん悪性度マップを更新することを容易化するように採用されることができる。例えば、１００のＲＦ時系列フレームが、悪性度マップを計算するのに使用される場合、前記悪性度マップ表示は、ＲＦフレーム＃１乃至＃１００を使用して、ｂモード画像＃１０１において開始することができる。次いで、ｂモード画像＃１０２において、ＲＦフレーム＃２乃至＃１０１から計算された悪性度マップが、表示され、以下同様である。したがって、（最初の１００のＲＦフレームを取得するための）がん悪性度マップ４２の表示を開始する際の最初の遅延の後に、がん悪性度マップ４２の後続の更新は、ｂモード画像３４の表示の更新と同じレートである。（超音波プローブ１２が移動された場合、がん悪性度マップ４２が再び同期される前に約１００のＲＦフレームの取得に対応する遅延が存在し、加えて、この重複技術は、十分に速い悪性度マップ推定に対して予測される）。
【００３０】
図３を参照すると、がん格付け分類器（又は複数の分類器）４６をトレーニングするのに機械学習を使用する例示的な方法が記載されている。この処理は、随意に、オフラインで、すなわち、超音波システム１０のマイクロプロセッサ、マイクロコントローラ、又は他の電子データ処理コンポーネント３０以外のコンピュータ６０により実行される。例えば、コンピュータ６０は、デスクトップコンピュータ、ノートブックコンピュータ、ネットワークベースのサーバコンピュータ、又はクラウドコンピュータシステム等でありうる。図３の処理は、トレーニングされた分類器４６を提供するために、図２を参照して記載された患者処置の前に実行される。
【００３１】
図３のトレーニングは、ラベル付けされたトレーニング標本６２に対して動作する。各ラベル付けされた標本は、（例えば前記ＲＦ時系列データから生成されたｂモード画像上の）同じ抽出が識別された生検の場所を持つ生検ＲＦ時系列超音波データを含む。各生検場所は、組織病理学的がん悪性度、すなわち、組織標本の病理学的分析により前記場所から抽出された組織標本に割り当てられるがん悪性度でラベル付けされる。ラベル付けされたトレーニング標本６２は、経直腸超音波ガイド生検の後に続いて標本の組織病理学的格付けを受けた過去の患者に対するデータであり、前記データに対して、前記生検中に取得されたＲＦ時系列超音波データが保存されている。トレーニング標本６２の各生検標本抽出に対して、医師は、前記場所の記録を提供するように前記ｂモード画像上の場所を適切にラベル付けする。トレーニング標本６２を作り上げるデータの過去の患者は、好ましくは、陽性の標本の統計的に代表的なサンプリングを提供するように選択され、すなわち、組織病理学的結果により示される様々な段階の前立腺がんを持つ患者である。トレーニング標本６２は、好ましくは、前立腺がんではない患者のサンプリングをも含む（陰性の標本、これらは、また又は代替的に、陰性の標本が、組織病理学ががんを示さない、すなわち１のグリソンスコアの前立腺器官の領域から取り出された生検標本を構成する場合に、前立腺がんの患者により提供されてもよい）。
【００３２】
動作Ｓ１２において、前記ＲＦ時系列データは、各生検場所を囲むマップ画素に対するフィーチャセット（すなわちフィーチャベクトル）を生成するように処理される。動作Ｓ１２は、図２の動作Ｓ３に適切に対応し、例えば同じマップ画素解像度及び同じフィーチャのセット、すなわち同じフィーチャベクトルである。代替的なアプローチにおいて、前記フィーチャのセットは、図３の機械学習トレーニング処理の一部として選択され、この場合、前記処理は、動作Ｓ３により抽出されたフィーチャベクトルを作り上げるローカルフィーチャを選択するオプションの動作Ｓ１４を含む。このようなフィーチャ選択は、動作Ｓ３において使用されたフィーチャベクトルを形成する最終的なフィーチャセットを形成するように最初のフィーチャセットの冗長フィーチャを識別し、取り除くように、手動で又は自動的に、例えば相互情報、相関、又は同様の統計を使用して、実行されることができる。他の適切なフィーチャ選択アルゴリズムは、全数検索、遺伝アルゴリズム、又は前進若しくは後進消去等を含む。
【００３３】
ＲＦ時系列超音波撮像データから抽出されたローカルフィーチャの場合、通常の経直腸超音波撮像ガイド生検処置は、典型的には、（ｂモード画像がＲＦ時系列データから生成されるために）ｂモード撮像が実行されるので、そのうちに必要なＲＦ時系列超音波撮像データを取得する。他方で、動作Ｓ３が、エラストグラフィ撮像データのような他のタイプの超音波撮像データからフィーチャを抽出する場合、ラベル付けされたトレーニング標本６２の超音波データは、対応するローカルフィーチャのトレーニングセットが前記トレーニング超音波撮像データから抽出されることを可能にするために必要なタイプの超音波データ（例えばエラストグラフィ撮像データ）を含むことを必要とすると理解される。
【００３４】
動作Ｓ１２及びオプションの動作Ｓ１４の出力は、生検場所に対応する各マップ画素を表すフィーチャベクトルである。（前記生検場所が識別される解像度に依存して、前記生検場所に広がる複数のマップ画素が存在しうる。）各々が対応する抽出された組織標本に対する組織病理学的がん悪性度でラベル付けされたこれらのフィーチャベクトルは、ラベル付けされたトレーニングセット６４を形成する。
【００３５】
動作Ｓ１６において、がん格付け分類器４６は、このトレーニングセット６４に対してトレーニングされる。前記トレーニングは、セット６４の入力トレーニングフィーチャベクトルに対するがん格付け分類器６４の出力と、対応する組織病理学的がん悪性度ラベルとの間の誤差を最小化するようにがん格付け分類器４６のパラメータを最適化する。がん格付け分類器４６は、例えば５つのグリソンスコアに対応する離散化された出力１乃至５を持つ、単一のマルチラベル分類器を有してもよい。代替的には、がん格付け分類器４６は、各々が異なるがん悪性度に対するバイナリ分類器のセットを有してもよく、例えば、グリソンスコア４のバイナリ分類器は、最適に、グリソンスコア４のラベルを持つトレーニングフィーチャベクトルに対して「１」を出力し、他のラベルを持つトレーニングベクトルに対して「０」を出力するようにトレーニングされる。一部の実施例において、分類器４６は、決定木のアンサンブル（ときどきランダムフォレストと称される）のような分類器のアンサンブルである。ある適切な分類器は、線形回帰、ロジスティック回帰、サポートベクトルマシン、及び決定木分類器等を含むが、これらに限定されない。アンサンブル分類器の使用の場合、各分類器の悪性決定の大多数のようなマップ画素の悪性度値が、算出されることができる。
【００３６】
多くのこのような分類器は、連続的な値を出力する。グリソンスコアのような離散的ながん悪性度を生成するために、閾値動作が、前記分類器の連続値出力に対して実行されることができ、これにより、前記マップ画素値が、離散値である。代替的には、閾値処理は実行されず、前記マップ画素は、連続値の分類器出力を直接的に割り当てられる。この場合、画像融合動作４８は、随意に、以前に記載されたような離散化された色ではなく、連続的な分類器出力にマッピングされる連続的な色のスペクトルを使用して色符号化を実行してもよい。
【００３７】
結果として生じるトレーニングされたがん格付け分類器４６（又はトレーニングされたパラメータ）は、がん悪性度分類動作Ｓ４を実行する際にマイクロプロセッサ、マイクロコントローラ、又は他の電子データ処理コンポーネント３０により使用するために超音波システム１０に適切にロードされる。
【００３８】
図１のシステムは、例えば、経直腸プローブ１２が前立腺器官の画像を取得するのに使用される、リアルタイム超音波撮像システム１０を含む。画像は、ｂモード撮像、ＲＦデータ、及びエラストグラフィ、又は後方散乱係数推定、若しくは減衰推定等のような他のＲＦデータに基づく方法を含むが、これらに限定されない。前記ＲＦデータは、従来のｂモード撮像に対してがん組織に関する追加の上方を提供する。生のＲＦ時系列データをｂモード画像に変換するのに伴う信号処理の様々なステップにより一部の情報が失われることが認識される。ここに開示されるように、前記超音波データ（例えばＲＦ時系列データ、及び／又はエラストグラフィデータ等）を使用して、がん悪性度の推定は、前立腺における各マップ画素又は領域の悪性度を推定するのにパターン認識及び機械学習技術を使用することにより実行される。各ボクセル又は領域（すなわちマップ画素）に対するがん悪性度が計算され、がん悪性度マップ４２が形成される。がん悪性度マップ４２は、前立腺のｂモード画像上に重ねられることができる、又は超音波装置１０が３Ｄ超音波撮像データを取得する場合に３Ｄにレンダリングされることができる。がん悪性度マップ４２は、プローブ１２又は生検装置１６をより良く配置するように前記超音波撮像及び生検システムにより使用されることができる。一度超音波プローブ１２が、特定の場所に移動されると、超音波撮像システム１０は、がん悪性度値を更新するように、がん悪性度マッパ４０により格付けされる更新された超音波画像を取得し、がん悪性度マップ４２は、これにより、適切に更新される。この処理は、がん悪性度マップ４２により示される高いがん悪性度の前立腺領域が識別されるまで、リアルタイムで繰り返されることができる。生検応用に関連して、前記識別された高いがん悪性度の前立腺領域は、生検標的として選択され、バイオプシーガン又はツール１６は、前記高い悪性度の領域から組織標本を取得するようにこの場所までガイドされる。
【００３９】
同様のワークフローは、標的療法に対しても考えられる。この応用において、前記高い悪性度のがんが識別され、治療ツール（例えば、小線源療法の場合に放射性シードを送達する針、又は無線周波数アブレーション針等）に対する標的として選択される。小線源療法の場合に、例えば、より多数のシードが、高い悪性度であるとがん悪性度マップ４２において示される場所に配置されてもよく、より少数のシードが、より低い悪性度として示される場所に配置されてもよい。ＩＭＲＴ計画応用において、がん悪性度マッパ４０は、計画画像（例えば、コンピュータ断層撮影、すなわちＣＴ、計画画像及び代替的には、計画ＣＴデータを増大する超音波ＲＦ時系列を持つ）の取得中に使用される。がん悪性度マップ４２は、基準マーカ又は解剖学的マーカ等を使用して前記計画画像と空間的に位置合わせされ、アラインされたがん悪性度マップ４２は、前記計画画像における高い悪性度のがん領域又は複数の領域をセグメント化する唯一の又は追加の上方を提供する。
【００４０】
がんを格付けする目的で、例示的な実施例は、超音波で生成されたがん悪性度マップ４２により示される最も高いがん悪性度の領域の標的サンプリングを実行するために生検処置をガイドするツールとしてがん悪性度マッパ４０を使用する。このアプローチにおいて、がん悪性度マップ４２は、生検標本収集をガイドするのに役立つが、前記生検標本の組織病理学的分析により生成されるがん格付けは、臨床的使用（すなわち、診断及び治療をガイドする）に対して受け入れられた格付けとして機能する。この例示的なアプローチは、臨床的格付けが、がん専門医により良好に受け入れられる組織病理学的格付けであるという利点を持つ。
【００４１】
代替的な実施例において、超音波で生成されたがん悪性度マップ４２は、臨床的使用に対する格付けとして機能する。すなわち、このような実施例において、生検は実行されず、代わりに、がん専門医は、がん格付けとしてがん悪性度マップ４２を信頼する。このアプローチは、がん悪性度マップ４２により提供されるがん格付けの特異性及び感度が、臨床的要件を満たすことを要求し、これは、がん悪性度マップ４２により生成される悪性度を記録し、これを組織病理学的悪性度と比較することにより経時的に決定されることができ、これらが、経時的にかつ十分な統計とともに満足のいく一致を示す場合に、がん悪性度マップ４２は、単独で合理的に信頼されうる。このアプローチは、生検標本収集と後の病理学的分析及び報告との間の遅延並びに侵襲的生検処置を消去する利点を持つ。
【００４２】
例示的な前立腺がんの例は、このようなアプローチが、前立腺の超音波撮像において一般的に及び有効に使用されるので、例示的な経直腸超音波プローブ１２を使用する。しかしながら、前述のように、開示された超音波ベースのがん格付けアプローチは、他のタイプのがんを格付けするのに有用に使用されうる。がんのタイプに依存して、異なるタイプの超音波プローブが使用されうる。例えば、乳がん撮像において、表面超音波プローブが、好ましいかもしれない。
【００４３】
例示的な実施例において、がん悪性度マッパ４０は、超音波装置１０のコンポーネントであるマイクロプロセッサ、マイクロコントローラ、又は他の電子データ処理コンポーネント３０により実現される。これは、マイクロプロセッサ又はマイクロコントローラ３０が、超音波装置１０と一体化され、例えば一部の実施例において電子コントローラとして機能し、したがって生のＲＦ時系列データを含む取得された超音波データに対して直接的なアクセスを持ち、例えば、前記ｂモード画像上のオーバレイとしてがん悪性度マップ４２を表示するために超音波装置１０の画像表示機能と一体化されることができるので、有利である。しかしながら、代替的には、がん悪性度マッパ４０が、前記ＲＦ時系列データを含む超音波撮像データを受信し、がん悪性度マップ４２を表示する表示コンポーネントを含む（又は超音波装置１０の表示コンポーネント２０にアクセスする）異なる電子データ処理装置上で実現されることが考えられる。例えば、がん悪性度マッパ４０は、ＵＳＢケーブル又は他のデータ接続により超音波装置１０と接続されたノートブックコンピュータ上で実現されてもよい。このような実施例において、がん悪性度マッパ４０は、前述のようにリアルタイムでがん悪性度マップ４２を更新するように前記超音波撮像と同時に実行してもよく、又は、代替的には、がん悪性度マッパ４０は、超音波撮像セッションが完了された後に実行されてもよく、保存されたＲＦ時系列超音波データに対して動作する。
【００４４】
図１及び２を参照してここに開示されたもののような様々な超音波ベースのがん格付けアプローチが、これらの動作を実行するようにマイクロプロセッサ、マイクロコントローラ、又は他の電子データ処理コンポーネント３０により読取可能かつ実行可能である命令を記憶する不揮発性記憶媒体により実現されてもよい。同様に、図３を参照してここに開示されたもののような様々な分類器トレーニングアプローチは、オフライン分類器トレーニングを実行するコンピュータ又は他の電子データ処理コンポーネントにより読取可能及び実行可能である命令を記憶する不揮発性記憶媒体により実現されてもよい。このような不揮発性記憶媒体は、非限定的な実例として、ハードディスクドライブ若しくは他の磁気記憶媒体、フラッシュメモリ、読取専用メモリ（ＲＯＭ）若しくは他の電子記憶媒体、光学ディスク若しくは他の光学記憶媒体、又はこれらの様々な組み合わせ等を含んでもよい。
【００４５】
本発明は、好適な実施例を三泡して記載されている。修正及び変更は、先行する詳細な記載を読み、理解すると他者が思いつきうる。本発明が、添付の請求項又はその同等物の範囲内に入る限り全てのこのような修正及び変更を含むと解釈されることが意図される。
【図１】
2018516135.tif 000003
【図２】
2018516135.tif 000004
【図３】
2018516135.tif 000005
【国際調査報告】
2018516135.tif 000006
2018516135.tif 000007
2018516135.tif 000008
2018516135.tif 000009
2018516135.tif 000010
2018516135.tif 000011