[go: up one dir, main page]

CN113842158A - Photoacoustic/ultrasonic endoscopic probe based on fixed-focus sound field and dynamic focus reconstruction algorithm - Google Patents

Photoacoustic/ultrasonic endoscopic probe based on fixed-focus sound field and dynamic focus reconstruction algorithm Download PDF

Info

Publication number
CN113842158A
CN113842158A CN202110910237.3A CN202110910237A CN113842158A CN 113842158 A CN113842158 A CN 113842158A CN 202110910237 A CN202110910237 A CN 202110910237A CN 113842158 A CN113842158 A CN 113842158A
Authority
CN
China
Prior art keywords
transducer
focus
sound field
photoacoustic
ultrasonic
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
Application number
CN202110910237.3A
Other languages
Chinese (zh)
Other versions
CN113842158B (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.)
SHENZHEN RESEARCH INSTITUTE CENTRAL SOUTH UNIVERSITY
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN202110910237.3A priority Critical patent/CN113842158B/en
Publication of CN113842158A publication Critical patent/CN113842158A/en
Application granted granted Critical
Publication of CN113842158B publication Critical patent/CN113842158B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
    • A61B5/0097Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying acoustic waves and detecting light, i.e. acousto-optic measurements
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Acoustics & Sound (AREA)
  • Gynecology & Obstetrics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

The invention relates to the field of endoscopic imaging, and discloses a photoacoustic/ultrasonic endoscopic probe based on a fixed-focus sound field and a dynamic focus reconstruction algorithm, which comprises the following steps of 1) determining an equivalent focus sound field and an equivalent focus of a transducer according to inherent parameters of the transducer; 2) determining the length of a path, which is F, of a signal fed back by the equivalent focus received by the transducer by taking the equivalent focus as a reference point, and determining the equivalent beam waist radius omega of a sound field focused by the transducer0And equivalent half depth of focus Z0(ii) a 3) And calculating the pixel value of a single pixel point obtained by scanning the transducer according to the parameters, and processing the pixel point by a delay superposition algorithm to obtain a final image. The invention optimizes the conventional delay superposition algorithm to obtain a novel delay superposition algorithm, and adopts a novel endoscopic probe which is combined with the algorithm for use, so that under the premise of ensuring the imaging speed,the imaging definition and accuracy of the final imaging are improved.

Description

基于定焦声场的光声/超声内窥探头及动态聚焦重建算法Photoacoustic/Ultrasonic Endoscopic Probe and Dynamic Focus Reconstruction Algorithm Based on Fixed-focus Sound Field

技术领域technical field

本发明涉及内窥成像领域,具体地涉及一种基于定焦声场的光声/超声内窥探头及动态聚焦重建算法。The invention relates to the field of endoscopic imaging, in particular to a photoacoustic/ultrasonic endoscopic probe based on a fixed-focus sound field and a dynamic focus reconstruction algorithm.

背景技术Background technique

声聚焦光声内窥成像是一种新兴的集高分辨,高穿透深度,高光学对比度为一体的生物医学内窥成像手段。在很多声聚焦光声内窥成像系统中,由于可以采用基于压电效应的换能器来接收超声,而这个换能器同时也可以发射超声,所以光声/超声成像由于同用一个换能器所以可以方便的集成在一起,形成声聚焦光声/超声成像系统。光声/超声成像在特定的场合中有很好的应用前景。Acoustic focused photoacoustic endoscopic imaging is an emerging biomedical endoscopic imaging method that integrates high resolution, high penetration depth, and high optical contrast. In many acoustic focusing photoacoustic endoscopic imaging systems, since a piezoelectric effect-based transducer can be used to receive ultrasound, and this transducer can also transmit ultrasound, the photoacoustic/ultrasonic imaging uses the same transducer. Therefore, the devices can be easily integrated together to form an acoustically focused photoacoustic/ultrasonic imaging system. Photoacoustic/ultrasound imaging has good application prospects in specific occasions.

在常规的光声/超声内窥成像探头中,其基本的结构可使其进行360度环形扫描成像。另外,也可以采用平面超声换能器,然后利用抛物面的超声及激光反射镜来形成聚焦声场。在此类系统成像时,光声和超声内窥成像算法都比较类似,其一般采用类似于B超中常用的超声束直线传播模型,直接按照像素相对于换能器的距离和声速计算得到的延迟时间,以及换能器内窥扫描时的旋转角度,来将所得数据进行扇形排列组合成像。其中,由于超声涉及到信号的往返,因此在超声成像时声速需要设置成实际声速的一半。In the conventional photoacoustic/ultrasonic endoscopic imaging probe, its basic structure enables it to perform 360-degree circular scanning imaging. In addition, a planar ultrasonic transducer can also be used, and then a parabolic ultrasonic and laser mirror can be used to form a focused sound field. In the imaging of such systems, the photoacoustic and endoscopic ultrasound imaging algorithms are relatively similar. Generally, they use a linear propagation model similar to the ultrasonic beam commonly used in B-mode ultrasound, and are calculated directly according to the distance of the pixel relative to the transducer and the speed of sound. The delay time, and the rotation angle of the transducer during endoscopic scanning are used to fan out the resulting data for imaging. Among them, since ultrasound involves the round-trip of the signal, the sound speed needs to be set to half of the actual sound speed during ultrasound imaging.

此种成像方法简单直接,但是由于定焦探头的应用,其只能在焦点位置得到最佳横向分辨率,包括切向上的以及沿着探头的轴向上的最佳分辨率,而在远离焦点区域,其横向分辨率将迅速降低。当超声传感器的数值孔径越大,其最佳横向分辨率最小,但是此种情况下其成像景深也越小,反之亦然。这样的现状,主要是由于当前的声聚焦光声/超声内窥成像算法没有考虑到实际换能器的聚焦声场造成的,因此迫切的需要一种新的算法来考虑换能器的声场,使得离焦区域也能有较高的横向分辨率,进而保证整体图像成像的清晰准确。This imaging method is simple and straightforward, but due to the application of the fixed-focus probe, it can only obtain the best lateral resolution at the focal position, including the best resolution in the tangential direction and along the axial direction of the probe, and the best resolution in the far-focus position. area, its lateral resolution will decrease rapidly. When the numerical aperture of the ultrasonic sensor is larger, its optimal lateral resolution is the smallest, but the imaging depth of field is also smaller in this case, and vice versa. This situation is mainly caused by the fact that the current acoustic focusing photoacoustic/ultrasonic endoscopic imaging algorithm does not consider the focused sound field of the actual transducer, so a new algorithm is urgently needed to consider the sound field of the transducer, so that The out-of-focus area can also have a higher lateral resolution, thereby ensuring the clarity and accuracy of the overall image.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的问题是提供一种基于定焦声场的光声/超声内窥动态聚焦重建算法,该方法能够提高离焦区域的横向分辨率,使得成像更加清晰准确。The problem to be solved by the present invention is to provide a photoacoustic/ultrasonic endoscopic dynamic focusing reconstruction algorithm based on a fixed-focus sound field, which can improve the lateral resolution of the out-of-focus region and make the imaging clearer and more accurate.

本发明进一步要解决的问题是提供一种基于定焦声场的光声/超声内窥探头,该探头能够提高离焦区域的横向分辨率,使得成像更加清晰准确。The further problem to be solved by the present invention is to provide a photoacoustic/ultrasonic endoscopic probe based on a fixed-focus sound field, which can improve the lateral resolution of the out-of-focus region and make the imaging clearer and more accurate.

为了解决上述技术问题,本发明提供了一种基于定焦声场的光声/超声内窥探头及动态聚焦重建算法,包括以下步骤:1)根据换能器的固有参数,确定所述换能器的等效聚焦声场和等效焦点;2)以等效焦点为基准点,确定所述换能器接收到的等效焦点反馈的信号所经过的路径长度为F,并确定出所述换能器聚集声场的等效束腰半径ω0和等效半焦深Z0;3)根据上述参数,计算出所述换能器扫描得到的单个像素点的像素值,并将原始的延时叠加算法进行优化,得到新型的延时叠加算法,通过该新型延时叠加算法,将扫描得到的图像像素点进行处理,以所述换能器接收到的该像素点所对应的超声信号的时间长短进行排序,从而获得最终的图像。In order to solve the above technical problems, the present invention provides a photoacoustic/ultrasonic endoscope probe based on a fixed-focus sound field and a dynamic focus reconstruction algorithm, including the following steps: 1) According to the inherent parameters of the transducer, determine the transducer 2) Taking the equivalent focus as the reference point, determine the length of the path traversed by the signal fed back by the equivalent focus received by the transducer as F, and determine the transducer The equivalent beam waist radius ω 0 and the equivalent semi-focal depth Z 0 of the focused sound field of the transducer; 3) According to the above parameters, calculate the pixel value of a single pixel point scanned by the transducer, and superimpose the original time delay The algorithm is optimized to obtain a new type of delay stacking algorithm. Through this new type of delay stacking algorithm, the scanned image pixels are processed, and the time length of the ultrasonic signal corresponding to the pixel received by the transducer is determined. sort to get the final image.

优选地,在步骤3中,通过如下公式得到最终的图像:Preferably, in step 3, the final image is obtained by the following formula:

Figure BDA0003203429890000021
Figure BDA0003203429890000021

其中,Nd为换能器扫描过的方向位置数,

Figure BDA0003203429890000022
为需要重建的像素坐标,S(i,t)为换能器在第i个位置时换能器接收到的信号,t为时间,
Figure BDA0003203429890000023
为第i个像素点所在位置的相位系数,
Figure BDA0003203429890000024
为第i个像素点所在位置的权重系数。通过该优选技术方案,将传统的延时叠加算法进行优化,得到新的延时叠加算法,延时叠加算法能够通过简单的计算推导,从而得到最终的图像数据,并且通过该叠加算法的计算,所需要的计算时间更短,效率更高,并且该算法,在保证了计算效率的同时,提高了成像的清晰度,使得最终的成像更加清晰准确。Among them, Nd is the number of directional positions scanned by the transducer,
Figure BDA0003203429890000022
is the pixel coordinate to be reconstructed, S(i,t) is the signal received by the transducer when the transducer is at the i-th position, t is the time,
Figure BDA0003203429890000023
is the phase coefficient at the position of the i-th pixel,
Figure BDA0003203429890000024
is the weight coefficient of the position of the i-th pixel. Through this preferred technical solution, the traditional time delay superposition algorithm is optimized to obtain a new time delay superposition algorithm. The time delay superposition algorithm can be deduced through simple calculation, so as to obtain the final image data, and through the calculation of the superposition algorithm, The required calculation time is shorter and the efficiency is higher, and the algorithm, while ensuring the calculation efficiency, improves the clarity of the imaging, making the final imaging clearer and more accurate.

进一步优选地,所述相位系数由像素点与换能器焦点的相对位置所决定,在光声成像中,相位系数为:Further preferably, the phase coefficient is determined by the relative position of the pixel point and the focus of the transducer. In photoacoustic imaging, the phase coefficient is:

Figure BDA0003203429890000031
Figure BDA0003203429890000031

在超声成像中,相位系数为:In ultrasound imaging, the phase coefficient is:

Figure BDA0003203429890000032
Figure BDA0003203429890000032

其中,v为超声波在介质中的传播速度,a为像素点相对于所述换能器焦点的轴向距离,b为像素点相对于所述换能器焦点的径向距离。通过该优选技术方案,通过对光声成像和超声成像的相位系数进行限定,使得该设备在进行光声成像和超声成像时能够更加准确的定位像素点与换能器焦点之间的相对位置关系,从而提高成像的准确性。Wherein, v is the propagation velocity of ultrasonic waves in the medium, a is the axial distance of the pixel point relative to the transducer focus, and b is the radial distance of the pixel point relative to the transducer focus. Through this preferred technical solution, by defining the phase coefficients of photoacoustic imaging and ultrasonic imaging, the device can more accurately locate the relative positional relationship between the pixel point and the transducer focus when performing photoacoustic imaging and ultrasonic imaging. , thereby improving the accuracy of imaging.

优选地,所述像素点距离换能器焦点的轴向距离a的表达式如下:Preferably, the expression of the axial distance a between the pixel point and the focus of the transducer is as follows:

Figure BDA0003203429890000033
Figure BDA0003203429890000033

其中,

Figure BDA0003203429890000034
为换能器在第i个位置时的焦点,
Figure BDA0003203429890000035
为换能器第i个位置时的旋转中心。通过该优选技术方案,确定了像素点与换能器焦点之间的轴向距离,从而能够提高图像像素点在轴向距离的方向上成像的准确性。in,
Figure BDA0003203429890000034
is the focus of the transducer at the i-th position,
Figure BDA0003203429890000035
is the rotation center of the transducer at the i-th position. With this preferred technical solution, the axial distance between the pixel point and the transducer focus is determined, so that the imaging accuracy of the image pixel point in the direction of the axial distance can be improved.

进一步优选地,所述像素点距离换能器焦点的径向距离b的表达式如下:Further preferably, the expression of the radial distance b of the pixel point from the transducer focus is as follows:

Figure BDA0003203429890000041
Figure BDA0003203429890000041

其中,

Figure BDA0003203429890000042
为换能器在第i个位置时的焦点,
Figure BDA0003203429890000043
为换能器第i个位置时的旋转中心。通过该优选技术方案,确定了像素点与换能器焦点之间的径向距离,从而能够提高图像像素点在径向距离的方向上成像的准确性。in,
Figure BDA0003203429890000042
is the focus of the transducer at the i-th position,
Figure BDA0003203429890000043
is the rotation center of the transducer at the i-th position. Through this preferred technical solution, the radial distance between the pixel point and the focus of the transducer is determined, so that the imaging accuracy of the image pixel point in the direction of the radial distance can be improved.

优选地,所述

Figure BDA0003203429890000044
为第i个像素所在位置的权重系数,所述权重系数可通过如下的公式计算得出:Preferably, the
Figure BDA0003203429890000044
is the weight coefficient of the position of the i-th pixel, and the weight coefficient can be calculated by the following formula:

Figure BDA0003203429890000045
Figure BDA0003203429890000045

其中,

Figure BDA0003203429890000046
为一个由当前像素点与换能器的距离所决定的位置限制函数。通过该优选技术方案,将权重系数通过位置限制函数进行限制,能够进一步对像素点的位置进行限定,从而能够更加准确的确定像素点的相对位置,进而保证了叠加后的图像的准确性。in,
Figure BDA0003203429890000046
is a position limit function determined by the distance between the current pixel and the transducer. Through the preferred technical solution, the weight coefficient is limited by the position limit function, and the position of the pixel point can be further limited, so that the relative position of the pixel point can be more accurately determined, thereby ensuring the accuracy of the superimposed image.

进一步优选地,所述位置限制函数通过以下方式来限制延时叠加的范围:Further preferably, the position limit function limits the range of the delay stacking in the following manner:

Figure BDA0003203429890000047
Figure BDA0003203429890000047

其中,a1、a2、b1和b2为位置限制函数的限制参数,所述限制参数由以下公式计算得出:Wherein, a 1 , a 2 , b 1 and b 2 are the limit parameters of the position limit function, and the limit parameters are calculated by the following formula:

Figure BDA0003203429890000048
Figure BDA0003203429890000048

a2=R,b1=-R,b2=R,a 2 =R, b 1 =-R, b 2 =R,

其中,D为内窥探头的直径,R为成像区域的半径。通过该优选技术方案,引入限制参数来进一步限定延时叠加的范围,并且将限制参数与探头的固有参数以及最终的成像参数相结合,从而提高了该限制函数的实用性,使得该限制函数能够准确的对延时叠加的范围进行限定,从而提高了经过延时叠加后图像的准确性。Among them, D is the diameter of the endoscopic probe, and R is the radius of the imaging area. Through this preferred technical solution, a limit parameter is introduced to further limit the range of time-lapse stacking, and the limit parameter is combined with the inherent parameters of the probe and the final imaging parameters, thereby improving the practicability of the limit function, making the limit function capable of The range of the time-lapse stacking is accurately limited, thereby improving the accuracy of the image after the time-lapse stacking.

本发明第二方面提供了一种基于定焦声场的光声/超声内窥探头,包括:透明外壳、氟化钙镜片、反射镜、保护结构和换能器,所述内窥探头由旋转电机驱动其运动,以完成侧向的环形扫描成像作业,所述透明外壳包裹在所述内窥探头的外表面,所述换能器发射的声波经过氟化钙镜片和反射镜的折射和聚焦,照射到探测物上,所述内窥探头能够实现本发明第一方面所提供的一种基于定焦声场的光声/超声内窥动态聚焦重建算法。A second aspect of the present invention provides a photoacoustic/ultrasonic endoscopic probe based on a fixed-focus sound field, comprising: a transparent casing, a calcium fluoride lens, a reflector, a protective structure and a transducer, the endoscopic probe is composed of a rotating motor Drive its movement to complete the lateral annular scanning imaging operation, the transparent shell wraps the outer surface of the endoscope probe, the sound wave emitted by the transducer is refracted and focused by the calcium fluoride mirror and the reflector, When irradiated on the probe, the endoscopic probe can implement the photoacoustic/ultrasonic endoscopic dynamic focus reconstruction algorithm based on the fixed-focus sound field provided in the first aspect of the present invention.

优选地,所述换能器采用平场超声换能器。通过该优选技术方案,采用平场超声换能器,其声场更为均匀,且制作难度相对于聚焦超声换能器的制作难度更低,并且接收面为平面便于接收从各个位置反馈的声波信号,从而提高最终图像成像的清晰度。Preferably, the transducer is a flat-field ultrasonic transducer. Through this preferred technical solution, a flat-field ultrasonic transducer is used, and its sound field is more uniform, and its manufacturing difficulty is lower than that of a focused ultrasonic transducer, and the receiving surface is flat to facilitate the reception of sound wave signals fed back from various positions , thereby improving the clarity of the final image imaging.

进一步优选地,所述反射镜可采用抛物面反射镜,且所述抛物面反射镜以弧高为基准线,所述基准线与所述氟化钙镜片的中心相垂直。通过该优选技术方案,通过抛物面反射镜,使得换能器与聚焦点相对于抛物面反射镜的中心呈镜像对称,并且该反射镜能够实现光声和超声的反射,从而使得该设备普遍适用于各种光声和超声的检查。Further preferably, the reflecting mirror can be a parabolic reflecting mirror, and the parabolic reflecting mirror takes the arc height as a reference line, and the reference line is perpendicular to the center of the calcium fluoride mirror. Through this preferred technical solution, the parabolic reflector is used to make the transducer and the focal point mirror-symmetrical with respect to the center of the parabolic reflector, and the reflector can realize photoacoustic and ultrasonic reflection, so that the device is generally applicable to various photoacoustic and ultrasound examinations.

通过上述技术方案,本发明提供的基于定焦声场的光声/超声内窥动态聚焦重建算法,将常规的延时叠加算法进行优化,对其中的参数进行优化限制,并引入的限制参数,从而得到一种新型的延时叠加算法,通过该新型延时叠加算法,能够在保证成像速度的前提下,提高最终成像的成像清晰度和准确性。Through the above technical solutions, the photoacoustic/ultrasonic endoscopic dynamic focusing reconstruction algorithm based on fixed-focus sound field provided by the present invention optimizes the conventional time delay stacking algorithm, optimizes and limits the parameters therein, and introduces the limiting parameters, thereby A new type of time-lapse stacking algorithm is obtained, through which the imaging clarity and accuracy of the final imaging can be improved on the premise of ensuring the imaging speed.

本发明还提供了基于定焦声场的光声/超声内窥探头,通过采用平面超声换能器与抛物面反射镜相结合进行工作,提高了该设备普遍的适用性,并且该设备能够接收从各个角度反馈回来的声波信号,从而增加了其最终成像的清晰度和准确性。The present invention also provides a photoacoustic/ultrasonic endoscope probe based on a fixed-focus sound field. By using a plane ultrasonic transducer in combination with a parabolic mirror to work, the general applicability of the device is improved, and the device can receive signals from various Angle feedback sound wave signal, which increases the clarity and accuracy of its final image.

本发明的其他特征和优点将在随后的具体实施方式部分予以详细说明。Other features and advantages of the present invention will be described in detail in the detailed description that follows.

附图说明Description of drawings

图1为现有技术中所使用的声聚焦光声内窥探头结构图以及扫描等效结构图;1 is a structural diagram of an acoustic focusing photoacoustic endoscope probe used in the prior art and a scanning equivalent structural diagram;

图2为本发明所使用的基于抛物面反射镜的声聚焦光声内窥探头结构图以及扫描等效结构图;Fig. 2 is the structure diagram and scanning equivalent structure diagram of the acoustic focusing photoacoustic endoscope probe based on the parabolic reflector used in the present invention;

图3为本发明所使用的反射镜结构示意图;Fig. 3 is the mirror structure schematic diagram used in the present invention;

图4为本发明所扫描的物体的结构图;Fig. 4 is the structure diagram of the object scanned by the present invention;

图5(a)和(b)分别为本发明的算法与原始算法对比的光声结果图;Fig. 5 (a) and (b) are respectively the photoacoustic result graph that the algorithm of the present invention is compared with the original algorithm;

图6(a)和(b)分别为本发明的算法与原始算法对比的超声结果图。Figures 6(a) and (b) are respectively the ultrasound result graphs comparing the algorithm of the present invention with the original algorithm.

附图标记reference number

1 透明外壳 2 中空声聚焦超声换能器1 Transparent housing 2 Hollow acoustically focused ultrasound transducer

3 多模光纤 4 顶部45度反射镜3 Multimode fiber 4 45 degree reflector on top

5 氟化钙镜片 6 反射镜5 calcium fluoride lens 6 reflector

7 保护结构 8 平场超声换能器7 Protective structure 8 Flat-field ultrasonic transducer

具体实施方式Detailed ways

以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and explain the present invention, but not to limit the present invention.

在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“设置”、“连接”应做广义理解,例如,术语“连接”可以是固定连接,也可以是可拆卸连接,或者是一体连接;可以是直接连接,也可以是通过中间媒介间接连接,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。In the description of the present invention, it should be noted that, unless otherwise expressly specified and limited, the terms "installation", "arrangement" and "connection" should be understood in a broad sense. For example, the term "connection" may be a fixed connection or a It can be a detachable connection or an integral connection; it can be a direct connection or an indirect connection through an intermediate medium, and it can be the internal communication of two elements or an interaction relationship between the two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific situations.

本发明的基于定焦声场的光声/超声内窥探头及动态聚焦重建算法的一种实施方式,包括以下步骤:1)根据换能器的固有参数,确定所述换能器的等效聚焦声场和等效焦点;2)以等效焦点为基准点,确定所述换能器接收到的等效焦点反馈的信号所经过的路径长度为F,并确定出所述换能器聚集声场的等效束腰半径ω0和等效半焦深Z0;3)根据上述参数,计算出所述换能器扫描得到的单个像素点的像素值,并将原始的延时叠加算法进行优化,得到新型的延时叠加算法,通过该新型延时叠加算法,将扫描得到的图像像素点进行处理,以所述换能器接收到的该像素点所对应的超声信号的时间长短进行排序,从而获得最终的图像。An embodiment of the photoacoustic/ultrasonic endoscope probe and dynamic focus reconstruction algorithm based on the fixed-focus sound field of the present invention includes the following steps: 1) According to the inherent parameters of the transducer, determine the equivalent focus of the transducer Sound field and equivalent focus; 2) Taking the equivalent focus as the reference point, determine the path length of the signal fed back by the equivalent focus received by the transducer as F, and determine the amount of sound field gathered by the transducer The equivalent beam waist radius ω 0 and the equivalent half focal depth Z 0 ; 3) According to the above parameters, calculate the pixel value of a single pixel point scanned by the transducer, and optimize the original delay stacking algorithm, A new type of delay stacking algorithm is obtained. Through the new type of delay stacking algorithm, the image pixels obtained by scanning are processed and sorted according to the time length of the ultrasonic signals corresponding to the pixels received by the transducer. Get the final image.

本发明相对于现有的算法不仅考虑了聚焦声场的实际形状,而且还通过对多个换能器位置采集的信号进行延时叠加和重建,实现了内窥探头在切向上的动态聚焦,从而有效的扩大其成像的景深,从而得到全方位的高分辨声聚焦光声内窥成像结果。Compared with the existing algorithm, the present invention not only considers the actual shape of the focused sound field, but also realizes the dynamic focusing of the endoscopic probe in the tangential direction by delaying, superimposing and reconstructing the signals collected by multiple transducer positions. Effectively expand the depth of field of its imaging, so as to obtain a full range of high-resolution acoustic focusing photoacoustic endoscopic imaging results.

本发明的基于定焦声场的光声/超声内窥探头及动态聚焦重建算法中的一些实施例,在步骤3中,利用上述步骤1和步骤2中得到的数据参数,通过以下公式获得最终的图像:In some embodiments of the photoacoustic/ultrasonic endoscope probe and dynamic focus reconstruction algorithm based on the fixed-focus sound field of the present invention, in step 3, using the data parameters obtained in the above steps 1 and 2, the final result is obtained by the following formula image:

Figure BDA0003203429890000071
Figure BDA0003203429890000071

其中,Nd为换能器扫描过的方向位置数,

Figure BDA0003203429890000072
为需要重建的像素坐标,S(i,t)为换能器在第i个位置时换能器接收到的信号,t为时间,
Figure BDA0003203429890000073
为第i个像素点所在位置的相位系数,
Figure BDA0003203429890000074
为第i个像素点所在位置的权重系数。Among them, Nd is the number of directional positions scanned by the transducer,
Figure BDA0003203429890000072
is the pixel coordinate to be reconstructed, S(i,t) is the signal received by the transducer when the transducer is at the i-th position, t is the time,
Figure BDA0003203429890000073
is the phase coefficient at the position of the i-th pixel,
Figure BDA0003203429890000074
is the weight coefficient of the position of the i-th pixel.

其中,相位系数由像素点与换能器焦点的相对位置所决定,在光声成像中,相位系数为:Among them, the phase coefficient is determined by the relative position of the pixel point and the transducer focus. In photoacoustic imaging, the phase coefficient is:

Figure BDA0003203429890000081
Figure BDA0003203429890000081

在超声成像中,相位系数为:In ultrasound imaging, the phase coefficient is:

Figure BDA0003203429890000082
Figure BDA0003203429890000082

其中,v为超声波在介质中的传播速度,a为像素点相对于所述换能器焦点的轴向距离,b为像素点相对于所述换能器焦点的径向距离,像素点距离换能器焦点的轴向距离a和径向距离b的表达式如下:Among them, v is the propagation speed of ultrasonic waves in the medium, a is the axial distance of the pixel point relative to the focal point of the transducer, b is the radial distance of the pixel point relative to the focal point of the transducer, and the pixel point distance The expressions for the axial distance a and the radial distance b of the energizer focus are as follows:

Figure BDA0003203429890000083
Figure BDA0003203429890000083

Figure BDA0003203429890000084
Figure BDA0003203429890000084

其中,

Figure BDA0003203429890000085
为换能器在第i个位置时的焦点,
Figure BDA0003203429890000086
为换能器第i个位置时的旋转中心,
Figure BDA0003203429890000087
为第i个像素所在位置的权重系数,所述权重系数可通过如下的公式计算得出:in,
Figure BDA0003203429890000085
is the focus of the transducer at the i-th position,
Figure BDA0003203429890000086
is the rotation center of the transducer at the i-th position,
Figure BDA0003203429890000087
is the weight coefficient of the position of the i-th pixel, and the weight coefficient can be calculated by the following formula:

Figure BDA0003203429890000088
Figure BDA0003203429890000088

其中,

Figure BDA0003203429890000089
为一个由当前像素点与换能器的距离所决定的位置限制函数,该位置限制函数通过以下方式来限制延时叠加的范围:in,
Figure BDA0003203429890000089
is a position limit function determined by the distance between the current pixel and the transducer. The position limit function limits the range of the delay stacking in the following ways:

Figure BDA00032034298900000810
Figure BDA00032034298900000810

其中,a1、a2、b1和b2为位置限制函数的限制参数,所述限制参数由以下公式计算得出:Wherein, a 1 , a 2 , b 1 and b 2 are the limit parameters of the position limit function, and the limit parameters are calculated by the following formula:

Figure BDA0003203429890000091
Figure BDA0003203429890000091

a2=R,b1=-R,b2=R,a 2 =R, b 1 =-R, b 2 =R,

其中,D为内窥探头的直径,R为成像区域的半径。通过上述的内窥成像算法可以看出,相对于现有正在使用的成像算法,本发明所采用的新型延时叠加算法在常规的延时叠加算法的基础上,对其进行优化,从而得到新型的延时叠加算法,该新型延时叠加算法不仅保留了传统延时叠加算法在成像的计算速度快,成像效率高的效果,还提高了最终的成像的清晰度,从而使得采用本发明的算法在内窥成像领域既能够保证成像效率,又能够提高成像的清晰准确。Among them, D is the diameter of the endoscopic probe, and R is the radius of the imaging area. It can be seen from the above-mentioned endoscopic imaging algorithm that, compared with the existing imaging algorithms in use, the novel time-lapse stacking algorithm adopted in the present invention is optimized on the basis of the conventional time-lapse stacking algorithm, thereby obtaining a new type of delay stacking algorithm. The new time-lapse superposition algorithm not only retains the effects of the traditional time-lapse superposition algorithm in imaging with fast calculation speed and high imaging efficiency, but also improves the final imaging clarity, so that the algorithm of the present invention is adopted. In the field of endoscopic imaging, it can not only ensure the imaging efficiency, but also improve the clarity and accuracy of imaging.

如图1所示,为现有的内窥成像算法所使用的内窥成像探头,包括透明外壳1,中空声聚焦超声传感器2,多模光纤3,以及顶部45度反射镜4部件,其用来进行光和超声的反射。在扫描时,随着顶部反射镜6的旋转,其可以进行侧向的环形扫描成像。As shown in Fig. 1, the endoscopic imaging probe used by the existing endoscopic imaging algorithm includes a transparent casing 1, a hollow acoustically focused ultrasonic sensor 2, a multi-mode optical fiber 3, and a top 45-degree reflector 4 components. to reflect light and ultrasound. During scanning, as the top mirror 6 rotates, it can perform lateral annular scanning imaging.

在图1中,等效束腰半径ω0和等效半焦深Z0采用如下公式得出:In Fig. 1, the equivalent beam waist radius ω 0 and the equivalent semi-focal depth Z 0 are obtained by the following formulas:

Figure BDA0003203429890000092
Figure BDA0003203429890000092

其中f0为内窥探头的中心频率,NA为换能器的数值孔径。where f 0 is the center frequency of the endoscopic probe, and NA is the numerical aperture of the transducer.

如图中所示,L1表示中空聚焦超声传感器与反射镜6中心之间的距离,L2表示反射镜6中心与中空聚焦超声传感器焦点之间的距离,在普通的声聚焦光声内窥重建算法中,其通常采用类似于B超成像中所用的超声束直线传播模型,在这种算法中,图像上用极坐标表示的任意一点(θ,r)的重建数值仅仅由单个探测器的信号决定:As shown in the figure, L1 represents the distance between the hollow focused ultrasonic sensor and the center of the mirror 6, and L2 represents the distance between the center of the mirror 6 and the focal point of the hollow focused ultrasonic sensor. In the ordinary acoustic focusing photoacoustic endoscopic reconstruction algorithm In this algorithm, the reconstruction value of any point (θ, r) represented by polar coordinates on the image is only determined by the signal of a single detector. :

I(θ,r)=Sθ(r/v))I(θ,r)=S θ (r/v))

其中Sθ(t)表示在θ角度的时候换能器收到的信号,v为超声在介质中的传播速度,然而,由于其没有考虑超声传感器的探测面的尺寸以及形状,造成了其只能在聚焦的位置得到最好的分辨率结果。where S θ (t) represents the signal received by the transducer at the angle θ, and v is the propagation speed of ultrasound in the medium. However, since it does not consider the size and shape of the detection surface of the ultrasonic sensor, it only The best resolution results are obtained where the focus is.

如图2至图3所示,为本发明中的内窥成像算法所使用的内窥成像探头,该探头包括透明外壳1,氟化钙镜片5,反射镜6,保护结构7以及使用的平场超声换能器8,探头整体在旋转电机驱动下进行运动,并完成侧向的环形扫描成像作业,其中,图中L3表示平场超声传感器8与氟化钙镜片5中心之间的距离,L4表示氟化钙镜片5中心与反射镜6中心之间的距离,L5表示抛物面反射镜6的焦距,在扫描时,由于反射镜6的存在,相当于平场超声传感器8处于图中虚线所标示的镜像位置,同时该反射镜6实现了声透镜的功能,因此在反射镜6后端的声场为聚焦声场,如图2所示,因此,在扫描过程中,与现有算法中所使用的内窥成像设备相对比,相当于现有的内窥成像设备中的中空聚焦超声传感器2围绕着旋转中心在以L3+L4为半径的扫描路径上进行旋转扫描。As shown in FIG. 2 to FIG. 3, it is an endoscopic imaging probe used by the endoscopic imaging algorithm in the present invention, the probe includes a transparent casing 1, a calcium fluoride lens 5, a reflector 6, a protective structure 7 and a flat surface used for In the field ultrasonic transducer 8, the probe moves as a whole under the drive of the rotating motor, and completes the lateral annular scanning imaging operation, wherein, L3 in the figure represents the distance between the flat field ultrasonic sensor 8 and the center of the calcium fluoride lens 5, L4 represents the distance between the center of the calcium fluoride lens 5 and the center of the reflector 6, and L5 represents the focal length of the parabolic reflector 6. During scanning, due to the presence of the reflector 6, it is equivalent to the plan ultrasonic sensor 8 being located in the dotted line in the figure. Marked mirror position, at the same time, the reflector 6 realizes the function of the acoustic lens, so the sound field at the rear end of the reflector 6 is a focused sound field, as shown in Figure 2. Therefore, during the scanning process, the same Compared with the endoscopic imaging device, the hollow focused ultrasonic sensor 2 in the existing endoscopic imaging device performs rotational scanning on a scanning path with a radius of L3+L4 around the center of rotation.

该氟化钙镜片5可透光反声,当进行光声检测时,采用脉冲激光对组织进行扫面探测时,脉冲激光可透过氟化钙镜片5照射到反射镜6上,该反射镜6为定制的抛物面反射镜6,该抛物面反射镜6可进行光声和超声的反射,当脉冲激光透过氟化钙镜片5照射到抛物面反射镜6上后,脉冲激光经过抛物面反射镜6的聚焦,将脉冲激光汇聚于聚焦声场的等效焦点,并以此确定聚集声场的等效束腰半径ω0和等效半焦深Z0,当脉冲激光作用在组织上时,被照射的组织迅速做出反应,反馈出的超声信号被平场超声换能器8所接收,从而具体的反馈信息的数据参数。The calcium fluoride lens 5 can transmit light and reflect sound. When performing photoacoustic detection, when a pulsed laser is used to scan the tissue for detection, the pulsed laser can pass through the calcium fluoride lens 5 and irradiate on the reflecting mirror 6. The reflecting mirror 6 is a customized parabolic reflector 6, which can reflect photoacoustic and ultrasonic waves. When the pulsed laser passes through the calcium fluoride lens 5 and irradiates the parabolic reflector 6, the pulsed laser passes through the parabolic reflector 6. Focusing, the pulsed laser is focused on the equivalent focus of the focused sound field, and the equivalent beam waist radius ω 0 and the equivalent half focal depth Z 0 of the focused sound field are determined accordingly. When the pulsed laser acts on the tissue, the irradiated tissue Respond quickly, and the feedback ultrasonic signal is received by the flat-field ultrasonic transducer 8, so as to specify the data parameters of the feedback information.

当进行超声探测时,平场超声换能器8发出超声进行探测时,发出的超声被透光反声的氟化钙镜片5所折射,经过折射的超声再次由抛物面反射镜6所聚焦,最终作用在组织上,组织受到超声的作用反馈出超声信号,该超声信号最终由平场超声换能器8所接收,并得到具体的反馈信息参数。When performing ultrasonic detection, when the flat-field ultrasonic transducer 8 emits ultrasonic waves for detection, the emitted ultrasonic waves are refracted by the light-transmitting and reflective calcium fluoride lens 5, and the refracted ultrasonic waves are again focused by the parabolic reflector 6, and finally Acting on the tissue, the tissue is subjected to the action of ultrasound to feed back an ultrasound signal, which is finally received by the flat-field ultrasound transducer 8, and specific feedback information parameters are obtained.

以图4所示的情况为例,此探测背景的中空区域直径为5mm,整体外围直径为23mm,在x方向上,在3mm到11mm的位置均匀分布有9个点目标。采集时所用的换能器直径为10mm,焦距为20mm,焦点到旋转中心的距离L5为7mm。聚焦超声传感器中心频率为15MHz,带宽为75%。Taking the situation shown in Figure 4 as an example, the diameter of the hollow area of the detection background is 5mm, and the overall outer diameter is 23mm. In the x-direction, there are 9 point targets evenly distributed at positions of 3mm to 11mm. The diameter of the transducer used for acquisition is 10mm, the focal length is 20mm, and the distance L5 from the focal point to the center of rotation is 7mm. The focused ultrasound sensor has a center frequency of 15MHz and a bandwidth of 75%.

图5(a)即用原始声聚焦光声内窥算法得到的光声图像,其在7mm位置的目标的横向分辨率最高,离开焦点区域越远,目标的横向分辨率扩展越大,系统的景深受限很大。Figure 5(a) is the photoacoustic image obtained by the original acoustic focusing photoacoustic endoscopic algorithm. The lateral resolution of the target at the 7mm position is the highest. The depth of field is very limited.

图5(b)为采用本发明中的基于定焦声场的光声内窥动态聚焦重建算法得到的光声图像,对比图5(a)可以看出,采用本发明中算法除7mm位置目标的横向分辨率和图5(a)相差不大以外,图中其余目标的横向分辨率都得到了显著的提高,由此可以看出本发明所使用的算法在光声成像中能够产生较为明显的提高。Fig. 5(b) is a photoacoustic image obtained by using the photoacoustic endoscopic dynamic focusing reconstruction algorithm based on the fixed-focus sound field in the present invention. Compared with Fig. 5(a), it can be seen that the algorithm in the present invention is used to divide the 7mm position target. Except that the lateral resolution is not much different from Figure 5(a), the lateral resolutions of other targets in the figure have been significantly improved. It can be seen that the algorithm used in the present invention can produce relatively obvious effects in photoacoustic imaging. improve.

图6(a)即用原始声聚焦光声内窥算法得到的超声图像,其在7mm位置的目标的横向分辨率最高,离开焦点区域越远,目标的横向分辨率扩展越大,系统的景深受限很大。Figure 6(a) is the ultrasound image obtained by the original acoustic focusing photoacoustic endoscopic algorithm. The lateral resolution of the target at the 7mm position is the highest. Very limited.

图6(b)为采用本发明中的基于定焦声场的超声内窥动态聚焦重建算法得到的超声图像,对比图6(a)可以看出,采用本发明中算法除7mm位置目标的横向分辨率和图6(a)相差不大以外,图中其余目标的横向分辨率都得到了显著的提高,由此可以看出本发明所使用的算法在超声成像的领域能够产生较为明显的提高。Fig. 6(b) is an ultrasonic image obtained by using the ultrasonic endoscopic dynamic focusing reconstruction algorithm based on the fixed-focus sound field in the present invention. Compared with Fig. 6(a), it can be seen that the algorithm in the present invention is used to divide the lateral resolution of the 7mm position target. 6(a), the lateral resolutions of other targets in the figure have been significantly improved. It can be seen that the algorithm used in the present invention can significantly improve the field of ultrasound imaging.

通过上述的基于定焦声场的光声/超声内窥动态聚焦重建算法,可以明显看出该算法在对于不同的系统结构中具有普遍适用性,并且在光声与超声内窥实验中,在保留焦区横向分辨率的前提下,显著提升了离焦区域的横向分辨率,由此得到最终的成像更加清晰准确。Through the above-mentioned dynamic focus reconstruction algorithm of photoacoustic/ultrasonic endoscopy based on fixed-focus sound field, it can be clearly seen that the algorithm has universal applicability in different system structures, and in the photoacoustic and endoscopic ultrasound experiments, the On the premise of the lateral resolution of the focal area, the lateral resolution of the out-of-focus area is significantly improved, resulting in a clearer and more accurate final image.

在本发明的描述中,参考术语“一个实施例”、“一些实施例”、“一种实施方式”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本发明中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。In the description of the present invention, reference to the description of the terms "one embodiment", "some embodiments", "one embodiment", etc. means that a particular feature, structure, material or characteristic described in connection with the embodiment or example is contained in in at least one embodiment or example of the present invention. In the present invention, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

以上结合附图详细描述了本发明的优选实施方式,但是,本发明并不限于此。在本发明的技术构思范围内,可以对本发明的技术方案进行多种简单变型,包括各个具体技术特征以任何合适的方式进行组合,为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。但这些简单变型和组合同样应当视为本发明所公开的内容,均属于本发明的保护范围。The preferred embodiments of the present invention have been described above in detail with reference to the accompanying drawings, however, the present invention is not limited thereto. Within the scope of the technical concept of the present invention, a variety of simple modifications can be made to the technical solutions of the present invention, including combining various specific technical features in any suitable manner. No further explanation is required. However, these simple modifications and combinations should also be regarded as the contents disclosed in the present invention, and all belong to the protection scope of the present invention.

Claims (10)

1.一种基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,包括以下步骤:1. a photoacoustic/ultrasonic endoscopic dynamic focusing reconstruction algorithm based on fixed-focus sound field, is characterized in that, comprises the following steps: 1)根据换能器的固有参数,确定所述换能器的等效聚焦声场和等效焦点;1) According to the inherent parameters of the transducer, determine the equivalent focused sound field and equivalent focus of the transducer; 2)以等效焦点为基准点,确定所述换能器接收到的等效焦点反馈的信号所经过的路径长度为F,并确定出所述换能器聚集声场的等效束腰半径ω0和等效半焦深Z02) Taking the equivalent focus as the reference point, determine the path length of the signal fed back by the equivalent focus received by the transducer as F, and determine the equivalent beam waist radius ω of the focused sound field of the transducer 0 and the equivalent semi-focal depth Z 0 ; 3)根据上述参数,计算出所述换能器扫描得到的单个像素点的像素值,并将原始的延时叠加算法进行优化,得到新型的延时叠加算法,通过该新型延时叠加算法,将扫描得到的图像像素点进行处理,以所述换能器接收到的该像素点所对应的超声信号的时间长短进行排序,从而获得最终的图像。3) According to the above parameters, calculate the pixel value of a single pixel point scanned by the transducer, and optimize the original delay stacking algorithm to obtain a new type of delay stacking algorithm. Through this new type of delay stacking algorithm, The image pixels obtained by scanning are processed and sorted according to the time length of the ultrasonic signals corresponding to the pixels received by the transducer, so as to obtain a final image. 2.根据权利要求1所述的基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,在步骤3中,通过如下公式得到最终的图像:2. The photoacoustic/ultrasonic endoscopic dynamic focusing reconstruction algorithm based on fixed-focus sound field according to claim 1, is characterized in that, in step 3, obtain the final image by following formula:
Figure FDA0003203429880000011
Figure FDA0003203429880000011
其中,Nd为换能器扫描过的方向位置数,
Figure FDA0003203429880000012
为需要重建的像素坐标,S(i,t)为换能器在第i个位置时换能器接收到的信号,t为时间,
Figure FDA0003203429880000013
为第i个像素点所在位置的相位系数,
Figure FDA0003203429880000014
为第i个像素点所在位置的权重系数。
Among them, Nd is the number of directional positions scanned by the transducer,
Figure FDA0003203429880000012
is the pixel coordinate to be reconstructed, S(i,t) is the signal received by the transducer when the transducer is at the i-th position, t is the time,
Figure FDA0003203429880000013
is the phase coefficient at the position of the i-th pixel,
Figure FDA0003203429880000014
is the weight coefficient of the position of the i-th pixel.
3.根据权利要求2所述的基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,所述相位系数由像素点与换能器焦点的相对位置所决定,在光声成像中,相位系数为:3. The photoacoustic/ultrasonic endoscopy dynamic focusing reconstruction algorithm based on fixed-focus sound field according to claim 2, wherein the phase coefficient is determined by the relative position of the pixel point and the focus of the transducer. In imaging, the phase coefficient is:
Figure FDA0003203429880000015
Figure FDA0003203429880000015
在超声成像中,相位系数为:In ultrasound imaging, the phase coefficient is:
Figure FDA0003203429880000021
Figure FDA0003203429880000021
其中,v为超声波在介质中的传播速度,a为像素点相对于所述换能器焦点的轴向距离,b为像素点相对于所述换能器焦点的径向距离。Wherein, v is the propagation velocity of ultrasonic waves in the medium, a is the axial distance of the pixel point relative to the transducer focus, and b is the radial distance of the pixel point relative to the transducer focus.
4.根据权利要求3所述的基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,所述像素点距离换能器焦点的轴向距离a的表达式如下:4. The photoacoustic/ultrasonic endoscopy dynamic focus reconstruction algorithm based on fixed-focus sound field according to claim 3, wherein the expression of the axial distance a of the pixel point from the transducer focus is as follows:
Figure FDA0003203429880000022
Figure FDA0003203429880000022
其中,
Figure FDA0003203429880000023
为换能器在第i个位置时的焦点,
Figure FDA0003203429880000024
为换能器第i个位置时的旋转中心。
in,
Figure FDA0003203429880000023
is the focus of the transducer at the i-th position,
Figure FDA0003203429880000024
is the rotation center of the transducer at the i-th position.
5.根据权利要求3所述的基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,所述像素点距离换能器焦点的径向距离b的表达式如下:5. the photoacoustic/ultrasonic endoscopy dynamic focusing reconstruction algorithm based on fixed-focus sound field according to claim 3, is characterized in that, the expression of the radial distance b of described pixel point distance transducer focus is as follows:
Figure FDA0003203429880000025
Figure FDA0003203429880000025
其中,
Figure FDA0003203429880000026
为换能器在第i个位置时的焦点,
Figure FDA0003203429880000027
为换能器第i个位置时的旋转中心。
in,
Figure FDA0003203429880000026
is the focus of the transducer at the i-th position,
Figure FDA0003203429880000027
is the rotation center of the transducer at the i-th position.
6.根据权利要求3所述的基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,所述
Figure FDA0003203429880000028
为第i个像素所在位置的权重系数,所述权重系数可通过如下的公式计算得出:
6. The photoacoustic/ultrasonic endoscopic dynamic focusing reconstruction algorithm based on fixed-focus sound field according to claim 3, wherein the
Figure FDA0003203429880000028
is the weight coefficient of the position of the i-th pixel, and the weight coefficient can be calculated by the following formula:
Figure FDA0003203429880000031
Figure FDA0003203429880000031
其中,
Figure FDA0003203429880000032
为一个由当前像素点与换能器的距离所决定的位置限制函数。
in,
Figure FDA0003203429880000032
is a position limit function determined by the distance between the current pixel and the transducer.
7.根据权利要求6所述的基于定焦声场的光声/超声内窥动态聚焦重建算法,其特征在于,所述位置限制函数通过以下方式来限制延时叠加的范围:7. The photoacoustic/ultrasonic endoscopic dynamic focus reconstruction algorithm based on fixed-focus sound field according to claim 6, wherein the position limit function limits the range of time delay stacking in the following manner:
Figure FDA0003203429880000033
Figure FDA0003203429880000033
其中,a1、a2、b1和b2为位置限制函数的限制参数,所述限制参数由以下公式计算得出:Wherein, a 1 , a 2 , b 1 and b 2 are the limit parameters of the position limit function, and the limit parameters are calculated by the following formula:
Figure FDA0003203429880000034
Figure FDA0003203429880000034
a2=R,b1=-R,b2=R,a 2 =R, b 1 =-R, b 2 =R, 其中,D为内窥探头的直径,R为成像区域的半径。Among them, D is the diameter of the endoscopic probe, and R is the radius of the imaging area.
8.一种基于定焦声场的光声/超声内窥探头,其特征在于,包括:透明外壳、氟化钙镜片、反射镜、保护结构和换能器,所述内窥探头由旋转电机驱动其运动,所述透明外壳包裹在所述内窥探头的外表面,所述换能器发射的信号波经过氟化钙镜片照射到反射镜上,并通过反射镜的聚焦,照射到探测物上,且所述内窥探头能够实现权利要求1~7中任一项算法。8. A photoacoustic/ultrasonic endoscopic probe based on a fixed-focus sound field, characterized in that it comprises: a transparent casing, a calcium fluoride lens, a reflector, a protective structure and a transducer, and the endoscopic probe is driven by a rotating motor When it moves, the transparent shell is wrapped around the outer surface of the endoscope probe, and the signal wave emitted by the transducer is irradiated to the reflector through the calcium fluoride lens, and is focused on the probe by the reflector. , and the endoscopic probe can implement the algorithm of any one of claims 1 to 7 . 9.根据权利要求8所述的基于定焦声场的光声/超声内窥探头,其特征在于,所述换能器采用平场超声换能器。9 . The photoacoustic/ultrasonic endoscopic probe based on a fixed-focus sound field according to claim 8 , wherein the transducer is a flat-field ultrasonic transducer. 10 . 10.根据权利要求8所述的基于定焦声场的光声/超声内窥探头,其特征在于,所述反射镜可采用抛物面反射镜,且所述抛物面反射镜以弧高为基准线,所述基准线与所述氟化钙镜片的中心相垂直。10. The photoacoustic/ultrasonic endoscope probe based on a fixed-focus sound field according to claim 8, wherein the reflector can be a parabolic reflector, and the parabolic reflector takes the arc height as a reference line, so The reference line is perpendicular to the center of the calcium fluoride lens.
CN202110910237.3A 2021-08-09 2021-08-09 Photoacoustic/ultrasound endoscopic probe and dynamic focus reconstruction algorithm based on fixed-focus sound field Active CN113842158B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110910237.3A CN113842158B (en) 2021-08-09 2021-08-09 Photoacoustic/ultrasound endoscopic probe and dynamic focus reconstruction algorithm based on fixed-focus sound field

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110910237.3A CN113842158B (en) 2021-08-09 2021-08-09 Photoacoustic/ultrasound endoscopic probe and dynamic focus reconstruction algorithm based on fixed-focus sound field

Publications (2)

Publication Number Publication Date
CN113842158A true CN113842158A (en) 2021-12-28
CN113842158B CN113842158B (en) 2023-11-21

Family

ID=78975688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110910237.3A Active CN113842158B (en) 2021-08-09 2021-08-09 Photoacoustic/ultrasound endoscopic probe and dynamic focus reconstruction algorithm based on fixed-focus sound field

Country Status (1)

Country Link
CN (1) CN113842158B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118986257A (en) * 2024-07-26 2024-11-22 江苏省苏北人民医院 Ultrasonic/photoacoustic imaging diagnosis and treatment integrated uterine cavity endoscope based on micro-array probe and working method thereof
CN119184628A (en) * 2024-11-27 2024-12-27 之江实验室 Photoacoustic endoscopic imaging assembly, computer equipment and endoscopic imaging method

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4208916A (en) * 1978-09-13 1980-06-24 Picker Corporation Electronic ultrasonic sector scanning apparatus and method
US20110230766A1 (en) * 2008-10-09 2011-09-22 Signostics Limited Ultrasound imaging modality improvement
KR20160048256A (en) * 2014-10-23 2016-05-04 포항공과대학교 산학협력단 Catheter and system for detecting the ultrasound signal and the photoacoustic signal
WO2017202362A1 (en) * 2016-05-25 2017-11-30 深圳开立生物医疗科技股份有限公司 Ultrasonic imaging method, system, and device
CN108261209A (en) * 2018-01-23 2018-07-10 湖南大学 The method of follow-on high-resolution sound focusing optoacoustic endoscopy imaging back projection imaging
CN110251093A (en) * 2019-07-24 2019-09-20 中南大学 A photoacoustic/ultrasound probe and scanning imaging method for acoustic focusing endoscopy
US20190290939A1 (en) * 2015-10-13 2019-09-26 Arcscan, Inc. Ultrasound therapeutic and scanning apparatus
CN111772581A (en) * 2020-06-11 2020-10-16 华南师范大学 High-sensitivity photoacoustic/ultrasonic dual-mode imaging device and method based on double-curvature linear array detector
CN112450983A (en) * 2020-11-19 2021-03-09 深圳先进技术研究院 Super-resolution imaging method and system

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4208916A (en) * 1978-09-13 1980-06-24 Picker Corporation Electronic ultrasonic sector scanning apparatus and method
US20110230766A1 (en) * 2008-10-09 2011-09-22 Signostics Limited Ultrasound imaging modality improvement
KR20160048256A (en) * 2014-10-23 2016-05-04 포항공과대학교 산학협력단 Catheter and system for detecting the ultrasound signal and the photoacoustic signal
US20190290939A1 (en) * 2015-10-13 2019-09-26 Arcscan, Inc. Ultrasound therapeutic and scanning apparatus
WO2017202362A1 (en) * 2016-05-25 2017-11-30 深圳开立生物医疗科技股份有限公司 Ultrasonic imaging method, system, and device
CN108261209A (en) * 2018-01-23 2018-07-10 湖南大学 The method of follow-on high-resolution sound focusing optoacoustic endoscopy imaging back projection imaging
CN110251093A (en) * 2019-07-24 2019-09-20 中南大学 A photoacoustic/ultrasound probe and scanning imaging method for acoustic focusing endoscopy
CN111772581A (en) * 2020-06-11 2020-10-16 华南师范大学 High-sensitivity photoacoustic/ultrasonic dual-mode imaging device and method based on double-curvature linear array detector
CN112450983A (en) * 2020-11-19 2021-03-09 深圳先进技术研究院 Super-resolution imaging method and system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
WANG, BO: "3D acoustic resolution-based photoacoustic endoscopy with dynamic focusing", 《QUANTITATIVE IMAGING IN MEDICINE AND SURGERY》, pages 685 - 696 *
XIAO, JIAYING: "Lithography of Aluminum Coated PVDF Annular Array for Photoacoustic Endoscopy", 《SENSORS AND ACTUATORS A-PHYSICAL》, pages 561 - 571 *
杨凯旋: "磁光声联合内窥图像仿真的研究", 《中国优秀硕士学位论文全文数据库》 *
魏宁宁: "基于动态聚焦的高分辨率光声内窥成像算法研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118986257A (en) * 2024-07-26 2024-11-22 江苏省苏北人民医院 Ultrasonic/photoacoustic imaging diagnosis and treatment integrated uterine cavity endoscope based on micro-array probe and working method thereof
CN119184628A (en) * 2024-11-27 2024-12-27 之江实验室 Photoacoustic endoscopic imaging assembly, computer equipment and endoscopic imaging method

Also Published As

Publication number Publication date
CN113842158B (en) 2023-11-21

Similar Documents

Publication Publication Date Title
JP4991860B2 (en) Resonant ultrasonic transducer
CN108261209B (en) Improved high-resolution acoustic focusing photoacoustic endoscopic imaging back-projection imaging method
CN104188625B (en) A kind of multi-modal micro imaging system
CN106419803B (en) Intravascular photoacoustic absorption, elasticity, sticky multi-modality imaging integration endoscope and its imaging method
CN105380586B (en) A combined stereoscopic scanning optical and acoustic endoscopic imaging device and method thereof
CN110251093B (en) Acoustic focusing endoscopic photoacoustic/ultrasonic probe and scanning imaging method
CN113842158B (en) Photoacoustic/ultrasound endoscopic probe and dynamic focus reconstruction algorithm based on fixed-focus sound field
CN109431543A (en) Sound wave condenser lens, supersonic imaging device and method
CN104977298A (en) Rotary type optical-acoustic combined imaging probe and conduit
JPH07184898A (en) Ultrasonic probe
JP2004089311A (en) Ultrasonic transmission/reception device
US20180168549A1 (en) Ultrasonic scanner with a multiple faceted mirror
US9161742B2 (en) Ultrasound imaging apparatus
CN108717045A (en) A kind of ultrasonic, three modality imaging system of optoacoustic and fluorescence
CN112450882B (en) Ultrasound probe, endoscope, endoscopic imaging system, and endoscopic imaging method
JP2012516182A (en) Acoustic device for ultrasound imaging
WO2022104701A1 (en) Ultrasound probe, endoscope, endoscopic imaging system, and endoscopic imaging method
JP6925815B2 (en) Ultrasonic diagnostic equipment
CN102824185B (en) Photoacoustic tomography system combined with acoustical transmission reflector and imaging method thereof
CN112444492A (en) Photoacoustic imaging method, computer device, system and storage medium
CN115363528A (en) Reconstruction method for elevation angle resolution of annular array photoacoustic computed tomography system
CN106691377A (en) Photoacoustic microimaging adaptive scanning system and method
CN112244889A (en) Method for determining vibration element array, puncture needle imaging method and ultrasonic equipment
JP2014061124A (en) Photoacoustic measuring apparatus, detection method of light scanning state of photoacoustic measuring apparatus, and sheet-like member for use in method
KR101797910B1 (en) Rotatory linear probe

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wu Hualin

Inventor after: Xiao Jiaying

Inventor after: Wang Yongjun

Inventor after: Jiang Jinsheng

Inventor after: Wang Bo

Inventor after: Ye Tong

Inventor before: Wang Bo

Inventor before: Xiao Jiaying

Inventor before: Wang Yongjun

Inventor before: Jiang Jinsheng

Inventor before: Wu Hualin

Inventor before: Ye Tong

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240527

Address after: Room A510, Virtual University Park Building, No. 19 Gaoxin South Fourth Road, Nanshan District, Shenzhen, Guangdong Province, 518000

Patentee after: SHENZHEN RESEARCH INSTITUTE, CENTRAL SOUTH University

Country or region after: China

Address before: 410083 Central South University headquarters, Yuelu District, Changsha City, Hunan Province

Patentee before: CENTRAL SOUTH University

Country or region before: China