CN103261909A - 用于成像非线性相互作用散射的方法 - Google Patents
用于成像非线性相互作用散射的方法 Download PDFInfo
- Publication number
- CN103261909A CN103261909A CN2011800507044A CN201180050704A CN103261909A CN 103261909 A CN103261909 A CN 103261909A CN 2011800507044 A CN2011800507044 A CN 2011800507044A CN 201180050704 A CN201180050704 A CN 201180050704A CN 103261909 A CN103261909 A CN 103261909A
- Authority
- CN
- China
- Prior art keywords
- signal
- measurement
- scattering
- image curve
- nonlinear interaction
- 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
- 230000003993 interaction Effects 0.000 title claims abstract description 137
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000003384 imaging method Methods 0.000 title claims description 53
- 238000005259 measurement Methods 0.000 claims abstract description 91
- 239000000463 material Substances 0.000 claims abstract description 51
- 230000008859 change Effects 0.000 claims description 42
- 230000035945 sensitivity Effects 0.000 claims description 28
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 11
- 230000002123 temporal effect Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 8
- 238000003491 array Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 4
- 230000009022 nonlinear effect Effects 0.000 claims description 4
- ASNHGEVAWNWCRQ-UHFFFAOYSA-N 4-(hydroxymethyl)oxolane-2,3,4-triol Chemical compound OCC1(O)COC(O)C1O ASNHGEVAWNWCRQ-UHFFFAOYSA-N 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 230000009021 linear effect Effects 0.000 description 44
- 238000001228 spectrum Methods 0.000 description 36
- 238000006073 displacement reaction Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 6
- 230000001902 propagating effect Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 5
- 230000005764 inhibitory process Effects 0.000 description 5
- 230000000644 propagated effect Effects 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 4
- 238000002059 diagnostic imaging Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 241000251468 Actinopterygii Species 0.000 description 2
- 208000004434 Calcinosis Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000000295 emission spectrum Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000004304 visual acuity Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000931526 Acer campestre Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 210000004712 air sac Anatomy 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000009514 concussion Effects 0.000 description 1
- 239000002872 contrast media Substances 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
- 230000006378 damage Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000008713 feedback mechanism Effects 0.000 description 1
- 239000004833 fish glue Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000010895 photoacoustic effect Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 239000007779 soft material Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- 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
-
- 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/96—Sonar systems specially adapted for specific applications for locating fish
-
- 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/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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Nonlinear Science (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
第一和第二脉冲波(105,105)沿着第一和第二发射波束(101,102)被发射,其中至少一个波束在至少一个方向上是宽的,并且第一和第二脉冲波之间的发射定时被选择成使得脉冲波前重叠在沿着材料物体中的至少一个测量或像曲线Γ(r)(107)传播的重叠区R(r,t)(106)。至少接收并处理由重叠区(106)中第一和第二波之间的非线性相互作用产生的散射信号,从而沿着Γ(r)形成非线性相互作用散射图像信号。通过改变第一和第二脉冲波(104,105)之间的相对发射定时或改变第一和第二发射波束(101,102)的至少一个的方向或两者,能够横向扫描测量或像曲线Γ(r)(107)。该方法适用于为电磁波和弹性波及其组合成像非线性散射源。
Description
技术领域
本发明解决了使用弹性和电磁波及其组合成像非线性散射。其应用在医学和技术领域两者中。
背景技术
用于电磁(EM)波和弹性(EL)波两者的材料性质通常示出非线性特性,其中材料参数取决于波的场变量的振幅。非线性材料特性的空间改变提供了EM和EL波两者的非线性散射,这种非线性散射源的成像在许多情况中对识别材料性质有用处。EM和EL波的正向波传播和局部散射两者都具有数学相似性,并且因此用于成像的方法和仪器具有相似的结构。EL波的使用的示例是用横波/剪切波和压缩波两者的材料试验、用压缩波的超声医学成像、和海底声纳和地质测量。EM波具有类似使用,其中尤其为提供二维信息给超声图像的医学成像研发波长在0.1-1mm范围的高GHz和THz范围中的EM技术。在红外线和光频范围中的EM成像也为材料试验和医学成像两者提供了有用信息。
非线性散射能够用于被分为参数的和共振散射类型的EM和EL波两者。对于EL波,参数散射来源于随局部弹性波场的振幅的局部弹性参数的非线性改变,其中非线性改变的空间改变产生非线性散射。对于EM波,参数散射来源于随局部EM波场的振幅的局部介电常数或磁导率的非线性改变,其中非线性改变的空间改变产生非线性散射。用被称为声波的弹性压缩波,例如用从软质肿瘤组织中微型钙化的超声非线性散射发现的,或用从像矿井的土层中的硬质物体或其他物体的声散射发现的,一个在软质材料和硬质材料之间的界面上得到强烈的非线性参数散射。例如,如用从血液中的微气泡或水中鱼鳔等等填充的气体的超声扩散得到的,在较硬材料和更加软的材料之间的界面得到强烈的非线性散射。
用单频率波段入射波,该参数非线性散射在散射波中产生入射频率波段的谐波分量。用局部相互作用的双波段入射波,该参数非线性散射围绕入射频率波段的回旋的波段,其提供围绕入射频率的和和差的波段。然而,材料参数的非线性改变也产生正向传播波的聚集的非线性失真。当高频脉冲的脉冲宽度增加超过低频脉冲的近似一半周期时,来自非线性正向传播失真的线性散射具有与局部非线性散射类似的子波,并且在这个情形中很难辨别从入射波的非线性传播失真的线性散射发生的信号分量和从局部非线性散射发生的信号分量。本发明以方法和仪器的方式呈现解决方法,这方法和仪器抑制来源于由非线性正向传播失真产生的分量的强线性散射的分量,并且提取局部非线性散射部件,从而产生局部非线性散射源的空间成像。
共振非线性散射具有相关的时间滞后,在一些情形中其能够被用于将来自入射波的局部非线性散射和正向传播失真的信号分量。然而,本发明进一步提供用于成像局部共振图像非线性散射源的优点。
对于声波,微气泡示出共振散射,例如其中该共振来源于具有壳层和气体的气泡的非线性弹性和体积近似气泡体积三倍的围绕气泡的共摆动流体质量之间的能量交换。由于弹性和质量随气泡压缩改变,所以共振频率受入射声波场的非线性影响,在散射场中产生具有大量入射频率的谐波分量(入射频率的n倍)以及甚至入射频率的子谐波分量(一小部分入射频率)的极其强的非线性散射。然而,为以大大高于气泡共振频率的频率成像,非线性散射是非常低的,并且本发明为以超过共振频率的频率加强成像微气泡提供了解决方法。
在波场和原子和分子之间的相互作用中发生共振非线性EM散射,这在量子物理学领域中最好地描述。EM共振散射的示例是与子谐波声散射具有相似性的荧光。双光子量子散射类似于第二谐波参数散射,但是包括随着过程中的时间滞后的详细动力。
在材料的EM和EL波中也得到非线性相互作用,其中例如在所谓的声光效应的过程中EL压缩波改变EM材料参数。材料中EM波的吸收产生材料的辐射力和局部加热,这在所谓的光声效应的过程中生成声波。因此,本发明解决了EM和EL波两者及其组合,其中本说明书和权利要求中涉及的波能够是EM和/或EL波。
发明内容
本发明内容给出本发明的内容的简要概述,并且不存在关于本发明范围的任何限制,其中本发明仅由附属权利要求限定。
根据本发明,提供了用于沿着实物/材料物体(material object)中至少一个测量或像曲线T(r)测量或成像两个波之间的非线性相互作用散射的方法,其中非线性材料参数产生波散射,该方法包括的步骤有a)在至少一个发射事件中发射,情况如下i)至少一个第一发射脉冲波,其沿着具有第一发射方向的第一发射波束具有第一发射时间和第一发射中心频率,和ii)至少一个第二发射脉冲波,其沿着具有第二发射方向的第二发射波束具有第二发射时间和第二发射中心频率,其中第一和第二发射方向和第一和第二发射时间被设置成事的来自第一和第二脉冲波的脉冲空间重叠在重叠区R(r,t)中,其沿着测量或像曲线Γ(r)在实物中随时间传播,和b)选择下列的一个或两者i)所述第一和第二发射时间之间的时间关系,和ii)所述第一和第二发射方向的至少一个,从而控制实物内的测量或像曲线Γ(r)的位置,并且c)用至少一个接收器检测从重叠区R(r,t)散射的波,并且生成至少一个接收信号,该信号包括至少非线性相互作用散射信号分量,这些分量涉及在所述重叠区R(r,t)中从第一和第二脉冲波之间的非线性相互作用散射的波,和d)处理所接收的信号,从而提取非线性相互作用散射信号分量并从而相应于所述测量或像曲线Γ(r)形成非线性相互作用测量或图像信号。
优选地,提取非线性相互作用信号分量的过程包括接收信号的时间频率过滤,其中所述过滤通过所述非线性相互作用散射信号分量并且阻断其他时间频率分量。
在优选实施例中,a)所述发射步骤包括至少两个发射事件,并且所述第一和第二发射方向和所述第一和第二发射时间之间的时间滞后对于所述至少两个发射事件是相同的,b)所述第一和第二脉冲波的一个在所述至少两个发射事件之间极性、振幅、相位和频率的至少一个改变,包括所述第一和第二脉冲波的至少一个对于至少一个所述发射事件具有零振幅,并且c)提取非线性信号分量的过程包括组合/结合来自所述至少两个发射事件的接收信号。
优选地在来自所述至少两个发射事件的信号在步骤c)中结合之前,来自所述至少两个发射事件的至少一个的接收信号经下列一个或两个操作,i)延迟修正和ii)脉冲失真修正。
在优选的实施例中,a)通过选择下列一项或两项在实物中选择多个测量或像曲线Γk(r),k=1...K,i)所述第一和第二发射时间之间的时间关系,和ii)所述第一和第二发射方向的至少一个,和b)如上所述获得相应于每个具体的测量或像曲线Γ(r)的测量或图像信号。
优选地,所述发射步骤包括发射多个发射事件,其中所述发射事件之间的发射时间足够长,以便至少两个不同的重叠区在物体内同时传播,并且所述第一和第二发射脉冲和所述第一和第二发射方向之间的相对正时在发射事件之间改变,以便为了测量或成像的改进率,所述至少两个不同的重叠区沿着至少两个不同的测量或像曲线传播,并且并行地接收和处理来自所述不同重叠区的散射信号。
在优选的实施例中,a)所述第一和第二脉冲波的一个的脉冲宽度比所述第一和第二波的另一个的时间振荡周期的一半短,以便非线性相互作用重叠区被分为高灵敏度的子区域,其中所述第一和第二波的所述另一个的波场的数值/大小(magnitude)足够高,以便并行地从多于一个测量或像曲线的一组Γk(r),k=1...K获得具有高灵敏度的散射信号,和b)并行地/同时(in parallel)接收和处理来自多于一个测量或像曲线的接收信号,从而为所述多于一个测量或像曲线同时获得图像信号。
优选地用于同一测量或像曲线Γ(r)的来自多个发射事件的接收信号被平均从而减少噪声并改进灵敏度。
优选地通过在至少两个测量或像曲线上的图像信号之间插值获得用于在所述测量或像曲线Γk(r)之间的图像点的图像信号。
优选地所述测量或像曲线Γk(r)跨表面空间分布,从而为跨所述表面散射的非线性相互作用的二维图像获得图像信号。
优选地所述测量或像曲线Γk(r)在三维区域内空间分布,从而为在所述区域中散射的非线性相互作用的三维图像获得图像信号。
根据另一方面,本发明提供了用于测量或成像在沿着实物/材料物体中至少一个测量或像曲线Γ(r)的两个波之间非线性相互作用散射的仪器,其中非线性材料参数产生波散射,该仪器包括a)发射装置,其经设置从而发射由下列组成的至少一种发射事件i)沿着具有第一发射方向的第一发射波束具有第一发射时间和第一发射中心频率的至少一个第一发射脉冲波,和ii)沿着具有第二发射方向的第二发射波束具有第二发射时间和第二发射中心频率的至少一个第二发射脉冲波,其中发射装置被设置成使得所述第一和第二波束在实物中的重叠区R(r,t)中重叠,其中所述重叠区在实物中沿着测量或像曲线T(r)随时间传播,b)控制装置,其经设置从而选择下列的一个或两个i)第一和第二发射时间之间的时间关系,和ii)所述第一和第二发射方向的至少一个,从而控制测量或像曲线Γ(r)在实物内的空间中的位置,c)接收装置,其被设置成至少从所述重叠区R(r,t)中所述第一和第二脉冲波之间的非线性相互作用接收来自实物的散射波,并且经设置成产生至少一个接收信号,该信号包括来自至少装置重叠区R(r,t)的至少非线性相互作用散射的分量,和d)处理装置,经设置从而处理所述接收信号,从而从所述重叠区R(r,t)提取非线性相互作用散射信号分量,从而沿着至少一个测量或像曲线Γ(r)形成非线性相互作用测量或图像信号。
优选地,所述处理装置包括用于接收信号的时间频率过滤的装置,其中所述过滤通过所述非线性相互作用散射信号分量并且阻断其他时间频率分量。
在优选的实施例中,a)所述发射装置经设置从而为每个测量或像曲线Γ(r)发射至少两个发射事件,b)所述发射装置包括用于在所述至少两个发射事件之间改变所述第一和第二脉冲波的一个的极性、振幅、相位和频率的至少一个的装置,其允许所述第一和第二脉冲波的至少一个为所述至少两个发射事件的至少一个具有零振幅,和c)所述处理装置包括用于在沿着所述测量或像曲线Γ(r)形成图像信号的过程中组合来自至少两个发射事件的接收信号的装置。
优选地所述处理装置包括在结合来自至少两个发射事件的信号之前,用i)延迟修正和ii)脉冲失真修正的一个或两个修正来自所述至少两个发射事件的至少一个的至少一个接收信号。
优选地,所述处理装置包括用于估算i)延迟修正和ii)脉冲失真修正的一个或两个的装置。
优选地所述用于接收的装置包括聚焦接收器***。优选地,所述聚焦接收器***被设置成随时间动态聚焦,以便该焦点在重叠区R(r,t)随时间移动时跟随重叠区R(r,t)。优选地所述聚焦接收器***是聚焦摄像机。
在优选的实施例中,所述聚焦接收器***包括a)多个接收器元件的阵列,其经设置从而并行地/同时(in parallel)接收散射信号,和b)从多组接收器元件并行地/同时生成多个接收信号的装置,以便不同的接收信号被并行地/同时分配给不同的测量或像曲线Γm(r),m=1...M。
优选地所述控制装置包括发射多组发射事件的装置和用于将来自多组发射事件的所述不同的接收信号平均的装置,从而产生不同的平均接收信号并将所述不同的平均接收信号分配给所述不同的测量或像曲线。
一组发射事件可以是单个发射事件或多个发射事件,其提供多个测量或像曲线。如其他地方更详细所述,单个发射事件可以引起多个的子区域的高灵敏度,其依次引起多个测量或像曲线。类似地,多个发射事件(其每个都可能引起或不引起多个子区域)可以引起多个平行测量或像曲线。然后,该组发射事件可以被发射多次,以便从多个测量或像曲线的每个上升的信号能够被平均(即为每个测量或像曲线产生平均的信号)。
优选地所述控制装置包括发射多组发射事件从而为每个测量或像曲线Γ(r)接收多个接收信号的装置和用于处理所述多个接收信号的装置,其中所述处理通过下列的至少一项i)将接收信号平均从而增加信噪比由此增加灵敏度,和ii)比较所述多个接收信号之间的差值,从而确定散射体的移动。
优选地所述控制装置包括这样的装置,其通过改变下列的一项或两项i)所述第一和第二发射时间之间的时间关系和ii)所述第一和第二发射方向的至少一个,从而在实物的区域中通过一组不同的测量或像曲线Γk(r),k=1...K扫描所述测量或像曲线Γ(r)。
优选地所述用于处理的装置包括用于在邻近测量或像曲线Γk(r)之间插值图像信号的装置,从而为所述邻近测量或像曲线之间的图像点生成图像信号。
在一些优选的实施例中,所述控制装置包括用于跨通过所述实物的表面扫描所述测量或像曲线的装置,从而为跨所述表面的非线性相互作用散射的二维图像提供图像信号。
所述控制装置可以包括用于在所述实物的三维区域内扫描所述测量或像曲线的装置,从而为所述区域中的非线性相互作用散射的三维图像提供图像信号。
优选地所述第一和第二发射波束的至少一个在至少方位方向上是宽的。更优选地,所述第一和第二发射波束的另一个在至少方位方向上时窄的,并且接收装置包括轴沿着方位方向上窄的发射波束的轴的聚焦接收波束,并且扫描测量或像曲线是通过扫描接收波束和方位方向上是窄的发射波束两者获得的。
优选地所述发射装置和接收装置包括元件阵列,其中相同阵列表面能够被用于发射所述第一和第二脉冲波两者并且至少接收非线性相互作用散射信号。
优选地,所述发射装置包括发射所述第一脉冲波的元件阵列和发射所述第二脉冲波的另一个元件阵列。
优选地,所述发射装置包含环形元件阵列。
所述第一和第二脉冲波的至少一个可以是弹性压缩波。所述第一和第二脉冲波的至少一个可以是弹性横波。所述第一和第二脉冲电磁波的至少一个可以是电磁波。
该方法和仪器适用于为电磁和弹性波及其组合的非线性材料参数成像。沿着第一和第二发射波束发射第一和第二脉冲波,其中在优选的实施例中至少一个波束在至少一个方向上是宽的。选择所述第一和第二脉冲波之间的发射定时,以便波脉冲在沿着实物中的测量或像曲线Γ(r)传播的重叠区R(r,t)中重叠。至少接收并处理在重叠区中所述第一和第二波之间非线性相互作用散射产生的散射信号,从而沿着Γ(r)形成非线性相互作用散射图像信号。通过改变第一和第二脉冲波之间的相对发射定时或改变第一和第二发射波束的方向或两者,可以二维或三维横向扫描测量或像曲线Γ(r),从而形成物体的二维或三维图像。该方法增强了散射体的测量或成像,这些散射体比周围介质更软或更硬,并且该方法例如对医学成像中的微钙化和造影剂微气泡成像、土层中像矿井的硬式结构的成像、潜水员血液中的气泡形成的成像、和水中的鱼等的鱼鳔的成像有用,但不限制于此。
附图说明
现在将仅通过示例并参考附图描述本发明的某些优选的实施例,其中:
图1a示出沿着测量或像曲线Γ(r)获得非线性相互作用测量或图像信号的基本原理;
图1b示出沿着测量或像曲线Γ(r)的重叠区R(r,t)的有效传播速度的计算;
图2a和2b示出入射第一和第二发射时间频率波段的和和差频率波段的形成和通过在时间频率域中过滤的非线性相互作用散射信号的提取;
图2c示出当第一中心频率小于第二发射波的时间频率带宽时,入射第一和第二发射时间频率波段的和和差频率波段的形成;
图3示出当一个脉冲的频率比另一个脉冲的带宽低很多时非线性相互作用源项被如何分解为子区域;
图4示出从一组第一和发射时间和第一和第二发射波束方向上升的灵敏度子区域的传播;
图5示出第一和第二发射时间之间的发射时间差的改变能够被怎样用于将灵敏度子区域移动到先前的低灵敏度区域;
图6示出根据本发明的仪器的方框图;
图7示出用于检测超声波的非线性相互作用散射信号的第一和第二波发射装置和接收装置的示例实施例;
图8示出用于检测超声波的非线性相互作用散射信号的第一和第二波发射装置和接收装置的另一个示例实施例;
图9示出用于检测超声波的非线性相互作用散射信号的第一和第二波发射装置和接收装置的还另一个示例实施例;
图10示出其中第一和第二发射波束的一个是窄波束的实施例;和
图11示出其中接收器***是摄像机的实施例。
具体实施方式
在图1a中示出由线109和110定边界的第一脉冲波波束101和由线111和112定边界的第二脉冲波波束102,其中第一和第二波束两者都入射到物体103上。在这个示例中,两个波束在为纸平面的方位面中都是宽的,但是在这个示例中优选在垂直于方位面的竖直方向上是窄的。图10中示出其中在方位方向上一个波束是窄的而另一个波束是宽的实施例。继续图1,104表示沿着具有垂直于脉冲波前的单位法线n1(r)的第一波束传播的第一脉冲波场,其中r是结合的方位方向和竖直方向的三维空间中的坐标向量。第二脉冲波场105沿着第二波束传播,其具有垂直于脉冲波前的单位法线n2(r)。第一和第二脉冲波具有第一和第二中心频率ω1=2πf1和ω2=2πf2,并且例如其能够用波场表示。
脉冲波104:
脉冲波105:
第一和第二脉冲波104和105在被示为106的区域R(r,t)中重叠。重叠区沿着示为107的图像或测量曲线Γ(r)随时间传播。通过变换第一和第二脉冲波之间的发射时间的差,测量或像曲线Γ(r)能够被如箭头108所示的横向移动。根据本发明的优选实施例,Γ(r)的这个移动允许为成像扫描物体。用第一(101)和第二(102)发射波束的固定选择,Γ(r)的尾端界线113和114由波束界线109-112限制。
在这个示例中,重叠区以有效传播速度沿着Γ(r)传播(移动),其中有效传播速度为
其中c是两个波束中的波的传播速度,并且是两个波束法线之间的角。参考图1b导出这个表达式,其中104和105示出在被示为106的重叠区R(r,t)的中心在时间t的第一和第二脉波前。在时间Δt之后,波前垂直于波前移动距离cΔt到分别用于第一和第二脉冲波的波前124和125,即沿着法线n1和n2移动。现在重叠区的中心已经传播到126,并且传播线127的长度是其中是波前124和125的法线之间的角。在120所示的细节示出 其证明是等式(2)。
当该非线性材料参数能够近似于波场中的二阶时,例如能够为超声波气压波计算波动方程,其包括非线性的正向传播和散射现象,该方程为:
其中p(r,t)是脉冲波场变量,其用于超声波压力波是声压。c0(r)是用于低场振幅的线性波传播速度,βp(r)是非线性传播参数,hp(r,t)是表示加热的波能吸收的卷积核,σ1(r)和γ(r)是线性散射参数,和σn(r)是非线性散射参数。左侧传播参数在>近似波长的范围上随r改变,而右侧散射参数在<近似波长的范围上随r改变。能够计算用于电磁波的类似的等式,其表示用于EM波的类似的非线性传播和散射现象。
线性材料参数中的空间改变产生沿着波的整体传播路径的独立于彼此的波p1(r,t)和p2(r,t)两者的线性散射,并且产生频率在入射波的波段内的散射波。对于相对于波中得到的场参数的振幅材料参数中非线性足够高的材料,非线性影响波的传播和散射两者。
等式(3)中的不同项对波的传播和散射具有不同影响:线性传播项(1)引导入射波的正向空间传播,而不增加新的频率分量。线性散射源项(4)产生入射波的局部三色,而不增加新的频率分量,也就是说线性散射波具有与入射波相同的频率分量。然而,对于相对于波场振幅在材料参数中非线性足够高的材料,非线性波的影响传播和局部散射两者。非线性材料参数的缓慢的改变(在>~波长的范围上接近于常数)将提供入射波的非线性正向传播失真,其在振幅中通过项(2)随传播距离积聚/增加。非线性材料参数的快速改变(在<~波长的范围上)通过项(5)产生入射波的局部非线性散射。
非线性传播(2)和散射(5)现象是在两者都与p(r,t)2成比例的这个近似中(二阶)。对于位两个分量的和p=p1+p2的波,非线性传播和散射两者都由下列给出:
在时间域中两个函数的放大在时间频率域中产生函数时间的傅里叶变换的卷积(即时间频率谱)。这个卷积引入功能的乘积中的频率分量,其为放大的因子的频率分量的和和差。我们在散射信号中定义两个组非线性失真项:
A组来源于入射波中的正向累积非线性传播失真分量的线性散射,即为等式(3)中项(2)和项(4)的组合,其中线性散射即等式(3)的项(4),和
B组直接来源于用线性传播的入射波中原始的频率分量的局部的非线性散射,即项(1)和项(5)直接的相互作用,其中局部的非线性散射即为项(5)。
也存在C组,其辱来自入射波中正向累积非线性传播失真分量的项(5),及等式(3)中的项(2)和项(5)之间的相互作用,但是典型的非线性材料参数过低以至于该组是可忽略的。
对于散射信号的A组和B组两者,对于自失真项,入射频谱之间的卷积引入谐波波段和子谐波波段,谐波波段中和频率给出大约是入射频率的整数倍的波段,并且子谐波波段中差频率给出低于由入射脉冲的包络线确定的入射频率以下的波段。对于非线性自失真,在入射波束的整个传播路径上积聚非线性传播失真,并且因此A组分量通常比局部散射中来源的B组分量强很多,除非在像用超声造影剂微气泡的共振散射或光学荧光和双光子散射的一些情况中,其中自失真散射(B组)是同等高的。非线性自失真的A组和B组分量也具有相同波形特征,以便其通常难以提取局部散射的非线性自失真分量(B组),除非在共振散射的一些情况中。
与自失真项相反,在沿着测量或像曲线Γ(r)的两个入射脉冲波之间的重叠区R(r,t)中只得到非线性相互作用项。因此能够使得提供入射波的非线性相互作用正向传播失真的正向传播距离比用于非线性自失真项的短很多,限制A组分量的振幅。通常想将重叠区R(r,t)设计得足够小,以便B组分量的振幅比A组分量的振幅高,这允许测量或成像局部非线性相互作用散射(B组)。
如关于图2a讨论的,这个现象用于本发明的第一实施例,从而为测量或像曲线Γ(r)获得散射波的非线性相互作用分量。在这个图201和202中示出入射脉冲波104和105的时间频率谱,其中f1≈f2/3用于插图的目的。相互作用项中的两个波场的乘积等于频率域内的卷积,其以入射波段的和和差频率产生非直线型相互作用项(等式(4))中的频率波段。然后,非线性相互作用项~2p1(r,t)p2(r,t)的时间频率谱是如203和204所示的201和202的卷积。A组和B组分量提供散射信号的类似的波形特征和从而A组分量的数值为局部非线性相互作用散射(B组)的灵敏度提供下限。然而,重叠区R(r,t)的尺寸足够小,存在许多情况,其中A组分量与局部非线性相互作用散射分量(B组)相比是弱的。
沿着Γ(r)接收或测量来自R(r,t)的非线性相互作用散射信号提供沿着测量或像曲线Γ(r)的基本的测量或图像信号。例如能够用时间灵敏探测器***获得沿着提取具有空间分辨率的散射波,由于沿着Γ(r)从R(r,t)的散射发生具有从图1a中Γ(r)的顶点113到底点114的累积时间滞后。例如图7-11中所述,也能够用聚焦检测***获得接收或测量空间分辨率沿着Γ(r)的非线性相互作用散射信号,该***有选择地观察沿着Γ(r)的聚焦区域,例如摄像机***(例如,CCD摄像机)或动态聚焦的接收器元件阵列(例如使用波束形成技术),或为有利地改进空间和对比度分辨度其能够结合时间滞后和聚焦探测/拾取***。
当入射频率波段201和202被设置成以使得非线性相互作用散射波段203和204不存在重叠时,如图2a所示,例如通过用图2a所示的带通滤波器205和206带通过滤时间频率域中的接收信号,能够线性地抑制散射波段201和202提取非线性相互作用散射分量。为了增加非线性相互作用散射的灵敏度,能够形式合并以图像信号的波段203和204两者的输出。
例如通过横向改变第一和第二脉冲波的发射之间的定时,或改变第一和第二波束的至少一个的方向,或其组合能够横向扫描测量或像曲线Γ(r),并且通过这个扫描,能够为一组多次测量或像曲线Γk(r),k=1...K获得图像数据,这些曲线在物体的二维或三维空间中分布,从而形成物体的二维或三维图像。通过当曲线之间的距离关于***的基本空间分辨率足够短时插值,能够获得这些测量或像曲线Γk(r)之间的图像数据,其中该基本空间分辨率由R(r,t)结合聚焦接收器***的空间分辨率的细节给出。
能够通过使用在第一和第二脉冲波的发射之间具有相同时间关系并且第一和第二波束的方向相同的至少两个发射事件测量散射体的移动或位移,以便来自第二发射事件的重叠区沿着与用于第一发射事件相同的测量或像曲线Γ(r)传播。然后随着来自两个发射事件的接收信号之间的相位和包络线的改变,会得到两个发射事件之间散射体的移动。相位改变通常被称作多普勒效应,并且引起对移动的最灵敏的检测。其只示出沿着两个入射波的和相位的坡度的移动或位移的分量,其中在图1a、1b的示例中其沿着n1+n2。随测量或像曲线的横向扫描获得的图像中的斑点的改变,并且分析区域(即沿着每个Γ(r)的至少两个发射事件)的两个图像之间的斑点图样的改变会为3D成像给出完全的位移矢量并且为2D成像给出平面中或表面上的位移矢量的分量。
在许多情形中,通过正向传播失真或局部非线性散射或其两者,能够从入射波段201或202的谐波自失真分量得到具有非线性散射波段203和204的干涉波段,减少对非线性相互作用散射的灵敏度。改进这个情况的一个方式使发射两个脉冲组,为第二组边和p1和P2的一个的极性。然后用于第二发射的散射非线性相互作用项~的极性会改变极性,而用于A组和B组两者散射的偶数的(第二、第四...)谐波自失真分量~p1 2和p2 2不会改变极性。由此,在该方法中从这两个发射时间中减去接收信号通常被称作脉冲倒置,然后会增强入射波段的偶数谐波分量上的非线性相互作用散射项。能够通过在时间域中过滤抑制入射波段201和202的线性散射分量。
在图2a中的示例中,主要是入射波段201的谐波分量干涉非线性相互作用光谱,并且在第二发射事件中能够有利地改变P2的信号,其中减法接收信号会抵消201及其全部谐波,而检索非线性相互作用光谱203和204。因此,入射光谱202能够通过过滤抑制。在这个示例中,也能够通过在第二发射事件中只发射第一脉冲p1并且从来自包括P2的第一发射事件的接收信号减去来自这个第二发射事件的接收信号,抑制来自入射光谱201及其谐波分量的干涉,其中第一和第二发射时间都用205(和选择性地206)带通过滤。忽略第二发射事件中的第二发射波段减少物体对波的暴露,并且减少物体应力或破坏的加热或其他潜在源。本领域的任何技术人员能够研发这些方案的改变。
图2b示出图2a中光谱的变形例,其中围绕f1的发射光谱221和围绕f2的222现在这样闭合以便在很好地低于较低发射波段221的频率的围绕f1-f2的非线性相互作用卷积光谱,而在接近较高的发射波段222频率两倍的频率(即接近于222的第二谐波分量)得到围绕f1+f2的非线性相互作用卷积光谱224。通过用过滤器225带通过滤能够分离低频相互作用散布范围223,并且用过滤器226能够分离高频相互作用散布范围224。也能够使用脉冲逆转技术,其中例如高频发射光谱222的极性在第二发射事件中改变。因此围绕f2的带阻过滤能够用于去除接收信号中的222的光谱。低非线性相互作用散射波段223的减少频率提供了接收时这个波段的低吸收率并且增加了对非线性相互作用散射的灵敏度。本领域的任何技术人员能够研发这些方案的改变。
为了进一步改进用于物体中大深度的非线性相互作用项~2p1(r,t)p2(r,t)的灵敏度,其对为低吸收选择尽可能低的低频(LF)有用,但是为了得到充分对准波束选择足够高频率。选择高的高频(HF)f2,从而得到足够的空间分辨率用于到物体中的给定成像深度,例如频率比f1:f2~1:5-1:30。在尤其优选的实施例中该比值大约1:10。为了超声波成像,例如在一个应用中能够选择f1~10MHz,从而用f1~1MHz成像低至40mm,或在另一个应用中选择f2~3.5MHz从而用f1-0.3MHz成像低至150mm,也就是说频率比大约1:10。为EM波的散射能够得到类似示例。
在图2c中,对于这种情况示出用于p1的示例入射LF光谱210和用于p2的示例入射HF光谱211,其卷积和和差频率光谱为来自等式(4)中的p1和p2的非线性相互作用乘积的212和213。212和213的中心频率分别为f2+f1和f2–f1。在这个示例中和和差频率光谱212和213和入射HF光谱211彼此重叠。在频率光谱210和211之间的卷积中,在LF光谱210的波段宽度B1上得到HF光谱211的邻近频率分量的平均。这个平均化也产生光谱212和213的边缘偏斜。当HF脉冲比LF脉冲短得多时,HF光谱211中的频率分辨率比210的波段宽度B1宽,并且这个平均对和和差光谱212和213的影响可忽略,并且在等式(5-9)中的分析中允许连续LF振荡的假设(即,无限长LF脉冲)。
图2c示出这样的示例,其中在和和差光谱212和213和原始HF光谱211之间存在相当大的重叠。这给出信号处理的机会,其抑制接收信号中A组散射分量(上面定义)。由此这个抑制改进了对B组散射分量的灵敏度,并且下面进一步讨论。然而,由于重叠区的尺寸小,A组分量通常具有非常低的振幅,以便在许多情况中得到对B组散射分量足够的灵敏度,而不需要增加对A组分量的抑制。在这个情形中能够通过使用组合LF和HF脉冲发射的两个发射事件通过脉冲逆转技术从接收信号检索B组分量,改变用于第二发射事件的LF脉冲p1的极性。然后,散射非线性相互作用项(212和213)~p1p2的极性会改变用于第二发射事件的极性,而线性散射HF信号(211)不改变。然后,减去来自两个发射事件的接收信号会强烈地已知线性散射HF信号(211)并且从R(r,t)提取非线性散射HF信号(212和213)。在这个过程中也会抑制LF脉冲的偶数谐波分量。
通过在时间频率域中过滤能够在接收信号中抑制来自LF信号的频率分量,包括基本和谐波分量两者,并且从两个发射事件减去过滤的接收信号会高度抑制来自高频波段211的线性散射频率分量,并为进一步处理提取非线性相互作用散射频率波段212和213,从而沿着测量或像曲线Γ(r)形成图像信号。也能够通过使用其中只发射第一脉冲波(即LF)的第一发射事件抑制来自第一频率范围(即,LF波段)的散射信号及其谐波分量,并且从来自其中发射第一和第二脉冲波两者的第二发射事件的接收信号减去来自这个第一发射事件的接收信号。然后需要为局部非线性相互作用散射信号的脉冲逆转提取分开用于LF脉冲p1的两个极性的每个的第一和第二发射时间,这意思是全部4个发射事件提取非线性相互作用散射信号。
也能够通过发射其中p1是零而p2非零的第一事件,随后发射具有相同p2和非零p1的第二发射事件,并从用于第二发射事件的接收信号减去用于第一发射事件的接收信号,提取局部非线性相互作用。甚至能够改变两个事件之间的P2的发射振幅,来自两个事件的一个具有接收信号的类似定标,从而在减去之后最小化来自p2的线性散射信号的任何剩余。也能够为手动和自动调节两者调节用于一个接收信号的定标变量,从而最小化得到信号中的保持分量210和211,从而补偿两个发射事件之间的发射和接收通道。这也能够为脉冲逆转发射完成。当该物体移动时,能够使用多于两个发射事件并按照本领域技术人员所已知的方法估算和补偿接收信号中的多普勒位移。根据上述讨论,本领域技术人员能够由此应用许多方案从接收信号提取局部非线性相互作用散射信号分量B组。
当散射体在发射事件之间移动时,原始光谱会在上述为图2a-c讨论的两个发射事件的信号的组合中渗漏,从而抵销不需要的分量。为了在如上所述的发射脉冲波的变量(例如极性和振幅)之间区分和发射事件之间的散射***移,必须使用增加数目的发射事件也估算散射***移,其中重叠区沿着相同发射测量或像曲线传播。根据已知方法,这允许散射***移的结合估算和不需要分量的抑制。
如上面所呈现的,当高频(HF第二脉冲波)的脉冲宽度TP2比低频(LF-第一脉冲波)时间周期T1短时,尤其是比T1/2短时,能够进一步抑制来自正向运送传播失真的线性散射的分量。图5a中示出这个情形的示例,其中HF脉冲宽度TP2<<T1/2=l/2f1。对于依靠LF振荡503的顶点的HF脉冲504,等式(3)的非线性传播项(2)产生由HF脉冲观察的传播速度的改变,HF脉冲产生用于HF脉冲的额外的非线性传播延迟。为了研究这个,应当注意:|2βp(r)p1(r,t)|=|x|<<1对于|x|<<1其允许近似值1-x≈l/(l+x),并且等式(3)左侧的传播项能够近似用于HF脉冲,其包括与LF脉冲的相互作用,如下
在这个传播算子中的时间导数的前面的分子是平方传播速度,并因此看出LF脉冲p1使传播速度变形用于共同传播HF脉冲p2,如下:
其中,p1是沿着共同传播HF脉冲的实际LF场变量。因此降至HF脉冲的深度r的传播时间是
其中p1(s)是与深度有关的沿着共同传播HF脉冲的平均LF场,积分变量s其中LAF场采取这个平均值的位置。在许多实际情况中,这接近传播HF脉冲的重心,并且t0(r)是沿着HF脉冲的传播降至范围r的零LF脉冲得出的传播滞后,并且τ(r)是降至深度r的非线性传播延迟。对于βpp1>0,LF脉冲产生HF脉冲(即负的非线性传播延迟)前进的传播速度的增加,而对于βpp1<0,LF脉冲产生HF脉冲的真实延迟(即正的非线性传播延迟)的传播速度的降低。
等式(3)中的非线性散射项(5)能够为非线性相互作用项是近似的,在与下列相同范围内
由此,可以看出在这个情形中(TP2<<T1/2)LF脉冲对B组分量上的影响是散射波的振幅调制,其与2σnp1成比例。由此在这个情形中在A组和B组散射分量的波形特征之间存在大的差异,其能够被用于通过在接收信号上信号处理进一步抑制接收信号中的A组分量。
例如这种处理,假设LF脉冲的振幅在第一散射中下降非常多,以至于在第一散射之后非线性相互作用可以忽略。然后A组和B组散射分量的相同的非线性传播延迟,并且能够使用下来模型用于来自具有反向极性的LF脉冲的两个发射事件的非线性相互作用接收信号。
其中t是从深度r的散射信号的到达时间,并且a(t)源于A类散射并且b(t)源于B类散射。在文献中存在许多估算信号之间延迟的方法,并且能够通过估算非线性传播延迟为y+和y-之间延迟的一半获得B组信号的估算并且校正估算的延迟并且将信号组合为
根据已知的方法本领域技术人员能够研发更精密的估算方案。
LF场变量围绕其平均值沿着HF脉冲的改变产生传播速度(等式(6)),其沿着HF脉冲改变并产生随传播距离积累的HF脉冲的正向传播失真。对于位于LF脉冲503的顶点的HF脉冲504,LF振荡的曲率产生失真,其中LF脉冲的前部和尾部得到相同的传播速度,而HF脉冲的中心在共同传播HF脉冲的中心为βpp1>0得到较高的传播速度并且为βpp1<0得到较低的传播速度。对于围绕LF脉冲502的零交点传播的HF脉冲501,沿着HF脉冲的平均LF场变量是零并且HF脉冲的非线性传播延迟是零。然而,在这个情形中,沿着HF脉冲存在LF脉冲振荡的强坡度,其产生HF脉冲的长度改变(脉冲宽度根据βpp1的信号的压缩或膨胀)。
沿着HF脉冲的LF场变量的改变也会产生局部非线性散射中的失真。因此A组和B组之间的差异能够被总结为
A组,正向非线性相互作用传播失真:
AI:沿着共同传播HF脉冲的平均LF压力产生HF脉冲的非线性传播延迟,所述延迟在相互作用区域中随着HF脉冲的传播距离积累。
AII:沿着共同传播HF脉冲的LF压力围绕其平均值的改变产生HF脉冲的非线性传播脉冲波形失真,所述失真在相互作用区域中随着HF脉冲的传播距离积累。
B组,局部非线性相互作用散射失真:
BI:沿着共同传播HF脉冲的平均LF压力产生与沿着HF脉冲的LF场变量的局部平均成比例的局部散射HF信号的非线性振幅调制。
BII:沿着共同传播HF脉冲的LF压力围绕其平均值的改变产生由沿着HF脉冲的LF场变量的局部改变给出的局部散射高频信号的脉冲波形失真。
在BII下的脉冲波形失真与AII下的不同,由于在AII下网络失真由沿着HF脉冲的传播路径的局部失真的累积给出,其中局部失真沿着传播路径改变,并且如下面讨论的甚至正负摆动,而在BII下该失真由LF脉冲场变量沿着在局部散射体的HF脉冲的局部改变给出。
该脉冲波形失真能够由操作人员在A组和B组散射分量上建模为
其中y+和y-是来自发射LF脉冲极性相反的两个发射事件的接收信号。Hq,q=±和b±是表示脉冲波形失真的算子。
为了缓和脉冲波形失真,能够在等式(11b)中由线性过滤器模型这些算子,脉冲响应hq,q=±和b±。例如能够从材料中的估算βp(r)p1(r)和材料中共同传播脉冲的非线性波传播的计算机模拟获得脉冲波形失真的估算。例如,能够从等式(7)获得βp(r)p1(r)的估算,其为估算的非线性传播延迟的空间梯度。
当得到失真算计的估算并且其具有逆转时,能够获得B组散射信号的估算,其具有对A组分量的强的抑制,如下
为了在图2c所示的情况中的A组和B组散射的波形特征的差异的更详细的分析,即Tp2<T1,假设LF和HF脉冲具有下列简单的平面波形
LF脉冲:p1(r,t)=a1(r,t-t1-n 1 r/c)cosω1(t-t1-n 1 r/c)
HF脉冲:p2(r,t)=p2(t-n 2 r/c) (13)
跨时间脉冲宽度TP2<<Tp1的HF脉冲,LF脉冲宽度,近似a1(r,t-t1-n 1 r/c)≈a1,即为常数,其中同意省略由相对于上面图2c讨论的入射LF光谱210的频率平均。然后,非线性相互作用散射源项源于
为第二脉冲(HF)引入迟至时间τ=t-n 2 r/c,其给出
并且估算τ上的傅里叶变换,如下
其中用ω=2πf定义相位θ(t1,r)=ω1(t1+(n 1-n 2)r/c)。注意P2(ω-ω1)是图2c中的212并且P2(ω+ω1)是213。进一步地注意HF脉冲的脉冲宽度Tp2与波段宽度B2有关,为Tp2=1/B2。然后具有等价式因为HF和LF脉冲的传播方向不同,所以正向传播延迟的累积和沿着HF脉冲传播方向的脉冲失真会具有最大延迟和失真围绕零振荡,并且具有最大和零的非线性相互作用散射信号的振幅的振荡取决于θ(t1,r)的值。右侧的第一相,即光谱的和项,对于θ(t1,r)=mπ,m=0,±1,...具有最大值,其中cosθ(t1,r)=±1。右侧的第二项,即光谱的差项对于θ的这些值是零,其中sinθ(t1,r)=0。对于θ(t1,r)=(2m+l)π/2,m=0,±1,...和项为零,其中cosθ(t1,r)=0,因为对于θ的这些值sinθ(t1,r)=±1用于差项最大。光谱P2(ω-ω1)和P2(ω+ω1)之间强烈重叠,差项被高度减少并且比和项提供更少的传播失真和散射。
对于B2<f1,光谱P2(ω+ω1)和P2(ω-ω1)具有与P2(ω)的最小重叠,并且和项P2(ω-ω1)+P2(ω+ω1)和差项P2(ω-ω1)-P2(ω+ω1)近似产生相同的散射量,除非两项的相对数值取决于θ(t1,r)。然而,因为HF脉冲传播方向n 2≠n 1-n 2,cosθ(t1,r)和sinθ(t1,r)项会在R(r,t)内随着r摆动,并且得到空间摆动并且近似等于从和和差项到R(r,t)内的非线性相互作用术语的贡献。在这个情形中,A组和B组项之间的波形特征的差异低。因此用信号处理难以抑制A组项,并且由此A组信号限制对B组信号的灵敏度,即局部非线性相互作用散射信号。因此对于B2<f1必须依靠小的重叠区R(r,t),从而能够测量局部非线性相互作用散射(B组)。
当B2增加超过f1时,光谱P2(ω+ω1)和P2(ω-ω1)可以与图2c中的P2(ω)重叠,并且进一步当B2增加超过2f1时,即TP2<T1/2,光谱P2(ω+ω1)和P2(ω-ω1)开始重叠。为了用强重叠谱进一步分析,重写等式(16)有用,如下
右侧上的第一项产生具有线性散射的正向传播失真的AI组(非线性传播滞延迟),并且BI组振幅调制具有局部非线性相互作用散射的a1。在正向传播失真(AII组)和局部非线性散射(BII组)两者中,第二项产生脉冲波形失真,其与沿着HF脉冲的LF脉冲的曲率有关。在正向传播失真(AII组)和局部非线性散射(BII组)两者中,第三项产生脉冲波形失真,其与沿着HF脉冲的LF脉冲的坡度有关。在光谱P2(ω-ω1)和P2(ω+ω1)之间具有强重叠,即第二和第三项(脉冲波形AII和BII)与第一项(用于AI组的传播延迟和用于BI组的散射振幅调制)相比被高度减少。
由此对于这个情形(图2c的TP2<T1/2)在B组(来自R(r,t)的局部非线性相互作用散射信号)和A组(从R(r,t)内的聚集的非线性相互作用传播失真分量的线性散射上升的信号)之间的存在强烈差异。这个波形特征的差异用于由等式(10、12)所述的信号处理,从而与如图2a和2b中举例说明的发射波相比,抑制接收信号中的A组分量,并增加对来自R(r,t)的局部非线性相互作用散射信号(B组)的灵敏度,其中Tp2>T1。
如上面关于图2c所讨论的,为了提取非线性相互作用散射信号,能够使用组合LF和HF脉冲发射的两个发射事件,为第二发射事件改变LF脉冲p1的极性。然而,在减去接收信号之前,用延迟和脉冲波形修正的一个或两者修正接收信号的一个或两者,例如在修正过滤器中,在减去信号之前。然后,校正过滤器产生下列的至少一个,i)延迟修正,其补偿非线性相互作用传播延迟(上面AI)的影响,和ii)脉冲波形失真修正,补偿线性散射接收信号中的非线性相互作用脉冲波形传播失真(上面AII)。
当HF脉冲波段宽度比LF大很多时(例如B2~5f1),重叠区分为子区域。为了更好地理解这个现象,参考图3,其中图3中的301示出灰度图中在固定时间在空间中随着其在n1的方向上传播的示例LF脉冲波场(104)的脉冲振荡。302示出LF脉冲沿着线303的空间振荡改变。304在灰度图中示出在固定时间在空间中随着在n2的方向上传播的示例HF脉冲波场(105)脉冲振荡。305示出HF脉冲沿着线306的空间改变的放大图例。图3b示出非线性相互作用源项2p1(r,f)p2(r,t),即图3a中的脉冲的乘积。可以看出为LF振荡的高值(两者正的和负的)得到源项的高值,并且在围绕LF振荡的零的源项中具有低值。这个产生非线性相互作用源项子区域308,其在中间具有非线性相互作用源项的低值子区域309。
非线性相互作用子区域308在n1+n2的方向上传播。应当注意,等式(8、9)中的相位θ(t1,r)=ω1(t1+(n 1-n 2)r/c)取决于沿着(n1-n2)的r的分量,并且由于可以看出θ(t1,r)在非线性相互作用子区域的传播方向上是常数。t1的改变(即,LF和HF脉冲之间的发射定时)引起非线性相互作用源项子区域的横向移动。当θ(t1,r)=mπ,m=0,±1,...时,cosθ(t1,r)=±1并且sinθ(t1,r)=0,并且HF脉冲取决于LF波的顶点或谷点。然后,得到等式(17)的延迟项和和项的最大值,即等式(17)右侧的第一和第二项,而等式(17)的右侧第三项(差项)是零。当θ(t1,r)=(2m+l)π/2,m=0,±1时,即cosθ(t1,r)=0并且sinθ(t1,r)=±1,HF脉冲沿着LF波的零并且位于低值子区域309中,其中差项是最大值并且延迟和和项是零。然而,如上面所讨论的,由于HF脉冲足够短,差项是可忽略的,并且非线性相互作用散射(B组)的灵敏度低。
由于HF带宽和LF之间具有这样大的差异(例如B2~5f1),从而能够为单个发射事件同时分析多次测量或像曲线,其中图4举例将三个测量或像曲线Γk-1(r)、Γk(r)、Γk+1(r)示为410、404、411,用于为三个不同的时间点所示的三个非线性相互作用子区域ta、tb、tc示为401、402、和403。由图中弯曲的箭头所示,非线性相互作用源项在n1+n2的方向上以等式(2)的速度cR向下传播。可以看出存在θ(t1,r)=mπ,m=0,±1,...的高非线性相互作用的区域(A区域),其中对成像非线性相互作用参数中空间改变的灵敏度高,用于足够短的HF脉冲。示出三个区域405、406、和407,但是如点412所示^-区域能够根据第一和第二发射波束101和102的LF脉冲的长度及其他细节进一步横向延伸。图4也举例示出用于θ(t1,r)=(2m+l)π/2,m=0,±1,...的低非线性相互作用的两个B区域408和409,其中对成像非线性相互作用参数中的空间改变的灵敏度低,用于足够短的HF脉冲。在竖直方向(垂直于纸平面)上重叠子区域(401、402、403)的厚度由在竖直方向上最薄的波束101和102的竖直厚度确定。
在典型的测量情况中能具有f2=10f1。图4中A类型的子区域的宽度是其中T1是LF震荡的时间周期,并且λ1和λ2分别是LF和HF脉冲波的宽度。例如用于图7-11中所述的***,用于HF的聚集接收器***的分辨率是Δr=2λ2FN,其中FN=F/D是接收器的光圈值,其作为焦距F和光圈宽度(例如,直径)D之间的比值给出。可以看出用的聚焦接收器***分析子区域的宽度,其对于给出FN<2.1。用许多EM和EL接收器***可以实现这些低光圈数。因此能够用平行接收器***分析几个A区域(例如405、406、407),并且用波束宽度例如由在测量或像曲线Γk-1(r)、Γk(r)、Γk+1(r)的每侧上的点垂直线显示的动态聚焦平行接收波束从这些区域接收信号。这些平行接收波束为几个曲线(例如,410、404、411)同时提供测量或图像信号,其中用相同的发射事件和根据已知方法在许多阵列元件上的接收信号上并行处理获得图像数据,并且参考图6的仪器方框图进一步地讨论。
为了增加对区域B中的非线性相互作用参数的成像灵敏度,能够发射第二组第一和第二脉冲波,其中如图5a所示改变第一和第二脉冲之间的定时关系量Δt1=π/2ω1,其中501示出用于位于低频脉冲502的零的对非线性相互作用参数成像灵敏度低的区域B的HF脉冲。在这个图中,HF脉冲501和LF脉冲502示出用于图4中低灵敏度子区域(区域B)的定时关系。在第二发射组中,改变LF和HF脉冲之间的发射定时的差异量Δt1=π/2ω1,以便LF脉冲到达相对位置503,其中现在HF脉冲504在对区域B中的非线性相互作用参数具有最大成像灵敏度的峰值LF场变量上。这与图4中的第一发射组的灵敏度子区域401、402、和403一起生成用于在图5b的区域B中的时间点ta、tb、和tc的灵敏度子区域510、511、和512。点515指示能够横向具有更多A区域和B区域,这取决于如上面参考图4讨论的波束的LF脉冲的长度及其他细节。由于如上面示例中所述的接收器***的足够高的焦点分辨率(其中FN<2.5),t1的改变提供了测量或新的像曲线,其中能够使用用接收信号的并行处理的许多平行接收波束沿着几个测量或像曲线与相同发射事件并行提供图像信号,其中图5b举例显示两个测量或像曲线Γm-1(r)、Γm(r)为用于5-区域的513和514。用于具有相同发射事件的A-和B-区域两者的平行测量或像曲线的总数取决于发射脉冲104和105的细节形状。
入射频率的大的差异也对提高共振非线性散射体的成像有用,其中选择第二发射频率达达高于共振频率,例如从而获得高的空间分辨率。然后,选择低于或等于共振频率的第一发射频率,并且使用参考图2c所述的方法从非线性相互作用共振散射体提取信号。例如,这对在大大高于共振频率的频率下高分辨率声成像造影剂微气泡有用,以及对高分辨率电磁波成像荧光分子有用。
使用发射事件之间的增加重复率,其中在第一和第二脉冲波之间的相对发射定时t改变,能够获得为许多重叠区Rt(r,t),i=2,...的进一步增加平行测量或成像,这些区域沿着物体中不同的测量或像曲线Γi(r),i=2,...同时传播。在这个情形中聚焦接收***适合于平行观察来自在物体中同时传播的不同重叠区的散射信号。该方法也为入射频率f1和f2提供平行测量,这些频率如此靠近以至于B2<f1,以便重叠区不分解为如图3-5所述的子区域,其中B2﹥﹥f1。该平行接收成像增加潜在帧速率,其打开平均接收信号或图像信号或图像信号的可能性,从而改进信号噪音比并增加用于成像从材料的局部非线性相互作用散射的灵敏度。
在上述讨论中,假设具有波场的非线性材料参数的第二顺序改变。这提供了等式(4)的简单的非线性相互作用项2p1(r,t)p2(r,t),其进一步提供图2a-c的卷积频谱。由于材料参数的更复杂非线性,能够具有有非线性波传播和非线性散射引入的频率分量。第二顺序非线性用于说明的目的,并且不是意图为本发明提供限制。本领域技术人员在本发明的保护范围中能够为每个特殊情况引入更复杂的非线性材料参数的细节,例如包含非线性共振散射体的波在材料中的传播和散射。
图6示出根据本发明的示例仪器的方框图,其中可以根据所述的方法操作仪器。该仪器包含用于分别将脉冲波束101和102发射到物体103中的装置601和602。至少一个波束在方位方向上(纸平面中)是宽的,并且在竖直方向上(垂直于纸平面)是窄的,其中这个图示出波束101和102两者在方位面上都是宽的示例。发射装置是由发射器603和604驱动的,其由控制单元605经由控制总线606在发射定时和发射振幅和发送频率和发射波束方向的至少一个上控制。根据已知方法,控制单元经由键盘或其他的用户界面从仪器操作员接收输入,从而为实际使用安装仪器。换句话说,控制单元605能够从软件或从反馈机构接收输入。接收器装置607从物体拾取散射波608,从而形成接收信号传送给接收信号预处理单元609,在这里这个单元的输出呗传送到图像处理单元610用于根据本发明进一步处理,以便单元607、609、和610中的信号调节和处理的最后结果是沿着至少一个测量或像曲线Γ(r)(107)提取局部非线性相互作用散射信号(B组),例如用参考图2-5所述的方法,并且形成测量或图像信号,其被传送到显示装置611用于显示。显示单元举例说明沿着一个测量或像曲线Γ(r)的非线性相互作用信号的图像,其为强度或颜色调制的显示线612。沿着测量或像曲线Γ(r)(107)的非线性散射信号的选择能够在接收器装置607、预处理单元609或图像处理单元610、或在结合的全部单元607、609和610中部分地或完全完成。通过变换第一和第二发射脉冲之间的发射时间关系和发射脉冲波束101和102的相对方向的至少一个,能够通过例如跨二维表面区域或在物体体积区域中的三维中横向扫描测量或像曲线107,从而提供物体的区域的二维(2D)或三维(3D)图像。
图6中的单元表示功能单元,并且在各种实施例中能够被疑相同物理单元(硬件)组合,这取决于可利用的优选技术和性能要求。例如,控制单元605和图像处理单元610能够被结合为单个计算机***上的软件模块,并且该单元603、604、和609的电子装置能够被制成到计算机总线606的插件。
即使也用接收器装置607,发射装置601和602能够类似地被组合为一个物理单元。在图7中得到这种情况的示例用于超声波成像,其中示出在物体701的表面702上放置的线性超声阵列。超声阵列根据已知的原理操作,并且能够沿着波束101和102同时发射第一和第二脉冲超声波104和105。对于大大分离的第一和第二频率,例如能够使用美国专利6,645,150或美国专利7,727,156中所述的阵列解决方法。能够用于接收波束703相同的阵列接收散射信号,其包含重叠区R(r,t)(106)。在这个示例中接收波束703包含整体测量或像曲线Γ(r)(107),其中通过接收信号的时间滞后获得沿着Γ(r)的空间分辨率。然后,应该注意散射信号的深度以等式(2)中重叠区R(r,t)的速度cR传播,而背散射信号以速度c传播,并且接收时间和沿着Γ(r)的深度之间的转化考虑此情况。如上所述并根据已知原理,能够通过使用的大孔径、动态聚焦接收波束改进与接收波束横向的空间分辨率。也能够使用几个并排平行接收波束,其同时从图4和5的许多A和B区域接收信号,根据已知方法并行处理来自全部区域的接收信号从而加速帧速率。在这个情况下,接收信号能够包含从重叠区R(r,t)的内侧和外侧两者线性散射的分量和来自内侧R(r,t)的非线性相互作用散射信号分量,特别是当在第一和第二频率之间存在大的差异时。然后在图6的单元609和610中通过参考图2-5所述的方法提取非线性相互作用散射信号。
图7中的非线性相互作用散射图像场的深度由波束101和102的宽度界定,并由此得到三角形式。图8中示出提供非线性相互作用散射图像区域的矩形形式的实施例,其中用于第一和第二波束101和102的发射装置是物理上分离的线性超声换能器阵列801和802。用于第一和第二脉冲波(104、105)、重叠区R(r,t)(106)和测量或像曲线Γ(r)(107)的标识与图1a和7中相同。来自第一阵列801的接收波束聚集到重叠区R(r,t)上,该波束被示为803。应该注意,在这个情形中所示的接收波束不覆盖整体测量或像曲线Γ(r),因此为了接收沿着Γ(r)的全部的非线性相互作用信号,波束位置/方向必须与焦点深度一起调向。根据已知方法,在一个示例中这能够通过选择不工作的有源接收元件,和/或为适当的方向调向调节延迟和聚焦波束实现,以便由于其实时移动,该聚焦按照重叠区R(r,t)。在另一个示例中,根据已知方法能够在801的全部接收元件上同时接收散射信号,并同时在处理单元中形成聚焦图像波束,以便随着重叠区R(r,t)获得接收波束交点。
图9示出根据本发明的用于超声波成像的集成发射与接收装置的再另一个实施例,其中超声阵列被设置为环形901。当物体能够由阵列围绕时,例如胸部成像时,这个实施例具有的优点在于其允许按照已知方法断层分析图像重构。该阵列由沿着整体环形901的独立元件903组成,该环形提供大的柔韧性,从而为发射和接收孔隙选择元件。在这个示例中,选择元件组904从而形成第一发射孔隙发射第一脉冲波波束101,并且选择另一元件组905从而形成第二发射孔隙发射第二脉冲波波束102。用于第一和第二脉冲波(104、105)、重叠区R(r,t)(106)和测量或像曲线Γ(r)(107)与图1a、7、和8中相同。该示例进一步示出接收元件组906,其给出接收波束907,以便接收波束的轴的全长沿着测量或像曲线Γ(r)(107)。例如,根据已知方法,扫描曲线F(r)形成图像能够通过改变第一和第二脉冲波的第一和第二发射时间之间的关系完成,或能够通过选择分享第一(904)和第二(905)发射孔隙的不同元件以及使用用于不同元件的发射脉冲的可变延迟,改变第一和第二脉冲波的第一和第二波束方向。然后完成用于接收孔隙(905)的元件的选择和来自元件的信号的电位延迟转向,从而匹配沿着每个实际测量或像曲线Γ(r)的接收波束的扫描。
在图1a和7-9中,脉冲波前被示为相当直的线,其产生也是直线的测量或像曲线。在实际情况中,具有有限的孔隙的折射使得波前变得有些弯曲,尤其是在波束的边缘。这使得测量或像曲线Γ(r)产生一些曲率,其能够被认为是电子接收波束转向。在相同的图中,已经说明第一和第二发射波束在相同平面中。在许多情况下这能够是有利的,由于曲线Γ(r)的横向扫描能够通过改变第一和第二脉冲波的发射定时关系获得。然而,从本发明的大体原理应当清楚,第一和第二发射波束能够在彼此交叉的不同平面中,其中发射定时关系必须被调节到波束彼此交叉,这在实际情况中能够由本领域技术人员改造。
按照本发明的发射和接收装置的变形例的另一个示例在图10中被示为物体1000表面上的1001,其中至少在方位面上第一和第二发射波束的一个是窄波束1002并且第一和第二发射波束的另一个是宽波束1003。然后,如1004所示,沿着波束1003的发射脉冲必须足够长,以便为想要的图像范围的整体长度沿着1002在两个脉冲波1004和1007之间得到相互作用。测量或像曲线Γ(r)是窄波束1002的轴1006,并且接收波束轴优选地沿着窄波束1002的轴对准。然后由宽波束(1003)的长脉冲1004和窄波束(1002)的短的发射脉冲1007的重叠给出重叠区R(r,t)1005。为了改进可能的分辨率,有利的是沿着窄波束发射高频脉冲。由于第一和第二波束之间的传播方向偏斜,类似于图4和5所示的,沿着波束1002会具有对非线性相互作用散射的灵敏度高和低的区域。类似于图5a和b中所示,能够由本领域技术人员改造的,通过改变第一和第二发射脉冲之间的发射定时的差异能够沿着1002移动高灵敏度区域。为了沿着整个接收波束得到对非线性相互作用散射的完全灵敏度,然后必须以所述第一和第二发射时间之间不同的时间滞后发射多于一个发射事件。通过横向扫描窄波束横向地扫描测量或像曲线Γ(r)。
在图7-10中,例如已经使用超声波成像,但是按照本发明发射器和接收器的类似几何设置也能够与EM波一起使用。为了以在GHz和THz范围中的频率EM成像,发射装置和接收装置能够是带状天线或微波激射器/激光二极管、和这些元件的阵列。为了在红外线光频范围中EM成像,用于发射装置的简单的解决方法是激光二极管系列,或机械方向转向激光二极管。用于接收探测装置的简单的解决方法能够是光敏二极管/晶体管或聚焦摄像机***(例如CCD摄像机),其提供来自R(r,t)的散射信号沿着Γ(r)的实时成像,其中可以直接观察到通过第一和第二脉冲波的相对的发射定时和/或第一和第二波束方向的改变横向扫描Γ(r)。图11中示出这种接收装置的示例实施例,其中1101和1102示出用于第一和第二波束101和102的发射这种,其脉冲波前104和105具有重叠区R(r,t)(106)和测量或像曲线Γ(r)(107)。接收器装置是EM摄像机***1109,其观察锥体1111聚焦到来自区域1110的非线性相互作用散射信号的观察区域上,该区域1110覆盖两个发射波束101和102直接的相互作用。
为了进一步增加对非线性相互作用散射信号的灵敏度,能够根据已知方法方便地将来自用于每个单独测量或像曲线Γ(r)的许多发射事件的接收信号或图像信号平均。
因此,尽管已经示出并描述并指出如应用于本发明的优选实施例的本发明的基本新颖特征,但是应当理解,本领域技术人员可以进行所示装置的形式和细节上及其操作的各种省略和代替和改变,而不偏离本发明的精神。
也明确地意图以大体上相同的方式实行大体上相同的功能的这些元件和/或方法步骤的全部组合在本发明的保护范围中。元素和/或方法台阶哪些实行实质上一样功能在实质上同样地到实现一样结果在...的区域之内发明的范围是也明确地打算让全部的综合。此外,应该认识到,结合本发明任何公开的形式或实施例所示和/或所述的结构和/或元件和/或方法步骤可以以任何其他所公开或所述或所提出的形式或实施例合并,作为设计选择的一般主题。因此,仅由权利要求的保护范围所述的限制本发明的意图。
Claims (34)
1.一种方法,其用于沿着具有产生波散射的非线性材料参数的材料物体中的至少一个测量或像曲线Γ(r)测量或成像两个波之间的非线性相互作用散射,所述方法包含的步骤:
a)在至少一个发射事件中发射,其中i)至少一个第一发射脉冲波,其沿着具有第一发射方向的第一发射波束具有第一发射时间和第一发射中心频率,和ii)至少一个第二发射脉冲波,其沿着具有第二发射方向的第二发射波束具有第二发射时间和第二发射中心频率,其中所述第一和第二发射方向和所述第一和第二发射时间被设置成使得所述第一和第二脉冲波的脉冲空间重叠在沿着测量或像曲线Γ(r)的材料物体中随时间传播的重叠区R(r,t)中,和
b)选择i)所述第一和第二发射时间之间的时间关系,和ii)所述第一和第二发射方向的至少一个中的的一项或两者,控制所述材料物体内的所述测量或像曲线Γ(r)的所述位置,和
c)用至少一个接收器检测从所述重叠区R(r,t)散射的波,并且生成至少一个接收信号,所述接收信号包括至少非线性相互作用散射信号分量,其涉及在所述重叠区R(r,t)中的所述第一和第二脉冲波之间的所述非线性相互作用散射的波,和
d)处理所述接收信号以提取所述非线性相互作用散射信号分量,并且沿着所述测量或像曲线Γ(r)形成非线性相互作用测量或图像信号。
2.根据权利要求1所述的用于测量或成像的方法,其中所述提取非线性相互作用信号分量的过程包括所述接收信号的时间频率过滤,其中所述过滤通过所述非线性相互作用散射信号分量并且阻断其他时间频率分量。
3.根据权利要求1或2所述的用于测量或成像的方法,其中
a)所述发射步骤包含至少两个发射事件,并且其中所述第一和第二发射方向和所述第一和第二发射时间之间的所述时间滞后对于所述至少两个发射事件是相同的,并且其中
b)所述第一和第二脉冲波中的一个在所述至少两个发射事件之间的极性、振幅、相位和频率中的至少一项改变,包括所述第一和第二脉冲波的至少一个对于所述发射事件中的至少一个具有零振幅的可能性,并且其中
c)提取非线性信号分量的所述过程包括组合来自所述至少两个发射事件的所述接收信号。
4.根据权利要求3所述的用于测量或成像的方法,其中在来自所述至少两个发射事件的信号在步骤c)中组合之前,来自所述至少两个发射事件的至少一个的所述接收信号是i)延迟修正和ii)脉冲失真修正中的一个或两者。
5.根据前述任一权利要求所述的用于测量或成像非线性相互作用散射的方法,其中
a)通过选择i)所述第一和第二发射时间之间的所述时间关系,和ii)所述第一和第二发射方向的至少一个中的一项或两者在所述材料物体中选择多个测量或像曲线Γk(r),k=1...K,和
b)根据前述任一权利要求获得相应于每个具体的测量或像曲线Γk(r)的所述测量或图像信号。
6.根据前述任一权利要求所述的用于测量或成像非线性相互作用散射的方法,其中所述发射步骤包含多个发射事件,并且所述发射事件之间的时间差是足够短的,以便在所述物体内同时传播至少两个不同的重叠区R(r,t),并且所述第一和第二发射脉冲之间的相对定时和所述第一和第二发射方向中的一个和两者在所述发射事件之间改变,以便所述至少两个不同的重叠区沿着至少两个不同的测量或像曲线Γ(r)传播,并且来自所述不同的重叠区的所述散射信号被并行地接收和处理以改进测量或成像率。
7.根据权利要求3所述的用于测量或成像的方法,其中
a)所述第一和第二脉冲波中的一个的脉冲宽度比所述第一和第二波中的另一个的时间振荡周期的一半短,以便所述非线性相互作用重叠区被分为高灵敏度的子区域,其中所述第一和第二波的所述另一个的波场的大小足够高,以便并行地从多于一个测量或像曲线的一组Γk(r),k=1..K获得具有高灵敏度的散射信号,和
b)来自所述多于一个测量或像曲线的所述接收信号被并行地接收和处理,从而为所述多于一个测量或像曲线并行地获得所述图像信号。
8.根据前述任一权利要求所述的方法,其中来自用于所述相同测量或像曲线Γ(r)的多个发射事件的接收信号被平均以减少噪声并改进灵敏度。
9.根据权利要求5所述的用于测量或成像非线性相互作用散射的方法,其中通过在至少两个测量或像曲线上的图像信号之间插值获得用于所述测量或像曲线Γk(r)的图像点的图像信号。
10.根据权利要求5或9所述的用于测量或成像非线性相互作用散射的方法,其中所述测量或像曲线Γk(r)被跨表面空间分布,从而获得用于跨所述表面的非线性相互作用散射的二维图像的图像信号。
11.根据权利要求5或9所述的用于测量或成像非线性相互作用散射的方法,其中所述测量或像曲线Γk(r)被空间分布在三维区域内,从而获得用于所述区域中的非线性相互作用散射的三维图像的图像信号。
12.一种仪器,其用于沿着具有产生波散射的非线性材料参数的材料物体中的至少一个测量或像曲线Γk(r)测量或成像两个波之间的非线性相互作用散射,所述仪器包括:
a)发射装置,其被设置成在由下面构成的至少一个发射事件中发射:i)至少一个第一发射脉冲波,其沿着具有第一发射方向的第一发射波束具有第一发射时间和第一发射中心频率,和ii)至少一个第二发射脉冲波,其沿着具有第二发射方向的第二发射波束具有第二发射时间和第二发射中心频率,其中所述发射装置被设置成使得所述第一和第二波束重叠在材料物体的重叠区R(r,t)中,其中所述重叠区沿着测量或像曲线Γ(r)在材料物体中随时间传播,和
b)控制装置,其被设置成选择i)所述第一和第二发射时间之间的时间关系,和ii)所述第一和第二发射方向中的至少一个中的一项或两者,从而在空间上控制所述材料物体内的所述测量或像曲线Γ(r)的所述位置,和
c)接收装置,其被设置成至少从所述重叠区R(r,t)中所述第一和第二脉冲波之间的非线性相互作用接收材料物体的散射波分量,并且被设置成生成至少一个接收信号,所述接收信号包括来自至少所述重叠区R(r,t)的至少非线性相互作用散射信号分量,和
d)处理装置,其被设置成处理所述接收信号以从所述重叠区R(r,t)提取所述非线性相互作用散射信号分量,以沿着所述至少一个测量或像曲线Γ(r)形成非线性相互作用测量或图像信号。
13.根据权利要求12所述的仪器,其中所述处理装置包括用于所述接收信号的时间频率过滤的装置,其中所述过滤通过所述非线性相互作用散射信号分量并且阻断其他时间频率分量。
14.根据权利要求12或13所述的仪器,其中
a)所述发射装置被设置成为每个测量或像曲线Γ(r)发射至少两个发射事件,和
b)所述发射装置包括用于改变所述至少两个发射事件之间的所述第一和第二脉冲波的一个的极性、振幅、相位和频率的至少一个的装置,包括所述第一和第二脉冲波的至少一个对所述至少两个发射事件的至少一个具有零振幅的可能性,并且其中
c)所述处理装置包括用于在沿着所述测量或像曲线Γ(r)形成图像信号的过程中结合来自至少两个发射事件的所述接收信号。
15.根据权利要求14所述的仪器,其中所述处理装置包括这样的装置,其在来自至少两个发射事件的所述信号被结合之前,用i)延迟修正和ii)脉冲失真修正中的一项或两者修正来自所述至少两个发射事件的至少一个的至少一个接收信号。
16.根据权利要求15所述的仪器,其中所述处理装置包括用于估算i)延迟修正和ii)脉冲失真修正中的一项或两者的装置。
17.根据权利要求12到16的任何一项所述的仪器,其中所述用于接收的装置包含聚焦接收器***。
18.根据权利要求17所述的仪器,其中所述聚焦接收器***被设置成随时间动态聚焦,以便焦点在重叠区R(r,t)随时间移动时跟随重叠区R(r,t)。
19.根据权利要求17或18的所述的仪器,其中所述聚焦接收器***是聚焦摄像机。
20.根据权利要求17、18或19所述的仪器,其中所述聚焦接收器***包括:
多个接收器元件的阵列,其被设置成并行接收所述散射信号,和
装置,其从多组接收器元件并行生成多个接收信号,以便不同的接收信号被并行地分配给不同的测量或像曲线Γm(r),m=1...M。
21.根据权利要求20所述的仪器,其中所述控制装置包含发射多组发射事件从而为每个测量或像曲线Γ(r)接收多个接收信号的装置和用于处理所述多个接收信号的装置,所述处理通过下列至少一个执行:
i)将所述接收信号平均以增加信噪比并由此增加灵敏度,和
ii)比较所述多个接收信号之间的差异,从而确定所述散射体的移动。
22.根据权利要求12到21中的任何一项所述的仪器,其中所述控制装置包括这样的装置,其通过改变i)所述第一和第二发射时间之间的时间关系和ii)所述第一和第二发射方向中的至少一个的一项和两者,通过所述材料物体的区域中的一组不同的测量或像曲线Γk(r),k=1...K扫描所述测量或像曲线Γ(r)。
23.根据权利要求22所述的仪器,其中所述用于处理的装置包括用于在邻近测量或像曲线Γk(r)之间插值图像信号从而为所述邻近测量或像曲线之间的图像点生成图像信号的装置。
24.根据权利要求21、22或23所述的仪器,其中所述控制装置包括用于跨贯穿所述材料物体的表面扫描所述测量或像曲线从而为跨所述表面的非线性相互作用散射的二维图像提供图像信号的装置。
25.根据权利要求21、22或23所述的仪器,其中所述控制装置包括用于在所述材料物体的三维区域内扫描所述测量或像曲线从而为所述区域中的非线性相互作用散射的三维图像提供图像信号的装置。
26.根据权利要求12到25的任何一项所述的仪器,其中所述第一和第二发射波束的至少一个在至少方位方向上是宽的。
27.根据权利要求12到26的任何一项所述的仪器,其中所述第一和第二发射波束的至少一个在竖直方向上是窄的。
28.根据权利要求27所述的仪器,其中所述第一和第二发射波束的另一个至少在方位方向上是窄的,并且所述接收装置包括聚焦接收波束,该聚焦接收波束的轴沿着在所述方位方向上是窄的所述发射波束的轴,并且扫描测量或像曲线是通过扫描所述接收波束和在方位方向上是窄的所述发射波束两者获得的。
29.根据权利要求12到28的任何一项所述的仪器,其中所述发射装置和接收装置包括元件阵列,其中相同的阵列表面能够被用于发射所述第一和第二脉冲波两者并至少接收所述非线性相互作用散射信号。
30.根据权利要求12到29中的任何一项所述的仪器,其中所述发射装置包括发射所述第一脉冲波的元件阵列和发射所述第二脉冲波的另一个元件阵列。
31.根据权利要求12到30中的任何一项所述的仪器,其中所述发射装置包含元件的环形阵列。
32.根据权利要求12到31中的任何一项所述的仪器,其中所述第一和第二脉冲波的至少一个是弹性压缩波。
33.根据权利要求12到32中的任何一项所述的仪器,其中所述第一和第二脉冲波的至少一个是弹性横波。
34.根据权利要求12到33中的任何一项所述的仪器,其中所述第一和第二脉冲波的至少一个是电磁波。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US37540310P | 2010-08-20 | 2010-08-20 | |
US61/375,403 | 2010-08-20 | ||
GB1018933.0A GB2484753B (en) | 2010-08-20 | 2010-11-09 | Method for imaging of nonlinear interaction scattering |
GB1018933.0 | 2010-11-09 | ||
PCT/GB2011/051564 WO2012035312A2 (en) | 2010-08-20 | 2011-08-19 | Method for imaging of nonlinear interaction scattering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103261909A true CN103261909A (zh) | 2013-08-21 |
CN103261909B CN103261909B (zh) | 2016-01-20 |
Family
ID=43414607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180050704.4A Expired - Fee Related CN103261909B (zh) | 2010-08-20 | 2011-08-19 | 用于成像非线性相互作用散射的方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8793079B2 (zh) |
EP (1) | EP2606374A2 (zh) |
CN (1) | CN103261909B (zh) |
GB (1) | GB2484753B (zh) |
WO (1) | WO2012035312A2 (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107231811A (zh) * | 2014-12-10 | 2017-10-03 | 瑟夫技术公司 | 用于非线***互散射的成像的方法 |
CN107530046A (zh) * | 2015-01-15 | 2018-01-02 | 罗德尼·赫林 | 漫射声学共焦成像器 |
CN107797160A (zh) * | 2017-09-01 | 2018-03-13 | 上海交通大学 | 弹性波和电磁波ct勘测数据联合分析***及方法 |
CN108338808A (zh) * | 2017-01-23 | 2018-07-31 | 美国西门子医疗解决公司 | 使用相干的剪切速度成像 |
CN109124680A (zh) * | 2017-06-16 | 2019-01-04 | 美国西门子医疗解决公司 | 用于形成超声探头的超声信号的方法和超声*** |
CN111050655A (zh) * | 2017-05-28 | 2020-04-21 | 利兰斯坦福初级大学董事会 | 通过非线性定位进行的超声成像 |
CN111819467A (zh) * | 2018-01-25 | 2020-10-23 | 海浪科技有限公司 | 波传播和散射参数的估计方法和仪器 |
CN113424073A (zh) * | 2018-12-17 | 2021-09-21 | 海浪科技有限公司 | 材料非线性体积弹性的超声估算 |
CN114089285A (zh) * | 2022-01-24 | 2022-02-25 | 安徽京淮健锐电子科技有限公司 | 一种基于一阶脉冲重复间隔pri的信号分选方法 |
US11885917B2 (en) | 2019-01-25 | 2024-01-30 | Surf Technology As | Methods and instrumentation for estimation of wave propagation and scattering parameters |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8923092B2 (en) | 2010-11-12 | 2014-12-30 | Chevron U.S.A. Inc. | System and method for investigating sub-surface features of a rock formation with acoustic sources generating coded signals |
EP2763585A1 (en) * | 2011-10-03 | 2014-08-13 | Surf Technology AS | Nonlinear imaging with dual band pulse complexes |
US10013767B2 (en) | 2013-11-01 | 2018-07-03 | The Research Foundation For The State University Of New York | Method for measuring the interior three-dimensional movement, stress and strain of an object |
CN105301582B (zh) * | 2015-10-14 | 2017-10-31 | 哈尔滨工程大学 | 一种水下声场中的弱信号增强方法 |
LT3544745T (lt) | 2016-11-22 | 2021-12-10 | Surf Technology As | Daugiajuostis ultragarso keitiklis |
DE112018003501T5 (de) | 2017-07-09 | 2020-04-23 | The Board Of Trustees Of The Leland Stanford Junior University | Ultraschallbildgebung mit spektraler compoundierung zur speckle-reduzierung |
CN111107791B (zh) * | 2017-09-01 | 2023-08-11 | 富士胶片索诺声公司 | 双频平面超声成像*** |
US10859695B2 (en) * | 2017-11-28 | 2020-12-08 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Acoustic system and method for characterizing granular media |
WO2020260944A1 (en) * | 2019-06-25 | 2020-12-30 | Surf Technology As | Suppression of multiple scattering noise in pulse echo imaging |
CN111024676B (zh) * | 2020-01-10 | 2022-03-04 | 河南工程学院 | 一种非线性z扫描测量方法及装置 |
CN113759354B (zh) * | 2020-06-02 | 2024-02-09 | 中国科学院声学研究所 | 一种适用于侧扫声呐的自适应底混响抑制方法 |
US20220105363A1 (en) | 2020-10-07 | 2022-04-07 | Surf Technology As | Methods and System for Stimulating Immune Response Against an Existing Cancer in a Patient |
US20220132240A1 (en) * | 2020-10-23 | 2022-04-28 | Alien Sandbox, LLC | Nonlinear Mixing of Sound Beams for Focal Point Determination |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5601086A (en) * | 1995-05-12 | 1997-02-11 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Beat frequency ultrasonic microsphere contrast agent detection system |
US6123669A (en) * | 1998-05-13 | 2000-09-26 | Kabushiki Kaisha Toshiba | 3D ultrasound imaging using 2D array |
WO2002029433A2 (en) * | 2000-10-02 | 2002-04-11 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging of nonlinearly intermodulated and harmonic frequency components |
US7301852B2 (en) * | 2003-08-13 | 2007-11-27 | Baker Hughes Incorporated | Methods of generating directional low frequency acoustic signals and reflected signal detection enhancements for seismic while drilling applications |
US20100036244A1 (en) * | 2008-01-09 | 2010-02-11 | Angelsen Bjoern A J | Nonlinear Elastic Wave Measurement and Imaging with Two-Frequency Elastic Wave Pulse Complexes |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5521882A (en) * | 1993-11-19 | 1996-05-28 | Schlumberger Technology Corporation | Measurement of formation characteristics using acoustic borehole tool having sources of different frequencies |
US5736642A (en) * | 1997-01-08 | 1998-04-07 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Nonlinear ultrasonic scanning to detect material defects |
CN1666568A (zh) | 2001-01-05 | 2005-09-07 | 比约恩·A·J·安杰尔森 | 宽带换能器 |
US7727156B2 (en) | 2005-07-26 | 2010-06-01 | Angelsen Bjoern A J | Dual frequency band ultrasound transducer arrays |
EA014167B1 (ru) * | 2005-09-08 | 2010-10-29 | Бьорн А.Й. Ангельсен | Способ построения изображения акустических свойств путем нелинейного низкочастотного управления свойствами высокочастотного рассеяния и распространения (варианты) и акустическое устройство для осуществления этого способа |
US8345509B2 (en) * | 2009-04-16 | 2013-01-01 | Chevron U.S.A., Inc. | System and method to create three-dimensional images of non-linear acoustic properties in a region remote from a borehole |
-
2010
- 2010-11-09 GB GB1018933.0A patent/GB2484753B/en not_active Expired - Fee Related
-
2011
- 2011-08-19 EP EP11817271.7A patent/EP2606374A2/en not_active Withdrawn
- 2011-08-19 US US13/213,965 patent/US8793079B2/en active Active
- 2011-08-19 CN CN201180050704.4A patent/CN103261909B/zh not_active Expired - Fee Related
- 2011-08-19 WO PCT/GB2011/051564 patent/WO2012035312A2/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5601086A (en) * | 1995-05-12 | 1997-02-11 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Beat frequency ultrasonic microsphere contrast agent detection system |
US6123669A (en) * | 1998-05-13 | 2000-09-26 | Kabushiki Kaisha Toshiba | 3D ultrasound imaging using 2D array |
WO2002029433A2 (en) * | 2000-10-02 | 2002-04-11 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging of nonlinearly intermodulated and harmonic frequency components |
US7301852B2 (en) * | 2003-08-13 | 2007-11-27 | Baker Hughes Incorporated | Methods of generating directional low frequency acoustic signals and reflected signal detection enhancements for seismic while drilling applications |
US20100036244A1 (en) * | 2008-01-09 | 2010-02-11 | Angelsen Bjoern A J | Nonlinear Elastic Wave Measurement and Imaging with Two-Frequency Elastic Wave Pulse Complexes |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107231811B (zh) * | 2014-12-10 | 2021-07-16 | 瑟夫技术公司 | 用于非线***互散射的成像的方法 |
CN107231811A (zh) * | 2014-12-10 | 2017-10-03 | 瑟夫技术公司 | 用于非线***互散射的成像的方法 |
CN107530046A (zh) * | 2015-01-15 | 2018-01-02 | 罗德尼·赫林 | 漫射声学共焦成像器 |
CN107530046B (zh) * | 2015-01-15 | 2021-07-13 | 罗德尼·赫林 | 漫射声学共焦成像器 |
CN108338808B (zh) * | 2017-01-23 | 2024-05-24 | 美国西门子医疗解决公司 | 使用相干的剪切速度成像 |
CN108338808A (zh) * | 2017-01-23 | 2018-07-31 | 美国西门子医疗解决公司 | 使用相干的剪切速度成像 |
CN111050655A (zh) * | 2017-05-28 | 2020-04-21 | 利兰斯坦福初级大学董事会 | 通过非线性定位进行的超声成像 |
US11369343B2 (en) | 2017-05-28 | 2022-06-28 | The Board Of Trustees Of The Leland Stanford Junior University | Ultrasound imaging by nonlinear localization |
CN109124680A (zh) * | 2017-06-16 | 2019-01-04 | 美国西门子医疗解决公司 | 用于形成超声探头的超声信号的方法和超声*** |
CN109124680B (zh) * | 2017-06-16 | 2021-09-07 | 美国西门子医疗解决公司 | 用于形成超声探头的超声信号的方法和超声*** |
US11303270B2 (en) | 2017-06-16 | 2022-04-12 | Siemens Medical Solutions Usa, Inc. | Pulse cancellation compensation method for contrast agent imaging |
CN107797160A (zh) * | 2017-09-01 | 2018-03-13 | 上海交通大学 | 弹性波和电磁波ct勘测数据联合分析***及方法 |
CN111819467A (zh) * | 2018-01-25 | 2020-10-23 | 海浪科技有限公司 | 波传播和散射参数的估计方法和仪器 |
CN111819467B (zh) * | 2018-01-25 | 2024-05-07 | 海浪科技有限公司 | 波传播和散射参数的估计方法和仪器 |
CN113424073A (zh) * | 2018-12-17 | 2021-09-21 | 海浪科技有限公司 | 材料非线性体积弹性的超声估算 |
US11885917B2 (en) | 2019-01-25 | 2024-01-30 | Surf Technology As | Methods and instrumentation for estimation of wave propagation and scattering parameters |
CN114089285A (zh) * | 2022-01-24 | 2022-02-25 | 安徽京淮健锐电子科技有限公司 | 一种基于一阶脉冲重复间隔pri的信号分选方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2606374A2 (en) | 2013-06-26 |
WO2012035312A3 (en) | 2012-09-27 |
US20120095699A1 (en) | 2012-04-19 |
GB2484753B (en) | 2013-01-02 |
US8793079B2 (en) | 2014-07-29 |
WO2012035312A2 (en) | 2012-03-22 |
CN103261909B (zh) | 2016-01-20 |
GB201018933D0 (en) | 2010-12-22 |
GB2484753A (en) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103261909B (zh) | 用于成像非线性相互作用散射的方法 | |
US9250323B2 (en) | Target detection utilizing image array comparison | |
JP6408297B2 (ja) | ビームフォーミング方法、計測イメージング装置、及び、通信装置 | |
US10324063B2 (en) | Methods and systems for measuring properties with ultrasound | |
CN108226927B (zh) | 基于加权迭代最小稀疏贝叶斯重构算法的sar成像方法 | |
Zhou et al. | Subsurface cavity imaging by crosshole borehole radar measurements | |
CN103260525A (zh) | 针对背景运动影响调节声辐射力效应的测量结果 | |
WO2013051943A1 (en) | Nonlinear imaging with dual band pulse complexes | |
US11209530B2 (en) | Method for imaging of nonlinear interaction scattering | |
US11841265B2 (en) | Heterogeneous subsurface imaging systems and methods | |
Yan et al. | Ultra-wideband MIMO ambiguity function and its factorability | |
Klibanov et al. | Convexification inversion method for nonlinear SAR imaging with experimentally collected data | |
Duric et al. | Acoustic tomography: promise versus reality | |
CN102253371B (zh) | 一种用于探地雷达成像的散射强度加权处理方法 | |
JP2009186449A (ja) | 2次元画像データ生成方法及び探査装置 | |
CN110297237B (zh) | 考虑天线方向图的探地雷达绕射叠加成像方法及*** | |
CN104280737A (zh) | 加权宽带时反算子分解声成像方法 | |
EP4194893A1 (en) | Method and system for processing beamformed data | |
Lei | A time-domain beamformer for UWB through-wall imaging | |
Taipodia et al. | A review of active and passive MASW techniques | |
Muir et al. | Seismo-acoustic sonar for buried object detection | |
Aftanas | Through wall imaging using M-sequence UWB radar system | |
Carlson et al. | High resolution image reconstruction from full-matrix capture data using minimum mean square error deconvolution of the spatio-temporal system transfer function | |
Nguyen et al. | 3D SAR image formation for underground targets using ultra-wideband (UWB) radar | |
Muir et al. | Comparison of acoustic and seismic excitation, propagation, and scattering at an air-ground interface containing a mine-like inclusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20180819 |
|
CF01 | Termination of patent right due to non-payment of annual fee |