(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2024-03-07
(54)【発明の名称】媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法およびシステム
(51)【国際特許分類】
A61B 8/14 20060101AFI20240229BHJP
【FI】
A61B8/14
【審査請求】有
【予備審査請求】未請求
(21)【出願番号】P 2023556951
(86)(22)【出願日】2022-03-15
(85)【翻訳文提出日】2023-11-14
(86)【国際出願番号】 EP2022056742
(87)【国際公開番号】W WO2022194888
(87)【国際公開日】2022-09-22
(32)【優先日】2021-03-19
(33)【優先権主張国・地域又は機関】EP
(81)【指定国・地域】
(71)【出願人】
【識別番号】508291928
【氏名又は名称】スーパー ソニック イマジン
【氏名又は名称原語表記】SUPER SONIC IMAGINE
(74)【代理人】
【識別番号】100145403
【氏名又は名称】山尾 憲人
(74)【代理人】
【識別番号】100135703
【氏名又は名称】岡部 英隆
(74)【代理人】
【識別番号】100227927
【氏名又は名称】中村 拓
(72)【発明者】
【氏名】ジャン,ボゥ
【テーマコード(参考)】
4C601
【Fターム(参考)】
4C601BB03
4C601DD19
4C601DD23
4C601EE11
4C601FF03
4C601GA20
4C601GA21
4C601JC05
4C601JC26
4C601LL33
(57)【要約】
本発明は、媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法に関し、上記方法は、処理システムによって実施されてもよく、本方法は、媒体の複数の体積単位のうちの少なくとも1つについて、モダリティおよび/またはソースのデータに基づいて光学特性が決定される光学決定ステップと、体積単位のうちの少なくとも1つについて、第2のモダリティおよび/またはソースのデータに基づいて蛍光特性が決定される蛍光決定ステップと、体積単位の決定された光学特性および蛍光特性に基づいて媒体のデータ表現がレンダリングされるレンダリングステップと、を含む。本発明はまた、対応する処理システムに関する。
【特許請求の範囲】
【請求項1】
媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法であって、
前記媒体の複数の体積単位のうちの少なくとも1つについて、第1のモダリティおよび/またはソースの前記データに基づいて、光学特性が決定される光学決定ステップと、
前記体積単位の少なくとも1つについて、第2のモダリティおよび/またはソースの前記データに基づいて、蛍光特性が決定される蛍光決定ステップと、
前記媒体のデータ表現が、前記体積単位の前記決定された光学特性および蛍光特性に基づいてレンダリングされるレンダリングステップと、
を含む方法。
【請求項2】
体積測定の前記媒体の前記データ表現が、視覚的表現および/または2Dもしくは3D画像を含む、
請求項1に記載の方法。
【請求項3】
前記光学特性が、少なくとも体積単位の光吸収率を含み、ならびに/あるいは、
前記光学特性が、前記少なくとも1つの体積単位の前記第1のモダリティおよび/またはソースの前記データに適用される、第1の媒体-光の相互作用マッピングに従って決定され、
任意選択的に、前記第1の媒体-光の相互作用マッピングが、前記第1のモダリティおよび/またはソースの前記データの関数として反射率、指向性、および吸収率のうちの少なくとも1つを決定する、
請求項1または2に記載の方法。
【請求項4】
前記蛍光特性が、前記少なくとも1つの体積単位の前記第2のモダリティおよび/またはソースの前記データに適用される第1の発光および/または吸収マッピングに従って決定される、
請求項1から3のいずれか一項に記載の方法。
【請求項5】
前記光学決定ステップにおいて、第3のモダリティおよび/またはソースの前記データに基づいて、かつ第1の媒体-光の相互作用マッピングとは異なる第2の媒体-光の相互作用マッピングに従って、前記体積単位のうちの少なくとも1つについて、前記光学特性が決定され、ならびに/あるいは、
前記蛍光決定ステップにおいて、第4のモダリティおよび/またはソースの前記データに基づいて、かつ第1の発光および/または吸収マッピングとは異なる第2の発光および/または吸収マッピングに従って、前記体積単位の少なくとも1つについて、蛍光特性が決定される、
請求項1から4のいずれか一項に記載の方法。
【請求項6】
前記蛍光決定ステップおよび前記光学決定ステップが、同時におよび/または並行して実行される、
請求項1から5のいずれか一項に記載の方法。
【請求項7】
前記決定された光学特性および蛍光特性がデータストレージに記憶され、前記レンダリングステップがその後、前記記憶された特性に基づいて1回または数回実行されるか、または
前記レンダリングステップが、前記蛍光決定ステップおよび/または前記光学決定ステップと同時に行われる、
請求項1から6のいずれか一項に記載の方法。
【請求項8】
前記第1のモダリティおよび/またはソースの前記データが、Bモード超音波撮像からの反射率情報を含み、ならびに/あるいは、
前記第2のモダリティおよび/またはソースの前記データが、せん断波エラストグラフィ超音波撮像からの組織弾性情報を含む、
請求項1から7のいずれか一項に記載の方法。
【請求項9】
前記第2のモダリティおよび/またはソースの前記データが、前記媒体に対して配置された介入デバイスの位置および/または向きを示す、介入デバイス情報、ならびに/あるいは、前記介入デバイスの形状の事前設定された情報を含み、
少なくとも1つの体積単位の前記蛍光決定ステップにおいて、前記介入デバイスによって占有されている体積単位を識別することによる前記介入デバイス情報に基づいて、蛍光特性が決定される、
請求項1から8のいずれか一項に記載の方法。
【請求項10】
前記レンダリングステップにおいて、前記データ表現の少なくとも一部が、前記媒体内に配置された介入デバイスの位置および/または向きを示す介入デバイス情報と、前記介入デバイスの形状の事前設定された情報と、に基づいてレンダリングされる、請求項1から8のいずれか一項に記載の方法。
【請求項11】
前記レンダリングステップが、前記データ表現が光線追跡体積レンダリング方法に従ってレンダリングされる光線追跡レンダリングステップを含む、
請求項1から10のいずれか一項に記載の方法。
【請求項12】
前記光線追跡レンダリングステップが、
前記媒体に対して事前定義された幾何学的関係で仮想光源を配置するサブステップと、
被写界深度および/または開口形状を定義するサブステップと、
前記媒体に対して事前定義された幾何学的関係で仮想ビュー平面を配置するサブステップであって、前記データ表現が前記仮想ビュー平面の関数としてレンダリングされる、サブステップと、
のうちの少なくとも1つのサブステップを含む、請求項11に記載の方法。
【請求項13】
前記マルチモダリティおよび/またはマルチソースデータが、前記媒体の画像データおよび/または3Dデータを含み、ならびに/あるいは、体積単位が、少なくとも1つのボクセルを含む、
請求項1から12のいずれか一項に記載の方法。
【請求項14】
前記光学決定ステップ、前記蛍光決定ステップ、および前記レンダリングステップのうちの少なくとも1つが、少なくとも1つの第1の人工知能ベースのアルゴリズムおよび/または機械学習アルゴリズムを実施または使用する、
請求項1から13のいずれか一項に記載の方法。
【請求項15】
少なくとも1つの第2の人工知能ベースのアルゴリズムおよび/または予め訓練された機械学習アルゴリズムが、前記レンダリングされたデータ表現の関数として事前定義されたタスクを実行し、前記タスクが、任意選択的に、回帰タスク、分類タスク、およびセグメント化タスクのうちの少なくとも1つを含む、請求項1から14のいずれか一項に記載の方法。
【請求項16】
データ処理システムによって実行された場合に、前記データ処理システムに、請求項1から15のいずれか一項に記載の方法を実行させる、コンピュータ可読命令を含むコンピュータプログラム。
【請求項17】
媒体のマルチモダリティおよび/またはマルチソースデータを処理するためのシステムであって、
前記媒体の複数の体積単位のうちの少なくとも1つについて、第1のモダリティおよび/またはソースの前記データに基づいて、光学特性を決定し、
前記体積単位のうちの少なくとも1つについて、第2のモダリティおよび/またはソースの前記データに基づいて、蛍光特性を決定し、
前記体積単位の前記決定された光学特性および蛍光特性に基づいて、前記媒体のデータ表現をレンダリングする
ように構成された、処理ユニットを含む、
システム。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法およびシステムに関する。特に、本開示は、特に医療撮像のための、本方法を実施する画像処理方法およびシステムに関する。
【背景技術】
【0002】
検査、特に例えば医療検査は、多くの場合、コンピュータ実装撮像方法、例えば超音波撮像によって支援される。この目的のために、検査された媒体(例えば、石、動物、人体またはその一部)からの検査データが取得され、それを検査ユーザまたは医師などの別のユーザに提示するために処理される。
【0003】
例えば、超音波撮像は、トランスデューサによって送信される1つまたは複数の超音波パルス(または波)による媒体の音波処理からなる。これらのパルスのエコーに応答して、例えば同じトランスデューサを使用して、超音波信号データが取得される。超音波撮像は、異なるモダリティ(またはモード)、例えば、Bモード(輝度モード)およびShearWave(登録商標)(SWE、せん断波エラストグラフィ)を有することができ、これらは同じ超音波撮像システムによって提供できる。
【0004】
そのような画像の解釈、またはより一般的には任意の検査データの解釈には、高度な専門知識が必要である。これは、特に、3次元画像データ、例えば3D超音波画像の場合である。
【0005】
検査データは、異なるモダリティおよび/または異なるソースのデータが、例えばユーザに示される同じ3D(三次元)表現に組み合わされる場合、さらに複雑になる可能性がある。そのような異なるモダリティは、異なるモードで同じシステムによって取得されたデータ、例えば、超音波システムによって取得されたBモード画像データおよびSWE画像データを含んでもよい。言い換えれば、異なるモダリティのデータは、同じソース(例えば、超音波)であるが異なるモダリティ(例えば、異なる超音波取得モード)のデータを含むことができる。
【0006】
異なるソースはまた、異なるシステム、例えば超音波システムおよび別の撮像システム、例えばコンピュータ断層撮影、磁気共鳴または陽電子放出断層撮影などの医療撮像システムによって得られたデータを含んでもよい。
【0007】
結果として、組み合わせられたマルチモダリティおよび/またはマルチソースデータ(例えば、レンダリングされた3D画像)の視覚的表現は、ユーザに過負荷な情報を提供する可能性があり、複雑な場合にあまり明瞭ではなく、潜在的に誤解の原因となり、その結果、より関連性の高い情報が他のあまり関連性のない情報によって隠される場合がある。
【0008】
画像、特に3D画像または深度情報を有する画像をレンダリングするための様々な既知の技術が存在する。例えば、医療画像の視覚化のために、以下の文献に記載されているように、様々な体積レンダリングアルゴリズムが進化している。
Zhang Q, Eagleson R, Peters TM. Volume visualization: a technical overview with a focus on medical applications. J Digit Imaging. 2011 Aug; 24(4): 640-64. doi: 10.1007/sl0278-010-9321-6. PMID: 20714917; PMCID: PMC3138940.
【0009】
特に、多断面再構成(MPR)は、任意に配置された直交平面または傾斜平面を使用して3D体積から2次元(2D)スライスを抽出する既知の画像処理技術である。MPRは、3D取得を伴う市販の超音波システムの標準的な解決策である。
【0010】
別の例は、オブジェクト表面が明示的または暗黙的にモデル化されることを必要とする表面レンダリング(SR)である。次いで、表面を陰影付けし、例えば光線追跡法によってレンダリングする。この技術は、ビデオゲームおよびデジタル映画で広く使用されている。しかし、表面モデリングは、事前に未知の数の未知の構造の自動セグメント化を必要とし、これは実際にはほとんど堅牢ではなく、追加の計算コストを必要とするため、医療用超音波の状況にはあまり適応していない。
【0011】
さらに別の例は、中間幾何学的表現を計算することなく3Dデータ全体を2D(2次元)画像として表示する体積レンダリング(VR)である。
【0012】
しかしながら、これらの技術は、上記の問題を解決または低減しない。
【0013】
さらに、医療検査(例えば、超音波撮像法を使用する)は、特に例えば超音波撮像によって導かれる医療介入のために、リアルタイムまたは準リアルタイムのモニタリングを必要とすることが多い。このため、既知のレンダリング方法の演算コストの増加とユーザの潜在的な情報過負荷の両方が、リアルタイム要件に関して不利になる可能性がある。
【発明の概要】
【0014】
現在、上述の問題を克服すること、特に、マルチモダリティおよび/またはマルチソースデータの視覚化、理解および解釈を容易にする、体積測定媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法およびシステムを提供することが依然として望ましい。さらに、方法およびシステムは、望ましくは、それぞれのデータ表現、例えば深度情報を有する3D画像または2D画像における上記マルチモダリティおよび/またはマルチソースデータの組み合わせ表現またはマージ表現を提供する。
【0015】
したがって、本開示の実施形態によれば、(例えば、体積測定)媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法が提供される。上記方法は、処理システムによって実施されてもよい。本方法は、以下のステップを含む:
・媒体の複数の体積単位のうちの少なくとも1つについて、第1のモダリティおよび/またはソースのデータに基づいて、光学特性が決定される光学決定ステップ、
・体積単位の少なくとも1つについて、第2のモダリティおよび/またはソースのデータに基づいて、蛍光特性が決定される蛍光決定ステップ、
・体積単位の決定された光学特性および蛍光特性に基づいて、体積測定媒体のデータ表現がレンダリングされるレンダリングステップ。
【0016】
そのような方法を提供することにより、異なるモダリティおよび/またはソースからのデータ(例えば、医療撮像データ)の視覚化、理解および/または解釈が容易になり、したがって、専門家ユーザ(例えば、医療専門家)と非専門家ユーザ(例えば、患者または訓練されていない専門家などの非専門家)の両方にとって、光学的に妥当な方法で、より直感的になる。
【0017】
本開示は、従来のレンダリングアルゴリズムを、多色蛍光機構のモデリングと組み合わせることによって適合させてもよい。このようにして、異なるモダリティおよび/またはソースを有するデータ体積単位を組み合わせ、統一されたレンダリングシーンにおいて解釈可能にできる。
【0018】
したがって、提案された技術は、複数の撮像モードおよび/または撮像ソース(例えば、第1のモダリティとしての3D Bモードおよび第2のモダリティとしての3D SWEモード)が、異なるカスタマイズ可能な蛍光色によってレンダリングされることを可能にし、その結果、それらが統一されたシーンで視覚化可能となる。この例では、Bモード撮像からの媒体のエコー輝度情報およびSWEからの組織弾性情報を、同時にレンダリングできる。
【0019】
言い換えれば、媒体特性は、同じ視覚化可能シーンで組み合わせることができる。これは、ユーザガイダンス、病理または治療のモニタリングのための視覚的な手がかりにつながる可能性があり、また、マルチモダリティおよび/またはマルチソースデータのより容易なユーザおよび患者の理解および解釈につながる可能性がある。
【0020】
さらなる結果として、結果として得られる統一されたシーンは、任意の機械学習または他のAI関連アプリケーションのための改善された基盤を提供できる。例えば、シーンは、CNN(畳み込みニューラルネットワーク)または事前定義されたサイズおよび/またはフォーマットの入力データを期待する別の機械学習アルゴリズムによって処理することができる、事前定義されたサイズを有する2Dまたは3D画像の形態であってもよい。結果として、2Dまたは3D画像は、さらなる事前処理を必要とせずに、CNN(または任意の他の機械学習アルゴリズム)のための単一の標準化された入力として使用できる。対照的に、元のマルチモダリティおよび/またはマルチソースデータは、それぞれの場合に特定の事前処理を必要とする、異なるまたは変化する解像度を有してもよい。
【0021】
さらに、光学特性および蛍光特性の決定に起因して、結果として得られるレンダリングされた表現は、元のマルチモダリティおよび/またはマルチソースデータと比較して、対照的に強化され得、より少ないノイズを含み得る。この状況は、人間がレンダリングされた表現を正しく解釈するのを助けるだけではない。それはまた、機械学習アルゴリズムによって実行される任意の分類および/または回帰タスクの結果を改善できる。
【0022】
例えば、レンダリングされたデータ表現は、元のデータのより現実的な重み付けがAIベースのアルゴリズムに提供されるようなものであってもよい。言い換えれば、(AIベースの)アルゴリズムは、より重み付けされた、すなわちより関連性の高いデータに対して、より敏感であり得る。一例では、SWEモダリティのデータは、蛍光決定ステップを通して増加した重み付けまたは注意(すなわち、強調表示および/または色付けできる)を取得できる。言い換えれば、上記データ(例えば、媒体内の病変をマーキングする)はまた、AIベースのアルゴリズムによって、注意を高めて処理されてもよい。
【0023】
本開示のステップはまた、1つまたは複数の機械学習アルゴリズムまたは任意の他のAI(人工知能)ベースのアルゴリズムによって、実行または少なくとも支援されてもよい。
【0024】
提案された方法のさらなる利点は、計算コストの可能な低減である。例えば、提案された方法では、オブジェクトまたは表面セグメント化なしでマルチモード3D画像の融合を可能にする。さらに、計算コストは、シーン内のオブジェクトの数とは無関係であることが有利であり、したがって、複数のオブジェクトに対して複数のレンダリングを行う必要はない。計算コストは、管理される体積単位の数にのみ依存し得る。さらに、蛍光特性の決定は、蛍光の物理的原理に依存し、任意の陰影付けには依存しない。これは、より物理的な直感を提供し、それ自体より自明であることによって、実用的なパラメータの最適化を容易にする。
【0025】
本方法は、媒体を複数の体積単位に分割するステップをさらに含み得る。上記体積単位は、いずれにせよ、例えば3次元マルチモダリティおよび/またはマルチソース画像データ内のボクセルの形態で、マルチモダリティおよび/またはマルチソースデータによって既に提供されてもよい。媒体の外側にも体積単位が存在することがさらに可能である。媒体の外側のこれらの体積単位のうちの少なくとも1つは、レンダリングされたデータ表現において考慮される別の手段、例えば介入デバイスによって占有されてもよい。例えば、介入デバイスの媒体および/または本体のサンプリングされた位置を表すボクセルは、座標系、例えばデカルトもしくは極座標系、または任意の他の事前定義された系に配置されてもよい。
【0026】
一例では、モダリティおよび/またはソースのデータは、媒体の3D画像情報を含んでもよい。上記情報は、デカルトもしくは極座標系または他の事前定義されたシステムなどの事前定義された座標系に従って、ボクセルモードによってボクセル上で取得されてもよい。そのような3D画像情報は、例えば、超音波マトリックスプローブ、すなわち、マトリックス形式で配置された複数のトランスデューサを有するか、またはいくつかのトランスデューサタイプをマージするプローブ、を使用して取得することができる。
【0027】
3D画像情報はまた、2D画像スライスの積み重ねから取得されてもよい。これは、機械的な1D線形超音波プローブから達成されてもよい。プローブは、所与の方向に機械的に回転させることができ、異なる角度で1Dプローブは2D画像スライスを取得する。
【0028】
光学特性が決定される少なくとも1つの体積単位は、蛍光特性が決定される体積単位と同じまたは異なる体積であってもよい。したがって、1つの体積単位は、決定された光学特性および/または決定された蛍光特性を有してもよく、またはその両方を有しなくてもよい。さらに、少なくとも1つの体積単位について、光学特性および蛍光特性のいずれも決定しないことが可能である。
【0029】
例示的な一実施形態では、いくつかの体積単位について光学特性が決定され、および/または、いくつかの体積単位について蛍光特性が決定される。別の例では、複数の体積単位のそれぞれについて光学特性が決定され、および/または複数の体積単位のそれぞれについて蛍光特性が決定される。
【0030】
レンダリングステップは、計算ステップまたは演算ステップとして理解されてもよい。それはまた、グラフィカルレンダリングステップおよび/または3Dレンダリングステップを含んでもよい。
【0031】
体積測定媒体のデータ表現は、複数の深さレベルを有する視覚的またはグラフィック表現を含んでもよい。
【0032】
一例では、データ表現、特に上記視覚的表現は、1つまたは複数の画像、例えば2次元(2D)または3次元(3D)画像を含んでもよく、またはその形態であってもよい。
【0033】
光学特性は、受動光学特性として理解することができる。
【0034】
蛍光特性は、(二次)光源を放射する特性として理解することができる。したがって、「光学特性」と「蛍光特性」との違いは、前者が受動的であり、後者が能動的であることであり得る。「受動的」は、それらが光源ではないことを意味することができ、「能動的」は、それらが光を放射することができることを意味することができる。
【0035】
光学決定ステップは、第1のモダリティおよび/またはソースのデータに対して実行されるだけでなく、異なるモダリティ(例えば、BモードおよびSWE)および/またはソース(例えば、超音波およびCT)のデータに基づいて、光学特性が決定されることも可能である。
【0036】
同様に、蛍光決定ステップにおいて、異なるモダリティおよび/またはソースのデータに基づいて蛍光特性が決定されることも可能である。例えば、各モダリティおよび/またはソースについて、異なるモダリティおよび/またはソースの区別を容易にするために、異なる蛍光色を決定することができる(例えば、SWEとして赤、認識された介入デバイス(針、マーカーまたは媒体中の関心領域など)のマーキングとして青)。
【0037】
一例では、光学特性は、体積単位の光吸収率を含んでもよく、または体積単位の光吸収率であってもよい。
【0038】
さらなる例では、光学特性は、少なくとも1つの体積単位の第1のモダリティおよび/またはソースのデータに適用される(例えば、事前定義された)媒体-光の相互作用マッピングに従って決定されてもよい。したがって、第1のモダリティおよび/またはソースのデータは、媒体-光の相互作用マッピングの入力として使用することができる。光学特性は、媒体-光の相互作用マッピングの出力であってもよい。
【0039】
媒体-光の相互作用マッピングは、第1のモダリティおよび/またはソースのデータの関数として反射率、指向性、および吸収率のうちの少なくとも1つを決定してもよい。
【0040】
このマッピングは、事前定義されたマッピング関数および/または事前定義されたマッピング規則であってもよく、あるいは、事前定義されたマッピング関数および/または事前定義されたマッピング規則を含んでもよい。例えば、マッピングは、光学特性を決定するためにどの種類のモダリティおよび/またはソースが使用されるかを(マッピング規則に基づいて)決定できる。例えば、光学特性は、Bモードデータ(この例では第1のモダリティおよび/またはソースである)に基づいて決定できる。さらなる例として、マッピングは、第1のモダリティおよび/またはソースのパラメータの値または値の範囲(例えば、Bモードデータ内の輝度値)の関数として、光学特性のパラメータの値または値の範囲(例えば、反射率、指向性、および吸収率)を、(マッピング規則および/または関数に基づいて)決定できる。
【0041】
蛍光特性は、少なくとも1つの体積単位の第2のモダリティおよび/またはソースのデータに適用される(例えば、事前定義された)発光および/または吸収マッピングに従って決定できる。上記マッピングは、例えば、発光マッピング、または発光および吸収マッピングであってもよい。したがって、上記マッピングは、少なくとも1つの体積単位の発光率を定義できる。
【0042】
したがって、この発光および/または吸収マッピングもまた、事前定義されたマッピング関数および/または事前定義されたマッピング規則であってもよく、またはそれを含んでもよい。例えば、マッピングは、どの種類のモダリティおよび/またはソースが蛍光特性を決定するために使用されるかを(マッピング規則に基づいて)決定できる。例えば、蛍光特性は、SWEデータ(この例では第2のモダリティおよび/またはソースである)に基づいて決定できる。さらなる例として、マッピングは、第2のモダリティおよび/またはソースのパラメータの値または値の範囲(例えば、SWEデータ内の輝度/色の値)の関数として、蛍光特性のパラメータの値または値の範囲(例えば、発光速度、光吸収率および/または色)を(マッピング規則および/または関数に基づいて)決定できる。
【0043】
別の例によれば、本開示によるマッピングは、所与のボクセル値の光学特性および/または蛍光特性をマッピングしてもよい。マッピングがマップを決定することも可能である。上記マップは、ソースおよび/またはモダリティデータのデータ構造に対応するデータ構造を提供してもよい。例えば、マップは、ソースおよび/またはモダリティデータを形成する画像と同じサイズであってもよい。したがって、マップは、各ピクセル/ボクセルにおける光学/蛍光特性を記憶してもよい。
【0044】
蛍光特性は、体積単位によって形成された二次光源を定義することができる。言い換えれば、体積単位は、蛍光特性に応じた2次光を放射してもよい。したがって、上記二次光は、第1のスペクトルに吸収され、体積単位によって第2のスペクトルの二次光として(再)放射される一次光に基づいて、決定されてもよい。
【0045】
体積単位の少なくとも1つの光学決定ステップでは、第1の媒体-光の相互作用マッピングとは異なる第2の媒体-光の相互作用マッピングに従って、第3のモダリティおよび/またはソースのデータに基づいて、光学特性を決定できる。
【0046】
体積単位の少なくとも1つの蛍光決定ステップでは、第1の発光および/または吸収マッピングとは異なる第2の発光および/または吸収マッピングに従って、第4のモダリティおよび/またはソースのデータに基づいて、蛍光特性を決定できる。
【0047】
言い換えれば、異なるモダリティおよび/またはソースのデータは、例えば、それぞれのデータモダリティおよび/またはソースについて異なる蛍光特性および/または光学特性を得るために、異なるマッピングに従って処理されてもよい。このようにして、異なるデータモダリティおよび/またはソースは、例えば人間または任意のAIツールによって、より容易かつより直感的に区別できる。例えば、
【0048】
発光および吸収マッピングはまた、処理システムのユーザによってカスタマイズ可能であってもよい。このようにして、ユーザは、特定のデータモダリティおよび/またはソースを、例えば強調表示および/または色付けすることができる方法を、自由に選択できる。
【0049】
蛍光決定ステップおよび光学決定ステップは、同時におよび/または並行して行ってもよい。このようにして、本方法の総処理時間を短縮できる。
【0050】
レンダリングステップは、蛍光決定ステップおよび/または光学決定ステップの後に行われてもよい。このように、レンダリングステップは、蛍光決定ステップおよび光学決定ステップで決定されたデータを用いてもよい。例えば、決定された光学特性および蛍光特性は、データストレージ(例えば、ローカルまたはクラウドデータストレージ)に記憶できる。レンダリングステップは、記憶された特性に基づいて1回または数回実行できる。このようにして、例えば3Dモデルのようなデータモデルに決定を記憶し、その後の任意の時点で異なるデータ表現のレンダリングを実行することが可能である。
【0051】
あるいは、レンダリングステップは、蛍光決定ステップおよび/または光学決定ステップと同時に(すなわち、同時期に)行われてもよい。例えば、レンダリングステップにおいて光線追跡を使用する場合、各光線(例えば、データ表現を形成する2D画像内の画素に対応する)について、上記光線によって満たされる体積単位の光学特性および蛍光特性を決定することが可能である。
【0052】
第1のモダリティおよび/またはソースのデータは、Bモード(輝度モード)超音波撮像からの反射率情報を含むことができる。例えば、上記第1のモダリティを使用して、レンダリングされたデータ表現の背景、例えばレンダリングされたデータ表現を形成する2Dまたは3D画像の背景を取得することができる。
【0053】
第2のモダリティおよび/またはソースのデータは、せん断波エラストグラフィ超音波撮像からの組織弾性情報を含むことができる。例えば、第2のモダリティを使用して、レンダリングされたデータ表現内の特定の区域を強調表示することができる。
【0054】
さらなる実施形態によれば、第2または第4のモダリティおよび/またはソースのデータは、介入デバイス情報を含むことができる。上記介入デバイス情報は、媒体に配置された介入デバイス(例えば、その6D姿勢)の位置および/もしくは向き、ならびに/または介入デバイスの形状の事前設定された情報を示すことができる。
【0055】
一例では、介入デバイス情報は、媒体の座標系におけるデバイスの位置および/または向きを示す座標を含むことができる。これらのデータは、同じディスプレイに表示されてもよいし、別のディスプレイに表示されてもよいし、統合シーン内に表示されなくてもよい。したがって、介入デバイス情報は、必ずしもデバイスの画像データを含む必要はないが、そうすることができる。
【0056】
したがって、画像の形態ではないモダリティおよび/またはソースのデータを含むことが可能である。上述したように、介入デバイスの画像情報は決定されず、その座標のみが決定されてもよい。この場合、蛍光特性は、これらの座標に関連付けられてもよく、それらは、媒体および介入デバイスを含む統合シーン内にレンダリングされてもよい。
【0057】
例えば、例示的な介入デバイスとして針を考えると、それをレンダリングするいくつかの可能性がある。
【0058】
針上の点は、例えば、Bモード画像の同じ座標系の下の体積内のボクセルによってボクセルで示され得る。次に、針ソースは、Bモード画像と位置合わせされた別の画像源であってもよく、例えば異なる蛍光色を有する以外は任意の他のモダリティおよび/またはソース(例えばSWE)のデータと同じ方法で処理されてもよい。
【0059】
あるいは、針上の点は、所与の画像、例えばBモード画像の座標系における座標のセットとして表されてもよい。次いで、蛍光色をこれらの座標に関連付けることができる。このように、針先はBモード画像の同じ座標系で表される。この代替案によれば、追加の針画像は使用されない。
【0060】
さらに別の代替形態では、針上の点は、グローバル基準座標系における座標のセットとして表されてもよい。例えば、介入デバイスのトランスデューサまたは送信機(針の先端など)は、その位置を知らせるためにシステムの信号捕捉器に信号を送信することができる。走査プローブはまた、その位置を知らせるために同じ捕捉器に信号を送信することができる。次に、この基準座標系は、一部またはすべてのソースによって共有されてもよい。このようにして、レンダリングシーンは、このグローバル座標系から構築することができる。
【0061】
さらに別の代替形態では、介入デバイスは、本開示の方法のステップにおいて媒体のように、すなわち媒体の一部であるかのように取り扱われる。言い換えれば、介入デバイスおよび媒体に適用される開示された方法にいかなる違いも存在し得ない。特に、同じデータモダリティおよび/またはソースを使用することができ、介入デバイスと媒体とを区別することなく、1つの共通の方法でデータを処理することができる。
【0062】
少なくとも1つの体積単位の蛍光決定ステップでは、介入デバイス情報に基づいて蛍光特性を決定することができる。この決定は、特にその位置および/または向きならびにその形状に従って、介入デバイス(すなわち、介入デバイスが配置されている場所)によって占有される媒体の体積単位を識別することによって行うことができる。
【0063】
したがって、蛍光決定ステップでは、介入デバイスによって占有される媒体の体積単位は、検出された介入デバイスについて決定された蛍光特性に従って着色され得る。このようにして、以下でより詳細に説明するように、例えば視覚的表現において、媒体内の介入デバイスを視覚化することができる。
【0064】
レンダリングステップにおいて、データ表現の少なくとも一部は、介入デバイス情報に基づいてレンダリングされてもよい。上記介入デバイス情報は、媒体内に配置された介入デバイスの位置および/または向きと、介入デバイスの形状の事前設定された情報とを示してもよい。
【0065】
介入デバイスの位置および/または向きは、介入デバイスの事前設定された情報に基づいて構成された走査波を使用して介入デバイスを走査することによって決定することができる。介入デバイスの形状の事前設定された情報は、介入デバイスの3D表面情報を含むことができる。それはまた、介入デバイスの異なる視点からの画像のセット(例えば、2Dまたは3D)を含むことができる。上記情報は、事前走査ステップ(例えば、光学決定ステップおよび/または蛍光決定ステップの前に実行される)で決定されてもよく、またはデータストレージ、例えば処理システムのローカルデータストレージまたはリモートデータストレージ、例えばクラウドシステムから取得されてもよい。データストレージはまた、異なる介入デバイスの介入デバイス情報を記憶することができる。さらに、それは、介入デバイスのさらなる情報、例えば技術的特性(その材料、その空間寸法、その重量など)および/または製造業者情報(製品ID、メーカID等)を記憶することができる。
【0066】
例えば、走査ステップで取得されたBモードデータまたは任意の他の走査データを使用して、介入デバイスの製品IDを決定することができる。この製品IDに基づいて、介入デバイスの形状(潜在的に3D形状)に関する正しい情報をデータストレージから取得することができる。
【0067】
蛍光決定ステップおよび/またはレンダリングステップにおいて介入デバイス情報を使用する結果として、介入デバイス(例えば、生検針)と生検によってアドレス指定された媒体(例えば、解剖学的構造)の両方を視覚化することが有利に可能になる。
【0068】
さらに、介入デバイスが媒体内のどこに移動すると予測されるか(例えば、その移動方向および/または移動経路)を決定することが可能である。上記予測は、それぞれ訓練された(AIベースの)アルゴリズムによって行われ得る。
【0069】
レンダリングステップは、データ表現が光線追跡体積レンダリング方法に従ってレンダリングされる光線追跡レンダリングステップを含むことができる。例えば、蛍光決定ステップおよび光学決定ステップは、それぞれ体積単位ごとに行われてもよい。次に、レンダリングされたデータ表現(例えば、2Dまたは3D画像)の各ピクセルのレンダリングステップにおいて、光線が決定され、ピクセル/ボクセルの色が、光線が満たす体積単位(例えば、ボクセル)の関数として演算される。
【0070】
フローバー、光学特性および蛍光特性は、レンダリングステップの前に決定される必要はない。蛍光決定ステップおよび光学決定ステップは、レンダリングステップと同時に行われてもよい。例えば、決定は、体積単位の値(例えば、ボクセル値)と光学特性および/または蛍光特性との間の事前定義されたマッピング規則であってもよい。そのような場合、それらの決定は、レンダリング自体と同時に行われてもよい。
【0071】
光線追跡レンダリングステップは、以下のサブステップの少なくとも1つを含んでもよい。
・媒体に対して事前定義された幾何学的関係(例えば、媒体の内側、媒体の近くもしくは媒体に隣接している、または媒体から離れている)で仮想光源を配置するステップ。
・被写界深度および/または開口形状を定義するステップ。
・媒体に対して事前定義された幾何学的関係で仮想ビュー平面を配置するステップであって、データ表現が仮想ビュー平面の関数としてレンダリングされる、ステップ。仮想ビュー平面の位置は、後者の実距離に対応するように、媒体からのユーザの実距離の関数として定義できる。これは、拡張現実レンズを使用して、視覚的表現を実際の媒体が位置する領域に一致させる場合に有用であり得る。
【0072】
マルチモダリティおよび/またはマルチソースデータは、媒体の画像データおよび/または3D情報データを含んでもよい。
体積単位は、少なくとも1つのボクセルを含むことができる。ボクセルは、デカルト座標系を形成することができる3D空間内の規則的なグリッド上の値を表すことができる(例えば、モダリティおよび/またはソースのデータがトランスデューサのマトリックス配列を有するプローブによって取得される場合)。したがって、上記3D空間は、媒体および任意選択で1つまたは複数の介入デバイスを含んでもよい。しかしながら、ボクセルは、極座標系に従って配置されてもよい(特に、モダリティおよび/またはソースのデータが、曲線に沿って整列された複数のトランスデューサを含む凸状アレイを有するプローブによって取得される場合)。
【0073】
一般に、本開示の方法で使用されるモダリティおよび/またはソースが異なるボクセル解像度および/または異なる座標系を有する場合、それらは基準システム(体積単位の基準解像度および基準座標系を含む)に標準化され得る。これにより、光学および蛍光決定ステップおよびレンダリングステップにおいて、それらのデータを融合できる。例えば、モダリティおよび/または最高解像度のソースのデータは、基準として使用されてもよく、他のモダリティおよび/またはソースのデータは、それぞれ補間されてもよい。あるいは、モダリティおよび/または最低解像度のソースのデータは、基準として使用されてもよく、他のモダリティおよび/またはソースのデータは、それぞれダウンサンプリングされてもよい。しかしながら、異なる座標系および/またはモダリティおよび/またはソースの異なる解像度を(事前に)調和させることは必須ではない。本開示によれば、それらはまた、レンダリングステップにおいてマージ(したがって、調和)されてもよい。
【0074】
例えば、レンダリングステップにおける各光線が光線追跡法に従って計算されるとき、補間によって異なる解像度のデータを処理できる。実際には、レンダリング中に、各光線を一連の点に離散化できる。したがって、これらの「光線点」の位置は、各データ体積内で既知である。これらの「光線点」でデータ値をサンプリングするために、その元の解像度からモダリティおよび/またはソースのデータを補間できる。値がサンプリングされると、光学/蛍光特性は、例えば、値ー特性マッピングから決定される。異なるモダリティおよび/またはソースのデータ解像度を事前に体系的に統一する必要はない。そのような技術の利点は、例えば、別の低解像度ソースが存在するために、高解像度データを劣化させる必要がないことである。同様に、別の高解像度ソースが存在するために、低解像度データをアップサンプリングする必要はない。したがって、本方法の計算コストを削減することができる。
【0075】
異なる座標系の場合、処理は同様であってもよい。基準座標系を定義または選択することができ、これは望ましくはレンダリングの座標系であり、光線を表すために使用される。そして、データ座標を基準座標に変換してもよい。そのため、各データボクセルは、基準座標系に配置され得る。変換は、事前に行われる必要はない。これは、レンダリング中にオンザフライで行うことができる。要約すると、ボクセル表現データを処理するために、データ解像度を事前に調和させる必要はない。代わりに、各光線がサンプリングされてもよく、データ値は、レンダリング中にオンザフライで補間によってサンプリングされた点で演算されてもよい。
【0076】
光学決定ステップ、蛍光決定ステップ、およびレンダリングステップのうちの少なくとも1つは、実装されてもよく、または少なくとも1つの第1の人工知能ベース(AIベース)アルゴリズムおよび/または(事前訓練された)機械学習アルゴリズムを使用してもよい。したがって、本開示の方法のステップのいずれか1つは、機械学習ベースまたは任意の他のAIアルゴリズムによって実行されてもよい。
【0077】
少なくとも1つの第2の人工知能ベースのアルゴリズムおよび/または予め訓練された機械学習アルゴリズムは、レンダリングされたデータ表現の関数として事前定義されたタスクを実行できる。例えば、上記タスクは、以下の非限定的なリスト、すなわち、回帰タスク、分類タスク、およびセグメント化タスク、または任意の他の所定のタスクのうちの少なくとも1つを含み得る。
【0078】
したがって、本開示のデータ表現は、必ずしも視覚的データ表現である必要はなく、(AIベースの)アルゴリズム、例えば上述のタスクの1つによって処理することができる任意の種類のデータセットであってもよい。
【0079】
これに加えて、データ表現はまた、レンダリングステップにおいて決定された任意のパラメータを含むか、または、レンダリングステップにおいて決定された任意のパラメータからなってもよい。例えば、そのようなパラメータは、介入デバイスの事前設定された点(例えば、生検針の先端)と媒体の事前設定された関心領域(これは、例えば、ユーザによって選択され、および/またはモダリティおよび/またはソースのデータによって、例えば、SWEデータによって、および/または人工知能専用モジュールを使用して決定される)との間の距離を含んでもよい。
【0080】
さらに、レンダリングされたデータ表現は、AIベースのアルゴリズム(すなわち、第2の人工知能ベースのアルゴリズムおよび/または予め訓練された機械学習アルゴリズム)のための入力として使用してもよい。次いで、上記AIベースのアルゴリズムは、元のデータ(すなわち、マルチモダリティおよび/またはマルチソースデータ)を理解または解釈する際にユーザを支援することができる任意の情報を決定してもよい。
【0081】
本開示はさらに、データ処理システムによって実行されると、データ処理システムに本開示による方法を実行させるコンピュータ可読命令を含むコンピュータプログラムに関する。
【0082】
本開示はさらに、体積測定媒体のマルチモダリティおよび/またはマルチソースデータを処理するためのシステム(すなわち、処理システム)に関する。上記システムは、
・媒体の複数の体積単位のうちの少なくとも1つについて、少なくとも第1のモダリティおよび/またはソースのデータに基づいて、光学特性を決定し、
・体積単位の少なくとも1つについて、少なくとも第2のモダリティおよび/またはソースのデータに基づいて、蛍光特性を決定し、
・体積単位の決定された光学特性および蛍光特性に基づいて、体積測定媒体のデータ表現をレンダリングする、
ように構成されている、処理ユニットを含んでもよい。
【0083】
システムは、医療システム、例えば超音波撮像システムであってもよく、または超音波撮像システムを含んでもよい。
【0084】
本開示はまた、コンピュータによって実行されるときに、本開示による方法のステップを実行するための命令を含むコンピュータプログラムを記録した、コンピュータによって読み取り可能な記録媒体に関してもよい。
【0085】
本開示およびその実施形態は、人間、動物、ならびに金属片、砂利、小石などの考慮されるべき任意の材料に、専用の医療デバイスの文脈で使用することができる。
【0086】
他の点で矛盾しない限り、上記の要素と本明細書内の要素との組み合わせを行うことができることを意図している。
【0087】
前述の一般的な説明および以下の詳細な説明の両方は例示的かつ説明的なものにすぎず、例示目的のために提供されており、特許請求の範囲に記載の本開示を限定するものではないことを理解されたい。
【0088】
本明細書に組み込まれ、本明細書の一部を構成する添付の図面は、例示目的のために提供され、本開示の実施形態を説明と共に示し、その原理を支持および例示するのに役立つ。
【図面の簡単な説明】
【0089】
【
図1】本開示の実施形態による、体積測定媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法のフローチャートを示す。
【
図2】
図1の方法の第1の例示的な実施形態を示す。
【
図3】
図1の方法の第2の例示的な実施形態を示す。
【
図4b】
図4aの光線追跡法によって得られた概略画像を示す。
【発明を実施するための形態】
【0090】
ここで、本開示の例示的な実施形態を詳細に参照し、その例を添付の図面に示す。可能な限り、同じまたは同様の部分を指すために図面全体を通して同じ参照番号が使用される。
【0091】
図1は、本開示の実施形態による、体積測定媒体のマルチモダリティおよび/またはマルチソースデータを処理するための方法のフローチャートを示す。
【0092】
任意選択の取得ステップでは、体積測定媒体の異なるモダリティおよび/または異なるソースのSIデータが取得される。一例では、体積測定(すなわち3D)媒体は、生体組織、特に患者または動物または植物のヒト組織であってもよい。
【0093】
媒体の複数の体積単位のうちの少なくとも1つの光学決定ステップS2aにおいて、少なくとも第1のモダリティおよび/またはソースのデータに基づいて光学特性が決定される。言い換えれば、上記ステップは、ステップS1の少なくとも第1のモダリティおよび/またはソース(または例えば、第1および第2のモダリティおよび/またはソース)の入力データを受信し、それらを処理して体積単位の光学特性(例えば、光吸収率)を決定する。
【0094】
体積単位の少なくとも1つの蛍光決定ステップS2bでは、少なくとも第2のモダリティおよび/またはソースのデータに基づいて、蛍光特性が決定される。言い換えれば、上記ステップは、ステップS1の少なくとも第2のモダリティおよび/またはソースの入力データを受信し、それらを処理して体積単位の蛍光特性(例えば、光度および/または色)を決定する。
【0095】
レンダリングステップS3において、体積単位の決定された光学特性および蛍光特性に基づいて、体積測定媒体のデータ表現がレンダリングされる。言い換えれば、上記ステップは、ステップS2aおよびS2bで決定された各体積単位の決定された蛍光および/または光学特性を入力として受信し、それらを処理してデータ表現をレンダリングする。
【0096】
例えば、データ表現は、視覚的表現を含んでもよく、視覚的表現であってもよく、または1つもしくは複数の画像、例えば2次元(2D)もしくは3次元(3D)画像の形態であってもよい。データ表現はまた、複数の深度レベルを有する視覚的またはグラフィック表現を含んでもよい。例えば、それは、媒体の異なる深度レベルの2D画像のセットを含んでもよい。データ表現は、ホログラムもしくは拡張現実シーン、またはユーザがシーン内をナビゲートできる任意の他のシーンをさらに含んでもよい。しかしながら、本開示のデータ表現は、必ずしも視覚的データ表現ではなく、例えば(AIベースの)アルゴリズム、例えば上述のタスクのうちの1つによって処理することができる任意の種類のデータセットであってもよい。
【0097】
レンダリングステップは、従来のレンダリング技術、例えば最大強度投影または光線追跡を使用してもよい。しかしながら、本開示の方法では、蛍光モデリングの追加のステップが、例えば従来の光線追跡フレームワークに追加される。これにより、異なるデータが個別化された蛍光色を数値的に放射し、次いで光学モデルを通って光線追跡され、最終的に統一されたシーンで提示されることが可能になる。
【0098】
本開示において適用可能なレンダリング技術のさらなる例は、多断面再構成(MPR)である。MPRは、任意に配置された直交平面または傾斜平面を使用して、3D体積から2次元(2D)スライスを抽出する画像処理技術である。別の可能な例は、オブジェクト表面が明示的または暗黙的にモデル化されることを必要とする表面レンダリング(SR)である。SRは、媒体内に介入デバイスをレンダリングするために有利に使用されてもよい。次いで、表面を陰影付けし、例えば光線追跡法によってレンダリングする。
【0099】
さらに別の例は、中間幾何学的表現を計算することなく3Dデータ全体を2D画像として表示する、本開示のレンダリングステップに特に有用な体積レンダリング(VR)である。
【0100】
本開示の方法は、いくつかの利点を示唆する。
【0101】
例えば、複数の3D撮像モード(例えば、3D Bモードおよび3D SWE)を視覚化する統一された方法が提供される。さらに、マルチモード3Dデータを視覚化するユーザ体験を向上させることができる。特に、SWEなどの機能情報を用いた超音波体積の直感的な解釈である。また、カスタマイズ可能な視点での体積内ナビゲーションも可能である。それは、例えば仮想眼鏡などを使用して、視覚的ガイダンスのために拡張現実デバイスにさらに統合できる。また、専門家および非専門家の教育または訓練、ならびに同僚間の彼らのコミュニケーションも容易にする。また、自由視点による情報探索が可能となり、撮像方法で強いられた制約が軽減される。したがって、例えば切除などの外科的または治療的使用のために、定性的な視覚的モニタリング能力を提供することができる。
【0102】
本方法は、GPU(すなわち、追加のハードウェアデバイスは不要)などの従来のハードウェアでさらに実施可能であってもよい。さらに、この技術は、オブジェクトまたは表面のセグメント化なしでマルチモード3D画像の融合を可能にする。計算コストはオブジェクトの数に依存せず、したがって、複数のオブジェクトに対して複数のレンダリングを行う必要はない。計算コストは、体積単位の数にのみ依存し、現在の計算能力に適合する。着色は、蛍光の物理的原理に依存し、ユーザの習慣および好みに応じて、もちろんカスタマイズできるとしても、任意の陰影付けには依存しない。これは、より物理的な直感を提供することによって、実用的なパラメータの最適化が容易になる。
【0103】
本方法は、処理システムによって実行されてもよい。そのようなシステムは、処理ユニットと、任意選択的に取得システム(例えば、トランスデューサ)および/または視覚化システム(例えば、ディスプレイ、眼鏡)とを含んでもよい。
【0104】
ただし、取得システムおよび/または視覚化システムが処理システムの外部にあることも可能である。例えば、トランスデューサ(またはトランスデューサのグループ)は、処理システムに遠隔接続可能であってもよい。例示的な一実施形態では、トランスデューサはIOTデバイスであり、かつ/あるいは、IOTデバイスに、ならびに/または処理ユニットおよび/もしくは視覚化システムを形成するスマートフォンに、接続可能である。トランスデューサは、インターネット、「クラウド」、4Gもしくは5Gプロトコル、WIFI、任意のローカルネットワーク、または任意の他のデータ接点もしくはリモート接続を介して、処理システムに接続可能であってもよい。
【0105】
処理システムおよび視覚化システムは、例えば、インターネット、「クラウド」、4Gもしくは5Gプロトコル、WIFI、任意のローカルネットワーク、または任意の他のデータ接点もしくはリモート接続を介して、リモート接続可能であることがさらに可能である。
【0106】
処理ユニットは、例えば、任意選択のバッファメモリと通信する中央処理装置(CPU)および/またはグラフィック処理装置(GPU)、任意選択的に中央処理システムにリンクされたメモリ(MEM)、ならびに、任意選択的に中央処理システムにリンクされたデジタル信号プロセッサ(DSP)、を含んでもよい。
【0107】
任意選択の取得システムは、少なくとも1つのトランスデューサ、例えば、パルスを送信し媒体(すなわち、組織)応答を受信するように構成された単一のトランスデューサ、を含んでもよい。また、複数のトランスデューサおよび/またはトランスデューサアレイを使用することも可能である。アレイは、媒体の二次元(2D)撮像を実行するように適合されてもよいが、アレイはまた、媒体の3D撮像を実行するように適合された二次元アレイであってもよい。トランスデューサアレイはまた、曲線に沿って整列した複数のトランスデューサを含む凸状アレイであってもよい。パルスを送信して応答を受信するために同じトランスデューサが使用されてもよく、または異なるトランスデューサが送信および受信に使用されてもよい。
【0108】
任意選択の視覚化システムは、1つまたは複数のディスプレイ、仮想現実レンズ、または2Dもしくは3D画像、ホログラムもしくは拡張現実シーンを示すための任意の他の電子デバイスを含むことができる。
【0109】
例えば、拡張現実レンズを使用することにより、ユーザは、媒体の視覚的表現を媒体の実際の位置に配置することによって、媒体の内側を仮想的に見ることができる。これは、介入デバイスが媒体に挿入され、関心点(例えば、病変)に向かうその内部の動きが制御される場合に特に有用であり得る。
【0110】
さらに、拡張現実レンズを使用することにより、外科医または任意の他の医師またはユーザは、フリーハンドで媒体の視覚的表現をナビゲートする(例えば、視点および/または視野角を変更する)ことができ、これは、彼/彼女が同時に、(例えば、本開示による介入デバイスを使用して)任意の介入を実行することができることを意味する。この目的のために、拡張現実レンズの動きを感知することができ、対応するセンサ信号は、視覚提示の視点および/または視野角を制御できる。
【0111】
ディスプレイおよび/または拡張現実レンズなどの異なるデバイス上にデータ表現(特に視覚的表現)を提示することも可能である。そのような視覚化デバイスはまた、必要性および/または特定の場合におけるそれらの有用性に応じて(例えば、検査または介入中のユーザの視線方向、視覚的表現の解像度などに応じて)、組み合わせてもよい。
【0112】
処理システムは、超音波撮像のためのデバイスを含んでもよく、または超音波撮像のためのデバイスに接続可能であってもよい。したがって、トランスデューサは超音波トランスデューサであってもよい。しかしながら、処理システムは、超音波以外の波(超音波波長とは異なる波長を有する波)を使用して、任意の撮像デバイスまたは医療システムを含んでもよく、またはそれらに接続可能であってもよい。
【0113】
したがって、マルチモダリティおよび/またはマルチソースデータは、異なるソース、例えば超音波システムおよび別のシステム、例えばコンピュータ断層撮影、磁気共鳴もしくは陽電子放射断層撮影のような医療撮像システム、または電磁(EM)追跡装置のような追跡システムから生じてもよい。
【0114】
さらなる実施形態によれば、モダリティおよび/またはソースデータは、介入デバイス情報を含んでもよい。上記介入デバイス情報は、例えば、以下に説明するように、介入デバイスの撮像および/または追跡のためのシステムから生じてもよい。上記介入デバイス情報は、媒体に配置された介入デバイス(例えば、その6D姿勢)の位置および/または向きを示してもよい。上記介入デバイス情報は、任意選択的に、介入デバイスの形状の事前設定された情報を含んでもよい。位置および/または向きは、上述したように、走査された媒体の座標系において介入デバイスにそれぞれの座標を提供することによって定義されてもよい。座標系はまた、体積単位の構造を事前定義してもよい。
【0115】
異なるデータモダリティの例は、異なる取得モード、例えば、輝度モード(Bモード)、せん断波エラストグラフィ(SWE)、および超音波撮像の場合の歪みエラストグラフィでシステム(例えば、1つの医療システム)によって取得されたデータであり得る。
【0116】
図2は、
図1の方法の第1の例示的な実施形態を示す。言い換えれば、
図2の方法は、特に記載のない限り、
図1の方法に対応する。特に、任意選択のステップ1で取得された(すなわち、走査された)マルチモダリティおよび/またはマルチソースデータは、2つの異なるモダリティ、すなわち、取得された第1のモダリティとしての(3D)Bモード(S1b参照)、および取得された第2のモダリティとしての(3D)SWE(S1a参照)、であってもよい。
【0117】
媒体(すなわち、組織)の複数の体積単位(例えば、ボクセル)の光学決定ステップS2aでは、ステップS1bの第1のモダリティ、およびステップS1aの第2のモダリティのデータに基づいて、光学特性が決定される。
【0118】
特に、光学特性は、体積単位の第1および第2のモダリティおよび/またはソースのデータに適用される組織-光の相互作用マッピングに従って決定される。組織-光の相互作用マッピングは、任意選択的に、反射率、指向性および吸収率のうちの少なくとも1つを決定する。体積単位の決定された光学特性は、光線追跡技術(これに関しては、
図4aおよび
図4bを参照)に従ってステップS3で画像をレンダリングするのに適し得る。
【0119】
媒体の複数の体積単位についての蛍光決定ステップS2bにおいて、ステップS1aの第2のモダリティのデータに基づいて蛍光特性が決定される。言い換えれば、上記ステップは、SWEデータを入力として受信し、それらを処理して体積単位の蛍光特性(例えば、光度および/または色)を決定する。結果として、SWE走査によって検出された媒体内の区域(例えば、病変)は、特定の蛍光特性によってマークされてもよい。蛍光特性は、第2のモダリティのデータに適用される発光および/または吸収マッピングに従って決定されてもよい。
【0120】
さらに、ステップS2Cにおいて、仮想光源は、媒体に対して事前定義された(幾何学的な)位置および/または姿勢で定義および配置されてもよい。上記光源は、ステップS3における光線追跡技術に従って媒体を示すレンダリングされたシーンを照らしてもよい。この目的のために、2Dビュー平面は、媒体に関して事前定義された(幾何学的)位置および/または姿勢で定義および配置されてもよい。
【0121】
レンダリングステップS3では、ステップS2aおよびS2bの体積単位の決定された光学特性および蛍光特性、ならびにステップS2Cで定義された照明特性およびビュー平面に基づいて、体積測定媒体のデータ表現がレンダリングされる。レンダリングステップは、
図4aおよび
図4bの文脈でさらに説明するように、光線追跡技術を使用してもよい。
【0122】
図3は、
図1の方法の第2の例示的な実施形態を示す図である。
図3の方法は、特に記載のない限り、
図2の方法に対応する。
【0123】
特に、
図2の方法に加えて、ステップS1c1およびS1c2が追加されており、例えば生検針からの介入デバイス情報が取得される。上記介入デバイス情報は、媒体に配置された介入デバイス(例えば、その6D姿勢)の位置および/または向きを示してもよい。上記介入デバイス情報は、介入デバイスの形状の事前設定された情報をさらに含んでもよい。
【0124】
例示的な介入デバイスは、針、特に生検針、マーカー、カテーテル、プロテーゼ、マイクロメカニカルシステム(MEMS)MEMS、ステント、弁、バイパス、ペースメーカー、または媒体に挿入され得るもしくは挿入された任意の他のデバイス、を含む。本開示による介入デバイスは、外科的介入に限定されず、外部デバイスが媒体(例えばマーカー)に挿入および/または埋め込まれる任意の介入に限定されることに留意されたい。したがって、介入デバイスは、任意の外部デバイスであってもよい。
【0125】
例示的な一実施形態では、上記介入デバイス情報は、介入デバイス(例えば、本開示による第4のモダリティおよび/またはソースを表す)を追跡および/または撮像するためのシステムによって取得されてもよい。上記システムは、超音波撮像システム、例えば、SWEデータおよび/またはBモードデータを取得するために使用されるシステムと同じであってもよい。あるいは、外部システムであってもよいし、後述するように、さらなる手段(例えば、センサ、トランスデューサまたはマーカー)を含んでもよい。
【0126】
例えば、その先端から信号を発するように構成された針の形態の介入デバイスを使用してもよい。針先の信号を捕捉することにより、先端の位置を特定することができる。したがって、この先端点は、媒体を示すデータ表現のシーン内にレンダリングすることができる。このタイプの視覚化は、生検位置の貴重な情報を、リアルタイムで有利に医師へ提供できる。
【0127】
ステップS1C1(
図3ではモード「ニードルプラス」と呼ばれる)において、媒体内(または媒体に関して、媒体の外側にある場合)の介入デバイスの位置および/または向きは、介入デバイスの事前設定された情報の関数として適合された第1の走査波を使用して介入デバイスを走査することによって、検出または追跡することができる。介入デバイスの事前設定された情報は、例えば、事前設定された波反射および吸収特性を含んでもよい。したがって、走査波は、特に周囲媒体内のその位置および向きを含む介入デバイスをより確実に検出するように適合され得る。この場合、介入デバイスは受動デバイスと見なすことができる。
【0128】
走査は、処理システムに関連付けられたプローブ、および任意選択でBモード走査および/またはSWEデータに使用される超音波システムによって行われてもよい。あるいは、プローブは、超音波以外を用いてもよい。
【0129】
さらなる代替形態によれば、介入デバイスは、その位置および/または向きを検出できるように構成された手段を備えてもよい。この場合、介入デバイスは能動デバイスと見なすことができる。例えば、介入デバイスは、針/マーカー位置を測定および/もしくはシグナリングするためのセンサ、マーカーならびに/またはトランスデューサを備えてもよい。さらなる例では、介入デバイスを追跡および/または撮像するためのシステムは、電磁(EM)追跡装置を含んでもよい。EM追跡装置は、介入デバイスのEMマイクロセンサが追跡されるEM場を生成することができる。
【0130】
ステップS1C2において、介入デバイスによって占有されている媒体内の(3D)領域、すなわち介入デバイスが配置されている場所、が決定される。この目的のために、ステップS1C1で検出された介入デバイスの位置および/または向きが使用される。さらに、介入デバイスの(3D)形状に関する情報が使用される。介入デバイスの形状の情報は、介入デバイスの3D表面情報を含んでもよい。それはまた、介入デバイスの異なる視点からの画像のセット(例えば、2Dまたは3D)を含むことができる。
【0131】
上記情報は、事前の走査ステップ(例えば、光学決定ステップおよび/または蛍光決定ステップの前に実行される)で(事前に)決定されてもよい。
図3の例示的な実施形態に示すように、ステップS1bで取得されたBモードデータを使用して形状を決定してもよい。しかし、他の走査方法で形状を決定してもよい。
【0132】
あるいは、介入デバイスの形状は、データストレージ、例えば処理システムのローカルデータストレージ、または、遠隔データストレージ、例えばクラウドシステム、から取得されてもよい。1つの例示的な利用では、形状は、後でいつでも得ることができるように、一度事前設定され、その後に保存されてもよい。
【0133】
しかしながら、介入デバイスは、場合によっては、介入中に、および/または介入自体の間に、進化し得ることに留意されたい。例えば、機械的な力によって介入中に変形、サイズ縮小、またはねじれる可能性がある。介入はまた、例えばマーカーである場合、放射線または組み込まれた流体(例えば放射性流体)を放射してもよい。したがって、システムは、事前記憶されたデータに基づいて、さらに介入中に得られる「ライブデータ」に基づいて、介入デバイスをレンダリングするように構成されることが望ましい。例えば、そのような「ライブデータ」は、上述したように、ステップS1bにおいて取得されたBモードデータからなることができる。
【0134】
さらに、データストレージは、異なる介入デバイスの介入デバイス情報を記憶してもよい。1つの介入デバイスの情報は、その形状だけでなく、そのさらなるパラメータ、例えば技術的特性(その材料、その空間寸法、その重量など)および/または製造業者情報(製品ID、メーカID等)を含んでもよい。
【0135】
例えば、走査ステップで取得されたBモードデータまたは任意の他の走査データを使用して、介入デバイスの製品IDを決定してもよい。この製品IDに基づいて、介入デバイスの(3D)形状に関する正しい情報をデータストレージから取得することができる。
【0136】
結果として、介入デバイスの3Dモデルが取得され、媒体に対する幾何学的関係に設定される。
【0137】
少なくとも1つの体積単位の蛍光決定ステップS2bでは、介入デバイス情報に基づいて、(介入デバイス固有の)蛍光特性を決定してもよい。この決定は、介入デバイスの位置および/または向きならびにその形状に従って、介入デバイスによって占有されているすべての体積単位を識別することによって行うことができる。特に、これらの体積単位はすべて、検出された介入デバイスについて決定された蛍光特性に従って着色されてもよい。
【0138】
結果として、介入デバイスによって占有される媒体内の区域は、ステップS1aのSWEデータに使用されるものとは異なる特定の蛍光特性によってマークされ得る。さらに、介入デバイスは、媒体および介入デバイスを含むそれぞれの3Dモデルにおいて視覚化されてもよい。
【0139】
さらなる例示的な実施形態によれば、レンダリングステップS2cにおいて、介入デバイス情報に基づいてデータ表現の少なくとも一部がレンダリングされてもよい。上記介入デバイス情報は、媒体内に配置された介入デバイスの位置および/または向きと、介入デバイスの形状の事前設定された情報とを示してもよい。したがって、体積単位の蛍光特性を決定し、その後にレンダリングする代わりに、最終的なデータ表現をレンダリングするために、介入デバイス情報を直接使用することも可能である。結果として、例えば介入デバイスが動くたびにシーン全体(すなわち、完全なデータ表現)をレンダリングする必要はなく、介入デバイスのみをレンダリングすることで、シーンのより効率的なリアルタイム更新が可能になる。
【0140】
例えば、表面レンダリング(SR)方法を使用して介入デバイスがレンダリングされてもよい。この場合、蛍光特性を介入デバイスの表面に適用することができ、SRレンダリングは、例えばステップS1bのBモードデータと共に実行される。したがって、介入デバイスの明示的な表面情報を有するデータと、bモードデータからのセグメント化情報のないデータとを同時に扱うことが可能になる。
【0141】
蛍光決定ステップおよび/またはレンダリングステップにおいて介入デバイス情報を使用する結果として、介入デバイス(例えば、生検針)および媒体(例えば、解剖学的構造)の両方を、視覚化することが可能になる。
【0142】
さらに、介入デバイスが媒体内のどこに移動すると予測されるか(例えば、その移動方向および/または移動経路)を決定することが可能である。上記予測は、それぞれ訓練された(AIベースの)アルゴリズムによって行われ得る。
【0143】
例えば、介入デバイスは、媒体中の関心領域に移動されるバイオプシーガンおよび生検針を含んでもよい。
【0144】
さらに、訓練されたAIベースのアルゴリズムは、ステップS1c1および/またはS1c2で得られたデータに基づいて、介入デバイスのタイプおよび/または形状を予測してもよい。上記予測に基づいて、介入デバイスの画像は、媒体内の介入デバイスの予測位置でレンダリングステップにおいてレンダリングされてもよい。画像は、処理システムのローカルデータストレージに記憶された、または例えばクラウドシステムに遠隔記憶された、介入デバイスの事前設定された画像データに基づいて生成されてもよい。
【0145】
一実施形態では、レンダリングされたデータ表現は、介入中および/または患者の検査(任意選択的に介入を含む)後に、介入の異なる段階で、記憶および/または提示されてもよい。また、検査時に更新されてもよい。例えば、腫瘍または任意の他の病変などの任意の進化の検証を可能にするために、以前の検査および現在の検査のデータ表現を(視覚的におよび/またはパラメータの形態で)提示してもよい。これは、埋め込まれたデバイスの任意の潜在的な進化にも適用されてもよい。例えば、デバイスは、より長い期間、例えば少なくとも2回の検査を伴う期間にわたって、媒体に挿入される1つまたは複数のマーカーまたは任意の他の要素を含んでもよい。そのような場合、上記埋込デバイスの進化、例えば、視覚的提示もしくは数字の形式(例えば、2つ以上の検査間の移動距離)における、マーカーの移動および/または変形をモニタリングすることができる。
【0146】
図4aは、本開示による光線追跡法の概略図を示す。光線追跡レンダリング方法(例えば、
図1~
図3の文脈で説明した方法で使用される)は、以下のサブステップを有してもよい。
1.ステップ1b(
図1~
図3を参照)のBモードデータから、ステップ2aにおいて、Bモード画像のすべての体積単位、例えばすべてのボクセルにおける減衰係数に、Bモードデータをマッピングすることによって、減衰マッピングが演算される。計算上効率的な代替案では、Bモードデータは、光線に沿った体積単位で(のみ)減衰係数にマッピングされ、これは次に後述する光線追跡法に使用される。言い換えれば、マッピングは「光線ごと」であってもよい。「関与する」とは、それぞれのモダリティおよび/またはソース(ここではBモード)におけるデータが提供される体積単位を意味する。
2.ステップS2cの光源データから、光線追跡法をステップ3で適用して、光減衰、および任意選択で局所勾配に関して入射光線の指向性を考慮することによって、関与する各体積単位における入射光強度を導出する。
3.ステップ1aのSWE体積の場合、ユーザは、ステップ2bで3×3の吸収マトリックスMを定義することによって、蛍光色をカスタマイズしてもよい。この行列は、RGB蛍光色成分のそれぞれについて、入射RGB強度の関数として光子吸収率を符号化する。RGBは、可能なカラーコーディングの一例にすぎない。RGBの可能な代替案は、例えばHSV(色相、彩度、明度)である。
4.SWE体積の場合、ユーザは、ステップ2bで3×3の蛍光発光マトリックスEをカスタマイズしてもよい。Eの主対角線上には、RGB成分ごとに発光率が設けられている。
5.次に、関与する各体積単位xにおける入射光強度、すなわち[I
R(x)I
G(x)I
B(x)]
Tが、SWE体積内で蛍光に変換される。
【数1】
ここで、
F
R(x)は、体積単位xが放射する蛍光の赤色成分であり、
F
G(x)は、体積単位xが放射する蛍光の緑色成分であり、
F
B(x)は、体積単位xが放射する蛍光の青色成分であり、
Eは、蛍光発光率を記録する3×3発光マトリックスであり、
Mは、光吸収率を記録する3×3吸収マトリックスであり、
I
R(x)は、体積単位xにおける入射光の赤色成分であり、
I(x)は、体積単位xにおける入射光の緑色成分であり、
I
B(x)は、体積単位xにおける入射光の青色成分である。
6.データ表現でレンダリングするために所与のピクセルを横切る視線Lに沿って、そのRGB強度成分[V
RV
GV
B]
Tは、ステップS3において、光線Lによって満たされるBモード体積単位の内側に到達した入射光源からの光強度と、蛍光発光からの蛍光強度とを組み合わせることによって計算される。
7.これらの2つの強度部分は、ステップS3において、この光線に沿って減衰を加えることによって、視線Lに沿って従来の光線追跡法によって再び変調することができる。強度寄与の2つの部分の減衰係数は異なっていてもよく、カスタマイズ可能であってもよい。
【0147】
図4bは、
図4aの光線追跡法によって得られた概略画像を示す。
図4bは、レンダリングされたデータ表現の例示的なシナリオを示す。この例では、表現は2つの超音波モダリティ、例えば、正常組織を示すための3D Bモードからのモダリティ、および病理組織(すなわち、病変)をユーザ定義または事前定義された蛍光色で強調するための3D SWEモードからのモダリティ、に基づく。
【0148】
合成ボールを含む体積(Bモード体積と見なされる)、および異なる組織弾性の小ボールを含む第2の体積(SWE体積と見なされる)のレンダリングの例が示されている。
【0149】
Bモード体積内のグレー値は、光減衰係数に比例すると考えることができる。したがって、仮想光源はボールの背後であって、むしろ右側に位置するように定義されているため、ボールの左側の領域は減衰の増加により暗く見える。上部小領域は、融合されたSWE情報、例えば、ユーザもしくはAIモジュールによって定義された、またはSWE走査を介して検出された、関心領域を表す。融合されたSWE情報を表す上記領域は、事前定義されたまたはユーザ定義の蛍光色を有してもよい。
【0150】
特許請求の範囲を含む説明全体を通して、「を含む」という用語は、特に明記しない限り、「少なくとも1つを含む」と同義であると理解されるべきである。さらに、特許請求の範囲を含む明細書に記載された任意の範囲は、特に明記しない限り、その最終値を含むと理解されるべきである。記載された要素の特定の値は、当業者に知られている許容された製造公差または工業公差内にあると理解されるべきであり、「実質的に」および/または「およそ」および/または「一般的に」という用語の使用は、そのような許容公差内にあることを意味すると理解されるべきである。
【0151】
本明細書における本開示は、特定の実施形態を参照して説明されているが、これらの実施形態は、本開示の原理および用途の単なる例示であることを理解されたい。
【0152】
本明細書および実施例は、例示としてのみ考慮されることが意図されており、本開示の真の範囲は、以下の特許請求の範囲によって示される。
【国際調査報告】