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 PDFInfo
- 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
Links
- 239000000523 sample Substances 0.000 title claims abstract description 36
- 238000003384 imaging method Methods 0.000 claims abstract description 60
- WUKWITHWXAAZEY-UHFFFAOYSA-L calcium difluoride Chemical compound [F-].[F-].[Ca+2] WUKWITHWXAAZEY-UHFFFAOYSA-L 0.000 claims description 14
- 229910001634 calcium fluoride Inorganic materials 0.000 claims description 14
- 238000012285 ultrasound imaging Methods 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 5
- 238000001839 endoscopy Methods 0.000 claims description 4
- 230000001681 protective effect Effects 0.000 claims description 4
- 238000002604 ultrasonography Methods 0.000 description 14
- 238000001514 detection method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000009558 endoscopic ultrasound Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000000034 method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0097—Detecting, 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
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
Description
技术领域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
其中,Nd为换能器扫描过的方向位置数,为需要重建的像素坐标,S(i,t)为换能器在第i个位置时换能器接收到的信号,t为时间,为第i个像素点所在位置的相位系数,为第i个像素点所在位置的权重系数。通过该优选技术方案,将传统的延时叠加算法进行优化,得到新的延时叠加算法,延时叠加算法能够通过简单的计算推导,从而得到最终的图像数据,并且通过该叠加算法的计算,所需要的计算时间更短,效率更高,并且该算法,在保证了计算效率的同时,提高了成像的清晰度,使得最终的成像更加清晰准确。Among them, Nd is the number of directional positions scanned by the transducer, 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, is the phase coefficient at the position of the i-th pixel, 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:
在超声成像中,相位系数为:In ultrasound imaging, the phase coefficient is:
其中,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:
其中,为换能器在第i个位置时的焦点,为换能器第i个位置时的旋转中心。通过该优选技术方案,确定了像素点与换能器焦点之间的轴向距离,从而能够提高图像像素点在轴向距离的方向上成像的准确性。in, is the focus of the transducer at the i-th position, 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:
其中,为换能器在第i个位置时的焦点,为换能器第i个位置时的旋转中心。通过该优选技术方案,确定了像素点与换能器焦点之间的径向距离,从而能够提高图像像素点在径向距离的方向上成像的准确性。in, is the focus of the transducer at the i-th position, 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.
优选地,所述为第i个像素所在位置的权重系数,所述权重系数可通过如下的公式计算得出:Preferably, the is the weight coefficient of the position of the i-th pixel, and the weight coefficient can be calculated by the following formula:
其中,为一个由当前像素点与换能器的距离所决定的位置限制函数。通过该优选技术方案,将权重系数通过位置限制函数进行限制,能够进一步对像素点的位置进行限定,从而能够更加准确的确定像素点的相对位置,进而保证了叠加后的图像的准确性。in, 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:
其中,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:
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
3 多模光纤 4 顶部45度反射镜3
5 氟化钙镜片 6 反射镜5
7 保护结构 8 平场超声换能器7
具体实施方式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
其中,Nd为换能器扫描过的方向位置数,为需要重建的像素坐标,S(i,t)为换能器在第i个位置时换能器接收到的信号,t为时间,为第i个像素点所在位置的相位系数,为第i个像素点所在位置的权重系数。Among them, Nd is the number of directional positions scanned by the transducer, 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, is the phase coefficient at the position of the i-th pixel, 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:
在超声成像中,相位系数为:In ultrasound imaging, the phase coefficient is:
其中,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:
其中,为换能器在第i个位置时的焦点,为换能器第i个位置时的旋转中心,为第i个像素所在位置的权重系数,所述权重系数可通过如下的公式计算得出:in, is the focus of the transducer at the i-th position, is the rotation center of the transducer at the i-th position, is the weight coefficient of the position of the i-th pixel, and the weight coefficient can be calculated by the following formula:
其中,为一个由当前像素点与换能器的距离所决定的位置限制函数,该位置限制函数通过以下方式来限制延时叠加的范围:in, 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:
其中,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:
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
在图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:
其中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
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
该氟化钙镜片5可透光反声,当进行光声检测时,采用脉冲激光对组织进行扫面探测时,脉冲激光可透过氟化钙镜片5照射到反射镜6上,该反射镜6为定制的抛物面反射镜6,该抛物面反射镜6可进行光声和超声的反射,当脉冲激光透过氟化钙镜片5照射到抛物面反射镜6上后,脉冲激光经过抛物面反射镜6的聚焦,将脉冲激光汇聚于聚焦声场的等效焦点,并以此确定聚集声场的等效束腰半径ω0和等效半焦深Z0,当脉冲激光作用在组织上时,被照射的组织迅速做出反应,反馈出的超声信号被平场超声换能器8所接收,从而具体的反馈信息的数据参数。The
当进行超声探测时,平场超声换能器8发出超声进行探测时,发出的超声被透光反声的氟化钙镜片5所折射,经过折射的超声再次由抛物面反射镜6所聚焦,最终作用在组织上,组织受到超声的作用反馈出超声信号,该超声信号最终由平场超声换能器8所接收,并得到具体的反馈信息参数。When performing ultrasonic detection, when the flat-field
以图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)
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)
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)
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 |
-
2021
- 2021-08-09 CN CN202110910237.3A patent/CN113842158B/en active Active
Patent Citations (9)
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)
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)
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 |