CN117064336A - 一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法 - Google Patents
一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法 Download PDFInfo
- Publication number
- CN117064336A CN117064336A CN202311041135.8A CN202311041135A CN117064336A CN 117064336 A CN117064336 A CN 117064336A CN 202311041135 A CN202311041135 A CN 202311041135A CN 117064336 A CN117064336 A CN 117064336A
- Authority
- CN
- China
- Prior art keywords
- raman
- contrast agent
- scanning
- transmission
- positioning
- 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.)
- Pending
Links
- 238000001069 Raman spectroscopy Methods 0.000 title claims abstract description 139
- 239000002872 contrast media Substances 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000004929 transmission Raman spectroscopy Methods 0.000 title claims abstract description 34
- 238000001228 spectrum Methods 0.000 title claims abstract description 21
- 238000001514 detection method Methods 0.000 claims abstract description 32
- 230000003595 spectral effect Effects 0.000 claims abstract description 28
- 239000000523 sample Substances 0.000 claims description 102
- 230000005540 biological transmission Effects 0.000 claims description 26
- 238000001237 Raman spectrum Methods 0.000 claims description 25
- 239000002105 nanoparticle Substances 0.000 claims description 8
- 238000012986 modification Methods 0.000 claims description 7
- 230000004048 modification Effects 0.000 claims description 7
- 239000002202 Polyethylene glycol Substances 0.000 claims description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 3
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 claims description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 3
- 239000010931 gold Substances 0.000 claims description 3
- 229910052737 gold Inorganic materials 0.000 claims description 3
- 229920001223 polyethylene glycol Polymers 0.000 claims description 3
- 229910052814 silicon oxide Inorganic materials 0.000 claims description 3
- 229910052709 silver Inorganic materials 0.000 claims description 3
- 239000004332 silver Substances 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 abstract description 6
- 206010028980 Neoplasm Diseases 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 2
- 210000001519 tissue Anatomy 0.000 description 37
- 239000002245 particle Substances 0.000 description 13
- 230000003902 lesion Effects 0.000 description 12
- 235000015277 pork Nutrition 0.000 description 7
- 239000011258 core-shell material Substances 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000004807 localization Effects 0.000 description 5
- 238000001356 surgical procedure Methods 0.000 description 5
- 229920000936 Agarose Polymers 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 238000013519 translation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004416 surface enhanced Raman spectroscopy Methods 0.000 description 2
- 238000000479 surface-enhanced Raman spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010052428 Wound Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000010171 animal model Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000002925 chemical effect Effects 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000002923 metal particle Substances 0.000 description 1
- 230000001394 metastastic effect Effects 0.000 description 1
- 206010061289 metastatic neoplasm Diseases 0.000 description 1
- 239000002086 nanomaterial Substances 0.000 description 1
- 229910000510 noble metal Inorganic materials 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000001766 physiological effect Effects 0.000 description 1
- 238000010882 preoperative diagnosis Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Classifications
-
- 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
- A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
-
- 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
- A61B5/0062—Arrangements for scanning
-
- 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/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- Veterinary Medicine (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明公开了一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,涉及医疗技术领域。可以完成厚组织中拉曼造影剂的三维实时定位;开发“两步法”扫描策略,首先沿圆柱组织轴向扫描一次,确定造影剂所在平面,然后对于造影剂所在平面进行环形扫描,即可完成病灶的定位;采用基于“谱峰比对数值”包含深度信息的信号,结合最小误差原则,进行病灶位置乃至图像的重构。本发明提供了一种可以进行深层病灶定位的拉曼光学技术,为临床的肿瘤检测以及定位提供更多选择,并为精准医疗提供新的检测定位方法。
Description
技术领域
本发明涉及医疗技术领域,尤其涉及一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法。
背景技术
无创检测体内深部病灶并获得其深度信息,对于临床医学诊断和治疗至关重要。人体的病灶如实体瘤等分布深浅不一,大多数病灶的解剖位置往往较深,从数厘米到十几厘米不等。外科手术是目前治疗恶性实体瘤和清扫转移***的主要方法。发展高精准的深层检测方法对于外科手术具有重要意义,不仅可以实现术前经皮定位、帮助规划手术减小创口,又能在术中实时引导、发现被掩盖的肉眼不可见的微小病灶,从而彻底切除病灶又尽可能保护正常组织,辅助术周导航,实现精准诊疗。
当前,X-CT、PET和MRI等多种技术在临床均有了广泛的应用。这些技术能够实现无创深层诊断,但受限于采集时间长、仪器昂贵或辐射剂量大等原因,更常用于初始诊断或术前观察,如在术前三维重建得到解剖或生理信息,为医生手术规划提供参考。但是,病灶的解剖位置可能会随着人体姿势改变或自然生理活动(如胸廓呼吸运动)发生移动,导致术中病灶位置与术前检测结果相比存在偏移,这为基于术前建模的虚拟或电磁术中引导等技术带来困难。因此,当前临床常用的检测手段在术前诊断和术中应用之间的衔接仍存在不足,难以同时实现无创诊断深层病灶与实时引导病灶的手术治疗。
表面增强拉曼光谱(Surface-Enhanced Raman Spectroscopy,SERS)是指由于表面等离激元的激发以及表面化学效应使得吸附在金属纳米结构表面的分子的拉曼散射产生极大增强的效应,可实现单分子和单颗粒检测,具有非常高的灵敏性。通过将拉曼分子吸附在粗糙贵金属颗粒表面,可以制备得到特异性好、灵敏度高的表面增强拉曼光谱探针。基于表面增强拉曼光谱探针(简称拉曼探针)的造影剂不易发生光漂白、不受生物组织自发荧光的干扰、且样品制备简单。此外,其特有的指纹式图谱使其具有超高的特异性,超窄的半宽峰有利于光学编码的多指标检测,且光学稳定性强。基于以上突出的优势,表面增强拉曼探针在生物医学领域具有非常广阔的应用前景。
但是,虽然现有的拉曼探针造影剂可以主动或被动地靶向病灶,但是如何确定造影剂在厚组织中的位置是一个挑战。当前研究存在:难以在临床安全激光照射条件下实现高穿透深度、并未解决病灶信号深度溯源,即“测不深”和“测不准”两大挑战。
拉曼光谱技术面临的瓶颈之一无法在安全激光照射下实现深层病灶的检测。增加激光功率是提升检测深度的直接方法,但也会导致单位面积中能量释放的速度提升,引起治疗副作用,导致生物组织表面迅速升温引起灼伤。目前几乎所有应用于动物模型的拉曼光谱检测都采用了显著高于MPE的激光光照,特别地,为了提升信噪比,往往采取大幅提高激光功率或延长采集时间的手段。
拉曼光谱技术面临的瓶颈之二是如何确定病灶深度。由于深层检测过程中获得的信号包含整个截面的光子信息,不具备纵深分辨率,所测得光谱是组织与颗粒信号的混合光谱。如何从混合光谱中提取颗粒的有效信息并判断信号深度,是亟需解决的问题。以往的透射拉曼光谱研究工作大多数只实现检测颗粒信号就戛然而止,但如果无法通过信号判断颗粒的大致位置,实现检测也仅仅意味着知晓这团厚组织中存在颗粒或病灶,这是远远不够,也无法解决实际临床问题的。
因此,本领域的技术人员致力于开发一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法。提供一种可以进行深层病灶定位的拉曼光学技术,为临床的肿瘤检测以及定位提供更多选择,并为精准医疗提供新的检测定位方法。
发明内容
有鉴于现有技术的上述缺陷,本发明所要解决的技术问题是在安全激光照射下实现深层病灶的检测,确定病灶深度,为临床的肿瘤检测以及定位提供更多选择,并为精准医疗提供新的检测定位方法。
为实现上述目的,本发明提供了一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,包括以下步骤:
步骤1、将含有拉曼探针造影剂的柱状生物组织样品置于样品台上,采用拉曼光谱仪在横向进行透射检测,即光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号;
步骤2、沿圆柱体的轴向扫描,即同时调整光源和拉曼采集探头的高度,并且保持处于同一高度以构成透射检测装置,根据扫描结果拟合得到曲线,并取信号强度最大位置为造影剂所在高度位置;
步骤3、对于所述步骤2选取的造影剂所在高度的平面进行环形扫描,即保持光源和拉曼采集探头位于对侧构成透射装置,通过同时旋转和拉曼采集探头或者仅旋转样品的方法,检测得到360°各个角度下的拉曼光谱;
步骤4、计算所述步骤3所得的拉曼光谱的谱峰比的对数值,并且根据谱峰比对数值-组织厚度的线性关系计算在各角度下造影剂和透射拉曼装置中拉曼采集探头的距离,最后综合所有角度的拉曼光谱,利用误差最小的原则完成造影剂的准确定位。
进一步地,所述拉曼探针造影剂包括金纳米颗粒内核、拉曼报告分子中间层、银外壳、表面修饰层。
进一步地,所述表面修饰层可以为聚乙二醇层或氧化硅层。
进一步地,所述拉曼采集探头置于入射激光的对侧,接收的是透射拉曼光谱信号。
进一步地,所述步骤1,光源和样品之间的距离可调,以调整入射激光的直径,并选择大直径光斑以达到安全激光照射剂量。
进一步地,所述步骤2,沿圆柱体轴向扫描步长可调整。
进一步地,所述步骤2,对于扫描结果采用高斯曲线进行拟合,得到信号最强的位置,对应于造影剂所在的平面。
进一步地,所述步骤3,采用圆形转盘作为样品台,旋转柱状生物组织进行透射拉曼信号采集。
进一步地,所述步骤3,选择环形扫描的角度间隔,角度间隔的范围包括10度至30度。
进一步地,所述谱峰比对数值采用的是单个谱峰对,即两个谱峰的强度的比值取对数的数据。
在本发明的较佳实施方式中,本发明开发了一种方法,可以完成厚组织中拉曼造影剂的三维实时定位。本方法结合“两步法”扫描策略,利用采集到的透射拉曼光谱信号,实现了生物组织中拉曼造影剂的三维定位。定位步骤如下:
步骤1:将含有拉曼探针造影剂的柱状生物组织样品置于样品台上,采用拉曼光谱仪在横向进行透射检测,即光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号。
步骤2:沿圆柱体的轴向扫描,即同时调整光源和拉曼采集探头的高度,并且保持两者处于同一高度以构成透射检测装置,根据扫描结果拟合得到曲线,并取信号强度最大位置为造影剂所在高度位置z;
步骤3:对于步骤2选取的造影剂所在高度的平面进行环形扫描,即保持光源和拉曼采集探头位于对侧构成透射装置,通过同时旋转光源和拉曼采集探头或者仅旋转样品的方法,检测得到360°各个角度下的拉曼光谱。
步骤4:计算步骤3所得的拉曼光谱的谱峰比的对数值,并且根据“谱峰比对数值”-组织厚度的线性关系计算在某角度下造影剂和透射拉曼装置中拉曼采集探头的距离,最后综合所有角度的拉曼光谱,利用误差最小的原则完成造影剂的准确定位。
当前存在的病灶三维定位或者成像方法,一般采用的是螺旋扫描或者逐层环形扫描的方法,扫描过程较为繁琐。本发明开发了“两步法”扫描策略,仅需通过两步,即沿圆柱组织轴向扫描一次,确定造影剂所在平面,然后对于造影剂所在平面进行环形扫描,即可完成病灶的定位。
当前的断层扫描(tomography)方法重建图像,一般是基于信号强度信息,结合滤波反投影等方法完成图像的重构,整体算法较为复杂,需要更长的重建时间。本发明采用了基于“谱峰比对数值”这一包含深度信息的信号,结合最小误差原则,进行病灶位置乃至图像的重构。本发明整体算法较为简单,仅需通过很少的环形扫描角度(对于实施例1,仅需36个扫描角度的数据)就可以完成二维平面(x,y)中的病灶定位
当前的拉曼造影剂定位方法,多数只能在很高的激光照射剂量下测量信号,无法实现安全激光剂量下的造影剂定位。本发明采用了大光斑入射的策略,可以通过调整激光和样品之间的间距,让入射的激光光斑面积较大,从而在总入射功率较大的基础上,仍旧满足安全激光照射剂量要求(因为面积大,所以密度小)。光源和样品之间的距离可调,以此调整入射激光的直径,并选择大直径光斑以达到安全激光照射剂量。在实施例1中,本发明选择的入射激光功率为300mW,并选择激光光斑直径为0.5cm,计算得到激光的功率密度为1.528W/cm2,符合安全激光照射剂量。
拉曼探针造影剂包括金纳米颗粒内核、拉曼报告分子中间层、银外壳、表面修饰层。
表面修饰层可以为聚乙二醇层或氧化硅层。
透射检测方法即将所述拉曼采集探头置于入射激光的对侧,接收的是透射拉曼光谱信号。
步骤1包括:所述光源和样品之间的距离可调,以此调整入射激光的直径,并选择大直径光斑以达到安全激光照射剂量。
步骤2包括:沿圆柱体轴向扫描步长可调整。
步骤2还包括:对于扫描结果采用高斯曲线进行拟合,得到信号最强的位置,对应于造影剂所在的平面。
步骤3包括:采用圆形转盘作为所述样品台,旋转柱状生物组织进行透射拉曼信号采集。
步骤3还包括:选择环形扫描的角度间隔,角度间隔的范围包括10度至30度。
步骤4通过使用环形扫描采集得到的光谱数据,通过数据处理得到“谱峰比对数值”,由此得到在角度i下的造影剂距拉曼采集探头预测距离为计算公式为:
“谱峰比对数值”采用的是单个谱峰对,即两个谱峰的强度的比值取对数的数据。可以选择不同的拉曼谱峰分析所述的“谱峰比对数值”。
步骤4的最小误差的原则如下,设定在在角度i下的造影剂距拉曼采集探头预测距离为对于平面每一点进行误差计算,平面某点距拉曼采集探头实际距离为di,综合所有环形扫描角度的预测结果,平面某点总误差为:/>利用最小误差原则即:寻找造影剂最有可能的平面位置。
本发明与现有技术相比较,具有如下显而易见的实质性特点和显著优点:
1.本发明完成了生物组织内拉曼造影剂的精确三维定位,填补了定位厚生物组织中拉曼造影剂这一需求的空白,且达到了很高的精确度,在直径5cm的柱状猪肉组织中定位误差不超过2mm。
2.本发明大大节约了扫描需要的步骤,对于生物组织中的造影剂定位,仅需通过沿轴向的扫描,以及环形扫描两步即可完成造影剂的三维定位。
3.本发明整体算法较为简单,仅需通过很少的环形扫描角度就可以完成二维平面中的病灶定位。
4.本发明符合安全激光照射剂量。
以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明,以充分地了解本发明的目的、特征和效果。
附图说明
图1是本发明的一个较佳实施例的基于透射拉曼光谱和“两步法”扫描的造影剂定位方法的过程示意图;
图2是本发明的一个较佳实施例的造影剂定位方法的流程图;
图3是本发明的一个较佳实施例的拉曼探针造影剂的透射电镜表征图;
图4是本发明的一个较佳实施例的拉曼探针造影剂的拉曼光谱图;
图5是本发明的一个较佳实施例的基于透射拉曼光谱和“两步法”扫描的造影剂定位方法对厚组织内造影剂凝胶样品进行(x,y)平面空间定位的结果图;
图6是本发明的一个较佳实施例的厚组织内造影剂凝胶样品分布的实际(x,y)位置的照片。
具体实施方式
以下参考说明书附图介绍本发明的多个优选实施例,使其技术内容更加清楚和便于理解。本发明可以通过许多不同形式的实施例来得以体现,本发明的保护范围并非仅限于文中提到的实施例。
在附图中,结构相同的部件以相同数字标号表示,各处结构或功能相似的组件以相似数字标号表示。附图所示的每一组件的尺寸和厚度是任意示出的,本发明并没有限定每个组件的尺寸和厚度。为了使图示更清晰,附图中有些地方适当夸大了部件的厚度。
如图1所示,造影剂定位方法的原理具体说明如下:
对于拉曼“谱峰比对数值”和造影剂包埋深度的线性关系推导如下:
首先是激光入射过程。在透射拉曼装置中,假设激光的传播距离大于在生物组织中的传输平均自由路径,那么可以在三维笛卡尔坐标系中得到激光透射穿透过程中的稳态扩散近似模型:
其中I(l)是激光传播距离为l后的光强,I0是激光初始强度,a是激光光斑入射半径,l是入射激光传播的距离,λ0是入射激光波长,μeff(λ0)是激光传播过程中的有效衰减系数,D(λ0)是激光光子扩散系数。
以上过程描述了激光入射组织后发生的衰减。其次是拉曼颗粒经由激光激发后产生拉曼光子并传输的过程。拉曼颗粒产生的拉曼光由强度为I(l)的激光激发,并与拉曼颗粒的浓度C、拉曼散射截面σ共同决定。产生的拉曼光子再经由组织衰减后,强度可以利用稳态传输方程描述如下:
其中λi是拉曼光子的波长,不同拉曼位移的光子具有不同的波长,μeff(λi)是拉曼光子的有效衰减系数,D(λi)是拉曼光子的扩散系数,d是拉曼颗粒(信号源)和拉曼采集探头(接收端)之间的距离,定义为造影剂包埋的深度d。
结合公式(1)和公式(2),并对于公式进行化简,检测到的拉曼光强为:
其中A定义如下,是拉曼光子波长无关项:
对于不同波长的拉曼光子,均可以用公式(3)描述其在透射***中光子的传输过程,而不同拉曼位移的拉曼光子也会有不同的有效衰减系数和扩散系数。当选定拉曼光子的波长为λ1和λ2时,其有效衰减系数和光子扩散系数都能够确定,检测到的拉曼光强为分别利用公式(3)描述,并作商得:
定义两个波长拉曼光子的有效衰减系数差为Δμeff:
Δμeff=μeff(λ2)-μeff(λ1) (6)
公式(5)两边取对数,结合公式(6),不同波长拉曼光子的“谱峰比对数值”最终可以简化为如下线性模型:
其中与激光传播相关的项都被约去,“谱峰比对数值”与造影剂包埋深度d成线性关系。且该线性模型斜率Δμeff仅与拉曼光子有效衰减系数差有关,该模型截距可以认为是拉曼颗粒不经过组织散射得到的拉曼信号的“谱峰比对数值”,是拉曼颗粒的一个特性,可以在颗粒完成开发并合成后确定,不会随测试条件变化而变化。
对于如何利用最小误差原则,结合线性关系实现造影剂的三维定位,具体说明如下:
首先,对于每个扫描角度来说,都可以确定一个拉曼采集探头所在位置,这在扫描过程中是已知的。定义拉曼采集探头的位置坐标为 其中r为柱状组织的半径,θi为某扫描角度对应的拉曼采集探头所处极角。利用所述的线性关系,可以得到某扫描角度下,拉曼采集探头和造影剂之间的距离,定义为造影剂预测深度/>这些数据均可通过扫描过程直接或者间接获取。
其次,对于造影剂所在组织断层平面内任意一点,位置坐标t=(xt,yt),定义误差函数E(t,θ)如下:
其中代表某一扫描角度下组织平面点和拉曼采集探头的实际距离(取欧氏距离),/>表示拉曼采集探头-组织平面点的预测距离(通过所述线性关系得到),作差即为误差。n表示扫描的角度总数,即对于多角度光谱信号误差取平均。该误差函数的含义是:组织平面内某点对于多角度光谱信号预测深度误差情况的平均。为了方便理解,假定如下理想场景:选择的组织平面点即为造影剂所在位置t0=(x0,y0),且所有的/>预测都绝对准确,即对于任意θi,有/> 那么此时得到的误差函数值应该为0。实际情况下,因为造影剂具有大小且光谱具有噪声,实际误差不可能为0,那么定位造影剂中心位置采用最小化误差的方法,即取误差函数最小值的组织平面位置,认为该位置gpre是造影剂中心位置的可能性最大,计算公式如下:
由此实现造影剂的定位。
基于所述原理,本发明提供了一种基于透射拉曼检测方式以及“两步法”扫描策略的造影剂定位方法,如图2所示,包括以下步骤:
步骤一:制备得到含有IR-780拉曼信号分子的核壳结构表面增强拉曼探针5mL0.6nmol/L。将上述获得的表面增强拉曼探针纳米颗粒与琼脂糖溶液混合成凝胶后,埋在柱状猪肉组织内部。将包含有拉曼凝胶的组织样品置于样品台上,采用透射拉曼装置在横向进行检测,光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号。
步骤二:首先沿圆柱体的轴向扫描,即同时调整光源和拉曼采集探头的高度,并且保持两者处于同一高度以构成透射装置,根据扫描结果拟合得到拉曼“谱峰比对数值”和扫描高度z的关系曲线,并取信号强度最大位置为造影剂所在高度位置z;
步骤三:根据步骤三选取的造影剂所在高度z进行环形扫描即保持光源和拉曼采集探头位于对侧构成透射装置,通过同时旋转光源和拉曼采集探头或者仅旋转样品的方法,检测不同角度下的拉曼光谱;
步骤四:计算所有光谱的谱峰比的对数值,获得该位置与造影剂的距离,综合计算出造影剂的(x,y)二维平面中误差函数的值,并通过最小误差原则,选择误差最小的点作为造影剂最有可能所处的点,最终完成造影剂定位。
实施例1
在该实施例中,如图3所示是核壳结构表面增强拉曼探针的透射电镜表征图,如图4所示是核壳结构表面增强拉曼探针的拉曼光谱图,如图5所示是本例中造影剂的平面空间定位的结果图,如图6所示是本例中造影剂的真实平面空间位置。
步骤一:制备得到含有IR-780拉曼信号分子的核壳结构表面增强拉曼探针5mL0.6nmol/L,将上述获得的表面增强拉曼探针纳米颗粒与琼脂糖溶液混合成凝胶后,埋在直径5cm的柱状猪肉组织内部。将包含有拉曼凝胶的组织样品置于样品台上,采用透射拉曼装置在横向进行检测,光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号。选择的入射激光功率为300mW,并选择激光光斑直径为0.5cm,计算得到激光的功率密度为1.528W/cm2,符合安全激光照射剂量。
步骤二:首先沿圆柱体的轴向扫描,即利用平移台以0.5cm的步长同时垂直调整光源和拉曼采集探头的高度,并且保持两者处于同一高度以构成透射检测装置,根据扫描结果拟合得到拉曼“谱峰比对数值”和扫描高度z的关系曲线,并取信号强度最大位置为造影剂所在高度位置z;
步骤三:根据步骤三选取的造影剂所在高度z进行环形扫描,即保持光源和拉曼采集探头位于对侧构成透射装置,同时旋转猪肉样品改变检测角度,共扫描36个点,间隔10度。
步骤四:计算所有光谱的谱峰比的对数值,获得该位置与造影剂的距离,综合计算出造影剂的(x,y)二维平面中误差函数的值。计算完成后,取最小误差函数点的105%作为阈值,并取阈值以下的值,绘制热图如图5所示,以此展示造影剂出现的可能性比较大的位置。并通过最小误差原则,选择误差最小的点作为造影剂最有可能所处的点,最终完成造影剂定位。
根据实验测试的光谱数据,使用表面增强拉曼探针内的拉曼信号分子IR-780的两个拉曼特征峰(525和1203cm-1)的强度,对于误差函数进行计算后,重构出误差函数二维热图如图5所示,可实现对拉曼探针的(x,y)分布的定位。实际的(x,y)分布如图6照片所示。选择误差最小的点作为造影剂最有可能所处的点,最终得到的预测误差,即预测位置和造影剂真实位置之间的距离,仅为1.3mm。这种方法提升了成像深度,可以获得无创非侵入性的检测和定位效果。
实施例2
步骤一:制备得到含有IR-780拉曼信号分子的核壳结构的表面增强拉曼探针5mL0.6nmol/L,将上述获得的表面增强拉曼探针纳米颗粒与琼脂糖溶液混合成凝胶后,分别取两块埋在直径7cm的柱状猪肉组织内部。将包含有拉曼凝胶的组织样品置于样品台上,采用透射拉曼装置在横向进行检测,光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号。选择的入射激光功率为300mW,并选择激光光斑直径为0.5cm,计算得到激光的功率密度为1.528W/cm2。
步骤二:首先沿圆柱体的轴向扫描,即利用平移台以0.5cm的步长同时垂直调整光源和拉曼采集探头的高度,并且保持两者处于同一高度以构成透射检测装置,根据扫描结果拟合得到拉曼“谱峰比对数值”和扫描高度z的关系曲线,并取信号强度极值的两个位置为造影剂所在高度位置z1,z2;
步骤三:根据步骤三选取的造影剂所在高度z1,z2进行环形扫描,即保持光源和拉曼采集探头位于对侧构成透射装置,同时旋转猪肉样品改变检测角度,各扫描18个点,间隔20度。
步骤四:分别计算两个平面环形扫描的所有光谱的“谱峰比对数值”,由此获得该位置与造影剂的距离,综合计算出造影剂的(x,y)二维平面中误差函数的值,并通过最小误差原则,选择误差最小的点作为造影剂最有可能所处的点,最终完成两块造影剂定位。
实施例3
步骤一:制备得到含有IR-780拉曼信号分子的核壳结构的表面增强拉曼探针5mL0.6nmol/L,将上述获得的表面增强拉曼探针纳米颗粒与琼脂糖溶液混合成凝胶后,埋在腹部直径约5cm的大鼠体内。将包含有拉曼凝胶的大鼠置于样品台上,采用透射拉曼装置在横向进行检测,光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号。选择的入射激光功率为300mW,并选择激光光斑直径为0.5cm,计算得到激光的功率密度为1.528W/cm2。
步骤二:首先沿圆柱体的轴向扫描,即利用平移台以0.5cm的步长同时垂直调整光源和拉曼采集探头的高度,并且保持两者处于同一高度以构成透射检测装置,根据扫描结果拟合得到拉曼“谱峰比对数值”和扫描高度z的关系曲线,并取信号强度最大位置为造影剂所在高度位置z;
步骤三:根据步骤三选取的造影剂所在高度z进行环形扫描,即保持光源和拉曼采集探头位于对侧构成透射装置,同时旋转猪肉样品改变检测角度,共扫描12个点,间隔30度。
步骤四:计算所有光谱的“谱峰比对数值”,并由此获得该位置与造影剂的距离,综合计算出造影剂的(x,y)二维平面中误差函数的值,并通过最小误差原则,选择误差最小的点作为造影剂最有可能所处的点,最终完成造影剂在大鼠中的定位。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (10)
1.一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,包括以下步骤:
步骤1、将含有拉曼探针造影剂的柱状生物组织样品置于样品台上,采用拉曼光谱仪在横向进行透射检测,即光源从左入射,采用拉曼采集探头在右侧收集拉曼光谱信号;
步骤2、沿圆柱体的轴向扫描,即同时调整光源和拉曼采集探头的高度,并且保持处于同一高度以构成透射检测装置,根据扫描结果拟合得到曲线,并取信号强度最大位置为造影剂所在高度位置;
步骤3、对于所述步骤2选取的造影剂所在高度的平面进行环形扫描,即保持光源和拉曼采集探头位于对侧构成透射装置,通过同时旋转光源和拉曼采集探头或者仅旋转样品的方法,检测得到360°各个角度下的拉曼光谱;
步骤4、计算所述步骤3所得的拉曼光谱的谱峰比的对数值,并且根据谱峰比对数值-组织厚度的线性关系计算在各角度下造影剂和透射拉曼装置中拉曼采集探头的距离,最后综合所有角度的拉曼光谱,利用误差最小的原则完成造影剂的准确定位。
2.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述拉曼探针造影剂包括金纳米颗粒内核、拉曼报告分子中间层、银外壳、表面修饰层。
3.如权利要求2所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述表面修饰层可以为聚乙二醇层或氧化硅层。
4.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述拉曼采集探头置于入射激光的对侧,接收的是透射拉曼光谱信号。
5.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述步骤1,光源和样品之间的距离可调,以调整入射激光的直径,并选择大直径光斑以达到安全激光照射剂量。
6.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述步骤2,沿圆柱体轴向扫描步长可调整。
7.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述步骤2,对于扫描结果采用高斯曲线进行拟合,得到信号最强的位置,对应于造影剂所在的平面。
8.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述步骤3,采用圆形转盘作为样品台,旋转柱状生物组织进行透射拉曼信号采集。
9.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述步骤3,选择环形扫描的角度间隔,角度间隔的范围包括10度至30度。
10.如权利要求1所述的基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法,其特征在于,所述谱峰比对数值采用的是单个谱峰对,即两个谱峰的强度的比值取对数的数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311041135.8A CN117064336A (zh) | 2023-08-17 | 2023-08-17 | 一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311041135.8A CN117064336A (zh) | 2023-08-17 | 2023-08-17 | 一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117064336A true CN117064336A (zh) | 2023-11-17 |
Family
ID=88718998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311041135.8A Pending CN117064336A (zh) | 2023-08-17 | 2023-08-17 | 一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117064336A (zh) |
-
2023
- 2023-08-17 CN CN202311041135.8A patent/CN117064336A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lin et al. | The emerging role of photoacoustic imaging in clinical oncology | |
Maloney et al. | Review of methods for intraoperative margin detection for breast conserving surgery | |
Kepshire et al. | A microcomputed tomography guided fluorescence tomography system for small animal molecular imaging | |
US8314406B2 (en) | Systems and methods for optical imaging using early arriving photons | |
JP5701615B2 (ja) | 電磁トラッキング及び光針による生検誘導 | |
US6503478B2 (en) | Chemically specific imaging of tissue | |
US8243269B2 (en) | Raman spectral analysis of sub-surface tissues and fluids | |
CN101959450B (zh) | 通过基于图像的x射线引导*** | |
JP5543212B2 (ja) | 前立腺癌のイメージング用システム、前立腺癌をイメージングするシステムの作動方法及びコンピュータ読取可能媒体 | |
US20050070803A1 (en) | Multiphoton photoacoustic spectroscopy system and method | |
CN107645938A (zh) | 成像导引的敞开式电离质谱测定 | |
CN104067313B (zh) | 成像装置 | |
WO2007111669A2 (en) | Combined x-ray and optical tomographic imaging system | |
US20180042483A1 (en) | Systems and methods for hyperspectral imaging | |
Carpenter et al. | Limited-angle x-ray luminescence tomography: methodology and feasibility study | |
US20150320319A1 (en) | Near infrared photonic prostatoscopy analyzer | |
Fixler et al. | Subcutaneous gold nanorods detection with diffusion reflection measurement | |
Chernomordik et al. | Inverse method 3-D reconstruction of localized in vivo fluorescence-application to Sjogren syndrome | |
US7809431B2 (en) | Method of optically imaging biological tissues by using fluorescence, in particular for defining regions of interest in tissues to be analyzed by tomography | |
CN117064336A (zh) | 一种基于透射拉曼光谱及两步法扫描拉曼造影剂定位方法 | |
EP1797818A2 (en) | Method and system for tomographic imaging using fluorescent proteins | |
Shiryaev et al. | Investigated spectral-fluorescent properties of endogenous porphyrins of the wild boar hepatobiliary system optimize the diagnostics and treatment of cholangiocarcinoma with FD and PDT | |
Grabtchak et al. | Radiance spectroscopy tool box for characterizing Au nanoparticles in tissue mimicking phantoms as applied to prostate | |
El-Sharkawy | Breast cancer detection by determination of optical properties of non-malignant and malignant breast tissues | |
Wallace et al. | Biomedical applications of THz imaging |
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 |