CN112105945A - 用于高分辨率宽带谐波成像的超声成像系统 - Google Patents
用于高分辨率宽带谐波成像的超声成像系统 Download PDFInfo
- Publication number
- CN112105945A CN112105945A CN201980027776.3A CN201980027776A CN112105945A CN 112105945 A CN112105945 A CN 112105945A CN 201980027776 A CN201980027776 A CN 201980027776A CN 112105945 A CN112105945 A CN 112105945A
- Authority
- CN
- China
- Prior art keywords
- waveform
- ultrasound system
- transmit
- echo signals
- frequency
- 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
- 238000012285 ultrasound imaging Methods 0.000 title abstract description 5
- 238000003384 imaging method Methods 0.000 title description 18
- 230000004044 response Effects 0.000 claims abstract description 14
- 230000005540 biological transmission Effects 0.000 claims abstract description 10
- 238000002604 ultrasonography Methods 0.000 claims description 40
- 239000000872 buffer Substances 0.000 claims description 8
- 239000000523 sample Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000002592 echocardiography Methods 0.000 abstract description 22
- 238000012545 processing Methods 0.000 abstract description 10
- 230000009467 reduction Effects 0.000 abstract description 6
- 230000035515 penetration Effects 0.000 abstract description 4
- 230000001419 dependent effect Effects 0.000 description 13
- 238000006243 chemical reaction Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 238000000034 method Methods 0.000 description 8
- 239000002131 composite material Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 238000013329 compounding Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 230000001427 coherent effect Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000002961 echo contrast media Substances 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 210000003813 thumb Anatomy 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 210000005166 vasculature Anatomy 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/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
- A61B8/145—Echo-tomography characterised by scanning multiple planes
-
- 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/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- 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
-
- 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
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- 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/895—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum
- G01S15/8952—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum using discrete, multiple frequencies
-
- 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/8959—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
- G01S15/8963—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes using pulse inversion
-
- 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/52019—Details of transmitters
- G01S7/5202—Details of transmitters for pulse systems
-
- 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/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
-
- 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/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
-
- 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/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
-
- 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/52077—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 with means for elimination of unwanted signals, e.g. noise or interference
Abstract
超声成像系统在相同的波束方向上相继发射两个波束，第二波束的波形是第一波形的相位反转形式，并且每个波形包含第一和第二主要频率分量。响应于每个波束发射而接收来自沿着波束方向的回波，并且通过在公共深度基础上以加法方式和减法方式组合两个回波序列来处理所述回波。由这种处理产生的回波分量从针对良好近场分辨率的高频率延伸到针对良好深度穿透性的低频率，并且包括基本、谐波、以及和或差频率分量。线性和非线性回波信号分量可以是为了散斑减少而复合的频率。
Description
技术领域
本发明涉及医学超声成像系统，并且具体涉及执行高分辨率宽带谐波成像的超声系统。
背景技术
谐波成像目前在超声成像中发现两个主要应用。一个是超声造影剂的成像。造影剂的微泡能够被设计为在由超声进行声穿透时振荡或破裂。这种振荡或破坏将使得从微泡返回的回波在谐波频率分量中是丰富的。对接收到的回波进行滤波以将谐波分量与基本频率分量分离。优选分离技术被已知为脉冲反转，如美国专利US 5951478(Hwang等人)中描述的。利用这些回波产生的图像能够锐利地分割包含造影剂的血流和脉管系统。
谐波成像的第二显著应用是组织谐波成像，其中，允许大致正弦发射波形当其穿过身体时经历自然失真。失真导致谐波频率分量的发展，其中，最显著的谐波频率分量通常在基本发射频率的二次谐波处。回波被接收并且以与组织谐波信号类似的方式被处理，以分离从身体返回的谐波分量。在美国专利US 5879303(Averkiou等人)中描述了组织谐波成像。对于由多路径散射引起的其低杂乱水平，从谐波分量产生的图像是期望的。
美国专利US 6440075(Averkiou)描述了增强非线性信号分量的产生的另一超声成像技术。这通过发射具有两个主要频率的波形来完成。当波形穿过组织或遇到微泡时，每个发射频率的谐波分量将如上面描述的那样发展。此外，两个发射频率分量将相互调制，由此发展和和差频率分量。这两种类型的非线性信号被接收，并且被用于形成通过使用两种非线性机制来增强的图像。该专利给出了和(sum)和差(difference)频率能够被形成并且被定位的若干方式的范例，诸如通过针对主要发射频率使用换能器通带的侧面并且针对差频率和谐波频率分量使用中心。
组织谐波成像具有两个显著的限制，其在极近场中成像的能力和其在远场中成像的能力。由于组织谐波成像的谐波分量仅在发射波穿过组织之后开始发展，因此对于谐波回波的返回而言很少或没有能量存在于极近场中。对于身体中的更大深度处的成像(其对于深腹部成像(诸如对肝脏进行成像)常常是必要的)，需要更低频率来抵抗深度相关的频率衰减的效应。由于二次谐波频率依据定义是基本发射频率的两倍，其依据性质是高频的，并且因此由深度相关的衰减更大地影响。该衰减能够降低接收到的回波的信噪比特性并且因此图像的诊断质量。此外，回波信号的谐波分量仅包含在回波中反射的能量的小部分，还减小其信噪比特性。因此期望能够以将在图像的近和远场两者中改进图像的分辨率的方式执行谐波成像。
发明内容
根据本发明的原理，描述了一种通过使用回波信号的宽带宽来执行高分辨率谐波成像的超声系统。发射波形被设计为包含两个主要频率分量。发射波形被发射两次，每次具有不同的发射调制。接收到的来自两个发射的回波通过加法和加法以两种方式组合，以提取图像信息的全谱。这分离了两个主要频率分量的基本、谐波和相互调制的分量。这些不同的分量包括提供良好近场分辨率的高频率相互调制和/或谐波分量、以及提供远场中的良好分辨率的强低频率基本分量。这些不同的频率分量的内容能够通过在信号分量分离之后使用跟踪滤波器来强调。
附图说明
在附图中：
图1以框图形式图示了根据本发明的原理构造的超声诊断成像系统。
图2图示了不同频率的波形的组合以产生多频率发射波形和其反转形式。
图3a和3b图示了产生超声系统中的多频率发射波形的两种方式。
图4图示了可以根据本发明的原理通过将互补回波信号相加和相减而提取的不同频率分量。
具体实施方式
首先参考图1，示出了根据本发明的原理构造的超声诊断成像系统。图1的超声系统利用发射器16，发射器16发射用于正被成像的对象内的谐波频率以及和和差频率信号的生成的多频率波束。发射器通过发射/接收(T/R)开关14被耦合到超声探头10的阵列换能器12的元件，这保护波束形成器免受高电压发射脉冲。发射器对确定发射波束的特性的多个控制参数作出响应，如附图中示出的。多频率发射脉冲的两个主要频率f1和f2被控制，其确定和(f1+f2)和差(f1-f2)频率分量将落在的频率。还控制两个发射的频率分量的幅度或强度a和b，其使得发射脉冲具有形式(asin(2πf1t)+bsin(2πf2t))。然而，由于和和差信号由相互调制的非线性效果产生，接收到的和和差信号分量将具有不是a和b强度的线性积的幅度c。两个发射频率的谐波也通过非线性效果来产生。
在图1中，换能器阵列12接收来自身体的包含不同频率(包括基本频率f1和f2、这两个频率的谐波、以及在换能器通带内的和和差频率分量)的全谱的回波。这些回波信号通过T/R开关14被耦合到波束形成器18，波束形成器18适当地延迟来自不同元件的回波信号，然后对它们进行组合以形成沿着接收到的回波方向的从浅到更深深度的回波信号的序列。优选地，波束形成器是对数字化回波信号操作以产生从场的近到远深度的离散相干数字回波信号的序列的数字波束形成器。波束形成器可以是响应于单个发射脉冲而产生沿着多个空间上不同的接收扫描线的回波信号的两个或更多个序列的多线波束形成器。来自沿着波束方向的相干波束形成回波被存储在线A缓冲器20a中。
在回波信号的线的这种第一接收之后，第二波束在与第一波束相同的波束方向上发射。发射器16具有为发射脉冲的相位或极性的另一可变发射参数，如附图中示出的。超声系统在相同的波束方向上发射不同发射极性、幅度和/或相位的两个或更多个波束。对于双脉冲实施例，响应于第二发射脉冲而接收的回波被存储在线B缓冲器20b中。响应于两个发射脉冲而接收的扫描线回波然后在空间(深度)基础上以通过加法和加法的两种方式被组合。备选地，接收到的回波的第二波束可以直接与存储的第一波束的回波组合而不缓冲。
当两个波束的回波通过求和器或加法器22被组合时，被称为脉冲反转分离的过程被执行。脉冲反转分离非线性信号分量，包括发射脉冲的两个频率的和、差和二次谐波分量，如在下面更充分地描述的。作为发射脉冲的不同的相位或极性的结果，异相基本(线性)回波分量将抵消，并且同相的非线性分量将组合以彼此加强，从而产生增强的且隔离的非线性谐波、和和差频率信号。非线性信号分离能够备选地通过功率调制的差分来执行，如美国专利US 5577505(Brock-Fisher等人)中描述的。根据本发明的原理，接收到的回波也在空间基础上通过减法器24以减法方式组合。减法组合抵消两个回波信号的非线性分量，并且加强了来自两个回波信号的基本线性频率分量。这种减法组合过程的结果也在下面更充分地进行描述。通过加法和减法组合分离的信号分量然后在每个深度位置处通过包括求和器或加法器26的组合器被组合成单个回波信号。
由求和器或加法器26产生的回波信号通过滤波器30来滤波，以移除不期望的信号，诸如带外噪声。优选地，滤波器30是当从场的增加深度接收回波信号时将其通带从高改变到低频率跟踪滤波器。跟踪滤波器能够是具有随着增加的深度而从高频带移动到低频带的通带的带通滤波器，或它能够是具有当从增加的深度接收回波时在频率上被移动更低的高频率切断的带通滤波器。当高频带由于深度相关的衰减而没有有用的或未污染的信号信息时，深度相关的跟踪滤波器通过排除高频率被进一步处理来改进分辨率。
经滤波的回波信号由检测器32检测，检测器32可以是幅度或相位检测器。在图1的图示的实施方式中，当回波要被处理以形成B模式(灰度)图像时，检测器是检测回波信号的包络的幅度检测器。当回波信号要被多普勒处理时，使用相位检测。幅度检测的回波信号由执行诸如扫描转换的操作以在图像场中形成组织的灰度图像的B模式处理器34处理，然后还由图像处理器36处理以显示在超声图像显示器38上。
图2概念地图示了在发射器16的控制下发射的多频率波束如何形成。第一频率的波形80与第二频率的波形82组合。这两个波形的多频率组合被示为复合波形84。诸如波形84的多频率波形被用于由换能器阵列12发射的第一波束。波形80被反转(180°相位改变)以产生第二波形86，第一波束的相位反转复制。波形86被用于由换能器阵列在与第一波束相同的波束方向上发射的第二波束。
图3a和3b详细地图示了适于在本发明的超声系统中使用的发射器16的两种实施方式。在图3a的实施方式中，发射脉冲的不同频率分量f1和f2的波形在数字操作中分离地形成，然后被组合以形成应用于换能器元件的复合多频率发射信号。f1生成器42产生f1发射信号分量，并且f2生成器44产生f2发射信号分量。生成器例如能够算法地生成信号分量。生成器响应于输入控制参数(诸如在附图中示出的选择针对发射波束的f1和f2频率分量的f1Sel.和f2Sel.)而产生其相应的发射波形。其他可变输入参数(未示出)可以是强度或切趾参数a和b、以及针对反转发射信号的相位或极性参数。备选地，由生成器42和44产生的输出波形可以在由求和器46组合成包含多个发射频率分量的复合发射波形之前或之后在幅度和相位或极性上改变。在图3a中，由生成器产生的波形由数字加权处理器电路43和45进行加权，数字加权处理器电路43和45将幅度或切趾加权因子a和b应用于所生成的波形。加权电路能够采取数字乘法器的形式，并且加权因子的符号(+1，-1)能够用于控制输出波形的极性。复合发射波形被应用于D/A转换器48以转换为模拟信号，其可以根据需要被进一步放大并滤波并且被用于驱动换能器阵列的换能器元件12'。
图3b图示了第二发射器实施方式，其中，复合多频率发射信号被事先产生，然后被存储在可以采取数字存储器形式的波形库50中。当特定多频率发射波形被期望时，其通过选择信号Sel.从库50中选择，通过加权电路51进行加权，并且被存储在发射寄存器52中。当发射器被触发以发射波束时，复合多频率波形通过时钟信号CLK被移出发射寄存器52，通过D/A转换器48被转换为模拟信号并且被应用于换能器元件12'。发射脉冲的幅度可以通过在A/D转换器之前的数字乘法器(诸如在加权电路中使用的数字乘法器)或通过在A/D转换器之后的模拟放大器改变，并且可以根据需要在模拟或数字域中被滤波。个体频率分量可以在发射寄存器52之后不在幅度上个体地调节，因为在该实施方式中波形在该点处已经是复合波形。
图4图示了超声系统如何在跨图像场的每个波束位置处使用图2的用于双脉冲发射的发射波形采集回波信号频率的宽普。在附图的顶部部分4a中，超声系统以加法方式组合已经被采集并且被存储在两个线缓冲器20a和20b中的回波信号84e和86e。通过求和器22进行的加法组合分别分离非线性差、二次谐波、以及和信号分量f2-f1、2f1、2f2和f2+f1。求和器22也将产生未在附图中示出的另一发射的频率分量的更高频率的二次谐波2f2。对于在该范例中被选择用于发射的f1和f2分量，非线性分量落在换能器通带中，如求和器22右侧的谱中示出的。为了图示的清楚起见，在附图中仅示出了非线性分量带的中心频率，但是将意识到，每个分量是实际上在中心频率周围的频率的带。应看到谱包括5.8MHz周围的和信号分量f2+f1的高频率、产生良好近场分辨率的高频率处的相互调制。这些分量由在发射之后不久的场中的发射信号分量f1和f2的相互调制产生。附图的中间部分4b示出了接收到的回波84e和86e通过减法器24以减法方式组合，从而实现线性基本频率分量f1和f2的回波信号，如减法器右侧示出的。频率f1是来自更深图像深度的回波中的能够在可检测幅度处返回的低频率，其中，低频率抵抗深度相关的衰减的效果。该信号也是基本频率，并且因此包含从更大深度处的组织反射的大部分能量。该信号还是响应于两个发射而被接收的回波的组合，并且因此，由于减法组合的原因，来自两个回波的分量彼此加强以产生来自更大深度的强接收回波。基本频率回波分量f2类似地从两次接收被加强，但是更易受深度相关的衰减的效果影响。在该范例中，f1和f2频率被审慎地选择，使得f1是被接收的最低频率并且能够进行来自更大深度的可接受返回，而f2频率落在非线性分量2f1与加法组合的f2+f1之间。
在加法和减法组合通过求和器26被组合之后，整个回波谱出现，如附图的底部部分4c中示出的。在该范例中，谱从小于2MHz的低频率延伸到几乎6MHz的高频率。因此，这些回波信号将由于相互调制的和频率分量f2+f1而在极近场中产生良好的高频率分辨率，并且由于强低频率基本信号分量f1而在远场中产生良好的分辨率。当从近到远场处理回波时从更高频率移动到更低频率的跟踪滤波器30(如图1所示)将甚至更多增强信噪比性能，并且强调该宽回波信号谱的益处。
图5以框图形式示出了本发明的超声系统的第二实施方式。图5的系统以与图1的系统相同的方式通过加法和减法组合器22和24来构造。这些组合器中的每个的输出被耦合到两个时变滤波器30a和30b。这些滤波器呈现随着从场的增加的深度接收回波信号而改变的通带特性，并且提供在图像的相应深度处提供最佳成像信号的频率分量的强调。经滤波的线性和非线性信号被耦合到检测器62a和62b，检测器62a和62b执行相应信号分量的幅度(包络)检测。检测到的信号然后通过时变加权电路64a和64b进行加权，时变加权电路64a和64b可以被构造为将所应用的检测到的信号乘以随着时间(深度)而改变的权重的乘法器电路。加权信号然后通过求和电路66进行组合。由于来自减法器24的线性信号和来自加法器22的非线性信号产生图像中的不同散斑图案，相干信号成像系统的伪影(通过求和电路66的两个信号路径的组合)将通过被称为频率复合的过程来实现组合的散斑图案的减少。为了描述用于散斑减少的频率复合的过程和效果，参见美国专利Re.35148(Lizzi等人)。求和电路66的组合的信号然后被耦合到B模式处理器34，以用于图像形成以及如图1的超声系统中的后续图像处理和显示。
图6a和6b分别图示了滤波器30b和30a的时变通带特性的范例。在浅深度处，滤波器30b的通带70被成形为通过由加热组合器22产生的更高频率的非线性分量，所述更高频率的非线性分量是相互调制的和信号分量f1+f2以及二次谐波分量2f1。在浅深度处，这些更高频率分量尚未显著地受深度相关的衰减影响，并且由于更高频率产生良好的分辨率，它们对于图像的浅(近场)深度是优选的。当从增加的深度(例如，中等范围)接收回波信号时，滤波器特性70转变为包括更低差频率分量f2-f1的通带，因为这些更低频率的非线性分量比更高频率的非线性分量更不受增加的深度相关的衰减影响。非线性回波信号分量在图像的中等范围内是强调的。当从场(深)的更大深度接收回波信号时，通带70继续转变为通过更低非线性信号分量f2-f1和2f1的更低通带，以相对排除最高频率的非线性和信号分量。因此，在考虑高频率信号的良好分辨率和高频率的深度相关的衰减的效果时，在每个深度处，通带70裁剪为通过将产生最佳图像的那些信号。
图6b提供了可以被用于线性信号分量的滤波器30a的时变滤波器特性72的范例。在浅深度处，滤波器特性72的通带在更高频率处是更大的，以强调f2频带的更高基本频率。这为图像的近场区域提供了具有良好分辨率的更高频率，并且所发射的f2频率将从近场返回回波信号，其中，非线性信号还未充分形成，从而进一步改进近场图像质量。在中等范围中，对于基本频带f1和f2两者的通过，滤波器特性72变得更平衡。在更深的深度处，滤波器特性72的响应转变到更低带以强调来自最不易受深度相关的衰减的效果影响的低基本频率信号f1的回波。因此，滤波器响应特性(诸如图6b中示出的特性)将利用来自f2频率回波的贡献来增强浅深度处的图像，并且将利用来自f1频带的回波的贡献来增强更深深度处的具有良好穿透性的图像。
图7a提供了通过加权电路64b应用于非线性回波信号分量的权重的变化的范例。曲线74示出了权重在最浅深度处时相对小，然后逐渐增加至中等范围内的更大加权值。在非线性信号分量已经有时间在声场中发展之后并且在深度相关的衰减的效果已经变得过度显著之前，这在中等范围内提供来自非线性回波信号的显著贡献。在更大深度处，曲线74示出了权重随着非线性回波信号的信噪比变得更极大地受深度相关的衰减影响而减小。
图7b图示了通过加权电路64b应用于线性回波信号分量的权重的变化的范例。曲线76示出了被应用在极近场(浅深度)中的更大权重，增加基本频率回波信号对最浅深度处的图像的贡献，其中，非线性信号还未显著发展。曲线76随着从中等范围接收回波以提供来自具有良好分辨率的非线性回波信号的相对更大的贡献而降低(较小加权)。在更大深度处，权重再次增加，如由上升曲线76指示的，其增加具有更好穿透性的更低基本频率回波的贡献，因为深度相关的衰减的效果更大地影响更高频率非线性回波信号。
加权电路的权重的值优选地还考虑两个其他因子。一个是对使线性和非线性信号路径中信号的相对幅度相等以用于更好频率复合的需要。当两个信号路径的信号适当地平衡时，实现散斑图案伪影的更好减少。在这方面，参见美国专利US 5957852(Hossack等人)。另一个是图像数据的对数压缩。用于B模式图像的灰度值通常被转换为对数改变的值，因为已经发现对数转换导致在诊断上更有用的灰度值的范围。尽管对数转换能够在图像信号路径的多种点处被执行，但是将转换函数合并到在加权电路64a和64b中使用的权重的选择中是高效的。因此，加权电路64a和64b能够执行三种功能：不同图像深度处的线性和非线性信号分量的最佳组合的选择；使用于有效频率复合和散斑减少的标称信号幅度相等；以及灰度图像值的对数转换。备选地，对数转换可以通过在每个加权电路之前或之后使用单独的对数转换存储器查找表来实现。
从上述内容看出，本发明的构思在图5的系统中的实施不仅通过线性和非线性回波信号分量两者的选择性组合使用来增强图像的整个深度上的图像分辨率和穿透性，图1的系统也是如此，其额外地提供了由散斑减少增强的图像。
应当注意，适于在本发明的实施方式中使用的超声系统以及特别地图1、3a和3b的超声系统的部件结构可以用硬件、软件或其组合来实施。超声系统的各种实施例和/或部件(例如，发射器、B模式处理器、图像处理器、加权电路、和波形生成器、或其中的部件和控制器)也可以被实施为一个或多个计算机或微处理器的部分。计算机或处理器可以包括计算设备、输入设备、显示单元和例如用于访问因特网的接口。计算机或处理器可以包括微处理器。微处理器可以连接到通信总线，例如，以访问PACS系统或用于输入训练图像的数据网络。计算机或处理器还可以包括存储器。存储器设备(诸如波形库和线缓冲器)可以包括随机存取存储器(RAM)和只读存储器(ROM)。计算机或处理器还可以包括存储设备，其可以是硬盘驱动器或可移除存储驱动器，诸如软盘驱动器、光盘驱动器、固态拇指驱动器等。存储设备也可以是用于将计算机程序或其他指令加载到计算机或处理器中的其他类似模块。
如本文所使用的，术语“计算机”或“模块”或“处理器”或“工作站”可以包括任何基于处理器或基于微处理器的系统，包括使用微控制器、简化指令集计算机(RISC)、ASIC、逻辑电路和能够执行本文所描述的功能的任何其它电路或处理器的系统。上述范例仅是示范性的，并且因此不旨在以任何方式限制这些术语的定义和/或含义。
计算机或处理器执行存储在一个或多个存储元件中的指令集，以便处理输入数据。存储元件还可以根据需要或期望来存储数据或其他信息。存储元件可以采取处理机器内的信息源或物理存储器元件的形式。
如上所述的超声系统的指令集(包括控制超声图像的采集、处理和显示的那些指令)可以包括指示计算机或处理器作为处理机器以执行特定操作(诸如，本发明的各种实施例的方法和过程，并且具体地，如上所述的接收的回波信号的相加和相减)的各种命令。指令集可以采取软件程序的形式。软件可以采取诸如系统软件或应用软件的各种形式，并且其可以被实现为有形和非瞬态计算机可读介质。此外，软件可以采取单独的程序或模块的集合的形式，诸如，波束形成器控制模块、较大程序内的程序模块或程序模块的部分。软件还可以包括面向对象编程形式的模块化编程。由处理机器对输入数据的处理可以响应于操作者命令，或者响应于先前处理的结果，或者响应于由另一处理机器进行的请求。
此外，对以下权利要求的限制不以模块加功能的格式写出，并且不旨在基于35U.S.C.112第六段来进行解释，除非和直到这样的权利要求限制明确地使用短语“用于…的模块”，之后是没有进一步结构的对功能的说明。
Claims (15)
1.一种用于产生高分辨率宽带谐波图像的超声系统，包括：
超声探头，其包括阵列换能器，所述阵列换能器适于发射超声波束并且接收来自沿着所发射的超声波束的方向的回波信号；
发射器，其被耦合到所述超声探头，并且适于控制所述探头在公共波束方向上发射以不同方式调制的第一和第二发射波形，每个波形包括至少两个主要频率分量；
缓冲器，其适于存储响应于至少所述第一发射波形的发射而接收到的回波信号；
求和器，其被耦合到所述缓冲器，并且适于在空间基础上以加法方式组合响应于所述第一发射波形和所述第二发射波形的发射而接收到的回波信号；
减法器，其被耦合到所述缓冲器，并且适于在空间基础上以减法方式组合响应于所述第一发射波形和所述第二发射波形的发射而接收到的回波信号；以及
图像处理器，其适于使用以加法方式和减法方式组合的回波信号来产生超声图像。
2.根据权利要求1所述的超声系统，其中，所述缓冲器还适于存储响应于所述第二发射波形的发射而接收到的回波信号。
3.根据权利要求1所述的超声系统，其中，所述发射器还适于控制所述探头发射第一波形和第二波形，所述第二波形是所述第一波形的反转形式。
4.根据权利要求2所述的超声系统，其中，所述发射器还适于控制所述探头发射第一波形和第二波形，所述第二波形处于与所述第一波形的180°相位关系。
5.根据权利要求1所述的超声系统，其中，所述阵列换能器还适于接收在所述两个分量频率的基本频率、谐波频率、以及和频率或差频率处的回波信号。
6.根据权利要求5所述的超声系统，其中，所述求和器还适于产生谐波频率分量、和频率分量和差频率分量的回波信号。
7.根据权利要求1所述的超声系统，其中，所述减法器还适于产生基本发射频率分量的回波信号。
8.根据权利要求7所述的超声系统，还包括组合器，所述组合器被耦合到所述求和器和所述减法器，并且适于产生基本频率分量、谐波频率分量、以及和频率分量或差频率分量的回波信号。
9.根据权利要求8所述的超声系统，还包括跟踪滤波器，所述跟踪滤波器被耦合到所述组合器和所述图像处理器。
10.根据权利要求1所述的超声系统，其中，所述求和器还适于在输出端处产生非线性回波信号分量，并且所述减法器还适于在输出端处产生线性回波信号分量；并且所述超声系统还包括：
第一时变滤波器，其被耦合到所述求和器的输出端，并且适于当从增加的深度接收非线性回波信号时从高通滤波器特性改变到更低频带滤波器特性；以及
第二时变滤波器，其被耦合到所述减法器的输出端，并且适于当从增加的深度接收线性回波信号时从高通滤波器特性改变到更低频带滤波器特性。
11.根据权利要求10所述的超声系统，其中，所述第二时变滤波器还适于强调高频线性回波信号对浅深度处的图像的贡献，并且强调低频线性回波信号对更深深度处的所述图像的贡献。
12.根据权利要求11所述的超声系统，其中，所述第一时变滤波器还适于强调非线性回波信号对所述中等范围内的所述图像的贡献。
13.根据权利要求10所述的超声系统，还包括：
第一检测器，其被耦合到所述第一时变滤波器的输出端；
第二检测器，其被耦合到所述第二时变滤波器的输出端；以及
求和电路，其被耦合到所述第一检测器和所述第二检测器的输出端，
其中，所述求和电路适于产生散斑减少的回波信号。
14.根据权利要求13所述的超声系统，还包括：
第一加权电路，其被耦合在所述第一检测器与所述求和电路之间，并且适于将根据深度而改变的权重应用于非线性回波信号；以及
第二加权电路，其被耦合在所述第二检测器与所述求和电路之间，并且适于将根据深度而改变的权重应用于线性回波信号。
15.根据权利要求14所述的超声系统，其中，所述第一加权电路还适于将最大权重应用于场的中等深度处的回波信号；并且
其中，所述第二加权电路还适于将最大权重应用于场的浅深度和更深深度处的回波信号。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201862661775P | 2018-04-24 | 2018-04-24 | |
| US62/661775 | 2018-04-24 | ||
| PCT/EP2019/059636 WO2019206709A1 (en) | 2018-04-24 | 2019-04-15 | Ultrasound imaging system for high resolution wideband harmonic imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112105945A true CN112105945A (zh) | 2020-12-18 |
Family
ID=66240116
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201980027776.3A Pending CN112105945A (zh) | 2018-04-24 | 2019-04-15 | 用于高分辨率宽带谐波成像的超声成像系统 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20210077078A1 (zh) |
| CN (1) | CN112105945A (zh) |
| WO (1) | WO2019206709A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113180725A (zh) * | 2021-03-18 | 2021-07-30 | 聚融医疗科技(杭州)有限公司 | 一种基于透镜回波的超声脉冲反转性能评估方法及系统 |
| CN113295777A (zh) * | 2021-04-07 | 2021-08-24 | 聚融医疗科技(杭州)有限公司 | 一种基于透镜回波的提升谐波成像性能的方法及系统 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7280713B2 (ja) * | 2019-02-27 | 2023-05-24 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10179589A (ja) * | 1996-11-08 | 1998-07-07 | Atl Ultrasound Inc | 高調波による超音波画像処理方法および装置 |
| US5879303A (en) * | 1996-09-27 | 1999-03-09 | Atl Ultrasound | Ultrasonic diagnostic imaging of response frequency differing from transmit frequency |
| US5908389A (en) * | 1996-09-27 | 1999-06-01 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging of harmonic frequencies with speckle reduction processing |
| US6251074B1 (en) * | 1996-11-26 | 2001-06-26 | Atl Ultrasound | Ultrasonic tissue harmonic imaging |
| US20010016685A1 (en) * | 2000-02-17 | 2001-08-23 | Jing-Wen Tsao | Ultrasonic diagnosis apparatus |
| US20020040188A1 (en) * | 2000-10-02 | 2002-04-04 | Michalakis Averkiou | Ultrasonic diagnostic imaging of nonlinearly intermodulated and harmonic frequency components |
| US20040015079A1 (en) * | 1999-06-22 | 2004-01-22 | Teratech Corporation | Ultrasound probe with integrated electronics |
| CN107205722A (zh) * | 2015-01-29 | 2017-09-26 | 皇家飞利浦有限公司 | 宽带混合的基波和谐波频率超声诊断成像 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4561019A (en) | 1983-05-16 | 1985-12-24 | Riverside Research Institute | Frequency diversity for image enhancement |
| US5706819A (en) | 1995-10-10 | 1998-01-13 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging with harmonic contrast agents |
| US5577505A (en) | 1996-02-06 | 1996-11-26 | Hewlett-Packard Company | Means for increasing sensitivity in non-linear ultrasound imaging systems |
| US5957852A (en) | 1998-06-02 | 1999-09-28 | Acuson Corporation | Ultrasonic harmonic imaging system and method |
| US11564659B2 (en) * | 2014-12-15 | 2023-01-31 | Canon Medical Systems Corporation | Ultrasonic diagnostic and image processing apparatus for tissue harmonic imaging by extracting nonlinear components from three signals via addition after phase rotation |
-
2019
- 2019-04-15 WO PCT/EP2019/059636 patent/WO2019206709A1/en active Application Filing
- 2019-04-15 US US17/050,235 patent/US20210077078A1/en active Pending
- 2019-04-15 CN CN201980027776.3A patent/CN112105945A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5879303A (en) * | 1996-09-27 | 1999-03-09 | Atl Ultrasound | Ultrasonic diagnostic imaging of response frequency differing from transmit frequency |
| US5908389A (en) * | 1996-09-27 | 1999-06-01 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging of harmonic frequencies with speckle reduction processing |
| JPH10179589A (ja) * | 1996-11-08 | 1998-07-07 | Atl Ultrasound Inc | 高調波による超音波画像処理方法および装置 |
| US6251074B1 (en) * | 1996-11-26 | 2001-06-26 | Atl Ultrasound | Ultrasonic tissue harmonic imaging |
| US20040015079A1 (en) * | 1999-06-22 | 2004-01-22 | Teratech Corporation | Ultrasound probe with integrated electronics |
| US20010016685A1 (en) * | 2000-02-17 | 2001-08-23 | Jing-Wen Tsao | Ultrasonic diagnosis apparatus |
| US20020040188A1 (en) * | 2000-10-02 | 2002-04-04 | Michalakis Averkiou | Ultrasonic diagnostic imaging of nonlinearly intermodulated and harmonic frequency components |
| CN107205722A (zh) * | 2015-01-29 | 2017-09-26 | 皇家飞利浦有限公司 | 宽带混合的基波和谐波频率超声诊断成像 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113180725A (zh) * | 2021-03-18 | 2021-07-30 | 聚融医疗科技(杭州)有限公司 | 一种基于透镜回波的超声脉冲反转性能评估方法及系统 |
| CN113180725B (zh) * | 2021-03-18 | 2023-03-24 | 聚融医疗科技(杭州)有限公司 | 一种基于透镜回波的超声脉冲反转性能评估方法及系统 |
| CN113295777A (zh) * | 2021-04-07 | 2021-08-24 | 聚融医疗科技(杭州)有限公司 | 一种基于透镜回波的提升谐波成像性能的方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2019206709A1 (en) | 2019-10-31 |
| US20210077078A1 (en) | 2021-03-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6656123B2 (en) | Combined fundamental and harmonic ultrasonic imaging at low MI or deeper depths | |
| US8454516B1 (en) | System and method for three dimensional harmonic ultrasound imaging | |
| US6458083B1 (en) | Ultrasonic harmonic imaging with adaptive image formation | |
| JP4433427B2 (ja) | 超音波散乱体をイメージングするためのシステム及び方法 | |
| JP4620261B2 (ja) | 高次の非線形成分による超音波撮像 | |
| EP1088244B1 (en) | Method and apparatus for harmonic imaging using multiple focal zones | |
| US20080249417A1 (en) | Non-Linear Ultrasonic Diagnostic Imaging Using Intermodulation Product Signals | |
| US20080275338A1 (en) | Nonlinear Ultrasonic Diagnostic Imaging Using Intermodulation Product Signals | |
| US6960169B2 (en) | Spread spectrum coding for ultrasound contrast agent imaging | |
| JP6978316B2 (ja) | ブロードバンドの混合された基本波及び高調波周波数超音波診断イメージング | |
| US8932225B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic method | |
| JPH11178824A (ja) | パルス反転ドップラー超音波診断画像処理方法及び装置 | |
| JP2003512914A (ja) | 動き効果を低減させた超音波パルス反転高調波分離 | |
| JP2004512857A (ja) | コントラストパルスシーケンスイメージングを使用する医療診断用超音波システム | |
| US6210334B1 (en) | Medical diagnostic ultrasound method and apparatus for harmonic detection using doppler processing | |
| CN112105945A (zh) | 用于高分辨率宽带谐波成像的超声成像系统 | |
| JP4642977B2 (ja) | 超音波診断装置及び超音波イメージング方法 | |
| US20100036249A1 (en) | Method and apparatus for multiline color flow and angio ultrasound imaging | |
| US6599248B1 (en) | Method and apparatus for ultrasound diagnostic imaging | |
| JP2009022462A (ja) | 超音波診断装置および超音波診断装置の制御プログラム | |
| US20060122510A1 (en) | Ultrasonic diagnostic apparatus and data processing method therefor | |
| US20050124895A1 (en) | Ultrasonic speckle reduction using nonlinear echo combinations | |
| JP4445255B2 (ja) | 組織で発生される高調波との広帯域周波数合成を用いた超音波スペックル低減の方法及び装置 | |
| CN111513755B (zh) | 一种超声成像设备、超声成像方法 | |
| EP3652560B1 (en) | An ultrasound probe and processing method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |