CN112089404A - 一种内窥OCT-Raman双模成像装置及成像方法 - Google Patents
一种内窥OCT-Raman双模成像装置及成像方法 Download PDFInfo
- Publication number
- CN112089404A CN112089404A CN202011075960.6A CN202011075960A CN112089404A CN 112089404 A CN112089404 A CN 112089404A CN 202011075960 A CN202011075960 A CN 202011075960A CN 112089404 A CN112089404 A CN 112089404A
- Authority
- CN
- China
- Prior art keywords
- oct
- raman
- imaging
- endoscopic
- cavity
- 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
- 238000001069 Raman spectroscopy Methods 0.000 title claims abstract description 132
- 238000003384 imaging method Methods 0.000 title claims abstract description 123
- 238000003332 Raman imaging Methods 0.000 claims abstract description 39
- 238000001237 Raman spectrum Methods 0.000 claims abstract description 16
- 230000003287 optical effect Effects 0.000 claims abstract description 13
- 239000000523 sample Substances 0.000 claims description 65
- 239000013307 optical fiber Substances 0.000 claims description 56
- 230000033001 locomotion Effects 0.000 claims description 22
- 230000005284 excitation Effects 0.000 claims description 20
- 230000001681 protective effect Effects 0.000 claims description 19
- 210000004712 air sac Anatomy 0.000 claims description 11
- 230000004927 fusion Effects 0.000 claims description 8
- 238000000034 method Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000009877 rendering Methods 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000002902 bimodal effect Effects 0.000 claims 3
- 230000009977 dual effect Effects 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000003595 spectral effect Effects 0.000 abstract description 8
- 230000003902 lesion Effects 0.000 abstract description 5
- 238000012800 visualization Methods 0.000 abstract description 4
- 238000012014 optical coherence tomography Methods 0.000 description 108
- 239000000835 fiber Substances 0.000 description 14
- 210000004204 blood vessel Anatomy 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004195 computer-aided diagnosis Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 210000001035 gastrointestinal tract Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- FGYHLIHEFHCMPP-FVWOBLMASA-N rhosin Chemical group N1=CC=NC2=CC(/C=N/NC(=O)[C@@H](CC=3C4=CC=CC=C4NC=3)N)=CC=C21 FGYHLIHEFHCMPP-FVWOBLMASA-N 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000002485 urinary effect Effects 0.000 description 1
- 230000000007 visual effect Effects 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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B1/00—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
- A61B1/06—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor with illuminating arrangements
- A61B1/0646—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor with illuminating arrangements with illumination filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B1/00—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
- A61B1/06—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor with illuminating arrangements
- A61B1/07—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor with illuminating arrangements using light-conductive means, e.g. optical fibres
-
- 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Signal Processing (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Optics & Photonics (AREA)
- Power Engineering (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Endoscopes (AREA)
Abstract
本发明公开了一种内窥OCT‑Raman双模成像装置及成像方法，包括拉曼成像系统，OCT成像系统、内窥镜系统、图像重建系统和显示系统。内窥镜系统分别与拉曼成像系统和OCT成像系统相连，拉曼成像系统和OCT成像系统皆与图像重建系统相连，图像重建系统与显示系统相连。本发明能够获得较高光学分辨率的三维图像，改变了拉曼光谱的谱线形式存在方式，将拉曼光谱与OCT三维成像相关联，实现了OCT腔体平铺三维图像及其对应的拉曼三维图像的重建，使病变位置对应更加明显，提高了成像的可视化程度，实现OCT成像数据和拉曼信号准确快速的采集。
Description
技术领域
本发明属于内窥镜技术领域，尤其涉及一种内窥OCT-Raman双模成像装置及成像方法。
背景技术
内窥镜作为一种无创或微创的成像方法，能够深入生物体内部体腔成像，已经广泛应用于血管，消化道，泌尿系统以及呼吸道等多个腔体的影像诊断。OCT是一种利用生物组织散射光相干原理成像的介观(微米尺度)活体组织高分辨率成像和观测手段，OCT内窥镜是内窥镜技术的重要分支。OCT技术有较高的分辨率的同时，又具有较大的层析能力，并且可以很容易地与光纤导管及内窥镜联用，对人体内器官进行成像。尤其是内窥频域OCT的能够实现高速实时监测、高分辨率、无损断层成像，同时能够对生物组织的形态结构、化学成分以及生理功能信息进行三维成像，在腔体类疾病和癌症的早期诊断方面具有重要的价值。
然而单纯的OCT内窥成像，其分析诊断难度高，需要大量的先验知识和临床经验。为了进一步提高诊断依据，可以加入从分子级上提供生物组织特征信息的拉曼检测。实际应用中，在OCT内窥探头内，加入拉曼激发和采集光纤，可以在获取OCT图像的同时，获得腔体表面的拉曼数据，该数据不仅可以作为诊断依据，还提升了计算机辅助诊断能力。
现有的OCT-Raman内窥探头，存在以下问题：首先，现有的拉曼内窥系统尽管能够获得内窥图像及其对应位置的拉曼光谱，但其拉曼光谱本身以谱线形式存在，需要复杂的数据分析才能够实现对病变的分析，与成像的内窥图像关联程度低，位置对应不明显，可视化程度较低，不能明显判断异常和病变部位。其次，在腔体内窥系统中，内窥探头的成像位置受成像导管在腔体内的位置所限制，其成像中心并不在血管中心，虽然这对内窥成像没有太大影响，但是对于拉曼检测来说，检测距离很大程度上影响了拉曼信号的强度。最后，大多数的内窥导管的机械旋转扫描是通过扭矩线缆实现的，在较大长度的腔体检测中，存在转速不均导致分均匀图像失真的问题，其转速和稳定性都受到一定程度的限制。
发明内容
本发明提供了一种内窥OCT-Raman双模成像装置及成像方法，针对腔体疾病的诊断，能够获得较高光学分辨率的三维图像的同时，改变了拉曼光谱的谱线形式存在方式，将拉曼光谱与OCT三维成像相关联，实现OCT腔体平铺三维图像及其对应的拉曼三维图像的重建，使病变位置对应更加明显，提高了成像的可视化程度；同时采用了微电机控制的探头，辅以外部气囊，实现OCT成像数据和拉曼信号准确快速的采集。
本发明提供的技术方案为：
一种内窥OCT-Raman双模成像装置，包括：拉曼成像系统、OCT成像系统、内窥镜系统、图像重建系统和显示系统；所述内窥镜系统分别与拉曼成像系统和OCT成像系统相连，所述拉曼成像系统和OCT成像系统皆与图像重建系统相连，所述图像重建系统与显示系统相连。
进一步的，所述拉曼成像系统还包括拉曼激发光光源装置、拉曼光谱仪，所述拉曼激发光光源装置和拉曼光谱仪分别通过光纤与内窥镜系统连接，拉曼光谱仪和图像重建系统连接。
进一步的，所述的OCT成像系统还包括激光光源装置、OCT干涉成像模块和数据采集模块；所述激光光源和OCT干涉成像模块分别通过光纤与内窥镜系统相连，数据采集模块与OCT干涉成像模块相连，数据采集模块与图像重建系统连接。
进一步的，所述的内窥镜系统包括运动控制模块和内窥探头，所述运动控制模块包括依次相连的光纤连接器、运动控制装置和成像导管，内窥探头包括光纤束、滤光片、反射镜、反射镜座、微型电机、连接管，光纤保护套、气囊导管、外保护套和气囊；
其中，光纤束和滤光片由光纤保护套相连，光纤束包括拉曼成像激发光纤、拉曼信号收集光纤、OCT光源光纤和OCT反射光纤；其中拉曼激发光纤、OCT光源光纤和OCT反射光纤呈三角形结构紧密排列组成中心光纤束，拉曼信号收集光纤围绕中心光纤束环形布置；反射镜固定在反射镜座上，反射镜座与连接管连接，连接管与光纤保护套由微型电机相连，所述连接管通微型电机作用，能够围绕光纤束旋转。微型电机由运动控制装置控制，微型电机包括空心定子和转子，空心定子与光纤保护套共轴连接，转子与连接管共轴连接，转子转动带动与连接管连接的反射镜转动；气囊导管与气囊相接，气囊导管固定在光纤保护套上，气囊随导管移动，气囊导管可调整气囊内气压。
进一步的，所述图像重建系统还包括OCT重建系统和拉曼重建系统；所述OCT重建系统与数据采集模块连接，用于将OCT原始腔体二维图像转换为OCT腔体平铺二维图像，再利用OCT三维重建算法，将一系列内窥探头轴向移动获得的数据，生成OCT腔体平铺三维图像，所述拉曼重建系统与拉曼光谱仪连接，用于生成与OCT腔体平铺三维图像对应的拉曼三维图像。
一种内窥OCT-Raman双模成像方法：
步骤1：通过拉曼成像系统输出拉曼激发光源，经过内窥镜系统照射管腔样本，OCT成像系统的输出OCT成像光源照射经过内窥镜系统照射管腔样本；
步骤2：管腔样本受到激发光照射，拉曼信号经内窥镜系统采集输入拉曼成像系统，OCT信号被内窥镜系统反射进入OCT成像系统；
步骤3：拉曼成像系统根据管腔样本的拉曼信号生成拉曼光谱送入图像重建系统，OCT成像系统根据OCT干涉成像模块生成OCT图像数据送入图像重建系统；
步骤4：每采集完一次信号，内窥探头转动一定角度重新采集信号，重复至旋转一圈；每旋转扫描一圈，内窥探头轴向移动一定距离重新采集信号，重复至完成轴向扫描。
步骤5：OCT成像系统通过极坐标展开算法，使OCT原始腔体二维图像转换为OCT腔体平铺二维图像。
步骤6：利用OCT三维重建算法，将一系列探头轴向移动获得的OCT原始腔体二维图像，生成OCT腔体平铺三维图像。
步骤7：拉曼成像系统根据拉曼三维成像算法生成与OCT腔体平铺三维图像对应的拉曼三维图像
进一步的，步骤5所述的极坐标展开算法为，利用OCT原始腔体二维图像寻找血管中心，重新建立以血管中心为原点，原X轴为极轴的极坐标系，通过图像坐标变换，将极坐标系的像素点变换到直角坐标系下，生成OCT腔体平铺二维图像；
进一步的，步骤6所述的三维重建算法包括体绘制算法和光线追踪算法；
进一步的，步骤7的所述拉曼三维重建算法为：拉曼成像系统根据管腔样本不同位置处不同拉曼位移的信号强度，获得表征腔体内壁不同物质的拉曼峰值，根据设置的拉曼位移，获取与OCT腔体平铺三维图像中对应腔体内壁位置的拉曼强度值，将该拉曼强度值作为纵坐标，XY轴与OCT腔体平铺三维图像对应，构成拉曼三维图像；
进一步的，步骤6所述的OCT腔体平铺三维图像与步骤7所述的拉曼三维图像通过图像融合算法可实现融合显示。
本发明具有如下优点或有益效果：
本发明所提供的内窥OCT-Raman双模成像装置设计了有可变气囊的、可高速旋转的内窥探头，使内窥探头的成像中心始终处在血管中心，保证了拉曼信号的准确性，解决了因探头转速不均导致的图像失真问题，通过同时将拉曼激发光和激光光源经内窥镜系统定位照射到样本腔体，采集拉曼信号和OCT成像反射信号的方法实现了双模数据的采集，通过极坐标展开算法形成腔体平铺图像，实现了OCT极坐标三维成像，利用拉曼三维重建算法实现拉曼三维成像，将两图像相互关联，提高了拉曼信息的可视化程度，使病变位置对应更加明显。
附图说明
图1是本发明实施例的内窥OCT-Raman双模成像装置的整体结构示意图。
图2是本发明实施例的内窥OCT-Raman双模成像装置的连接示意图。
图3是本发明实施例的运动控制模块的结构示意图。
图4是本发明实施例的内窥探头的结构示意图。
图5是本发明实施例的内窥探头多模光纤结构的示意图。
图6是本发明实施例的内窥OCT-Raman双模成像方法的流程图。
图7是本发明实施例的内窥OCT干涉成像模块结构示意图
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
图1是本发明实施例的内窥OCT-Raman双模成像装置的整体结构示意图，图2是本发明实施例的内窥OCT-Raman双模成像装置的连接示意图。在本实施方式中，内窥OCT-Raman双模成像装置包括拉曼成像系统10、OCT成像系统20、内窥镜系统30、图像重建系统40和显示系统50。内窥镜系统30分别与拉曼成像系统10和OCT成像系统20相连，拉曼成像系统10和OCT成像系统20皆与图像重建系统40相连，图像重建系统40与显示系统50相连；其中拉曼成像系统10包括拉曼激发光光源装置110和拉曼光谱仪120，OCT成像系统包括激光光源装置210和OCT干涉成像模块220和数据采集模块230，内窥镜系统30包括运动控制模块310和内窥探头320，成像导管313，图像重建系统40包括拉曼重建系统410和OCT重建系统420；
图3是本发明实施例的运动控制模块的结构示意图；图4是本发明实施例的内窥探头的结构示意图；图5是本发明实施例的内窥探头多模光纤结构的示意图。其中，运动控制模块310包括光纤连接器311、运动控制装置312和成像导管313，内窥探头320包括光纤束3201、滤光片3202、反射镜3203、反射镜座3204、微型电机3205、连接管3206，光纤保护套3207、气囊导管3208、外保护套3209和气囊3210；微型电机3205由空心定子3211和转子3212组成；光纤束3201包括拉曼成像激发光纤32011，OCT成像光源光纤32012，OCT成像反射光纤32013，拉曼信号收集光纤32014；光纤束3201和滤光片3202由光纤保护套3207相连，反射镜3203固定在所述反射镜座3204上与所述连接管3206连接，连接管3206与光纤保护套3207由微型电机3205相连，连接3206管通过微型电机3205作用，可围绕光纤束3201高速旋转，气囊导管3208与气囊3210相接，气囊导管3208固定在光纤保护套3207上，气囊3210随导管3206移动，可调整气囊3210内气压。
本发明实施例的内窥OCT-Raman双模成像装置的工作原理为：拉曼激发光光源装置110向内窥镜系统30输出激发光进入拉曼成像激发光纤32011，激光光源装置210向OCT干涉成像模块220输出扫频光源，OCT干涉成像模块220将光分为两路，一路进入干涉成像模块220的参考臂223，经反射后，作为参考光，另一路进入内窥镜系统30的OCT成像光源光纤32012；内窥镜系统中的运动控制装置312控制成像导管轴向运动到样本腔体的成像部位，内窥探头320的拉曼成像激发光纤32011和OCT光源光纤32012照射官腔样本成像部分，通过拉曼信号收集光纤32014收集该部位的拉曼信号输入到拉曼光谱仪120，生成拉曼光谱信号后，输入拉曼重建系统410中，与此同时，通过OCT成像反射光纤32013反射OCT成像的光信号，反射管腔样本的光信号进入OCT干涉成像模块220，经光纤耦合器222耦合后，进入探测器221，探测器将输入光信号转变为电信号输入数据采集模块230，经过傅里叶变换，输入OCT重建系统420；通过运动控制装置312发出控制信号，控制内窥探头320内部旋转设备进行360度的旋转，再发出控制信号，控制内窥探头320轴向移动扫描，每采集完一次信号，内窥探头320转动一定角度重新采集信号，重复至旋转一圈；每旋转扫描一圈，内窥探头320轴向移动一定距离重新采集信号，重复至完成轴向扫描。图像重建系统40根据接收的信号，利用OCT的极坐标三维重建算法和拉曼三位重建算法，完成OCT腔体平铺三维图像和拉曼三维图像。
上述过程中，内窥探头320受到运动控制装置312的轴向移动信号后，气囊导管3208首先减小气囊3210压力，运动控制装置312轴向拉动光纤保护套3207，气囊导管3208固定在光纤保护套3207上，一同拖动气囊3210轴向运动，运动结束后增大气囊3210压力，使成像导管313位于血管横切面的几何中心；当内窥探头320需要转动时，运动控制装置312给与微型电机3205一定频率的转动信号，转子3212带动连接管3206上的反射镜3203转动，从而采集管腔360度范围的信号。
(其中OCT干涉成像模块220的干涉光谱信号可以表示为：
其中，Pr为探测器探测到的参考光功率，Po为管腔样本臂射出的照射到管腔样本的光功率，k(t)为波数，是时间的因变量，(z)表示扫频光源的相干函数。其中第三积分项是扫频OCT系统干涉光谱信号中的有用信息，对该项进行傅立叶变换后得到的信息就是从波数空间转换到深度坐标的结果，从而获得管腔样本轴向各层的散射强度信息，再对结果进行数据处理和图像重建就能够实现对管腔样本深度结构成像。
对上式表示的干涉光谱信号做傅立叶变换，其结果为：
其中，S(k)为扫频光源的功率谱密度函数，δ(z)为零光程差附近直流项，a(z)为管腔样本轴向坐标z处返回的光的幅度)。
图6是本发明实施例的内窥OCT-Raman双模成像方法的流程图，本发明实施例的内窥OCT-Raman双模成像方法包括以下步骤：
步骤1：通过拉曼成像系统10输出拉曼激发光源，经过内窥镜系统30照射管腔样本，OCT成像系统20输出OCT成像光源经过内窥镜系统30照射管腔样本；
在步骤1中，拉曼成像系统10和OCT成像系统20包括产生输出光源的激光器及滤波器组成的光路系统；
步骤2：管腔样本受到激发光照射，拉曼信号经内窥镜系统30采集输入拉曼成像系统10，OCT信号被内窥镜系统30反射进入OCT成像系统20；
步骤3：拉曼成像系统10根据管腔样本的拉曼信号生成拉曼光谱送入图像重建系统40，OCT成像系统20根据OCT干涉成像模块220生成OCT图像数据送入图像重建系统40；
在步骤3中，拉曼成像系统10使用拉曼光谱仪120产生拉曼光谱，OCT干涉成像模块220使用扫频OCT成像装置，OCT信号进入光纤耦合器222耦合后进入探测器221，被探测器中的电荷耦合元件捕获光信号转变为电信号送入数据采集模块230，数据采集模块230经过傅里叶变换算法生成OCT原始腔体二维图像；
步骤4：每采集完一次信号，内窥探头320转动一定角度重新采集信号，重复至旋转一圈；每旋转扫描一圈，内窥探头320轴向移动一定距离重新采集信号，重复至完成轴向扫描；
步骤5：OCT成像系统20通过极坐标展开算法，使OCT原始腔体图像转换为OCT腔体平铺二维图像；
在步骤5中，极坐标展开算法为：首先对图像进行去噪预处理后，利用霍夫变换算法，将OCT原始腔体二维图像的血管找出，并拟合为圆形，确定该圆心位置；重新建立以血管中心为原点，原X轴正方向为极轴的极坐标系，通过图像坐标变换公式x＝ρcosθ，y＝ρsinθ，将极坐标系的像素点变换到直角坐标系下，生成OCT腔体平铺二维图像。
步骤6：利用OCT三维重建算法，将一系列内窥探头320轴向移动获得的OCT原始腔体二维图像，生成OCT腔体平铺三维图像；
在步骤6中，三维重建算法包括体绘制算法和光线追踪算法，基本原理为通过内窥探头的OCT二维成像数据及其轴向移动大小将轴向移动生成的一系列OCT原始腔体二维图像构造成三维体数据，再根据体绘制算法或更为复杂的光线追踪算法在显示器中实现三维重建。
步骤7：拉曼成像系统10根据拉曼三维成像算法生成与OCT腔体平铺三维图像对应的拉曼三维图像。
在步骤7中，拉曼三维重建通过采集样本某点处的拉曼信号，选择指定谱峰的强度、谱峰面积或者谱峰宽等参数并以灰度图或伪色图的形式显示，图像明暗、色彩的变化与指定谱峰所代表的物质成分相对应，可提供物质的浓度、分布状况等生物信息，具体实现步骤包括：
步骤7.1展开：将腔体一周的拉曼信号，与腔体内壁位置一一对应，使用极坐标转换算法，使拉曼信号与腔体平铺图像各点一一对应。
步骤7.2预处理：使用滤波和拟合算法，完成拉曼光谱的去尖峰、去噪和基线校正，解决实际拉曼信号采集中受仪器、环境以及荧光影响导致采集的拉曼信号有较大的噪声与背景漂移等问题。
步骤7.3二维成像：根据腔体成分的观察需要，设定固定的拉曼位移，获取与OCT腔体平铺三维图像中对应腔体内壁位置的拉曼强度值，该拉曼强度值作为拉曼二维成像的Z轴坐标，拉曼三维成像的X轴与OCT腔体平铺三维图像X轴对应，形成表征物质成分变化的拉曼二维图像；并且拉曼强度值不仅可以与作为坐标值，还可以匹配对应颜色表，使其高度变换更加明显。
步骤7.4插值：根据已知的拉曼二维成像像素点，选择合适的插值函数进行插值，增加拉曼成像的像素点，改善其视觉效果；对于拉曼成像，需考虑激发光束特点和采集点处采集信息的真实分布情况，合理使用最近邻域法、双线性内插法和基于小波变换的图像插值法等方法。
步骤7.5三维成像：将一系列探头轴向移动获得的拉曼信号，重建为二维拉曼图像后，根据其与OCT腔体平铺三维图像中对应Y轴的位置，进行三维重建，生成拉曼三维图像。
在步骤6所述的OCT腔体平铺三维图像与步骤7所述的拉曼三维图像通过一种图像融合算法可实现融合显示：在拉曼三维图像中，计算峰值变化的梯度，其梯度的绝对值作为伪彩的强度值；伪彩设为绿色，将OCT腔体平铺三维图像变为三通道图像，将OCT腔体平铺三维图像的绿色通道替换为伪彩的绿色值，实现OCT腔体平铺三维图像与拉曼三维图像的融合。
在步骤6所述的OCT腔体平铺三维图像与步骤7所述的拉曼三维图像通过另一种图像融合算法可实现融合显示：在需要对腔体中固定组成或特殊成分定位或分析时，可利用与OCT腔体平铺三维图像中对应腔体内壁位置的拉曼信号，利用各位置拉曼信号的特征峰与目标成分的标准图谱的进行模板匹配，通过相似度阈值判断该位置是否含有目标成分，确定含有目标成分的位置后，对OCT腔体平铺三维图像颜色通道进行变换覆盖，实现OCT腔体平铺三维图像与拉曼三维图像的融合。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种内窥OCT-Raman双模成像装置，其特征在于，包括：拉曼成像系统、OCT成像系统、内窥镜系统、图像重建系统和显示系统；所述内窥镜系统分别与拉曼成像系统和OCT成像系统相连，所述拉曼成像系统和OCT成像系统皆与图像重建系统相连，所述图像重建系统与显示系统相连。
2.根据权利要求1所述的内窥OCT-Raman双模成像装置，其特征在于，所述拉曼成像系统包括拉曼激发光光源装置、拉曼光谱仪，所述拉曼激发光光源装置和拉曼光谱仪分别通过光纤与内窥镜系统连接，拉曼光谱仪和图像重建系统连接。
3.根据权利要求1所述的内窥OCT-Raman双模成像装置，其特征在于，所述OCT成像系统包括激光光源装置、OCT干涉成像模块和数据采集模块；所述激光光源和OCT干涉成像模块分别通过光纤与内窥镜系统相连，所述数据采集模块与OCT干涉成像模块相连，数据采集模块与图像重建系统连接。
4.根据权利要求1所述的内窥OCT-Raman双模成像装置，其特征在于，所述内窥镜系统包括运动控制模块和内窥探头；所述运动控制模块包括依次相连的光纤连接器、运动控制装置和成像导管，所述内窥探头包括光纤束、滤光片、反射镜、反射镜座、微型电机、连接管、光纤保护套、气囊导管、外保护套和气囊；
其中，所述光纤束和滤光片由光纤保护套相连，所述光纤束包括拉曼成像激发光纤、拉曼信号收集光纤、OCT成像光源光纤和OCT成像反射光纤；其中拉曼成像激发光纤、OCT成像光源光纤和OCT成像反射光纤呈三角形结构紧密排列组成中心光纤束，拉曼信号收集光纤围绕中心光纤束环形布置；
所述反射镜固定在反射镜座上，反射镜座与连接管连接，连接管与光纤保护套由微型电机相连，所述连接管通过微型电机作用，能够围绕光纤束旋转；所述微型电机由运动控制装置控制，微型电机包括空心定子和转子，空心定子与光纤保护套共轴连接，转子与连接管共轴连接，所述转子转动带动与连接管连接的反射镜转动；
所述气囊导管与气囊相接，气囊导管固定在光纤保护套上，气囊随气囊导管移动。
5.根据权利要求1所述的内窥OCT-Raman双模成像装置，其特征在于，所述图像重建系统包括OCT重建系统和拉曼重建系统；所述OCT重建系统与数据采集模块连接，用于将OCT原始腔体二维图像转换为OCT腔体平铺二维图像，再利用OCT三维重建算法，将一系列内窥探头轴向移动获得的数据，生成OCT腔体平铺三维图像，拉曼重建系统与拉曼光谱仪连接，用于生成与OCT腔体平铺三维图像对应的拉曼三维图像。
6.一种如权利要求1所述的内窥OCT-Raman双模成像方法，其特征在于，包括以下步骤：
步骤1：通过拉曼成像系统输出拉曼激发光源，经过内窥镜系统照射管腔样本；通过OCT成像系统输出OCT成像光源，经过内窥镜系统照射管腔样本；
步骤2：管腔样本受到激发光照射，拉曼信号经内窥镜系统采集输入拉曼成像系统；OCT信号被内窥镜系统反射进入OCT成像系统；
步骤3：拉曼成像系统根据管腔样本的拉曼信号生成拉曼光谱送入图像重建系统；OCT成像系统根据OCT干涉成像模块生成OCT图像数据送入图像重建系统；
步骤4：每采集完一次信号，内窥探头转动重新采集信号，重复至旋转一圈；获得OCT原始腔体二维图像，每旋转扫描一圈，内窥探头轴向移动重新采集信号，重复至完成轴向扫描；
步骤5：OCT成像系统通过极坐标展开算法，使OCT原始腔体二维图像转换为OCT腔体平铺二维图像；
步骤6：利用OCT三维重建算法，将一系列内窥探头轴向移动获得的OCT原始腔体二维图像，生成OCT腔体平铺三维图像；
步骤7：拉曼成像系统根据拉曼三维成像算法生成与OCT腔体平铺三维图像对应的拉曼三维图像。
7.根据权利要求6所述的内窥OCT-Raman双模成像方法，其特征在于，步骤5的所述极坐标展开算法为：利用OCT原始腔体二维图像寻找血管中心，重新建立以血管中心为原点，原X轴为极轴的极坐标系，通过图像坐标变换，将极坐标系的像素点变换到直角坐标系下，生成OCT腔体平铺二维图像。
8.根据权利要求6所述的内窥OCT-Raman双模成像方法，其特征在于，步骤6所述的三维重建算法包括体绘制算法和光线追踪算法。
9.根据权利要求6所述的内窥OCT-Raman双模成像方法，其特征在于，步骤7的所述拉曼三维重建算法为：拉曼成像系统根据管腔样本不同位置处不同拉曼位移的信号强度，获得表征腔体内壁不同物质的拉曼峰值，根据设置的拉曼位移，获取与OCT腔体平铺三维图像中对应腔体内壁位置的拉曼强度值，将该拉曼强度值作为纵坐标，XY轴与OCT腔体平铺三维图像对应，构成拉曼三维图像。
10.根据权利要求6所述的内窥OCT-Raman双模成像方法，其特征在于，步骤6所述的OCT腔体平铺三维图像与步骤7所述的拉曼三维图像通过图像融合算法能够实现融合显示。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011075960.6A CN112089404B (zh) | 2020-10-10 | 2020-10-10 | 一种内窥OCT-Raman双模成像装置及成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011075960.6A CN112089404B (zh) | 2020-10-10 | 2020-10-10 | 一种内窥OCT-Raman双模成像装置及成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112089404A true CN112089404A (zh) | 2020-12-18 |
| CN112089404B CN112089404B (zh) | 2022-04-05 |
Family
ID=73782628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011075960.6A Active CN112089404B (zh) | 2020-10-10 | 2020-10-10 | 一种内窥OCT-Raman双模成像装置及成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112089404B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113160040A (zh) * | 2021-05-07 | 2021-07-23 | 上海大学 | 用于医疗内窥镜成像的管腔图像展开方法 |
| CN113812928A (zh) * | 2021-11-22 | 2021-12-21 | 北京航空航天大学 | 基于拉曼光谱和光学相干断层成像的多模态成像装置 |
| CN113812929A (zh) * | 2021-11-22 | 2021-12-21 | 北京航空航天大学 | 一种多模态成像装置 |
| CN116990236A (zh) * | 2023-09-25 | 2023-11-03 | 中国科学院空天信息创新研究院 | 采样控制方法、装置及系统 |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102056534A (zh) * | 2008-04-03 | 2011-05-11 | 英弗拉雷德克斯公司 | 用于进行血管内结构分析以补偿化学分析模态的系统和方法 |
| US20110282166A1 (en) * | 2009-12-18 | 2011-11-17 | The Regents Of The University Of California | System and Method for Efficient Coherence Anti-Stokes Raman Scattering Endoscopic and Intravascular Imaging and Multimodal Imaging |
| CN102256530A (zh) * | 2008-12-18 | 2011-11-23 | 富士胶片株式会社 | 光学探头和光学观察装置 |
| CN102525382A (zh) * | 2012-02-08 | 2012-07-04 | 中国科学院光电技术研究所 | 螺旋扫描共路干涉型内窥扫频oct实时成像方法及系统 |
| CN104541153A (zh) * | 2012-07-02 | 2015-04-22 | 新加坡国立大学 | 涉及利用光纤拉曼光谱的内窥镜进行实时癌症诊断的方法 |
| CN106793982A (zh) * | 2014-09-15 | 2017-05-31 | 圣纳普医疗(巴巴多斯)公司 | 使用组合模态光学探头的系统和方法 |
| CN107518862A (zh) * | 2016-06-20 | 2017-12-29 | 魏晋 | 集成光学相干与非相干成像系统的一种光纤及成像方法 |
| CN107576645A (zh) * | 2017-08-14 | 2018-01-12 | 中国农业大学 | 一种拉曼光谱提取系统及方法 |
| CN107941782A (zh) * | 2017-12-11 | 2018-04-20 | 南京航空航天大学 | 可内窥光纤拉曼探针及检测装置 |
| CN108040243A (zh) * | 2017-12-04 | 2018-05-15 | 南京航空航天大学 | 多光谱立体视觉内窥镜装置及图像融合方法 |
| CN108618758A (zh) * | 2018-04-27 | 2018-10-09 | 华南师范大学 | 血管内光声-光学相干断层成像-近红外光多模态成像装置与方法 |
| CN109069156A (zh) * | 2016-03-16 | 2018-12-21 | Gi视窗公司 | 用于为磁吻合设备提供精确放置的靶向系统 |
| CN109085119A (zh) * | 2018-09-04 | 2018-12-25 | 席鹏 | 一种拉曼层析光谱检测的共聚焦三维成像系统和实现方法 |
| CN109387501A (zh) * | 2018-11-29 | 2019-02-26 | 中国政法大学 | 利用快速三维拉曼光谱重构技术鉴定交叉笔画顺序的方法 |
| CN208973833U (zh) * | 2018-05-02 | 2019-06-14 | 南开大学 | 基于自发荧光技术和光学相干层析术的甲状腺颈部组织分类和识别系统 |
| CN109982657A (zh) * | 2016-11-11 | 2019-07-05 | 直观外科手术操作公司 | 具有多模态图像显示的外科手术系统 |
| US10365219B1 (en) * | 2018-02-06 | 2019-07-30 | Xidian University | Two-mode Raman optical projection tomography system |
| CN110140020A (zh) * | 2016-12-31 | 2019-08-16 | 诺基亚美国公司 | 混合式拉曼光谱与光学相干性断层扫描成像 |
| CN110367941A (zh) * | 2019-08-20 | 2019-10-25 | 东北大学秦皇岛分校 | 一种探测光融合非接触光声-光学相干层析双模成像系统 |
-
2020
- 2020-10-10 CN CN202011075960.6A patent/CN112089404B/zh active Active
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102056534A (zh) * | 2008-04-03 | 2011-05-11 | 英弗拉雷德克斯公司 | 用于进行血管内结构分析以补偿化学分析模态的系统和方法 |
| CN102256530A (zh) * | 2008-12-18 | 2011-11-23 | 富士胶片株式会社 | 光学探头和光学观察装置 |
| US20110282166A1 (en) * | 2009-12-18 | 2011-11-17 | The Regents Of The University Of California | System and Method for Efficient Coherence Anti-Stokes Raman Scattering Endoscopic and Intravascular Imaging and Multimodal Imaging |
| CN102525382A (zh) * | 2012-02-08 | 2012-07-04 | 中国科学院光电技术研究所 | 螺旋扫描共路干涉型内窥扫频oct实时成像方法及系统 |
| CN104541153A (zh) * | 2012-07-02 | 2015-04-22 | 新加坡国立大学 | 涉及利用光纤拉曼光谱的内窥镜进行实时癌症诊断的方法 |
| CN106793982A (zh) * | 2014-09-15 | 2017-05-31 | 圣纳普医疗(巴巴多斯)公司 | 使用组合模态光学探头的系统和方法 |
| CN109069156A (zh) * | 2016-03-16 | 2018-12-21 | Gi视窗公司 | 用于为磁吻合设备提供精确放置的靶向系统 |
| CN107518862A (zh) * | 2016-06-20 | 2017-12-29 | 魏晋 | 集成光学相干与非相干成像系统的一种光纤及成像方法 |
| CN109982657A (zh) * | 2016-11-11 | 2019-07-05 | 直观外科手术操作公司 | 具有多模态图像显示的外科手术系统 |
| CN110140020A (zh) * | 2016-12-31 | 2019-08-16 | 诺基亚美国公司 | 混合式拉曼光谱与光学相干性断层扫描成像 |
| CN107576645A (zh) * | 2017-08-14 | 2018-01-12 | 中国农业大学 | 一种拉曼光谱提取系统及方法 |
| CN108040243A (zh) * | 2017-12-04 | 2018-05-15 | 南京航空航天大学 | 多光谱立体视觉内窥镜装置及图像融合方法 |
| CN107941782A (zh) * | 2017-12-11 | 2018-04-20 | 南京航空航天大学 | 可内窥光纤拉曼探针及检测装置 |
| US10365219B1 (en) * | 2018-02-06 | 2019-07-30 | Xidian University | Two-mode Raman optical projection tomography system |
| CN108618758A (zh) * | 2018-04-27 | 2018-10-09 | 华南师范大学 | 血管内光声-光学相干断层成像-近红外光多模态成像装置与方法 |
| CN208973833U (zh) * | 2018-05-02 | 2019-06-14 | 南开大学 | 基于自发荧光技术和光学相干层析术的甲状腺颈部组织分类和识别系统 |
| CN109085119A (zh) * | 2018-09-04 | 2018-12-25 | 席鹏 | 一种拉曼层析光谱检测的共聚焦三维成像系统和实现方法 |
| CN109387501A (zh) * | 2018-11-29 | 2019-02-26 | 中国政法大学 | 利用快速三维拉曼光谱重构技术鉴定交叉笔画顺序的方法 |
| CN110367941A (zh) * | 2019-08-20 | 2019-10-25 | 东北大学秦皇岛分校 | 一种探测光融合非接触光声-光学相干层析双模成像系统 |
Non-Patent Citations (2)
| Title |
|---|
| CHETAN A.PATIL: "Combined Raman spectroscopy and optical coherence tomography device for tissue characterization", 《OPTICS LETTERS》 * |
| 王磊: "新型膀胱肿瘤光学诊断技术研究进展", 《现代泌尿外科杂志》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113160040A (zh) * | 2021-05-07 | 2021-07-23 | 上海大学 | 用于医疗内窥镜成像的管腔图像展开方法 |
| CN113812928A (zh) * | 2021-11-22 | 2021-12-21 | 北京航空航天大学 | 基于拉曼光谱和光学相干断层成像的多模态成像装置 |
| CN113812929A (zh) * | 2021-11-22 | 2021-12-21 | 北京航空航天大学 | 一种多模态成像装置 |
| CN113812929B (zh) * | 2021-11-22 | 2022-05-13 | 北京航空航天大学 | 一种多模态成像装置 |
| WO2023088110A1 (zh) * | 2021-11-22 | 2023-05-25 | 北京航空航天大学 | 一种多模态成像装置 |
| WO2023088100A1 (zh) * | 2021-11-22 | 2023-05-25 | 北京航空航天大学 | 基于拉曼光谱和光学相干断层成像的多模态成像装置 |
| CN116990236A (zh) * | 2023-09-25 | 2023-11-03 | 中国科学院空天信息创新研究院 | 采样控制方法、装置及系统 |
| CN116990236B (zh) * | 2023-09-25 | 2023-12-01 | 中国科学院空天信息创新研究院 | 采样控制方法、装置及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112089404B (zh) | 2022-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112089404B (zh) | 一种内窥OCT-Raman双模成像装置及成像方法 | |
| US6552796B2 (en) | Apparatus and method for selective data collection and signal to noise ratio enhancement using optical coherence tomography | |
| EP2579777B1 (en) | Second generation hand-held optical imager | |
| JP4619803B2 (ja) | 蛍光断層画像取得装置 | |
| US9332942B2 (en) | Systems, processes and computer-accessible medium for providing hybrid flourescence and optical coherence tomography imaging | |
| US8712504B2 (en) | Hand-held optical probe based imaging system with 3D tracking facilities | |
| JP2002200037A (ja) | 内視鏡装置 | |
| CN102499645A (zh) | 光声与荧光双模一体断层成像系统及成像方法 | |
| CN114145719B (zh) | 双模冠脉血管图像三维融合的方法和融合系统 | |
| JP2010179042A (ja) | 光構造観察装置及びその構造情報処理方法、光構造観察装置を備えた内視鏡システム | |
| CN113520272B (zh) | 一种内窥导管-多模态光学成像耦合检测系统 | |
| CN109745074B (zh) | 一种三维超声成像的系统及方法 | |
| WO2012002302A1 (ja) | 光断層画像化装置及び光断層画像化方法 | |
| CN209847228U (zh) | 一种三维超声成像装置 | |
| JP2012013432A (ja) | 断層画像処理装置及び方法、並びに光干渉断層画像診断装置 | |
| WO2011158849A1 (ja) | 断層画像処理装置及び方法、並びに光干渉断層画像診断装置 | |
| JP5748281B2 (ja) | 光干渉断層画像処理方法及びその装置 | |
| JP2013066559A (ja) | 光干渉断層画像処理方法及びその装置 | |
| JP2012010776A (ja) | 断層画像処理装置及び方法、並びに光干渉断層画像診断装置 | |
| CN108986084B (zh) | 一种新型oct图像显示方法 | |
| JP5752538B2 (ja) | 光干渉断層画像処理方法及びその装置 | |
| Jiang et al. | Hand-held 3D photoacoustic imager with GPS | |
| JP5812785B2 (ja) | 光断層画像処理装置及び光断層画像処理装置の作動方法 | |
| CN114468998B (zh) | 单视角反射式近红外二区荧光动态断层成像系统及方法 | |
| EP2201889A1 (en) | Optical apparatus for acquiring structure information and its processing method of optical interference signal |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230517 Address after: Room 401, Building 2, Shengtai Technology Building, Jiangning District, Nanjing City, Jiangsu Province, 211106 Patentee after: Nanjing Yingyan Medical Technology Co.,Ltd. Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Patentee before: Nanjing University of Aeronautics and Astronautics |
|
| TR01 | Transfer of patent right |