(57)【要約】  本開示は、超音波画像を生成するシステムおよび方法を提供する。この方法は、超音波装置から出力される生の超音波信号を受信することと、生の超音波信号を強化超音波画像に変換する動作を実行することとを含む。強化超音波画像は、放射状対称フィルタリングによりさらに処理され、動径対称画像を生成することができる。強化画像と放射対称画像の両方は、医師が超音波データのみに基づいて肝臓がんまたは腎臓がんの診断を行うためのものである。【選択図】図１
2020508174.tif 000002
【特許請求の範囲】
【請求項１】
一時的ではない、コンピュータ読み取り可能メモリと、
１つ以上のプロセッサと、
プログラミング命令を含むコンピュータ読み取り可能な媒体とを備えた超音波画像を生成するシステムであって、前記プログラミング命令は、前記１つ以上のプロセッサで実行されたとき、前記システムが、
超音波装置から出力される生の超音波信号を取得し、
前記生の超音波信号をフィルタリングすることによりフィルタリングされた超音波画像データを生成して、信号ドメインから周波数ドメインに情報を変換し、
前記フィルタリングされた超音波画像データから画像の位相とエネルギーの特徴を抽出し、
前記画像の位相とエネルギーの特徴を周波数領域から信号領域に変換することにより、フィルタリングされた超音波信号を生成し、
組織界面からの後方散乱超音波信号に基づいて透過マップ推定を生成し、
前記透過マップ推定に基づいて前記フィルタリングされた超音波信号の前記画像位相とエネルギー特性を強化し、それに関連する軟組織の部分が強化超音波画像でより暗くなったり、より明るくなるようにする
システム。
【請求項２】
前記１つ以上のプロセッサによって実行されたときに、前記システムが、
強化超音波画像の暗い球面形状を分析することにより、強化超音波画像に対して放射状対称フィルタリングを実行し、
放射状対称画像を生成して、強化超音波画像の関心領域を識別する
プログラミング命令をさらに含む、請求項１に記載のシステム。
【請求項３】
前記１つ以上のプロセッサによって実行されたときに、前記システムが、
前記放射対称性画像に基づいて表面トポグラフィマップを出力させる
プログラミング命令をさらに含む、請求項２に記載のシステム。
【請求項４】
前記画像の位相およびエネルギーの特徴を抽出するステップがアルファスケールフィルタリングによって実行される、請求項１に記載のシステム。
【請求項５】
前記アルファスケールフィルタリングが周波数領域で実行される、請求項４に記載のシステム。
【請求項６】
前記アルファスケールフィルタリングが、以下の式によって定義される、請求項５に記載のシステム。
【数１】
2020508174.tif 000028
αは定数の微分パラメータ、σはフィルタのアルファスケールパラメータ、ｎｃはフィルタαから計算された単位正規化定数である。
【請求項７】
ｎｃが以下の式に従って計算される、請求項６に記載のシステム。
【数２】
2020508174.tif 000029
ｓはスケールパラメータで、αは微分パラメータである。
【請求項８】
前記透過マップ推定が、以下の式に基づいて軟組織における散乱および減衰効果を組み合わせることを含む、請求項１に記載のシステム。
【数３】
2020508174.tif 000030
ＵＳ（ｘ、ｙ）はローカルエネルギー画像、ＵＳＡ（ｘ、ｙ）は信号伝送マップ、ＵＳＥ（ｘ、ｙ）は強化超音波画像、αは、局所領域の軟組織のエコー源性を表す定数値である。
【請求項９】
減衰関数をイメージング深度の関数としてモデル化するランベルト・ベールの法則を使用して信号伝送マップＵＳＡ（ｘ、ｙ）を推定することにより、ＵＳＥ（ｘ、ｙ）が抽出される請求項８に記載のシステム。
【請求項１０】
前記ランベルト・ベールの法則が以下の式によって定義される、請求項９に記載のシステム。
【数４】
2020508174.tif 000031
ＵＳｏ（ｘ、ｙ）は初期強度画像、ＵＳＴ（ｘ、ｙ）は減衰強度画像、ηは減衰係数、ｄ（ｘ、ｙ）は超音波トランスデューサの表面からの距離である。
【請求項１１】
ＵＳＡ（ｘ、ｙ）が、以下の式によって定義される目的関数を最小化することによって取得される、請求項１０に記載のシステム。
【数５】
2020508174.tif 000032
ωはインデックスセット、ｏは要素単位の乗算、＊は畳み込み演算、Ｄｊは８つのキルシュフィルタと１つのラプラシアンフィルタで構成される高次微分フィルタのバンクを使用して取得され、Ｗｊは、以下の数式に従って計算された重み行列である。
【数６】
2020508174.tif 000033
【請求項１２】
ＵＳＥ（ｘ、ｙ）が以下の式に従って計算される、請求項１１に記載のシステム。
【数７】
2020508174.tif 000034
【請求項１３】
前記放射状対称フィルタリングが、以下の式によって指定される、請求項２に記載のシステム。
【数８】
2020508174.tif 000035
「＊」は畳み込み演算を示し、Ａｎは等方性ガウス関数であり、Ｆｎは次のように定義される。
【数９】
2020508174.tif 000036
【数１０】
2020508174.tif 000037
また、αは放射状の厳密性パラメータであり、ｋｎは異なる半径にわたるスケーリング係数である。
【請求項１４】
前記放射状対称画像が以下の式によって定義される、請求項２に記載のシステム。
【数１１】
2020508174.tif 000038
Ｓは考えられるすべての範囲にわたるすべての対称性の寄与の合計である。
【請求項１５】
コンピューティングデバイスにより、超音波装置から出力された生の超音波信号を取得し、
前記生の超音波信号をフィルタリングすることによりフィルタリングされた超音波画像データを生成して、情報を信号領域から周波数領域に変換し、
前記フィルタリングされた超音波画像データから画像位相とエネルギー特徴を抽出し、
前記画像の位相とエネルギーの特徴を周波数領域から信号領域に変換することにより、フィルタリングされた超音波信号を生成し、
組織界面からの後方散乱超音波信号に基づいて透過マップ推定を生成し、
前記透過マップ推定に基づいて前記フィルタリングされた超音波信号の画像位相とエネルギー特性を強化し、それに関連する軟組織の部分が強化超音波画像でより暗くなったり、より明るくなるようにすることを含む、超音波画像を生成する方法。
【請求項１６】
強化超音波画像内の暗い球面形状を分析することにより、前記強化超音波画像上で、コンピューティングデバイスにより放射状対称フィルタリングを実行し、
コンピューティングデバイスによって、放射状対称画像を生成して、前記強化超音波画像内の関心領域を識別することを更に含む、請求項１５に記載の方法。
【請求項１７】
前記放射状対称画像に基づいて表面トポグラフィマップを出力することをさらに含む、請求項１６に記載の方法。
【請求項１８】
前記後方散乱超音波信号が、組織との散乱および減衰を含む相互作用によって変調される、請求項１５に記載の方法。
【請求項１９】
前記フィルタリングされた超音波画像データを生成するステップがフーリエ変換によって実行される、請求項１５に記載の方法。
【請求項２０】
前記フィルタリングされた超音波信号を生成するステップが逆フーリエ変換によって実行される、請求項１５に記載の方法。
【請求項２１】
前記画像位相およびエネルギー特徴を抽出するステップがアルファスケールフィルタリングにより実行される、請求項１５に記載の方法。
【請求項２２】
前記アルファスケールフィルタリングが周波数領域で実行される、請求項２１に記載の方法。
【請求項２３】
前記アルファスケールフィルタリングが以下の式によって定義される、請求項２２に記載の方法。
【数１２】
2020508174.tif 000039
αは定数の微分パラメータ、σはフィルタのアルファスケールパラメータ、ｎｃはフィルタαから計算された単位正規化定数である。
【請求項２４】
ｎｃが以下の式に従って計算される、請求項２３に記載の方法。
【数１３】
2020508174.tif 000040
ｓはスケールパラメータで、αは微分パラメータである。
【請求項２５】
前記透過マップ推定が以下の式に基づいて軟組織における散乱効果と減衰効果を組み合わせることを含む、請求項１５に記載の方法。
【数１４】
2020508174.tif 000041
ＵＳ（ｘ、ｙ）はローカルエネルギー画像、ＵＳ  Ａ（ｘ、ｙ）は信号伝送マップ、ＵＳ  Ｅ（ｘ、ｙ）は強化超音波画像、αは局所領域の軟組織のエコー源性を表す定数値である。
【請求項２６】
減衰関数をイメージング深度の関数としてモデル化するランベルト・ベールの法則を使用して前記信号伝送マップＵＳ  Ａ（ｘ、ｙ）を推定することによりＵＳ  Ｅ（ｘ、ｙ）が抽出される、請求項２５に記載の方法。
【請求項２７】
前記ランベルト・ベールの法則が以下の式によって定義される、請求項２６に記載の方法。
【数１５】
2020508174.tif 000042
ＵＳｏ（ｘ、ｙ）は初期強度画像、ＵＳＴ（ｘ、ｙ）は減衰強度画像、ηは減衰係数、ｄ（ｘ、ｙ）は超音トランスデューサの表面からの距離である。
【請求項２８】
ＵＳＡ（ｘ、ｙ）が、以下の式によって定義される目的関数を最小化することによって得られる、請求項２７に記載の方法。
【数１６】
2020508174.tif 000043
ωはインデックスセット、ｏは要素単位の乗算、＊は畳み込み演算、Ｄｊは８つのキルシュフィルタと１つのラプラシアンフィルタで構成される高次微分フィルタのバンクを使用して取得され、Ｗｊは、以下の式に従って計算された重み行列である。
【数１７】
2020508174.tif 000044
【請求項２９】
ＵＳＥ（ｘ、ｙ）が以下の式に従って計算される、請求項２８に記載の方法。
【数１８】
2020508174.tif 000045
【請求項３０】
前記放射状対称フィルタリングが、以下の式によって指定される、請求項１６に記載の方法。
【数１９】
2020508174.tif 000046
「＊」は畳み込み演算を示し、Ａｎは等方性ガウス関数であり、Ｆｎは次のように定義される。
【数２０】
2020508174.tif 000047
【数２１】
2020508174.tif 000048
αは放射状の厳密性パラメータであり、ｋｎは異なる半径にわたるスケーリング係数である。
【請求項３１】
前記放射状対称画像が以下の式によって定義される、請求項１６に記載の方法。
【数２２】
2020508174.tif 000049
Ｓは考えられるすべての範囲にわたるすべての対称性の寄与の合計である。
【発明の詳細な説明】
【技術分野】
【０００１】
本出願は３５Ｕ．Ｓ．Ｃ．§１１９（ｅ）の下で２０１７年２月２７日に提出された米国仮特許出願第６２／４６４，０５８号に対する優先権を主張する。前述の出願は参照により本明細書に組み込まれる。
【０００２】
この文書は一般に画像処理に関する。より詳細には、この文書は、改善された癌診断および治療監視を促進する超音波画像処理のためのシステムおよび方法に関する。
【背景技術】
【０００３】
Ｃ型肝炎ウイルス（「ＨＣＶ」）および非アルコール性脂肪肝疾患（「ＮＡＦＬＤ」）は、北米における慢性肝疾患の２つの最も一般的な原因である。これはおそらく、米国で最も一般的なタイプの肝疾患である。臨床診療では、経腹的超音波は、その可用性、低コスト、および放射線被曝がないため、最初の画像診断法として最も広く使用されている。ＮＡＦＬＤは、線維症、肝硬変、肝がん、肝移植を必要とする肝不全、および死亡につながる可能性がある。腹部超音波検査では軽度の肝脂肪症を検出できず、単純な脂肪症、非アルコール性脂肪性肝炎（ＮＡＳＨ）、および肝線維症を区別できず、オペレータに依存し、腹腔内ガスの干渉を受け、肥満患者では画質が悪く技術的に困難である。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本開示は、強化超音波画像を生成するためのシステムおよび方法を提供する。本システムおよび方法によって処理された超音波画像は、医師が生検サンプルを必要とせずに超音波データのみに基づいて肝臓または腎臓病の診断を行うことを可能にする視覚的明瞭度のある軟組織パターンを示す。
【課題を解決するための手段】
【０００５】
一態様では、システムは、非一時的なコンピュータ可読メモリ、１つ以上のプロセッサ、およびプログラミング命令を含むコンピュータ可読媒体を含むことができる。プログラミング命令は、１つ以上のプロセッサで実行されると、システムが超音波デバイスから出力される生の超音波信号を取得し、生の超音波信号をフィルタリングして信号ドメインから周波数ドメインに情報を変換することにより、フィルタリングされた超音波画像データを生成する。システムはまた、フィルタリングされた超音波画像データから画像位相およびエネルギー特徴を抽出し、画像位相およびエネルギー特性を周波数領域から信号領域に変換することにより、フィルタリングされた超音波信号を生成し得る。システムは、組織インターフェースからの後方散乱超音波信号に基づいて透過マップ推定をさらに生成し、透過マップ推定に基づいてフィルタリングされた超音波信号の画像位相およびエネルギー特性を強化して、強化超音波画像では関連する軟組織の部分をより暗くしたり、より明るくしたりすることができる。
【０００６】
いくつかの実施形態では、システムは、強化超音波画像内の暗い球形を分析することにより強化超音波画像上で放射状対称フィルタリングを実行し、放射状対称画像を生成して強化超音波画像内の関心領域を識別することができる。いくつかの実施形態では、システムは、放射状対称画像に基づいて表面トポグラフィマップをさらに出力し得る。いくつかの実施形態では、フィルタリングされた超音波画像データを生成するステップは、フーリエ変換によって実行される。いくつかの実施形態では、フィルタリングされた超音波信号を生成するステップは、逆フーリエ変換によって実行される。
【０００７】
いくつかの実施形態では、画像位相およびエネルギー特徴を抽出するステップは、アルファスケールフィルタリングにより実行される。いくつかの実施形態では、アルファスケールフィルタリングは周波数領域で実行され、アルファスケールフィルタリングは以下の式で定義できる。
【０００８】
【数１】
2020508174.tif 000003
【０００９】
αは定数の微分パラメータ、σはフィルタのアルファスケールパラメータ、ｎｃはフィルタαから計算された単位正規化定数で、ｎｃは以下の数式に従って計算される。
【００１０】
【数２】
2020508174.tif 000004
【００１１】
ｓはスケールパラメータで、αは微分パラメータである。
【００１２】
いくつかの実施形態では、透過マップ推定は、以下の式に基づいて軟組織の散乱および減衰効果を組み合わせることを含む。
【００１３】
【数３】
2020508174.tif 000005
【００１４】
ＵＳ（ｘ、ｙ）はローカルエネルギー画像、ＵＳＡ（ｘ、ｙ）は信号伝送マップ、ＵＳＥ（ｘ、ｙ）は強化超音波画像、αは局所領域の軟部組織におけるエコー源性の定数値の代表値である。ＵＳＥ（ｘ、ｙ）は、減衰関数をイメージング深度の関数としてモデル化するランベルト・ベールの法則を使用して信号伝送マップＵＳＡ（ｘ、ｙ）を推定することにより抽出できる。ランベルト・ベールの法則は、以下の式で定義される。
【００１５】
【数４】
2020508174.tif 000006
【００１６】
ＵＳｏ（ｘ、ｙ）は初期強度画像、ＵＳＴ（ｘ、ｙ）は減衰強度画像、ηは減衰係数、ｄ（ｘ、ｙ）は超音波トランスデューサ表面からの距離である。ＵＳＡ（ｘ、ｙ）は、以下の式で定義される目的関数を最小化することにより得られる。
【００１７】
【数５】
2020508174.tif 000007
【００１８】
ωはインデックスセット、ｏは要素単位の乗算、＊は畳み込み演算、Ｄｊは８つのキルシュフィルタと１つのラプラシアンフィルタで構成される高次微分フィルタのバンクを使用して取得され、Ｗｊは、以下の式に従って計算された重み行列である。
【００１９】
【数６】
2020508174.tif 000008
【００２０】
ＵＳＥ（ｘ、ｙ）は、以下の式に従って計算される。
【００２１】
【数７】
2020508174.tif 000009
【００２２】
いくつかの実施形態では、放射状対称フィルタリングは、以下の式によって指定される。
【００２３】
【数８】
2020508174.tif 000010
【００２４】
「＊」は畳み込み演算を示し、Ａｎは等方性ガウス関数であり、Ｆｎは次のように定義される。
【００２５】
【数９】
2020508174.tif 000011
【００２６】
【数１０】
2020508174.tif 000012
【００２７】
また、αは放射状の厳密性パラメータであり、ｋｎは異なる半径にわたるスケーリング係数である。放射対称画像は、以下の式によって定義される。
【００２８】
【数１１】
2020508174.tif 000013
【００２９】
Ｓは考えられるすべての範囲にわたるすべての対称性の寄与の合計である。
【００３０】
別の態様では、超音波画像を生成する方法も提供される。この方法は、コンピューティングデバイスによって、超音波デバイスから出力される生の超音波信号を取得することと、生の超音波信号をフィルタリングして情報を信号領域から周波数領域に変換することによって、フィルタリング超音波画像データを生成することを含む。この方法は、フィルタリングされた超音波画像データから画像位相およびエネルギー特徴を抽出し、画像領域およびエネルギー特性を周波数領域から信号領域に変換することによりフィルタリングされた超音波信号を生成することも含む。この方法はさらに、組織インターフェースからの後方散乱超音波信号に基づいて透過マップ推定を生成し、それに関連する軟組織の部分が強化超音波画像ではより暗くなったり明るくなったりできるように、透過マップ推定に基づいてフィルタリング超音波信号の画像相およびエネルギー特性を強化することを含む。
【００３１】
いくつかの実施形態では、方法は、強化超音波画像内の暗球形状を分析し、強化超音波画像内の関心領域を識別する放射状対称画像を生成することにより強化超音波画像上で放射状対称フィルタリングを実行することをさらに含む。いくつかの実施形態では、この方法は、放射対称性画像に基づいて表面トポグラフィマップを出力することをさらに含む。
【図面の簡単な説明】
【００３２】
本解決策は、以下の図面を参照して説明され、図面全体を通して、同様の数字は同様のアイテムを表す。
【図１】強化超音波画像を生成するプロセスの例を示す。
【図２】放射対称画像を生成するプロセスの例を示す。
【図３】脂肪肝疾患を示す非健康な肝臓の生検（組織病理学）画像を示す。
【図４】図４Ａ〜Ｄ（集合的に「図４」）  正常／健康な肝臓のＢモード（元の）超音波画像（図４Ａ）および正常／健康な肝臓の強化超音波画像（図４Ｂ）、病変肝臓のＢモード（元の）超音波画像（図４Ｃ）と病変肝臓の強化超音波画像（図４Ｄ）を示す。
【図５】図５Ａ〜Ｅ（集合的に「図５」）  異なる画像タイプの比較を示す。図５Ａは正常／健康な肝臓のＢモード（元の）超音波画像を示し、図５Ｂは正常／健康な肝臓組織に対応する放射状対称画像の表面トポグラフィ視覚化を示し、図５Ｃは図５Ｂと比較して異なる放射状対称フィルタパラメータが適用されている正常／健康な肝臓組織に対応する放射状対称画像の表面トポグラフィ視覚化を示し、図５Ｄは正常／健康な組織に対応する強化超音波画像の表面トポグラフィ視覚画を示し、図５Ｅは図５Ｄと比較して異なるフィルタパラメータが適用されている正常／健康な組織に対応する強化超音波画像の表面トポグラフィ視覚化を示す。
【図６】図６Ａ〜Ｅ（集合的に「図６」）  異なる画像タイプの比較を示す。図６Ａは病変肝臓のＢモード（元の）超音波画像を示し、図６Ｂは病変肝臓組織に対応する放射状対称画像の表面トポグラフィ視覚化を示し、図６Ｃは図６Ｂと比較して異なる放射状対称フィルタパラメータが適用された病変肝臓組織に対応する放射状対称画像の表面トポグラフィ視覚化を示し、図６Ｄは病変組織に対応する強化超音波画像の表面トポグラフィの視覚化を示し、図６Ｅは図６Ｄと比較して異なるフィルタパラメータが使用された病気組織に対応する強化超音波画像の表面トポグラフィ視覚化を示す。
【図７】本解決策を実施することができる例示的な超音波装置の概略図である。
【図８】コンピューティングデバイスの例示的なアーキテクチャの概略図である。
【発明を実施するための形態】
【００３３】
本明細書に一般的に記載され、添付の図面に示される実施形態の構成要素は、多種多様な異なる構成で配置および設計され得ることが容易に理解されるであろう。したがって、図に表されているさまざまな実施形態の以下のより詳細な説明は、本開示の範囲を限定するものではなく、さまざまな実施形態の単なる代表例である。実施形態の様々な態様が図面に提示されているが、特に示されない限り、図面は必ずしも縮尺通りに描かれていない。
【００３４】
本解決策は、腎臓および肝臓に特に焦点を合わせた癌診断および治療監視という観点における超音波画像処理および次世代超音波画像化プラットフォームの設計のためのシステムおよび方法に関する。本解決策の超音波画像処理は、生の超音波信号を使用して、フィルタリングされた超音波データから画像特徴を抽出する。以前のプラットフォームは画像データを使用するが、これらの以前のプラットフォームには多くの欠点がある。たとえば、これらの以前のプラットフォームでは、腎臓および肝臓がんの検出に対する感度と特異性が低く、がんの管理と治療に大きな問題を引き起こすため、問題が発生する。一般に、肝転移の画像の外観は非特異的であり、組織診断には生検標本が必要である。本解決策は、本明細書で説明されるような新しい計算アルゴリズムを超音波装置に実装することによりこれらの欠点を解決する。
【００３５】
脂肪肝疾患を示す非健康肝臓の生検（組織病理学）画像が図３に示されている。矢印で強調表示されている生検画像の円形パターンは、肝臓に脂肪肝疾患があることを示している。しかしながら、円形パターンは、超音波画像を調べることでは表示または検出できない。したがって、超音波画像に基づいた病変診断は分かりにくい。超音波データを強化し、診断を行うのに有用なより多くの情報を提供できる画像処理装置が必要である。このような方法は、本解決策によって実施される。
【００３６】
いくつかのシナリオでは、本解決策は、脂肪肝、肝臓癌、および腎臓癌の早期検出に使用される。生検は現在、このような早期発見の標準的な手法である。本解決策は、生検を実施することによるサンプルを取得する必要はない。したがって、本解決策は、早期の検出を可能にし、患者に過度のストレスをかけることなく検出を可能にする。
【００３７】
本発明は、肝臓および／または腎臓の超音波データの強化をもたらすフレームワークを作成するために、異なるステップを組み合わせる。各ステップの最適化も主要な貢献の１つである。アルゴリズムの主要部分は、肝臓および腎臓の超音波データからの局所位相画像特徴の抽出に基づいている。さまざまな医療画像データを処理するために、位相ベースの画像強化と処理が以前から提案されている。画像位相情報は、大きさ情報よりも画像の視覚的外観に大きく寄与することが長い間知られているシーンの解釈における重要な要素である。位相特徴は強度不変であり、特に超音波画像の処理に重要な特性であるノイズに対してより確固たるものである。位相情報の抽出は、Ｂモード超音波画像（元の画像）がフーリエ変換操作によって周波数領域に変換され、バンドパス直交フィルタで乗算される周波数領域で実行される。
【００３８】
図１を参照すると、本解決策を理解するのに有用な画像強調のための例示的な方法１００の図が提供されている。ステップ１０２〜１０８の動作の組み合わせは、肝臓および／または腎臓超音波データの強化をもたらすフレームワークを作成する。
【００３９】
プロセスの主要部分は、超音波信号からの局所位相画像特徴の抽出に基づいている。画像位相情報は、大きさ情報よりも画像の視覚的外観に大きく寄与するため、シーンの解釈における重要な要素である。位相特徴は強度不変であり、特に超音波画像の処理に重要な特性であるノイズに対してより確固たるものである。
【００４０】
位相情報の抽出は、矢印１０１で示すように、生の超音波信号またはＢモード超音波画像が超音波装置から受信される周波数領域で実行される。次に、生の超音波画像信号は、１０２に示すように、信号領域から周波数領域に変換される。この変換を実現するためにフーリエ変換が使用される。その後、１０４で示されるように、フィルタリングされた超音波画像データは、バンドパス直交フィルタで乗算される。フーリエ変換およびバンドパス直交フィルタは、当技術分野では周知である。本明細書では、任意の既知のまたは既知のフーリエ変換および／またはバンドパス直交フィルタを制限なく使用することができる。たとえば、１０４では最適化されたバンドパス直交フィルタが採用されている。最適化はフィルタパラメータの選択に関連している。通常、これは試行錯誤によって実行される。適切なパラメータのセットが見つかると、これらは一定に保たれ、すべての超音波画像はこの最適化されたフィルタでフィルタリングされる。しかしながら、本明細書では自動フィルタパラメータ最適化が採用され得る。自動フィルタパラメータ最適化の技術は、当技術分野で知られている。本明細書では、自動フィルタパラメータ最適化のための任意の既知のまたは既知の技術を制限なく使用することができる。
【００４１】
いくつかのシナリオでは、バンドパス直交フィルタはアルファスケールフィルタを含む。そのため、１０４ではアルファスケールフィルタリングを採用できる。アルファスケールフィルタリングは、フィルタリングされた生の超音波信号から位相およびエネルギー特徴を抽出するために実行される。この特徴抽出により、医師が超音波データのみに基づいて肝臓がんまたは腎臓がんの診断を行うことができる視覚的明瞭度を有する軟組織パターンを示す強化超音波画像の生成が容易になる。特に、他の周波数領域直交フィルタ（ログガボールフィルタなど）などのフィルタリングは、軟部組織の特徴を強化しないため肝臓や腎臓の診断アプリケーションでの単独使用には適していない。
【００４２】
ログガボールフィルタは、アルファスケールフィルタに比べて計算量が少なく、リソース集約的であることを理解されたい。したがって、当業者は、超音波画像データを処理するときにアルファスケールフィルタを使用する動機付けを受けないであろう。しかしながら、発明者らは、アルファスケールフィルタとプロセス１０２の他の２つのフーリエ変換との組み合わせ、およびプロセス１０６−１０８の透過マップフィルタリングは、特定の軟組織特徴の視覚化に関心があるとき、特定の非自明な利点を提供することを認識した。これらの利点には、医師が超音波データのみに基づいて肝臓癌または腎臓癌の診断を行うことができる視覚的明瞭度を備えた軟組織特徴の強化が含まれる。
【００４３】
いくつかのシナリオでは、次の式（１）で定義されるように、周波数領域のアルファスケールフィルタが構築される。
【００４４】
【数１２】
2020508174.tif 000014
【００４５】
αはフィルタがＤＣ条件を満たすためにα＝０．２になるように選択された定数導関数パラメータである。σはフィルタアルファスケールパラメータ（たとえば、いくつかの実施形態では、フィルタアルファスケールパラメータは２５である）で、ｎｃは、次の式（２）を使用してフィルタα値から計算された単位正規化定数である。
【００４６】
【数１３】
2020508174.tif 000015
【００４７】
ｓはスケールパラメータで、αは微分パラメータである。いくつかの実施形態では、スケールパラメータは２であり、微分パラメータは１．８３である。
【００４８】
構築されたアルファスケールフィルタでフーリエ変換された超音波画像をフィルタリングし、逆フーリエ変換操作を使用して、画像の位相とエネルギー特徴を超音波画像から抽出する。局所エネルギー画像は、肝臓または腎臓の超音波画像の基礎となる構造情報をエンコードする。
【００４９】
アルファスケールフィルタリングの完了後、方法１００は１０５に進み、超音波画像データが逆フーリエ変換プロセスを介して周波数領域から信号領域に変換される。次に、１０６で１０５の結果を使用して、透過マップ推定を取得する。組織内の超音波信号の相互作用は、２つの主要なカテゴリー（つまり、散乱と減衰）に特徴付けることができる。後方散乱された超音波信号の情報（組織界面から超音波トランスデューサまで）は、これら２つの相互作用によって変調されるため、構造情報コーディングのメカニズムと見なすことができる。これに基づいて、モデルが開発された。このモデルは、本明細書では超音波信号送信マップと呼ばれる。超音波信号伝送マップは、超音波画像から適切な肝臓および／または腎臓組織構造を取り戻すのに役立つ。これを実現するために、組織内の散乱効果と減衰効果を組み合わせた線形補間モデルが採用されている。線形補間モデルは、次の式（３）で定義される。
【００５０】
【数１４】
2020508174.tif 000016
【００５１】
ＵＳ（ｘ、ｙ）は１０４で計算された局所エネルギー画像、ＵＳＡ（ｘ、ｙ）は信号伝送マップ、ＵＳＥ（ｘ、ｙ）は強化された肝臓／腎臓超音波画像、αは、局所領域の組織のエコー源性を表す定数値である。現在の解決策では、３つの異なる強化結果を得るために、αの３つの異なる値が使用されている。
【００５２】
１０６および１０８で示されるように、１０４で生成されたローカルエネルギー画像は、後続のフィルタリングプロセスで使用される。１０６と１０８のフィルタリングプロセスは、通常、抽出された位相とエネルギー特性を強化するために実行され、それにより、それに関連する軟組織の部分は、強化超音波画像においてより暗くなったり、より明るくなったりする（例えば、図３に示されるように）。
【００５３】
１０６では、減衰関数をイメージング深度の関数としてモデル化するよく知られたランベルト・ベールの法則を使用して信号伝送マップＵＳＡ（ｘ、ｙ）を推定することにより、ＵＳＥ（ｘ、ｙ）が抽出される。ランベルト・ベールの法則は、次の式（４）で定義される。
【００５４】
【数１５】
2020508174.tif 000017
【００５５】
ＵＳｏ（ｘ、ｙ）は初期強度画像（アルファスケールフィルタリングステップから取得したフィルタリングされた画像）、ＵＳＴ（ｘ、ｙ）は減衰強度画像、ηは減衰係数（例えば、いくつかの実施形態では、使用される減衰係数は２である）、ｄ（ｘ、ｙ）は超音波トランスデューサ表面からの距離である。
【００５６】
ＵＳＴ（ｘ、ｙ）が得られると、式（５）で定義される目的関数を最小化することによりＵＳＡ（ｘ、ｙ）が得られる。
【００５７】
【数１６】
2020508174.tif 000018
【００５８】
ωはインデックスセット、ｏは要素ごとの乗算、＊は畳み込み演算である。Ｄｊは、８つのキルシュフィルタと１つのラプラシアンフィルタで構成される高次微分フィルタのバンクを使用して取得される。ＷＪは、次の式（６）を使用して計算された重み行列である。
【００５９】
【数１７】
2020508174.tif 000019
【００６０】
ＵＳＡ（ｘ、ｙ）が推定されると、ＵＳＥ（ｘ、ｙ）は次の式（７）を使用して計算される。
【００６１】
【数１８】
2020508174.tif 000020
【００６２】
図３は、脂肪肝疾患を示す非健康な肝臓の生検（組織病理学）画像を示す。生検画像の円形パターンは、矢印で強調されているように、肝臓に脂肪肝疾患があることを示している。しかしながら、円形パターンは、超音波画像を調査することでは見えず、または検出できない（図４Ａおよび４Ｃ）。円形パターンは、Ｂモード超音波画像の図４Ｃから生成された図４Ｄに示されるように、本解決策を通じて生成された強化超音波画像で見ることができる。したがって、本解決策は、医師が超音波データのみに基づいて肝臓または腎臓疾患の診断を行い、生検標本を取得する必要性を回避するための方法を提供する。
【００６３】
ここで図２を参照すると、強化超音波画像に基づいて放射状対称画像を生成する方法が提供されている。このプロセスは、２０１で強化超音波画像の放射状対称フィルタリングを実行することから始まる。強調された超音波画像は、局所的な放射対称性を利用して、強調された画像内の関心領域を識別することによってさらに処理される。分析は、病変肝臓が健康な肝臓と比較して高度な放射状対称性を持っているという事実に依存している。解析は、高速放射状特徴検出アルゴリズムを使用して、強化された画像内の暗い球形を検索することにより実行される。２０３で、プロセスは、図５〜６に示されるように、放射状対称画像の生成を含む。図５Ｂは、正常／健康な肝臓組織に対応する放射状対称画像の表面トポグラフィ視覚化を示し、図５Ｃは、図５Ｂと比較して異なる放射状対称フィルタパラメータが適用された、正常／健康な肝臓組織に対応する放射状対称画像の表面トポグラフィ視覚化を示す。図６Ｂは、病変肝組織に対応する放射状対称画像の表面トポグラフィ視覚化を示し、図６Ｃは、図６Ｂと比較して異なる放射状対称フィルタパラメータが適用された、病変肝組織に対応する放射状対称画像の表面トポグラフィ視覚化を示す。２０５で、プロセスは、放射対称性画像に基づいて表面トポグラフィマップを出力することを任意で含む。高速放射対称アルゴリズムについては、以下で詳しく説明する。
【００６４】
マルチメディア画像の処理には、高速放射対称性が使用された。各半径ｎに対して、アルゴリズムは画像の勾配を使用して、プラスとマイナスの両方の影響を受けるピクセルに投じられる。これらのピクセルは、以下の式（８）および（９）を使用して計算される。
【００６５】
【数１９】
2020508174.tif 000021
【００６６】
【数２０】
2020508174.tif 000022
【００６７】
上記の式では、「ラウンド」は各ベクトル要素を最も近い整数に四捨五入し、「ｇ」は画像の勾配（たとえば、強化画像の勾配）で、「ｎ」は、画像内で検索された球面構造の半径値を表わす。ｐ−ｖｅとｐ＋ｖｅは、それぞれ、中心に向かって、または中心から遠ざかる方向に勾配ｇ（ｐ）を持つピクセルに対応する。これらのピクセルを使用して、それぞれＯｎおよびＭｎとして示される方向および大きさの投影画像が計算される。影響を受けるピクセルのそれぞれについて、ＯｎおよびＭｎの対応するポイントｐ＋ｖｅが１と｜｜ｇ（ｐ）｜｜それぞれ増加する。同様に、負の影響を受けるピクセルの場合、式（１０）〜（１３）で定義されているように、対応するポイントは各画像で同じ量だけ減少する。
【００６８】
【数２１】
2020508174.tif 000023
【００６９】
これらの画像を使用して、放射状対称応答画像は次のように定義される。
【００７０】
【数２２】
2020508174.tif 000024
【００７１】
ここで、「＊」は畳み込み演算を示す。Ａｎは等方性ガウス関数であり、Ｆｎは式（１４）および（１５）で定義される。
【００７２】
【数２３】
2020508174.tif 000025
【００７３】
【数２４】
2020508174.tif 000026
【００７４】
αは放射状の厳密性パラメータであり、ｋｎは異なる半径にわたるスケーリング係数である。最終的な完全な放射状対称変換は、さまざまな半径値に対してこの操作を実行し、結果の機能マップを式（１６）のように合計することによって定義される。
【００７５】
【数２５】
2020508174.tif 000027
【００７６】
最終的なトポグラフィ画像は、放射状対称および位相画像上の画像強度値に基づいてトポグラフィマップを表示することによって取得される。病変肝臓は、健康な肝臓と比較して、より高いトポグラフィを持っている。表面トポグラフィ画像の例は、図５Ｄ−Ｅに健康な肝臓、図６Ｄ−Ｅに病変肝臓が示されている。具体的には、図５Ｄは、正常／健康組織に対応する強化超音波画像の表面トポグラフィ視覚化を示し、図５Ｅは、図５Ｄと比較して異なるフィルタパラメータが使用された、正常／健康組織に対応する強化超音波画像の表面トポグラフィ視覚化を示す。図６Ｄは病変組織に対応する強化超音波画像の表面トポグラフィの視覚化を示し、図６Ｅは図６Ｄと比較して異なるフィルタパラメータが使用されている病変組織に対応する強化超音波画像の表面トポグラフィ視覚化を示す。
【００７７】
ここで図７を参照すると、本解決策を実施することができる例示的な超音波装置７００の概略図が提供されている。超音波装置７００は、計算装置７０４が配置される装置スタンド７０２を含む。コンピューティングデバイス７０４は、一般に、超音波装置７００の動作を制御するように構成される。そのような制御は、入力デバイス７１６を介した、および／または所定の規則に従って、ユーザとソフトウェアの相互作用に応答し得る。ロボットアーム７０６は、デバイススタンド７０２に移動可能に取り付けられている。ロボットアーム７０６は、複数の関節７０８を備えた関節アームを含むが、これに限定されない。超音波トランスデューサ７１０は、ロボットアーム７０６の遠位端に配置される。超音波トランスデューサ７１０は、検査されるオブジェクトの表面７１２上で動かされ、生の超音波信号を生成する。オブジェクトは、患者位置決めテーブル７１４上に配置される。生の超音波信号は、本解決策に従って処理するためにローカルコンピューティングデバイス７０４またはリモートコンピューティングデバイスに提供される。この処理の結果は、生の超音波データの強化超音波画像（例えば、図４に示されているものなど）への変換である。
【００７８】
ここで図８を参照すると、コンピューティングデバイス８００の例示的なアーキテクチャの図が提供されている。コンピューティングデバイス８００は、超音波デバイス（例えば、図７の超音波装置７００）に対してローカルであり得る。この場合、コンピューティングデバイス７０４は、コンピューティングデバイス８００と同一または実質的に類似している。他のシナリオでは、コンピューティングデバイス８００は、超音波装置から離れて配置されている。この場合、超音波装置のコンピューティングデバイス７０４は、ネットワーク（例えば、イントラネットまたはインターネット）を介した２つのコンピューティングデバイス間の通信を促進するネットワークインターフェースを含む。
【００７９】
コンピューティングデバイス８００は、図８に示されるものよりも多いまたは少ないコンポーネントを含むことができる。しかしながら、示されているコンポーネントは、本解決策を実施する例示的な実施形態を開示するのに十分である。図８のハードウェアアーキテクチャは、生の超音波信号に基づいて強化超音波画像の生成を容易にするように構成された代表的なコンピューティングデバイスの一実施形態を表す。そのようなものとして、図８のコンピューティングデバイス８００は、本解決策に従ってそのような強化超音波画像を提供する方法の少なくとも一部を実施する。
【００８０】
コンピューティングデバイス８００のいくつかまたはすべてのコンポーネントは、ハードウェア、ソフトウェア、および／またはハードウェアとソフトウェアの組み合わせとして実装することができる。ハードウェアには、限定されないが、１つ以上の電子回路が含まれる。電子回路は、限定されないが、受動部品（たとえば、抵抗器およびキャパシタ）および／または能動部品（たとえば、増幅器および／またはマイクロプロセッサ）を含むことができる。受動的および／または能動的構成要素は、本明細書に記載の方法論、手順、または機能のうちの１つ以上を実行するように適合、配置、および／またはプログラムすることができる。
【００８１】
図８に示すように、コンピューティングデバイス８００は、ユーザインターフェース８０２、中央処理装置（「ＣＰＵ」）８０６、システムバス８１０、システムバス８１０を介してコンピューティングデバイス８００の他の部分に接続されアクセス可能なメモリ８１２、およびシステムバス８１０に接続されたハードウェアエンティティ８１４を備える。ユーザインターフェースは、入力デバイス（例えば、キーパッド８５０、カーソル制御デバイス８５８、および／またはカメラ８６０）および出力デバイス（例えば、スピーカー８５２、ディスプレイ８５４、および／または発光ダイオード８５６）を含むことができ、コンピューティングデバイス８００の動作を制御するためのユーザとソフトウェアの相互作用を容易にする。
【００８２】
ハードウェアエンティティ８１４の少なくともいくつかは、ランダムアクセスメモリ（「ＲＡＭ」）、ディスクドライバ、および／またはコンパクトディスク読み取り専用メモリ（「ＣＤ−ＲＯＭ」）であり得るメモリ８１２へのアクセスおよび使用を含むアクションを実行する。ハードウェアエンティティ８１４は、本明細書に記載の方法論、手順、または機能のうちの１つ以上を実行するように構成された命令８２０（例えば、ソフトウェアコード）の１つ以上のセットが記憶されるコンピュータ可読記憶媒体８１８を備えるディスクドライブユニット８１６を含むことができる。命令８２０は、コンピューティングデバイス８００による実行中に、メモリ８１２内および／またはＣＰＵ８０６内に完全にまたは少なくとも部分的に存在することもできる。メモリ８１２およびＣＰＵ８０６はまた、機械可読媒体を構成できる。本明細書で使用される「機械可読媒体」という用語は、命令８２０の１つ以上のセットを格納する単一の媒体または複数の媒体（例えば、集中型または分散型データベース、および／または関連キャッシュおよびＲＰＵ）を指す。本明細書で使用される「機械可読媒体」という用語は、コンピューティングデバイス８００による実行のための命令セット８２０を格納、エンコード、または搬送することができ、これにより、コンピューティングデバイス８００は、本開示の方法論のいずれか１つまたは複数を実行する。
【００８３】
いくつかのシナリオでは、ハードウェアエンティティ８１４は、強化超音波画像の提供を容易にするためにプログラムされた電子回路（例えば、プロセッサ）を含む。これに関して、電子回路は、コンピューティングデバイス８００にインストールされたソフトウェアアプリケーション８６４にアクセスして実行できることを理解されたい。ソフトウェアアプリケーション８２４は、一般に、強化超音波画像の生成を促進するように動作する。議論が進むにつれて、ソフトウェアアプリケーション８６４のその他の機能が明らかになるだろう。
【００８４】
本明細書で説明される本解決策を実装するソフトウェアアプリケーション８６４は、コンピュータ可読記憶媒体にソフトウェアプログラムとして格納され、ＣＰＵ８０６上で実行するように構成される。さらに、本解決策のソフトウェア実装は、限定されないが、分散処理、コンポーネント／オブジェクト分散処理、並列処理、仮想マシン処理を含むことができる。様々なシナリオにおいて、ネットワーク環境に接続されたネットワークインターフェースデバイス８６２は、命令８２０を使用してネットワーク上で通信する。
【００８５】
本解決策は、その精神または本質的な特徴から逸脱することなく、他の特定の形態で具現化され得る。説明された実施形態は、すべての点で例示的であり、限定的ではないと見なされるべきである。したがって、本発明の範囲は、この詳細な説明より添付の特許請求の範囲によって示される。特許請求の範囲と同等の意味および範囲内にあるすべての変更は、その範囲内に含まれるものとする。
【００８６】
本明細書全体を通して、特徴、利点、または類似の用語への言及は、本発明で実現され得る特徴および利点のすべてが本発明の任意の単一の実施形態であるべきであることを意味しない。むしろ、特徴および利点に言及する言葉は、実施形態に関連して説明される特定の特徴、利点、または特性が本発明の少なくとも１つの実施形態に含まれることを意味すると理解される。したがって、本明細書全体にわたる特徴および利点、ならびに類似の言語の議論は、必ずしもそうとは限らないが、同じ実施形態を指すことがあり得る。
【００８７】
さらに、説明された本発明の特徴、利点、および特性は、一つ以上の実施形態において任意の適切な方法で組み合わせることができる。当業者は、本明細書の説明に照らして、特定の実施形態の一つ以上の特定の特徴または利点なしで本発明を実施できることを認識するであろう。他の例では、本発明のすべての実施形態に存在するわけではない特定の実施形態で、追加の特徴および利点が認識される場合がある。
【００８８】
本明細書を通して「一実施形態」、「実施形態」、または類似の用語への言及は、示された実施形態に関連して説明される特定の特徴、構造、または特性が本発明の少なくとも一つの実施形態に含まれることを意味する。したがって、本明細書全体にわたる「一実施形態では」、「実施形態では」という語句、および類似の文言は、必ずしもそうではないが、すべて同じ実施形態を指す場合がある。
【００８９】
本明細書で使用される場合、単数形「ａ」、「ａｎ」、および「ｔｈｅ」は、文脈がそうでないことを明確に指示しない限り、複数の参照を含む。別に定義されない限り、本明細書で使用されるすべての技術用語および科学用語は、当業者によって一般に理解されるのと同じ意味を有する。本明細書で使用されている「含む」という用語は、「含むが、これに限定されない」ことを意味する。
【００９０】
本明細書に開示および特許請求される装置、方法、およびアルゴリズムのすべては、本開示に照らして過度の実験を行うことなく作成および実行することができる。本発明を好ましい実施形態に関して説明してきたが、概念、精神、および本発明の範囲から逸脱することなく、装置、方法、および方法のステップの順序に変更を適用できることは当業者には明らかであろう。より具体的には、本明細書に記載の構成要素に特定の構成要素を追加、組み合わせ、または置換しても、同じまたは同様の結果が得られることは明らかであろう。当業者に明らかなそのような類似の置換および修正はすべて、定義された本発明の精神、範囲、および概念内にあるとみなされる。
【００９１】
上記で開示された特徴および機能、ならびに代替は、多くの他の異なるシステムまたは用途に組み合わされてもよい。様々な現在予測できないまたは予期しない代替、修正、変形、または改善が当業者によって行われ得、それらのそれぞれも開示された実施形態に含まれることが意図される。
【図１】
2020508174.tif 000050
【図２】
2020508174.tif 000051
【図３】
2020508174.tif 000052
【図４】
2020508174.tif 000053
【図５】
2020508174.tif 000054
【図６】
2020508174.tif 000055
【図７】
2020508174.tif 000056
【図８】
2020508174.tif 000057
【国際調査報告】
2020508174.tif 000058