JP2022542005A - 超音波を使用して不均一な媒体を非侵襲的に特性化するための方法及びシステム - Google Patents
超音波を使用して不均一な媒体を非侵襲的に特性化するための方法及びシステム Download PDFInfo
- Publication number
- JP2022542005A JP2022542005A JP2022500603A JP2022500603A JP2022542005A JP 2022542005 A JP2022542005 A JP 2022542005A JP 2022500603 A JP2022500603 A JP 2022500603A JP 2022500603 A JP2022500603 A JP 2022500603A JP 2022542005 A JP2022542005 A JP 2022542005A
- Authority
- JP
- Japan
- Prior art keywords
- point
- medium
- input
- focusing
- calculated
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 94
- 238000002604 ultrasonography Methods 0.000 title claims abstract description 78
- 230000004044 response Effects 0.000 claims abstract description 91
- 239000011159 matrix material Substances 0.000 claims abstract description 67
- 238000004364 calculation method Methods 0.000 claims description 40
- 230000005540 biological transmission Effects 0.000 claims description 30
- 238000005259 measurement Methods 0.000 claims description 29
- 238000012512 characterization method Methods 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 4
- 101710084218 Master replication protein Proteins 0.000 claims description 3
- 101710112083 Para-Rep C1 Proteins 0.000 claims description 3
- 101710112078 Para-Rep C2 Proteins 0.000 claims description 3
- 102100022881 Rab proteins geranylgeranyltransferase component A 1 Human genes 0.000 claims description 3
- 102100022880 Rab proteins geranylgeranyltransferase component A 2 Human genes 0.000 claims description 3
- 101710119887 Trans-acting factor B Proteins 0.000 claims description 3
- 101710119961 Trans-acting factor C Proteins 0.000 claims description 3
- 230000000875 corresponding effect Effects 0.000 description 31
- 239000010410 layer Substances 0.000 description 17
- 230000008569 process Effects 0.000 description 14
- 239000000523 sample Substances 0.000 description 14
- 238000012360 testing method Methods 0.000 description 13
- 238000002592 echocardiography Methods 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 8
- 238000012285 ultrasound imaging Methods 0.000 description 7
- 230000004075 alteration Effects 0.000 description 6
- 210000004185 liver Anatomy 0.000 description 6
- 230000001934 delay Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000003491 array Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000002123 temporal effect Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 235000013372 meat Nutrition 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 210000000577 adipose tissue Anatomy 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- TVEXGJYMHHTVKP-UHFFFAOYSA-N 6-oxabicyclo[3.2.1]oct-3-en-7-one Chemical compound C1C2C(=O)OC1C=CC2 TVEXGJYMHHTVKP-UHFFFAOYSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 230000001594 aberrant effect Effects 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000012623 in vivo measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/043—Analysing solids in the interior, e.g. by shear waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/26—Arrangements for orientation or scanning by relative movement of the head and the sensor
- G01N29/262—Arrangements for orientation or scanning by relative movement of the head and the sensor by electronic orientation or focusing, e.g. with phased arrays
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02483—Other human or animal parts, e.g. bones
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02491—Materials with nonlinear acoustic properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/106—Number of transducers one or more transducer arrays
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Gynecology & Obstetrics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
- トランスデューサによって生成された不均一な媒体の領域に一連の入射超音波を生成するステップであって、一連の入射超音波は送信基底(i)を構成する、ステップと;
- 入力送信基底(i)と出力受信基底(u)との間で定義された実験的な反射行列Rui(t)を記録するステップと;
- 第1の点(P1)を中心にした入力フォーカルスポットに対応する実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第2の点(P2)を中心にした出力フォーカルスポットに対応する実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、間の媒体の応答REP(r,Δr)を決定するステップであって、前記応答は、媒体中の空間位置rの中心点(PC)の関数として表され、前記中心点(PC)は、第1の点と第2の点と(P1、P2)の中間に位置し、第1及び第2の点を通過する測定軸(AXm)の原点であり、前記測定軸は、媒体の第1の軸(X)に対して角度βを形成し、第1の点(P1)は、測定軸上の距離座標+Δrにあり、第2の点(P2)は、測定軸上の距離座標-Δrにある、ステップと;
を含む。
Ninは、送信基底iの要素数であり。
Noutは、出力受信基底uの要素数であり、
Rui(t)は実験的反射行列であり、Rui(uout,iin,τ(rin,rout,uout, iin))は、時間τでの送信iinに続いてトランスデューサuoutによって記録された実験的反射行列Rui(t)の要素であり、
τは、送信基底iと第1の点P1との間の超音波の外向き飛行時間τinと、第2の点と受信基底uとの間の超音波の戻り飛行時間τoutと、の合計であり、次の式によって表れる。
τ(rin,rout,uout,iin)=τin(rin,iin)+τout(rout,uout)
- 距離座標Δrの複数の値と同じ中心点(PC)とに対して、かつ、所定の角度βに対応している同じ測定軸(AXm)に対して、計算された複数の応答REP(r,Δr)である応答プロファイルPR(δr)を算出するステップであって、δrは中心点から第2の点の距離である、すなわち、Δr=δr.uβになる値であり、uβは角度βで定義される測定軸AXmの方向の単位ベクトルである、ステップ、
をさらに含む。
- 応答プロファイルの係数に基づいて中心点(PC)の解像度w(r)を算出するステップであって、ここで、解像度w(r)は、ゼロ距離座標(|Δr|=0)の辺りを中心とする、応答プロファイルPR(δr)の前記ピークの幅である、ステップ、
をさらに含む。
-解像度w(r)と理論解像度w0(r)に基づいて中心点(PC)のフォーカシング基準F(r)を決定するステップであって、理論解像度w0(r)は、入力送信基底(i)と出力受信基底(u)とに基づいて算出される、ステップ、
をさらに以下を含む。
-フォーカシング基準F(r)が媒体の複数の点で計算され、入力フォーカシング及び/または出力フォーカシングのための少なくとも1つの計算パラメータが、前記複数の点に対するフォーカシング基準F(r)の平均を最小化または最大化することによって最適化される画像最適化ステップ、
をさらに以下を含む。
- パルス(ω1)、送信基底(i)、及び、受信基底(u)について、中心点(PC)で第1の分析計算、ここで、理論解像度は、トランスデューサアレイを中心点(PC)から見た場合のトランスデューサアレイが離れる角度によって計算される。
- パルス範囲(Δω)、送信基底(i)、及び、受信基底(u)について、中心点(PC)で第2の分析計算、ここで、理論解像度は、実験的反射行列Rui(t)からの信号の周波数スペクトルによって重み付けされた中心点(PC)からトランスデューサアレイを見た角度の前記パルス範囲に亘った積分計算である、及び
- 先ず、入力仮想トランスデューサ(TVin)に対応する媒体の第1の点と、送信基底(i)と、の間、次に、出力仮想トランスデューサ(TVout)に対応する媒体の第2の点と、受信基底(u)と、の間の波動伝播シミュレーションの第3の計算であって、前記シミュレーションは、応答REP(r,Δr)と媒体内の波動伝播のモデルを使用している。
- 2つの相互応答間の平均相関係数である対称性のレベルα(r)を算出するステップであって、前記平均は、所定の解像度wd(r)より大きい係数の距離座標値について計算され、及び/または、角度値の範囲βについて、または、所定の角度値βdについて計算され、
第1の点(P1)を中心とした入力フォーカルスポットに対応する実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第2の点(P2)を中心とした出力フォーカルスポットに対応する実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の媒体の第1の応答REP1(r,Δr)と、
第2の点(P2)を中心とした入力フォーカルスポットに対応する実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第1の点(P1)を中心とした出力フォーカルスポットに対応する実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の媒体の第2の応答REP2(r,-Δr)。
Re[.]は実数部の数学演算子であり、
|.|は係数数学演算子であり、
<.>は平均数学演算子であり、
*は複素共役演算子である。
- 応答REP(r,Δr)の係数の2乗の平均であるアフォーカル強度Ioff(r)を算出するステップであって、平均は、所定の解像度wd(r)よりも大きい係数の距離値に対して計算される、及び/または、角度値の範囲βまたは所定の角度値βdについて計算され、例えば、以下の式によって計算される。
IM(r)=α(r)・Ioff(r)
である多重散乱IM(r)を算出するステップ、
- 1-対称性のレベルα(r)とアフォーカル強度Ioff(r)との積、すなわち、
IN(r)=(1-α(r))・Ioff(r)
であるノイズ強度IN(r)を算出するステップ、
をさらに含み、これにより以下の関係式を有する。
Ioff(r)=IM(r)+IN(r)
- ゼロ距離座標値(|Δr|=0)の場合、すなわち、第1の点(P1)、第2の点(P2)、及び、中心点(PC)が一致する媒体の点の場合の応答REP(r,Δr=0)の2乗係数の値である共焦点強度Ion(r)を算出するステップと、
- 次の式に基づいて計算された単一散乱強度IS(r)を算出するステップと、
Ion(r)=IS(r)+2IM(r)+IN(r)
をさらに含む。
- 媒体の局所特性パラメータの画像を決定するステップであって、前記局所特性パラメータは、応答REP(r,Δr)に基づいて決定される、ステップ、
をさらに含む。
- 不均一な媒体の領域で一連の入射超音波を生成するのに適したトランスデューサの第1のアレイ(10)であって、この一連の入射超音波は送信基底(i)を構成し、前記領域によって後方散乱された超音波の時間の関数として記録するのに適している、トランスデューサの第1のアレイ(10)と、
- トランスデューサの第1のアレイに結合されたコンピューティングユニット(42)であって、
- 入力送信基底(i)と出力受信基底(u)との間で定義された実験的反射行列Rui(t)を記録することと、
- 第1の点を中心とした入力フォーカルスポットに対応する実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第2の点(P2)を中心とした出力フォーカルスポットに対応する実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の媒体のREP(r,Δr)を算出することであって、前記応答は、媒体中の空間位置rの中心点(PC)の関数として表され、前記中心点(PC)は、第1の点と第2の点と(P1、P2)の中間に位置し、第1及び第2の点を通過する測定軸(AXm)の原点であり、前記測定軸は、媒体の第1の軸(X)に対して角度βを形成し、第1の点(P1)は、測定軸上の距離座標+Δrにあり、第2の点(P2)は、測定軸上の距離座標-Δrにある、ことと、
に適しているコンピューティングユニット(42)と、
を備える。
本明細書では、不均一なサンプルの非侵襲的超音波特性評価の方法とシステムについて説明する。これらの方法及びシステムは、図4に示される定義に基づく。空間位置rの中心点PCは、媒体の空間参照系で定義され、第1の点P1と第2の点P2との間の中間に位置する。測定軸AXmが定義され、第1の点P1及び第2の点P2を通過し、トランスデューサ11のアレイの第1の軸Xと角度βを形成する。中心点PCは、測定軸の原点AXm(測定軸上のゼロの距離座標)に位置する。第1の点P1は距離座標-Δrにあり、2第2の点P2は測定軸の原点である中心点PCから距離座標+Δrにある。
- トランスデューサ11のアレイ10を用いて、不均一な媒体の領域に一連の入射超音波USinを生成するステップであって、一連の入射超音波は、送信基底iである、ステップと、
- 入力送信基底iと出力受信基底uとの間で定義された実験的反射行列Rui(t)を記録するステップと、
- 第1の点P1の空間位置rinの入力仮想トランスデューサTVinと第2の点P2の空間位置routの出力仮想トランスデューサTVoutの間の媒体の応答REP(r,Δr)を決定するステップであって、応答は、空間位置rの中心点PCの関数として表されている、ステップと、
を含む。
r=(rout+rin)/2と、Δr=(rout-rin)/2
- 第1の点P1を中心にした入力フォーカルスポットを生成する実験的反射行列Rui(t)に基づく入力フォーカシングプロセスであって、入力フォーカルスポットは、入力仮想トランスデューサTVinに対応する、入力フォーカシングプロセスと、
- 第2の点P2を中心にした出力フォーカルスポットを生成する実験的反射行列Rui(t)に基づく出力フォーカシングプロセスであって、出力フォーカルスポットは、出力仮想トランスデューサTVoutに対応する、出力フォーカシングプロセスと、
が適用される。
Ninは、送信基準iの要素数であり、
Noutは、出力受信基底uの要素数であり、
Rui(t)は実験的反射行列であり、Rui(uout,iin,τ(rin,rout,uout,iin))は、時間τでの送信iinに続くトランスデューサuoutによって記録された実験的反射行列Rui(t)の要素である。
τ(rin,rout,uout,iin)=τin(rin,iin)+τout(rout,uout) (式2)
上記の方法で決定された媒体の応答を得た後、応答プロファイルPR(δr)を決定するステップを実行することができ、応答プロファイルは、複数の距離座標Δrの値について計算された複数の応答REP(r,Δr)である。この応答プロファイルPR(δr)は、空間位置rの同じ中心点PCに対して、同じ角度βの方向に対応する同じ測定軸AXmに沿って考慮される。従って、応答プロファイルPR(δr)は、複数の距離δrに対して決定され、距離δrは、中心点PCに対する第2の点P2の横座標、すなわちΔr=δr.uβとなる値であり、uβは角度βによって定義される、測定軸AXmの方向の単位ベクトルである。言い換えると、応答プロファイルPR(δr)は、図5の画像の垂直スライスであり、この応答プロファイルは、画像スライスの曲線である。
- 幅がフォーカルスポットに関連する単一散乱に対応する値δrが低い場合の過剰強度(曲線の最大値)、及び
- 多重散乱現象とノイズから生じるエコーが原因の、δrのすべての値に存在するインコヒーレントな背景、
を示す。
- 図7のパートA(左側)では、「ファントム」と呼ばれる試験媒体(媒体A)の画像を見ることができ、これは、異なる剛性の2つのシリンダを含む。
- 図7のパートB(中央)では、上記の試験媒体(媒体B)上に肉の層が配置された媒体の画像を見ることができる。
- 図7のパートC(右側)では、「イン・ビボ」肝臓である媒体(媒体C)の画像を見ることができる。
次に、応答プロファイルの係数に基づいて、角度βの測定軸AXmの方向の中心点PCの解像度w(r)を決定するステップを実行できる。従って、この解像度は、超音波画像の解像度の局所的な推定値である。
- 角度βがπ/2の場合、すなわち、角度β=π/2に対応する応答プロファイルに対して取得された値w(r)による軸方向の解像度、及び
- ゼロ角度の場合、すなわち、すなわち、角度β=0に対応する応答プロファイルに対して取得された値w(r)による横方向の解像度、
を取得することを可能にする。
第1の変形によれば、理論的解像度w0(r)は、パルス(ω1)、送信基底(i)、及び受信基底(u)について中心点(PC)での第1の分析計算によって算出される。それはトランスデューサアレイを中心点(PC)から見た角度で計算される。それは、空間位置rの中心点を超音波処理するために送信中に、またはこの中心点から来るエコーを収集するために受信中に使用される開口の最大半角に依存する。
解像度w(r)及び理論的解像度w0(r)に基づいて、中心点PCのフォーカシング基準F(r)を算出するステップを実行することができる。理論上の解像度は、例えば、入力送信基底i、出力受信基底u、及び媒体内の超音波の伝播のモデリングに基づいて算出される。
F(r)=w(r)/w0(r) (式3)
または
F(r)=w0(r)/w(r) (式4)
次いで、対称性レベルα(r)を決定するステップを実行でき、この対称性レベルは、2つの相互応答間、すなわち、入力仮想トランスデューサ(TVin)及び出力仮想トランスデューサ(TVout)の位置rin及びroutを交換することによる平均相関係数である。従って、空間位置rの中心点PCでの媒体の音響相互関係の基準がテストされる。多重散乱現象から生じる信号は相互的であるが、電子ノイズから生じる信号は相関していないことに留意されたい。
第1の点P1を中心とした入力フォーカルスポットを生成する実験的反射行列の入力フォーカルスポットに基づいて計算された入力仮想トランスデューサTVin(空間位置rin)と、第2の点P2を中心とした出力フォーカルスポットを生成する実験的反射行列の出力フォーカルスポットに基づいて計算された出力仮想トランスデューサTVout(空間位置rout)と、の間の媒体の第1の応答REP1(r,Δr)、及び、
第2の点P2を中心とした入力フォーカルスポットを生成する実験的反射行列の入力フォーカルスポットに基づいて計算された入力仮想トランスデューサTVin(空間位置rin)と、第1の点P1を中心とした出力フォーカルスポットを生成する実験的反射行列の出力フォーカルスポットに基づいて計算された出力仮想トランスデューサTVout(空間位置rout)と、の間の媒体の第2の応答REP2(r,-Δr)、
が算出される。
Re[.]は実数部の数学演算子であり、
|.|は係数数学演算子であり、
<.>は平均数学演算子であり、この演算子は、1つまたは複数の変数(上記では、例えば、所定の解像度および角度値βより大きい距離座標値)に従って適用され得、
*は複素共役演算子である。
次いで、中心点PCの第1の多重散乱インジケータε(r)を算出するステップを実行することができ、この多重散乱インジケータは、例えば、以下の式によって計算される。
α(r)は、媒体の中心点PCで定義される対称性レベルである。
次いで、中心点PCのアフォーカル強度Ioff(r)を算出するステップを実行することが可能であり、このアフォーカル強度は、応答REP(r,Δr)の二乗係数の平均であり、平均が計算される。所定の解像度wd(r)(すなわち|Δr|>wd(r))より大きい係数の距離値の場合、及び/または、角度値βの範囲または所定の角度値βdに対して計算された場合、例えば次の式で計算される。
I(r,Δr)=|REP(r,Δr)|2 (式9)
ここで、
I(r,Δr)は点の強度、すなわちこの点での応答の2乗係数であり、
|.|は関数数学演算子であり、
<.>は平均数学演算子であり、この演算子は、1つまたは複数の変数(上記では、例えば、所定の解像度および角度値βより大きい距離座標値)に従って適用され得る。
次いで、IM(r)を算出するステップを実行でき、この多重散乱強度は、対称性レベルα(r)とアフォーカル強度Ioff(r)の積である。すなわち、
IM(r)=α(r)・Ioff(r) (式10)
次いで、ノイズ強度IN(r)を算出するステップを実行できる。このノイズ強度は、1から対称レベルα(r)を引いたものとアフォーカル強度Ioff(r)の積である。すなわち、
IN(r)=(1-α(r))・Ioff(r) (式11)
これにより、以下の関係を有する。
Ioff(r)=IM(r)+IN(r) (式12)
次いで、共焦点強度Ion(r)を算出するステップを実行できる。この共焦点強度は、距離座標値Δrがゼロ、すなわち|Δr|=0、の場合の応答REP(r,Δr=0)の2乗係数の値になり、第1の点P1、第2の点P2、及び中心点PCが一致する媒体内の点を意味する。
単一散乱強度IS(r)を算出するステップを実行できる。この単一散乱強度は、次の式に基づいて計算される。
Ion(r)=IS(r)+2IM(r)+IN(r) (式14)
上記の計算の1つの可能な利用法は、超音波画像の計算を最適化することである。
a)最適な音速c1optの曲線B1の単層モデル
b)最適な音速c2optの曲線B2の2層モデル
c)最適な音速c3optの曲線B1の3層モデル
11 トランスデューサ
20 媒体
21 散乱体
40 システム
41 プロービングデバイス
42 コンピューティングユニット
43 ディスプレイデバイス
Claims (20)
- 不均一な媒体の非侵襲的超音波特性評価のための方法であって、当該方法は、
- トランスデューサによって生成された、不均一な媒体の領域に一連の入射超音波を生成するステップであって、前記一連の入射超音波は、送信基底(i)を構成する、ステップと;
- 入力送信基底(i)と出力受信基底(u)の間で定義された実験的反射行列Rui(t)を記録するステップと;
- 第1の点を中心とした入力フォーカルスポットに対応する前記実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第2の点(P2)を中心とした出力フォーカルスポットに対応する前記実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の前記媒体の応答REP(r,Δr)を算出するステップであって、前記応答は、前記媒体中の空間位置rの中心点(PC)の関数として表され、前記中心点(PC)は、前記第1の点(P1)と前記第2の点(P2)との中間に位置し、前記第1及び第2の点を通過する測定軸AXmの原点であり、前記測定軸は、前記媒体の第1の軸(X)に対して角度βを形成し、前記第1の点(P1)は、前記測定軸上の距離座標+Δrにあり、前記第2の点(P2)は、前記測定軸上の距離座標-Δrにある、ステップと、
含むことを特徴とする方法。 - 前記応答REP(r,Δr)を算出するステップにおいて、前記実験的反射行列の前記入力フォーカシングは、前記送信基底と前記入力仮想トランスデューサとの間の波の外向きの飛行時間を使用し、前記出力フォーカシングは、前記出力仮想トランスデューサと前記受信基底の間の波の戻り飛行時間を使用することを特徴とする請求項1に記載の方法。
- 前記媒体の前記応答REP(r,Δr)は、以下の式によって計算され、
Ninは、送信基底iの要素数であり。
Noutは、出力受信基底uの要素数であり、
Rui(t)は実験的反射行列であり、Rui(uout,iin,τ(rin,rout,uout, iin))は、時間τで送信iinに続いてトランスデューサuoutによって記録された前記実験的反射行列Rui(t)の要素であり、
τは、前記送信基底iと前記第1の点(P1)との間の超音波の外向き飛行時間τinと、前記第2の点(P2)と前記受信基底uとの間の超音波の戻り飛行時間τoutと、の合計であり、
τ(rin,rout,uout,iin)=τin(rin,iin)+τout(rout,uout)
の式によって表れることを特徴とする請求項1または2に記載の方法。 - - 距離座標Δrの複数の値と同じ中心点(PC)とに対して、かつ、所定の角度βに対応している同じ測定軸(AXm)に対して、計算された複数の応答REP(r,Δr)である応答プロファイルPR(δr)を算出するステップであって、δrは中心点から前記第2の点の距離である、すなわち、Δr=δr.uβになる値であり、uβは角度βで定義される測定軸AXmの方向の単位ベクトルである、ステップ、
をさらに含むことを特徴とする請求項1~3のいずれか一項に記載の方法。 - - 前記応答プロファイルPR(δr)の係数に基づいて前記中心点(PC)の解像度w(r)を算出するステップであって、ここで、解像度w(r)は、ゼロ距離座標(|Δr|=0)の辺りを中心とする、前記応答プロファイルPR(δr)のピークの幅である、ステップ、
をさらに含むことを特徴とする請求項4に記載の方法。 - 前記ピークの幅は、前記ピークの最大高さの一部である高さで推定され、前記一部は、例えば、前記ピークの前記最大高さの半分であることを特徴とする請求項5に記載の方法。
- - 解像度w(r)及び理論解像度w0(r)に基づいて前記中心点(PC)のフォーカシング基準F(r)を算出するステップであって、前記理論解像度w0(r)は、前記入力送信基底(i)と前記出力受信基底(u)とに基づいて算出される、ステップ、
をさらに含むことを特徴とする請求項4に記載の方法。 - - 前記フォーカシング基準F(r)が前記媒体の複数の点で計算され、前記入力フォーカシング及び/または前記出力フォーカシングのための少なくとも1つの計算パラメータが、前記複数の点に対する前記フォーカシング基準F(r)の平均値を最小化または最大化することによって最適化される画像最適化ステップ、
をさらに含む請求項7に記載の方法。 - 前記少なくとも1つの計算パラメータは、前記媒体中の音速を含むことを特徴とする請求項8に記載の方法。
- 理論解像度は、
- パルス(ω1)、前記送信基底(i)、及び、前記受信基底(u)について、前記中心点(PC)における第1の分析計算であって、前記理論解像度は、トランスデューサアレイを前記中心点(PC)から見た角度によって計算される、第1の分析計算と、
- パルス範囲(Δω)、前記送信基底(i)、及び、前記受信基底(u)について、前記中心点(PC)における第2の分析計算であって、前記理論解像度は、前記実験的反射行列Rui(t)からの信号の周波数スペクトルによって重み付けされた前記中心点(PC)からトランスデューサアレイを見た角度の前記パルス範囲に亘った積分計算である、第2の分析計算と、
- 先ず、前記入力仮想トランスデューサ(TVin)に対応する媒体の前記第1の点と、前記送信基底(i)と、の間、次に、前記出力仮想トランスデューサ(TVout)に対応する媒体の前記第2の点と、前記受信基底(u)と、の間の波動伝播シミュレーションの第3の計算であって、前記シミュレーションは、前記応答REP(r,Δr)と前記媒体内の波動伝播のモデルを使用している、第3の計算と、
のリストに含まれる技術によって算出されることを特徴とする請求項7に記載の方法。 - - 2つの相互応答間の平均相関係数である対称性のレベルα(r)を算出するステップであって、平均値は、所定の解像度wd(r)より大きい係数の距離座標値について計算され、及び/または、角度値の範囲βについて、または、所定の角度値βdについて計算され、
第1の点(P1)を中心とした入力フォーカルスポットに対応する前記実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第2の点(P2)を中心とした出力フォーカルスポットに対応する前記実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の前記媒体の第1の応答REP1(r,Δr)と、
第2の点(P2)を中心とした入力フォーカルスポットに対応する実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第1の点(P1)を中心とした出力フォーカルスポットに対応する前記実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の前記媒体の第2の応答REP2(r,-Δr)、であるステップ、
をさらに備えることを特徴とする請求項1~3のいずれか一項に記載の方法。 - - 前記応答REP(r,Δr)の係数の2乗の平均であるアフォーカル強度Ioff(r)を算出するステップであって、平均は、所定の解像度wd(r)よりも大きい係数の距離値に対して計算される、及び/または、角度値の範囲βまたは所定の角度値βdについて計算され、例えば、
- 対称性のレベルα(r)とアフォーカル強度Ioff(r)との積、すなわち、
IM(r)=α(r)・Ioff(r)
である多重散乱IM(r)を算出するステップと、
- 1-対称性のレベルα(r)とアフォーカル強度Ioff(r)との積、すなわち、
IN(r)=(1-α(r))・Ioff(r)
であるノイズ強度IN(r)を算出ステップと、
をさらに含み、これにより、
Ioff(r)=IM(r)+IN(r)
の関係式を有する、ステップと、
をさらに含むことを特徴とする請求項11に記載の方法。 - - ゼロ距離座標値(|Δr|=0)の場合、すなわち、第1の点(P1)、第2の点(P2)、及び、中心点(PC)が一致する媒体の点の場合の応答REP(r,Δr=0)の2乗係数の値である共焦点強度Ion(r)を算出するステップと、
- 次の式
Ion(r)=IS(r)+2IM(r)+IN(r)
に基づいて計算された単一散乱強度IS(r)を算出するステップと、
をさらに含むことを特徴とする請求項14に記載の方法。 - - 前記媒体の局所特性パラメータの画像を決定するステップであって、前記局所特性パラメータは、応答REP(r,Δr)に基づいて決定される、ステップ、
をさらに含むことを特徴とする請求項1~17のいずれか一項に記載方法。 - 前記局所特性パラメータは、解像度w(r)、フォーカシング基準F(r)、対称性のレベルα(r)、第1の多重散乱インジケータε(r)、第2の多重散乱インジケータγ(r)、アフォーカル強度Ioff(r)、共焦点強度Ion(r)、多重散乱強度IM(r)、単一散乱強度IS(r)、及び、ノイズ強度IN(r)を含むリストから選択されることを特徴とする請求項18に記載の方法。
- 不均一な媒体の非侵襲的超音波特性評価のためのシステム(40)であって、当該システムは、
- 前記不均一な媒体の領域で一連の入射超音波を生成するのに適したトランスデューサの第1のアレイ(10)であって、前記一連の入射超音波は送信基底(i)を構成し、前記領域によって後方散乱された超音波の時間の関数として記録するのに適している、トランスデューサの第1のアレイ(10)と、
- 前記トランスデューサの前記第1のアレイに結合されたコンピューティングユニット(42)であって、
- 前記入力送信基底(i)と出力受信基底(u)との間で定義された実験的反射行列Rui(t)を記録することと、
- 第1の点を中心とした入力フォーカルスポットに対応する前記実験的反射行列の入力フォーカシングに基づいて計算された空間位置rinの入力仮想トランスデューサ(TVin)と、第2の点(P2)を中心とした出力フォーカルスポットに対応する前記実験的反射行列の出力フォーカシングに基づいて計算された空間位置routの出力仮想トランスデューサ(TVout)と、の間の前記媒体のREP(r,Δr)を算出することであって、前記応答は、前記媒体中の空間位置rの中心点(PC)の関数として表され、前記中心点(PC)は、第1の点と第2の点と(P1、P2)の中間に位置し、第1及び第2の点を通過する測定軸(AXm)の原点であり、前記測定軸は、媒体の第1の軸(X)に対して角度βを形成し、第1の点(P1)は、測定軸上の距離座標+Δrにあり、第2の点(P2)は、測定軸上の距離座標-Δrにある、ことと、
に適しているコンピューティングユニット(42)と、
を備えることを特徴とするシステム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1908904A FR3099580A1 (fr) | 2019-08-02 | 2019-08-02 | Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène |
FR1908904 | 2019-08-02 | ||
PCT/FR2020/051416 WO2021023933A1 (fr) | 2019-08-02 | 2020-07-31 | Procédé et système de caractérisation ultrasonore non invasive d'un milieu hétérogène |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022542005A true JP2022542005A (ja) | 2022-09-29 |
Family
ID=68807030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022500603A Pending JP2022542005A (ja) | 2019-08-02 | 2020-07-31 | 超音波を使用して不均一な媒体を非侵襲的に特性化するための方法及びシステム |
Country Status (8)
Country | Link |
---|---|
US (1) | US20240032889A1 (ja) |
EP (1) | EP4007914A1 (ja) |
JP (1) | JP2022542005A (ja) |
KR (1) | KR20220038678A (ja) |
CN (1) | CN114144118A (ja) |
AU (1) | AU2020324570A1 (ja) |
FR (1) | FR3099580A1 (ja) |
WO (1) | WO2021023933A1 (ja) |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0256481B1 (de) * | 1986-08-20 | 1994-01-12 | Siemens Aktiengesellschaft | Verfahren und Einrichtung zur adaptiven Fokussierung bei einem medizinischen Ultraschall-Bildgabegerät |
US5551433A (en) * | 1994-08-05 | 1996-09-03 | Acuson Corporation | Method and apparatus for a geometric aberration transform in an adaptive focusing ultrasound beamformer system |
FR2815717B1 (fr) * | 2000-10-20 | 2003-01-10 | Centre Nat Rech Scient | Procede et dispositif non invasif de focalisation d'ondes acoustiques |
WO2003017843A1 (en) * | 2001-08-24 | 2003-03-06 | The Brigham And Women's Hospital, Inc. | Ultrasound therapy |
US9117439B2 (en) * | 2008-03-13 | 2015-08-25 | Supersonic Imagine | Method and apparatus for ultrasound synthetic imagining |
FR2932339B1 (fr) | 2008-06-09 | 2012-11-23 | Centre Nat Rech Scient | Procede et dispositif de sondage par propagation d'ondes |
DE102009025463A1 (de) * | 2009-06-12 | 2011-03-10 | Technische Universität Dresden | Anordnung und Verfahren zur Bestimmung von Schichtdicken und Schallgeschwindigkeiten in Medien mit Hilfe von Ultraschall |
ITGE20090070A1 (it) * | 2009-08-31 | 2011-03-01 | Esaote Spa | Metodo e dispositivo per il rilevamento e la visualizzazione di informazioni emodinamiche in particolare del flusso ematico nelle vene, mediante ultrasoni |
CN103969651A (zh) * | 2014-04-24 | 2014-08-06 | 中国科学院声学研究所 | 自适应声学成像方法 |
-
2019
- 2019-08-02 FR FR1908904A patent/FR3099580A1/fr active Pending
-
2020
- 2020-07-31 JP JP2022500603A patent/JP2022542005A/ja active Pending
- 2020-07-31 US US17/631,929 patent/US20240032889A1/en active Pending
- 2020-07-31 EP EP20820244.0A patent/EP4007914A1/fr active Pending
- 2020-07-31 KR KR1020227002139A patent/KR20220038678A/ko active Search and Examination
- 2020-07-31 WO PCT/FR2020/051416 patent/WO2021023933A1/fr active Application Filing
- 2020-07-31 AU AU2020324570A patent/AU2020324570A1/en active Pending
- 2020-07-31 CN CN202080052379.4A patent/CN114144118A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
CN114144118A (zh) | 2022-03-04 |
AU2020324570A1 (en) | 2022-01-27 |
KR20220038678A (ko) | 2022-03-29 |
WO2021023933A1 (fr) | 2021-02-11 |
FR3099580A1 (fr) | 2021-02-05 |
US20240032889A1 (en) | 2024-02-01 |
EP4007914A1 (fr) | 2022-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11346819B2 (en) | Methods and systems for non-invasively characterizing a heterogeneous medium using ultrasound | |
US11776526B2 (en) | Method and system for ultrasonic characterization of a medium | |
US11761928B2 (en) | Method and system for ultrasonic characterization of a medium | |
Bottenus et al. | Acoustic reciprocity of spatial coherence in ultrasound imaging | |
US11768181B2 (en) | Method and system for ultrasonic characterization of a medium | |
Lambert et al. | Ultrasound matrix imaging—Part I: The focused reflection matrix, the F-factor and the role of multiple scattering | |
US20220082693A1 (en) | Method and system for ultrasonic characterization of a medium | |
US20240036004A1 (en) | Method and system for ultrasonic characterization of a medium | |
Szasz | Advanced beamforming techniques in ultrasound imaging and the associated inverse problems | |
JP2022542005A (ja) | 超音波を使用して不均一な媒体を非侵襲的に特性化するための方法及びシステム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230614 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240213 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240315 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240401 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240516 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240603 |