JP2013503655A5 - - Google Patents

Download PDF

Info

Publication number
JP2013503655A5
JP2013503655A5 JP2012526083A JP2012526083A JP2013503655A5 JP 2013503655 A5 JP2013503655 A5 JP 2013503655A5 JP 2012526083 A JP2012526083 A JP 2012526083A JP 2012526083 A JP2012526083 A JP 2012526083A JP 2013503655 A5 JP2013503655 A5 JP 2013503655A5
Authority
JP
Japan
Prior art keywords
function
selected position
difference
data processing
positive
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2012526083A
Other languages
English (en)
Other versions
JP5628314B2 (ja
JP2013503655A (ja
Filing date
Publication date
Application filed filed Critical
Priority claimed from PCT/EP2010/062816 external-priority patent/WO2011026866A1/en
Publication of JP2013503655A publication Critical patent/JP2013503655A/ja
Publication of JP2013503655A5 publication Critical patent/JP2013503655A5/ja
Application granted granted Critical
Publication of JP5628314B2 publication Critical patent/JP5628314B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Description

特に、本発明の態様は、患者の臓器を分析するデータ処理方法を提示する。データ処理方法は、以下のステップを含む。分析期間における臓器を表す一連のインプット画像が、提供される。個々のインプット画像は、分析期間内に含まれる対応する取得の瞬間における臓器の対応する位置の呼び掛け信号(例えば、超音波)に対する応答を、各々が示す(例えば、ピクセルのレベルでの、若しくはピクセルのグループでの)インプット値のセットを含む。時間に関する分析関数が、(例えば、関心領域の、若しくはピクセルのグループの)選択された位置のセットのうちの各々と、関連付けられる。分析関数は、(例えば、適合処理により決定されるものとしての)一連のインプット画像内の選択された位置のインプット値の傾向をモデル化するものである。時間に関する参照関数が、(例えば、健康な組織を含む臓器の領域のインプット値を適合させることにより)分析関数に対して提供される。本発明の実施形態に係る解決手段では、個々の選択された位置の分析関数が、参照関数と比較される。この動作は、選択された位置の分析関数と参照関数との間の相違の極性の分析期間における傾向を表す極性傾向を判別する(言い換えれば、それは、個々の瞬間で分析関数が参照関数より大きいか若しくは小さいかを示す相違の極性が、分析期間において経時的にどのように変化するかを表す。例えば、相違が、分析関数と参照関数の間の差異として決定されるならば、極性の傾向は、この差異がヌルのとき、常に正のとき、常に負のとき、正から負に変化するとき、若しくはその逆のときを、示す。一方、相違が、分析関数と参照関数の間の比率として決定されるならば、極性の傾向は、この比率が1のとき、常に1より大きいとき、常に1より小さいとき、1より大きいものからいより小さいものに変化するとき、若しくはその逆のときを、示す。)。続いてパラメトリック画像が生成される。パラメトリック画像は、個々の選択された位置に対するパラメトリック値を含む。パラメトリック値は、(例えば、極性傾向のクラスを表す値、及び/又はその計測により)選択された位置の極性傾向を示す。
本発明の一つの実施形態では、選択された位置を分類するステップが、差異関数が正である差異関数に関する、分析期間における積分に従って、正のエネルギを計算するステップと、差異関数が負である差異関数の絶対値に関する、分析期間における積分に従って、負のエネルギを計算するステップとを含む。(例えば、差異サンプル若しくは差異関数が夫々正若しくは負であるとき、分析期間における、差異サンプルの合計に若しくは差異関数の積分に設定される正のエネルギ若しくは負のエネルギによるものである。)続いて選択された位置は、正のエネルギと負のエネルギの合計が閾値より小さい場合ヌルクラスに分類され、正のエネルギ(例えば、相対的正のエネルギ)が更なる閾値を超える場合、正の単極性クラスに分類され、負のエネルギ(例えば、相対的負のエネルギ)が更なる閾値を超える場合、負の単極性クラスに分類され、及び/又は、それら以外の場合、正−負の両極性クラス若しくは負−正の両極性クラスに分類される。
本発明の一つの実施形態では、選択された位置を正−負の両極性クラス若しくは負−正の両極性クラスに分類するステップが、差異関数が分析期間内に最大値に到達する正のピークの瞬間を計算するステップと、差異関数が分析期間内に最小値に到達する負のピークの瞬間を計算するステップとを含む。続いて選択された位置は、最大ピークの瞬間が最小ピークの瞬間を超える場合、正−負の両極性クラスに分類され、最小ピークの瞬間が最大ピークの瞬間を超える場合、負−正の両極性クラスに分類される。
本発明の一つの実施形態では、パラメトリック画像を生成するステップが、個々の選択された位置に対して、強度値を計算するステップを含む。強度値は、分析期間における選択された位置の分析関数と参照関数との間の相違を計測する。個々の選択された位置のパラメトリック値は、対応する強度値に従って設定される。
本発明の別の態様は、患者の臓器を分析する治療方法を提示する。患者の臓器を分析する治療方法は以下のステップを含む。呼び掛け信号が、分析期間において臓器に加えられる。方法は、分析期間において、臓器を表す一連のインプット画像を取得することにより、継続する。個々のインプット画像は、分析期間内に含まれる対応する取得の瞬間における臓器の対応する位置の呼びかけ信号に対する応答を、各々が示すインプット値のセットを含む。(インプット画像は、時間に関する分析関数を、選択された位置のセットの各々と関連付けるように処理され、分析関数は、一連のインプット画像内の選択された位置のインプット値の傾向をモデル化し、個々の選択された位置の分析関数を、選択された位置の分析関数と参照関数との間の相違の極性に関する、分析期間における傾向を表す極性傾向を判別する、時間に関する参照関数と比較し、個々の選択された位置に対するパラメトリック値を含むパラメトリック画像を生成し、パラメトリック値は選択された位置の極性傾向を示す。)臓器の状況は、パラメトリック画像に従って評価される。
本発明の実施形態に係る解決手段では、分析処理における、対応するモデル関数と参照関数との間の比較に従って、個々のピクセルに対するパラメトリック値が計算される。対応するパラメトリック値に基づく値を個々のピクセルに割り当てることにより、パラメトリック画像が生成される。特に、この比較は、極性傾向に基づくものであり、この極性傾向は、モデル関数と参照関数との間の相違の極性の分析処理における、傾向を表す(例えば、対応する差異関数の符号が、時間に対してどのように変化するかを表す)。特に、本発明の実施形態では、この極性傾向の定性的インジケータ及び/又は定量的インジケータが計算される。
別途、分析関数を参照関数と比較(し、それらの相違の極性傾向を判別)することも、可能である。例えば、この動作は、(差異サンプルを計算すること無く)個々の分析サンプルを、対応する参照サンプルと比較するのみでも実行され得る。
同様に、強度値は、モデル関数と参照関数との間の相違(例えば、差異)に関する、他のどの手段に従って、計算されてもよい。

Claims (17)

  1. 患者の臓器を分析するデータ処理方法(A1−A14)において、
    前記データ処理方法は、
    分析期間における臓器を表す一連のインプット画像を提供するステップ(A1)であって、個々のインプット画像は、前記分析期間内に含まれる対応する取得の瞬間における臓器の対応する位置の呼び掛け信号に対する応答を、各々が示すインプット値のセットを含む、ステップと、
    選択された位置のセットのうちの各々と、時間に関する分析関数を関連付けるステップ(A2、A3、A4、A5)であって、前記分析関数は、一連のインプット画像内の選択された位置のインプット値の傾向をモデル化するものである、ステップと、
    前記分析関数に対して時間に関する参照関数を提供するステップ(A2’、A4’、A5、A6’)とを含み、
    個々の選択された位置の分析関数を、選択された位置の分析関数と参照関数との間の相違の極性の分析期間における傾向を表す極性傾向を判別する参照関数と、比較し、前記相違の極性は、選択された位置の分析関数の特徴が参照関数の特徴よりも大きい、等しい、または小さいかを示す、ステップ(A7)と、
    個々の選択された位置に対するパラメトリック値を含むパラメトリック画像を生成するステップ(A8−A13)であって、パラメトリック値は、選択された位置の極性傾向を示す、ステップとにより特徴付けられる
    方法。
  2. 極性傾向が、選択された位置の分析関数と参照関数との間の差異の符号に関する、分析期間における傾向を表す
    請求項1に記載のデータ処理方法(A1−A14)。
  3. 個々の選択された位置の分析関数を参照関数と比較するステップ(A7)が、
    分析期間における、選択された位置の分析関数と参照関数との間の差異に基づいて、差異関数の表現を計算するステップであって、個々の選択された位置に対するパラメトリック値が、選択された位置の差異関数の符号に関する、分析期間における傾向に従って計算される(400)、ステップを含む
    請求項2に記載のデータ処理方法(A1−A14)。
  4. 差異関数の表現を計算する(A7)ステップが、
    差異関数の離散表現を計算する(A’)ステップであって、差異関数の離散表現は一連の差異サンプルを含み、それらの各々は、対応するサンプリングの瞬間における分析関数の値と参照関数の値との間の差異に等しい、ステップを含む
    請求項3に記載のデータ処理方法(A1−A14)。
  5. 差異関数の表現を計算する(A7)ステップが、
    前記一連の差異サンプルを適合させることにより差異関数のパラメトリック表現を計算するステップ(415)を含む
    請求項4に記載のデータ処理方法(A1−A14)。
  6. パラメトリック画像を生成する(A8−A13)ステップが、
    選択された位置の極性傾向に従って、複数の所定のクラスの間の一つに、個々の選択された位置を分類する(A8)ステップと、
    対応するクラスに従って、個々の選択された位置のパラメトリック値を設定する(A10)ステップとを含む
    請求項1乃至5のうちのいずれか一に記載のデータ処理方法(A1−A14)。
  7. 個々の選択された位置を分類する(A8)ステップが、
    差異関数が分析期間内にて主として正である場合、選択された位置を正の単極性クラスに分類し、
    差異関数が分析期間内にて主として負である場合、選択された位置を負の単極性クラスに分類し、
    差異関数が分析期間の第1の部分内にて主として正であり、第1の部分に続く分析期間の残りの第2の部分内にて主として負である場合、選択された位置を正−負の両極性クラスに分類し、
    差異関数が分析期間の更なる第1の部分内にて主として負であり、更なる第1の部分に続く分析期間の残りの更なる第2の部分内にて主として正である場合、選択された位置を負−正の両極性クラスに分類し、及び/又は、
    差異関数が分析期間内にて実質的にゼロである場合、選択された位置をヌルクラスに分類するステップ(406−469)を含む
    請求項6に記載のデータ処理方法(A1−A14)。
  8. 選択された位置を分類するステップ(406−469)が、
    差異関数が正である差異関数に関する、分析期間における積分に従って、正のエネルギを計算するステップ(409、430;415−418、430)と、
    差異関数が負である差異関数の絶対値に関する、分析期間における積分に従って、負のエネルギを計算するステップ(409、430;415−418、430)と、
    正のエネルギと負のエネルギの合計が閾値より小さい場合、選択された位置をヌルクラスに分類し、正のエネルギが更なる閾値を超える場合、選択された位置を正の単極性クラスに分類し、負のエネルギが更なる閾値を超える場合、選択された位置を負の単極性クラスに分類し、及び/又は、それら以外の場合、選択された位置を正−負の両極性クラス若しくは負−正の両極性クラスに分類するステップ(436−469)とを含む
    請求項7に記載のデータ処理方法(A1−A14)。
  9. 選択された位置を正−負の両極性クラス若しくは負−正の両極性クラスに分類するステップ(436−469)が、
    差異関数が分析期間内に最大値に到達する正のピークの瞬間を計算するステップ(451;457)と、
    差異関数が分析期間内に最小値に到達する負のピークの瞬間を計算するステップ(451;457)と、及び、
    最大ピークの瞬間が最小ピークの瞬間を超える場合、選択された位置を正−負の両極性クラスに分類し、最小ピークの瞬間が最大ピークの瞬間を超える場合、選択された位置を負−正の両極性クラスに分類するステップ(463−469)とを含む
    請求項8に記載のデータ処理方法(A1−A14)。
  10. パラメトリック画像を生成する(A8−A13)ステップが、
    個々の選択された位置に対して、分析期間における選択された位置の分析関数と参照関数との間の相違を計測する強度値を計算(A8)し、対応する強度値に従って個々の選択された位置のパラメトリック値を設定するステップを含む
    請求項1乃至9のうちのいずれか一に記載のデータ処理方法(A1−A14)。
  11. 個々の選択された位置に対する強度値を計算(A8)するステップが、
    差異関数の絶対値に関する、分析期間における、積分に従って、強度値を設定するステップを含む
    請求項3を直接若しくは間接に引用する際の請求項10に記載のデータ処理方法(A1−A14)。
  12. パラメトリック画像をディスプレイ(A14)するステップであって、パラメトリック画像の個々のパラメトリック値が、対応するクラスのためのクラス可視化ディメンジョン及び/又は対応する強度値のための強度可視化ディメンジョンを有するグラフ表現により、ディスプレイされる、ステップを更に含む
    請求項6乃至11のうちのいずれか一に記載のデータ処理方法(A1−A14)。
  13. クラス可視化ディメンジョンは、各々が対応するクラスのためのものである複数のカラーを含み、強度可視化ディメンジョンは、強度値に対応する輝度を含む
    請求項12に記載のデータ処理方法(A1−A14)。
  14. 参照関数を提供する(A2’、A4’、A5、A6’)ステップが、
    参照位置のセットを選択するステップと、
    一連のインプット画像内の参照位置のインプット値の傾向をモデル化する参照関数を決定する(A5,A6’)ステップとを含む
    請求項1乃至13のうちのいずれか一に記載のデータ処理方法(A1−A14)。
  15. 参照関数を決定する(A5,A6’)ステップが、
    個々のインプット画像内の参照位置のインプット値の平均値を計算するステップ(A5)と、
    一連のインプット画像内の平均値の傾向をモデル化する参照関数を決定するステップ(A6’)とを含む
    請求項14に記載のデータ処理方法(A1−A14)。
  16. データ処理システム上で稼働する際、請求項1乃至15のうちのいずれか一に記載のデータ処理方法(A1−A14)のステップをデータ処理システムに実行させるコード手段を含むコンピュータプログラム(300)。
  17. 請求項1乃至15のうちのいずれか一に記載のデータ処理方法(A1−A14)のステップを実行するための手段(300)を含む診断システム(100)。
JP2012526083A 2009-09-01 2010-09-01 経時的な動的挙動に基づくパラメトリック画像 Active JP5628314B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
EP09169189.9 2009-09-01
EP09169189 2009-09-01
EP09170731.5 2009-09-18
EP09170731 2009-09-18
PCT/EP2010/062816 WO2011026866A1 (en) 2009-09-01 2010-09-01 Parametric images based on dynamic behavior over time

Publications (3)

Publication Number Publication Date
JP2013503655A JP2013503655A (ja) 2013-02-04
JP2013503655A5 true JP2013503655A5 (ja) 2014-07-24
JP5628314B2 JP5628314B2 (ja) 2014-11-19

Family

ID=42797469

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012526083A Active JP5628314B2 (ja) 2009-09-01 2010-09-01 経時的な動的挙動に基づくパラメトリック画像

Country Status (7)

Country Link
US (1) US8929634B2 (ja)
EP (1) EP2473972B1 (ja)
JP (1) JP5628314B2 (ja)
CN (1) CN102483847B (ja)
AU (1) AU2010291254B2 (ja)
CA (1) CA2769164C (ja)
WO (1) WO2011026866A1 (ja)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10130342B2 (en) 2007-12-28 2018-11-20 Bracco Suisse Sa Initialization of fitting parameters for perfusion assessment based on bolus administration
CN103228214B (zh) * 2010-11-29 2015-04-22 株式会社日立医疗器械 图像处理装置及图像处理方法
WO2014024453A1 (ja) * 2012-08-07 2014-02-13 パナソニック株式会社 医用データ処理装置、医用データ処理方法、及び超音波診断装置
WO2014054810A1 (ja) 2012-10-04 2014-04-10 株式会社 東芝 超音波診断装置
US20140172459A1 (en) * 2012-12-14 2014-06-19 Koninklijke Philips Electronics N.V. Clinical support system and method
EP2936433B1 (en) * 2012-12-21 2018-09-19 Bracco Suisse SA Segmentation in diagnostic imaging applications based on statistical analysis over time
US10368842B2 (en) 2014-04-07 2019-08-06 Bracco Suisse S.A. Estimation of acoustic level in-situ with non-fundamental analysis
CN107106109B (zh) * 2014-11-06 2020-11-06 皇家飞利浦有限公司 计算机断层扫描系统
US10631806B2 (en) * 2015-02-19 2020-04-28 Canon Medical Systems Corporation X-ray diagnostic apparatus, image processing apparatus, and image processing method
CN108289654B (zh) 2015-12-10 2021-03-30 博莱科瑞士股份有限公司 通过动态阈值化来检测固定化造影剂
JP7308196B2 (ja) * 2017-11-22 2023-07-13 コーニンクレッカ フィリップス エヌ ヴェ 超音波肺評価
DE102018215958A1 (de) * 2018-09-19 2020-03-19 Siemens Healthcare Gmbh Verfahren zur Ermittlung eines Befähigungsvorgangs, diesen nutzenden Verfahren zur Einstellung eines dosisbezogenen Aufnahmeparameters, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5445813A (en) 1992-11-02 1995-08-29 Bracco International B.V. Stable microbubble suspensions as enhancement agents for ultrasound echography
IN172208B (ja) 1990-04-02 1993-05-01 Sint Sa
AU636481B2 (en) 1990-05-18 1993-04-29 Bracco International B.V. Polymeric gas or air filled microballoons usable as suspensions in liquid carriers for ultrasonic echography
IL104084A (en) 1992-01-24 1996-09-12 Bracco Int Bv Sustainable aqueous suspensions of pressure-resistant and gas-filled blisters, their preparation, and contrast agents containing them
CN1068229C (zh) 1993-12-15 2001-07-11 勃勒柯研究有限公司 超声对比介质、含该介质的对比剂及方法
JPH08336531A (ja) 1995-06-12 1996-12-24 Hitachi Medical Corp 超音波診断装置
WO2004110279A1 (en) 2003-06-12 2004-12-23 Bracco Research Sa Blood flow estimates through replenishment curve fitting in ultrasound contrast imaging
JP4494744B2 (ja) 2003-09-25 2010-06-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
CN101027693B (zh) 2004-08-09 2010-05-12 伯拉考国际股份公司 基于多个掩码的用于医学成像的像对准方法及设备
EP1855596B1 (en) * 2005-02-23 2015-07-01 Koninklijke Philips N.V. Ultrasonic diagnostic imaging system for detecting lesions of the liver
AU2006234326B2 (en) 2005-04-14 2011-06-16 Bracco Suisse S.A. Perfusion assessment based on animated perfusion imaging
EP2474268B1 (en) * 2005-05-06 2013-07-10 Vasonova, Inc. Apparatus for endovascular device guiding and positioning
US9198639B2 (en) 2005-11-10 2015-12-01 Bracco Suisse S.A. Detection of immobilized contrast agent in medical imaging applications based on flow dynamics analysis
JP2009028194A (ja) * 2007-07-26 2009-02-12 Ge Medical Systems Global Technology Co Llc 超音波撮像装置
JP2009100971A (ja) 2007-10-24 2009-05-14 Ge Medical Systems Global Technology Co Llc 超音波撮像装置
EP2234543B1 (en) 2007-12-28 2016-11-02 Bracco Suisse SA Quantification analisys of immobilized contrast agent in medical imaging applications
JP5455379B2 (ja) * 2009-01-07 2014-03-26 株式会社東芝 医用画像処理装置、超音波診断装置、及び医用画像処理プログラム
JP5593382B2 (ja) 2009-06-08 2014-09-24 ブラッコ・シュイス・ソシエテ・アノニム パラメトリック画像の自動スケーリング

Similar Documents

Publication Publication Date Title
JP2013503655A5 (ja)
JP5672144B2 (ja) 心拍数・呼吸数検出装置,方法およびプログラム
JP5847735B2 (ja) 視覚機能検査方法、視覚機能検査システム及びソフトウェア
CN103260526B (zh) 具有峰值强度检测功能的超声成像系统和方法
JP2004174254A5 (ja)
JP5576775B2 (ja) 画像処理装置、画像処理方法、及び画像処理プログラム
JP5828649B2 (ja) 画像処理装置、画像処理方法、及びコンピュータプログラム
RU2011130813A (ru) Система визуализации акустических изображений и способ визуализации акустических изображений с квантификацией контраста
Haase et al. Automated and objective action coding of facial expressions in patients with acute facial palsy
US20220309668A1 (en) Using a Set of Machine Learning Diagnostic Models to Determine a Diagnosis Based on a Skin Tone of a Patient
Liu et al. Simple model for encoding natural images by retinal ganglion cells with nonlinear spatial integration
Zheng et al. Corneal‐Smart Phone: A novel method to intelligently estimate postmortem interval
US20150310254A1 (en) System and method for automated identification of abnormal ciliary motion
Balkenius et al. Comparing analysis methods in functional calcium imaging of the insect brain
US11645859B2 (en) Analysis device, analysis method, analysis program and display device
Penttinen et al. Deducing self-interaction in eye movement data using sequential spatial point processes
TW202111617A (zh) 區域性腦年齡預測之方法、非暫時性電腦可讀媒體及設備
JP7385046B2 (ja) 色斑予測方法、装置、設備及び記憶媒体
JP5524589B2 (ja) 診断支援システム、方法及びコンピュータプログラム
EP2980757B1 (en) Quantification and imaging methods of the echo-texture feature
JP5924406B2 (ja) 画像計測装置、画像計測方法、および、画像計測用プログラム
WO2015125759A1 (ja) 加齢分析方法及び加齢分析装置
Ahmed et al. Automatic Region of Interest Extraction from Finger Nail Images for Measuring Blood Hemoglobin Level
JP2012115319A (ja) 器官運動解析装置及び器官運動解析プログラム
WO2013041977A1 (en) Method of retinal image enhancement and tool therefor