CN110392842A - 用于测量表面运动的方法和声学装置 - Google Patents
用于测量表面运动的方法和声学装置 Download PDFInfo
- Publication number
- CN110392842A CN110392842A CN201780057675.1A CN201780057675A CN110392842A CN 110392842 A CN110392842 A CN 110392842A CN 201780057675 A CN201780057675 A CN 201780057675A CN 110392842 A CN110392842 A CN 110392842A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- wave
- energy converter
- measuring process
- point
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 99
- 230000033001 locomotion Effects 0.000 title claims abstract description 50
- 238000002604 ultrasonography Methods 0.000 claims abstract description 72
- 238000005259 measurement Methods 0.000 claims abstract description 56
- 230000008569 process Effects 0.000 claims abstract description 47
- 238000013508 migration Methods 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims description 22
- 230000004044 response Effects 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 11
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000007493 shaping process Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 238000013507 mapping Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 230000002463 transducing effect Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 210000004209 hair Anatomy 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 206010002091 Anaesthesia Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000037005 anaesthesia Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000000763 evoking effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000001562 sternum Anatomy 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- 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/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- 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
-
- 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/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S15/36—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated with phase comparison between the received signal and the contemporaneously transmitted signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1102—Ballistocardiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
- A61B5/1135—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing by monitoring thoracic expansion
-
- 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/15—Transmission-tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4209—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
-
- 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/4477—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device using several separate ultrasound transducers or probes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Gynecology & Obstetrics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种用于检测表面(21)上多个点(P)的运动的方法，包括测量步骤，在该测量步骤中，将入射超声波朝向表面发射到空气中，并检测在空气中经表面(21)反射得到的超声波。在该测量步骤中，用入射超声波以多个入射角照射每个测量点，并且由多个超声接收换能器(3a)组成的接收换能器网络(3)来检测反射超声波。通过确定所述测量点的两个波束成形信号之间的延迟和/或相移来确定表面在该测量点处的运动。
Description
技术领域
本发明涉及用于检测表面运动的方法和声学装置。
本发明尤其涉及用于检测表面运动的方法，包含如下测量步骤，在所述测量步骤中，使用超声波发射装置将至少一个入射超声波朝向表面发射到空气中，并检测由所述至少一个入射超声波在空气中经所述表面反射得到的至少一个超声波所代表的反射信号。
背景技术
文献US4122427描述了这种方法的示例，其中，通过向测量通道的表面发射频率为40Hz数量级的超声波来测量位于测量通道的表面运动。超声波由单个换能器发射。
发明内容
本发明的目的是进一步改进这种类型的方法，尤其为了实现更高的表面运动检测效率。
为此目的，根据本发明，所述类型的方法在每次测量时，其特征在于：
-用所述至少一个入射超声波以多个入射角照射至少属于所述表面的每个测量点来测量多个测量点的运动，
-由包括多个超声接收换能器的接收换能器网络来检测反射信号，并且通过将所述反射信号至少进行接收的波束成形来确定每个测量点的波束成形信号，
并且特征还在于，还包括至少一个运动确定步骤，在该运动确定步骤中，通过确定所考虑的测量点的两个波束成形信号之间的延迟或相移的至少一个来确定表面在该测量点处的所述运动。
由于这些条件，可以测量各种液体或固体的表面运动，无论光滑还是粗糙，无论平坦与否，也无需考虑它们对光的透明度。表面运动可以在大面积上成像，例如几十个平方厘米，频率高达一千赫兹或更多。对于最小可检测速度为每秒零点几毫米的量级，本发明所述测量方法的灵敏度达到1微米。最后，可使用低声波功率，例如60-70dB SPL(声压级，soundpressure level)的量级。
根据本发明所述方法的优选实施例，还可以使用以下条件中的一个和/或另一个：
-在每个测量步骤中，用所述至少一个入射波照射所述表面的每个测量点，入射角在超过至少20度的入射角范围内扩展；
-在每个测量步骤中，测量在面积大于10cm2的所述表面(21)上的基本上任何点的运动；
-超声波发射装置和接收换能器网络是二维的；
-在基本垂直的两个方向上，发射换能器网络(2)的孔径至少等于超声波接收装置的孔径；
-在至少一个方向上，发射换能器网络的孔径等于超声波接收装置的孔径的至少三倍；
-超声波发射装置的孔径在每个方向上至少等于20cm；
-超声波发射装置包括发射换能器网络，发射换能器网络包括多个超声发射换能器；
-超声发射换能器被分成几组，并且在所述测量步骤中，属于同一组的超声发射换能器同时发射同一信号；
-超声波发射装置包括至少一个超声发射换能器，并被布置成发射到混合腔中，该混合腔适配于将所述至少一个入射超声波在发送到表面之前引起多次反射；
-超声波的频率小于100kHz；
-超声波以每秒500次以上的频率发射；
-在每个测量步骤k中，在至少在表面的不同点处，在所述点处计算接收的波束成形信号Sk，
并且在每个运动确定步骤中，对于两个不同的k值，通过确定在所述点处接收的波束成形信号Sk之间的延迟或相移来确定所述表面上每个点的运动。
-在每个测量步骤k中，超声波发射装置向所述表面发射非聚焦的入射超声波；
-在每个测量步骤k中，超声波发射装置连续地向所述表面的不同点发射聚焦的入射超声波，并且由对应于在所述点处的聚焦的入射超声波的反射信号来确定在所述点处接收的波束成形信号Sk；
-接收的波束成形信号Sk由以下公式确定：
其中：
-rj是由指数为j的超声接收换能器(3a)检测到的信号，
-t是时间，
-dj是点P与指数为j的超声接收换能器(3a)之间的距离，
-c是超声波在空气中的速度。
-超声波发射装置包括发射换能器网络，发射换能器网络包括多个超声发射换能器，
在每个测量步骤k中，测量每个超声发射换能器和每个超声接收换能器之间的各自脉冲响应，然后至少在所述表面的不同点处，计算所述点处发射和接收的波束成形信号S'k，
并且在每个运动确定步骤中，对于两个不同的k值，通过确定在所述点处发射和接收的波束成形信号S'k之间的延迟或相移来确定所述表面的每个点的运动。
-在所述点处每个发射和接收的波束成形信号由以下公式确定：
其中：
-i是1～M之间的指数，指的是超声发射换能器，
-j是1～N之间的指数，指的是超声接收换能器，
-hijk是指数为i的超声发射换能器和指数为j的超声接收换能器之间的脉冲响应，
-t是时间，
-dij是超声波从指数为i的超声发射换能器经所述表面上的所考虑的点(P)处反射，到达指数为j的超声接收换能器而行进的距离，
-c是超声波在空气中的速度；
-在至少一些测量步骤中，确定预定观察区域内每个测量点的波束成形信号，并且确定属于所述表面的测量点为那些将所述波束成形信号最大化的测量点；
-在运动确定步骤中，确定同一测量点处的两个测量步骤之间的波束成形信号的延迟dt，以及：
-与dt·c成正比的测量点的行程δ
-和/或与dt·c/Δt成正比的测量点的速度，
其中c是超声波在空气中的速度，Δt是所述两个测量步骤之间的时间间隔；
-在运动确定步骤中，确定同一测量点处的两个测量步骤之间的波束成形信号的相移以及：
-与成正比的测量点的行程δ
-和/或与成正比的测量点的速度，
其中c是超声波在空气中的速度，f是超声波的频率，Δt是所述两个测量步骤之间的时间间隔。
此外，本发明还涉及一种用于检测反射超声波的表面的运动的装置，该装置包括超声波发射装置，超声波接收装置，控制超声波发射装置和接收由超声波接收装置接收的信号的控制装置，其中控制装置适配于执行若干连续的测量步骤，在每个测量步骤中，超声波发射装置将至少一个入射超声波朝向所述表面发射到空气中，并且由超声波接收装置接收由所述至少一个入射超声波在空气中经所述表面反射得到的至少一个超声波所代表的反射信号。
其特征在于：
超声波发射装置适配于由所述至少一个入射超声波以多个入射角照射至少属于所述表面(21)的多个测量点(P)，其特征还在于：
超声波接收装置是包括多个超声接收换能器的接收换能器网络，其特征还在于：
控制装置适配于在每个测量步骤中通过至少将所述反射信号进行接收的波束成形来确定每个测量点的波束成形信号，并且其特征还在于:
控制装置适配于通过确定所考虑的测量点的两个波束成形信号之间的延迟或相移的至少一个来确定该测量点处的所述表面运动。
本发明的其他特征和优势将在以下对其实施例之一的参考附图的描述中显示，实施例作为示例，但并不局限于此。
附图说明
在附图中：
-图1是实施本发明的超声波发生装置示例的示意图，
-图2示出了可用于本发明的发射换能器网络的示例，
-图3是示出用于一个或多个超声换能器的混合腔示例的剖面的示意图，
-图4是图3所示混合腔的内部透视图；
-图5是示出了组装在载体上的发射换能器网络和接收换能器网络的透视示意图，
-图6是由图1至5的装置测量表面行进速度的映射的示例图，
-图7是由图1至5的装置测量表面行程的映射的示例图，
-图8是由图1至5的装置测量表面上的点随时间的行进速度曲线的示例图，以及
-图9是由图1至5的装置测量表面上的点随时间的行程曲线的示例图。
在不同的附图中，相同的标记指代相同或相似的元件。
具体实施方式
图1示出了用于测量表面21的运动的超声装置1，表面21可以是任何种类，不论是固体或者液体。例如，表面21可以是人或动物的皮肤。
超声装置1包括发射换能器网络2和接收换能器网络3，发射换能器网络2可以包括一个或多个(M个)超声发射换能器2a(TE1，TE2，......TEM)，接收换能器网络3可以包括多个(N个)超声接收换能器3a(TR1，TR2，......TRN)。
换能器网络2和3都可以优选为二维网络。
换能器网络2和3都可以布置在同一载体上，并且在这种情况下，换能器2a和3a可以互相间隔布置。
换能器2a和3a可以是任何已知类型。在一个实施例中，超声发射换能器2a可以是高带宽扬声器，超声接收换能器3a可以是高带宽麦克风。
如下所述，超声发射换能器2a可选地布置成在将超声波发送到所述表面21之前将超声波发射到混合腔4中。
发射换能器网络2和接收换能器网络3旨在测量属于所述表面21的多个测量点P的运动，例如测量在面积超过10cm2以及可能几十cm2的所述表面21上的基本上所有点的运动。
为此目的，发射换能器网络2和接收换能器网络3各自设计成具有相对大的孔径，使得所述表面21的每个测量点P，被来自发射换能器网络2的入射超声波以大角度入射角照射，并且使得接收换能器3a以大角度视角观察所述表面21上的每个点P。
发射换能器网络2的孔径可以用网络的横向尺寸OEY、OEZ来定义，例如沿着垂直Z轴和横向水平Y轴(参见图1，2，5)。该孔径对应于从每个测量点P观察发射换能器网络2的发射孔径角α。因此，每个测量点P由入射超声波以多个入射角照射，形成宽度为α的入射角范围。例如，孔径角α可以至少等于20度。
接收换能器网络3的孔径可以用网络的ORY'、ORZ横向尺寸来定义，例如沿垂直轴Z和水平横向轴Y'(参见图1，5)。该孔径对应于从每个测量点P观察接收换能器网络3的接收孔径角β。因此，每个测量点P由超声接收换能器3a以多个视角观察，形成宽度为β的视角范围。例如，接收孔径角β可以至少等于20度。
发射孔径角α可以优选为大于接收孔径角β。
特别地，发射换能器网络2沿轴Y、Z的至少一条或两条轴的孔径，优选大于接收换能器网络沿相应轴Y'、Z的孔径，这确保了至少部分反射超声波总是返回到接收换能器网络3。
可选地，发射换能器网络2沿轴Y、Z的至少一条或两条轴的孔径，可以是接收换能器网络沿着相应的轴Y'、Z的孔径的3倍以上(或甚至5倍，或甚至7倍)。
超声装置1可以由电子电路11控制，电子电路11又可以连接到计算机12或类似设备。
电子电路11可包括：
-连接到每个发射换能器2a(E1-EM)的数模转换器6(D/A1-D/AM)，如有必要，该数模转换器6与缓冲存储器9(B1-BN)相关联，
-由计算机12控制的电子中央单元5(CPU)，用来控制每个模数转换器6，
-至少一个连接到电子中央单元5的中央存储器7(MEM)，
-可能地，特定的信号处理电路8(DSP)
-连接到每个接收换能器3a(TR1-TRN)、并分别与缓冲存储器9(B'1-B'N)相关联的模数转换器10(A/D1-A/DN)。
可以将超声发射换能器2a发射的入射超声波频率控制在100kHz以下，例如在40～70kHz之间。
优选地，电子中央单元5适于以大于每秒500次的频率发射入射超声波，使得可以随时间追踪所述表面21的运动。
如图2在特定实施例中所示，发射换能器网络2可以是由刚性板13承载的超声发射换能器2a的二维网络。
可以控制这些超声发射换能器2a同时或顺序地发射入射超声波。
如果合适的话，超声发射换能器2a可以分成几组，并且控制同一组的所有超声发射换能器2a同时发射相同的信号。这导致了倍增网络，人为地增加了发射换能器网络2的发射孔径。
替代地或另外地，如上所述并且如图1和3所示，发射换能器网络2可以与混合腔4耦合。在这种情况下，发射换能器网络2可以缩减到单个超声发射换能器2a，因为此时发射换能器网络2的孔径OE'Y、OE'Z不是由换能器的空间分布限定，而是由混合腔的发射面尺寸限定，入射超声波通过该发射面朝向所述表面21发射到空气中。
混合腔4可以是例如板或任何其他的固体物体，或者它可以是安装有发射换能器网络2的壳体，如图3所示。在这种情况下，混合腔4可以具有底部14、侧壁15、和平行于底部14并且通过侧壁15保持远离底部14的发射壁。发射换能器网络2可以安装在混合腔4中的底部14上。发射壁16可以具有孔17，使得入射超声波能够输出到表面21。这些孔17构成尽可能多的超声波源，并限定发射换能器网络2的孔径。
如图4所示，反射壁18可以设置在混合腔内，例如附接到发射壁16，优选为基本垂直于底部14。反射壁18以各种方式取向，以促进超声波在混合腔4中的多次反射。
上述装置用于测量表面21的运动，包括任何点的行程和速度。
通过处理由超声接收换能器3a接收的反射信号来执行该测量，信号处理方法可以例如为以下面描述的三种方法之一。
方法1
在这第一信号处理方法中，通过同时激起M个超声发射换能器2a来实现入射超声波(入射超声波发射)对表面21的照射。
在每个对应于入射波发射的测量步骤k中，由超声接收换能器3a接收的信号通常由接收的波束成形来处理。
因此，对于处于观察区域(即可能被表面21占据的区域)的每个点P，接收处的波束成形信号Sk(t)计算如下：
其中：
-rj是接收网络3中换能器Rj接收的信号，
-t是时间，
-dj是点P和换能器Rj之间的距离，
-c是超声波在空气中的速度。
当点P属于反射的表面21时这个信号是最大的。
因此可以确定所述表面21的形貌，进而确定所述表面21的图像，尽管该图像通常是近似的。
另一方面，通过比较对应于不同连续发射的信号Sk(t)，可以非常精确地确定在两次发射k和k+1之间，所述表面上每个点P处的行程δ。
对于所述表面21上的给定点P，在至少一个运动确定步骤中比较连续信号Sk(t)和Sk+1(t)，其示出了相移相移与P在发射k和k+1之间在垂直于所述表面21方向行进的行程δ相关，关系式为：
其中：
-c是超声波在空气中的速度，
-f是超声波的频率，
-θ是超声波在点P处入射和反射的平均角度，即表面法线与反射波束之间的平均角度(根据表面21的形貌确定)。
考虑到表面21对超声波反射，表面21的行程d相当于位于表面21后面的虚源(virtual source)的行程2d·cosθ，因此得到上述公式(1)。
因此，通过以下公式计算出点P在发射k和k+1之间的行程δ：
以类似的方式，可以相对于Sk(t)(t由每次入射波发射进行计数)测量信号Sk+1(t)的延迟的不是相移通过以下公式计算发射k和k+1之间的行程δ：
在这两种情况下，可以容易地追踪表面21的正常速度：V＝δ/Δt，其中Δt是发射k和k+1之间的时间间隔。
方法1非常快速且特别有效，特别是在以宽角度范围照射表面21、而因此在发射换能器网络2具有宽发射孔径的情况下。
因此在这种情况下，发射换能器网络2特别适于采用倍增网络或混合腔。
方法2
在这第二信号处理方法中，混合腔对于发射换能器网络2不是优选的。
每个超声发射换能器2a连续发射脉冲波，发射被认为是准时的，并且发射的是球面波。由所有超声发射换能器2a进行的每一连续发射系列构成测量步骤k，然后由所有超声发射换能器2a在连续发射的后续系列k+1中进行重复发射。
对于每个发射系列k，确定每个超声发射换能器Ei和每个超声接收换能器TRj之间的脉冲响应hijk(t)。
对于待观察区域上的给定点P，其从发射器Ei到接收器Rj的超声波覆盖距离(记作dij)已知。且每一对(Ei，Rj)的距离都已知。
然后为超声接收换能器TRj计算点P发射和接收的波束成形信号。信号rjk(t)通过对所有延迟的脉冲响应求和获得，以便补偿行进时间dij/c。
该信号表示如下：
对应于点P和发射系列k、对于所有超声接收换能器3a的信号S'k(t)是所有接收器的波束成形信号的总和：
如果P是表面21上的质点(material point)，那么信号S将是振幅最大的短脉冲，因为所有接收器计算的所有脉冲将具有相同的延迟。否则，信号将在时间上更加错列并且振幅更低。
该方法使得可以通过第一方式确定表面21的形貌，从而确定表面21的(通常不是非常精确的)图像。
在至少一个运动确定步骤中，对于表面上的相同点P，采用对应于不同测量步骤k的连续信号S'k(t)，可以非常精确地确定点P在两个连续发射系列k和k+1之间的行程δ。
在本发明的所有实施例中，可以在执行测量步骤时执行运动确定步骤。当然，这种计算也可以在延迟的基础上进行。
如果点P在t和t+dt之间向垂直于表面21的方向运动了距离δ，则新信号S'k+1(t)将是相对于S'k(t)、来自dt＝2δcosθ(符号定义同上)的时间延迟脉冲(相对于入射波发射计数)。
时移dt给出了表面21在点P处的行程：
δ＝dt·c/(2cosθ) (5)
可以容易地回溯表面的法向速度：V＝δ/Δt，其中Δt是发射系列k和k+1之间的时间间隔。
如在方法1中那样，由信号S'k(t)和S'k+1(t)之间的相移可以以等效的方式确定表面21的行程和速度。
方法2非常精确，可以测量振幅远低于波长的表面运动。然后可以使用抛物线或更高阶插值来计算距离声学信号采样周期最近部分的时移。
方法3
在第三方法中，发射换能器网络2可以是倍增孔径类型或具有混合腔4。
发射波在表面21上反射并由N个超声接收换能器3a收集。从表面21到接收器的波路径将被认为是直线，而没有任何反射。
方法3包括至少一个校正步骤(一个预备步骤，有重新校正需要时可以按定期间隔进行重复)，在校正步骤时确定超声发射换能器TEi和超声接收换能器TRj之间的脉冲响应hij(t)。在傅立叶变换之后，这些脉冲响应被标记为Hij(省略角脉动ω以简化计算)。如上述方法2中那样，该预备步骤可以通过使每个超声发射换能器TEi连续发射脉冲来执行。
该预备步骤之后是若干连续的测量步骤，以指数k进行标记。在每个测量步骤k中，入射波连续聚焦在表面21的不同点P上，例如这些点P在预备步骤中确定。
这些连续的聚焦操作如下进行。
根据时间反演法和可逆性原理，使发射器i发射反相时间脉冲响应(即在频域Hij *中，其中上标*表示共轭复数)，应获得以下响应：
Rij＝HijHij*＝|Hij|2
Rij的解是实数，对应于以初始时间为中心的脉冲。因此根据该观察，通过简单地发射时间上反向和延迟的脉冲响应，可以使得以一定的延迟向任何一个接收器j发送脉冲。
实际上，可以从单个发射器i向所有接收器发射任何波阵面。该波阵面仅采用在接收器j处的到达延迟(记作Tj)来表示。为此，在超声发射换能器TEi重新发射信号之前组合信号就已足够：
其中Ei是超声发射换能器TEi发射的、在频域中的信号。
因此，可以合成到达超声接收换能器TRj的球面波阵面，该球面波阵面由表面21上的点P构成的发散中心产生。换句话说，入射超声波可聚焦于表面21上的点P处。
当所有超声发射换能器TEi发射合成相同波阵面所需的信号时，该波阵面合成得到进一步增强。所有超声发射换能器TEi发射的信号可以以下述矩阵表示：
在超声接收换能器TRj处，接收的信号是：
可以以矩阵形式总结如下：
由各个超声接收换能器TRj接收的所有信号可以在聚焦于表面21的点P之后进行组合，方法为在对接收的信号加和之前补偿其延迟来实行波束成形：
指数k代表测量步骤，即连续地聚焦在表面21的各个考虑点P上的入射波发射系列。
在矩阵表示法中，发射时在点P处的聚焦和接收时的波束成形被简化为以下单个信号：
Sk＝T R＝T H E＝T H HT*TT* (11)
上标T*表示矩阵是共轭转置的。
信号Sk主要表示在点P处收集的信息。
当该点在测量步骤k和测量步骤k+1之间运动时，结果是Sk和Sk+1之间的相移
该相移使得可以在至少一个运动确定步骤中，使用上述公式(2)回溯表面21在点P处的行程δ，或者如方法1中所解释的回溯行进速度V。
应当注意，上述第三方法的计算也可以在时域而不是在频域中进行。相反，方法1和2的计算也可以在频域中进行。
在设想的不同方法中，根据应用，δ和V的值可以是与上述公式成比例的值。
可选方案
接收的波束成形可由混频器实现，如以下文献所述：
-Nicolas Quieffin,Stefan Catheline,Ros Kiri Ing and Mathias Fink,Acoustic source localization model using in-skull reverberation and timereversal,Applied Physics Letters vol.90,063902(2007)；
-N.Etaix,M.Fink and R.K.Ing,Acoustic imaging device with onetransducer,J.Acoust.Soc.Am.131(5),pp.EL395-EL399,2012；
-N.Etaix,J.Dubois,M.Fink and R.K.Ing,Increasing the modal density inplates for mono-element focusing in air,J.Acoust.Soc.Am.,Vol.134(2),pp.1049-1054,2013.
图6和图7示出了根据本发明的方法测量人体下胸骨的表面21的行进速度和行程的映射示例，所述表面21由心跳引发振动。图8和9示出了根据本发明的方法测量相同固体表面21的行进速度和行程随时间的变化。
因此，本发明特别适用于人体表面运动的分析，特别适用于呼吸运动或心血管系统的研究。
本发明特别有利，尤其是在肺通气的研究中，因为它可以追踪肺容量的变化，从而无需干扰受试者就能追溯吸入/呼出空气的流速，这与例如涉及到口罩或面罩的方法不同。
本发明还可用于自动监测患者的呼吸，例如在重症监护中，从麻醉中醒来，在急诊室的候诊室中，或者是处于危险中的婴儿等。
本发明还可用于医疗领域之外，例如快速检测人的入睡，检测人的压力情况等。
最后，在医疗领域之外，本发明还可用于检测固定在皮肤上和被组织隐藏的物体。通过激发主体(具体来说是观察区域)，皮肤的表面运动将被隐藏的对象打断，这可以通过本申请文件要求保护的装置观察到。
Claims (22)
1.一种用于检测反射超声波的表面(21)的运动的方法，包括若干连续的测量步骤，在每个测量步骤中，使用超声波发射装置(2)将至少一个入射超声波朝向所述表面(21)发射到空气中，并检测由所述至少一个入射超声波在空气中经所述表面(21)反射得到的至少一个超声波所代表的反射信号，
在每个测量步骤中，其特征在于：
-至少属于所述表面的多个测量点(P)的运动由所述至少一个入射超声波以多个入射角照射每个测量点(P)来测量，
-由多个超声接收换能器(3a)组成的接收换能器网络(3)来检测所述反射信号，并且通过至少将所述反射信号进行接收的波束成形来确定每个测量点(P)的波束成形信号，
并且其特征还在于，还包括至少一个运动确定步骤，在所述运动确定步骤中，通过确定所考虑的测量点(P)的两个波束成形信号之间的延迟或相移的至少一个，来确定所述表面(21)在该所考虑的测量点(P)处的所述运动。
2.根据权利要求1所述的方法，其特征在于：在每个测量步骤中，用所述至少一个入射波、以在至少20度的入射角范围内扩展的入射角照射所述表面(21)的每个测量点(P)。
3.根据权利要求1或2所述的方法，其特征在于：在每个测量步骤中，测量在面积大于10cm2的所述表面(21)上的基本上任何点的运动。
4.根据前述权利要求中任一项所述的方法，其特征在于：所述超声波发射装置(2)和所述接收换能器网络(3)是二维的。
5.根据权利要求4所述的方法，其特征在于：所述发射换能器网络(2)的孔径(OEY，OEZ，OE'Y，OE'Z)在两个基本垂直的方向上，至少等于所述超声波接收装置(3)的孔径(ORY'，ORZ)。
6.根据权利要求5所述的方法，其特征在于：所述发射换能器网络(2)的孔径(OEY，OEZ，OE'Y，OE'Z)在至少在一个方向上，等于所述超声波接收装置(3)的孔径(ORY，ORZ)的至少三倍。
7.根据权利要求4至6中任一项所述的方法，其特征在于：所述超声波发射装置(2)的孔径(OEY，OEZ，OE'Y，OE'Z)在每个方向上至少等于20cm。
8.根据前述权利要求中任一项所述的方法，其特征在于：所述超声波发射装置(2)包括发射换能器网络(2)，所述发射换能器(2)包括多个超声发射换能器(2a)。
9.根据权利要求8所述的方法，其特征在于：所述超声发射换能器(2a)被分成几组，并且在所述测量步骤中，由属于同一组的超声发射换能器(2a)同时发射相同的信号。
10.根据前述权利要求中任一项所述的方法，其特征在于：所述超声波发射装置(2)包括至少一个超声发射换能器(2a)，所述超声发射换能器(2a)被布置成发射到混合腔(4)中，所述混合腔(4)适配于将所述至少一个入射超声波在发送到所述表面(21)之前引起多次反射。
11.根据前述权利要求中任一项所述的方法，其特征在于：所述超声波的频率小于100kHz。
12.根据前述权利要求中任一项所述的方法，其特征在于：所述超声波以高于每秒500次的频率发射。
13.根据前述权利要求中任一项所述的方法，其特征在于：在每个测量步骤k中，在至少在所述表面(21)的不同点(P)处的所述点(P)处，计算接收的波束成形信号Sk，
并且在每个运动确定步骤中，对于两个不同的k值，通过确定在所述点(P)处接收的波束成形信号Sk之间的延迟或相移来确定所述表面(21)上的每个点(P)的运动。
14.根据权利要求13所述的方法，其特征在于：在每个测量步骤k中，所述超声波发射装置(2)朝向所述表面(21)发射非聚焦的入射超声波。
15.根据权利要求13所述的方法，其特征在于：在每个测量步骤k中，所述超声波发射装置(2)连续地向所述表面(21)的不同点(P)发射聚焦的入射超声波，并且通过在所述点(P)处与所述聚焦的入射超声波相关的所述反射信号来确定在所述点(P)处接收的所述波束成形信号Sk。
16.根据权利要求13至15中任一项所述的方法，其特征在于：所述接收的波束成形信号Sk由以下公式确定：
其中：
-rj是由指数为j的超声接收换能器(3a)检测到的信号，
-t是时间，
-dj是点P与指数为j的超声接收换能器(3a)之间的距离，
-c是超声波在空气中的速度。
17.根据权利要求1至12中任一项所述的方法，其特征在于：所述超声波发射装置(2)包括发射换能器网络(2)，所述发射换能器网络(2)包括多个超声发射换能器(2a)，
在每个测量步骤k中，在每个超声发射换能器(2a)和每个超声接收换能器(3a)之间测量各自的脉冲响应，然后至少在所述表面(21)的不同点(P)处，计算所述点(P)发射和接收的波束成形信号S'k，
并且在每个运动确定步骤中，对于两个不同的k值，通过确定在所述点(P)的发射和接收的波束成形信号S'k之间的延迟或相移来确定所述表面(21)的每个点(P)的运动。
18.根据权利要求17所述的方法，其特征在于：所述点(P)的每个发射和接收的波束成形信号由以下公式确定：
其中：
-i是1和M之间的指数，指的是超声发射换能器(2a)，
-j是1和N之间的指数，指的是超声接收换能器(3a)，
-hijk是指数为i的超声发射换能器(2a)和指数为j的超声接收换能器(3a)之间的脉冲响应，
-t是时间，
-dij是超声波从指数为i的超声发射换能器(2a)经所述表面(21)的考虑点(P)反射到达指数为j的超声接收换能器(3a)而行进的距离，
-c是超声波在空气中的速度。
19.根据前述权利要求中任一项所述的方法，其特征在于：在至少一些测量步骤中，为预定观察区域的每个测量点(P)确定波束成形信号，并且确定属于表面(21)的测量点(P)为那些将所述波束成形信号最大化的测量点(P)。
20.根据前述权利要求中任一项所述的方法，其特征在于：在所述运动确定步骤中，确定同一测量点(P)处的两个测量步骤之间的波束成形信号延迟dt，并且：
-所述测量点(P)的行程δ与dt·c成正比
-和/或所述测量点(P)的速度与dt·c/Δt成正比，
其中c是超声波在空气中的速度，Δt是所述两个测量步骤之间的时间间隔。
21.根据权利要求1至19中任一项所述的方法，其特征在于：在所述运动确定步骤中，确定同一测量点(P)处的两个测量步骤之间的波束成形信号相移并且：
-所述测量点(P)的行程δ与成正比
-和/或所述测量点(P)的速度与成正比，
其中c是超声波在空气中的速度，f是超声波的频率，Δt是所述两个测量步骤之间的时间间隔。
22.一种用于检测反射超声波的表面(21)的运动的装置，包括超声波发射装置(2)，超声波接收装置(3)，控制超声波发射装置(2)和接收由超声波接收装置(3)检测到的信号的控制装置(5)，其中所述控制装置(5)适配于执行若干连续的测量步骤，在每个测量步骤中，所述超声波发射装置(2)将至少一个入射超声波朝向所述表面(21)发射到空气中，并且由所述超声波接收装置(3)接收由所述至少一个入射超声波在空气中经所述表面(21)反射得到的至少一个超声波所代表的反射信号，其特征在于：
所述超声波发射装置(2)适配于由所述至少一个入射超声波以多个入射角照射至少属于所述表面(21)的多个测量点(P)，其特征还在于：
所述超声波接收装置(3)是由多个超声接收换能器(3a)组成的接收换能器网络(3)，其特征还在于：
控制装置(5)适配于在每个测量步骤中通过将所述反射信号至少进行接收的波束成形来确定每个测量点(P)的波束成形信号，
并且控制装置(5)适配于通过确定所考虑的测量点(P)的两个波束成形信号之间的延迟或相移的至少一个来确定所述表面(21)在所述所考虑的测量点(P)处的所述运动。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1656935A FR3054325B1 (fr) | 2016-07-21 | 2016-07-21 | Procede et dispositif acoustique pour mesurer des mouvements de surface. |
| FR1656935 | 2016-07-21 | ||
| PCT/FR2017/051883 WO2018015638A1 (fr) | 2016-07-21 | 2017-07-10 | Procédé et dispositif acoustique pour mesurer des mouvements de surface. |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110392842A true CN110392842A (zh) | 2019-10-29 |
| CN110392842B CN110392842B (zh) | 2023-04-04 |
Family
ID=57137094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201780057675.1A Active CN110392842B (zh) | 2016-07-21 | 2017-07-10 | 用于测量表面运动的方法和声学装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US11249187B2 (zh) |
| EP (1) | EP3488261B1 (zh) |
| CN (1) | CN110392842B (zh) |
| ES (1) | ES2836267T3 (zh) |
| FR (1) | FR3054325B1 (zh) |
| WO (1) | WO2018015638A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| IL267774B2 (en) * | 2019-07-01 | 2023-11-01 | Omer Eshet | Ultrasonic system and method for monitoring vital indicators |
| FR3121830A1 (fr) | 2021-04-16 | 2022-10-21 | Centre National De La Recherche Scientifique | Méthode de caractérisation de la vibration d’une surface |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4625555A (en) * | 1984-02-07 | 1986-12-02 | Terumo Kabushiki Kaisha | Ultrasonic measurement method, and apparatus therefor |
| CN1343310A (zh) * | 1999-03-15 | 2002-04-03 | 弹性测量体系弹性推动公司 | 使用切变波的成像方法和装置 |
| CN101936756A (zh) * | 2010-08-31 | 2011-01-05 | 华南理工大学 | 一种多频相控阵超声多普勒流量检测系统及方法 |
| CN102027440A (zh) * | 2008-03-18 | 2011-04-20 | 艾利普提克实验室股份有限公司 | 对象与运动检测 |
| US20110208060A1 (en) * | 2010-02-24 | 2011-08-25 | Haase Wayne C | Non-contact Biometric Monitor |
| CN102549451A (zh) * | 2009-07-20 | 2012-07-04 | 罗伯特·博世有限公司 | 超声测量装置和用于分析超声信号的方法 |
| CN102590814A (zh) * | 2012-03-02 | 2012-07-18 | 华南理工大学 | 一种超声波探头空间位置和三维姿态的检测装置及方法 |
| US20130301391A1 (en) * | 2010-11-16 | 2013-11-14 | Qualcomm Incorporated | System and method for object position estimation based on ultrasonic reflected signals |
| CN104105449A (zh) * | 2011-12-01 | 2014-10-15 | 毛伊图像公司 | 使用基于声脉冲和多孔多普勒超声的运动检测 |
| CN105518451A (zh) * | 2013-07-22 | 2016-04-20 | 爱莎莉 | 用于表征超声波在线性可变厚度的波导中的传播的超声方法和设备 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4122427A (en) * | 1976-06-24 | 1978-10-24 | Herbert Karsh | Motion monitor |
| US4197856A (en) * | 1978-04-10 | 1980-04-15 | Northrop Robert B | Ultrasonic respiration/convulsion monitoring apparatus and method for its use |
| US8550998B2 (en) * | 2008-01-09 | 2013-10-08 | Bjørn A. J. Angelsen | Nonlinear elastic wave measurement and imaging with two-frequency elastic wave pulse complexes |
| WO2009088306A2 (en) * | 2008-01-09 | 2009-07-16 | Angelsen Bjoern A J | Nonlinear elastic imaging with two-frequency elastic pulse complexes |
| US9939413B2 (en) * | 2008-01-09 | 2018-04-10 | Surf Technology As | Measurement and imaging of scatterers with memory of scatterer parameters using at least two-frequency elastic wave pulse complexes |
| CN101836110B (zh) * | 2008-03-31 | 2013-05-01 | 日本克劳特克雷默尔株式会社 | 超声波探伤方法及其装置 |
| GB2463257A (en) * | 2008-09-04 | 2010-03-10 | Elio Poggiagliolmi | Stabilized platform for an ultrasonic Doppler vibrometer for sensing motion of a surface |
| JP5468589B2 (ja) * | 2011-09-27 | 2014-04-09 | 富士フイルム株式会社 | 超音波診断装置および超音波画像生成方法 |
| JP5702326B2 (ja) * | 2012-03-30 | 2015-04-15 | 富士フイルム株式会社 | 超音波プローブおよびそれを備える超音波診断装置 |
| KR20150118493A (ko) * | 2014-04-14 | 2015-10-22 | 삼성전자주식회사 | 초음파 장치 및 그 제어 방법 |
| KR20170060852A (ko) * | 2015-11-25 | 2017-06-02 | 삼성메디슨 주식회사 | 초음파 영상을 제공하는 방법 및 이를 위한 초음파 장치 |
| JP6933038B2 (ja) * | 2017-08-10 | 2021-09-08 | コニカミノルタ株式会社 | 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法 |
-
2016
- 2016-07-21 FR FR1656935A patent/FR3054325B1/fr not_active Expired - Fee Related
-
2017
- 2017-07-10 WO PCT/FR2017/051883 patent/WO2018015638A1/fr unknown
- 2017-07-10 CN CN201780057675.1A patent/CN110392842B/zh active Active
- 2017-07-10 US US16/319,445 patent/US11249187B2/en active Active
- 2017-07-10 EP EP17748529.9A patent/EP3488261B1/fr active Active
- 2017-07-10 ES ES17748529T patent/ES2836267T3/es active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4625555A (en) * | 1984-02-07 | 1986-12-02 | Terumo Kabushiki Kaisha | Ultrasonic measurement method, and apparatus therefor |
| CN1343310A (zh) * | 1999-03-15 | 2002-04-03 | 弹性测量体系弹性推动公司 | 使用切变波的成像方法和装置 |
| CN102027440A (zh) * | 2008-03-18 | 2011-04-20 | 艾利普提克实验室股份有限公司 | 对象与运动检测 |
| CN102549451A (zh) * | 2009-07-20 | 2012-07-04 | 罗伯特·博世有限公司 | 超声测量装置和用于分析超声信号的方法 |
| US20110208060A1 (en) * | 2010-02-24 | 2011-08-25 | Haase Wayne C | Non-contact Biometric Monitor |
| CN101936756A (zh) * | 2010-08-31 | 2011-01-05 | 华南理工大学 | 一种多频相控阵超声多普勒流量检测系统及方法 |
| US20130301391A1 (en) * | 2010-11-16 | 2013-11-14 | Qualcomm Incorporated | System and method for object position estimation based on ultrasonic reflected signals |
| CN104105449A (zh) * | 2011-12-01 | 2014-10-15 | 毛伊图像公司 | 使用基于声脉冲和多孔多普勒超声的运动检测 |
| CN102590814A (zh) * | 2012-03-02 | 2012-07-18 | 华南理工大学 | 一种超声波探头空间位置和三维姿态的检测装置及方法 |
| CN105518451A (zh) * | 2013-07-22 | 2016-04-20 | 爱莎莉 | 用于表征超声波在线性可变厚度的波导中的传播的超声方法和设备 |
Non-Patent Citations (1)
| Title |
|---|
| 孔涛等: "空气耦合超声换能器声场计算与测量研究", 《机械工程学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EP3488261B1 (fr) | 2020-08-12 |
| FR3054325B1 (fr) | 2019-08-16 |
| US11249187B2 (en) | 2022-02-15 |
| EP3488261A1 (fr) | 2019-05-29 |
| CN110392842B (zh) | 2023-04-04 |
| WO2018015638A1 (fr) | 2018-01-25 |
| FR3054325A1 (fr) | 2018-01-26 |
| ES2836267T3 (es) | 2021-06-24 |
| US20190331793A1 (en) | 2019-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210085292A1 (en) | Ultrasound imaging with sparse array probes | |
| US7470232B2 (en) | Method and apparatus for non-invasive ultrasonic fetal heart rate monitoring | |
| US7066888B2 (en) | Method and apparatus for determining an ultrasound fluid flow centerline | |
| Koblitz et al. | Asymmetry and dynamics of a narrow sonar beam in an echolocating harbor porpoise | |
| US11744553B2 (en) | Ultrasound system, method and computer program product | |
| KR101501479B1 (ko) | 초음파 최적화 방법과 그를 위한 초음파 의료 장치 | |
| Au et al. | Dolphin biosonar signals measured at extreme off-axis angles: Insights to sound propagation in the head | |
| JP2003164453A (ja) | 三次元ドップラー超音波撮像システム及び方法 | |
| Sandrin et al. | Ultrafast two-dimensional ultrasonic speckle velocimetry: A tool in flow imaging | |
| JPH0120899B2 (zh) | ||
| CN102551791A (zh) | 一种超声成像方法和装置 | |
| CN106901777A (zh) | 一种多功能超声探头与微血管成像及血流速度检测方法 | |
| EP2766738B1 (en) | Three dimensional (3d) transverse oscillation vector velocity ultrasound imaging | |
| CN110392842A (zh) | 用于测量表面运动的方法和声学装置 | |
| Bera et al. | Three-dimensional beamforming combining micro-beamformed RF datasets | |
| JPH01145043A (ja) | 超音波診断装置 | |
| CN205697831U (zh) | 一种带有a超和/或m超功能的胎儿监护装置 | |
| Leech et al. | Physical principles of ultrasound | |
| JP2856471B2 (ja) | 超音波診断装置 | |
| Sannnou et al. | 1P5-16 Examination of effectiveness of signal-to-noise ratio factor in estimation of sound speed of medium | |
| US20200237341A1 (en) | Medical Ultrasound system using transducer with circular formation elements placement and built-in MEMS motion sensors | |
| Nambu et al. | The frontal 3-D imaging for PTCA using ultrasonic ring array probe | |
| Leech | Physical Principles and the Basic Exam | |
| JPS60249947A (ja) | 超音波ｃｔ装置 | |
| JPS62204731A (ja) | 超音波減衰音速測定装置 |
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 | ||
| CB02 | Change of applicant information |
Address after: Paris France Applicant after: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE Applicant after: ECOLE SUPERIEURE DE PHYSIQUE ET DE CHIMIE INDUSTRIELLES DE LA VILLE DE PARIS Applicant after: SORBONNE UNIVERSITE Applicant after: INSERM (INSTITUT NATIONAL DE LA SANTE ET DE LA RECHERCHE MEDICALE) Applicant after: University of Western dais, Paris Address before: Paris France Applicant before: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE Applicant before: ECOLE SUPERIEURE DE PHYSIQUE ET DE CHIMIE INDUSTRIELLES DE LA VILLE DE PARIS Applicant before: SORBONNE UNIVERSITE Applicant before: INSERM (INSTITUT NATIONAL DE LA SANTE ET DE LA RECHERCHE MEDICALE) Applicant before: University of Paris |
|
| CB02 | Change of applicant information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20220701 Address after: Paris France Applicant after: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE Applicant after: ECOLE SUPERIEURE DE PHYSIQUE ET DE CHIMIE INDUSTRIELLES DE LA VILLE DE PARIS Applicant after: SORBONNE UNIVERSITE Applicant after: INSERM (INSTITUT NATIONAL DE LA SANTE ET DE LA RECHERCHE MEDICALE) Applicant after: University of Paris Address before: Paris France Applicant before: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE Applicant before: ECOLE SUPERIEURE DE PHYSIQUE ET DE CHIMIE INDUSTRIELLES DE LA VILLE DE PARIS Applicant before: SORBONNE UNIVERSITE Applicant before: INSERM (INSTITUT NATIONAL DE LA SANTE ET DE LA RECHERCHE MEDICALE) Applicant before: UNIVERSITE PARIS DIDEROT - PARIS 7 |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240207 Address after: Paris France Patentee after: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE Country or region after: France Patentee after: ECOLE SUPERIEURE DE PHYSIQUE ET DE CHIMIE INDUSTRIELLES DE LA VILLE DE PARIS Country or region after: Fa Guo Patentee after: SORBONNE UNIVERSITE Patentee after: INSERM (INSTITUT NATIONAL DE LA SANTE ET DE LA RECHERCHE MEDICALE) Patentee after: University of Western dais, Paris Patentee after: ASSISTANCE PUBLIQUE - HOPITAUX DE PARIS Address before: Paris France Patentee before: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE Country or region before: France Patentee before: ECOLE SUPERIEURE DE PHYSIQUE ET DE CHIMIE INDUSTRIELLES DE LA VILLE DE PARIS Country or region before: Fa Guo Patentee before: SORBONNE UNIVERSITE Patentee before: INSERM (INSTITUT NATIONAL DE LA SANTE ET DE LA RECHERCHE MEDICALE) Patentee before: University of Western dais, Paris |
|
| TR01 | Transfer of patent right |