[go: up one dir, main page]

CN117969677B - Acoustic field detection method, system, device and medium based on coded phase imaging - Google Patents

Acoustic field detection method, system, device and medium based on coded phase imaging Download PDF

Info

Publication number
CN117969677B
CN117969677B CN202410364077.0A CN202410364077A CN117969677B CN 117969677 B CN117969677 B CN 117969677B CN 202410364077 A CN202410364077 A CN 202410364077A CN 117969677 B CN117969677 B CN 117969677B
Authority
CN
China
Prior art keywords
light field
distribution
wavefront
sound field
plane
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202410364077.0A
Other languages
Chinese (zh)
Other versions
CN117969677A (en
Inventor
徐英明
施钧辉
孙明丽
陈睿黾
阮永都
李驰野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Lab
Original Assignee
Zhejiang Lab
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang Lab filed Critical Zhejiang Lab
Priority to CN202410364077.0A priority Critical patent/CN117969677B/en
Publication of CN117969677A publication Critical patent/CN117969677A/en
Application granted granted Critical
Publication of CN117969677B publication Critical patent/CN117969677B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/22Details, e.g. general constructional or apparatus details
    • G01N29/24Probes
    • G01N29/2418Probes using optoacoustic interaction with the material, e.g. laser radiation, photoacoustics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/22Details, e.g. general constructional or apparatus details

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Optics & Photonics (AREA)
  • Electrostatic, Electromagnetic, Magneto- Strictive, And Variable-Resistance Transducers (AREA)

Abstract

The invention discloses a sound field detection method, a system, equipment and a medium based on code phase imaging, wherein the system comprises the following components: the parallel light field transmitting module is used for generating a pulse parallel light field to illuminate a sound field medium, the pulse parallel light field carries out wave front modulation through the sound field medium, imaging is carried out through the imaging module, modulation is carried out through the coded plate, and light field intensity distribution is obtained through the light field detector; the system further comprises: the ultrasonic transducer module is used for generating a to-be-detected high-frequency sound field; the control and calculation module is used for controlling the synchronization of the parallel light field emission module and the light field detector, and storing the light field intensity distribution acquired by the light field detector to calculate the complex amplitude distribution of the high-frequency sound field to be detected. The invention can indirectly detect the parameter information of the sound field by utilizing the complex amplitude change of the light field after passing through the sound field medium, and the method has the advantages of simultaneously detecting the intensity and the phase of the light field, and has the characteristics of high stability, high spatial resolution and high phase resolution.

Description

基于编码相位成像的声场探测方法、系统、设备、介质Acoustic field detection method, system, device and medium based on coded phase imaging

技术领域Technical Field

本发明属于光场相位成像和声场探测领域,具体涉及一种基于编码相位成像的声场探测方法、系统、设备、介质。The present invention belongs to the field of light field phase imaging and sound field detection, and in particular relates to a sound field detection method, system, equipment and medium based on coded phase imaging.

背景技术Background technique

传统的声场探测方法一般采用超声探头逐点扫描方法进行探测,这种方法基于超声波在被检测材料中的传播特性,利用超声探头在材料表面或内部进行逐点扫描,通过分析超声波的反射、透射或散射等信息,来检测材料的内部结构和缺陷,然而该方法需要长时间扫描,并且干扰了待测声场的分布。Traditional sound field detection methods generally use an ultrasonic probe to scan point by point. This method is based on the propagation characteristics of ultrasonic waves in the material being detected. It uses an ultrasonic probe to scan point by point on the surface or inside of the material. By analyzing the reflection, transmission or scattering information of the ultrasonic wave, the internal structure and defects of the material are detected. However, this method requires a long scanning time and interferes with the distribution of the sound field to be measured.

目前还可以通过记录透射光场条纹强度分布判断声场的分布参数,该方法是一种光学与声学相结合的技术。这种方法通常涉及到声光相互作用的原理,其中声波会对光场产生调制作用,使得透射光场产生特定的条纹强度分布。具体来说,当声波传播通过介质时,它会引起介质密度的变化,这种变化进而会改变介质的折射率。当光场穿过这个受到声波调制的介质时,光的传播路径和相位会受到影响,从而导致透射光场的强度分布发生变化。但是该方法探测的为经过声场介质后衍射的光场强度分布,无法得到声场介质的密度分布,受光源强度及介质衍射影响,误差较大。At present, the distribution parameters of the acoustic field can also be determined by recording the intensity distribution of the transmitted light field stripes. This method is a technology that combines optics and acoustics. This method usually involves the principle of acoustic-optical interaction, in which the sound wave modulates the light field, causing the transmitted light field to produce a specific stripe intensity distribution. Specifically, when the sound wave propagates through the medium, it causes a change in the density of the medium, which in turn changes the refractive index of the medium. When the light field passes through this medium modulated by the sound wave, the propagation path and phase of the light will be affected, resulting in a change in the intensity distribution of the transmitted light field. However, this method detects the intensity distribution of the light field diffracted after passing through the acoustic field medium, and cannot obtain the density distribution of the acoustic field medium. It is affected by the intensity of the light source and the diffraction of the medium, and has a large error.

发明内容Summary of the invention

本发明针对上述声场探测技术的局限性,提出了一种基于编码相位成像的声场探测方法、系统、设备、介质。In view of the limitations of the above-mentioned sound field detection technology, the present invention proposes a sound field detection method, system, device and medium based on coded phase imaging.

根据本发明实施例第一方面,提供了一种基于编码相位成像的声场探测系统,所述系统包括:According to a first aspect of an embodiment of the present invention, a sound field detection system based on coded phase imaging is provided, the system comprising:

平行光场发射模块,所述平行光场发射模块用于产生脉冲平行光场以照明声场介质,脉冲平行光场经声场介质进行波前调制,经成像模块成像,经编码板调制,通过光场探测器获取光场强度分布;A parallel light field emission module, which is used to generate a pulse parallel light field to illuminate the acoustic field medium. The pulse parallel light field is wavefront modulated by the acoustic field medium, imaged by the imaging module, modulated by the encoding plate, and the light field intensity distribution is obtained by the light field detector;

所述系统还包括:The system further comprises:

超声换能器模块,用于产生待测高频声场;Ultrasonic transducer module, used to generate high-frequency sound field to be measured;

控制及计算模块,用于控制平行光场发射模块和光场探测器的同步,存储光场探测器获取的光场强度分布以计算待测的高频声场的复振幅分布。The control and calculation module is used to control the synchronization of the parallel light field emission module and the light field detector, and store the light field intensity distribution obtained by the light field detector to calculate the complex amplitude distribution of the high-frequency sound field to be measured.

根据本发明实施例第二方面,提供了一种基于编码相位成像的声场探测方法,基于上述的基于编码相位成像的声场探测系统实现,所述方法包括:According to a second aspect of an embodiment of the present invention, a sound field detection method based on coded phase imaging is provided, which is implemented based on the above-mentioned sound field detection system based on coded phase imaging. The method includes:

获取超声换能器模块产生的待测高频声场;Acquire the high-frequency sound field to be measured generated by the ultrasonic transducer module;

获取光场探测器记录的光场强度分布;Obtaining the light field intensity distribution recorded by the light field detector;

将光场探测器所在的平面设置为第一平面,编码板所在的平面为第二平面;设置编码板的前表面处初始的光场复振幅分布;初始的光场复振幅分布在第一平面、第二平面、约束面间迭代传播;基于光场强度分布计算误差,当误差小于阈值时完成迭代,得到编码板前表面处的波前分布;The plane where the light field detector is located is set as the first plane, and the plane where the encoding plate is located is set as the second plane; the initial light field complex amplitude distribution at the front surface of the encoding plate is set; the initial light field complex amplitude distribution is iteratively propagated between the first plane, the second plane, and the constraint surface; the error is calculated based on the light field intensity distribution, and the iteration is completed when the error is less than a threshold value to obtain the wavefront distribution at the front surface of the encoding plate;

将编码板前表面处的波前分布回传至超声换能器模块所在的平面,从而计算将编码板前表面处的波前分布在超声换能器模块处的波前振幅和相位分布,得到声场的复振幅分布。The wavefront distribution at the front surface of the coding plate is transmitted back to the plane where the ultrasonic transducer module is located, so as to calculate the wavefront amplitude and phase distribution of the wavefront distribution at the front surface of the coding plate at the ultrasonic transducer module to obtain the complex amplitude distribution of the sound field.

根据本发明实施例第三方面,提供了一种电子设备,包括存储器和处理器,所述存储器与所述处理器耦接;其中,所述存储器用于存储程序数据,所述处理器用于执行所述程序数据以实现上述的基于编码相位成像的声场探测方法。According to a third aspect of an embodiment of the present invention, an electronic device is provided, comprising a memory and a processor, wherein the memory is coupled to the processor; wherein the memory is used to store program data, and the processor is used to execute the program data to implement the above-mentioned sound field detection method based on coded phase imaging.

根据本发明实施例第四方面,提供了一种计算机可读存储介质,其上存储有计算机程序,所述程序被处理器执行时实现上述的基于编码相位成像的声场探测方法。According to a fourth aspect of an embodiment of the present invention, a computer-readable storage medium is provided, on which a computer program is stored. When the program is executed by a processor, the above-mentioned sound field detection method based on coded phase imaging is implemented.

与现有技术相比,本发明的技术效果为:Compared with the prior art, the technical effects of the present invention are:

本发明提供了一种基于编码相位成像的声场探测系统与方法,将超声换能器放置于脉冲平行光场中,基于相干调制成像单次重建光场的复振幅分布,结合超短脉冲探测光场的瞬时变化,同时可以实现连续的实时探测。同时,本发明通过重建声场位置的光场相位分布,间接得到声场介质的密度分布,根据声场介质的密度分布,可以判断超声换能器产生的待测高频声场的分布参数,从而实现无损声场探测,在不打乱原始声场的前提下实现声场分布参数的探测。并且,本发明系统光路简单,性能稳定,成本低。The present invention provides a sound field detection system and method based on coded phase imaging, which places an ultrasonic transducer in a pulsed parallel light field, reconstructs the complex amplitude distribution of the light field based on coherent modulation imaging, and combines ultrashort pulses to detect instantaneous changes in the light field, thereby achieving continuous real-time detection. At the same time, the present invention indirectly obtains the density distribution of the sound field medium by reconstructing the light field phase distribution at the sound field position. According to the density distribution of the sound field medium, the distribution parameters of the high-frequency sound field to be measured generated by the ultrasonic transducer can be determined, thereby achieving non-destructive sound field detection, and realizing the detection of sound field distribution parameters without disrupting the original sound field. In addition, the system of the present invention has a simple optical path, stable performance, and low cost.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

图1为本发明实施例提供的基于编码相位成像的声场探测系统的示意图;FIG1 is a schematic diagram of an acoustic field detection system based on coded phase imaging provided by an embodiment of the present invention;

图2为本发明实施例提供的超声换能器模块以及二维平移模块的示意图;FIG2 is a schematic diagram of an ultrasonic transducer module and a two-dimensional translation module provided in an embodiment of the present invention;

图3为本发明实施例提供的在空间脉冲激光器下基于编码相位成像的声场探测系统的示意图;3 is a schematic diagram of an acoustic field detection system based on coded phase imaging under a spatial pulse laser provided by an embodiment of the present invention;

图4为本发明实施例提供的在光纤脉冲激光器下基于编码相位成像的声场探测系统的示意图;FIG4 is a schematic diagram of an acoustic field detection system based on coded phase imaging under a fiber pulse laser provided by an embodiment of the present invention;

图5为本发明实施例提供的基于编码相位成像的声场探测方法的流程示意图;FIG5 is a schematic flow chart of a method for detecting an acoustic field based on coded phase imaging according to an embodiment of the present invention;

图6为本发明实施例提供的获取编码板前表面处的波前分布的流程示意图;6 is a schematic diagram of a process for obtaining the wavefront distribution at the front surface of the coding plate according to an embodiment of the present invention;

图7为本发明实施例提供的在光纤脉冲激光器下重建声场介质的振幅和相位分布;FIG7 is an amplitude and phase distribution of an acoustic field medium reconstructed under a fiber pulse laser provided by an embodiment of the present invention;

图8为本发明实施例提供的一种电子设备的示意图。FIG. 8 is a schematic diagram of an electronic device provided by an embodiment of the present invention.

图中,1-平行光场发射模块、2-声场介质模块、3-超声换能器模块、4-成像模块、5-编码板、6-光场探测器、7-控制及计算模块、8-二维平移模块、9-空间脉冲激光器、10-空间光滤波器、11-准直透镜、12-消色差成像透镜、13-光纤脉冲激光器、14-光纤跳线。In the figure, 1-parallel light field emission module, 2-acoustic field medium module, 3-ultrasonic transducer module, 4-imaging module, 5-encoding board, 6-light field detector, 7-control and calculation module, 8-two-dimensional translation module, 9-spatial pulse laser, 10-spatial light filter, 11-collimating lens, 12-achromatic imaging lens, 13-fiber pulse laser, 14-fiber jumper.

具体实施方式Detailed ways

这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present invention. Instead, they are merely examples of devices and methods consistent with some aspects of the present invention as detailed in the appended claims.

在本发明使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本发明。在本发明和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。The terms used in the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. The singular forms "a", "the" and "the" used in the present invention and the appended claims are also intended to include plural forms, unless the context clearly indicates other meanings. It should also be understood that the term "and/or" used herein refers to and includes any or all possible combinations of one or more associated listed items.

应当理解,尽管在本发明可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本发明范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。It should be understood that although the terms first, second, third, etc. may be used in the present invention to describe various information, these information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of the present invention, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining".

如图1所示,本发明提供了一种基于编码相位成像的声场探测系统,所述系统包括:平行光场发射模块1,所述平行光场发射模块1用于产生脉冲平行光场以照明声场介质,脉冲平行光场经声场介质2进行波前调制,经成像模块4成像,经编码板5调制,通过光场探测器6获取光场强度分布;As shown in FIG1 , the present invention provides an acoustic field detection system based on coded phase imaging, the system comprising: a parallel light field emission module 1, the parallel light field emission module 1 is used to generate a pulse parallel light field to illuminate an acoustic field medium, the pulse parallel light field is wavefront modulated by an acoustic field medium 2, imaged by an imaging module 4, modulated by a coding plate 5, and the light field intensity distribution is obtained by a light field detector 6;

所述系统还包括:The system further comprises:

超声换能器模块3,用于产生待测高频声场;Ultrasonic transducer module 3, used to generate a high-frequency sound field to be measured;

控制及计算模块7,用于控制平行光场发射模块1和光场探测器6的同步,存储光场探测器6获取的光场强度分布以计算待测的高频声场的复振幅分布。The control and calculation module 7 is used to control the synchronization of the parallel light field transmitting module 1 and the light field detector 6, and store the light field intensity distribution obtained by the light field detector 6 to calculate the complex amplitude distribution of the high-frequency sound field to be measured.

进一步地,所述系统还包括:如图2所示,二维平移模块8,用于在x方向和y方向上移动超声换能器模块3。Furthermore, the system further comprises: as shown in FIG. 2 , a two-dimensional translation module 8 for moving the ultrasonic transducer module 3 in the x-direction and the y-direction.

进一步地,所述平行光场发射模块1选用空间激光器产生平行光场系统或光纤激光器产生平行光场系统;所述平行光场发射模块1产生的脉冲平行光场的脉宽小于10ns,重复频率小于1kHz。Furthermore, the parallel light field emission module 1 uses a space laser to generate a parallel light field system or a fiber laser to generate a parallel light field system; the pulse width of the pulse parallel light field generated by the parallel light field emission module 1 is less than 10ns, and the repetition frequency is less than 1kHz.

进一步地,所述成像模块4采用缩束成像系统,用于低频大口径声场探测;或,所述成像模块4采用成像放大系统,用于高频小口径声场探测。Furthermore, the imaging module 4 adopts a reduced beam imaging system for low-frequency and large-caliber sound field detection; or, the imaging module 4 adopts an imaging amplification system for high-frequency and small-caliber sound field detection.

实施例1Example 1

如图3所示,本实例提供了一种基于编码相位成像的声场探测系统,其中平行光场发射模块1选用空间激光器10;具体地,所述系统包括:依次布置的空间脉冲激光器9、空间光滤波器10、准直透镜11、声场介质2、消色差成像透镜12、编码板5、光场探测器6;As shown in FIG3 , this example provides an acoustic field detection system based on coded phase imaging, wherein the parallel light field emission module 1 uses a spatial laser 10; specifically, the system comprises: a spatial pulse laser 9, a spatial light filter 10, a collimating lens 11, an acoustic field medium 2, an achromatic imaging lens 12, a coding plate 5, and an optical field detector 6, which are arranged in sequence;

其中,所述空间脉冲激光器9,用于产生空间脉冲光场,用于照明瞬态声场介质;在本实例中,所述空间脉冲激光器9的光源波长为1064nm,脉宽为10ns,重复频率为10Hz。The spatial pulse laser 9 is used to generate a spatial pulse light field for illuminating a transient acoustic field medium. In this example, the light source wavelength of the spatial pulse laser 9 is 1064 nm, the pulse width is 10 ns, and the repetition frequency is 10 Hz.

空间光滤波器10,用于滤除杂散光,产生球面脉冲光场Q;A spatial optical filter 10 is used to filter out stray light and generate a spherical pulse light field Q;

准直透镜11,用于将球面脉冲光场Q准直为平行光场,用于照明声场介质。The collimating lens 11 is used to collimate the spherical pulse light field Q into a parallel light field for illuminating the acoustic field medium.

声场介质模块2,用于传播声场,对平行光场产生波前调制;The acoustic field medium module 2 is used to propagate the acoustic field and generate wavefront modulation on the parallel light field;

超声换能器模块3,用于产生待测的高频声场;Ultrasonic transducer module 3, used to generate a high-frequency sound field to be measured;

所述成像模块4采用消色差成像透镜12,用于对声场介质进行成像。在本实例中,所述消色差成像透镜12,焦距为100mm,口径为50.8mm,成像放大倍率为3。The imaging module 4 uses an achromatic imaging lens 12 for imaging the sound field medium. In this example, the achromatic imaging lens 12 has a focal length of 100 mm, an aperture of 50.8 mm, and an imaging magnification of 3.

在本实例中,所述编码板5的单位像元尺寸为9μm×9μm,对于波长1064nm,相位变化为0-π分布;光场探测器6的单位像素尺寸为9μm×9μm,动态范围为16bit。In this example, the unit pixel size of the encoding plate 5 is 9 μm×9 μm, and for a wavelength of 1064 nm, the phase change is a 0-π distribution; the unit pixel size of the light field detector 6 is 9 μm×9 μm, and the dynamic range is 16 bits.

实施例2Example 2

如图4所示,本实例提供了一种基于编码相位成像的声场探测系统,其中平行光场发射模块1选用光纤脉冲激光器13;具体地,所述系统包括:依次布置的光纤脉冲激光器13、光纤跳线14、准直透镜11、声场介质2、消色差成像透镜12、编码板5、光场探测器6。As shown in FIG4 , this example provides an acoustic field detection system based on coded phase imaging, wherein the parallel light field emission module 1 uses a fiber pulse laser 13; specifically, the system includes: a fiber pulse laser 13, a fiber jumper 14, a collimating lens 11, an acoustic field medium 2, an achromatic imaging lens 12, a coding plate 5, and a light field detector 6 arranged in sequence.

其中,所述光纤脉冲激光器13,用于产生光纤脉冲光场,用于照明瞬态变化声场介质;在本实例中,所述光纤脉冲激光器13中光源波长为532nm,脉宽为10ps,重复频率为1kHz。The fiber pulse laser 13 is used to generate a fiber pulse light field for illuminating a transiently changing acoustic field medium. In this example, the fiber pulse laser 13 has a light source wavelength of 532 nm, a pulse width of 10 ps, and a repetition frequency of 1 kHz.

光纤跳线14,用于产生传输光场,并最终产生球面光场W;The optical fiber jumper 14 is used to generate a transmission light field and finally generate a spherical light field W;

准直透镜11,用于将球面脉冲光场W准直为平行光场,用于照明声场介质;A collimating lens 11, used to collimate the spherical pulse light field W into a parallel light field for illuminating the acoustic field medium;

声场介质模块2,用于传播声场,对平行光场产生波前调制;The acoustic field medium module 2 is used to propagate the acoustic field and generate wavefront modulation on the parallel light field;

超声换能器模块3,用于产生待测的高频声场;Ultrasonic transducer module 3, used to generate a high-frequency sound field to be measured;

所述成像模块4采用消色差成像透镜12,用于对声场介质进行成像。在本实例中,所述的消色差成像透镜12,焦距为150mm,口径为50.8mm,缩小倍率为3。The imaging module 4 uses an achromatic imaging lens 12 for imaging the sound field medium. In this example, the achromatic imaging lens 12 has a focal length of 150 mm, an aperture of 50.8 mm, and a reduction ratio of 3.

编码板5,用于对光场波前进行调制,算法迭代时其复振幅透过率需要已知,本实施例中采用单位像元尺寸为18μm×18μm,对于波长532nm,相位变化为0-π分布;The coding plate 5 is used to modulate the wavefront of the light field. Its complex amplitude transmittance needs to be known during algorithm iteration. In this embodiment, the unit pixel size is 18 μm×18 μm. For a wavelength of 532 nm, the phase change is a 0-π distribution.

光场探测器6,用于记录光场的强度分布,记录该强度矩阵,本实施例中采用的单位像素尺寸为9μm×9μm,动态范围为16bit。The light field detector 6 is used to record the intensity distribution of the light field and the intensity matrix. In this embodiment, the unit pixel size used is 9 μm×9 μm and the dynamic range is 16 bits.

如图5所示,本发明实施例还提供了一种基于编码相位成像的声场探测方法,基于上述的基于编码相位成像的声场探测系统实现,所述方法包括:As shown in FIG5 , an embodiment of the present invention further provides a sound field detection method based on coded phase imaging, which is implemented based on the above-mentioned sound field detection system based on coded phase imaging. The method includes:

步骤S1,超声换能器模块3产生待测高频声场。Step S1, the ultrasonic transducer module 3 generates a high-frequency sound field to be measured.

具体地,利用二维平移模块8将超声换能器模块3平移至光场可视范围内,启动超声换能器模块3发出待测高频声场。Specifically, the two-dimensional translation module 8 is used to translate the ultrasonic transducer module 3 into the visible range of the light field, and the ultrasonic transducer module 3 is started to emit a high-frequency sound field to be measured.

步骤S2,获取光场探测器6记录的光场强度分布。Step S2, obtaining the light field intensity distribution recorded by the light field detector 6.

具体地,利用控制及计算模块7控制空间脉冲激光器9和光场探测器6时间同步,并 且存储光场探测器6记录的光场强度分布矩阵Specifically, the control and calculation module 7 is used to control the time synchronization of the spatial pulse laser 9 and the light field detector 6, and store the light field intensity distribution matrix recorded by the light field detector 6. .

步骤S3,将光场探测器6所在的平面设置为第一平面,编码板5所在的平面为第二平面;设置编码板5的前表面处初始的光场复振幅分布;初始的光场复振幅分布在第一平面、第二平面、约束面间迭代传播;基于光场强度分布计算误差,当误差小于阈值时完成迭代,得到编码板5前表面处的波前分布。Step S3, setting the plane where the light field detector 6 is located as the first plane, and the plane where the encoding plate 5 is located as the second plane; setting the initial light field complex amplitude distribution at the front surface of the encoding plate 5; iteratively propagating the initial light field complex amplitude distribution between the first plane, the second plane, and the constraint surface; calculating the error based on the light field intensity distribution, completing the iteration when the error is less than a threshold, and obtaining the wavefront distribution at the front surface of the encoding plate 5.

具体地,如图6所示,所述步骤S3具体包括以下子步骤:Specifically, as shown in FIG6 , step S3 specifically includes the following sub-steps:

步骤S301,设置编码板5的前表面处初始的光场复振幅分布为,其中A和分别为随机猜测的振幅和相位分布,n为对 应时间序号,为对应光场探测器6的记录时间。 Step S301, setting the initial light field complex amplitude distribution at the front surface of the coding plate 5 to , where A and are the amplitude and phase distributions of random guesses, n is the corresponding time sequence number, is the recording time of the corresponding light field detector 6.

步骤S302,将编码板5的前表面处初始的光场复振幅分布作为第1次迭代值Step S302: taking the initial light field complex amplitude distribution at the front surface of the coding plate 5 as the first iteration value .

步骤S303,光场复振幅分布经编码板5调制后,得到第一波前,表达式如下:Step S303, after the complex amplitude distribution of the light field is modulated by the encoding plate 5, a first wavefront is obtained, which is expressed as follows:

;

式中,k为迭代次数,k=1,2,3…,为编码板5的透过率复振幅分布; Where k is the number of iterations, k=1,2,3…, is the transmittance complex amplitude distribution of the coding plate 5;

步骤S304,第一波前在自由空间中传播至第一平面时,利用光场强度分布进行约束,得到第二波前,表达式如下:Step S304, when the first wavefront propagates to the first plane in the free space, the light field intensity distribution is used for constraint to obtain the second wavefront, which is expressed as follows:

;

;

式中,L1表示第二平面与第一平面之间的距离,表示光场自由空间传播L距 离;表示对波前取相位角; Where L1 represents the distance between the second plane and the first plane, represents the free space propagation distance L of the light field; It means taking the phase angle of the wavefront;

步骤S305,将第二波前回传至第二平面,并利用透过率复振幅分布进行更新,得到编码板5前表面处的第三波前,表达式如下:Step S305, the second wavefront is transmitted back to the second plane, and updated using the transmittance complex amplitude distribution to obtain the third wavefront at the front surface of the encoding plate 5, which is expressed as follows:

;

;

步骤S306,将第三波前传输至约束面,所述约束面为第三波前的频谱面或焦点面;构造孔径约束函数,利用孔径约束函数对时间和波长进行约束,得到第四波前,表达式如下:Step S306, transmitting the third wavefront to a constraint surface, which is a spectrum surface or a focal surface of the third wavefront; constructing an aperture constraint function, and using the aperture constraint function to constrain time and wavelength to obtain a fourth wavefront, which is expressed as follows:

;

式中,L2表示第二平面与约束面之间的距离,为孔径函数约束; Where L2 represents the distance between the second plane and the constraint surface, is the aperture function constraint;

步骤S307,将第四波前传输至第二平面,得到第五波前,表达式如下:Step S307, transmitting the fourth wavefront to the second plane to obtain the fifth wavefront, which is expressed as follows:

;

步骤S308,按上述过程进行迭代,并基于第二波前与脉冲衍射光斑强度矩阵计算误差,当误差小于阈值时完成迭代,将迭代得到的第五波前作为编码板5前表面处的波前分布。Step S308, iterate according to the above process, and calculate the error based on the second wavefront and the pulse diffraction spot intensity matrix. When the error is less than the threshold, the iteration is completed, and the fifth wavefront obtained by the iteration is used as the wavefront distribution at the front surface of the encoding plate 5.

具体地,误差的表达式如下:Specifically, the error expression is as follows:

;

式中,表示光场探测器6记录的光场强度分布;表示第二波 前。 In the formula, represents the light field intensity distribution recorded by the light field detector 6; Indicates before the second wave.

步骤S4,将编码板5前表面处的波前分布回传至超声换能器模块3所在的平面,从而计算将编码板5前表面处的波前分布在超声换能器模块3处的波前振幅和相位分布,得到声场的复振幅分布。Step S4, transmit the wavefront distribution at the front surface of the coding plate 5 back to the plane where the ultrasonic transducer module 3 is located, so as to calculate the wavefront amplitude and phase distribution of the wavefront distribution at the front surface of the coding plate 5 at the ultrasonic transducer module 3, and obtain the complex amplitude distribution of the sound field.

具体地,获取编码板5前表面处的波前分布,记为Specifically, the wavefront distribution at the front surface of the coding plate 5 is obtained, which is recorded as ;

将编码板5前表面处的波前分布回传至超声换能器模块3所在的平面,得到超声换能器模块3所在的平面处的波前分布,表达式如下:The wavefront distribution at the front surface of the coding plate 5 is transmitted back to the plane where the ultrasonic transducer module 3 is located, and the wavefront distribution at the plane where the ultrasonic transducer module 3 is located is obtained, and the expression is as follows:

;

式中,其中为编码板5所在的平面为第二平面至超声换能器模块3所在的平面之 间的距离; In the formula, is the distance between the plane where the encoding plate 5 is located and the second plane and the plane where the ultrasonic transducer module 3 is located;

计算将编码板5前表面处的波前分布在超声换能器模块3处的波前振幅和相位分布;即于声场介质2的密度,得到声场的复振 幅分布。 Calculate the wavefront amplitude of the wavefront at the front surface of the encoder plate 5 distributed at the ultrasonic transducer module 3 and phase distribution ; That is, the density of the sound field medium 2 is used to obtain the complex amplitude distribution of the sound field.

如图7所示,图7中的(a)为声场位置光场的振幅和图7中的(b)为相位分布,其中振幅中无法分辨声场分布信息,但是该位置的相位条纹分布与声场的分布相同,通过计算条纹间隔可以得到该声场为1.96MHz。As shown in Figure 7, (a) in Figure 7 is the amplitude of the light field at the sound field position and (b) in Figure 7 is the phase distribution, where the sound field distribution information cannot be distinguished from the amplitude, but the phase fringe distribution at this position is the same as the distribution of the sound field. By calculating the fringe interval, it can be obtained that the sound field is 1.96MHz.

综上所述,本发明提供了一种基于编码相位成像的声场探测系统与方法,将超声换能器放置于脉冲平行光场中,基于相干调制成像单次重建光场的复振幅分布,结合超短脉冲探测光场的瞬时变化,同时可以实现连续的实时探测。同时,本发明通过重建声场位置的光场相位分布,间接得到声场介质的密度分布,根据声场介质的密度分布,可以判断超声换能器产生的待测高频声场的分布参数,从而实现无损声场探测,在不打乱原始声场的前提下实现声场分布参数的探测。并且,本发明系统光路简单,性能稳定,成本低。In summary, the present invention provides a sound field detection system and method based on coded phase imaging, which places an ultrasonic transducer in a pulsed parallel light field, reconstructs the complex amplitude distribution of the light field based on single coherent modulation imaging, and combines ultrashort pulses to detect instantaneous changes in the light field, while achieving continuous real-time detection. At the same time, the present invention indirectly obtains the density distribution of the sound field medium by reconstructing the light field phase distribution at the sound field position. According to the density distribution of the sound field medium, the distribution parameters of the high-frequency sound field to be measured generated by the ultrasonic transducer can be determined, thereby achieving non-destructive sound field detection and detecting the sound field distribution parameters without disrupting the original sound field. In addition, the system of the present invention has a simple optical path, stable performance, and low cost.

本说明书还提供了一种计算机可读存储介质,该存储介质存储有计算机程序,计算机程序可用于执行上述数据同步的方法。The present specification also provides a computer-readable storage medium, which stores a computer program, and the computer program can be used to execute the above-mentioned data synchronization method.

本说明书还提供了图8所示的电子设备的示意结构图。如图8所述,在硬件层面,该电子设备包括处理器、内部总线、网络接口、内存以及非易失性存储器,当然还可能包括其他业务所需要的硬件。处理器从非易失性存储器中读取对应的计算机程序到内存中然后运行,以实现上述数据同步的方法。This specification also provides a schematic structural diagram of an electronic device as shown in FIG8. As shown in FIG8, at the hardware level, the electronic device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory, and may also include hardware required for other services. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to implement the above-mentioned data synchronization method.

当然,除了软件实现方式之外,本说明书并不排除其他实现方式,比如逻辑器件抑或软硬件结合的方式等等,也就是说以下处理流程的执行主体并不限定于各个逻辑单元,也可以是硬件或逻辑器件。Of course, in addition to software implementation, this specification does not exclude other implementation methods, such as logic devices or a combination of software and hardware, etc., that is to say, the executor of the following processing flow is not limited to each logic unit, but can also be hardware or logic devices.

在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device, PLD)(例如现场可编程门阵列(Field ProgrammableGateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera HardwareDescriptionLanguage)、Confluence、CUPL(Cornell University ProgrammingLanguage)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。In the 1990s, it was very clear whether the improvement of a technology was hardware improvement (for example, improvement of the circuit structure of diodes, transistors, switches, etc.) or software improvement (improvement of the method flow). However, with the development of technology, many improvements of the method flow today can be regarded as direct improvements of the hardware circuit structure. Designers almost always obtain the corresponding hardware circuit structure by programming the improved method flow into the hardware circuit. Therefore, it cannot be said that the improvement of a method flow cannot be implemented with hardware entity modules. For example, a programmable logic device (PLD) (such as a field programmable gate array (FPGA)) is such an integrated circuit whose logical function is determined by the user's programming of the device. Designers can "integrate" a digital system on a PLD by programming themselves, without having to ask chip manufacturers to design and make dedicated integrated circuit chips. Moreover, nowadays, instead of manually making integrated circuit chips, this kind of programming is mostly implemented by "logic compiler" software, which is similar to the software compiler used when developing and writing programs, and the original code before compilation must also be written in a specific programming language, which is called hardware description language (HDL). There is not only one kind of HDL, but many kinds, such as ABEL (Advanced Boolean Expression Language), AHDL (Altera Hardware Description Language), Confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java Hardware Description Language), Lava, Lola, MyHDL, PALASM, RHDL (Ruby Hardware Description Language), etc. The most commonly used ones are VHDL (Very-High-Speed Integrated Circuit Hardware Description Language) and Verilog. Those skilled in the art should also know that it is only necessary to program the method flow slightly in the above-mentioned hardware description languages and program it into the integrated circuit, and then it is easy to obtain the hardware circuit that implements the logic method flow.

控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20 以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。The controller may be implemented in any suitable manner, for example, the controller may take the form of a microprocessor or processor and a computer-readable medium storing a computer-readable program code (e.g., software or firmware) executable by the (micro)processor, a logic gate, a switch, an application-specific integrated circuit (ASIC), a programmable logic controller, and an embedded microcontroller. Examples of controllers include, but are not limited to, the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320. The memory controller may also be implemented as part of the control logic of the memory. It is also known to those skilled in the art that, in addition to implementing the controller in a purely computer-readable program code manner, the controller may be implemented in the form of a logic gate, a switch, an application-specific integrated circuit, a programmable logic controller, and an embedded microcontroller by logically programming the method steps. Therefore, such a controller may be considered as a hardware component, and the devices for implementing various functions included therein may also be considered as structures within the hardware component. Or even, the devices for implementing various functions may be considered as both software modules for implementing the method and structures within the hardware component.

上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。The systems, devices, modules or units described in the above embodiments may be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is a computer. Specifically, the computer may be, for example, a personal computer, a laptop computer, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本说明书时可以把各单元的功能在同一个或多个软件和/或硬件中实现。For the convenience of description, the above device is described in various units according to their functions. Of course, when implementing this specification, the functions of each unit can be implemented in the same or multiple software and/or hardware.

本领域内的技术人员应明白,本说明书的实施例可提供为方法、系统、或计算机程序产品。因此,本说明书可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本说明书可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems, or computer program products. Therefore, this specification may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

本说明书是参照根据本说明书实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。This specification is described with reference to the flowcharts and/or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of this specification. It should be understood that each process and/or box in the flowchart and/or block diagram, as well as the combination of the processes and/or boxes in the flowchart and/or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。In a typical configuration, a computing device includes one or more processors (CPU), input/output interfaces, network interfaces, and memory.

内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and/or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本说明书可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems or computer program products. Therefore, this specification may take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware. Furthermore, this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。This specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. This specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

以上所述仅为本说明书的实施例而已,并不用于限制本说明书。对于本领域技术人员来说,本说明书可以有各种更改和变化。凡在本说明书的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本说明书的权利要求范围之内。The above description is only an embodiment of the present specification and is not intended to limit the present specification. For those skilled in the art, the present specification may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present specification shall be included in the scope of the claims of the present specification.

Claims (9)

1. A method of sound field detection based on code phase imaging, characterized in that it is implemented by a sound field detection system based on code phase imaging, said system comprising:
The parallel light field transmitting module (1) is used for generating a pulse parallel light field to illuminate a sound field medium, the pulse parallel light field carries out wave front modulation through the sound field medium (2), imaging is carried out through the imaging module (4), modulation is carried out through the encoding plate (5), and light field intensity distribution is obtained through the light field detector (6);
The system further comprises:
An ultrasonic transducer module (3) for generating a sound field of high frequency to be detected;
the control and calculation module (7) is used for controlling the synchronization of the parallel light field emission module (1) and the light field detector (6), and storing the light field intensity distribution acquired by the light field detector (6) so as to calculate the complex amplitude distribution of the high-frequency sound field to be detected;
the method comprises the following steps:
The ultrasonic transducer module (3) generates a sound field of high frequency to be detected;
acquiring light field intensity distribution recorded by a light field detector (6);
Setting the plane of the light field detector (6) as a first plane, and setting the plane of the coding plate (5) as a second plane; setting an initial light field complex amplitude distribution at the front surface of the code plate (5); the initial optical field complex amplitude distribution is iteratively propagated among the first plane, the second plane and the constraint plane; calculating errors based on the light field intensity distribution, and finishing iteration when the errors are smaller than a threshold value to obtain wavefront distribution at the front surface of the coding plate (5);
and transmitting the wave front distribution at the front surface of the coding plate (5) back to the plane where the ultrasonic transducer module (3) is located, so as to calculate wave front amplitude and phase distribution of the wave front at the front surface of the coding plate (5) at the ultrasonic transducer module (3) and obtain complex amplitude distribution of the sound field.
2. A sound field detection method based on code phase imaging according to claim 1, characterized in that the process of acquiring the wavefront distribution at the front surface of the code plate (5) comprises:
setting the initial light field complex amplitude distribution at the front surface of the code plate (5) as Wherein A andThe amplitude and phase distribution of random guess are respectively, n is the corresponding time sequence number, and t n is the recording time of the corresponding light field detector (6);
Taking the initial light field complex amplitude distribution at the front surface of the code plate (5) as a1 st iteration value E 1(x,y,tn)=E0(x,y,tn);
after the light field complex amplitude distribution is modulated by the encoding plate (5), a first wavefront is obtained, and the expression is as follows:
ψk(x,y,tn)=Ek(x,y,tn)·T(x,y);
wherein k is the iteration number, k=1, 2,3 …, and T (x, y) is the transmittance complex amplitude distribution of the coding plate (5);
when the first wavefront propagates to the first plane in the free space, the first wavefront is constrained by using the light field intensity distribution to obtain a second wavefront, and the expression is as follows:
Wherein L 1 denotes a distance between the second plane and the first plane, Representing the light field free space propagation L distance; angle (G (x, y, t n)) represents the phase angle taken for the wavefront;
And (3) transmitting the second wavefront back to the second plane, and updating by using the transmittance complex amplitude distribution to obtain a third wavefront at the front surface of the coding plate (5), wherein the expression is as follows:
E'(x,y,tn)=E(x,y,tn)+T*(x,y)/max(T*(x,y)·T(x,y))·[ψ'(x,y,tn)-ψ(x,y,tn)];
transmitting the third wavefront to a constraint surface, wherein the constraint surface is a spectrum surface or a focus surface of the third wavefront; constructing an aperture constraint function, and constraining time and wavelength by using the aperture constraint function to obtain a fourth wavefront, wherein the expression is as follows:
wherein L 2 represents the distance between the second plane and the constraint surface, and H (x, y) is the aperture function constraint;
transmitting the fourth wavefront to the second plane to obtain a fifth wavefront, where the expression is as follows:
And iterating according to the process, calculating an error based on the second wavefront and the pulse diffraction light spot intensity matrix, finishing iteration when the error is smaller than a threshold value, and taking a fifth wavefront obtained by iteration as wavefront distribution at the front surface of the coding plate (5).
3. The method of claim 2, wherein calculating the error based on the light field intensity distribution comprises:
The expression of the error is as follows:
Wherein I (x, y, t n) represents the light field intensity distribution recorded by the light field detector (6); g (x, y, t n) represents the second wavefront.
4. A sound field detection method based on coded phase imaging according to claim 1, characterized in that the step of transmitting back the wave front distribution at the front surface of the code plate (5) to the plane in which the ultrasound transducer module (3) is located, thereby calculating the wave front amplitude and phase distribution of the wave front distribution at the front surface of the code plate (5) at the ultrasound transducer module (3), resulting in a complex amplitude distribution of the sound field comprises:
Acquiring a wavefront distribution at the front surface of the code plate (5), denoted as E (x, y, t n);
the wave front distribution at the front surface of the coding plate (5) is transmitted back to the plane where the ultrasonic transducer module (3) is located, and the wave front distribution at the plane where the ultrasonic transducer module (3) is located is obtained, and the expression is as follows:
Wherein L 3 is the distance between the second plane and the plane of the ultrasonic transducer module (3);
Calculating a wavefront amplitude abs (U (x, y, t n)) and a phase distribution angle (U (x, y, t n)) that distribute a wavefront at the front surface of the encoding plate (5) at the ultrasound transducer module (3); based on the density of the sound field medium (2), a complex amplitude distribution of the sound field is obtained.
5. The encoding phase imaging based sound field detection method of claim 1, wherein the system further comprises: a two-dimensional translation module (8) for moving the ultrasound transducer module (3) in x-direction and y-direction.
6. The sound field detection method based on code phase imaging according to claim 1, wherein the parallel light field emission module (1) selects a space laser to generate a parallel light field system or a fiber laser to generate a parallel light field system; the pulse width of the pulse parallel light field generated by the parallel light field transmitting module (1) is smaller than 10ns, and the repetition frequency is smaller than 1kHz.
7. The sound field detection method based on code phase imaging according to claim 1, wherein the imaging module (4) adopts a beam-shrinking imaging system for low-frequency large-caliber sound field detection; or the imaging module (4) adopts an imaging amplification system for detecting the high-frequency small-caliber sound field.
8. An electronic device comprising a memory and a processor, wherein the memory is coupled to the processor; wherein the memory is configured to store program data and the processor is configured to execute the program data to implement the encoding phase imaging based sound field detection method of any one of claims 1 to 7.
9. A computer-readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the encoding phase imaging based sound field detection method according to any one of claims 1 to 7.
CN202410364077.0A 2024-03-28 2024-03-28 Acoustic field detection method, system, device and medium based on coded phase imaging Active CN117969677B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410364077.0A CN117969677B (en) 2024-03-28 2024-03-28 Acoustic field detection method, system, device and medium based on coded phase imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410364077.0A CN117969677B (en) 2024-03-28 2024-03-28 Acoustic field detection method, system, device and medium based on coded phase imaging

Publications (2)

Publication Number Publication Date
CN117969677A CN117969677A (en) 2024-05-03
CN117969677B true CN117969677B (en) 2024-07-26

Family

ID=90846369

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410364077.0A Active CN117969677B (en) 2024-03-28 2024-03-28 Acoustic field detection method, system, device and medium based on coded phase imaging

Country Status (1)

Country Link
CN (1) CN117969677B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118655090B (en) * 2024-08-16 2025-02-07 之江实验室 A method, device, storage medium and electronic device for photoacoustic imaging demixing
CN118795034B (en) * 2024-09-13 2025-02-14 之江实验室 Sound field detection system, method, electronic device, and medium based on dual-detector phase imaging

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112683794A (en) * 2020-12-11 2021-04-20 中国科学院上海光学精密机械研究所 Phase imaging and element detection device and method based on wavefront modulation

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL137447A (en) * 2000-07-23 2007-03-08 Israel Atomic Energy Comm Apparatus and method for probing light absorbing agents in biological tissues
WO2005069997A2 (en) * 2004-01-20 2005-08-04 Northeastern University Enhanced detection of acousto-photonic emissions in optically turbid media using a photo-refractive crystal-based detection system
US9354166B2 (en) * 2006-10-06 2016-05-31 California Institute Of Technology Time-reversal of variance-encoded light (TROVE)
US10203274B2 (en) * 2006-10-06 2019-02-12 California Institute Of Technology Optical focusing inside scattering media with time-reversed ultrasound microbubble encoded (TRUME) light
EP2110703B1 (en) * 2008-04-15 2010-08-18 JDS Uniphase Corporation Retarder-based despeckle device and method for laser illumination systems
US8930145B2 (en) * 2010-07-28 2015-01-06 Covidien Lp Light focusing continuous wave photoacoustic spectroscopy and its applications to patient monitoring
US9679360B2 (en) * 2013-05-10 2017-06-13 Trustees Of Princeton University High-resolution light-field imaging
CN106018286B (en) * 2016-05-18 2018-12-11 哈尔滨工业大学 Reconstruction detection device and method for reconstructing based on the distribution of the dispersive medium optical parameter of light-field camera and laser of frequency modulation
IL251636B (en) * 2017-04-06 2018-02-28 Yoav Berlatzky Coherence camera system and method thereof
CN107913056B (en) * 2017-11-07 2019-11-15 华南理工大学 A method for realizing arbitrary ultrasonic fields
US11030776B2 (en) * 2019-02-01 2021-06-08 Molecular Devices (Austria) GmbH Calibration of a light-field imaging system
CN114111860B (en) * 2021-12-03 2022-08-30 北京科技大学 Distributed phi-OTDR sensing method and system based on multi-frequency pulse coding
CN114486812B (en) * 2022-01-11 2024-07-23 清华大学 A complex amplitude imaging method and system
CN114544776B (en) * 2022-02-15 2025-06-06 南京华秦光声科技有限责任公司 Super-resolution ultrasonic holographic ghost imaging device and imaging method based on photoacoustic waves
CN115685565B (en) * 2022-08-05 2024-07-26 华侨大学 A method and device for controlling a complex light field passing through a scattering medium
CN116087111B (en) * 2023-03-29 2023-06-30 之江实验室 A photoacoustic microscopy imaging system and method based on arrayed frequency field modulated light beams
CN116399244A (en) * 2023-05-08 2023-07-07 浙江大学 High-resolution surface measurement method and device based on broad-spectrum laser and wavefront coding
CN116763257A (en) * 2023-06-25 2023-09-19 南京理工大学 A multi-pixel joint encoding DMD high-speed skin light field control method
CN117539068A (en) * 2023-11-20 2024-02-09 合肥工业大学 Devices and methods for controlling multi-degree-of-freedom information in three-dimensional light fields through scattering media

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112683794A (en) * 2020-12-11 2021-04-20 中国科学院上海光学精密机械研究所 Phase imaging and element detection device and method based on wavefront modulation

Also Published As

Publication number Publication date
CN117969677A (en) 2024-05-03

Similar Documents

Publication Publication Date Title
CN117969677B (en) Acoustic field detection method, system, device and medium based on coded phase imaging
US9947118B2 (en) Non-interferometric phase measurement
JP5364233B2 (en) Electromagnetic field simulator and electromagnetic field simulation program
CN110095968B (en) Off-axis digital holographic image reconstruction device and method and microscopic imaging system
JP5066307B1 (en) Optical microphone
US11412118B2 (en) Electromagnetic wave phase/amplitude generation device, electromagnetic wave phase/amplitude generation method, and electromagnetic wave phase/amplitude generation program
CN113970367B (en) Free field laser absolute calibration method of interference type optical fiber hydrophone
JP2009139134A (en) Large area invariant light spot imaging apparatus and method
CN109357615A (en) Common device for speckle interference and shear speckle interference
CN105352915B (en) A kind of dynamic measurement method of refractive index Two dimensional Distribution
Karatodorov et al. Multi-pass probing for high-sensitivity tomographic interferometry
CN106950811A (en) A kind of digital composite holographic imaging method and device
TW202232475A (en) Optical measurement system and optical measurement method
CN115824400A (en) Laser power detection system, welding system and detection method for laser welding workpiece
CN115494005A (en) Semiconductor defect detection device and method based on infrared micro-digital holography
CN112741599B (en) A Time-Reversal Ultrasonic Encoded Light Focusing Method Guided by Photoacoustic Signal Inversion
CN117968866B (en) A single ultrashort pulse spatiotemporal coupling measurement method, system, device, and medium
CN118795034B (en) Sound field detection system, method, electronic device, and medium based on dual-detector phase imaging
CN115836209A (en) Observation device and observation method
KR101995122B1 (en) Method and apparatus for recording and projecting 3 dimensional holographic image using scattering layer
JP7246543B1 (en) Communication device, program, and communication method
CN116718126A (en) A quality inspection system and method for gene chip assembly
CN108508094B (en) A test device and test method for sound absorption capacity of materials based on schlieren method
CN114723616A (en) A kind of restoration method, device and equipment of blurred image
JP2018036561A (en) Hologram reproduction device and program

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
GR01 Patent grant
GR01 Patent grant