JP6698511B2 - Ultrasonic signal processing device and ultrasonic transmitting/receiving device using the same - Google Patents
Ultrasonic signal processing device and ultrasonic transmitting/receiving device using the same Download PDFInfo
- Publication number
- JP6698511B2 JP6698511B2 JP2016245761A JP2016245761A JP6698511B2 JP 6698511 B2 JP6698511 B2 JP 6698511B2 JP 2016245761 A JP2016245761 A JP 2016245761A JP 2016245761 A JP2016245761 A JP 2016245761A JP 6698511 B2 JP6698511 B2 JP 6698511B2
- Authority
- JP
- Japan
- Prior art keywords
- wave
- velocity
- signal processing
- ultrasonic signal
- wave component
- 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
Links
- 238000009826 distribution Methods 0.000 claims description 76
- 238000005259 measurement Methods 0.000 claims description 60
- 238000006073 displacement reaction Methods 0.000 claims description 53
- 238000000034 method Methods 0.000 claims description 34
- 230000005540 biological transmission Effects 0.000 claims description 32
- 238000007689 inspection Methods 0.000 claims description 31
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000002123 temporal effect Effects 0.000 claims description 20
- 230000001902 propagating effect Effects 0.000 claims description 17
- 239000000523 sample Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000005855 radiation Effects 0.000 description 8
- 210000004204 blood vessel Anatomy 0.000 description 7
- 238000013016 damping Methods 0.000 description 7
- 206010028980 Neoplasm Diseases 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000010297 mechanical methods and process Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 201000007270 liver cancer Diseases 0.000 description 1
- 208000019423 liver disease Diseases 0.000 description 1
- 208000014018 liver neoplasm Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
本発明は、超音波信号処理装置、および、それを用いた超音波送受信装置に係り、被検体内においてせん断波を発生させ、その伝搬速度を計測することにより、生体組織の性状を評価する技術に関する。 The present invention relates to an ultrasonic signal processing device and an ultrasonic transmitting/receiving device using the same, and a technique for evaluating the properties of a biological tissue by generating a shear wave in a subject and measuring its propagation velocity. Regarding
超音波やMRI（Magnetic Resonance Imaging）、X線CT（Computed Tomography）に代表される医療用の画像表示装置は、目視できない生体内の情報を数値または画像の形態で提示する装置として広く利用されている。中でも超音波を利用して画像を表示する超音波撮像装置は、他の装置と比較して高い時間分解能を備えており、例えば拍動下の心臓を滲みなく画像化できる性能を持つ。 Medical image display devices represented by ultrasound, MRI (Magnetic Resonance Imaging), and X-ray CT (Computed Tomography) are widely used as devices that present invisible information in the body in the form of numerical values or images. There is. Above all, an ultrasonic imaging apparatus that displays an image using ultrasonic waves has a higher time resolution than other apparatuses, and has an ability to image a pulsating heart without blurring.
被検体である生体内を伝搬する波は、主に縦波と横波に区別される。超音波撮像装置の製品に搭載されている組織形状を映像化する技術や、血流速度を計測する技術は、主に縦波（音速約1540m/s）の情報を利用している。 The waves propagating in the living body, which is the subject, are mainly classified into longitudinal waves and transverse waves. The technology for visualizing the tissue shape and the technology for measuring blood flow velocity, which are installed in the products of the ultrasonic imaging device, mainly use the information of the longitudinal wave (sound velocity about 1540 m/s).
近年、生体内を伝搬する横波（以降、せん断波）を利用して組織の弾性率を評価する技術が注目されており、慢性肝疾患や癌に対する臨床利用が進められている。この技術では、計測対象となる組織内部にせん断波を発生させ、その伝搬速度から弾性率等の弾性を表す評価指標を算出する。せん断波を発生させる手法は、機械方式と放射圧方式に大別される。機械方式は、バイブレータ等を利用して体表面に1kHz程度の振動を与えてせん断波を発生させる方式で、振動源となる駆動装置が必要である。一方、放射圧方式は、超音波を組織内の局所に集中させる集束超音波を利用して生体内に音響放射圧を加え、瞬時的に発生する組織変位を利用してせん断波を発生させる。いずれの方式も、発生したせん断波による組織変位を超音波により計測することにより、伝搬速度を算出する。せん断波の伝搬速度から、組織性状を表す弾性率等の特性値を計算によって求める。 In recent years, a technique for evaluating the elastic modulus of a tissue using a transverse wave (hereinafter, shear wave) propagating in a living body has attracted attention, and clinical use for chronic liver disease and cancer has been advanced. In this technique, a shear wave is generated inside the tissue to be measured, and an evaluation index representing elasticity such as elastic modulus is calculated from its propagation velocity. Methods of generating shear waves are roughly classified into mechanical methods and radiation pressure methods. The mechanical method is a method of generating a shear wave by applying a vibration of about 1 kHz to the body surface using a vibrator or the like, and requires a driving device as a vibration source. On the other hand, in the radiation pressure method, acoustic radiation pressure is applied to the inside of a living body by using focused ultrasonic waves that concentrate ultrasonic waves locally in the tissue, and shear waves are generated by using tissue displacement that occurs instantaneously. In either method, the propagation velocity is calculated by measuring the tissue displacement due to the generated shear wave with ultrasonic waves. From the propagation velocity of the shear wave, the characteristic values such as elastic modulus representing the tissue property are calculated.
非特許文献１には、検査対象を伝搬するせん断波の波面振幅の分布の計測データを、フーリエ変換した後、せん断波の伝搬する方向により、いくつかのデータに分割し、速度マップを生成する方法が開示されている。具体的には、非特許文献１の技術では、伝搬するせん断波の波面振幅の２次元空間（ｘ−ｚ面）における分布を、フーリエ変換により２次元の空間周波数（ｋｘ−ｋｚ面）における強度分布に変換する。この２次元空間周波数（ｋｘ−ｋｚ面）の強度分布に、２次元空間周波数（ｋｘ−ｋｚ面）の中心から所定の角度範囲を取り出すフィルタを適用することにより、８つに等分した角度範囲ごとの空間周波数（ｋｘ−ｋｚ面）の強度分布に分割する。分割した角度範囲ごとの空間周波数（ｋｘ−ｋｚ面）の強度分布を、それぞれ実空間（ｘ−ｚ面）の波面振幅の分布に再変換する。そして、各角度範囲の振幅の分布についてせん断波の速度の分布を算出し、最後に、角度範囲ごとのせん断波の速度分布を加算することにより全体の速度分布の画像を生成する。
In
特許文献１には、検査対象の機械的特性を計測する目的でせん断波を励起し、多方向に伝搬する波面を検出し、最低波速度を判定し、判定結果に基づいて媒体の波速度や材料特性を算出する方法が開示されている。また、異なる方向に伝搬する波成分を算出し、各伝搬方向についての波速度マップを生成する方法が開示されている。特許文献１の技術では、せん断波の２次元空間（ｘ−ｚ面）における分布を、２つの空間周波数（ｋｘ、ｋｚ）と時間周波数（ｆ）の空間（ｋ−ｆ空間）の分布に変換し、波の速度ｃをｃ＝ｆ／ｋにより算出している。
In
せん断波を利用して組織の弾性を評価する方法は、弾性を定量的に計測できるため、腫瘍診断において極めて重要であり、臨床価値が高い。せん断波を利用して組織の弾性を精度よく求めるためには、伝搬速度を正確に計測する必要がある。例えば、予め伝搬方向を仮定して２以上の計測点を設け、両計測点の通過に要する時間を正確に計測することが望まれる。しかし、せん断波の伝搬経路に検査対象以外の構造物が存在する場合、または背景組織が線維や脂肪などの複雑な組織構造を持つ場合には、せん断波は、反射、屈折、回折などの影響により伝搬方向が変化する。伝搬方向の変化は、伝搬速度の計測の主たる誤差要因となっている。 The method of evaluating elasticity of tissue using shear waves is extremely important in tumor diagnosis and has high clinical value because elasticity can be quantitatively measured. In order to obtain the elasticity of tissue accurately using shear waves, it is necessary to measure the propagation velocity accurately. For example, it is desirable to provide two or more measurement points assuming the propagation direction in advance and to accurately measure the time required for passing through both measurement points. However, when there is a structure other than the inspection target in the propagation path of the shear wave, or when the background tissue has a complicated tissue structure such as fiber or fat, the shear wave is affected by reflection, refraction, diffraction, etc. Causes the propagation direction to change. The change in the propagation direction is the main error factor in measuring the propagation velocity.
非特許文献１に記載の技術では、角度範囲ごとのせん断波の速度分布を最終的には加算するため、すべての方向のせん断波および反射波が含まれた画像が生成される。このため、反射波の影響を低減することはできるが、屈折波や回折波の影響は生成した画像に含まれている。一方、特許文献１に記載の技術は、最低波速度を求めたり、異なる２以上の伝搬方向についての波速度マップを生成することはできるが、伝搬方向が変化した屈折波や回折波の影響を低減して速度を計測することはできない。
In the technique described in Non-Patent
本発明の目的は、せん断波の速度を精度よく計測することのできる超音波信号処理装置を提供することにある。 An object of the present invention is to provide an ultrasonic signal processing device capable of accurately measuring the velocity of shear waves.
上記目的を達成するために、本発明によれば、検査対象に超音波を送受信して得られた超音波信号を受け取って、検査対象を伝搬するせん断波による検査対象の複数の計測点の変位を時系列に計測する計測部と、計測部の計測結果に基づいて、変位の周波数分布を求め、周波数ごとに、所定の強度以上の波成分を選択する選択部と、選択部が選択した波成分の速度を求める速度算出部とを有する超音波信号処理装置が提供される。 In order to achieve the above object, according to the present invention, an ultrasonic signal obtained by transmitting and receiving ultrasonic waves to and from an inspection target is received, and displacements of a plurality of measurement points of the inspection target due to shear waves propagating through the inspection target are received. A time-series measurement unit, a frequency distribution of displacement is obtained based on the measurement result of the measurement unit, a selection unit that selects a wave component having a predetermined intensity or more for each frequency, and a wave selected by the selection unit. There is provided an ultrasonic signal processing device having a speed calculation unit for calculating the speed of a component.
本発明によれば、せん断波の速度を精度よく計測することができる。 According to the present invention, the velocity of shear waves can be accurately measured.
以下、本発明の実施形態を図面に従い説明する。 Embodiments of the present invention will be described below with reference to the drawings.
図１に実施形態の超音波送受信装置（超音波撮像装置）の一構成例のブロック図を示す。本実施形態の超音波送受信装置は、図１のように、超音波信号処理装置（制御部３０）として、計測部３１と選択部３２と速度算出部３３とを備えている。計測部３１は、超音波信号を用いて、検査対象１００を伝搬するせん断波による検査対象の複数の計測点の変位を時系列に計測する。選択部３２は、計測部３１の計測結果に基づいて、上記変位の周波数分布を求め、さらに、周波数ごとに、所定の強度（振幅）以上の波成分を選択する。速度算出部３３は、選択部３２が選択した波成分の速度を求める。
FIG. 1 shows a block diagram of an example of the configuration of an ultrasonic transceiver (ultrasonic imaging device) according to an embodiment. As shown in FIG. 1, the ultrasonic transmission/reception device of the present embodiment includes a
なお、ここでいう波成分は、複数の計測点の変位の時間変化で表される波（すなわちせん断波）を構成する波成分であって、時間周波数と空間周波数とその強度（振幅）の組み合わせによって定義される。 The wave component mentioned here is a wave component that constitutes a wave (that is, a shear wave) represented by the time change of displacement of a plurality of measurement points, and is a combination of time frequency, spatial frequency, and its intensity (amplitude). Defined by
このように、本実施形態では、強度（振幅）が所定値以上の波成分を選択部３２が選択することにより、せん断波が検査対象の組織によって屈折や回折して生じた波成分を除去し、本来のせん断波である主成分を選択することができる。よって、選択した波成分の速度を速度算出部３３により求めることにより、屈折波や回折波の影響を低減し、せん断波の速度を精度よく求めることができる。せん断波の速度を用いて、弾性等の組織性状を表す特性値を求めることができる。
As described above, in the present embodiment, the
選択部３２は、周波数ごとに所定の強度（振幅）以上の波成分を選択するために、例えば次のような処理を行う。まず、計測部３１が計測した複数の計測点の変位の時間変化を、時間周波数と空間周波数を２軸とする空間の波成分の強度分布に変換する。例えば、選択部３２は、時間と複数の計測点の配置方向の距離とを２軸とする空間（平面）における変位の振幅分布を、２次元フーリエ変換することにより、時間周波数と空間周波数を２軸とする空間（平面）における波成分の強度分布に変換する。つぎに、選択部３２は、時間周波数と空間周波数を２軸とする空間における波成分の強度分布から、時間周波数ごとまたは空間周波数ごとに所定の強度以上の波成分を選択する。これにより、強度が比較的小さい屈折波や回折波を除去し、強度が大きい波成分、すなわち主成分のせん断波の波成分を選択することができる。
The
波成分の選択方法としては、どのような方法を用いてもよいが、例えば、選択部３２は、時間周波数と空間周波数を２軸とする空間の所定の範囲を選択するフィルタを生成し、生成したフィルタを時間周波数と空間周波数を２軸とする空間の波成分の強度分布に適用することにより、所定の強度以上の波成分を選択する構成とする。フィルタの選択する範囲は、時間周波数と空間周波数を２軸とする空間の波成分の強度分布における、最大強度の波成分を含む範囲を選択するように設定することが望ましい。これにより、時間周波数と空間周波数を２軸とする空間にフィルタを適用して所定の範囲を選択することで、所定の強度以上の波成分を選択できる。
Although any method may be used as a method of selecting the wave component, for example, the selecting
なお、フィルタの選択範囲を、操作者から受け付ける入力部１３をさらに備える構成としてもよい。これにより、操作者の所望する速度以上の波成分を選択することができる。
Note that the filter selection range may further include the
なお、選択部３２は、波成分の伝搬方向が所定の方向範囲の波成分を選択する機能を備えていることが望ましい。これにより、反射波を除去することができる。
The
一方、速度算出部３３は、選択部３２が選択した波成分の速度を求めるために、例えば以下のような処理を行う。すなわち、速度算出部３３は、選択部３２が選択した波成分の、時間周波数と空間周波数を２軸とする空間における強度分布を、複数の計測点の変位の時間変化に再変換する。具体的には例えば、時間周波数と空間周波数を２軸とする空間（平面）における、選択した波成分の強度分布を、２次元逆フーリエ変換することにより、時間と距離とを２軸とする空間（平面）における複数の計測点の変位の振幅分布に再変換することができる。速度算出部３３は、再変換より得た、時間と距離とを２軸とする空間（平面）における複数の計測点の変位の振幅分布に、例えば任意の２点を設定し、波成分が２点を通過するのに要する時間と、２点の距離とに基づいて、その２点間の波の速度を求めることができる。
On the other hand, the
また、本実施形態の超音波送受信装置は、検査対象１００の複数の点について、速度算出部３３が求めた速度の分布を示す速度分布画像を生成する画像生成部３４をさらに有する構成であることが望ましい。
Further, the ultrasonic transmission/reception device of the present embodiment is configured to further include an
なお、本実施形態では、検査対象１００内でせん断波を発生させ、伝搬させる方法としては、検査対象１００に振動を与える方法であればどのようなものでもよく、音響放射圧や加振装置などの機械的な加振方法の他、心拍動など検査対象自身による加振も含む。
In the present embodiment, as a method of generating and propagating a shear wave in the
以下、本実施形態の超音波送受信装置の構成についてさらに詳しく説明する。以下の説明では、集束超音波を検査対象１００に照射して、音響放射圧によりせん断波を発生させる場合を例に説明する。
Hereinafter, the configuration of the ultrasonic transmission/reception device of this embodiment will be described in more detail. In the following description, a case will be described as an example in which focused ultrasonic waves are applied to the
図１に示したように、本実施形態の超音波送受信装置１には、送受信制御部２０と、制御部３０が配置されている。超音波送受信装置１には、探触子１０と、入力部３と、表示部１５が接続されている。送受信制御部２０は、探触子１０を構成する各振動子に受け渡す送信信号を生成する送信ビームフォーマ２１と、探触子１０の各振動子の出力から、検査対象１００内の所定の点についての受信信号を生成する受信ビームフォーマ２２とを備えている。
As shown in FIG. 1, a transmission/
制御部３０は、上述した計測部３１、選択部３２、速度算出部３３および画像生成部３４を含む。制御部３０には、メモリ１６が接続されている。
The
制御部３０の計測部３１、選択部３２、速度算出部３３および画像生成部３４の機能は、ソフトウエアによって実現することも可能であるし、その一部または全部をハードウエアによって実現することも可能である。ソフトウエアによって実現する場合、制御部３０を、ＣＰＵ（Central Processing Unit）やＧＰＵ（Graphics Processing Unit）等のプロセッサにより構成し、メモリ１６に予め格納されたプログラムを読み込んで実行することにより、計測部３１、選択部３２、速度算出部３３および画像生成部３４の機能を実現する。また、ハードウエアによって実現する場合には、ＡＳＩＣ（Application Specific Integrated Circuit）のようなカスタムＩＣや、ＦＰＧＡ（Field-Programmable Gate Array）のようなプログラマブルＩＣを用い、計測部３１、選択部３２、速度算出部３３および画像生成部３４の動作を少なくとも実現するように回路設計を行えばよい。
The functions of the
以下、各部の動作を図２、図３を用いて具体的に説明する。ここでは、制御部３０は、ソフトウエアにより実現する場合を例に説明する。図２は、装置全体の動作を示すフローチャートであり、図３（ａ）、（ｂ）は、探触子１０において送受信される超音波を示す説明図である。
The operation of each unit will be specifically described below with reference to FIGS. 2 and 3. Here, the case where the
まず、制御部３０は、図３（ａ）、（ｂ）のように、送受信制御部２０に検査対象１００のＲＯＩ（Ｒｅｇｉｏｎ ｏｆ Ｉｎｔｅｒｅｓｔ）３００の位置を指示し、第１超音波２３を探触子１０から送信させて、検査対象１００内で音響放射圧を生成することにより、検査対象１００のＲＯＩ３００に局所的に圧力を加える（図２のステップ１０１）。これにより、検査対象１００にせん断波が発生する。具体的には、送受信制御部２０は、音響放射圧を生じさせる第１超音波２３をＲＯＩ３００内の所定の位置に照射するための送信信号を送信ビームフォーマ２１に生成させる。送信ビームフォーマ２１が生成した送信信号は、探触子１０を構成する振動子にそれぞれ受け渡され、探触子１０は、所定の音響強度を有し、所定の深さの焦点に収束する第１超音波２３を所定の送信方向に送信する（図２のステップ１０１）。
First, as shown in FIGS. 3A and 3B, the
第１超音波２３が照射された対象物１００には、第１超音波２３による音響放射圧が加わり、第１超音波２３の照射を停止すると、圧力負荷が排除されて復元力が働き、せん断波が伝搬する。せん断波は、第１超音波２３が照射された位置を起点に放射状に伝搬する。図３（ａ）、（ｂ）の図面では、右方向に伝搬する波面を計測する場合を例示している。
Acoustic radiation pressure by the first
制御部３０は、送受信制御部２０に、第２超音波２４を送信するように指示して、第２超音波２４を探触子１０から送信させ（ステップ１０２）、検査対象１００のＲＯＩ３００内のせん断波が伝搬する方向（例えばｘ方向）に等間隔に設定した複数の計測点３０１の組織の変位を計測する（ステップ１０３）。具体的には、制御部３０の制御下で送受信制御部２０は、送信ビームフォーマ２１に送信信号を生成させ、図３（ａ），（ｂ）のように、探触子１０から複数の計測点３０１に所定のタイミングで１以上の第２超音波２４を照射させる（ステップ１０２）。第２超音波２４を照射された計測点３０１から反射等した超音波は、探触子１０の振動子によって受信される。送受信制御部２０は、第２超音波２４の複数の計測点３０１をそれぞれ通って深さ方向（ｚ方向）に延びる複数の受信走査線を設定し、受信ビームフォーマ２２を動作させて、受信走査線上の複数の点について振動子の出力信号を整相加算させて受信信号（ＲＦ（Ｒａｄｉｏ Ｆｒｅｑｕｅｎｃｙ）信号）を生成させる。したがって、ＲＦ信号は、計測点３０１の受信走査線上の異なる深さｚの点における受信信号が、受信走査線方向に連結された信号である。送受信制御部２０は、第２超音波２４の送信および反射等した超音波の受信を、時間経過とともに、所定の時間間隔で繰り返し、時間経過ごとのＲＦ信号を複数の計測点３０１についてそれぞれ生成する。（ステップ１０３）。
The
なお、図３（ａ）の例は、複数の計測点３０１ごとに第２超音波２４を照射し、その都度、送信した計測点３０１に受信走査線を設定してＲＦ信号を得る構成を示している。一方、図３（ｂ）の例は、複数の計測点３０１にわたって１つの第２超音波２４を照射し、複数の計測点３０１にそれぞれ受信走査線を設定してＲＦ信号を得る構成を示している。図３（ａ）の送受信パターンは、感度に優れるが、何回も第２超音波２４を送信する必要があるため時間分解能はあまり高くできない。一方、図３（ｂ）の送受信パターンは、１回の第２超音波２４で複数本の受信走査線について受信信号を得ることができるため、時間分解能が高いが、第２超音波２４の照射範囲が広いため、感度はあまり高くできない。検査対象１００の部位や、必要な時間分解能に応じて、図３（ａ），（ｂ）のいずれかの送受信パターンを用いればよい。
Note that the example of FIG. 3A shows a configuration in which the second
つぎに、制御部３０は、計測部３１に指示して、複数の計測点３０１およびその深さ方向（ｚ方向）の各点について変位を計測させる（ステップ１０３）。図４（ａ）、図５（ａ）、図６は、計測部３１が求めたせん断波の伝搬方向（ｘ方向）と深さ方向（ｚ方向）とを２軸とする空間（平面）における変位の分布を示す説明図である。計測部３１は、受信ビームフォーマ２２からＲＦ信号を受け取り、同一の計測点３０１について時系列に得たＲＦ信号同士の相互相関演算により、複数の計測点３０１およびその深さ方向（ｚ方向）の各点について組織変位を求める。これにより、図４（ａ）、図５（ａ）、図６に示すように、せん断波の伝搬方向(ｘ方向、距離方向）に設定した複数の計測点３０１およびその深さ方向（ｚ方向）の各点（図４（ａ）の画素）についての変位（せん断波の振幅）の時間変化が得られる。図４（ａ）、図５（ａ）、図６に示した例では、せん断波の伝搬方向（ｘ方向）と深さ方向（ｚ方向）とを２軸とする空間（平面）における変位を示している。図４（ａ）においては変位の大きさは、画素の輝度で表され、図５（ａ）、図６においては、変位の大きい領域がハッチングで示されている。また、図５（ａ）では、時系列の組織変位の分布を求めた結果が、時系列方向（ｔ方向）についても示されている。
Next, the
図６に示すように、発生したせん断波は、着目方向に伝搬する主成分（透過波）６１の他に、検査対象１００の生体組織の線維や脂肪などの構造物で反射され着目方向とは逆方向に伝搬する反射波６２や、構造物で屈折、回折、散乱等した屈折波・回折波・散乱波６３を含んでいる。そのため、せん断波全体の速度を計測した場合、反射波６２や屈折波・回折波・散乱波６３が主成分６１の速度を過小評価させる原因となる。そこで、本実施形態では、この後述べる選択部３２が、主成分６１を抽出する処理を行うことにより、反射波６２や屈折波・回折波・散乱波６３の影響を低減し、主成分６１の速度を精度よく計測する。
As shown in FIG. 6, in addition to the main component (transmitted wave) 61 propagating in the direction of interest, the generated shear wave is reflected by structures such as fibers and fat of the biological tissue of the
すなわち、制御部３０は、選択部３２に指示して、変位の周波数分布を求め、さらに、周波数ごとに、所定の強度（振幅）以上の波成分を選択する（ステップ１０４）。これにより、選択部３２は、主成分６１である透過波を中心とする速度成分を反射波６２や屈折波・回折波・散乱波６３から分離して抽出する。
That is, the
図７は、図２のステップ１０４における選択部３２の動作を示すフローチャートである。選択部３２は、ステップ１０３において計測部３１が算出した、せん断波の伝搬方向(ｘ方向）と深さ方向（ｚ方向）を２軸とする空間（平面）の各点についての変位の時間変化（すなわち、この平面を伝搬するせん断波の振幅の時間変化）を計測部３１から受け取り、特定の深度ｚ＝ｚ０の時空間データを抽出する（ステップ７０１）。抽出されたデータの例を図４（ｂ）、図５（ｂ）に示す。図４（ｂ）、図５（ｂ）は、選択部３２が算出した時間ｔと伝搬方向（ｘ方向）を２軸とする平面の変位（振幅）を示す説明図である。図４（ｂ）および図５（ｂ）のように、深度ｚ＝ｚ０の時空間データは、時間ｔと伝搬方向（ｘ方向）を２軸とする平面の変位（振幅）として表される。なお、図４（ｂ）は、着目方向の逆方向に伝搬する波がなく、主成分６１と屈折波・回折波・散乱波６３のみが伝搬している例であり、図５（ｂ）は、主成分６１と屈折波・回折波・散乱波６３が着目方向に伝搬する他に、反射波６２とその回折波や散乱波６３が着目方向の逆方向に伝搬している例を示している。
FIG. 7 is a flowchart showing the operation of the
つぎに、選択部３２は、２Ｄ−ＦＦＴ（２次元フーリエ変換）を実行し、時間ｔと伝搬方向（ｘ方向）を２軸とする平面の変位（振幅）を、時間周波数ｆと空間周波数ｋｘと強度（振幅）で表される波成分に分離する。図４（ｃ）、図５（ｃ）は、時間周波数ｆと空間周波数ｋｘを２軸とする空間（平面）の波成分の強度（振幅）分布（以下周波数パワー分布とも呼ぶ）である。選択部３２は、図４（ｃ）、図５（ｃ）のように、周波数パワー分布を得る（ステップ７０２）。
Next, the
つぎに、選択部３２は、周波数パワー分布（時間周波数ｆと空間周波数ｋｘを２軸とする空間における波成分の強度分布）から、時間周波数ｆごとまたは空間周波数ｋｘごとに所定の強度以上の波成分を選択する。図４（ｄ）、図５（ｅ）は、時間周波数ｆと空間周波数ｋｘを２軸とする空間に設定されたフィルタを示す説明図である。図４（ｄ）、図５（ｅ）のように、選択部３２は、時間周波数ｆまたは空間周波数ｋｘごとに波成分の強度（振幅）が最も大きい点を含むように範囲８０１を設定し、この範囲８０１のみを選択するフィルタを生成する。そのため、ここでは、選択部３２は、まず時間周波数ｆと空間周波数ｋｘを２軸とする空間の、波成分の強度（振幅）の最大値の点を、時間周波数ｆまたは空間周波数ｋｘごとにプロットする（ステップ７０３）。図８は、図４（ｃ）の波成分の強度の最大値をプロットしたグラフである。プロットされた点を含む所定の範囲を設定してもよいが、図８の例では、プロットした点に直線をフィッティングすることにより、その傾きによって表される、強度が最大の波成分の位相速度（Ｖ＝ｆ／ｋ）、すなわち実態としては群速度に相当する平均速度Ｖ０を求める（ステップ７０４）。つぎに、強度が最大の波成分の平均速度Ｖ０を含み、周波数空間（時間周波数ｆと空間周波数ｋｘを２軸とする空間）上で平均速度Ｖ０から所定の速度だけ小さい速度Ｖ１から平均速度Ｖ０から速度だけ大きい速度Ｖ２の範囲の波成分を抽出するための範囲８０１を周波数空間上に設定し、フィルタを作成する（ステップ７０５）。具体的には、図４（ｄ）、図５（ｅ）、図８のように、平均速度Ｖ０を示す直線から角度α１、α２の範囲８０１を設定し、フィルタを生成する。
Next, the
次に、選択部３２は、ステップ７０５で生成したフィルタ（図４（ｄ）、図５（ｅ））を、周波数パワー分布（図４（ｃ）、図５（ｃ））に適用し、範囲８０１の周波数パワー分布のみを抽出する（ステップ７０６）。これにより、選択部３２は、主成分６１である透過波を中心とする速度の波成分を反射波６２や屈折波・回折波・散乱波６３から分離して抽出することができる。図５（ｆ）は、フィルタ適用後の周波数パワー分布を示す。
Next, the
つぎに、選択部３２は、フィルタで抽出した主成分６１を中心とする速度の波成分に、２Ｄ−ＩＦＦＴ（２次元逆フーリエ変換）を施し、再び時間ｔと伝搬方向（ｘ方向）を２軸とする平面の変位（振幅）データに戻す（ステップ７０７）。
Next, the
つぎに、選択部３２は、ステップ７０１に戻り、ステップ深さｚをγだけ変えて、ステップ７０１〜７０７の処理を繰り返す（ステップ７０８、７０９）。これを全ての深度について実行する。図５（ｇ）は、深さ方向（ｚ方向）と伝搬方向（ｘ方向）と時間方向（ｔ方向）について、主成分６１の変位（振幅）を示す説明図である。ステップ７０１〜７０７の処理を繰り返すことにより、図５（ｇ）のように、深さ方向（ｚ方向）と伝搬方向（ｘ方向）と時間方向（ｔ方向）について、主成分６１の変位（振幅）を示す波面データが生成される。
Next, the
なお、図５（ｃ）に示したように、反射波６２の強度が大きい場合、周波数パワー分布には、反射波６２の波成分の強度も大きく表れる。そのため、選択部３２は、波成分の伝搬方向が、主成分６１方向の波成分を選択する処理をステップ７０２の後で行うことが望ましい。図５（ｄ）は、波方向の伝搬方向を選択する方向性フィルタである。例えば、図５（ｄ）のような方向性フィルタを用いて、時間周波数ｆと空間周波数ｋｘを２軸とする空間において、反射波６２の波成分が現れる領域（空間の第１象限と第３象限）のデータを除去し、主成分６１の波成分が現れる領域（第2象限と第４象限）のデータのみを選択する方向性フィルタを、図５（ｃ）の時間周波数ｆと空間周波数ｋｘを２軸とする空間の変位の分布に積分する。これにより、反射波を除去した周波数パワー分布が得られる。この反射波を除去した周波数パワー分布を用いて、ステップ７０３〜７０５においてフィルタを生成することができる。なお、方向性フィルタ５０１を、ステップ７０６の処理の後で積分してもよい。
Note that, as shown in FIG. 5C, when the intensity of the reflected
速度算出部３３は、選択部３２が生成した深さ方向（ｚ方向）と伝搬方向（ｘ方向）と時間方向（ｔ方向）について主成分６１の変位（振幅）を示す波面データに基づいて、速度を算出する（ステップ１０５）。
Based on the wavefront data indicating the displacement (amplitude) of the
図９は、図２のステップ１０５における速度算出部３３の処理を示すフローチャートである。まず、速度算出部３３は、図５（ｇ）のように、深さ方向（ｚ方向）と伝搬方向（ｘ方向）を２軸とする空間上に、一定間隔（ｄ）離れた２つの計測点５０１、５０２を設定する（ステップ９０１）。つぎに、速度算出部３３は、各計測点５０１，５０２における変位データの時間変化を選択部３２が生成したデータから抽出する。そして、計測点５０１，５０２ごとに、変位データの時間変化から波面特徴量（例えば、最大振幅値、最小振幅値、ゼロクロス点など）が現れる時間ｔを求め、計測点間の時間ｔの差Δｔを求めることにより、せん断波の主成分の到達時間差Δｔ（距離ｄを伝搬するのに要する時間（Δｔ）：Time of flight）を算出する（ステップ９０２）。この時間差Δｔで距離ｄを除することにより、せん断波の主成分の速度Ｖを求める（Ｖ＝ｄ／Δｔ）（ステップ９０３）。以上の処理を、深さ方向（ｚ方向）と伝搬方向（ｘ方向）を２軸とする空間全体で実行することにより、速度分布データを生成する（ステップ９０４）。
FIG. 9 is a flowchart showing the processing of the
画像生成部３４は、速度算出部３３が生成した深さ方向（ｚ方向）と伝搬方向（ｘ方向）を２軸とする空間の速度分布を用いて弾性率画像を生成し、表示部１５に表示させる（ステップ１０６）。図１０は、図２のステップ１０６における画像生成部３４の処理を詳しく示すフローチャートである。画像生成部３４は、深さ方向（ｚ方向）と伝搬方向（ｘ方向）を２軸とする空間の各点の速度データを速度算出部３３から受け取り（ステップ１００１）、予め定めた式（１）により、各点の弾性率Ｅを算出する（ステップ１００２）。
The
Ｅ＝３ρＶ２ ・・・（１）
ただし、Ｖは速度であり、ρは予め定めて定めておいた密度である。ここで速度Ｖは、周波数依存性の有無により位相速度と群速度に対応する。
E=3ρV 2 (1)
However, V is a velocity, and ρ is a density that is determined in advance. Here, the velocity V corresponds to the phase velocity and the group velocity depending on the presence or absence of frequency dependence.
画像生成部３４は、求めた弾性率Ｅの分布を、深さ方向（ｚ方向）と伝搬方向（ｘ方向）を２軸とする空間において示す弾性率画像を生成し（ステップ１００３）、表示部１５に表示させる（ステップ１００４）。図１１は、表示部１５の表示画面の領域１１１に、弾性率画像を示した画面例である。操作者は、主成分６１から得た弾性率画像から、組織の弾性を正確に把握することができる。
The
また、画像生成部３４は、受信ビームフォーマ２２の出力するＲＦ信号を受け取って、ＲＦ信号の振幅を輝度情報に変換し、検査対象１００の生体組織の形状を示す画像（Ｂモード像）を生成する。画像生成部３４は、生成したＢモード像を、図１１の領域１１２に表示する。これにより、操作者は、並べて表示された領域１１２のＢモード像と領域１１１の弾性率画像を見比べることができ、組織の形状と、弾性率との関係を把握することができる。また、画像生成部３４は、Ｂモード像と弾性率画像の一方を半透明化して重畳表示させてもよい。
Further, the
また、表示部１５の表示画面の領域１１３には、図１１のように、選択部３２がステップ７０２で生成した周波数パワー分布の画像に、平均速度Ｖ０を示す直線と、フィルタの範囲８０１を示す下限速度Ｖ１と上限速度Ｖ２を示す直線を示してもよい。これにより、操作者は、表示されている弾性率画像がどのようなフィルタによって選択された波成分から得た速度に基づく弾性率かを把握することができる。
Further, in the
さらに、画像生成部３４は、領域１１３の周波数パワー分布の近傍に、フィルタの範囲８０１の下限速度Ｖ１と上限速度Ｖ２が矢印で示された表示バー１１４を表示させ、入力部１３を介して、操作者が表示バー１１４上で、下限速度Ｖ１と上限速度Ｖ２を変更できるように構成してもよい。入力部１３を介して、操作者が下限速度Ｖ１と上限速度Ｖ２を変更した場合、選択部３２は、変更された下限速度Ｖ１と上限速度Ｖ２によりフィルタを生成して、図７のステップ７０１〜７０８を再度実行し、主成分６１の変位分布を再計算する。速度算出部３３および画像生成部３４は、再計算された変位分布を用いて、弾性率を求め、弾性率画像を生成し、領域１１１に再表示させる。これにより、操作者は、変更後の下限速度Ｖ１と上限速度Ｖ２で抽出された主成分６１で求めた弾性率を把握することができる。
Further, the
なお、ここでは弾性率を算出し、弾性率画像を表示する例について説明したが、表示する特性値は、弾性率に限られるものではなく、速度を用いて算出できる生体組織の特性値であればどのような値であってもよい。 It should be noted that although the example in which the elastic modulus is calculated and the elastic modulus image is displayed has been described here, the displayed characteristic value is not limited to the elastic modulus, and may be a characteristic value of a biological tissue that can be calculated using velocity. Any value will do.
また、表示画面の領域１１１には、弾性率画像に重畳させて、操作者が所望する画像を表示させる構成としてもよい。図１２(ａ)、（ｂ）は、操作が重畳表示させる画像を選択するボタン１２２を表示画面に表示させた画面例である。図１２（ａ）、（ｂ）では、Ｂモード像、速度画像、弾性率画像、粘性率画像、勾配画像、減衰画像、血管画像のボタンが表示されている。操作者が、入力部１３を介して、これらボタンのうち２以上を選択することにより、画像生成部３４は、それらの画像の重畳画像を生成して領域１１１に表示させる。図１２（ａ）では、Ｂモード像と弾性率画像が領域１１１に重畳表示され、図１２（ｂ）では、Ｂモード像と弾性率画像と血管画像が領域１１１に重畳表示されている。速度画像は、速度算出部３３が求めた速度の分布を画像にしたものである。粘性率画像は、粘性率の分布を示す画像であり、勾配画像は、位相速度の勾配を示す画像である。また、減衰画像は、超音波の減衰率を示す画像であり、血管画像は、ドプラー撮像や造影剤を注入して得た血流速度や血管分布等の血管情報を示す画像である。粘性率、勾配、減衰、血管の各画像は、いずれも公知の方法によって生成可能であるため、生成方法についての詳細な説明は省略するが、勾配画像と粘性率画像と減衰画像の生成方法についてのみ、以下簡単に説明する。
Further, in the
図１３は、勾配画像の生成処理を示すフローチャートである。図１３において、図７のフローと同様の処理を行うステップには、図７と同じ符号を付している。制御部３０は、計測部３１が算出した、せん断波の伝搬方向(ｘ方向）と深さ方向（ｚ方向）の空間（平面）の各点についての変位の時間変化から、特定の深度ｚ＝ｚ０の時空間データを抽出し（ステップ７０１）、さらに、特定の方位方向の範囲（ｘ方向）のデータを抽出する（ステップ１３１）。制御部３０は、抽出したデータに対して、２次元フーリエ変換を施した後、変位が最大値を取る点をプロットする（ステップ７０２，７０３）。プロットした点ごとに、位相速度Ｖを、Ｖ＝ｆ／ｋ（ｆ：時間周波数、ｋ：空間周波数）にＶより算出し、図１４のように、時間周波数ｆと位相速度Ｖの空間にプロットする。制御部３０は、プロットした点に直線または曲線をフィッティングし、その速度勾配を求める（ステップ１３２）。これを、全方位方向（ｘ方向）および全深さ方向（ｚ方向）について行うことにより、方位方向（ｘ方向）と深さ方向（ｚ方向）を２軸とする空間における速度勾配分布が得られる（ステップ１３３、７０８）。画像生成部３４は、得られた速度勾配分布から勾配画像を生成する（ステップ７０９）。
FIG. 13 is a flowchart showing the process of generating a gradient image. 13, steps that perform the same processing as in the flow of FIG. 7 are assigned the same reference numerals as in FIG. 7. The
一方、図１５は、粘性画像および減衰率画像の生成処理を示すフローチャートである。図１５において、図１３のフローと同様の処理を行うステップには、図１３と同じ符号を付している。制御部３０は、ステップ７０１、ステップ１３１、ステップ７０２および７０３までは、勾配画像を生成する処理と同様であり、図１４のように、変位が最大値を取る点の位相速度Ｖを時間周波数ｆにごとにプロットする。このプロットした点に対して、公知のVoigtモデルやMaxwellモデルによる機械モデルフィッティングを行って、粘性値および減衰率を推定する（ステップ１４１）。これを、全方位方向（ｘ方向）および全深さ方向（ｚ方向）について行うことにより、方位方向（ｘ方向）と深さ方向（ｚ方向）を２軸とする空間における粘性分布および減衰率分布が得られる。画像生成部３４は、得られた粘性分布および減衰率分布からそれぞれ粘性画像および減衰率画像を生成する（ステップ１４２）。
On the other hand, FIG. 15 is a flowchart showing the generation processing of the viscosity image and the attenuation rate image. In FIG. 15, steps that perform the same processing as in the flow of FIG. 13 are assigned the same reference numerals as in FIG. The
上述してきた実施形態では、波成分の強度分布から所定の強度以上の波成分をフィルタを用いて選択部３２が選択する構成であったが、選択部３２が、所定の速度範囲の波成分を選択する構成とすることも可能である。
In the above-described embodiment, the
上述してきたように、本実施形態によれば、せん断波の屈折波・回折波・散乱波６３等を除去して主成分６１の速度を精度よく測定できるため、弾性率等の組織の特性を精度よく算出できる。よって、超音波を用いて腫瘍を診断する場合、組織形状を映像化した画像（例えばＢモード像）の輝度情報を利用して腫瘍の位置、サイズおよび形状の評価を行うだけでなく、せん断波の速度に基づいて得られる腫瘍の弾性率等の特性値を画像に基づいて、腫瘍の性状の評価を行うことができる。また、得られる血管分布や、粘性率等の他の情報も同時に参照することもできる。このように、形状と性状の両情報を活用することで、より確定診断に近い超音波診断が実現し、CTやMRIなどの追加検査が回避できる。その結果、被ばくなどの患者への身体的負担だけでなく、医療費に関する経済的負担の軽減効果が期待できる。
As described above, according to the present embodiment, since the refraction wave/diffraction wave/scattered
１…超音波送受信装置、３…入力部、１０…探触子、１３…入力部、１５…表示部、１６…メモリ、２０…送受信制御部、２１…送信ビームフォーマ、２２…受信ビームフォーマ、３０…制御部、３１…計測部、３２…選択部、３３…速度算出部、３４…画像生成部、１００…検査対象
DESCRIPTION OF
Claims (10)
前記計測部の計測結果に基づいて、前記変位の周波数分布を求め、周波数ごとに、所定の強度以上の波成分を選択する選択部と、
前記選択部が選択した前記波成分の速度を求める速度算出部とを有し、
前記選択部は、前記計測部が計測した前記複数の計測点の変位の時間変化を、時間周波数と空間周波数を２軸とする空間における前記波成分の強度分布に変換し、前記時間周波数と空間周波数の空間の前記波成分の強度分布から、前記時間周波数ごとに前記所定の強度以上の波成分を選択する
ことを特徴とする超音波信号処理装置。 An ultrasonic signal obtained by transmitting and receiving ultrasonic waves to and from the inspection object is received, and a measuring unit that measures displacements of a plurality of measurement points of the inspection object due to shear waves propagating through the inspection object in time series,
Based on the measurement result of the measuring unit, obtain the frequency distribution of the displacement, for each frequency, a selection unit for selecting a wave component of a predetermined intensity or more,
A velocity calculation unit that obtains the velocity of the wave component selected by the selection unit,
The selecting unit converts a temporal change in displacement of the plurality of measurement points measured by the measuring unit into an intensity distribution of the wave component in a space having two axes of a temporal frequency and a spatial frequency, and the temporal frequency and the space. An ultrasonic signal processing device , wherein a wave component having a predetermined intensity or more is selected for each time frequency from the intensity distribution of the wave component in a frequency space .
前記選択部は、前記入力部を介して前記操作者から受け付けた範囲に前記フィルタの範囲を設定することを特徴とする超音波信号処理装置。 The ultrasonic signal processing apparatus according to claim 3 , further comprising an input unit that receives a setting of a range selected by the filter from an operator,
The ultrasonic signal processing device, wherein the selection unit sets the range of the filter to a range received from the operator via the input unit.
前記計測部の計測結果に基づいて、前記変位の時間変化の波を構成する波成分のうち、所定の速度範囲の波成分を選択する選択部と、
前記選択部が選択した前記波成分の速度を求める速度算出部とを有し、
前記選択部は、前記計測部が計測した前記複数の計測点の変位の時間変化を、時間周波数と空間周波数を２軸とする空間における前記波成分の強度分布に変換し、前記空間の前記波成分の強度の最大値の点に直線をフィッティングし、前記直線の傾きから前記強度が最大の波成分の位相速度を求め、求めた前記位相速度を含む所定の位相速度範囲の波成分を抽出するための範囲を前記空間上で選択する
ことを特徴とする超音波信号処理装置。 An ultrasonic signal obtained by transmitting and receiving ultrasonic waves to and from the inspection object is received, and a measurement unit that measures displacements of a plurality of measurement points of the inspection object by a shear wave propagating through the inspection object in time series,
Based on the measurement result of the measurement unit, a selection unit that selects a wave component in a predetermined velocity range among the wave components that form the wave of the displacement with time.
A velocity calculation unit that obtains the velocity of the wave component selected by the selection unit,
The selecting unit converts a temporal change in displacement of the plurality of measurement points measured by the measuring unit into an intensity distribution of the wave component in a space having a temporal frequency and a spatial frequency as two axes, and converts the wave in the space. A straight line is fitted to the point of the maximum value of the intensity of the component, the phase velocity of the wave component having the maximum intensity is obtained from the slope of the straight line, and the wave component in a predetermined phase velocity range including the obtained phase velocity is extracted. An ultrasonic signal processing device, characterized in that a range for selecting is selected in the space .
前記超音波信号処理装置は、請求項１ないし９のいずれか１項に記載の超音波信号処理装置であることを特徴とする超音波送受信装置。 A transmission/reception control unit that transmits and receives ultrasonic waves to and from the inspection target, and an ultrasonic signal processing device that processes the ultrasonic signals received by the transmission and reception control unit,
The ultrasonic signal processing apparatus is the ultrasonic signal processing apparatus according to any one of claims 1 to 9 .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016245761A JP6698511B2 (en) | 2016-12-19 | 2016-12-19 | Ultrasonic signal processing device and ultrasonic transmitting/receiving device using the same |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016245761A JP6698511B2 (en) | 2016-12-19 | 2016-12-19 | Ultrasonic signal processing device and ultrasonic transmitting/receiving device using the same |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018099180A JP2018099180A (en) | 2018-06-28 |
| JP6698511B2 true JP6698511B2 (en) | 2020-05-27 |
Family
ID=62714809
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016245761A Active JP6698511B2 (en) | 2016-12-19 | 2016-12-19 | Ultrasonic signal processing device and ultrasonic transmitting/receiving device using the same |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6698511B2 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7490387B2 (en) * | 2019-02-26 | 2024-05-27 | キヤノンメディカルシステムズ株式会社 | Ultrasound diagnostic equipment |
| US11850099B2 (en) | 2019-02-26 | 2023-12-26 | Canon Medical Systems Corporation | Ultrasonic diagnostic apparatus |
| JP7236312B2 (en) * | 2019-04-04 | 2023-03-09 | 富士フイルムヘルスケア株式会社 | Ultrasound diagnostic device, signal processing device, and program |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6393502B2 (en) * | 2014-04-09 | 2018-09-19 | 国立大学法人京都大学 | Method for obtaining elastic properties of biological tissue |
-
2016
- 2016-12-19 JP JP2016245761A patent/JP6698511B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2018099180A (en) | 2018-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5735718B2 (en) | Ultrasonic diagnostic apparatus and elasticity evaluation method | |
| US9072493B1 (en) | Ultrasonic diagnostic apparatus and elastic evaluation method | |
| JP5730978B2 (en) | Ultrasonic diagnostic apparatus and method | |
| JP6216736B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic method | |
| JP6353929B2 (en) | Ultrasonic diagnostic apparatus and elasticity evaluation method | |
| JP6987496B2 (en) | Analyst | |
| WO2015040710A1 (en) | Diagnostic ultrasound apparatus, medical image-processing device and medical image-processing method | |
| KR20150070859A (en) | Method and apparatus for obtaining elasticity information of interest of region using shear wave | |
| JP2008136860A (en) | Ultrasonic diagnostic apparatus and image processing program for it | |
| JP6698511B2 (en) | Ultrasonic signal processing device and ultrasonic transmitting/receiving device using the same | |
| KR20140036650A (en) | Method, apparatus and system for analysing elastography of tissue using 1-dimensional ultrasound probe | |
| US11490876B2 (en) | Ultrasonic diagnostic device and method for evaluating physical properties of biological tissue | |
| JP6861624B2 (en) | Ultrasonic transmitter / receiver and ultrasonic transmitter / receiver method | |
| KR102545007B1 (en) | Ultrasound imaging apparatus and controlling method for the same | |
| CN111789632B (en) | Ultrasonic diagnostic apparatus, signal processing apparatus, and storage medium | |
| JP7354632B2 (en) | Ultrasonic diagnostic device and method of controlling the ultrasonic diagnostic device | |
| JP2015128554A (en) | Ultrasonic diagnostic equipment | |
| Urban et al. | Harmonic pulsed excitation and motion detection of a vibrating reflective target | |
| JP6537706B2 (en) | ULTRASONIC IMAGING DEVICE AND METHOD OF PROCESSING ULTRASONIC RECEIVING SIGNAL | |
| Soler López et al. | Application of ultrasound in medicine part ii: the ultrasonic transducer and its associated electronics | |
| JP7347445B2 (en) | Ultrasonic signal processing device, ultrasonic diagnostic device, and ultrasonic signal processing method | |
| Marsico et al. | Elastography: Technical Aspects | |
| Silvestri et al. | US Basic Principles |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190304 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200109 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200121 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200304 |
|
| 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: 20200421 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200428 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6698511 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| 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 |