CN105572049B - Optoacoustic quantifies elastograph imaging method and device - Google Patents
Optoacoustic quantifies elastograph imaging method and device Download PDFInfo
- Publication number
- CN105572049B CN105572049B CN201510944323.0A CN201510944323A CN105572049B CN 105572049 B CN105572049 B CN 105572049B CN 201510944323 A CN201510944323 A CN 201510944323A CN 105572049 B CN105572049 B CN 105572049B
- Authority
- CN
- China
- Prior art keywords
- tissue sample
- point
- computer
- signal
- quantifies
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 19
- 238000010168 coupling process Methods 0.000 claims abstract description 57
- 238000005859 coupling reaction Methods 0.000 claims abstract description 57
- 230000008878 coupling Effects 0.000 claims abstract description 56
- 239000007788 liquid Substances 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000005284 excitation Effects 0.000 claims abstract description 11
- 230000005540 biological transmission Effects 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims description 23
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000002033 PVDF binder Substances 0.000 claims description 5
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 5
- 239000010931 gold Substances 0.000 claims description 5
- 229910052737 gold Inorganic materials 0.000 claims description 5
- 229920002981 polyvinylidene fluoride Polymers 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000002604 ultrasonography Methods 0.000 claims description 3
- 239000012528 membrane Substances 0.000 claims description 2
- 230000008014 freezing Effects 0.000 claims 3
- 238000007710 freezing Methods 0.000 claims 3
- 230000005611 electricity Effects 0.000 claims 1
- 238000009434 installation Methods 0.000 claims 1
- 239000007787 solid Substances 0.000 claims 1
- 238000002091 elastography Methods 0.000 abstract description 14
- 238000005259 measurement Methods 0.000 abstract description 7
- 230000001066 destructive effect Effects 0.000 abstract description 3
- 230000003760 hair shine Effects 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 70
- 229920001817 Agar Polymers 0.000 description 19
- 239000008272 agar Substances 0.000 description 19
- 238000001514 detection method Methods 0.000 description 19
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000003902 lesion Effects 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000002099 shear wave elastography Methods 0.000 description 2
- 208000026062 Tissue disease Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 208000019425 cirrhosis of liver Diseases 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000001766 physiological effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011896 sensitive detection Methods 0.000 description 1
- 210000001685 thyroid gland Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
- G01N2021/1706—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids in solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
- G01N2021/1708—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids with piezotransducers
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种光声定量弹性成像方法及装置,所述方法包括:激光器发出脉冲激光,该脉冲激光通过聚焦透镜聚焦,照在组织样品上,激发出光声信号,光声信号经过耦合槽的耦合液后被超声探测器接收;超声探测器接收的光声信号经放大器放大后,被示波器采集,将信号信息存储到计算机中,计算机控制步进电机逐点移动组织样品,对应的组织样品上形成X‑Y二维平面扫描区域;示波器采集完全部信号后,计算机计算各点的组织样品定量弹性模量;根据计算的定量弹性模量,重建出组织样品的定量弹性二维图像;所述装置包括光声激发源、信号采集/传输/重建组件、耦合槽、步进电机和X‑Y二维扫描平台。本发明可以实现无损、高分辨的组织弹性定量测量和成像。
The invention discloses a method and device for photoacoustic quantitative elastography. The method comprises: a laser emits pulsed laser light, the pulsed laser light is focused through a focusing lens, and shines on a tissue sample to excite a photoacoustic signal, and the photoacoustic signal passes through a coupling groove The coupling liquid is received by the ultrasonic detector; the photoacoustic signal received by the ultrasonic detector is amplified by the amplifier, collected by the oscilloscope, and the signal information is stored in the computer. The computer controls the stepping motor to move the tissue sample point by point, and the corresponding tissue sample Form an X-Y two-dimensional plane scanning area; after the oscilloscope collects all the signals, the computer calculates the quantitative elastic modulus of the tissue sample at each point; according to the calculated quantitative elastic modulus, the quantitative elastic two-dimensional image of the tissue sample is reconstructed; The device includes a photoacoustic excitation source, a signal acquisition/transmission/reconstruction component, a coupling slot, a stepper motor and an X-Y two-dimensional scanning platform. The invention can realize non-destructive, high-resolution quantitative measurement and imaging of tissue elasticity.
Description
技术领域technical field
本发明涉及一种弹性成像方法及装置,尤其是一种光声定量弹性成像方法及装置,属于生物医学检测技术领域。The invention relates to an elastography method and device, in particular to a photoacoustic quantitative elastography method and device, belonging to the technical field of biomedical detection.
背景技术Background technique
组织的病理变化往往伴随其机械性质特别是弹性模量的变化,因此可以通过分析组织的物理参数实现组织的病变检测。目前,通过检测组织的弹性预测组织疾病的方法已经被医生及影像学研究者广泛地接收,如利用超声弹性检测的方法实现无损肝硬化的检测方法已经应用于临床疾病的参数检测。The pathological changes of tissues are often accompanied by the changes of their mechanical properties, especially the elastic modulus, so the detection of tissue lesions can be realized by analyzing the physical parameters of tissues. At present, the method of predicting tissue diseases by detecting tissue elasticity has been widely accepted by doctors and imaging researchers. For example, the detection method of non-destructive liver cirrhosis using ultrasound elasticity detection method has been applied to the parameter detection of clinical diseases.
现有的定量弹性检测方法主要为静态超声弹性检测和超声剪切波弹性成像。静态超声弹性检测方法采用静态或准静态的组织激励方法,人为利用已知大小的力挤压组织是组织产生形变,利用超声成像获得组织中的位移大小,根据弹性模量等于(力/位移大小),可以获得定量弹性模量。该方法的缺点是,将组织建模为纯弹性体,忽略组织具有粘弹性质的事实,不可避免的会导致测量结果不准确。同时该方法中施加在组织上的力在组织中的分布受到组织样品形状及激励手段的影响,导致测量结果与真实值偏差较大。超声剪切波弹性成像利用聚焦超声在焦点处引起的非线性效应激励剪切波,通过测量剪切波在组织中的传播速度获得组织的弹性模量。该方法主要应用于如乳腺及甲状腺等大组织器官的弹性检测,存在激发效率较低对小病灶区检测不准确等缺点,且聚焦区域能量过高导致温度升高对生物组织造成损伤或者探测区域生理性质的改变不可避免。同时。目前的弹性定量检测和成像主要能够在大器官能够实现,对于血管疾病的弹性成像还停留在成像应变图而无法实现定量弹性成像的程度。本文提出的组织弹性模量的光声定量测量与成像方法能够实现无损小病灶区的定量弹性测量和成像,与已有弹性测量和成像技术相比具有无法比拟的优越性。The existing quantitative elastic testing methods are mainly static ultrasonic elastic testing and ultrasonic shear wave elastography. The static ultrasonic elastic detection method adopts a static or quasi-static tissue excitation method, artificially uses a known force to squeeze the tissue to cause deformation of the tissue, and uses ultrasonic imaging to obtain the displacement in the tissue. According to the elastic modulus equal to (force/displacement ), the quantitative elastic modulus can be obtained. The disadvantage of this approach is that modeling the tissue as a pure elastic body, ignoring the fact that the tissue has viscoelastic properties, inevitably leads to inaccurate measurements. At the same time, the distribution of the force exerted on the tissue in this method is affected by the shape of the tissue sample and the excitation means, resulting in a large deviation between the measurement result and the true value. Ultrasonic shear wave elastography uses the nonlinear effect induced by focused ultrasound at the focal point to excite shear waves, and obtains the elastic modulus of tissues by measuring the propagation velocity of shear waves in tissues. This method is mainly used in the elastic detection of large tissues and organs such as breast and thyroid gland. It has the disadvantages of low excitation efficiency and inaccurate detection of small lesions, and the high energy of the focused area will cause damage to biological tissues or detection areas due to temperature rise. Changes in physiological properties are inevitable. at the same time. The current quantitative elasticity detection and imaging can mainly be realized in large organs, and the elasticity imaging of vascular diseases still stays at the level of imaging strain maps and cannot achieve quantitative elasticity imaging. The photoacoustic quantitative measurement and imaging method of tissue elastic modulus proposed in this paper can realize quantitative elastic measurement and imaging of small lesions without damage, and has incomparable advantages compared with existing elastic measurement and imaging techniques.
目前,利用光声方法检测组织弹性在我国已有报道,如2011年9月7日公开的发明专利:光声弹性成像方法与装置,申请人:华南师范大学,申请日:2011年1月14日,申请号:201110008213.5,其采用了强度调制的连续光源激发产生光声信号,通过测量该信号与调制信号之间的相位差并逐点扫描可以重建出检测组织样品的弹性分布图像。但在上述方法中,未能够实现组织弹性模量的定量测量,也就是不能获得弹性模量的绝对值,只能提供相对值,降低了该方法在实际应用中测量结果的准确性。At present, the use of photoacoustic method to detect tissue elasticity has been reported in my country, such as the invention patent published on September 7, 2011: photoacoustic elastography method and device, applicant: South China Normal University, application date: January 14, 2011 Japan, application number: 201110008213.5, which uses an intensity-modulated continuous light source to excite a photoacoustic signal. By measuring the phase difference between the signal and the modulated signal and scanning point by point, the elastic distribution image of the tissue sample can be reconstructed. However, in the above method, the quantitative measurement of the tissue elastic modulus cannot be realized, that is, the absolute value of the elastic modulus cannot be obtained, and only relative values can be provided, which reduces the accuracy of the measurement results of this method in practical applications.
发明内容Contents of the invention
本发明的目的是为了解决上述现有技术的缺陷,提供了一种光声定量弹性成像方法,该方法可以实现无损、高分辨的组织弹性定量测量和成像。The object of the present invention is to solve the above-mentioned defects in the prior art, and provide a photoacoustic quantitative elastography method, which can realize non-destructive, high-resolution quantitative measurement and imaging of tissue elasticity.
本发明的另一目的在于提供一种实现上述方法的光声定量弹性成像装置。Another object of the present invention is to provide a photoacoustic quantitative elastography device for realizing the above method.
本发明的目的可以通过采取如下技术方案达到:The purpose of the present invention can be achieved by taking the following technical solutions:
光声定量弹性成像方法,所述方法包括以下步骤:A photoacoustic quantitative elastography method, the method comprising the steps of:
1)将组织样品放置在X-Y二维扫描平台上,并浸没在耦合槽的耦合液中;将聚焦透镜设置在组织样品正上方,并调整聚焦透镜的高度使聚焦透镜的焦点不偏离组织样品表面;将超声探测器对准组织样品,并使超声探测器的下端进入耦合槽的耦合液中;1) Place the tissue sample on the X-Y two-dimensional scanning platform and immerse it in the coupling liquid of the coupling tank; set the focusing lens directly above the tissue sample, and adjust the height of the focusing lens so that the focus of the focusing lens does not deviate from the surface of the tissue sample ;Aim the ultrasonic probe at the tissue sample, and make the lower end of the ultrasonic probe enter the coupling liquid in the coupling tank;
2)激光器发出脉冲激光,该脉冲激光通过聚焦透镜聚焦,照在组织样品上,激发出光声信号,光声信号经过耦合槽的耦合液后被超声探测器接收;2) The laser emits pulsed laser light, which is focused by the focusing lens and irradiated on the tissue sample to excite a photoacoustic signal, which is received by the ultrasonic detector after passing through the coupling liquid in the coupling tank;
3)超声探测器接收的光声信号经放大器放大后,被示波器采集,示波器将采集的信号信息存储到计算机中,计算机控制步进电机逐点移动组织样品,对应的组织样品上形成X-Y二维平面扫描区域,步进电机每移动一次,示波器就进行一次信号采集;3) The photoacoustic signal received by the ultrasonic detector is amplified by the amplifier and collected by the oscilloscope. The oscilloscope stores the collected signal information in the computer. The computer controls the stepping motor to move the tissue sample point by point, and forms an X-Y two-dimensional image on the corresponding tissue sample. In the plane scanning area, every time the stepper motor moves, the oscilloscope performs a signal acquisition;
4)示波器采集完全部信号后,计算机将各点的信号对时间进行一次积分,得到组织样品表面振动位移的时间函数;获取组织样品表面振动位移从零上升至其最大值时所需时间,利用该时间计算各点的组织样品定量弹性模量;根据计算的定量弹性模量,重建出组织样品的定量弹性二维图像。4) After the oscilloscope has collected all the signals, the computer integrates the signals at each point with respect to time to obtain the time function of the vibration displacement of the tissue sample surface; obtain the time required for the surface vibration displacement of the tissue sample to rise from zero to its maximum value, and use At this time, the quantitative elastic modulus of the tissue sample at each point is calculated; according to the calculated quantitative elastic modulus, the quantitative elastic two-dimensional image of the tissue sample is reconstructed.
作为一种优选方案,步骤1)中,所述使超声探测器的下端进入耦合槽的耦合液中具体为:使超声探测器的下端进入耦合槽的耦合液深度5-8mm处。As a preferred solution, in step 1), the step of allowing the lower end of the ultrasonic probe to enter the coupling liquid in the coupling tank specifically includes: allowing the lower end of the ultrasonic probe to enter the coupling liquid in the coupling tank at a depth of 5-8mm.
作为一种优选方案,步骤3)中,所述计算机控制步进电机逐点移动组织样品是指:计算机利用Labview程序控制步进电机逐点移动组织样品;步骤4)中,所述计算机将各点的信号对时间进行一次积分是指:计算机利用Matlab程序将各点的信号对时间进行一次积分。As a preferred solution, in step 3), the computer controls the stepping motor to move the tissue sample point by point refers to: the computer uses the Labview program to control the stepping motor to move the tissue sample point by point; in step 4), the computer moves each Integrating the signal of each point with respect to time means that the computer uses the Matlab program to integrate the signal of each point with respect to time once.
作为一种优选方案,步骤4)中,所述计算各点的组织样品定量弹性模量,采用下式:As a preferred solution, in step 4), the quantitative elastic modulus of the tissue sample at each point of the calculation adopts the following formula:
其中,ρ为生物组织密度,R为激光光斑半径,tmax为组织样品表面振动位移从零上升至其最大值时所需时间。Among them, ρ is the density of biological tissue, R is the radius of the laser spot, and t max is the time required for the surface vibration displacement of the tissue sample to rise from zero to its maximum value.
本发明的另一目的可以通过采取如下技术方案达到:Another object of the present invention can be achieved by taking the following technical solutions:
光声定量弹性成像装置,所述装置包括光声激发源、信号采集/传输/重建组件、耦合槽、步进电机和X-Y二维扫描平台,所述光声激发源包括激光器和聚焦透镜;所述信号采集/传输/重建组件包括超声探测器、放大器、示波器和计算机,所述超声探测器、放大器、示波器和计算机依次相连,所述计算机安装有采集控制及信号处理系统;所述步进电机与计算机相连,所述X-Y二维扫描平台放置在耦合槽中,所述耦合槽中充满耦合液;A photoacoustic quantitative elastography device, the device includes a photoacoustic excitation source, a signal acquisition/transmission/reconstruction component, a coupling groove, a stepper motor and an X-Y two-dimensional scanning platform, and the photoacoustic excitation source includes a laser and a focusing lens; The signal acquisition/transmission/reconstruction assembly includes an ultrasonic detector, an amplifier, an oscilloscope and a computer, the ultrasonic detector, the amplifier, an oscilloscope and the computer are connected in sequence, and the computer is equipped with an acquisition control and signal processing system; the stepping motor Connected to a computer, the X-Y two-dimensional scanning platform is placed in a coupling tank, and the coupling tank is filled with coupling liquid;
测试时,组织样品放置在X-Y二维扫描平台上,并浸没在耦合槽的耦合液中;所述聚焦透镜设置在组织样品正上方,且聚焦透镜的焦点不偏离组织样品表面,所述激光器发出的脉冲激光通过聚焦透镜聚焦,照在组织样品上;所述超声探测器对准组织样品,且超声探测器的下端进入耦合槽的耦合液中,接收组织样品被激发出的光声信号;所述计算机控制步进电机逐点移动组织样品。During the test, the tissue sample is placed on the X-Y two-dimensional scanning platform and immersed in the coupling liquid in the coupling tank; the focusing lens is set directly above the tissue sample, and the focus of the focusing lens does not deviate from the surface of the tissue sample, and the laser emits The pulsed laser light is focused through the focusing lens and irradiated on the tissue sample; the ultrasonic probe is aimed at the tissue sample, and the lower end of the ultrasonic probe enters the coupling liquid in the coupling groove to receive the photoacoustic signal excited by the tissue sample; The computer controls the stepper motor to move the tissue sample point by point.
作为一种优选方案,所述装置还包括仪器固定/支撑器械组件,所述仪器固定/支撑器械组件用于固定/支撑X-Y二维扫描平台、聚焦透镜以及超声探测器。As a preferred solution, the device further includes an instrument fixing/supporting instrument assembly, and the instrument fixing/supporting instrument assembly is used for fixing/supporting the X-Y two-dimensional scanning platform, the focusing lens and the ultrasonic probe.
作为一种优选方案,所述超声探测器为水听器,其响应频率为200KHz-15MHz,直径为1mm;超声探测器的压电转换部件是一厚度为28μm的金电极聚偏二氟乙烯膜,接收光声信号时,该金电极聚偏二氟乙烯膜对准组织样品。As a preferred solution, the ultrasonic detector is a hydrophone with a response frequency of 200KHz-15MHz and a diameter of 1 mm; the piezoelectric conversion part of the ultrasonic detector is a gold electrode polyvinylidene fluoride film with a thickness of 28 μm , when receiving a photoacoustic signal, the polyvinylidene fluoride membrane of the gold electrode is aligned with the tissue sample.
作为一种优选方案,所述示波器的采样率为2.5GHz,所述计算机安装的采集控制及信号处理系统利用Labview和Matlab程序编写而成。As a preferred solution, the sampling rate of the oscilloscope is 2.5GHz, and the acquisition control and signal processing system installed in the computer is programmed by using Labview and Matlab programs.
作为一种优选方案,所述激光器发出的脉冲激光波长为400~2500nm,脉冲宽度为1~50ns,重复频率为1Hz~5KHz。As a preferred solution, the wavelength of the pulse laser emitted by the laser is 400-2500 nm, the pulse width is 1-50 ns, and the repetition frequency is 1 Hz-5 KHz.
作为一种优选方案,所述耦合槽中的耦合液为水,监测水温,使水温与组织样品的温度保持一致。As a preferred solution, the coupling liquid in the coupling tank is water, and the water temperature is monitored to keep the water temperature consistent with the temperature of the tissue sample.
本发明相对于现有技术具有如下的有益效果:Compared with the prior art, the present invention has the following beneficial effects:
1、本发明将示波器采集的信号对时间进行一次积分,得到组织样品表面振动位移的时间函数,获取组织样品表面振动位移从零上升至其最大值时所需时间,利用该时间计算组织定量弹性模量,与已有的相对弹性测量的方法相比,无需正常组织作为参考,具有更高的准确性。1. In the present invention, the signal collected by the oscilloscope is integrated once over time to obtain the time function of the surface vibration displacement of the tissue sample, and the time required for the surface vibration displacement of the tissue sample to rise from zero to its maximum value is obtained, and the quantitative elasticity of the tissue is calculated using this time Modulus, compared with existing relative elastic measurement methods, does not require normal tissue as a reference and has higher accuracy.
2、本发明利用激光器发出脉冲激光,该脉冲激光通过聚焦透镜聚焦,照在组织样品上,激发出光声信号,从而进行进行组织弹性检测,与传统的超声弹性检测方法相比,具有组织特异性和高分辨能力。2. The present invention uses a laser to emit pulsed laser light. The pulsed laser light is focused by a focusing lens and irradiated on a tissue sample to excite a photoacoustic signal, thereby performing tissue elasticity detection. Compared with the traditional ultrasonic elasticity detection method, it has tissue specificity. and high resolution capabilities.
3、本发明采用的超声探测器为水听器,其响应频率为200KHz-15MHz,直径为1mm,具有探测灵敏度高、无带宽限制的优点,从而保证了高灵敏检测的能力。3. The ultrasonic detector used in the present invention is a hydrophone with a response frequency of 200KHz-15MHz and a diameter of 1mm. It has the advantages of high detection sensitivity and no bandwidth limitation, thereby ensuring the ability of highly sensitive detection.
4、本发明的光声定量弹性检测的方法具备快速检测的能力,实现该方法的装置结构简单、使用方便,可以广泛应用于组织的弹性成像中,便于产业化。4. The photoacoustic quantitative elastic detection method of the present invention has the ability of rapid detection, and the device for realizing the method is simple in structure and easy to use, and can be widely used in tissue elastic imaging, which is convenient for industrialization.
附图说明Description of drawings
图1为本发明实施例1的光声定量弹性成像装置的结构示意图。FIG. 1 is a schematic structural diagram of a photoacoustic quantitative elastography device according to Embodiment 1 of the present invention.
图2为本发明实施例1的位移对时间的依赖关系曲线图。FIG. 2 is a graph showing the dependence of displacement on time in Example 1 of the present invention.
图3为本发明实施例2的琼脂样品a的示意图。Fig. 3 is a schematic diagram of agar sample a of Example 2 of the present invention.
图4为本发明实施例2的琼脂样品b的示意图。Fig. 4 is a schematic diagram of agar sample b of Example 2 of the present invention.
图5为本发明实施例2的琼脂样品a的光声图像。Fig. 5 is a photoacoustic image of agar sample a in Example 2 of the present invention.
图6为本发明实施例2的琼脂样品b的光声图像。Fig. 6 is a photoacoustic image of agar sample b in Example 2 of the present invention.
图7为本发明实施例2的琼脂样品a的光声弹性图像。Fig. 7 is a photoacoustic elastic image of agar sample a in Example 2 of the present invention.
图8为本发明实施例2的琼脂样品b的光声弹性图像。Fig. 8 is a photoacoustic elastic image of agar sample b in Example 2 of the present invention.
图9为图3和图4虚线处的弹性模量曲线图。FIG. 9 is a curve diagram of elastic modulus at the dotted line in FIG. 3 and FIG. 4 .
其中,1-耦合槽,2-步进电机,3-激光器,4-聚焦透镜,5-超声探测器,6-放大器,7-示波器,8-计算机,9-组织样品。Among them, 1-coupling tank, 2-stepping motor, 3-laser, 4-focusing lens, 5-ultrasonic detector, 6-amplifier, 7-oscilloscope, 8-computer, 9-tissue sample.
具体实施方式Detailed ways
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。The present invention will be further described in detail below in conjunction with the embodiments and the accompanying drawings, but the embodiments of the present invention are not limited thereto.
实施例1:Example 1:
如图1所示,本实施例的光声定量弹性成像装置包括光声激发源、信号采集/传输/重建组件、耦合槽1、步进电机2、X-Y二维扫描平台和仪器固定/支撑器械组件(图中未示出),所述光声激发源包括激光器3和聚焦透镜4;所述信号采集/传输/重建组件包括超声探测器5、放大器6、示波器7和计算机8,所述超声探测器5、放大器6、示波器7和计算机8依次相连;所述示波器7的采样率为2.5GHz;所述计算机8安装有采集控制及信号处理系统,该系统利用Labview和Matlab程序编写而成;所述步进电机2与计算机8相连,所述X-Y二维扫描平台放置在耦合槽1中,所述耦合槽1中充满耦合液;As shown in Figure 1, the photoacoustic quantitative elastography device of this embodiment includes a photoacoustic excitation source, a signal acquisition/transmission/reconstruction component, a coupling tank 1, a stepping motor 2, an X-Y two-dimensional scanning platform and an instrument fixing/supporting device Components (not shown in the figure), the photoacoustic excitation source includes a laser 3 and a focusing lens 4; the signal acquisition/transmission/reconstruction assembly includes an ultrasonic probe 5, an amplifier 6, an oscilloscope 7 and a computer 8, and the ultrasonic Detector 5, amplifier 6, oscilloscope 7 and computer 8 are connected successively; The sampling rate of described oscilloscope 7 is 2.5GHz; Described computer 8 is equipped with acquisition control and signal processing system, and this system utilizes Labview and Matlab program to form; The stepper motor 2 is connected to the computer 8, the X-Y two-dimensional scanning platform is placed in the coupling tank 1, and the coupling tank 1 is filled with coupling liquid;
所述超声探测器5为水听器,其响应频率为200KHz-15MHz,直径为1mm;超声探测器的压电转换部件是一厚度为28μm的金电极聚偏二氟乙烯膜。The ultrasonic probe 5 is a hydrophone with a response frequency of 200 KHz-15 MHz and a diameter of 1 mm; the piezoelectric conversion part of the ultrasonic probe is a gold electrode polyvinylidene fluoride film with a thickness of 28 μm.
所述仪器固定/支撑器械组件用于固定/支撑X-Y二维扫描平台、聚焦透镜4以及超声探测器5。The instrument fixing/supporting instrument assembly is used to fix/support the X-Y two-dimensional scanning platform, the focusing lens 4 and the ultrasonic probe 5 .
实现本实施例的光声定量弹性成像方法,主要采用以下原理:To realize the photoacoustic quantitative elastography method of this embodiment, the following principles are mainly adopted:
组织被激光照射产生剪切波,其方程为:The tissue is irradiated by laser to generate shear wave, whose equation is:
其中,uz为在z方向即激发光线方向的位移,是剪切波位移也是组织被激励后产生光声信号时的位移。为剪切波波速,μ为剪切模量,v=μ/ρ为动态剪切粘性系数,Fz为光致热弹辐射力。在激光焦点位置,M=β0αI0,β0为热膨胀系数,α为组织的吸收系数,I0为激光强度。r为位置自变量,R为激光光斑半径,φ(t)为激光强度随时间变化关系,脉冲宽度为t0。Among them, u z is the displacement in the z direction, that is, the direction of the excitation light, which is the shear wave displacement and also the displacement when the photoacoustic signal is generated after the tissue is excited. is the shear wave velocity, μ is the shear modulus, v=μ/ρ is the dynamic shear viscosity coefficient, and F z is the photothermoelastic radiation force. At the laser focus position, M=β 0 αI 0 , β 0 is the thermal expansion coefficient, α is the tissue absorption coefficient, and I 0 is the laser intensity. r is the independent variable of the position, R is the radius of the laser spot, φ(t) is the relationship of the laser intensity with time, and the pulse width is t 0 .
按照常微分方程的解法,并进行汉克尔变换得到:According to the solution of the ordinary differential equation, and carry out the Hankel transformation to get:
上面的方程写成格林函数的形式:The above equation is written in the form of Green's function:
此剪切位移同时是产生光声信号时生物组织表面的位移。This shear displacement is simultaneously the displacement of the surface of the biological tissue when the photoacoustic signal is generated.
因为光声信号和位移之间满足:Because the relationship between the photoacoustic signal and the displacement satisfies:
考虑在激光焦点处温度变化为时间慢变化,通过对上式进行分析变量:Considering that the temperature change at the laser focus is a time-slow change, the variable is analyzed by the above formula:
光声信号幅值从零变化到其最大值所经历时间满足tmax=R/ct,剪切模量μ满足μ=ρ(R/tmax)2;其中,ρ为生物组织密度,可取为1100kg/m3,通过利用弹性模量和剪切模量关系E=2μ(1+η),其中,η为生物组织泊松比,由于生物组织不可压缩,取为0.499,从而获得组织的定量弹性模量大小:The time elapsed for the amplitude of the photoacoustic signal to change from zero to its maximum value satisfies t max =R/c t , and the shear modulus μ satisfies μ=ρ(R/t max ) 2 ; where ρ is the density of biological tissue, which can be is 1100kg/m 3 , by using the relationship between elastic modulus and shear modulus E=2μ(1+η), where η is the Poisson’s ratio of biological tissue, which is taken as 0.499 due to the incompressibility of biological tissue, thereby obtaining the Quantitative modulus of elasticity:
因此,本实施例的光声定量弹性成像方法包括以下步骤:Therefore, the photoacoustic quantitative elastography method of this embodiment includes the following steps:
1)将组织样品9放置在X-Y二维扫描平台上,并浸没在耦合槽1的耦合液中;将聚焦透镜4设置在组织样品9正上方,并调整聚焦透镜4的高度使聚焦透镜4的焦点不偏离组织样品9表面;将超声探测器5对准组织样品9,即金电极聚偏二氟乙烯膜对准组织样品9,并使超声探测器5的下端进入耦合槽1的耦合液深度约5-8mm处中,耦合液为水,监测水温,使水温与组织样品9的温度保持一致;1) The tissue sample 9 is placed on the X-Y two-dimensional scanning platform, and immersed in the coupling liquid of the coupling tank 1; the focusing lens 4 is arranged directly above the tissue sample 9, and the height of the focusing lens 4 is adjusted so that the focusing lens 4 The focus does not deviate from the surface of the tissue sample 9; the ultrasonic probe 5 is aligned with the tissue sample 9, that is, the polyvinylidene fluoride film of the gold electrode is aligned with the tissue sample 9, and the depth of the coupling liquid that makes the lower end of the ultrasonic probe 5 enter the coupling groove 1 is about In the 5-8mm position, the coupling liquid is water, and the water temperature is monitored to keep the water temperature consistent with the temperature of the tissue sample 9;
2)激光器3发出脉冲激光,发出的脉冲激光波长为400~2500nm,脉冲宽度为1~50ns,重复频率为1Hz~5KHz,该脉冲激光通过聚焦透镜4聚焦,照在组织样品9上,激发出光声信号,光声信号经过耦合槽1的耦合液后被超声探测器5接收;2) The laser 3 emits pulsed laser light with a wavelength of 400-2500 nm, a pulse width of 1-50 ns, and a repetition rate of 1 Hz-5 KHz. The pulsed laser light is focused by the focusing lens 4 and irradiated on the tissue sample 9 to excite light The acoustic signal, the photoacoustic signal is received by the ultrasonic detector 5 after passing through the coupling liquid in the coupling tank 1;
3)超声探测器5接收的光声信号经放大器6放大后,被示波器7采集,示波器7将采集的信号信息存储到计算机8中,计算机8利用Labview程序控制步进电机2逐点移动组织样品9,对应的组织样品9上形成X-Y二维平面扫描区域,步进电机2每移动一次,示波器7就进行一次信号采集;3) The photoacoustic signal received by the ultrasonic detector 5 is amplified by the amplifier 6 and collected by the oscilloscope 7. The oscilloscope 7 stores the collected signal information in the computer 8, and the computer 8 uses the Labview program to control the stepping motor 2 to move the tissue sample point by point 9. An X-Y two-dimensional plane scanning area is formed on the corresponding tissue sample 9, and each time the stepping motor 2 moves once, the oscilloscope 7 performs a signal acquisition;
4)示波器7采集完全部信号后,计算机8利用Matlab程序将各点的信号对时间进行一次积分,得到组织样品表面振动位移的时间函数;获取组织样品表面振动位移从零上升至其最大值时所需时间,如图2所示,图中给出了四种位移随着时间的变化关系,不同的弹性模量,位移到达最大值所需的时间是不相同的;通过上述式(6),利用该时间计算各点的组织样品定量弹性模量;根据计算的定量弹性模量,重建出组织样品的定量弹性二维图像。4) After the oscilloscope 7 has collected all the signals, the computer 8 uses the Matlab program to integrate the signals of each point with respect to time to obtain the time function of the vibration displacement of the tissue sample surface; when the vibration displacement of the tissue sample surface rises from zero to its maximum value, The required time, as shown in Figure 2, provides four kinds of displacement relationships with time in the figure, different elastic modulus, the time required for displacement to reach the maximum value is not the same; through the above formula (6) , using the time to calculate the quantitative elastic modulus of the tissue sample at each point; according to the calculated quantitative elastic modulus, the quantitative elastic two-dimensional image of the tissue sample is reconstructed.
实施例2:Example 2:
本实施例是利用琼脂样品进行的实验,主要包括以下步骤:The present embodiment is the experiment that utilizes agar sample to carry out, mainly comprises the following steps:
1)在浓度为20g/L的琼脂中加入10%的墨水做成方形样品,在它中间做成浓度为15g/L的琼脂并加入3%的墨水做成方圆形样品,这形成了图3中的琼脂样品a;在浓度为30g/L的琼脂中加入10%的墨水做成方形样品,在它中间做成浓度为25g/L的琼 脂并加入3%的墨水做成方圆形样品,这形成了图4中的琼脂样品b;1) Add 10% ink to make a square sample in the agar with a concentration of 20g/L, make agar with a concentration of 15g/L in the middle of it and add 3% ink to make a square sample, which forms a graph Agar sample a in 3; add 10% ink to agar with a concentration of 30g/L to make a square sample, make agar with a concentration of 25g/L in the middle of it and add 3% ink to make a square sample , which forms agar sample b in Figure 4;
2)启动激光器,输出脉冲激光波长为532nm,脉宽为10ns,重复频率是15Hz;该脉冲激光通过聚焦透镜聚焦后照射在琼脂样品a和b上,琼脂样品a和b被激发出光声信号,光声信号经过耦合槽中的耦合液后被超声探测器接收;2) start the laser, the output pulse laser has a wavelength of 532nm, a pulse width of 10ns, and a repetition rate of 15Hz; the pulse laser is irradiated on agar samples a and b after being focused by a focusing lens, and agar samples a and b are excited to emit photoacoustic signals, The photoacoustic signal is received by the ultrasonic detector after passing through the coupling liquid in the coupling tank;
3)超声探测器接收的光声信号经放大器放大后,将其传输到示波器进行数据采集,示波器再将数据传输并存储到计算机中,计算机控制步进电机对应的在样品上形成扫描区域;3) After the photoacoustic signal received by the ultrasonic detector is amplified by the amplifier, it is transmitted to the oscilloscope for data collection, and the oscilloscope then transmits and stores the data to the computer, and the computer controls the stepping motor to form a corresponding scanning area on the sample;
4)采集完全部信号后,对采集的数据先归一化,然后用最大值投影法和弹性投影法重建光声图像和光声弹性图像,图5和图6分别是琼脂样品a和b的光声图像,可以看出,显示的背景琼脂光声成像几乎没有对比度,中间圆形区域的光声成像也几乎没有对比度;图7和图8分别是琼脂样品a和b的光声弹性图像,可以看出,显示的光声弹性图像则显示大的对比度,这说明了光声弹性成像的必要性;图9中的a和b分别显示的是图3和图4虚线处琼脂样品的弹性值,可以看出弹性值的明显变化。4) After all the signals are collected, the collected data are first normalized, and then the photoacoustic image and the photoacoustic image are reconstructed by the maximum projection method and the elastic projection method. Figure 5 and Figure 6 are the photoacoustic images of agar samples a and b respectively. Acoustic image, it can be seen that the photoacoustic imaging of the displayed background agar has almost no contrast, and the photoacoustic imaging of the middle circular area has almost no contrast; Figure 7 and Figure 8 are the photoacoustic elastic images of agar samples a and b respectively, which can be It can be seen that the displayed photoacoustic elastography image shows a large contrast, which illustrates the necessity of photoacoustic elastography; a and b in Fig. 9 show the elasticity values of the agar samples at the dotted lines in Fig. 3 and Fig. 4 respectively, A clear change in the elasticity value can be seen.
综上所述,本发明提出了一种新的组织定量弹性模量测量方法,与已有的相对弹性测量的方法相比,无需正常组织作为参考,具有更高的准确性;利用聚焦激光进行组织弹性检测,与传统的超声弹性检测方法相比,具有组织特异性和高分辨能力;采用的超声探测器为水听器,具有探测灵敏度高、无带宽限制的优点,从而保证了高灵敏检测的能力。In summary, the present invention proposes a new method for measuring the quantitative elastic modulus of tissue, compared with the existing method for measuring relative elasticity, it does not require normal tissue as a reference, and has higher accuracy; Tissue elasticity detection, compared with the traditional ultrasonic elasticity detection method, has tissue specificity and high resolution; the ultrasonic detector used is a hydrophone, which has the advantages of high detection sensitivity and no bandwidth limitation, thus ensuring high sensitivity detection Ability.
以上所述,仅为本发明专利较佳的实施例,但本发明专利的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明专利所公开的范围内,根据本发明专利的技术方案及其发明构思加以等同替换或改变,都属于本发明专利的保护范围。The above is only a preferred embodiment of the patent of the present invention, but the scope of protection of the patent of the present invention is not limited thereto. Equivalent replacements or changes to the technical solutions and their inventive concepts all fall within the scope of protection of the invention patent.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510944323.0A CN105572049B (en) | 2015-12-16 | 2015-12-16 | Optoacoustic quantifies elastograph imaging method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510944323.0A CN105572049B (en) | 2015-12-16 | 2015-12-16 | Optoacoustic quantifies elastograph imaging method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105572049A CN105572049A (en) | 2016-05-11 |
CN105572049B true CN105572049B (en) | 2018-09-11 |
Family
ID=55882429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510944323.0A Active CN105572049B (en) | 2015-12-16 | 2015-12-16 | Optoacoustic quantifies elastograph imaging method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105572049B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106419803B (en) * | 2016-08-31 | 2018-06-19 | 华南师范大学 | Intravascular photoacoustic absorption, elasticity, sticky multi-modality imaging integration endoscope and its imaging method |
CN106501182B (en) * | 2016-09-22 | 2019-06-04 | 南京大学 | A method for non-destructive measurement of elasticity using photoacoustic eigenspectroscopy |
CN106769877A (en) * | 2016-11-29 | 2017-05-31 | 哈尔滨工业大学 | A kind of biological tissue's frequency domain photoacoustic imaging detection method and system |
CN106680214B (en) * | 2017-03-20 | 2019-08-27 | 华南师范大学 | A photoacoustic laparoscopic device and method for simultaneously detecting elasticity and viscosity |
CN107132187B (en) * | 2017-04-05 | 2020-08-28 | 清华大学 | Photoacoustic imaging system and imaging method |
CN109781241B (en) * | 2019-02-01 | 2020-11-03 | 华南师范大学 | A device and method for measuring ultrasonic field distribution based on photoacoustic effect |
CN109870413A (en) * | 2019-03-05 | 2019-06-11 | 上海市同济医院 | Experimental device and method for photoacoustic detection of varicocele model in adolescent rats |
CN113777045B (en) * | 2020-06-10 | 2022-10-18 | 复旦大学 | Super-resolution functional photoacoustic imaging method based on single-particle multilateral localization tracking |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1382983A (en) * | 2002-03-21 | 2002-12-04 | 江苏大学 | Nondestructive quantitative testing method and device for photoacoustic performance of film-based coupling structure of micro-opto-electro-mechanical system |
CN102175776A (en) * | 2011-01-14 | 2011-09-07 | 华南师范大学 | Photoacoustic elastic imaging method and device |
CN102297898A (en) * | 2011-05-17 | 2011-12-28 | 南京理工大学 | Laser ultrasonic measuring method for third order elastic constant of metal |
CN102879335A (en) * | 2012-09-26 | 2013-01-16 | 华南师范大学 | Portable noninvasive real-time photoacoustic viscoelastic detector |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004351153A (en) * | 2003-05-30 | 2004-12-16 | Olympus Corp | Optical imaging apparatus |
-
2015
- 2015-12-16 CN CN201510944323.0A patent/CN105572049B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1382983A (en) * | 2002-03-21 | 2002-12-04 | 江苏大学 | Nondestructive quantitative testing method and device for photoacoustic performance of film-based coupling structure of micro-opto-electro-mechanical system |
CN102175776A (en) * | 2011-01-14 | 2011-09-07 | 华南师范大学 | Photoacoustic elastic imaging method and device |
CN102297898A (en) * | 2011-05-17 | 2011-12-28 | 南京理工大学 | Laser ultrasonic measuring method for third order elastic constant of metal |
CN102879335A (en) * | 2012-09-26 | 2013-01-16 | 华南师范大学 | Portable noninvasive real-time photoacoustic viscoelastic detector |
Also Published As
Publication number | Publication date |
---|---|
CN105572049A (en) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105572049B (en) | Optoacoustic quantifies elastograph imaging method and device | |
JP5235586B2 (en) | Biological information processing apparatus and biological information processing method | |
CN102175776B (en) | Photoacoustic elastic imaging method and device | |
JP5201920B2 (en) | measuring device | |
Vaithilingam et al. | Three-dimensional photoacoustic imaging using a two-dimensional CMUT array | |
CN102579080B (en) | Integrated portable confocal opto-acoustic microscopy imaging device and method | |
JP5538856B2 (en) | Photoacoustic device | |
CN101813672B (en) | Rapid three-dimensional photoacoustic imaging system based on ultrasonic plane array detector and method thereof | |
CN102944521B (en) | Non-contact photoacoustic and optical coherence tomography dual-imaging device and detection method thereof | |
CN106560160A (en) | Object Information Acquiring Apparatus And Control Method Thereof | |
JP2010088873A (en) | Biological information imaging apparatus | |
CN110455719A (en) | Three-dimensional photoacoustic imaging system and method | |
CN109946236A (en) | A polarized photoacoustic microscopy imaging device and imaging method | |
CN104323762A (en) | Photoacoustic microscopy imaging-based quantitative detection device for nevus flammeus blood vessel | |
Rebling et al. | Optoacoustic characterization of broadband directivity patterns of capacitive micromachined ultrasonic transducers | |
CN106725348A (en) | A kind of absorber absorption coefficient measurement and the simultaneously lossless detection method and its device of photoacoustic imaging | |
CN106680214B (en) | A photoacoustic laparoscopic device and method for simultaneously detecting elasticity and viscosity | |
CN111528921B (en) | Information acquisition device | |
CN105572048B (en) | The quantitative elasticity of optoacoustic, the two-parameter imaging method of viscosity and device | |
JP5183406B2 (en) | Biological information processing apparatus and biological information processing method | |
Kim et al. | Measurements of acoustic radiation force of ultrahigh frequency ultrasonic transducers using model-based approach | |
CN1168980C (en) | Method and device for measuring photoacoustic signals in biological tissue using probe ultrasound beam | |
CN104161520A (en) | Epidermal melanin concentration determination method and device based on photoacoustic effect principle | |
CN110261315A (en) | A kind of scanning near-field opto-acoustic microscopic imaging instrument | |
CN106769877A (en) | A kind of biological tissue's frequency domain photoacoustic imaging detection method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |