CN113812926A - 一种基于激光多普勒测振的磁声耦合成像***及方法 - Google Patents
一种基于激光多普勒测振的磁声耦合成像***及方法 Download PDFInfo
- Publication number
- CN113812926A CN113812926A CN202111137246.XA CN202111137246A CN113812926A CN 113812926 A CN113812926 A CN 113812926A CN 202111137246 A CN202111137246 A CN 202111137246A CN 113812926 A CN113812926 A CN 113812926A
- Authority
- CN
- China
- Prior art keywords
- target body
- vibration
- laser doppler
- frequency
- excitation
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 57
- 238000005259 measurement Methods 0.000 title claims abstract description 32
- 230000008878 coupling Effects 0.000 title claims abstract description 29
- 238000010168 coupling process Methods 0.000 title claims abstract description 29
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000005284 excitation Effects 0.000 claims abstract description 63
- 230000003068 static effect Effects 0.000 claims abstract description 23
- 230000009471 action Effects 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 10
- 230000001939 inductive effect Effects 0.000 claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims description 39
- 239000013598 vector Substances 0.000 claims description 23
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000004907 flux Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000003321 amplification Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 2
- 239000000523 sample Substances 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000002245 particle Substances 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005670 electromagnetic radiation Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000035559 beat frequency Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011895 specific detection Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于激光多普勒测振的磁声耦合成像***及方法，***包括：激励源、激励线圈、磁铁、激光多普勒振动传感器、导轨和控制模块；激励源产生激励电流至激励线圈；激励线圈位于目标体的上方位置，其用于在激励电流的作用下，引发目标体中产生感应涡流；磁铁分别位于目标体的上下方，产生静磁场；目标体在感应涡流和静磁场的共同作用下，产生与激励源同频的振动；导轨围绕目标体设置；激光多普勒振动传感器在导轨上移动并360度采集目标体的振动信号；控制模块用于对振动信号进行处理，并根据振动信号，对目标体内部的电学参数分布进行三维重建。本发明通过检测目标体表面高频振动信息，改善成像方法的信噪比，提升成像空间分辨率。
Description
技术领域
本发明涉及生物医学成像技术领域，更具体的说是涉及一种基于激光多普勒测振的磁声耦合成像***及方法。
背景技术
磁声耦合成像(Magneto-acoustic Tomography,MAT)是一种无损生物组织电特性功能成像技术。该方法将目标体放在静磁场和脉冲磁场中，脉冲磁场在目标体中产生涡流，涡流在静磁场作用下生成洛伦兹力，引发目标体发生超声振动，通过对声信号的采集，MAT能够重建出目标体内的电导率参数空间分布。MAT技术能够先于生物组织结构变化，提供更有效的组织功能信息用于早期疾病诊断。
虽然MAT技术已在成像理论和技术上取得了不断发展，但仍存在许多问题亟待研究。当前，MAT***中普遍采用基于压电效应的超声探头对信号进行检测，此类压电超声探头技术成熟，但是，应用于MAT成像存在以下几方面的问题：
(1)MAT成像中超声信号微弱，声场信号检测方式具有较低信噪比。例如,0.1T磁场强度下，1μA mm偶极矩受到的洛伦兹力仅为10-10N量级；在1T强磁场中，直径10mm的弹性导电球体，探测位置为1mm远时，能够接收到的超声压强为1/6N/m2(1/6×10-5atm)，而且，声压幅值还会随探测距离增加而不断减小；在感应式MAT成像中，1T磁场强度、1000V/m电场和0.5μs脉宽激励信号条件下，经过计算得知，直径10mm弹性导电球体，探测距离为50mm时，能够探测到的超声压强仅为0.015Pa。
(2)MAT成像使用的超声探头中心频率较低，多为10MHz以下，限制了MAT的空间分辨率。根据MAT成像理论，输出超声信号与激励信号同频，激励信号频率是决定MAT空间分辨率的关键因素:生物组织中声速约为1500m/s，MAT实验中经常使用1MHz激励脉冲信号(脉宽为1μs)，因此，理论上MAT成像的空间分辨率约为1500m/s×1μs＝1.5mm。由于MAT中采用的超声探头中心频率较低，因此MAT成像空间分辨率被限制在mm量级。同时，超声探头测量频带较窄，检测时会错过一定量带有组织电特性的其他频率信息，进一步影响了成像重建效果。
(3)强电磁辐射对超声探头的干扰。MAT成像***普遍采用0.1T以上的静磁场和1MHz高频脉冲磁场，MAT成像的空间环境复杂，电磁辐射剧烈变化、干扰大。基于压电效应的超声探头属于电磁敏感设备，为了改善信噪比，MAT需要多次(1000次)测量后取平均信号波形作为最终数据，导致测量时间较长。
因此，如何克服目前磁声耦合成像方法中信号采集方面存在的局限性问题，是本领域技术人员亟需解决的问题。
发明内容
有鉴于此，本发明提供了一种基于激光多普勒测振的磁声耦合成像***，通过检测目标体表面高频(10MHz量级)振动信息(如位移、速度和加速度)，避免遭受传统磁声耦合成像方法中的电磁干扰，改善成像方法的信噪比，提升成像空间分辨率，凭借较宽的频率带宽来获得更丰富的组织电特性信息。
为了实现上述目的，本发明采用如下技术方案：
一种基于激光多普勒测振的磁声耦合成像***，包括：激励源、激励线圈、磁铁、激光多普勒振动传感器、导轨和控制模块；
其中，所述激励源在所述控制模块的控制下产生激励电流至所述激励线圈；所述激励线圈位于目标体的上方位置，其用于在所述激励电流的作用下，引发所述目标体中产生感应涡流；所述磁铁设置有一对，且分别位于所述目标体的上下方，其用于产生静磁场；所述目标体在所述感应涡流和所述静磁场的共同作用下，产生与所述激励源同频的振动；
所述导轨围绕所述目标体设置；所述激光多普勒振动传感器安装在所述导轨上，且在所述控制模块的控制下沿所述导轨移动并360度采集所述目标体的振动信号；
所述控制模块用于对所述振动信号进行处理，并根据所述振动信号，对所述目标体内部的电学参数分布进行三维重建。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像***中，所述控制模块包括滤波放大器、数据采集卡、主控单元和计算机；所述滤波放大器、所述数据采集卡和所述计算机依次与所述激光多普勒振动传感器电性连接；所述主控单元分别与所述激励源、所述激光多普勒振动传感器和所述计算机电性连接；
所述滤波放大器用于对所述振动信号进行滤波、除噪和幅值放大处理；
所述数据采集卡用于将处理后的所述振动信号收集并上传至所述计算机；
所述计算机用于根据处理后的所述振动信号对所述目标体内部的电学参数分布进行三维重建，并通过所述主控单元对所述激励源和所述激光多普勒振动传感器进行参数设定和运行管理。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像***中，所述激光多普勒振动传感器包括激光器、光分束器、反射光路、参考光路、光学探测器、速度解调器和伺服分析模块；
所述激光器用于发射激光束；所述光分束器将所述激光束分成两束，其中一束经所述反射光路照射至所述目标体上，生成反射光并返回至所述光学探测器；另一束作为参考光，并经所述参考光路入射至所述光学探测器；所述光学探测器用于将所述反射光和所述参考光混频成脉冲光，并对脉冲光进行光电转换，生成脉冲电信号；所述速度解调器用于对所述脉冲电信号进行解调，生成速度信号；所述伺服分析模块用于根据积分演算方法，将所述速度信号转换成位移信号。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像***中，所述反射光路包括第一反光镜、第二反光镜、第一半透半反镜、透镜和第三反光镜；
经所述光分束器分光后的其中一束激光束依次通过所述第一反光镜和所述第二反光镜反射后，通过所述第一半透半反镜和所述透镜后，照射至处于振动状态的所述目标体上，在多普勒效应下产生频率发生偏移的反射光；
生成的反射光经所述透镜汇聚光线后，再依次经所述第一半透半反镜和所述第三反光镜调整光轴后，照射至所述光学探测器中。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像***中，所述参考光路包括声光调制器和第二半透半反镜；所述参考光经所述声光调制器驱动，发生偏移，并经所述第二半透半反镜调整光轴后，照射至所述光学探测器中。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像***中，所述声光调制器为AOM布拉格盒。
本发明还提供一种基于激光多普勒测振的磁声耦合成像方法，其适用于上述基于激光多普勒测振的磁声耦合成像***，包括以下步骤：
利用所述控制模块对所述激励源和所述激光多普勒振动传感器进行参数设置；
利用所述磁铁产生磁通密度向量为B0的静磁场，使目标体置于静磁场空间内；
启动所述激励源，产生高频激励电流至所述激励线圈；
利用所述激励线圈产生的环形激励电流，在所述目标体中引发感应电流密度向量为J的感生涡流；
感生涡流在静磁场作用下，产生洛伦兹力F＝J×B0；在洛伦兹力F作用下，所述目标体内部产生与所述激励源同频的高频振动；
在所述控制模块的控制下，所述激光多普勒振动传感器沿所述轨道移动，对目标体进行360°扫描探测，得到所述目标体的振动信号；
利用所述控制模块对所述振动信号进行处理，并根据所述振动信号，对所述目标体内部的电学参数分布进行三维重建。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像方法中，所述对目标体进行360°扫描探测，得到所述目标体的振动信号，包括：
利用所述激光多普勒振动传感器发射频率为f0的激光束，并将所述激光束分束成两束，其中一束照射至处于振动状态的所述目标体上，在多普勒效应下，生成频率为fr的反射光；另一束作为参考光；
将参考光在驱动频率fB下进行调制，使其发生频移；
将反射光与频移后的参照光在同一空间位置上发生干涉，混频成频率大小f的脉冲光；
对所述脉冲光进行光电转换和频率解调，得到多普勒频率fD大小和与之成比例的目标体的振动速度；计算公式如下：
fD＝f-fB；
v＝λfD/2；
其中，v表示目标体的振动速度，λ表示波长；
利用下式计算所述目标体的振幅A：
A＝v/ω；
其中，ω表示目标体以角频率ω作简谐振动。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像方法中，所述目标体内洛伦兹力引发的表面振动位移的关系模型为：
其中，ρ为声波引起的密度变化，Cs为目标体内声速；∮V dr'为体积分符号；J(r′)为声源处的电流密度；B0(r′)为声源处的磁场强度；t为时间；r′为声源的位置矢量；rb为探测器在目标体表面的位置矢量。
优选的，在上述一种基于激光多普勒测振的磁声耦合成像方法中，所述目标体内部的电学参数的计算公式为：
其中，σ表示所述目标体内部组织电导率，J表示感应电流密度向量，B0表示静磁场磁通密度向量，B1表示J引起的交变磁场强度。
经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种基于激光多普勒测振的磁声耦合成像***及方法，能够克服目前磁声耦合成像方法中信号采集方面存在的局限性，通过检测目标体表面高频(10MHz量级)振动信息(如位移、速度和加速度)，避免遭受传统磁声耦合成像方法中的电磁干扰，改善成像方法的信噪比，提升成像空间分辨率，凭借较宽的频率带宽来获得更丰富的组织电特性信息。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的基于激光多普勒测振的磁声耦合成像***的结构示意图；
图2附图为本发明提供的激光多普勒振动传感器的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例公开了一种基于激光多普勒测振的磁声耦合成像***，包括：激励源1、激励线圈2、磁铁3、激光多普勒振动传感器4、导轨5和控制模块；
其中，激励源1在控制模块的控制下产生激励电流至激励线圈2；激励线圈2位于目标体的上方位置，其用于在激励电流的作用下，引发目标体中产生感应涡流；磁铁3设置有一对，且分别位于目标体的上下方，其用于产生静磁场；目标体在感应涡流和静磁场的共同作用下，产生与激励源同频的振动；
导轨5围绕目标体设置；激光多普勒振动传感器4安装在导轨5上，且在控制模块的控制下沿导轨5移动并360度采集目标体的振动信号；
控制模块用于对振动信号进行处理，并根据振动信号，对目标体内部的电学参数分布进行三维重建。
具体的，控制模块包括滤波放大器6、数据采集卡7、主控单元8和计算机9；滤波放大器6、数据采集卡7和计算机9依次与激光多普勒振动传感器4电性连接；主控单元8分别与激励源1、激光多普勒振动传感器4和计算机9电性连接；
滤波放大器6用于对振动信号进行滤波、除噪和幅值放大处理；
数据采集卡7用于将处理后的振动信号收集并上传至计算机9；
计算机9用于根据处理后的振动信号对目标体内部的电学参数分布进行三维重建，并通过主控单元8对激励源1和激光多普勒振动传感器4进行参数设定和运行管理。
本发明实施例对目标体振动进行检测的过程如下：
选取高重复频率的大功率电流源作为激励源1；采用圆形激励线圈2置于目标体上，激励线圈2中脉冲激励信号在目标体中引发感应电流密度向量为J的感生涡流；一对永磁铁3分别放置在目标体上下两侧，从而生成磁通密度向量为B0的静磁场，目标体置于静磁场中；涡流在静磁场作用下产生洛伦兹力F＝J×B0，引发目标体内部发生与激励源同频的高频振动，振动由振源向目标体表面进行传播；激光多普勒振动传感器4向目标体发射探测激光，同时接收获取目标体表面反射回来的激光信号并进行光电检测变换；激光多普勒振动传感器4安装于导轨5，对目标体表面进行360°扫描测量；通过滤波放大器6对振动传感器4输出的微弱信号进行滤波除噪声、幅值放大预处理；经过处理的信号被数据采集卡7收集并上传至计算机9；计算机9通过主控单元8对激励源1和激光多普勒振动传感器4进行参数设定和运行管理。
在一个具体实施例中，如图2所示，激光多普勒振动传感器4包括激光器A、光分束器BS1、反射光路、参考光路、光学探测器B、速度解调器C和伺服分析模块D；
激光器A用于发射激光束；光分束器BS1将激光束分成两束，其中一束经反射光路照射至目标体上，生成反射光并返回至光学探测器B；另一束作为参考光，并经参考光路入射至光学探测器B；光学探测器B用于将反射光和参考光混频成脉冲光，并对脉冲光进行光电转换，生成脉冲电信号；速度解调器C用于对脉冲电信号进行解调，生成速度信号；伺服分析模块D用于根据积分演算方法，将速度信号转换成位移信号。
其中，反射光路包括第一反光镜M1、第二反光镜M2、第一半透半反镜BS2、透镜E和第三反光镜M3；
经光分束器BSA分光后的其中一束激光束依次通过第一反光镜M1和第二反光镜M2反射后，通过第一半透半反镜BS2和透镜E后，照射至处于振动状态的目标体上，在多普勒效应下产生频率发生偏移的反射光；
生成的反射光经透镜E汇聚光线后，再依次经第一半透半反镜BS2和第三反光镜M3调整光轴后，照射至光学探测器B中。
参考光路包括声光调制器F和第二半透半反镜BS3；参考光经声光调制器F驱动，发生偏移，并经第二半透半反镜BS3调整光轴后，照射至光学探测器B中。
更有利的，声光调制器F为AOM布拉格盒。
本发明还公开一种基于激光多普勒测振的磁声耦合成像方法，包括以下步骤：
步骤一、利用控制模块对激励源和激光多普勒振动传感器进行参数设置；
步骤二、利用磁铁产生磁通密度向量为B0的静磁场，使目标体置于静磁场空间内；
步骤三、启动激励源，产生高频激励电流至激励线圈；
步骤四、利用激励线圈产生的环形激励电流，在目标体中引发感应电流密度向量为J的感生涡流；
步骤五、感生涡流在静磁场作用下，产生洛伦兹力F＝J×B0；在洛伦兹力F作用下，目标体内部产生与激励源同频的高频振动；
步骤六、在控制模块的控制下，激光多普勒振动传感器沿轨道移动，对目标体进行360°扫描探测，得到目标体的振动信号；
步骤七、利用控制模块对振动信号进行处理，并根据振动信号，对目标体内部的电学参数分布进行三维重建。
其中，步骤六中激光多普勒振动传感器对目标体进行360°扫描探测，得到目标体的振动信号，具体探测过程如下：
1、利用激光多普勒振动传感器发射频率为f0的激光束，并将激光束分束成两束，其中一束照射至处于振动状态的目标体上，在多普勒效应下，生成频率为fr的反射光；另一束作为参考光。
2、将参考光在驱动频率fB下进行调制，使其发生频移。
在BS1分出一束参考光是为了检出振动的方向性。在参考光路中设置了声光调制器(AOM布拉格盒)，其驱动频率为fB，这样，参考光发生了频移，然后经过半透半反镜BS3调整光轴，照在光学探测器上。
3、将反射光与频移后的参照光在同一空间位置上发生干涉，混频成频率大小f的脉冲光。
4、对脉冲光进行光电转换和频率解调，得到多普勒频率fD大小和与之成比例的目标体的振动速度；计算公式如下：
fD＝f-fB；
v＝λfD/2；
其中，v表示目标体的振动速度，λ表示波长。
fD的推导过程为：
在激光多普勒测振传感器中，多选用连续波He-Ne激光器作为光源，其输出的激光束通过光分束器BS1分成二束，一束输出照射在振动物体上，另一束作为参考光。输出的光通过反射镜M1和M2调整光轴，透过半透半反镜BS2，再穿过透镜，集中照射在振动目标体上。
频率为f0的激光束(波长为λ)照射在振动速度为v的物体表面上，由于产生了多普勒效应，则反射光的频率fr为:
c是光速，由于多普勒效应，照射光频率发生了变化。多普勒频率为:
由于c＞＞v，所以
因为c＝f0λ，所以:
fD＝2v/λ (3)。
多普勒频移反射光通过透镜会聚光线，经BS2、M3调整光轴。这束光在光学探测器上与已频移的参考光重叠，发生干涉。反射光频率为fr＝f0+fD，又因为光在光学探测器上进行光电转换，其拍频为f＝(f0+fD)-(f0-fB)＝fD+fB，所以fD＝f-fB。
5、利用下式计算目标体的振幅A：A＝v/ω；
其中，ω表示目标体以角频率ω作简谐振动。
振幅A推导过程具体为：
假设振动物体以角频率ω作简谐振动，振动物体表面速度为:
v＝Aω·cos(ωt) (4)
其中，A为振幅，又因为当t→0时，cos(ωt)＝1，因此:A＝v/ω。
这样，对于振动是近似简谐振动时，假设已知振动速度和频率，根据方程(5)，就可以计算出高频振动下的位移。
在一个具体实施例中，目标体内洛伦兹力引发的表面振动位移的关系模型为：
其中，ρ为声波引起的密度变化，Cs为目标体内声速；∮V dr′为体积分符号；J(r′)为声源处的电流密度；B0(r′)为声源处的磁场强度；t为时间；r′为声源的位置矢量；rb为探测器在目标体表面的位置矢量。
在一个具体实施例中，目标体内部的电学参数的计算公式为：
其中，σ表示目标体内部组织电导率，J表示感应电流密度向量，B0表示静磁场磁通密度向量，B1表示J引起的交变磁场强度。
目标体内洛伦兹力引发的表面振动位移的关系模型和目标体内部的电学参数的具体推导过程如下：
基于表面振动信息的磁声耦合成像算法。将磁声耦合成像技术应用于生物体时，生物组织的声学特性与流体非常相似，磁声耦合成像的声场方程可以从流体的声学理论出发，选取声学特性均匀的模型进行研究。记洛伦兹力为F＝J×B0，J为感应电流密度向量，B0为静磁场磁通密度向量，B1为J引起的交变磁场强度，cs为组织体内声速，静息时的介质密度为ρ0，静息时介质中的静压强为p0，声波引起的密度变化为ρ，介质中的质点速度为v，位移为u，声压为p。
磁声耦合成像中，介质内的声波幅度很小，声压和密度变化也很小，可以导出线性化耦合方程：
将方程(6)中的第一个方程求梯度，并乘以
有：
方程(6)中的第二个方程对t求导，有：
方程(7)、(8)相减，可消去变量p，有：
进一步，将质点位移定义为对质点速度的积分u＝∫vdt，代入方程(9)，有：
将
代入方程(6)中第一个方程，有：
因此，可通过方程(9)、(10)的求解，得到质点位移(或速度)，进一步通过方程(11)得到声压。同理，激光多普勒测振仪测量得到目标体表面位移(或速度)后，可以通过方程(11)得到目标体表面的声压信息。
在MAT成像中，描述声压信号分布的波动方程为：
使用格林函数方法，获得方程(12)的解析解为：
其中，r和r′分别为探测器和声源的位置矢量，R＝|r-r′|为源与探测器之间的距离。假设探测器在目标体表面rb，即r＝rb时，可将方程(11)代入方程(13)，得到组织内洛伦兹力引发的表面振动位移的关系模型为：
重建得到J(r′)×B0(r′)后，可根据方程(15)计算得到组织内电导率σ参数：
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.一种基于激光多普勒测振的磁声耦合成像***，其特征在于，包括：激励源、激励线圈、磁铁、激光多普勒振动传感器、导轨和控制模块；
其中，所述激励源在所述控制模块的控制下产生激励电流至所述激励线圈；所述激励线圈位于目标体的上方位置，其用于在所述激励电流的作用下，引发所述目标体中产生感应涡流；所述磁铁设置有一对，且分别位于所述目标体的上下方，其用于产生静磁场；所述目标体在所述感应涡流和所述静磁场的共同作用下，产生与所述激励源同频的振动；
所述导轨围绕所述目标体设置；所述激光多普勒振动传感器安装在所述导轨上，且在所述控制模块的控制下沿所述导轨移动并360度采集所述目标体的振动信号；
所述控制模块用于对所述振动信号进行处理，并根据所述振动信号，对所述目标体内部的电学参数分布进行三维重建。
2.根据权利要求1所述的一种基于激光多普勒测振的磁声耦合成像***，其特征在于，所述控制模块包括滤波放大器、数据采集卡、主控单元和计算机；所述滤波放大器、所述数据采集卡和所述计算机依次与所述激光多普勒振动传感器电性连接；所述主控单元分别与所述激励源、所述激光多普勒振动传感器和所述计算机电性连接；
所述滤波放大器用于对所述振动信号进行滤波、除噪和幅值放大处理；
所述数据采集卡用于将处理后的所述振动信号收集并上传至所述计算机；
所述计算机用于根据处理后的所述振动信号对所述目标体内部的电学参数分布进行三维重建，并通过所述主控单元对所述激励源和所述激光多普勒振动传感器进行参数设定和运行管理。
3.根据权利要求1所述的一种基于激光多普勒测振的磁声耦合成像***，其特征在于，所述激光多普勒振动传感器包括激光器、光分束器、反射光路、参考光路、光学探测器、速度解调器和伺服分析模块；
所述激光器用于发射激光束；所述光分束器将所述激光束分成两束，其中一束经所述反射光路照射至所述目标体上，生成反射光并返回至所述光学探测器；另一束作为参考光，并经所述参考光路入射至所述光学探测器；所述光学探测器用于将所述反射光和所述参考光混频成脉冲光，并对脉冲光进行光电转换，生成脉冲电信号；所述速度解调器用于对所述脉冲电信号进行解调，生成速度信号；所述伺服分析模块用于根据积分演算方法，将所述速度信号转换成位移信号。
4.根据权利要求3所述的一种基于激光多普勒测振的磁声耦合成像***，其特征在于，所述反射光路包括第一反光镜、第二反光镜、第一半透半反镜、透镜和第三反光镜；
经所述光分束器分光后的其中一束激光束依次通过所述第一反光镜和所述第二反光镜反射后，通过所述第一半透半反镜和所述透镜后，照射至处于振动状态的所述目标体上，在多普勒效应下产生频率发生偏移的反射光；
生成的反射光经所述透镜汇聚光线后，再依次经所述第一半透半反镜和所述第三反光镜调整光轴后，照射至所述光学探测器中。
5.根据权利要求3所述的一种基于激光多普勒测振的磁声耦合成像***，其特征在于，所述参考光路包括声光调制器和第二半透半反镜；所述参考光经所述声光调制器驱动，发生偏移，并经所述第二半透半反镜调整光轴后，照射至所述光学探测器中。
6.根据权利要求5所述的一种基于激光多普勒测振的磁声耦合成像***，其特征在于，所述声光调制器为AOM布拉格盒。
7.一种基于激光多普勒测振的磁声耦合成像方法，其适用于如权利要求-6任一项所述的基于激光多普勒测振的磁声耦合成像***，其特征在于，包括以下步骤：
利用所述控制模块对所述激励源和所述激光多普勒振动传感器进行参数设置；
利用所述磁铁产生磁通密度向量为B0的静磁场，使目标体置于静磁场空间内；
启动所述激励源，产生高频激励电流至所述激励线圈；
利用所述激励线圈产生的环形激励电流，在所述目标体中引发感应电流密度向量为J的感生涡流；
感生涡流在静磁场作用下，产生洛伦兹力F＝J×B0；在洛伦兹力F作用下，所述目标体内部产生与所述激励源同频的高频振动；
在所述控制模块的控制下，所述激光多普勒振动传感器沿所述轨道移动，对目标体进行360°扫描探测，得到所述目标体的振动信号；
利用所述控制模块对所述振动信号进行处理，并根据所述振动信号，对所述目标体内部的电学参数分布进行三维重建。
8.根据权利要求7所述的一种基于激光多普勒测振的磁声耦合成像方法，其特征在于，所述对目标体进行360°扫描探测，得到所述目标体的振动信号，包括：
利用所述激光多普勒振动传感器发射频率为f0的激光束，并将所述激光束分束成两束，其中一束照射至处于振动状态的所述目标体上，在多普勒效应下，生成频率为fr的反射光；另一束作为参考光；
将参考光在驱动频率fB下进行调制，使其发生频移；
将反射光与频移后的参照光在同一空间位置上发生干涉，混频成频率大小f的脉冲光；
对所述脉冲光进行光电转换和频率解调，得到多普勒频率fD大小和与之成比例的目标体的振动速度；计算公式如下：
fD＝f-fB；
v＝λfD/2；
其中，v表示目标体的振动速度，λ表示波长；
利用下式计算所述目标体的振幅A：
A＝v/ω；
其中，ω表示目标体以角频率ω作简谐振动。
9.根据权利要求7所述的一种基于激光多普勒测振的磁声耦合成像方法，其特征在于，所述目标体内洛伦兹力引发的表面振动位移的关系模型为：
其中，ρ为声波引起的密度变化，Cs为目标体内声速；∮Vdr'为体积分符号；J(r')为声源处的电流密度；B0(r')为声源处的磁场强度；t为时间；r'为声源的位置矢量；rb为探测器在目标体表面的位置矢量。
10.根据权利要求7所述的一种基于激光多普勒测振的磁声耦合成像方法，其特征在于，所述目标体内部的电学参数的计算公式为：
其中，σ表示所述目标体内部组织电导率，J表示感应电流密度向量，B0表示静磁场磁通密度向量，B1表示J引起的交变磁场强度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111137246.XA CN113812926B (zh) | 2021-09-27 | 2021-09-27 | 一种基于激光多普勒测振的磁声耦合成像***及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111137246.XA CN113812926B (zh) | 2021-09-27 | 2021-09-27 | 一种基于激光多普勒测振的磁声耦合成像***及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113812926A true CN113812926A (zh) | 2021-12-21 |
| CN113812926B CN113812926B (zh) | 2024-05-10 |
Family
ID=78921606
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111137246.XA Active CN113812926B (zh) | 2021-09-27 | 2021-09-27 | 一种基于激光多普勒测振的磁声耦合成像***及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113812926B (zh) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070055158A1 (en) * | 2005-08-22 | 2007-03-08 | Siemens Medical Solutions Usa, Inc. | Automated identification of cardiac events with medical ultrasound |
| CN102085096A (zh) * | 2010-12-13 | 2011-06-08 | 中国医学科学院生物医学工程研究所 | 注入电流式磁声耦合成像装置 |
| CN102357036A (zh) * | 2011-11-01 | 2012-02-22 | 东南大学 | 差频超声纵向力磁声电效应检测与成像***的测试台 |
| CN102788836A (zh) * | 2012-07-26 | 2012-11-21 | 中国科学院电工研究所 | 一种磁声显微成像方法及成像*** |
| US20150148675A1 (en) * | 2013-11-23 | 2015-05-28 | Massachusetts Institute Of Technology | System and method for analyzing tissue using shear waves |
| CN104688224A (zh) * | 2015-03-31 | 2015-06-10 | 中国医学科学院生物医学工程研究所 | 一种应用于声学非均匀媒介磁声耦合成像重建方法 |
| CN105816156A (zh) * | 2016-05-17 | 2016-08-03 | 中国医学科学院生物医学工程研究所 | 一种结合热声效应信号的电导率磁声成像装置及成像方法 |
| US20170107558A1 (en) * | 2015-10-14 | 2017-04-20 | The Regents Of The University Of California | Assessment Of Blood Coagulation Using An Acoustic Radiation Force Based Optical Coherence Elastography (ARF-OCE) |
| CN107174202A (zh) * | 2017-05-05 | 2017-09-19 | 深圳大学 | 一种基于主动探测的磁声成像方法及*** |
| KR20180061869A (ko) * | 2016-11-30 | 2018-06-08 | 한국과학기술원 | 텐던 긴장력 모니터링용 센서와 이를 이용한 텐던의 긴장력 진단 시스템 |
| CN109375230A (zh) * | 2018-11-08 | 2019-02-22 | 北京仿真中心 | 一种高精度激光回波频率调制***和方法 |
| CN112244813A (zh) * | 2020-10-22 | 2021-01-22 | 无锡鸣石峻致医疗科技有限公司 | 一种低场核磁共振弹性测量方法及*** |
| CN113080926A (zh) * | 2021-05-11 | 2021-07-09 | 中国科学院电工研究所 | 一种线圈检测式磁声电成像装置及方法 |
-
2021
- 2021-09-27 CN CN202111137246.XA patent/CN113812926B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070055158A1 (en) * | 2005-08-22 | 2007-03-08 | Siemens Medical Solutions Usa, Inc. | Automated identification of cardiac events with medical ultrasound |
| CN102085096A (zh) * | 2010-12-13 | 2011-06-08 | 中国医学科学院生物医学工程研究所 | 注入电流式磁声耦合成像装置 |
| CN102357036A (zh) * | 2011-11-01 | 2012-02-22 | 东南大学 | 差频超声纵向力磁声电效应检测与成像***的测试台 |
| CN102788836A (zh) * | 2012-07-26 | 2012-11-21 | 中国科学院电工研究所 | 一种磁声显微成像方法及成像*** |
| US20150148675A1 (en) * | 2013-11-23 | 2015-05-28 | Massachusetts Institute Of Technology | System and method for analyzing tissue using shear waves |
| CN104688224A (zh) * | 2015-03-31 | 2015-06-10 | 中国医学科学院生物医学工程研究所 | 一种应用于声学非均匀媒介磁声耦合成像重建方法 |
| US20170107558A1 (en) * | 2015-10-14 | 2017-04-20 | The Regents Of The University Of California | Assessment Of Blood Coagulation Using An Acoustic Radiation Force Based Optical Coherence Elastography (ARF-OCE) |
| CN105816156A (zh) * | 2016-05-17 | 2016-08-03 | 中国医学科学院生物医学工程研究所 | 一种结合热声效应信号的电导率磁声成像装置及成像方法 |
| KR20180061869A (ko) * | 2016-11-30 | 2018-06-08 | 한국과학기술원 | 텐던 긴장력 모니터링용 센서와 이를 이용한 텐던의 긴장력 진단 시스템 |
| CN107174202A (zh) * | 2017-05-05 | 2017-09-19 | 深圳大学 | 一种基于主动探测的磁声成像方法及*** |
| CN109375230A (zh) * | 2018-11-08 | 2019-02-22 | 北京仿真中心 | 一种高精度激光回波频率调制***和方法 |
| CN112244813A (zh) * | 2020-10-22 | 2021-01-22 | 无锡鸣石峻致医疗科技有限公司 | 一种低场核磁共振弹性测量方法及*** |
| CN113080926A (zh) * | 2021-05-11 | 2021-07-09 | 中国科学院电工研究所 | 一种线圈检测式磁声电成像装置及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113812926B (zh) | 2024-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Humphrey | Nonlinear propagation in ultrasonic fields: measurements, modelling and harmonic imaging | |
| JP6067590B2 (ja) | 非合焦超音波による超音波振動法 | |
| US6520911B1 (en) | Ultrasound-hall effect imaging system and method | |
| US5921928A (en) | Acoustic force generation by amplitude modulating a sonic beam | |
| US5903516A (en) | Acoustic force generator for detection, imaging and information transmission using the beat signal of multiple intersecting sonic beams | |
| US9488573B2 (en) | Acousto-electromagnetic investigation of physical properties of an object | |
| US9717471B2 (en) | Method and apparatus for multiple-wave doppler velocity meter | |
| CN105816156B (zh) | 一种结合热声效应信号的电导率磁声成像装置及成像方法 | |
| US20070038095A1 (en) | Ultrasound vibrometry | |
| CN107174202B (zh) | 一种基于主动探测的磁声成像方法及*** | |
| JP2002514756A (ja) | 熱音響組織スキャナ | |
| CN113406008B (zh) | 一种基于白光干涉的光声成像装置及方法 | |
| WO2005069997A2 (en) | Enhanced detection of acousto-photonic emissions in optically turbid media using a photo-refractive crystal-based detection system | |
| CN109730722A (zh) | 基于聚焦超声声振信号的弹性成像方法 | |
| US20140031684A1 (en) | System for transcranial ultrasound imaging | |
| CN206080489U (zh) | 基于f‑p腔光纤声传感器的磁声信号检测及成像*** | |
| CN106037638A (zh) | 一种去除热声效应影响的电导率磁声成像装置及成像方法 | |
| JP4874497B2 (ja) | 幅広ビーム映像化 | |
| CN113812926B (zh) | 一种基于激光多普勒测振的磁声耦合成像***及方法 | |
| US7018334B2 (en) | Imaging of a region in a scattering medium | |
| Wu et al. | Fiber optic ultrasonic sensor using Raman-Nath light diffraction | |
| Gunawan et al. | Characterizing acoustic impedance of several saline solution utilizing range finder acoustic sensor | |
| Weng et al. | Orbital-angular-momentum-based image sensor using high resolution photoacoustic tomography | |
| CN110927729A (zh) | 基于位移衰减特性的声辐射力脉冲弹性成像方法 | |
| CN102415900A (zh) | 一种生物医用非接触式电磁超声检测*** |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |