JP7432426B2 - Ultrasonic diagnostic equipment, signal processing equipment, and signal processing programs - Google Patents
Ultrasonic diagnostic equipment, signal processing equipment, and signal processing programs Download PDFInfo
- Publication number
- JP7432426B2 JP7432426B2 JP2020066134A JP2020066134A JP7432426B2 JP 7432426 B2 JP7432426 B2 JP 7432426B2 JP 2020066134 A JP2020066134 A JP 2020066134A JP 2020066134 A JP2020066134 A JP 2020066134A JP 7432426 B2 JP7432426 B2 JP 7432426B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- roi
- representative value
- dimensional
- representative
- 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
- 238000012545 processing Methods 0.000 title claims description 33
- 238000002604 ultrasonography Methods 0.000 claims description 39
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 21
- 239000000523 sample Substances 0.000 claims description 21
- 230000007423 decrease Effects 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 description 20
- 230000005540 biological transmission Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 11
- 238000012937 correction Methods 0.000 description 10
- 210000001519 tissue Anatomy 0.000 description 7
- 238000007689 inspection Methods 0.000 description 6
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 208000004930 Fatty Liver Diseases 0.000 description 2
- 206010019708 Hepatic steatosis Diseases 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 208000010706 fatty liver disease Diseases 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 208000019423 liver disease Diseases 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 231100000240 steatosis hepatitis Toxicity 0.000 description 2
- 206010012426 Dermal cyst Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 206010033675 panniculitis Diseases 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 210000004304 subcutaneous tissue Anatomy 0.000 description 1
- 230000002123 temporal 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/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
本発明は、超音波診断装置に係り、被検体内の音響減衰特性の計測手法に関わる。 The present invention relates to an ultrasonic diagnostic apparatus, and relates to a method for measuring acoustic attenuation characteristics within a subject.
超音波やMRI(Magnetic Resonance Imaging)装置、X線CT(Computed Tomography)装置に代表される医療用の画像表示装置は、目視できない生体内の情報を数値または画像の形態で提示する装置として広く利用されている。中でも超音波を利用して画像を表示する超音波撮像装置は、他の装置と比較して高い時間分解能を備えており、例えば拍動下の心臓を滲みなく画像化できる性能を持つ。 Medical image display devices, such as ultrasound, MRI (Magnetic Resonance Imaging) devices, and X-ray CT (Computed Tomography) devices, are widely used to present in-vivo information that cannot be seen with the naked eye in the form of numerical values or images. has been done. Among these, ultrasonic imaging devices that display images using ultrasound have higher temporal resolution than other devices, and have the ability, for example, to image a beating heart without blurring.
超音波診断装置は、検査対象に向けて超音波を送信し、散乱体からの反射信号(以下受信RF(Radio Frequency)信号と呼ぶ)を受信して画像を構成する。基本的には送受信に要する時間と音速から散乱体までの距離を計測し、超音波エネルギーの強度に基づく輝度の空間分布を構成することで超音波画像(以下、B像)が生成される。 An ultrasonic diagnostic apparatus transmits ultrasonic waves toward an object to be examined, receives reflected signals from scatterers (hereinafter referred to as received RF (Radio Frequency) signals), and forms an image. Basically, an ultrasound image (hereinafter referred to as B image) is generated by measuring the distance to the scatterer based on the time required for transmission and reception and the speed of sound, and constructing a spatial distribution of brightness based on the intensity of ultrasound energy.
また、超音波診断装置は、受信信号を用いて、血流速度や、組織の硬さや、減衰特性等の音響的特性を定量計測することにより、疾患の診断をする技術の開発も提案されている。 In addition, it has been proposed to develop technology for diagnosing diseases using ultrasound diagnostic equipment by quantitatively measuring acoustic characteristics such as blood flow velocity, tissue hardness, and attenuation characteristics using received signals. There is.
例えば、特許文献1には、生体組織の減衰特性を求める手法が開示されている。特許文献1では、まず、探触子から対象物に対して、焦点距離の異なる第1及び第2の集束させた超音波ビームを送信し、対象物の所定の計測点についての第1および第2の超音波ビームによる受信信号をそれぞれ得る。次に、減衰特性算出部は、所定の計測点について得た、第1の超音波ビームによる受信信号と第2の超音波ビームによる受信信号の超音波エネルギーの深度方向の変化を用いて対象物の減衰特性を求めている。
For example,
また、特許文献2には、超音波画像に含まれる各種のノイズやスペックルを低減し、被検体組織に関する情報を強調するために、画像のフィルタリング処理を行う技術であって、筋組織構造のように、局所領域でみた場合に組織構造が特定の方向へ揃っているときには、その構造をはっきりとさせるようにフィルタのパラメータを調整する方法が提案されている。すなわち、超音波画像における局所領域のエッジの方向と、前記局所領域におけるエッジの大きさと、局所領域におけるエッジの方向一様性を算出し、フィルタのパラメータを設定する。 Furthermore, Patent Document 2 discloses a technology that performs image filtering processing in order to reduce various noises and speckles included in ultrasound images and emphasize information regarding the subject's tissue. A method has been proposed in which when the tissue structure is aligned in a specific direction when viewed in a local region, the filter parameters are adjusted to make the structure clear. That is, the direction of the edge of the local region in the ultrasound image, the size of the edge in the local region, and the uniformity of the direction of the edge in the local region are calculated, and the parameters of the filter are set.
生体組織の音響特性の1つとして、組織内部を伝播する音響エネルギーの減衰特性がある。減衰特性とは、超音波が媒質中を伝播する際、超音波エネルギーがその伝播距離に応じて、吸収、散乱の現象により減衰する特性である。肝疾患の1つである脂肪肝は、正常肝と比較して音響エネルギーの減衰が大きくなることが知られている。そこで、減衰率を計測することにより、疾患の診断に役立てることが期待されている。 One of the acoustic properties of living tissue is the attenuation property of acoustic energy propagating inside the tissue. The attenuation characteristic is a characteristic in which when an ultrasonic wave propagates in a medium, the ultrasonic energy is attenuated due to absorption and scattering phenomena depending on the propagation distance. Fatty liver, a type of liver disease, is known to cause greater attenuation of acoustic energy compared to normal liver. Therefore, it is expected that measuring the attenuation rate will be useful in diagnosing diseases.
超音波の信号強度は、伝搬距離が大きくなるにつれ滑らかに減少するが、実際にはランダムな信号強度の変化(スペックル)が含まれている。減衰率の算出は、距離当たりの信号強度の低下率から算出するため、このスペックルが、計測値のバラつきを引き起こし、減衰率の算出精度を低下させるという問題がある。 Although the signal strength of ultrasound decreases smoothly as the propagation distance increases, it actually includes random changes in signal strength (speckles). Since the attenuation rate is calculated from the rate of decrease in signal strength per distance, there is a problem in that this speckle causes variations in the measured values and reduces the accuracy of the attenuation rate calculation.
また、超音波診断装置には、深度に応じて受信信号にゲインをかけるTime Gain Control (TGC)機能が実装されており、超音波エネルギーの強度分布を輝度分布として表すB像の輝度が、浅部から深部まで一様になるように調整する。しかし、実際の受信RF信号の深度方向の減衰は、深度に応じていない場合がある。例えば、送信超音波が集束されて照射される場合、受信RF信号は送信超音波の焦点付近で強度が大きくなる。 In addition, the ultrasound diagnostic equipment is equipped with a Time Gain Control (TGC) function that applies a gain to the received signal according to the depth, and the brightness of the B image, which represents the intensity distribution of ultrasound energy as a brightness distribution, is Adjust so that it is uniform from the top to the bottom. However, the actual attenuation of the received RF signal in the depth direction may not depend on the depth. For example, when transmitted ultrasound is focused and irradiated, the intensity of the received RF signal increases near the focal point of the transmitted ultrasound.
また、実際の受信RF信号の減衰は、深度方向だけでなく、深度方向と方位方向に2次元的に発生する。例えば、肝疾患の診断を想定した場合、肝実質と皮下組織やのう胞、横隔膜では音速、密度、減衰率が異なるため、構造物の境界で超音波が反射、屈折、散乱を起こし、深度方向と方位方向の音圧プロファイルに対して2次元的に影響を与える。 Furthermore, the actual attenuation of the received RF signal occurs not only in the depth direction but also two-dimensionally in the depth direction and the azimuth direction. For example, when diagnosing liver disease, the sound speed, density, and attenuation rate differ between the liver parenchyma, subcutaneous tissue, cyst, and diaphragm, so ultrasound waves are reflected, refracted, and scattered at the boundaries of structures, and Affects the sound pressure profile in the azimuth direction two-dimensionally.
そのため、深度方向のみのTGC機能では、深度方向のB像の輝度分布を一様にすることはできず、依然B像上で生体組織の形態が視認しにくくなり、診断が困難となる場合がある。 Therefore, with the TGC function only in the depth direction, it is not possible to make the brightness distribution of the B image uniform in the depth direction, and the morphology of the living tissue may still be difficult to see on the B image, making diagnosis difficult. be.
本発明の目的は、超音波の減衰率を方向にかかわらず精度よく算出可能にすることにある。 An object of the present invention is to enable accurate calculation of the attenuation rate of ultrasonic waves regardless of direction.
上記目的を達成するために、本発明の超音波診断装置は、超音波探触子から信号を受け取って受信ビームフォーミングを行うことにより、撮像領域に設定された複数の受信走査線についてそれぞれ受信信号を生成する受信ビームフォーマと、受信信号を処理する信号処理部とを有する。信号処理部は、複数の受信信号を受信走査線の順に並べた2次元信号空間に2次元のROIを複数設定する2次元ROI設定部と、ROI内の信号を用いて代表値を求め、求めた代表値をROI内に設定した代表点の処理後信号値とする代表値算出部とを有する。 In order to achieve the above object, the ultrasound diagnostic apparatus of the present invention receives signals from an ultrasound probe and performs reception beamforming, thereby providing a reception signal for each of a plurality of reception scanning lines set in an imaging area. and a signal processing unit that processes the received signal. The signal processing unit includes a two-dimensional ROI setting unit that sets a plurality of two-dimensional ROIs in a two-dimensional signal space in which a plurality of received signals are arranged in the order of received scanning lines, and a two-dimensional ROI setting unit that determines a representative value using the signal in the ROI. and a representative value calculation unit that uses the representative value as the processed signal value of the representative point set within the ROI.
本発明の超音波診断装置によれば、2次元のROIを複数設定し、それぞれの代表値を処理後信号値とするため、超音波の減衰率を方向にかかわらず精度よく算出することができる。 According to the ultrasound diagnostic apparatus of the present invention, since a plurality of two-dimensional ROIs are set and each representative value is used as the processed signal value, the attenuation rate of ultrasound can be calculated with high accuracy regardless of the direction. .
以下、本発明の実施形態を図面に従い説明する。 Embodiments of the present invention will be described below with reference to the drawings.
<<実施形態1>>
<超音波診断装置の全体構成>
図1に実施形態の超音波診断装置1の一構成例のブロック図を示す。本実施形態の超音波送受信装置1は、送受信制御部20と、信号処理部30を備えている。また、送受信制御部20には、探触子10が接続されている。信号処理部30には、外部入力デバイス13と、表示部16が接続されている。
<<
<Overall configuration of ultrasound diagnostic equipment>
FIG. 1 shows a block diagram of a configuration example of an ultrasonic
送受信制御部20は、送信ビームフォーマ21と受信ビームフォーマ22とを備えている。送信ビームフォーマ21は、送信信号を生成し、探触子10を構成する配列された振動子に出力する。探触子10の複数の振動子は、送信信号を超音波に変換し、検査対象100に向かって送信する。検査対象100において反射や散乱された超音波のうち探触子10に到達した超音波は、探触子10の各振動子によって電気信号に変換される。受信ビームフォーマ22は、各振動子の電気信号を受け取って、撮像領域に設定した受信走査線上の複数の点について焦点を合わせるように遅延させて加算(整相加算)することにより、受信信号(以下、受信RF信号と称する)を生成(受信ビームフォーミング)する。
The transmission/
信号処理部30は、2次元ROI設定部31と、ROI適合性分析部32と、代表値算出部33と、代表値適合性判定部34と、補間部35と、画像生成部36と、減衰率算出部37とを備えて構成される。
The
2次元ROI設定部31は、複数の受信RF信号を受信走査線の順に並べた2次元信号空間を形成し、2次元信号空間に2次元のROIを複数設定する。
The two-dimensional
代表値算出部33は、ROI内の信号を用いて代表値を求め、求めた代表値をROI内に設定した代表点の処理後信号値とする。
The representative
このように、2次元のROIを受信RF信号の2次元信号空間に複数設定し、それぞれの代表値を求め、ROIの代表値の処理後信号値とするためROI内の受信信号に含まれるスペックルを2次元空間において低減することができる。これにより、超音波の減衰率を方向にかかわらず精度よく算出することができる。 In this way, multiple two-dimensional ROIs are set in the two-dimensional signal space of the received RF signal, the representative value of each is determined, and the specifications included in the received signal within the ROI are used to obtain the processed signal value of the representative value of the ROI. can be reduced in two-dimensional space. Thereby, the attenuation rate of ultrasonic waves can be calculated with high precision regardless of the direction.
上記2次元ROI設定部31は、複数のROIを、2次元信号空間の受信走査線に沿った方向および受信走査線の配列方向の少なくとも一方について、互いにずれた位置に設定することが望ましい。
It is preferable that the two-dimensional
このとき、2次元ROI設定部31は、設定するROIのサイズや形状を、受信RF信号に基づいて適応的に変更してもよい。
At this time, the two-dimensional
ROI適合性分析部32は、ROI内の信号強度の分布が、予め定めた信号強度分布に適合しているか否かを判定し、適合しているROIを有効ROIとし、適合していないROIを無効ROIとする。
The ROI
代表値適合性判定部34は、代表値算出部33が算出した複数のROIについての代表値の分布が、予め定めた条件に適合しているかを判定する。代表値適合性判定部34は、条件に適合していない代表値を処理後信号値として採用しない。例えば、代表値適合性判定部34は、複数のROIの代表値の信号強度の分布が、音波伝播方向にそって徐々に減少している場合、条件に適合すると判定する。
The representative value
補間部35は、2次元信号空間において、複数のROIの代表値を用いて、ROIとROIの間の2次元信号空間の信号値を補間により求める。
The
画像生成部36は、処理後信号を用いて、Bモード像を生成する。画像生成部36は、処理後信号を用いて、Bモード像の輝度値を最適化する構成としてもよい。
The
減衰率算出部37は、処理後信号を用いて、深さ方向の減衰率を求める。
The attenuation
信号処理部30は、CPU(Central Processing Unit)やGPU(Graphics Processing Unit)等のプロセッサーと、メモリとを備えたコンピュータ等によって構成され、CPUが、メモリに格納されたプログラムを読み込んで実行することにより、上述した2次元ROI設定部31と、ROI適合性分析部32と、代表値算出部33と、代表値適合性判定部34と、補間部35と、画像生成部36と、減衰率算出部37の機能をソフトウエアにより実現する。なお、信号処理部30は、その一部または全部をハードウエアによって実現することも可能である。例えば、ASIC(Application Specific Integrated Circuit)のようなカスタムICや、FPGA(Field-Programmable Gate Array)のようなプログラマブルICを用いて信号処理部7を構成し、信号処理部7の各部の機能を実現するように回路設計を行えばよい。
The
図2に実施形態の超音波診断装置の一構成例のフローチャートを示す。全体の動作を、このフローチャートに沿って説明する。 FIG. 2 shows a flowchart of a configuration example of the ultrasonic diagnostic apparatus according to the embodiment. The entire operation will be explained along this flowchart.
まず、探触子10から図3のよういに超音波を送受信する。具体的には、送信ビームフォーマ21によって生成された送信信号が探触子10に出力され、探触子10から、超音波のビームである送信超音波301が送信され、検査対象100内部の散乱体303から散乱してきた反射超音波302を探触子10が受信する。受信ビームフォーマ22は、探触子10からの信号を整相加算することにより、設定した受信走査線に沿った受信RF信号を生成する。受信RF信号は、整相加算処理後、さらに常用対数で圧縮したものである。
First, ultrasonic waves are transmitted and received from the
なお、受信RF信号には、スペックルと呼ばれるランダムな信号の強度変化が含まれている。これは、媒質内の散乱体303から散乱された超音波がランダムに干渉しあうことで発生するもので、媒質の音響特性や、媒質内の散乱体の形状および密度によって変化する。また、検査対象の媒質は減衰があることを想定している。従って、この媒質内部を超音波が伝播すると、伝播距離に応じて、超音波エネルギーが減少する。 Note that the received RF signal includes random signal intensity changes called speckles. This occurs when the ultrasonic waves scattered from the scatterers 303 in the medium randomly interfere with each other, and changes depending on the acoustic characteristics of the medium and the shape and density of the scatterers in the medium. Furthermore, it is assumed that the medium to be inspected has attenuation. Therefore, when an ultrasonic wave propagates inside this medium, the ultrasonic energy decreases depending on the propagation distance.
ステップ201では、2次元ROI設定部31には、受信走査線に沿った受信RF信号が受信ビームフォーマ22から入力される。
In
2次元ROI設定部31は、複数の受信RF信号を受信走査線の順に横軸(方位)方向(探触子の振動子の並び方向)に並べることにより、2次元信号空間を形成する。したがって、この2次元信号空間は、横軸が方位方向、すなわち探触子10の振動子の並び方向、縦軸が深さ方向にあたる時間となっている。
The two-dimensional
次に、ステップ202では、送信した超音波の音圧プロファイルを用いて受信信号を補正し、補正後信号を得る。これにより、送信した超音波の音圧の強弱に起因する受信RF信号の信号強度の分布を補正する。音圧プロファイルは、均質なファントムもしくはシミュレーションによって予め計測や算出を行い、超音波診断装置内のメモリに保存したものを使用する。補正方法としては、受信RF信号をM(i,t),音圧プロファイルをR(i,t)、補正後の信号をC(i,t)として、
次に、ステップ203では、2次元ROI設定部31は、設定する複数のROIの仕様及び位置を決定する。ここでは、2次元ROI設定部31は、複数のROIを、2次元信号空間の受信走査線に沿った方向および受信走査線の配列方向の少なくとも一方について、互いにずれた位置に設定する。また、2次元ROI設定部31は、設定するROIのサイズや形状を、受信RF信号に基づいて適応的に変更する。
Next, in
図4を用いて、ROI仕様と、位置の設定方法の例を説明する。ROI401は、受信RF信号中に深さ方向のサイズがサンプル数h、方位方向のサイズがサンプル数wのサイズであり、その代表点(例えば中心)を(xn,zn)の位置にくるように設定されたROIである。ROI402は、ROI401とサイズは同じであるが、中心位置を深さ方向にサンプル数dh、方位方向にサンプル数dwだけ異なる位置に設定したものである。よって、ROI401の代表点(xn,zn)に対し、ROI402を代表点の位置は(xn+dw,zn+dhs)となる。
An example of ROI specifications and a position setting method will be explained using FIG. 4. The
2次元ROI設定部31によるROI401,402の深さ方向のサイズ(サンプル数h)と方位方向のサイズ(サンプル数w)の決定方法、ならびに、ROI401とROI402の位置ずれ量(深さ方向のサンプル数dh、方位方向のサンプル数dws)の決定方法の一例を図5を用いて説明する。
The two-dimensional
2次元ROI設定部31は、図5のステップ501において、ステップ201において、受信RF信号を並べて形成した2次元信号空間を受け取る。次に,ステップ502において2次元フーリエ変換し,受信RF信号の空間周波数の分布を求め,2次元ROI設定部31は,この空間周波数分布に基づき,スペックルの方位方向,深さ方向のサイズの分布を求める。
In
次にステップ503にて、2次元ROI設定部31は、スペックルの空間的なサイズ分布に基づき、ROIのサイズを決定する。具体的な決定方法としては、例えば、ROIの深さ方向のサンプル数h、方位方向のサンプル数w及び深さ方向の移動サンプル数dh、方位方向の移動サンプル数dwを、スペックルの空間サイズの平均の任意倍にする。
Next, in
次に、図2のステップ204では、ROI適合性分析部32は、2次元ROI設定部31が設定したROIが、ステップ205以降の処理に使用するかどうか、適合性を判断する。適合性の判断の方法の模式図を図6に示す。まず、ROI適合性分析部32は、2次元ROI設定部31が設定したROIがROI601、ROI602である場合、それぞれROI内の各点の受信RF信号の信号強度(輝度)分布603、信号強度(輝度)分布604を算出する。そして、ROI適合性分析部32は、各々の信号強度分布と、理想データ605との違いを定量化する。なお、理想データは、予め均質ファントムやシミュレーションなどで算出し、超音波診断装置のメモリに保存したものを使用する。
Next, in
図6において、ROI601の信号強度(輝度)分布603は、理想分布605と非常に似ている。一方、構造物606の含まれているROI602の信号強度(輝度)分布604は、理想分布605から乖離している。ROI適合性分析部32は、このような理想分布607からの乖離度から、ROIの適合性の有無を判定する。
In FIG. 6, a signal intensity (brightness) distribution 603 of an
ROI適合性分析部32によって、適合性があると判断されたROIに関しては、次のステップ205にて、代表値算出部33が代表値の算出を行う。
For the ROI that is determined to be compatible by the ROI
図7のフローを用いて、上記ステップ204、205のROI適合性分析部32による乖離度の算出と、代表値算出部33による代表値の算出の処理フローを具体的に説明する。
The process flow of calculating the degree of deviation by the ROI
図7のステップ701において、2次元ROI設定部31は、図4で示したようにひとつ前に設定したROIの位置に基づいて、ROIの設定位置を決定し、設定する。
In
ステップ702は、ROI適合性分析部32は、ROI中の信号強度(輝度)ヒストグラムRmeasuredを算出する。
In
ステップ703は、ROI適合性分析部32は、ROI中の信号強度(輝度
ヒストグラムRmeasuredと理想の輝度ヒストグラムRmodelとの乖離度Sを、
ステップ704にて、ROI適合性分析部32は、予め設定した閾値と、式(2)で算出した乖離度Sの比較によって適合性を判定する。すなわち、ROI適合性分析部32は、予め設定した閾値よりも乖離度Sが小さい場合、適合性があると判定する。
In
ステップ705において、ROI適合性分析部32は、適合性があると判定したROIに対しては、有効ROIとしてROIの代表点(xn,zn)の値μをROI中の値から算出する。代表点の値としては、ROI中の平均値や、中央値、最大値、最小値といった値を用いる。なお、本実施形態では、最初の代表点の位置をROIの中央の点としているが、中央の点に限られず重心やいずれかの辺の中点や四隅の点の予め定めた位置の点としてもよい。
In
また、ROI適合性分析部32は、適合性がないと判断したROIに関しては、無効ROIとしてステップ706において代表点の値を定数τと設定する。なお、τは、0など、受信RF信号上に現れない定数を決定する。
In addition, the ROI
ステップ207において、ROI適合性分析部32は、算出した代表値をメモリ(付図示)に保存する。
In
最後に、ステップ208(ステップ707)にて、ROI適合性分析部32は、全ての領域の分析が完了するまで上記処理を繰り返す。
Finally, in step 208 (step 707), the ROI
図8は、図7のステップ701から707の処理途中の様子を示した図である。受信RF信号の2次元信号空間805中に、適合すると判定された有効ROI803と、適合しないと判定された無効ROI804が位置することがわかる。各々のROIに対して、代表点802が設定されている。
FIG. 8 is a diagram showing the process during
図9は、図7のステップ701から707の処理を完了した際の様子を示した図である。処理が完了すると、無効ROIの代表点の集合体である無効領域の903と、有効ROIの代表点の集合体である有効領域904が受信RF信号805中に分布する形になる。
FIG. 9 is a diagram showing the state when the processing from
つぎに、代表値適合性判定部34は、代表値算出部33が算出した複数のROIについての代表値の分布が、予め定めた条件に適合しているかを判定する。
Next, the representative value
まず、図2のステップ208において、代表値適合性判定部34は代表値の選出を行う。図10は、代表値の選出のフローチャートである。
First, in
まず、ステップ1001で、代表値適合性判定部34は、代表点間を線分で結び、ステップ1002で、図11に示した、3点に囲まれた領域1102を作成する。
First, in
次に、ステップ1003にて、代表値適合性判定部34は、その3点に囲まれた領域1102内に、図9に示した無効領域(無効ROI)903が含まれていないか判定する。
Next, in
最後に、ステップ1004にて、代表値適合性判定部34は、3点の代表点の値を用い、3点に囲まれた領域中を通過する超音波ビーム方向の代表値の勾配を算出する。
Finally, in
ステップ1004の代表値適合性判定部34の操作を図示したものが図11である。図11にて、代表値適合性判定部34は、k1,k2,k3の3点に囲まれた領域1102の内部を通過するビームベクトル1101を設定し、そのビームベクトル1101の方向の代表値の勾配を算出する。
FIG. 11 illustrates the operation of the representative value
次に、ステップ1005にて、代表値適合性判定部34は、ビームベクトル1101の方向の勾配が負であれば、代表値を適正として選出する。
Next, in
ここで、代表値適合性判定部34が、ビームベクトル1101の方向の勾配が負であったときに適正として選出する理由としては、検査対象100を伝搬するにつれ超音波は減衰するものであると想定しており、超音波のビーム方向、つまり音波の伝播方向に沿って必ずエネルギーが減少する、すなわち代表値の勾配は負となる。
Here, the reason why the representative value
ただし、ROI適合性分析部32によって、そもそも輝度分布が適正でないと判断されている無効領域(無効ROI)は必ずしもエネルギーが減少するとは限らないため、ステップ1003にて予め除去している。
However, since the energy of an invalid region (invalid ROI) whose brightness distribution is determined to be inappropriate by the ROI
図2、ステップ209では、ステップ208で選出された代表点間に対し2次元の補間処理を行う。これにより、代表点の分布をより滑らかにすることができ,画像生成部36で生成される画像をより滑らかにすることができる。また,複数のROIを間隔をあけて配置した場合であっても、ROI間の信号値を算出することができるため、処理速度を向上させることができる。なお、補間方法は、2次元の線形補間、2次元のスプライン補間、2次元のBスプライン補間方法等の公知の方法を用いる。
In
以上により、2次元空間で深度方向および方位方向にスペックルを除去した信号(代表値および補間処理値)を生成することができる。 As described above, it is possible to generate a signal (representative value and interpolated value) from which speckles are removed in the depth direction and azimuth direction in a two-dimensional space.
最後に、図2のステップ210にて、補間部35は、画像の生成を行う。画像の生成方法に関しては、実施形態2、実施形態3で述べる。
Finally, in
<<実施形態2>>
図12は、図2のステップ210での、画像生成部36によるB像の画像の生成処理のフローチャートの例を示したものである。
<<Embodiment 2>>
FIG. 12 shows an example of a flowchart of the image generation process of the B image by the
まず、ステップ1201で、画像生成部36には、図2のステップ209から出力された補間処理後の信号が入力される。
First, in
次に、ステップ1202にて、画像生成部36は、音圧プロファイルを読み込む。なお、この音圧プロファイルは、図2のステップ202の音圧プロファイル補正で用いたものと同じものである。
Next, in
次に、ステップ1203にて、画像生成部36は、B像の描画領域全体のゲイン補正を行うための2次元補正ゲインを算出する。2次元補正ゲインG(I,t)は、補間処理後の信号をS(I,t)、音圧プロファイルをP(I,t)とすると、
次に、ステップ1204において、受信RF信号を読み込む。
Next, in
次に、ステップ1205では、ステップ1205で読み込んだ受信RF信号を2次元補正ゲインによって補正し、補正後の信号を作成する。補正後の信号をC(i,t)、受信RF信号をM(i,t)とすると、
ステップ1206では、画像生成部36は、式(4)で作成した補正後の信号C(i,t)を用いてB像を生成することにより、2次元補正ゲインで補正されたB像を得ることができる。画像生成部36は、生成した画像を表示部16に表示する。
In
図13は、実施形態2で作成された補正B像の表示例である。ゲインマップ1302と、補正B像1303が並べて表示されており、ゲインマップ1302中には、図2のステップ208で選出された代表点1202と、図2のステップ204で分析された無効領域(無効ROI)1201が示されている。ユーザは、B像単体はもちろんであるが、図13に示したようにゲインマップと比較することで、補正が適切に行われるか判断することができる。
FIG. 13 is a display example of the corrected B image created in the second embodiment. A gain map 1302 and a corrected B image 1303 are displayed side by side, and the gain map 1302 includes the
<<実施形態3>>
図14は、図2のステップ210での、画像生成部36による減衰率マップの画像の生成処理のフローチャートの例を示したものである。減衰率は、単位距離、単位周波数あたりの超音波エネルギーの減衰量である。
<<Embodiment 3>>
FIG. 14 shows an example of a flowchart of the attenuation rate map image generation process performed by the
実施形態3では、まず、ステップ1401で、画像生成部36には、図2のステップ209から出力された補間処理後の信号が入力される。
In the third embodiment, first, in
次に、ステップ1402で、画像生成部36は、補間処理後信号の時間軸に音速の値(生体内で約1540 m/s)をかけることで距離に換算する。
Next, in
次に、ステップ1403で、画像生成部36は、ビーム方向の勾配を算出する。このビーム方向の勾配は、超音波の伝播に伴う超音波エネルギーの減衰量そのものを示している。
Next, in
さらにステップ1404にて、画像生成部36は、ステップ1402で算出したビーム方向の勾配から、減衰率を算出する。減衰率α[dB/MHz/cm]は、ビーム方向の勾配をD[dB/cm]、送信に用いた周波数をf[MHz]とすると、
画像生成部36は、上記処理を、ユーザが外部入力デバイス13を介して予め設定した減衰率マップ表示領域内の全領域に対して行い、減衰率マップを作成する。画像生成部36は、作成した減衰率マップを表示部16に表示する。
The
図15は、減衰率マップの表示部16における表示形態の例である。超音波画像1500上に、減衰率マップ表示領域1503が重畳され減衰率マップが表示されている。この表示形態では、減衰率マップと、代表点1502、無効領域1501をユーザは比較することができる。よって、ユーザは、代表点や無効領域の選出のされ方から、減衰率の計測が適切に行われているか判断することができる。
FIG. 15 is an example of the display format of the attenuation rate map on the
図16は、減衰率マップの表示部16における表示形態のもう1つの例である。超音波画像1600上に、減衰率マップ表示領域1603が重畳されている。減衰率の大きさは、減衰率の表示レンジ1605で定められたカラーバーに対応するカラーによって示されている。この表示形態では、ユーザがB像1602と減衰率マップとの比較から、無効値1604が適切に検出され、減衰率の計測が適切に行われているか判断することができる。
FIG. 16 is another example of the display format of the attenuation rate map on the
上述してきたように、本実施形態1の超音波診断装置によれば、2次元空間で深度方向および方位方向にスペックルを除去した信号(代表値および補間処理値)を生成することができる。よって、この信号を用いて実施形態2のようにB像を生成してもよいし、実施形態3のように減衰率を算出してもよい。 As described above, the ultrasonic diagnostic apparatus of the first embodiment can generate signals (representative values and interpolated values) from which speckles are removed in the depth direction and azimuth direction in a two-dimensional space. Therefore, using this signal, the B image may be generated as in the second embodiment, or the attenuation rate may be calculated as in the third embodiment.
B像は、超音波の伝搬による減衰が精度よく補正されているため、高精細な画像となっている。よって、ユーザは、この高精細なB像を見て診断を行うことができる。 The B image is a high-definition image because the attenuation due to ultrasound propagation has been accurately corrected. Therefore, the user can make a diagnosis by looking at this high-definition B image.
また。実施形態3のように、減衰率を算出して表示することにより、減衰率が所定値以上に大きい場合は、脂肪肝等の減衰率に特徴が表れる疾患の判定に役立てることができる。 Also. By calculating and displaying the attenuation rate as in the third embodiment, if the attenuation rate is greater than a predetermined value, it can be useful for determining a disease characterized by the attenuation rate, such as fatty liver.
1:超音波送受信装置、10:探触子、13:外部入力デバイス、16:表示部、20:送受信制御部、21:送信ビームフォーマ、22:受信ビームフォーマ、21:送信ビームフォーマ、30:信号処理部、31:2次元ROI設定部、32:ROI適合性分析部、33:代表値算出部、34:代表値適合性判定部、35:補間部、36:画像生成部、37:減衰率算出部、100:検査対象 DESCRIPTION OF SYMBOLS 1: Ultrasonic transmitter/receiver, 10: Probe, 13: External input device, 16: Display section, 20: Transmission/reception control section, 21: Transmission beamformer, 22: Reception beamformer, 21: Transmission beamformer, 30: Signal processing unit, 31: Two-dimensional ROI setting unit, 32: ROI compatibility analysis unit, 33: Representative value calculation unit, 34: Representative value compatibility determination unit, 35: Interpolation unit, 36: Image generation unit, 37: Attenuation Rate calculation unit, 100: Inspection target
Claims (7)
前記信号処理部は、
複数の前記受信信号を前記受信走査線の順に並べた2次元信号空間に2次元のROIを複数設定する2次元ROI設定部と、
前記ROI内の信号を用いて代表値を求め、求めた代表値を前記ROI内に設定した代表点の処理後信号値とする代表値算出部と、
前記代表値算出部が算出した複数の前記ROIについての前記代表値の分布が、予め定めた条件に適合しているかを判定する代表値適合性判定部とを備え、
前記代表値適合性判定部は、前記条件に適合していない代表値を前記処理後信号値として採用しないことを特徴とする超音波診断装置。 a reception beamformer that receives signals from an ultrasound probe and performs reception beamforming to generate reception signals for each of a plurality of reception scanning lines set in an imaging region; and a signal processing unit that processes the reception signals. and has
The signal processing section includes:
a two-dimensional ROI setting unit that sets a plurality of two-dimensional ROIs in a two-dimensional signal space in which the plurality of received signals are arranged in the order of the received scanning lines;
a representative value calculation unit that obtains a representative value using the signal within the ROI and uses the obtained representative value as a processed signal value of a representative point set within the ROI;
a representative value suitability determination unit that determines whether the distribution of the representative values for the plurality of ROIs calculated by the representative value calculation unit conforms to predetermined conditions;
The ultrasound diagnostic apparatus is characterized in that the representative value suitability determination unit does not adopt a representative value that does not meet the conditions as the processed signal value.
前記信号処理部は、
複数の前記受信信号を前記受信走査線の順に並べた2次元信号空間に2次元のROIを複数設定する2次元ROI設定部と、
前記ROI内の信号を用いて代表値を求め、求めた代表値を前記ROI内に設定した代表点の処理後信号値とする代表値算出部と、
前記処理後信号値を用いて、Bモード像を生成する画像生成部とを有することを特徴とする超音波診断装置。 a reception beamformer that receives signals from an ultrasound probe and performs reception beamforming to generate reception signals for each of a plurality of reception scanning lines set in an imaging region; and a signal processing unit that processes the reception signals. and has
The signal processing section includes:
a two-dimensional ROI setting unit that sets a plurality of two-dimensional ROIs in a two-dimensional signal space in which the plurality of received signals are arranged in the order of the received scanning lines;
a representative value calculation unit that obtains a representative value using the signal within the ROI and uses the obtained representative value as a processed signal value of a representative point set within the ROI;
An ultrasonic diagnostic apparatus comprising: an image generation section that generates a B-mode image using the processed signal value.
前記信号処理部は、
複数の前記受信信号を前記受信走査線の順に並べた2次元信号空間に2次元のROIを複数設定する2次元ROI設定部と、
前記ROI内の信号を用いて代表値を求め、求めた代表値を前記ROI内に設定した代表点の処理後信号値とする代表値算出部と、
前記処理後信号値を用いて、Bモード像の輝度値を最適化する画像生成部とを有することを特徴とする超音波診断装置。 a reception beamformer that receives signals from an ultrasound probe and performs reception beamforming to generate reception signals for each of a plurality of reception scanning lines set in an imaging region; and a signal processing unit that processes the reception signals. and has
The signal processing section includes:
a two-dimensional ROI setting unit that sets a plurality of two-dimensional ROIs in a two-dimensional signal space in which the plurality of received signals are arranged in the order of the received scanning lines;
a representative value calculation unit that obtains a representative value using the signal within the ROI and uses the obtained representative value as a processed signal value of a representative point set within the ROI;
An ultrasonic diagnostic apparatus comprising: an image generation unit that optimizes a brightness value of a B-mode image using the processed signal value.
前記信号処理部は、
複数の前記受信信号を前記受信走査線の順に並べた2次元信号空間に2次元のROIを複数設定する2次元ROI設定部と、
前記ROI内の信号を用いて代表値を求め、求めた代表値を前記ROI内に設定した代表点の処理後信号値とする代表値算出部と、
前記2次元信号空間において、複数の前記ROIの前記代表値を用いて、前記ROIとROIの間の前記2次元信号空間の信号値を補間により求める補間部とを有することを特徴とする超音波診断装置。 a reception beamformer that receives signals from an ultrasound probe and performs reception beamforming to generate reception signals for each of a plurality of reception scanning lines set in an imaging region; and a signal processing unit that processes the reception signals. and has
The signal processing section includes:
a two-dimensional ROI setting unit that sets a plurality of two-dimensional ROIs in a two-dimensional signal space in which the plurality of received signals are arranged in the order of the received scanning lines;
a representative value calculation unit that obtains a representative value using the signal within the ROI and uses the obtained representative value as a processed signal value of a representative point set within the ROI;
an interpolation unit that obtains signal values in the two-dimensional signal space between the ROIs by interpolation using the representative values of the plurality of ROIs in the two-dimensional signal space. Diagnostic equipment.
前記ROI内の信号を用いて代表値を求め、求めた代表値を前記ROI内に設定した代表点の処理後信号値とする代表値算出部と、
前記代表値算出部が算出した複数の前記ROIについての前記代表値の分布が、予め定めた条件に適合しているかを判定する代表値適合性判定部とを備え、
前記代表値適合性判定部は、前記条件に適合していない代表値を前記処理後信号値として採用しないことを特徴とする信号処理装置。 a two-dimensional ROI setting unit that sets a plurality of two-dimensional ROIs in a two-dimensional signal space in which received signals after receiving beamforming are arranged in order of receiving scanning lines;
a representative value calculation unit that obtains a representative value using the signal within the ROI and uses the obtained representative value as a processed signal value of a representative point set within the ROI ;
a representative value suitability determination unit that determines whether the distribution of the representative values for the plurality of ROIs calculated by the representative value calculation unit conforms to predetermined conditions;
The signal processing device, wherein the representative value suitability determination unit does not adopt a representative value that does not meet the conditions as the processed signal value .
受信ビームフォーミング後の受信信号を受信走査線の順に並べ、形成された2次元信号空間に2次元のROIを複数設定する2次元ROI設定手段と、
前記ROI内の信号を用いて代表値を求め、求めた代表値を前記ROI内に設定した代表点の処理後信号値とする代表値算出手段と、
前記代表値算出手段が算出した複数の前記ROIについての前記代表値の分布が、予め定めた条件に適合しているかを判定する代表値適合性判定手段として機能させる信号処理プログラムであって、
前記代表値適合性判定手段は、前記条件に適合していない代表値を前記処理後信号値として採用しない
ことを特徴とする信号処理プログラム。 a two-dimensional ROI setting means for arranging the received signals after receiving beamforming by the computer in the order of the received scanning lines, and setting a plurality of two-dimensional ROIs in the formed two-dimensional signal space;
representative value calculating means for determining a representative value using the signal within the ROI, and determining the determined representative value as a processed signal value of a representative point set within the ROI ;
A signal processing program that functions as a representative value suitability determining means for determining whether the distribution of the representative values for the plurality of ROIs calculated by the representative value calculating means conforms to a predetermined condition. ,
The representative value suitability determination means does not adopt a representative value that does not meet the conditions as the processed signal value.
A signal processing program characterized by :
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020066134A JP7432426B2 (en) | 2020-04-01 | 2020-04-01 | Ultrasonic diagnostic equipment, signal processing equipment, and signal processing programs |
CN202110190790.4A CN113491535B (en) | 2020-04-01 | 2021-02-19 | Ultrasonic diagnostic apparatus, signal processing apparatus, and recording medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020066134A JP7432426B2 (en) | 2020-04-01 | 2020-04-01 | Ultrasonic diagnostic equipment, signal processing equipment, and signal processing programs |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021159510A JP2021159510A (en) | 2021-10-11 |
JP7432426B2 true JP7432426B2 (en) | 2024-02-16 |
Family
ID=77995008
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020066134A Active JP7432426B2 (en) | 2020-04-01 | 2020-04-01 | Ultrasonic diagnostic equipment, signal processing equipment, and signal processing programs |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7432426B2 (en) |
CN (1) | CN113491535B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010082230A (en) | 2008-09-30 | 2010-04-15 | Fujifilm Corp | Ultrasonic signal processing apparatus and method |
WO2017068892A1 (en) | 2015-10-22 | 2017-04-27 | 株式会社日立製作所 | Ultrasound diagnostic device and attenuation characteristic measurement method |
WO2017110756A1 (en) | 2015-12-24 | 2017-06-29 | オリンパス株式会社 | Ultrasonic observation device, method for operating ultrasonic observation device, and program for operating ultrasonic observation device |
JP2017158917A (en) | 2016-03-11 | 2017-09-14 | 東芝メディカルシステムズ株式会社 | Medical diagnostic apparatus and medical analysis program |
US20190328364A1 (en) | 2018-04-30 | 2019-10-31 | Esaote Spa | Method and ultrasound system for shear wave elasticity imaging |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101219063B (en) * | 2007-01-12 | 2011-03-23 | 深圳迈瑞生物医疗电子股份有限公司 | B image equalization method and system structure based on two-dimension analysis |
JP5730979B2 (en) * | 2013-11-08 | 2015-06-10 | 日立アロカメディカル株式会社 | Ultrasonic diagnostic apparatus and elasticity evaluation method |
JP6164652B2 (en) * | 2014-05-29 | 2017-07-19 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ultrasonic diagnostic equipment |
JP5948527B1 (en) * | 2014-12-22 | 2016-07-06 | オリンパス株式会社 | Ultrasonic observation apparatus, operation method of ultrasonic observation apparatus, and operation program of ultrasonic observation apparatus |
EP3278737A4 (en) * | 2015-03-31 | 2019-01-02 | Olympus Corporation | Ultrasonic observation apparatus, operating method of ultrasonic observation apparatus, and operating program for ultrasonic observation apparatus |
JP6216736B2 (en) * | 2015-04-08 | 2017-10-18 | 株式会社日立製作所 | Ultrasonic diagnostic apparatus and ultrasonic diagnostic method |
JP2017113145A (en) * | 2015-12-22 | 2017-06-29 | オリンパス株式会社 | Ultrasonic observation device, operation method of ultrasonic observation device, and operation program of ultrasonic observation device |
US11138723B2 (en) * | 2016-07-22 | 2021-10-05 | Canon Medical Systems Corporation | Analyzing apparatus and analyzing method |
-
2020
- 2020-04-01 JP JP2020066134A patent/JP7432426B2/en active Active
-
2021
- 2021-02-19 CN CN202110190790.4A patent/CN113491535B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010082230A (en) | 2008-09-30 | 2010-04-15 | Fujifilm Corp | Ultrasonic signal processing apparatus and method |
WO2017068892A1 (en) | 2015-10-22 | 2017-04-27 | 株式会社日立製作所 | Ultrasound diagnostic device and attenuation characteristic measurement method |
WO2017110756A1 (en) | 2015-12-24 | 2017-06-29 | オリンパス株式会社 | Ultrasonic observation device, method for operating ultrasonic observation device, and program for operating ultrasonic observation device |
JP2017158917A (en) | 2016-03-11 | 2017-09-14 | 東芝メディカルシステムズ株式会社 | Medical diagnostic apparatus and medical analysis program |
US20190328364A1 (en) | 2018-04-30 | 2019-10-31 | Esaote Spa | Method and ultrasound system for shear wave elasticity imaging |
Also Published As
Publication number | Publication date |
---|---|
CN113491535B (en) | 2024-04-26 |
CN113491535A (en) | 2021-10-12 |
JP2021159510A (en) | 2021-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5730979B2 (en) | Ultrasonic diagnostic apparatus and elasticity evaluation method | |
US7985182B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic image acquiring method | |
JP6216736B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic method | |
CN108685596B (en) | Tissue property estimation using ultrasound medical imaging | |
US8475382B2 (en) | Ultrasound diagnostic apparatus and method for tracing movement of tissue | |
WO2014162966A1 (en) | Diagnostic ultrasound apparatus and elasticity evaluation method | |
EP1685799B1 (en) | Ultrasonic diagnostic apparatus and ultrasonic image acquiring method | |
EP1747757B1 (en) | Ultrasound system for displaying an elastic image | |
CN102415902B (en) | Ultrasonic diagnostic apparatus and ultrasonic image processng apparatus | |
US20120108972A1 (en) | Ultrasound diagnostic apparatus and method for tracing movement of tissue | |
EP2423872A1 (en) | Ultrasonic image processing apparatus | |
JP6457107B2 (en) | Ultrasonic diagnostic apparatus and attenuation characteristic measuring method | |
JP4575737B2 (en) | Ultrasonic imaging device | |
EP1867285A1 (en) | An ultrasonic imaging apparatus, an image-processing apparatus, and an ultrasonic image-processing method | |
US11219429B2 (en) | Ultrasound imaging apparatus and controlling method for the same | |
EP2157442B1 (en) | Formation of an elastic image in an ultrasound system | |
US20090069684A1 (en) | Ultrasonic imaging apparatus and a method for generating an ultrasonic image | |
WO2023274763A1 (en) | Ultrasonic quantification of acoustic attenuation coefficient in the presence of elevation aperture blockage | |
JP7432426B2 (en) | Ultrasonic diagnostic equipment, signal processing equipment, and signal processing programs | |
US20240156441A1 (en) | Tissue property estimation with ultrasound medical imaging | |
Loizou et al. | Introduction to ultrasound imaging and speckle noise | |
WO2013002351A1 (en) | Ultrasound diagnostic device and medical image processing device | |
JP2024018506A (en) | Blood flow image forming device and blood flow image forming program | |
CN112294362A (en) | Ultrasonic CT apparatus, image processing method, and recording medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20211013 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230327 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231027 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231114 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240112 |
|
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: 20240130 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240205 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7432426 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 |