JP6502367B2 - 異方性柔軟媒体を特性評価するための方法及び超音波デバイス、並びにそのような特性評価デバイス用の超音波プローブのセット - Google Patents

異方性柔軟媒体を特性評価するための方法及び超音波デバイス、並びにそのような特性評価デバイス用の超音波プローブのセット Download PDF

Info

Publication number
JP6502367B2
JP6502367B2 JP2016549242A JP2016549242A JP6502367B2 JP 6502367 B2 JP6502367 B2 JP 6502367B2 JP 2016549242 A JP2016549242 A JP 2016549242A JP 2016549242 A JP2016549242 A JP 2016549242A JP 6502367 B2 JP6502367 B2 JP 6502367B2
Authority
JP
Japan
Prior art keywords
propagation
shear wave
ultrasonic
transducer
flexible medium
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.)
Active
Application number
JP2016549242A
Other languages
English (en)
Other versions
JP2017504435A (ja
Inventor
ミカエル・タンター
マチュー・ペルノー
マティアス・フィンク
ジャン−リュック・ジェニソン
Original Assignee
センター ナショナル ド ラ ルシェルシュ サイエンティフィーク
センター ナショナル ド ラ ルシェルシュ サイエンティフィーク
アンセルム(アンスティチュート・ナシオナル・ドゥ・ラ・サンテ・エ・ドゥ・ラ・ルシェルシュ・メディカル)
エコール・シュペリュール・ドゥ・フィシック・エ・シミー・アンデュストリエル・ドゥ・ラ・ヴィル・ドゥ・パリ
ユニヴェルシテ パリ ディドロ パリ 7
ユニヴェルシテ パリ ディドロ パリ 7
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by センター ナショナル ド ラ ルシェルシュ サイエンティフィーク, センター ナショナル ド ラ ルシェルシュ サイエンティフィーク, アンセルム(アンスティチュート・ナシオナル・ドゥ・ラ・サンテ・エ・ドゥ・ラ・ルシェルシュ・メディカル), エコール・シュペリュール・ドゥ・フィシック・エ・シミー・アンデュストリエル・ドゥ・ラ・ヴィル・ドゥ・パリ, ユニヴェルシテ パリ ディドロ パリ 7, ユニヴェルシテ パリ ディドロ パリ 7 filed Critical センター ナショナル ド ラ ルシェルシュ サイエンティフィーク
Publication of JP2017504435A publication Critical patent/JP2017504435A/ja
Application granted granted Critical
Publication of JP6502367B2 publication Critical patent/JP6502367B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4272Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue
    • A61B8/4281Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue characterised by sound-transmitting media or devices for coupling the transducer to the tissue

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Vascular Medicine (AREA)
  • Gynecology & Obstetrics (AREA)
  • Acoustics & Sound (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Description

本発明は、異方性柔軟媒体を特性評価するための方法及びデバイス、並びにそのような特性評価デバイス用の一組の超音波プローブのセットに関する。
骨格筋や心筋等の繊維を含む柔軟媒体の特性評価は、その媒体内のせん断波の伝播をイメージングすることによるものが既に提案されている。特に、繊維に沿った及び繊維に垂直な媒体の剛性率の値を測定することが、二回のせん断波放出の間に数度回転される超音波トランスデューサアレイを用いて、媒体内にせん断波を連続的に伝播させて、異なる方向でのその伝播を毎回イメージングすることによるものが提案されている(非特許文献1)。このようにして、せん断波の伝播速度が各方向において決定される。非特許文献2に教示されているように、速度の測定結果を特異値分解することによって、繊維の方向と、繊維に沿った及び繊維に垂直な剛性率とをどちらも得ることができる。
しかしながら、これら既知の方法は、比較的遅いという欠点を有するので、特に、患者の心筋や患者の骨格筋等の動く繊維媒体を特性評価することができない。
国際公開第2004/021038号 国際公開第2000/055616号
W.‐N.Lee,M.Pernot,M.Couade,E.Messas,P.Bruneval,A.Bel,A.A.Hagege,M.Fink,M.Tanter,"Mapping Myocardial Fiber Orientation Using Echocardiography‐Based Shear Wave Imaging",IEEE Trans.Med.Imaging,Vol.31,pp. 554−562,2012 W.‐N.Lee,B.Larrat,M.Pernot,M.Tanter,"Ultrasound Elastic Tensor Imaging: Comparison with MR Diffusion Tensor Imaging in the Myocardium",Physics in Medicine and Biology,Vol.57,pp.5075−5095,2012 Royer,Dieulesaint,2000; Elastic Waves in Solids I: Free and Guided Propagation: Springer−Verlag Berlin Heidelberg
本発明は、特に、こうした欠点を解消するという目的を有する。
このため、本発明は、繊維を含む少なくとも一つの部分と外面とを有する異方性柔軟媒体を特性評価するための方法を提案し、本方法は以下のステップを備える:
(a)測定ステップ。この測定ステップにおいては、異方性柔軟媒体内の中心領域から発散する波として伝播する少なくとも一つのせん断波を発生させて、前記異方性柔軟媒体の表面から複数の超音波観測トランスデューサを用いて、それら超音波観測トランスデューサを固定したままで、中心領域からの複数の所定の伝播方向における少なくとも一つのせん断波の伝播を観測する。複数の所定の伝播方向は、0度及び180度とは異なる角度を成す少なくとも二つの方向を有する。超音波観測トランスデューサは、少なくとも所定の伝播方向に沿って位置する。測定ステップは、50ms未満の期間内に行われる。
(b)少なくとも一つの計算ステップ。この計算ステップにおいては、測定ステップ(a)で収集されたデータから、複数の所定の伝播方向の各々についてせん断波の少なくとも一つの伝播パラメータを決定する。
(c)特性評価ステップ。この特性評価ステップにおいては、計算ステップ(b)で伝播方向の各々について決定されたせん断波の少なくとも一つの伝播パラメータから、異方性柔軟媒体の繊維の方向と、繊維に垂直な方向におけるレオロジー弾性パラメータと、繊維の方向におけるレオロジー弾性パラメータとのうちから選択された異方性柔軟媒体の少なくとも一つのレオロジー特性を決定する。
こうした構成を用いることで、測定プローブを回転させずに、繊維の方向及び/又は弾性パラメータを即座に得ることが可能になり、つまりは、in vivoで測定を容易に行うことが可能になる。
本発明に係る方法の多様な実施形態においては、任意で、以下の構成を更に適用することができる:
‐ 特性評価ステップ(c)で決定されるレオロジー弾性パラメータは弾性率である;
‐ 測定ステップ(a)において、異方性柔軟媒体内の特定の深さ範囲にわたってせん断波を発生させ、その深さ範囲内の複数の異なる深さにおいて少なくとも一つのせん断波の伝播を観測する。計算ステップ(b)において、複数の所定の伝播方向(P)の各々についてのせん断波の少なくとも一つの伝播パラメータを複数の異なる深さにおいて決定する。そして、特性評価ステップ(c)において、異なる複数の深さにおける少なくとも一つのレオロジー特性を決定する;
‐ 測定ステップ(a)において、せん断波の伝播が観測される所定の伝播方向の数は、3個から20個の間、有利には5個から10個の間である;
‐ 測定ステップ(a)において、複数の超音波観測トランスデューサを全て同時に用いて、複数の所定の伝播方向の全てについて、少なくとも一つのせん断波の伝播を同時に観測する;
‐ 測定ステップ(a)において、複数のせん断波を連続的に放出して、複数の超音波観測トランスデューサのうち一部を用いて、複数の所定の伝播方向の少なくとも一つにおける各せん断波の伝播を連続的に観測する;
‐ 測定ステップ(a)において、複数の所定の伝播方向に沿ってそれぞれ位置する複数の所定の測定点におけるせん断波の通過を検出し、計算ステップ(b)において決定されるせん断波の伝播パラメータは、せん断波の伝播速度、又は、中心領域からのせん断波の飛行時間のいずれかを表すパラメータである;
‐ 異方性柔軟媒体の表面に垂直な中心軸からせん断波を発生させて、せん断波の通過が検出される複数の測定点はそれぞれ中心軸から2cm未満、有利には0.5cm未満の距離に位置する;
‐ 測定ステップ(a)において、超音波観測トランスデューサによる少なくとも毎秒300ショットの超音波圧縮波の率での音響圧縮波の放出によって、所定の測定点におけるせん断波の通過を検出し、異方性柔軟媒体によって反射された超音波信号を超音波観測トランスデューサによって感知することによって、所定の測定点でのせん断波の通過の際の異方性柔軟媒体の内部変位(絶対変位若しくは変形、又は変位速度若しくは変形速度)を決定する;
‐ 測定ステップ(a)において、測定点毎に単一の超音波観測トランスデューサを使用する;
‐ 測定ステップ(a)において、超音波観測トランスデューサを含む超音波励起トランスデューサが、中心領域に集束させた励起超音波を放出して、中心軸に沿って異方性柔軟媒体を変位させることによって、せん断波を発生させる;
‐ 特性評価ステップ(c)において、特異値(異なる複数の所定の伝播方向における伝播パラメータの値)に対応する対応弾性テンソルを決定して、その弾性テンソルを特異値分解することによって、少なくとも一つのレオロジー特性を決定する。
‐ 測定ステップ(a)において、中心軸からせん断波を発生させて、特性評価ステップ(c)において、計算ステップ(b)で計算された所定の伝播方向における伝播パラメータの値に応じて、計算ステップ(b)で計算された所定の伝播方向における伝播パラメータの値の補間によって、略楕円形の曲線C(V(θ)・cosθ,V(θ)・sinθ)を決定し、ここで、V(θ)は、中心軸に垂直な平面内の伝播パラメータの値であり、θは、その平面に属する基準系に対する相対的な平面内の伝播方向を示す角度であり; V(θ)の最大値に対応する角度θに対応する繊維の方向、V(θ)に従う繊維の方向におけるレオロジー弾性パラメータ、V(θ+π/2)に従う繊維に垂直な方向におけるレオロジー弾性パラメータを決定する;
‐ 異方性柔軟媒体は、動いている人間又は動物の筋肉の少なくとも一部を有し、本方法は、測定ステップ(a)、計算ステップ(b)、及び特性評価ステップ(c)を複数回行って、そこから、筋肉の収縮に関係する生体パラメータを推定する。
本発明の更なる対象は、繊維を含む少なくとも一つの部分と外面とを有する異方性柔軟媒体を特性評価するためのデバイスであり、その特性評価デバイスは、励起プローブ及び複数の超音波観測トランスデューサを制御する電子制御デバイスを備え、励起プローブは、中心領域から異方性柔軟媒体内にせん断波を発生させるように構成され、超音波観測トランスデューサは、中心領域から複数の所定の伝播方向に沿って位置し、複数の所定の伝播方向は、0度及び180度とは異なる角度を成す少なくとも二つの方向を有し、励起プローブ及び超音波観測トランスデューサを異方性柔軟媒体の表面に配置した際に、電子制御デバイスは以下のことを行うように構成される:
(a)中心領域から異方性柔軟媒体内に発散して伝播するのに適した少なくとも一つのせん断波を励起プローブによって発生させて、50ms未満の総観測期間内に、中心領域からの所定の伝播方向におけるせん断波の伝播を超音波観測トランスデューサに観測させること;
(b)超音波観測トランスデューサによって収集されたデータから、複数の所定の伝播方向の各々におけるせん断波の少なくとも一つの伝播パラメータを同時に決定すること;
(c)複数の所定の伝播方向の各々におけるせん断波の少なくとも一つの伝播パラメータから、異方性柔軟媒体の繊維の方向と、繊維に垂直な方向におけるレオロジー弾性パラメータと、繊維の方向におけるレオロジー弾性パラメータとのうちから選択された異方性柔軟媒体の少なくとも一つのレオロジー特性を決定すること。
本発明に係る特性評価デバイスの多様な実施形態においては、任意で、以下の構成を更に適用することができる:
‐ 超音波観測トランスデューサの数は、3個から20個の間、有利には5個から10個の間である;
‐ 電子制御デバイスは、各超音波観測トランスデューサと向き合うせん断波の通過を検出するように構成され、中央ユニットによって決定されるせん断波の伝播パラメータは、せん断波の伝播速度、又はせん断波の飛行時間を表すパラメータである;
‐ 電子制御デバイスは、以下のことを行うことによって、せん断波の通過を検出するのに適している;
‐ 少なくとも毎秒300ショットの超音波圧縮波の率で超音波観測トランスデューサに音響圧縮波を放出させて、
‐ 異方性柔軟媒体によって反射された超音波信号を超音波観測トランスデューサによって感知して、
‐ 超音波観測トランスデューサに向き合うせん断波の通過の際の異方性柔軟媒体の内部変位を決定すること;
励起プローブは、超音波観測トランスデューサを含む略ディスク状の超音波トランスデューサである;
最後に、本発明の更なる対象は、上記デバイス用の超音波プローブのセットであり、中心軸に沿って超音波圧縮波を放出して、放射圧によって中心軸に沿って異方性柔軟媒体を変位させて、中心軸から媒体内に発散するせん断を伝播させるように構成された略ディスク状の超音波励起トランスデューサを含み、その超音波励起トランスデューサは、中心軸に対して発散する異なる複数の伝播方向にそれぞれ分配されていて且つ中心軸からe cm未満の距離に位置する複数の超音波観測トランスデューサを有し、それら伝播方向は、0度及び180度とは異なる角度を成す少なくとも二つの方向を有する。有利には、超音波観測トランスデューサは、中心軸に中心が有る円の上に等分配され、各超音波観測トランスデューサは、中心軸に平行なビームとして圧縮波を放出するように構成される。
本発明の他の特徴及び利点は、添付図面を参照して、非限定的な例として与えられる本発明の一実施形態の以下の説明から明らかになるものである。
本発明の一実施形態に係る特性評価方法の適用を示す概略図である。 図1に見て取れる超音波プローブのセットの平面図である。 図2の超音波プローブのセットが属する特性評価デバイスの概略図である。
異なる図面において、同じ参照符号は同一又は同様の要素を指称する。
図1及び図2に概略的に示されるように、本発明の対象は、人間や動物の体、特に、生体の一部等といった、繊維を含む少なくとも一つの部分3と外面1とを有する異方性柔軟媒体Cを特性評価することである。例えば、外面は、人間又は動物の皮膚であり得て、部分3は、心臓、特に心筋や、骨格筋の一部であり得る。任意で、部分3は、非繊維組織2によって外面1から離隔され得る。
特性評価は、励起プローブ5及び複数の観測プローブ6を含むプローブセット4を用いて行われる。
励起プローブ5は、特に、実質的にディスク状又はカップ状の超音波励起トランスデューサであり得て、励起プローブの使用時に異方性柔軟媒体Cの深さ方向に位置する中心軸Zを有する。任意で、この励起プローブ5は、外面1に当てられる凹面を含み、一般的に、この凹面には、超音波検査で従来使用されているゲル8が充填される。励起プローブ5は、中心軸Zに沿って異方性柔軟媒体Cに超音波圧縮波9を放出するように構成され、その集束点は、中心軸Zに沿って数センチメートル(例えば1から6cm)にわたって広がり且つ中心軸Zに垂直に数ミリメートル(例えば0.2から3mm)の幅を有する中心領域10内に存在する。
中心軸Zに沿った集束点10の位置及び長さは、統計的には、励起プローブが外面1に配置された際に、検査されている人間又は動物の特性評価される繊維部分3内に少なくとも部分的に集束点が見出されるように設定される。従って、人間の心筋の特性評価について、集束点10は、例えば2から4cmの深さzで始まり、また、例えば5から10cmの深さzで終わる。
例えば、励起プローブ5は、3cm未満、有利には2cm未満の半径Rを有し得る。
観測プローブ6の数は、2個以上、有利には3個以上、例えば3個から20個の間、有利には5個から10個の間であり得る。図示されている例では、観測プローブの数は8個である。
観測プローブ6は、中心軸Z周りに分配されていて、それぞれ、中心軸Zに対して発散している異なる伝播方向Pを向いている。各場合において、伝播方向Pは、0度及び180度とは異なる角度を成す少なくとも二つの伝播方向を有する。
観測プローブ6は、中心軸Zから、2cm未満、例えば0.5cm未満の距離に位置し得る。観測プローブ6は、有利には中軸Z上に中心がある円の上に等分配され得る。以下において、発散している伝播方向Pは、中心軸Zに垂直な平面内の直交基準系X、Y内において、X軸に対する角度θによって位置特定される。
観測プローブ6はそれぞれ、医療用イメージングで従来使用されているタイプの一つの超音波観測トランスデューサ、又は、任意で、一組の超音波観測トランスデューサであり得て、中心軸Zに平行な比較的長距離において集束されて、軸Zに平行なビームを形成する超音波圧縮波12をそれぞれ放出するようになっている。
有利には、各観測プローブ6は、単一の超音波観測トランスデューサによって形成される。複数の超音波観測トランスデューサ6は、Z軸に対して放射状(半径方向)に延伸するアレイとして構成され、媒体Cに向けてある程度の凹面を有し、各超音波観測トランスデューサ6の集束点12aが、媒体C内の集束点10と略同じ深さに位置するようにして、また、この集束点12aが、Z軸に対して半径方向に小さな寸法、例えば0.1から2mm程度の寸法を有し、直交方向に数ミリメートル程度の大きな寸法を有するようにする。異なる超音波観測トランスデューサに対応する複数の集束点12aは、好ましくは互いに離隔される。
有利には、観測プローブ6は励起プローブ5に含まれ得て、例えば、励起プローブ5の厚さ方向に設けられた凹部に含まれ得る。
励起プローブ5及び観測プローブ6は、例えばマルチワイヤケーブル7を介して制御デバイスと通信する。図3に示されるように、この制御デバイスは、特定の電子筐体20と通信するコンピュータ19等を備え得て、電子筐体20自体は、ケーブル7を介してプローブセット4に接続される。制御デバイスは、任意で、電子筐体20及びコンピュータ19の機能を全て統合した単一の装置となり得る点に留意されたい。
図示されている例では、電子筐体20は、トランスデューサの数の分だけのチャネル、例えば9つのチャネルを含み得て、各チャネルは、励起プローブ5(T0)と、観測プローブ6(T1〜T8)に接続される。各チャネルは、バッファメモリ15A(B0〜B8)が付随していて且つマイクロプロセッサ等の電子中央ユニット16(CPU)と通信するアナログ・デジタル変換器15(A/D0〜A/D8)を含み得て、電子中央ユニット16自体は、メモリ17(MEM)、信号処理回路18(DSP)、コンピュータ19と通信し得る。電子中央ユニット16は、任意で、心電図装置14(ECG)と更に通信し得る。
上記デバイスは以下のように動作する。
(a)測定ステップ
使用者が媒体Cの繊維部分3を特性評価したい場合、使用者は上述のようにプローブセット4を皮膚1に適用して、測定ステップ(a)を開始する。測定ステップ(a)中においては、電子中央ユニット16がまず、励起プローブ5によって集束された超音波を短期間にわたって放出して、例えば、特許文献1に一般的な原理が説明されているように、波圧の効果によって軸Zに沿った組織の変位を生じさせる。
せん断波11がZ軸に対して実質的に放射状(半径方向)に伝播して(図1を参照)、この伝播中において、電子中央ユニット16は、全ての観測プローブ6に同時に媒体Cを観測させて、せん断波11の伝播を観測する。このため、特許文献2及び特許文献1に特に説明されているように、電子中央ユニット16は、全ての観測プローブ6に超音波圧縮波を高い率(例えば、毎秒300ショット以上)で同時に放出させて、観測プローブ6が組織3によって反射された超音波信号を感知する。そして、それらの信号は、まずバッファメモリ15aに記憶され、その後メモリ17に記憶される。
この測定ステップは、例えば数ミリ秒にわたって続く。
任意で、電子筐体20は、トランスデューサの数よりも少ないチャネルを含み得る。例えば、電子筐体20は、超音波励起トランスデューサ5用のチャネルA/D0‐B0と、スイッチングデバイス(図示せず)を介して複数の超音波観測トランスデューサ6に連続して接続されるチャネルA/D1‐B1とを含み得る。この場合、測定ステップ(a)は、せん断波がそれぞれ連続して複数回放出されて、その伝播を複数の観測プローブ6のうち一つを介して観測するプロセスを含み、そのプロセスが各観測プローブ6について繰り返される。この場合においても、測定ステップは非常に短く、50ms未満である。
(b)計算ステップ
続く計算ステップ(b)は、実時間で又はオフラインで行われ、電子中央ユニット16又はコンピュータ19が、メモリに記憶されている反射信号を利用する。この処理中においては、特許文献2及び特許文献1に特に説明されているように、信号は、各観測プローブ6に向かい合っている組織3の経時的な内部変位又は変形を決定するのに用いられ、例えば、感知信号の相関計算によって決定される。
従って、弾性せん断波の通過を各観測プローブ6の前方で位置特定することができ、また、このことを組織3の各深さにおいて行うことができる。例えば、せん断波の通過を、組織3内の所定の数の深さ方向箇所、例えば10から30箇所の深さ方向箇所に、例えばミリメートル×ミリメートル単位で位置特定し得る。
そして、電子中央ユニット16又はコンピュータ19は、各深さにおいて、励起プローブ5と各観測プローブ6との間のせん断波の飛行時間(TOF,time of flight)t、又は飛行時間を表す他のパラメータ、又はせん断波の伝播を表す他の伝播データを決定して、そこから、せん断波の伝播速度Vが、各観測プローブ6に対応する放射方向Pにおける関連深さにおいて推定される。
(c)特性評価ステップ
そして、繊維の方向と、繊維に垂直な方向におけるレオロジー弾性パラメータ(特に、弾性率)と、繊維の方向に沿ったレオロジー弾性パラメータ(特に、弾性率)とのなかから選択された媒体3の少なくとも一つレオロジー特性を決定することができる。有利には、これらレオロジー特性は同時に決定される。
上記レオロジー特性は、多様な方法、例えば、補間(内挿)によって、より好ましくは、特異値分解によって決定され得る。
1.補間(内挿)
この方法では、各伝播方向Pについて計算ステップ(b)で計算された伝播速度の計算値に応じて、略楕円形の曲線C(V(θ)・cosθ,V(θ)・sinθ)が補間によって決定される。ここで、V(θ)は、XY平面内の伝播パラメータの値である。繊維の方向は、V(θ)の最大値に対応する角度θに対応していて、繊維の方向におけるレオロジー弾性パラメータがV(θ)に従って決定され、繊維に垂直な方向におけるレオロジー弾性パラメータがV(θ+π/2)に従って決定される。例えば、これらのパラメータが繊維に平行な弾性率Eparと、繊維に垂直な弾性率Eperpである場合、これらの弾性率は、以下の式によって決定可能である:
Figure 0006502367
ここで、ρは媒体3の密度である。
2.特異値分解
この方法では、まず、各深さについて、媒体3内の超音波の弾性伝播テンソルが、測定された伝播速度から決定される。この場合、伝播がXY平面内のみにおいて得られるものであるので、このテンソルは階数2(rank2)の行列Mである:
Figure 0006502367
ここで、成分Eは弾性率である。特に、この弾性テンソルは、非特許文献3で定義されている。
計算ステップ(b)においては、上記非特許文献2に記載されているように、以下の式を解くことによって、速度V(θ)からテンソルMを決定することができる:
Figure 0006502367
ここで、ρは媒体3の密度であり、Nはトランスデューサの数であり、θ…θは、観測プローブ6に対応する異なる伝播方向Pの角度である。
上記非特許文献2に記載されているように、テンソルMを特異値分解して、上記パラメータEpar、Eperpを与える対角行列を得ることができる:
Figure 0006502367
この決定は、関連する深さにおけるファイバーの方向θの決定を伴う。何故ならば以下の関係式があるからである:
M=RM
ここで、Rは、角度θに対応する回転行列である:
Figure 0006502367
測定ステップ(a)、計算ステップ(b)、及び特性評価ステップ(c)は、比較的高速で繰り返され得て、心臓又は他の筋肉の動作に従って、例えば、媒体3を形成する筋肉の完全な動作サイクルにわたる媒体3のレオロジー特性、パラメータの連続的な値を与えることができる。心臓の場合、心電計14によって電子中央ユニット16に通信されるデータを用いて、連続的な測定は心拍の周期と同期され得る。そこから、筋肉の収縮に関連する生理パラメータの値、例えば、収縮性、最大硬度や他のパラメータが推定される。
感知プローブが感知する信号は、媒体3を形成する筋肉の動作によるその媒体3の変形又は変位を位置特定することによって、互いに相対的な連続的測定を再調整することも可能にし得る点に留意されたい。
1 繊維を含む部分
2 非繊維組織
3 外面(皮膚)
4 プローブセット
5 励起プローブ
6 観測プローブ
7 ケーブル
8 ゲル
9 超音波圧縮波
10 中心領域
11 せん断波
12 超音波圧縮波

Claims (21)

  1. 繊維を含む少なくとも一つの部分(3)と外面(1)とを備える異方性柔軟媒体(C)を特性評価するための方法であって、
    測定ステップ(a)と、少なくとも一つの計算ステップ(b)と、特性評価ステップ(c)とを備え、
    前記測定ステップ(a)において、前記異方性柔軟媒体の中心領域(10)から発散して伝播する少なくとも一つのせん断波(11)を発生させて、前記異方性柔軟媒体の表面(1)から複数の超音波観測トランスデューサ(6)を用いて、前記複数の超音波観測トランスデューサを固定したままで、前記中心領域(10)からの複数の所定の伝播方向(P)における前記少なくとも一つのせん断波の伝播を観測し、前記複数の所定の伝播方向(P)が、0度及び180度とは異なる角度を成す少なくとも二つの方向を有し、前記超音波観測トランスデューサ(6)が少なくとも前記所定の伝播方向(P)に沿って配置され、前記測定ステップが50ms未満の期間内に行われ、
    前記計算ステップ(b)において、前記複数の所定の伝播方向(P)の各々について、前記測定ステップ(a)で収集されたデータから、前記せん断波の少なくとも一つの伝播パラメータを決定し、
    前記特性評価ステップ(c)において、前記計算ステップ(b)で前記複数の伝播方向(P)の各々について決定された前記せん断波の少なくとも一つの伝播パラメータから、前記異方性柔軟媒体の繊維の方向と、前記繊維に垂直な方向におけるレオロジー弾性パラメータと、前記繊維の方向におけるレオロジー弾性パラメータとのうちから選択された前記異方性柔軟媒体の少なくとも一つのレオロジー特性を決定する、方法。
  2. 前記特性評価ステップ(c)で決定されるレオロジー弾性パラメータが弾性率である、請求項1に記載の方法。
  3. 前記測定ステップ(a)において、前記異方性柔軟媒体内の特定の深さ範囲にわたって前記せん断波を発生させて、前記少なくとも一つのせん断波の伝播を前記深さ範囲内の複数の深さにおいて観測し、
    前記計算ステップ(b)において、前記複数の所定の伝播方向の各々について、前記複数の深さにおける前記せん断波の少なくとも一つの伝播パラメータを決定し、
    前記特性評価ステップ(c)において、前記複数の深さにおける前記少なくとも一つのレオロジー特性を決定する、請求項1又は2に記載の方法。
  4. 前記測定ステップ(a)において、前記せん断波の伝播が観測される所定の伝播方向(P)の数が3個から20個までの間である、請求項1から3のいずれか一項に記載の方法。
  5. 前記測定ステップ(a)において、前記複数の超音波観測トランスデューサ(6)を全て同時に用いて、前記少なくとも一つのせん断波の伝播を前記複数の所定の伝播方向(P)の全てにおいて同時に観測する、請求項1から4のいずれか一項に記載の方法。
  6. 前記測定ステップ(a)において、複数のせん断波を連続的に放出して、前記複数の超音波観測トランスデューサ(6)の一部を用いて、前記複数の所定の伝播方向(P)の少なくとも一つにおいて各せん断波の伝播を連続的に観測する、請求項1から4のいずれか一項に記載の方法。
  7. 前記測定ステップ(a)において、前記所定の伝播方向(P)に沿ってそれぞれ位置する複数の所定の測定点におけるせん断波の通過を検出し、
    前記計算ステップ(b)において決定される前記せん断波の伝播パラメータが、前記せん断波の伝播速度、又は前記中心領域(10)からの前記せん断波の飛行時間のいずれかを表すパラメータである、請求項1から6のいずれか一項に記載の方法。
  8. 前記異方性柔軟媒体の表面(1)に垂直な中心軸(Z)から前記せん断波を発生させて、前記せん断波の通過が検出される前記複数の測定点(6)の各々が、前記中心軸から2cm未満の距離に位置する、請求項7に記載の方法。
  9. 前記測定ステップ(a)において、
    前記超音波観測トランスデューサ(6)による少なくとも毎秒300ショットの超音波圧縮波の率での音響圧縮波(12)の放出によって、前記所定の測定点における前記せん断波の通過を検出し、
    前記異方性柔軟媒体によって反射された超音波信号を前記超音波観測トランスデューサ(6)によって感知し、
    前記所定の測定点における前記せん断波の通過の際の前記異方性柔軟媒体の内部変位を決定する、請求項7又は8に記載の方法。
  10. 前記測定ステップ(a)において、測定点毎に単一の超音波観測トランスデューサ(6)を使用する、請求項7から9のいずれか一項に記載の方法。
  11. 前記測定ステップ(a)において、前記超音波観測トランスデューサ(6)を有する超音波励起トランスデューサ(5)が、前記中心領域(10)に集束された超音波励起波(9)を放出して、中心軸(Z)に沿って前記異方性柔軟媒体を変位させることによって、前記せん断波を発生させる、請求項1から10のいずれか一項に記載の方法。
  12. 前記特性評価ステップ(c)において、複数の所定の伝播方向(P)における伝播パラメータの値である特異値に対応する弾性テンソルを決定して、前記弾性テンソルを特異値分解することによって、前記少なくとも一つのレオロジー特性を決定する、請求項1から11のいずれか一項に記載の方法。
  13. 前記測定ステップ(a)において、中心軸(Z)から前記せん断波を発生させ、前記特性評価ステップ(c)において、前記計算ステップ(b)で計算された前記所定の伝播方向(P)における伝播パラメータの値に従った補間によって、略楕円形の曲線C(V(θ)・cosθ,V(θ)・sinθ)を決定し、V(θ)が前記中心軸(Z)に垂直な平面内の伝播パラメータの値であり、θが、前記平面に属する基準系に対する前記平面内の伝播方向を示す角度であり、V(θ)の最大値に対応する角度θに対応する前記繊維の方向と、V(θ)に従う前記繊維の方向におけるレオロジー弾性パラメータと、V(θ+π/2)に従う前記繊維に垂直な方向におけるレオロジー弾性パラメータとを決定する、請求項1から11のいずれか一項に記載の方法。
  14. 前記異方性柔軟媒体が、動いている人間又は動物の筋肉の少なくとも一部分を備え、前記方法が前記測定ステップ(a)、前記計算ステップ(b)及び前記特性評価ステップ(c)を複数回連続して行い、前記筋肉の収縮に関係する生理パラメータを推定する、請求項1から13のいずれか一項に記載の方法。
  15. 繊維を含む少なくとも一つの部分と外面とを有する異方性柔軟媒体を特性評価するためのデバイスであって、励起プローブ(5)及び複数の超音波観測トランスデューサ(6)を制御する電子制御デバイス(19、20)を備え、前記励起プローブ(5)が、中心領域(10)から前記異方性柔軟媒体にせん断波を発生させるように構成されていて、前記複数の超音波観測トランスデューサ(6)が、前記中心領域(10)から複数の所定の伝播方向(P)に沿って位置し、前記複数の所定の伝播方向(P)が、0度及び180度とは異なる角度を成す少なくとも二つの方向を有し、
    前記異方性柔軟媒体の表面に前記励起プローブ(5)及び前記超音波観測トランスデューサ(6)が配置された際に、前記電子制御デバイス(19、20)が、
    (a)前記中心領域(10)から前記異方性柔軟媒体内に発散しながら伝播する少なくとも一つのせん断波を前記励起プローブ(5)に発生させ、50ms未満の総観測期間内に、前記中心領域(10)からの前記所定の伝播方向(P)における前記せん断波の伝播を前記超音波観測トランスデューサ(6)に観測させ、
    (b)前記超音波観測トランスデューサ(6)が収集したデータから、前記複数の所定の伝播方向(P)の各々について前記せん断波の少なくとも一つの伝播パラメータを同時に決定し、
    (c)前記複数の所定の伝播方向(P)の各々についての前記せん断波の少なくとも一つの伝播パラメータから、前記異方性柔軟媒体の繊維の方向と、前記繊維に垂直な方向におけるレオロジー弾性パラメータと、前記繊維の方向におけるレオロジー弾性パラメータとのうちから選択された前記異方性柔軟媒体の少なくとも一つのレオロジー特性を決定するように構成されている、デバイス。
  16. 前記超音波観測トランスデューサの数が3個から20個までの間である、請求項15に記載のデバイス。
  17. 前記電子制御デバイス(19、20)が、各超音波観測トランスデューサ(6)に向き合う前記せん断波の通過を検出するように構成されていて、中央ユニットによって決定される前記せん断波の伝播パラメータが、前記せん断波の伝播速度、又は前記せん断波の飛行時間のいずれかを表すパラメータである、請求項16に記載のデバイス。
  18. 前記電子制御デバイス(19、20)が、
    少なくとも毎秒300ショットの超音波圧縮波の率で音響圧縮波を前記超音波観測トランスデューサ(6)に放出させ、
    前記異方性柔軟媒体によって反射された超音波信号を前記超音波観測トランスデューサ(6)に観測させ、
    前記超音波観測トランスデューサ(6)に向き合う前記せん断波の通過の際の前記異方性柔軟媒体の内部変位を決定することによって、
    前記せん断波の通過を検出するように構成されている、請求項17に記載のデバイス。
  19. 前記励起プローブ(5)が、前記超音波観測トランスデューサ(6)を含む略ディスク状の超音波トランスデューサである、請求項18に記載のデバイス。
  20. 請求項15から19のいずれか一項に記載のデバイス用の超音波プローブセットであって、中心軸に沿って超音波圧縮波を放出して、放射圧によって前記中心軸に沿って前記異方性柔軟媒体を変位させ、前記中心軸から発散するせん断波の伝播を前記異方性柔軟媒体内に生じさせるように構成された略ディスク状の超音波励起トランスデューサ(5)を含み、前記超音波励起トランスデューサが、前記中心軸(Z)に対して発散する複数の伝播方向にそれぞれ分配され且つ前記中心軸から2cm未満の距離に配置される複数の超音波観測トランスデューサ(6)を含み、各超音波観測トランスデューサが、前記中心軸(Z)に平行なビーム(12)として圧縮波を放出するように構成されていて、前記伝播方向が、0度及び180度とは異なる角度を成す少なくとも二つの方向を有する、超音波プローブセット。
  21. 前記複数の超音波観測トランスデューサ(6)が、前記中心軸(Z)に中心がある円の上に等分配されている、請求項20に記載の超音波プローブセット。
JP2016549242A 2014-01-31 2015-01-12 異方性柔軟媒体を特性評価するための方法及び超音波デバイス、並びにそのような特性評価デバイス用の超音波プローブのセット Active JP6502367B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1400265A FR3017041B1 (fr) 2014-01-31 2014-01-31 Procede et dispositif ultrasonore de caracterisation des milieux mous anisotropes, et ensemble de sonde ultrasonore pour un tel dispositif de caracterisation
FR1400265 2014-01-31
PCT/FR2015/050058 WO2015114232A1 (fr) 2014-01-31 2015-01-12 Procede et dispositif ultrasonore de caracterisation des milieux mous anisotropes, et ensemble de sonde ultrasonore pour un tel dispositif de caracterisation

Publications (2)

Publication Number Publication Date
JP2017504435A JP2017504435A (ja) 2017-02-09
JP6502367B2 true JP6502367B2 (ja) 2019-04-17

Family

ID=50543210

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016549242A Active JP6502367B2 (ja) 2014-01-31 2015-01-12 異方性柔軟媒体を特性評価するための方法及び超音波デバイス、並びにそのような特性評価デバイス用の超音波プローブのセット

Country Status (8)

Country Link
US (1) US11399804B2 (ja)
EP (1) EP3099240B1 (ja)
JP (1) JP6502367B2 (ja)
KR (1) KR102369563B1 (ja)
CN (1) CN106232013B (ja)
CA (1) CA2939013C (ja)
FR (1) FR3017041B1 (ja)
WO (1) WO2015114232A1 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR112017013861B1 (pt) 2014-12-24 2022-07-12 Supersonic Imagine Método para elastrografia de onda de cisalhamento para geração de imagem de um campo de observação em um meio anisotrópico e aparelho de geração de imagem para a implementação do método
JP6259953B2 (ja) * 2015-08-03 2018-01-10 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. せん断波を用いる測定に関する超音波システム及び方法
JP7133480B2 (ja) * 2016-06-30 2022-09-08 ザ・ユニバーシティ・オブ・ブリティッシュ・コロンビア 腹部の超音波せん断波振動エラストグラフィ
CN107440740B (zh) * 2017-07-21 2021-06-25 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
WO2020043561A1 (en) * 2018-08-29 2020-03-05 Koninklijke Philips N.V. Ultrasound system and methods for smart shear wave elastography
US20220015741A1 (en) * 2019-01-09 2022-01-20 Koninklijke Philips N.V. Ultrasound system and method for shear wave characterization of anisotropic tissue
WO2023198719A1 (en) 2022-04-11 2023-10-19 Institut National de la Santé et de la Recherche Médicale Method and apparatus for ultrasonic evaluation of an isolated organ

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2791136B1 (fr) 1999-03-15 2001-06-08 Mathias Fink Procede et dispositif d'imagerie utilisant les ondes de cisaillement
FR2844058B1 (fr) 2002-09-02 2004-11-12 Centre Nat Rech Scient Procede et dispositif d'imagerie utilisant des ondes de cisaillement
FR2844178B1 (fr) 2002-09-06 2005-09-09 Dispositif et procede pour la mesure de l'elasticite d'un organe humain ou animal et l'etablissement d'une representation a deux ou trois dimensions de cette elasticite
WO2008139245A1 (en) * 2007-05-16 2008-11-20 Super Sonic Imagine Method and device for measuring a mean value of visco-elasticity of a region of interest
US20100228126A1 (en) * 2009-03-06 2010-09-09 Mirabilis Medica Inc. Ultrasound treatment and imaging applicator
EP2453800B1 (en) * 2009-07-17 2018-11-21 Koninklijke Philips N.V. Spatially-fine shear wave dispersion ultrasound vibrometry sampling
US10172527B2 (en) * 2009-07-31 2019-01-08 Supersonic Imagine Method and apparatus for measuring a physical parameter in mammal soft tissues by propagating shear waves
JP6148010B2 (ja) 2009-11-25 2017-06-14 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 焦束されたスキャンラインビーム形成での超音波剪断波撮像
WO2011138722A1 (en) * 2010-05-03 2011-11-10 Andrey Rybyanets Resonantly amplified shear waves
GB201009006D0 (en) * 2010-05-28 2010-07-14 Ntnu Technology Transfer As Ultrasound acquisition
US20130317361A1 (en) * 2011-02-04 2013-11-28 Hitachi Medical Corporation Ultrasound diagnostic apparatus and method
WO2013116809A1 (en) * 2012-02-03 2013-08-08 Los Alamos National Security, Llc Ultrasound waveform tomography with tv regularization
EP2816958B1 (en) * 2012-02-21 2020-03-25 Maui Imaging, Inc. Determining material stiffness using multiple aperture ultrasound
BR112014026753A2 (pt) * 2012-04-27 2017-06-27 Echosens dispositivo para a medida de parâmetro ultrassono ou biomecânico de meio viscoelástico
US10324065B2 (en) * 2014-01-06 2019-06-18 Samsung Electronics Co., Ltd. Ultrasound diagnostic apparatus, ultrasound image capturing method, and computer-readable recording medium

Also Published As

Publication number Publication date
JP2017504435A (ja) 2017-02-09
CA2939013C (fr) 2023-03-14
CA2939013A1 (fr) 2015-08-06
EP3099240A1 (fr) 2016-12-07
CN106232013B (zh) 2020-08-14
WO2015114232A1 (fr) 2015-08-06
EP3099240B1 (fr) 2024-02-28
US11399804B2 (en) 2022-08-02
CN106232013A (zh) 2016-12-14
CN111772678A (zh) 2020-10-16
US20160345938A1 (en) 2016-12-01
KR20160135167A (ko) 2016-11-25
FR3017041B1 (fr) 2016-03-04
FR3017041A1 (fr) 2015-08-07
KR102369563B1 (ko) 2022-03-03

Similar Documents

Publication Publication Date Title
JP6502367B2 (ja) 異方性柔軟媒体を特性評価するための方法及び超音波デバイス、並びにそのような特性評価デバイス用の超音波プローブのセット
JP6294548B2 (ja) 心臓壁の電気双極子密度の幾何学的測定のための装置
JP6838057B2 (ja) 種々のイメージング・モダリティを用いて得られる画像を位置合せし、画像レジストレーションを検証するためのシステムおよび方法
CN106999162B (zh) 超声波诊断装置以及弹性评价方法
JP5594361B2 (ja) トランスデューサのアレイを備えた装置により患者の心筋収縮性を測定する方法及び患者の心筋収縮性を測定するための装置
JP2006280532A (ja) 超音波診断装置
US10893805B2 (en) Method of and apparatus for characterizing spatial-temporal dynamics of media excitable for deformation
US20150099960A1 (en) Ultrasonic probe and medical apparatus including the same
US20140187940A1 (en) Method of calculating displacement of shear wave, method of calculating mechanical modulus of body, and system using the methods
JP2018520733A5 (ja)
US11490876B2 (en) Ultrasonic diagnostic device and method for evaluating physical properties of biological tissue
WO2020042020A1 (zh) 一种超声弹性检测设备及剪切波弹性成像方法、装置
JP6321162B2 (ja) 繊維状媒体のマッピングの方法及びデバイス
JP6933262B2 (ja) 超音波診断システム及び画像生成方法
KR102608323B1 (ko) 보강적인 전단파 초음파 이미징을 위한 방법, 시스템 및 컴퓨터 프로그램 제품
Urban et al. Investigation of the effects of myocardial anisotropy for shear wave elastography using acoustic radiation force and harmonic vibration
CN111772678B (en) Ultrasonic method and apparatus for characterizing weak anisotropic soft media, and ultrasonic probe assembly for such characterization apparatus
JP5491671B2 (ja) 関心領域の粘弾性の平均値を測定するための方法
JP2022543248A (ja) 心臓用装置
Zhao et al. Robust shear wave speed measurement using comb-push ultrasound radiation force
Jiao et al. A shear wave endoscopic elasticity imaging approach with micro focused piezoelectric transducer

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20171225

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180824

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180903

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20181203

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190225

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190320

R150 Certificate of patent or registration of utility model

Ref document number: 6502367

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313115

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250