CN108309298A - 一种基于激光超声的磁声电成像装置 - Google Patents
一种基于激光超声的磁声电成像装置 Download PDFInfo
- Publication number
- CN108309298A CN108309298A CN201810035254.5A CN201810035254A CN108309298A CN 108309298 A CN108309298 A CN 108309298A CN 201810035254 A CN201810035254 A CN 201810035254A CN 108309298 A CN108309298 A CN 108309298A
- Authority
- CN
- China
- Prior art keywords
- laser
- ultrasound
- module
- target imaging
- field
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
-
- 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
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
- A61B5/0536—Impedance imaging, e.g. by tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
-
- 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
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5261—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种基于激光超声的磁声电成像装置,包括激光超声激励模块、检测与重建模块和控制与同步模块。控制与同步模块分别与激光超声激励模块及检测与重建模块连接。激光超声激励模块产生超声信号;检测与重建模块获取目标成像体的电参数图像;控制与同步模块控制所述磁声电成像装置协同一致工作。脉冲激光器发射的脉冲激光经过滤光片进行激光的衰减或滤光后,进入激光光束调整***调整,然后进入激光超声激励***,通过光声转化转变为超声束。超声束与静磁场相结合,在目标成像体内产生局部的电场源,通过非接触的检测线圈检测感应的电信号,再经过微弱信号处理子***处理和图像重建,得到目标成像体的重建图像。
Description
技术领域
本发明涉及一种磁声电成像装置。
背景技术
利用组织电学特性的电阻抗成像,具备无损和功能成像的技术优点,是继形态、结构成像之后的新一代医学成像方法,对于生命科学研究和疾病的早期诊断具有重要的价值,可为生命科学研究和临床提供一种全新的诊断信息。一方面,当生物组织发生早期病变而且尚未在形态结构上产生改变时,生物组织内各类化学物质所带电荷量和电荷的空间分布将首先发生变化,宏观上表现为病变部位生物组织的电特性:电阻抗、电导率、介电常数发生改变,因此通过对组织电特性的成像便可以无损检测组织的病变情况,从而达到对病变组织早期诊断早期治疗的目的。美国癌症研究所的专家Layton早于2006年12月撰文,证实了电阻抗成像在乳腺癌和***癌的早期诊断上具有独特的优越性,预言电阻抗成像将以无创、快捷、高准确率、低成本的优势替代活检,成为临床癌症诊断的首选方法。另一方面,电阻抗成像方法可以显示人体各器官完成生理活动,如呼吸、心脏搏动时的变化,不仅可以对人体心血管、食道、胃部等生理活动过程进行连续动态的监测,而且可以发现人体某些组织的病理改变,如癌变、脑水肿等,因此该技术在研究人体生理功能和病变机理方面的具有重要的价值和广泛的应用前景。
磁声电成像方法是2008首次提出用于电导率成像的新方法,随着成像技术的快速发展,该成像方法得到越来越多的关注。目前在方法学的研究基础上,主要集中在如何满足临床需求的低成本早期诊断成像***方面。为提高成像的分辨率和对比度,学者们分别从重建算法、超声激励源、微弱信号检测方法、实验平台等方面进行研究,例如在CN201410773975.8一种磁-声多场耦合成像实验装置中提到的磁声电装置中需要超声屏蔽盒,激励线圈和接收线圈放置在超声屏蔽盒内,因此超声装置与电磁之间有电磁干扰问题,这影响成像的灵敏度和分辨率。
发明内容
本发明的目的是克服现有的磁声电成像方法传统超声激励和电磁检测之间的电磁干扰问题,提高成像的灵敏度和分辨率,提出一种基于激光超声的磁声电成像装置。
所述的基于激光超声的磁声电成像装置包括激光超声激励模块,检测与重建模块和控制与同步模块。控制与同步模块分别与激光超声激励模块及检测与重建模块连接。激光超声激励模块产生强脉冲、高频率、高频宽的超声信号,检测与重建模块获取目标成像体的电参数图像,控制与同步模块的作用是控制整个装置的协同一致工作。
所述的激光超声激励模块包括脉冲激光器、滤波片、小孔光阑、激光光束调整***、磁体、目标成像体和激光超声激励***。脉冲激光器、滤波片、小孔光阑、激光光束调整***和激光超声激励***从前向后依次顺序排列,即脉冲激光器的后面放置滤波片,滤波片的后面为小孔光阑、小孔光阑的后面布置激光光束调整***,激光光束调整***的后面是激光超声激励***。激光超声激励***作用于目标成像体。所述的磁***于目标成像体的周围。磁体作用于目标成像体,在目标成像体内产生一定强度的静磁场,所产生的静磁场大于0.1T。磁体可以是电磁体也可以是永磁体。磁体与检测线圈连接,组成一个具有激励和检测功能的电信号探测探头。所述的检测线圈位于磁体产生的磁场均匀区域,检测线圈的中心轴线与磁场均匀区域的磁场方向平行。脉冲激光器的出光口与滤波片的中心、小孔光阑的中心和激光光束调整***的中心位于同一直线上。脉冲激光器发射的激光经过滤光片滤光、然后经过小孔光阑和激光光束调整***,将激光束调整到需要的尺寸,再经过激光超声激励***后,激光转变为超声信号,到达目标成像体。所的脉冲激光器可以是Nd:YAG脉冲激光器,也可以是其它飞秒脉冲激光器。
激光光束调整***的作用是将激光束调整到所需要的尺寸,例如如果需要产生聚焦超声,则激光光束调整***将激光扩束,以助于后续得到聚焦的超声束;如果需要产生平面超声束,则激光光束调整***可以使其产生发散角比较小的准直光束。激光光束调整***可以是光纤,使其传输更远的距离,有助于在一些非电场合的应用。
所述的激光超声激励***是一个激光超声激励换能器,即为由激光产生高强超声的转换***。激光超声激励***与激光光束调整***配合使用。如果需要产生平面超声,则激光超声激励***的载体为平面载体。平面载体的表面镀有一层平面激光超声增强介质,紧贴在平面激光超声增强介质之上的是一层弹性膜。同理如果需要产生聚焦超声,则激光超声激励***的载体为聚焦载体,相对应的,在聚焦载体上涂有一层聚焦激光超声增强介质,在聚焦激光超声增强介质之上是弹性膜。所述的平面载体和聚焦载体的材质可以是玻璃、也可以是亚克力树脂,还可以是其它具有良好透光度的材质。所述的平面激光超声增强介质和聚焦激光超声增强介质可以是纳米金属介质,也可以是碳纳米颗粒、石墨烯薄膜或者烟灰。所述的弹性膜是具有弹性的保护薄膜,例如聚二甲基硅氧烷。所述的平面载体的形状是平面,聚焦载体的形状是凹面。所述的激光超声激励***可以通过耦合剂直接与目标成像体接触,也可以直接作用到目标成像体的表面;还可以直接把目标成像体的表面作为激光激励超声***,通过脉冲激光器调整后的光束直接激励目标成像体表面,表面产生的超声信号直接耦合到目标成像体的深处,因此也是一种非接触的激光超声激励***。所述的平面激光超声增强介质和聚焦激光超声增强介质可以是一种碳纳米管、石墨烯等一种材料的单结构膜也可以是混合膜结构,比如可以是金属沉积膜和碳纳米颗粒膜的混合膜。
所述的检测与重建模块包括检测线圈,微弱信号处理子***,图像重建算法和计算机显示子***,微弱信号处理子***连接计算机显示子***,图像重建算法和计算机显示子***采集微弱信号处理子***获取的数据,进行图像重建。检测线圈接收到的电信号后接入微弱信号处理子***,处理后的信号通过计算机显示子***进行显示和保存,利用图像重建算法对图像进行重建,具体的重建算法包括时间反转法和压缩感知重建算法。
所述的时间反转法在磁声电成像中借助互易定理变为波动场的逆问题,已知对于任意波动场源H(r),以δ(t)函数激励的场u(r,t),δ(t)为冲激函数,c0为声速,t为时间,若满足如下波动方程:
为拉普拉斯算符。
那么在无穷远边界条件下,该波动方程的解可用格林函数的积分表示为:
其中,r'表示源点,r表示场点,u(r,t)为t时刻在空间r处的激励场,H(r')表示场源分布函数。G(r,t|r',0)表示点源激励下场量u(r,t)的格林函数,它同样满足波动方程。
利用时间反转法求解声场的计算公式为:
其中trd=2T0-t+|r-r'|/c0,S为环绕目标区域的闭合曲面。n为扫描平面的法向分量,u'(r',trd)为激励场势函数的一阶导数,u(r',trd)为激励场势函数,u”(r',trd)为激励场势函数的二阶导数,为扫描平面S上的面积分符号、c0为声速、T0为时间场的反转时刻,T0需要足够大,使得对于t>T0的时刻u(r,t)=0。
因为波动场源与初始时刻场量的一阶导数有关,也就是所以可对声场的计算公式两端求导,并取0时刻的值,得到场源的表达式为:
由此即可重建出分布波动源H(r'),其中π为3.14,n为扫描平面的法向分量,u'(r',trd)为激励场量的一阶导数,u”(r',trd)为激励场量的二阶导数,为扫描平面S上的面积分符号。
所述的压缩感知重建算法原理如下:
首先借助互易定理,检测线圈检测到的感应电动势uAB为:
其中uAB为检测线圈检测到的感应电动势,表示在激光激励超声***和静磁场共同作用的区域Ω内的积分,为互易电流密度的旋度,为速度势函数,静磁场B0、目标成像体密度ρ0、作用区域Ω内,r'为超声换能器激励位置,r为超声换能器所在点的位置,J2为互易得电流密度、B0为静磁场磁感应强度、ρ0为目标成像体的密度、为速度势、为哈密顿(Hamilton)算子。
该公式说明,磁声电信号uAB的大小是由互易电流密度的旋度速度势函数静磁场磁感应强度B0和密度ρ0在作用区域Ω内的相互作用决定的。感应电动势uAB除了是时间的函数以外,它还是与超声换能器激励位置r'有关的函数。用r表示超声换能器所在位置,则uAB可改写为u(r,t)。
采用压缩感知的方法重建磁声电信号分布场源的方法如下:
令压缩感知矩阵KΦ=ACS,则由正交匹配追踪算法求解的矩阵方程为:
U=ACSθ (2)
其中,U为变换矩阵,Φ是小波基下的正交变换矩阵,例如傅立叶基或小波基,K是由矩量法求得的系数矩阵,同时满足高斯分布。因此可利用ACS=KΦ直接得到压缩感知矩阵ACS。θ为转换系数,正交匹配追踪算法的计算方法如下所述:
第一步:首先初始化残差r0=U,搜索列向量集合Λ=空集;
第二步:在第i次迭代中,计算压缩感知矩阵的每列与残差之间的相关系数<ri,aj>,将相关系数最大的列aj并入Λ,即Λ=[Λ,αj],其中j表示感知矩阵ACS的第j列向量,ri表示第i次迭代的残差;
第三步:利用最小二乘法逼近测量信号:其中为测量信号,||||为绝对值表达式,argmin表示使函数U-ΛH取得最小值的所有自变量H的集合;
第四步:利用公式更新残差;
第五步:如果没有达到迭代步数则返回第二步继续迭代,否则输出
其中,为磁声电信号重建分布场源。
为了得到电导率分布,借助赫姆霍兹定理,由互易电流密度旋度和边界条件重建互易电流密度矢量,然后再由电流密度矢量重建电导率。
所述的控制与同步模块包括三维控制台支架和控制子***,控制子***与三维控制台支架上的三个步进电机连接,首先控制子***发送一个脉冲信号作为外触发信号,外触发信号分别控制激光超声激励模块脉冲激光器、检测与重建模块的微弱电信号处理子***、图像重建算法和计算机显示子***,以及三维控制台支架。在控制子***的控制下,脉冲激光器发射激光信号,同时微弱电信号处理子***接收检测线圈的信号并处理。微弱电信号处理子***每采集一次信号后,在控制子***的控制下,三维控制台支架在步进电机的带动下移动一步,进行下一个位置的采集,位于三维控制台支架上的目标成像体由此实现其三维移动。
所述的激光超声非接触式磁声电成像装置可以为便携式装置。若为便携式装置,磁体和检测线圈可以为一体化的探测探头。在控制子***的控制下,通过控制磁体和检测线圈的移动实现三维扫描,也可以固定磁体和检测线圈,通过三维控制台支架控制目标成像体移动。
本发明工作过程如下:激光超声激励模块的脉冲激光器发射的脉冲激光首先经过滤光片进行激光的衰减或滤光,之后经过小孔光阑,进入激光光束的调整***,使激光源调整为所需尺寸的激光束。然后,激光束进入激光超声激励***,通过光声转化,使其成为转变为超声束,超声束与静磁场相结合,在目标成像体内产生局部的电场源,通过非接触的检测线圈检测感应的电信号,检测电信号经过微弱信号处理子***和重建算法,最后在计算机显示子***显示重建的图像。
附图说明
图1本发明基于激光超声的非接触式磁声电成像方法与装置示意图;
图2激光超声转换传感器示意图;
图中:A01脉冲激光器,A02激光器出光口,A03滤波片,A04小孔光阑,A05激光光束调整***,A06磁体,A07检测线圈,A08目标成像体,A09微弱电信号处理子***,A10计算机显示子***,A11激光超声激励***,A12三维控制台支架A13控制子***,B01平面载体,B02平面激光超声增强介质B03弹性膜,C01聚焦载体,C02聚焦激光超声增强介质,C03弹性膜。
具体实施方式
以下结合附图和具体实施方式进一步说明本发明。
如图1所示,本发明基于激光超声的磁声电成像装置包括激光超声激励模块、检测与重建模块和控制与同步模块。控制与同步模块分别与激光超声激励模块及检测与重建模块连接。
激光超声激励模块的作用是产生强脉冲、高频率、高频宽的超声信号;检测与重建模块的作用是获取目标成像体A08的电参数图像;控制与同步模块的作用是控制整个装置协同一致工作。
所述的激光超声激励模块包括脉冲激光器A01、目标成像体A08、滤波片A03、小孔光阑A04、激光光束调整***A05、磁体A06,以及激光超声激励***A11。脉冲激光器、滤波片、小孔光阑、激光光束调整***和激光超声激励***从前向后依次顺序排列,即脉冲激光器A01后面是滤波片A03,滤波片A03后面是小孔光阑A04、小孔光阑A04后面是激光光束调整***A05,激光光束调整A05***后面是激光超声激励***A11。激光超声激励***A11作用于目标成像体A08。所述的脉冲激光器A01的出光口与滤波片A03的中心、小孔光阑A04的中心、激光光束调整***A05的中心位于同一直线上。脉冲激光器A01的激光出射口A02发射的激光经过滤波片A03滤光,然后经过小孔光阑A04调整到需要的尺寸,调整光束尺寸后的激光束射入激光光束调整***A05,然后再经过激光超声激励***A11后激光变为超声信号。所述的磁体A06位于目标成像体A08的周围,磁体A06作用于目标成像体A08,在目标成像体A08内产生一定强度的静磁场,所产生的静磁场大于0.1T。磁体A06可以是电磁体也可以是永磁体,磁体A06与检测线圈A07连接,组成一个具有激励和检测功能的电信号探测探头。所述的检测线圈A07位于磁体产生的磁场均匀区域,检测线圈A07的中心轴线与均匀区域的磁场方向平行。
所的脉冲激光器可以是Nd:YAG脉冲激光器也可以是其它飞秒脉冲激光器。
所述的激光光束调整***A05对激光束进行调整。如果需要产生聚焦超声,则激光光束调整***A05为激光扩束***,使激光束达到一定尺寸,以助于后续得到聚焦的超声束。如果需要产生平面超声束,则激光光束调整***A05使激光束产生发散角比较小的准直光束。激光光束调整***可以是光纤,使激光束的传输距离更远,有助于在一些非电场合的应用。
所述的激光超声激励***A11是一个激光超声激励换能器,即由激光产生高强超声的转换***。激光超声激励***A11与激光光束调整***A05配合使用,如果需要产生平面超声,则激光超声激励***的载体为平面载体B01,平面载体B01的表面镀有一层平面激光超声增强介质B02,平面激光超声增强介质B02之上是一层弹性膜B03。如果需要产生聚焦超声,则激光超声激励***的载体为聚焦载体C01,相对应的,在聚焦载体C01的上表面涂覆有一层聚焦激光超声增强介质C02,聚焦激光超声增强介质C02之上是弹性膜C03。所述的平面载体B01和聚焦载体C01的材质是玻璃、亚克力树脂,或其它具有良好透光度的材质。所述的平面激光超声增强介质B02和聚焦激光超声增强介质C02可以是纳米金属介质,也可以是碳纳米颗粒、石墨烯薄膜或者烟灰。所述的弹性膜B03和聚焦激光超声增强介质C03可以是具有弹性的保护薄膜,例如聚二甲基硅氧烷。平面载体B01的形状是平面,聚焦载体C01的形状是凹面。所述的激光超声激励***A11可以通过耦合剂直接与目标成像体接触,也可以直接作用到目标成像体A08的表面,还可以直接把目标成像体A08的表面作为激光激励超声***,通过脉冲激光器A01调整后的光束直接激励目标成像体A08表面,目标成像体A08表面产生的超声信号直接耦合到目标成像体A08深处,因此也是一种非接触的激光超声激励***。
所述的检测与重建模块包括检测线圈A07,微弱信号处理子***A09,图像重建算法和计算机显示子***A10,微弱信号处理子***A09连接计算机显示子***A10,图像重建算法和计算机显示子***A10采集微弱信号处理子***A09获取的数据,进行图像重建。检测线圈A07接收到的电信号送入微弱信号处理子***A09,微弱信号处理子***A09处理后,信号由计算机显示子***A10进行显示和保存,利用图像重建算法对目标成像体A08进行图像重建。
所述的控制与同步模块包括三维控制台支架A12和控制子***A13。控制子***A13分别与三维控制台支架A12上的三个步进电机连接,控制三个步进电机的运动。首先控制子***A13发送一个脉冲信号作为外触发信号,分别作用于脉冲激光器A01、微弱电信号处理子***A09、计算机显示子***A10和三维控制台支架A12上的三个步进电机。在控制子***A13的控制下,脉冲激光器A01发射激光信号,同时微弱电信号处理子***A09开始接收检测线圈A07的信号并处理,每采集一次信号后,三维控制台支架A12在步进电机的带动下移动一步,进行下一个位置的采集,由此实现目标成像体A08的三维移动。
所述的激光超声非接触式磁声电成像装置可以为便携式装置,若为便携式装置,则可以把磁体A06和检测线圈A07做成一体化的探测探头,在控制子***的控制下,可以通过控制磁体A06和检测线圈A07的移动实现三维的扫描,也可以固定磁体A06和检测线圈A07,通过三维控制台支架A12控制目标成像体A08的移动。
所述的平面激光超声增强介质B02和聚焦激光超声增强介质C02可以是单层膜结构也可以是混合膜结构,比如可以是金属沉积膜和碳纳米颗粒膜的混合膜。
本发明基于激光超声的磁声电成像装置的工作过程如下:
首先脉冲激光器A01发射一束脉冲激光,脉冲激光经过激光器出光口A02射出后,首先经过滤光片A03进行激光的衰减或滤光,之后进入激光光束调整***A05,使其成为可使用的激光束。然后激光束进入激光超声激励***A11,通过光声转化,使激光束转变为可利用的超声束,所述的超声束与静磁场相结合,在目标成像体A08内产生局部的电场源。通过非接触的检测线圈A07检测磁声电信号,检测到的电信号经过微弱信号处理子***A09,通过图像重建算法重建图像,在计算机显示子***A10显示。
Claims (8)
1.一种基于激光超声的磁声电成像装置,其特征在于:所述的装置包括激光超声激励模块、检测与重建模块和控制与同步模块;控制与同步模块分别与激光超声激励模块及检测与重建模块连接;激光超声激励模块产生超声信号;检测与重建模块获取目标成像体(A08)的电参数图像;控制与同步模块控制所述磁声电成像装置协同一致工作。
2.如权利要求1所述的基于激光超声的磁声电成像装置,其特征在于:所述的激光超声激励模块包括脉冲激光器(A01)、目标成像体(A08)、滤波片(A03)、小孔光阑(A04)、激光光束调整***(A05)、磁体(A06),以及激光超声激励***(A11);脉冲激光器(A01)、滤波片(A03)、小孔光阑(A04)、激光光束调整***(A05)和激光超声激励***(A11)从前向后依次顺序排列,即脉冲激光器(A01)后面是滤波片(A03),滤波片(A03)后面是小孔光阑(A04)、小孔光阑(A04)后面是激光光束调整***(A05),激光光束调整***(A05)后面是激光超声激励***(A11);激光超声激励***(A11)作用于目标成像体(A08);所述的脉冲激光器(A01)的出光口与滤波片(A03)的中心、小孔光阑(A04)的中心、激光光束调整***(A05)的中心位于同一直线上;脉冲激光器A01的激光出射口(A02)发射的激光经过滤波片(A03)滤光,经过小孔光阑(A04)调整到需要的尺寸,然后射入激光光束调整***(A05),再经过激光超声激励***(A11)后转变为超声信号;所述的磁体(A06)位于目标成像体(A08)周围,磁体(A06)作用于目标成像体(A08),在目标成像体(A08)内产生静磁场;磁体(A06)与检测线圈(A07)连接,组成探测探头;所述的检测线圈(A07)位于磁体产生的磁场均匀区域,检测线圈(A07)的中心轴线与均匀区域的磁场方向平行。
3.如权利要求1所述的基于激光超声的磁声电成像装置,其特征在于:所述的激光超声激励***(A11)为激光超声激励换能器;需要产生平面超声时,激光超声激励***(A11)的载体为平面载体(B01),平面载体(B01)的表面镀有一层平面激光超声增强介质(B02),平面激光超声增强介质(B02)之上有一层弹性膜(B03);需要产生聚焦超声时,激光超声激励***的载体为聚焦载体(C01),在聚焦载体(C01)的上表面涂覆有一层聚焦激光超声增强介质(C02),聚焦激光超声增强介质(C02)之上是弹性膜(C03);平面载体(B01)为平面,聚焦载体(C01)为凹面;所述的激光超声激励***(A11)通过耦合剂直接与目标成像体接触,或者直接作用到目标成像体(A08)的表面。
4.如权利要求3所述的基于激光超声的磁声电成像装置,其特征在于:所述的目标成像体(A08)的表面作为激光激励超声***;通过脉冲激光器(A01)调整后的激光束直接作用于目标成像体(A08)表面,目标成像体(A08)的表面产生的超声信号直接耦合到目标成像体(A08)的深处。
5.如权利要求1所述的基于激光超声的磁声电成像装置,其特征在于:所述的控制与同步模块包括三维控制台支架(A12)和控制子***(A13);控制子***(A13)与三维控制台支架(A12)上的三个步进电机连接;首先控制子***(A13)发送脉冲信号作为外触发信号,分别控制脉冲激光器(A01)、微弱电信号处理子***(A09)、计算机显示子***(A10)和三维控制台支架(A12);在控制子***(A13)的控制下,脉冲激光器(A01)-发射激光信号,微弱电信号处理子***(A09)接收检测线圈(A07)的检测信号并处理,每采集一次信号后三维控制台支架(A12)在步进电机的带动下移动一步,进行下一个位置的采集,由此实现目标成像体(A08)的三维移动。
6.如权利要求1所述的基于激光超声的磁声电成像装置,其特征在于:所述的检测与重建模块包括检测线圈(A07),微弱信号处理子***(A09),图像重建算法和计算机显示子***(A10),微弱信号处理子***(A09)连接计算机显示子***(A10),图像重建算法和计算机显示子***(A10)采集微弱信号处理子***(A09)获取的数据,进行图像重建;检测线圈(A07)接收到的电信号送入微弱信号处理子***(A09)处理,之后由计算机显示子***(A10)保存,利用图像重建算法对目标成像体(A08)进行图像重建,所述的图像重建算法包括时间反转法和压缩感知重建算法。
7.如权利要求6所述的基于激光超声的磁声电成像装置,其特征在于:所述的时间反转法在磁声电成像中借助互易定理变为波动场的逆问题,
已知对于任意波动场源H(r),以δ(t)函数激励的场u(r,t),δ(t)为冲激函数,c0为声速,t为时间,若满足如下波动方程:
▽为拉普拉斯算符;
那么在无穷远边界条件下,该波动方程的解可用格林函数的积分表示为:
其中,r'表示源点,r表示场点,u(r,t)为t时刻在空间r处的激励场,H(r')表示场源分布函数。G(r,t|r',0)表示点源激励下场量u(r,t)的格林函数,它同样满足波动方程;
利用时间反转法求解声场的计算公式为:
其中trd=2T0-t+|r-r'|/c0,S为环绕目标区域的闭合曲面。n为扫描平面的法向分量,u'(r',trd)为激励场势函数的一阶导数,u(r',trd)为激励场势函数,u”(r',trd)为激励场势函数的二阶导数,为扫描平面S上的面积分符号、c0为声速、T0为时间场的反转时刻,T0需要足够大,使得对于t>T0的时刻u(r,t)=0;
因为波动场源与初始时刻场量的一阶导数有关,也就是所以对声场的计算公式两端求导,并取0时刻的值,得到场源的表达式为:
由此即能够重建出分布波动源H(r'),其中π为3.14,n为扫描平面的法向分量,u'(r',trd)为激励场量的一阶导数,u”(r',trd)为激励场量的二阶导数,为扫描平面S上的面积分符号。
8.如权利要求6所述的基于激光超声的磁声电成像装置,其特征在于:所述的压缩感知重建算法重建磁声电信号分布场源的方法如下:
令压缩感知矩阵KΦ=ACS,则由正交匹配追踪算法求解的矩阵方程为:
U=ACSθ (2)
其中,U为变换矩阵,Φ是小波基下的正交变换矩阵,在傅立叶基或小波基条件下,K是由矩量法求得的系数矩阵,同时满足高斯分布;因此利用ACS=KΦ直接得到压缩感知矩阵ACS;θ为转换系数,正交匹配追踪算法的计算方法如下所述:
第一步:首先初始化残差r0=U,搜索列向量集合Λ=空集;
第二步:在第i次迭代中,计算压缩感知矩阵的每列与残差之间的相关系数<ri,aj>,将相关系数最大的列aj并入Λ,即Λ=[Λ,αj],其中j表示感知矩阵ACS的第j列向量,ri表示第i次迭代的残差;
第三步:利用最小二乘法逼近测量信号:其中为测量信号,||||为绝对值表达式,argmin表示使函数U-ΛH取得最小值的所有自变量H的集合;
第四步:利用公式更新残差;
第五步:如果没有达到迭代步数则返回第二步继续迭代,否则输出其中,为磁声电信号重建分布场源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810035254.5A CN108309298B (zh) | 2018-01-15 | 2018-01-15 | 一种基于激光超声的磁声电成像装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810035254.5A CN108309298B (zh) | 2018-01-15 | 2018-01-15 | 一种基于激光超声的磁声电成像装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108309298A true CN108309298A (zh) | 2018-07-24 |
CN108309298B CN108309298B (zh) | 2021-04-09 |
Family
ID=62893798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810035254.5A Active CN108309298B (zh) | 2018-01-15 | 2018-01-15 | 一种基于激光超声的磁声电成像装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108309298B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107495965A (zh) * | 2017-08-25 | 2017-12-22 | 中国科学院电工研究所 | 一种激光聚焦超声激励的磁声电成像方法与装置 |
CN109141493A (zh) * | 2018-09-25 | 2019-01-04 | 中国科学院电工研究所 | 光驱动的超声探头及其超声成像*** |
CN110710972A (zh) * | 2019-10-21 | 2020-01-21 | 中国石油大学(华东) | 一种线圈检测式磁声电层析成像的信息处理方法 |
CN111239730A (zh) * | 2020-01-19 | 2020-06-05 | 浙江大学 | 一种基于时间反转和压缩感知的电磁非视线成像方法 |
CN111481198A (zh) * | 2020-04-24 | 2020-08-04 | 中国科学院电工研究所 | 磁声监控下的激光诱导超声激励方法与装置 |
CN113080926A (zh) * | 2021-05-11 | 2021-07-09 | 中国科学院电工研究所 | 一种线圈检测式磁声电成像装置及方法 |
CN115607112A (zh) * | 2022-11-29 | 2023-01-17 | 暨南大学附属第一医院(广州华侨医院) | 一种基于光磁声的一体化智能成像***及方法 |
CN115868956A (zh) * | 2023-03-01 | 2023-03-31 | 暨南大学附属第一医院(广州华侨医院) | 一种用于光磁声成像的激励源抗干扰方法 |
CN117547242A (zh) * | 2024-01-12 | 2024-02-13 | 杭州永川科技有限公司 | 磁感应断层成像设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101247758A (zh) * | 2005-05-11 | 2008-08-20 | 明尼苏达大学评议会 | 用于利用磁感应进行成像的方法和设备 |
CN102805621A (zh) * | 2012-07-26 | 2012-12-05 | 中国科学院电工研究所 | 一种磁声电成像***及成像方法 |
CN102860825A (zh) * | 2012-10-16 | 2013-01-09 | 中国科学院电工研究所 | 基于洛仑兹力力学效应的磁声电阻抗成像***及成像方法 |
CN102894974A (zh) * | 2012-10-16 | 2013-01-30 | 中国科学院电工研究所 | 一种磁声电成像***及成像方法 |
US20150366458A1 (en) * | 2012-05-07 | 2015-12-24 | Helmholtz Zentrum Munchen Deutsches Forschungszentrum Fur Gesundheit Und Umwelt (Gmbh) | Apparatus and method for frequency-domain thermo-acoustic tomographic imaging |
CN107495965A (zh) * | 2017-08-25 | 2017-12-22 | 中国科学院电工研究所 | 一种激光聚焦超声激励的磁声电成像方法与装置 |
-
2018
- 2018-01-15 CN CN201810035254.5A patent/CN108309298B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101247758A (zh) * | 2005-05-11 | 2008-08-20 | 明尼苏达大学评议会 | 用于利用磁感应进行成像的方法和设备 |
US20150366458A1 (en) * | 2012-05-07 | 2015-12-24 | Helmholtz Zentrum Munchen Deutsches Forschungszentrum Fur Gesundheit Und Umwelt (Gmbh) | Apparatus and method for frequency-domain thermo-acoustic tomographic imaging |
CN102805621A (zh) * | 2012-07-26 | 2012-12-05 | 中国科学院电工研究所 | 一种磁声电成像***及成像方法 |
CN102860825A (zh) * | 2012-10-16 | 2013-01-09 | 中国科学院电工研究所 | 基于洛仑兹力力学效应的磁声电阻抗成像***及成像方法 |
CN102894974A (zh) * | 2012-10-16 | 2013-01-30 | 中国科学院电工研究所 | 一种磁声电成像***及成像方法 |
CN107495965A (zh) * | 2017-08-25 | 2017-12-22 | 中国科学院电工研究所 | 一种激光聚焦超声激励的磁声电成像方法与装置 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107495965A (zh) * | 2017-08-25 | 2017-12-22 | 中国科学院电工研究所 | 一种激光聚焦超声激励的磁声电成像方法与装置 |
CN109141493A (zh) * | 2018-09-25 | 2019-01-04 | 中国科学院电工研究所 | 光驱动的超声探头及其超声成像*** |
CN109141493B (zh) * | 2018-09-25 | 2020-08-25 | 中国科学院电工研究所 | 光驱动的超声探头及其超声成像*** |
CN110710972A (zh) * | 2019-10-21 | 2020-01-21 | 中国石油大学(华东) | 一种线圈检测式磁声电层析成像的信息处理方法 |
CN111239730A (zh) * | 2020-01-19 | 2020-06-05 | 浙江大学 | 一种基于时间反转和压缩感知的电磁非视线成像方法 |
CN111481198A (zh) * | 2020-04-24 | 2020-08-04 | 中国科学院电工研究所 | 磁声监控下的激光诱导超声激励方法与装置 |
CN111481198B (zh) * | 2020-04-24 | 2023-05-02 | 中国科学院电工研究所 | 磁声监控下的激光诱导超声激励方法与装置 |
CN113080926A (zh) * | 2021-05-11 | 2021-07-09 | 中国科学院电工研究所 | 一种线圈检测式磁声电成像装置及方法 |
CN115607112A (zh) * | 2022-11-29 | 2023-01-17 | 暨南大学附属第一医院(广州华侨医院) | 一种基于光磁声的一体化智能成像***及方法 |
CN115868956A (zh) * | 2023-03-01 | 2023-03-31 | 暨南大学附属第一医院(广州华侨医院) | 一种用于光磁声成像的激励源抗干扰方法 |
CN117547242A (zh) * | 2024-01-12 | 2024-02-13 | 杭州永川科技有限公司 | 磁感应断层成像设备 |
CN117547242B (zh) * | 2024-01-12 | 2024-05-14 | 杭州永川科技有限公司 | 磁感应断层成像设备 |
Also Published As
Publication number | Publication date |
---|---|
CN108309298B (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108309298A (zh) | 一种基于激光超声的磁声电成像装置 | |
AU2017203545B2 (en) | Transmissive imaging and related apparatus and methods | |
US6535625B1 (en) | Magneto-acoustic imaging | |
US6974415B2 (en) | Electromagnetic-acoustic imaging | |
CN103356189B (zh) | 磁共振和超声参数化图像融合的方法 | |
US6490471B2 (en) | Electromagnetical imaging and therapeutic (EMIT) systems | |
CN103415258B (zh) | 用于***的超声检测的***和方法 | |
CN107669270B (zh) | 一种电磁聚焦超声激励的磁声电成像方法与装置 | |
CN108294751B (zh) | 一种磁声电-超声检测装置 | |
US9508154B2 (en) | Medical imaging apparatus and method of operating the same | |
CN103025248A (zh) | 图像信息获取装置、图像信息获取方法及图像信息获取程序 | |
JP2011183149A (ja) | 測定装置 | |
Clement et al. | Superresolution ultrasound imaging using back-projected reconstruction | |
US20160089117A1 (en) | Ultrasound imaging apparatus and method using synthetic aperture focusing | |
KR20150041723A (ko) | 초음파 프로브 및 이를 포함한 의료 장치 | |
US11272906B2 (en) | Ultrasonic imaging device and method for controlling same | |
KR102218308B1 (ko) | 초음파 영상 처리 장치 및 방법 | |
US10018598B2 (en) | Ultrasonic imaging apparatus and method of controlling the same | |
JP2017140092A (ja) | 被検体情報取得装置 | |
KR20160054360A (ko) | 영상 장치 및 영상화 방법 | |
CN107536608A (zh) | 基于灵敏度场理论的非侵入式声电成像方法 | |
KR101611449B1 (ko) | 프로브를 구비한 의료영상 장치의 제어 방법 및 그 의료영상 장치 | |
RU2662079C1 (ru) | Способ микроволновой томографии сверхвысокого разрешения | |
Sun et al. | Radiation theory comparison for magnetoacoustic tomography with magnetic induction (MAT-MI) | |
JP6012776B2 (ja) | 被検体情報取得装置および被検体情報取得装置の制御方法 |
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 |