CN112057044A - 一种oct图像反射和衰减系数提取方法 - Google Patents

一种oct图像反射和衰减系数提取方法 Download PDF

Info

Publication number
CN112057044A
CN112057044A CN202010861217.7A CN202010861217A CN112057044A CN 112057044 A CN112057044 A CN 112057044A CN 202010861217 A CN202010861217 A CN 202010861217A CN 112057044 A CN112057044 A CN 112057044A
Authority
CN
China
Prior art keywords
tissue
coefficient
reflection
attenuation
absorption
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
Application number
CN202010861217.7A
Other languages
English (en)
Inventor
徐亚伟
李俊峡
张毅
焦力群
黄进宇
丁明
陶魁园
李世强
张弢
刘子旭
鲁坤
匡皓
陆维
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Forssmann Medical Technology Co ltd
Original Assignee
Nanjing Forssmann Medical Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Forssmann Medical Technology Co ltd filed Critical Nanjing Forssmann Medical Technology Co ltd
Priority to CN202010861217.7A priority Critical patent/CN112057044A/zh
Publication of CN112057044A publication Critical patent/CN112057044A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062Arrangements for scanning
    • A61B5/0066Optical coherence imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0073Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Mathematical Analysis (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Mathematical Optimization (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

一种OCT图像反射和衰减系数提取方法,其特征是它包括:(一)衰减系数提取和(二)反射系数提取;通过对提取公式的简化处理,能很快得到需要的图像。本发明具有方法简单,可操作性强的优点。

Description

一种OCT图像反射和衰减系数提取方法
技术领域
本发明涉及一种OCT技术,尤其是一种OCT图像处理技术,具体地说是一种OCT图像反射和衰减系数提取方法。
背景技术
目前,OCT等技术手段已经广泛的用于各种疾病的诊断和治疗。OCT技术分为时域OCT和频域OCT。时域OCT采用宽带光源,通过参考臂调整z轴延迟,对待测物的不同深度成像,速度低,在临床实用度较低;频域OCT分为两种:一种是基于宽带光源加上分光光谱仪结构的频域OCT系统;另一种是基于扫频激光器和点探测器结构的频域OCT(OFDI)。由于OFDI的成像速度快,结构简单,在临床得到了广泛的应用和研究。基于OFDI的OCT信号与组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂的吸收等因素相关,干扰因素较多。其中组织的反射和衰减系数是真正有临床价值的信号。这里提出一种反射和衰减系数提取方法。
发明内容
本发明的目的是针对现有OCT图像提取中使用的OFDI方法因人体组织干扰较多,影响处理质量和速度的问题,发明一种OCT图像反射和衰减系数提取方法。
本发明的技术方案是:
一种OCT图像反射和衰减系数提取方法,其特征是它包括:
(一)衰减系数提取:
Figure BDA0002648204370000011
考虑组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂吸收等因素,OCT信号可以用上式表达;其中I0是镜头出射光功率,T(z)是镜头的点扩散函数,η是造影剂吸收系数,S(z)表征光源的相干长度,μb是组织的反射系数,μt是组织衰减系数;公式1可以写成:
Figure BDA0002648204370000012
两边取对数:
Figure BDA0002648204370000021
组织的衰减系数主要和组织的散射和吸收相关,当光入射到组织上时会发生散射,在所有的散射光中只有在镜头的数值孔径(NA)之内的散射光才能被镜头收集,并形成OCT信号;因此,组织的反射系数μb远远小于衰减系数μt;再加上对数变换和组织的吸收,因此,公式(3)的第二项和第三项相比可以忽略不计;由于相同组织的衰减系数相同,因此公式(3)可以简化成:
Figure BDA0002648204370000022
对公式(4)进行线性拟合即可得到不同组织的衰减系数。
(二)反射系数提取;
忽略组织的吸收并且假设反射系数和散射系数成线性关系,即:
μt=α*μb (5)
公式(1)可以写成:
Figure BDA0002648204370000023
两边积分可得:
Figure BDA0002648204370000024
I(∞)=0,则:
Figure BDA0002648204370000025
通过公式(8)即可得到组织的反射系数;
公式(8)忽略了组织的吸收并且假设衰减和反射系数呈线性关系;考虑到组织的吸收也是和组织特性相关,即和反射或者散射系数相关,并且在实际情况中,反射系数与散射系数不一定成线性关系;可以进一步的将公式(8)扩展成:
Figure BDA0002648204370000031
通过调整a、b和c得到更加符合实际情况的反射系数;进一步的还可以增加到三次项和更高次项。
本发明的有益效果是:
本发明能明显提高图像处理速度和质量,降低计算难度,提高处理速度。
附图说明
图1是本发明的组织的反射与散射示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示。
一种OCT图像反射和衰减系数提取方法,其特征是它包括:
(一)衰减系数提取:
Figure BDA0002648204370000032
考虑组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂吸收等因素,OCT信号可以用上式表达;其中I0是镜头出射光功率,T(z)是镜头的点扩散函数,η是造影剂吸收系数,S(z)表征光源的相干长度,μb是组织的反射系数,μt是组织衰减系数;公式1可以写成:
Figure BDA0002648204370000033
两边取对数:
Figure BDA0002648204370000041
组织的衰减系数主要和组织的散射和吸收相关。从图1可以看出,当光入射到组织上时会发生散射,在所有的散射光中只有在镜头的数值孔径(NA)之内的散射光才能被镜头收集,并形成OCT信号;因此,组织的反射系数μb远远小于衰减系数μt;再加上对数变换和组织的吸收,因此,公式(3)的第二项和第三项相比可以忽略不计;由于相同组织的衰减系数相同,因此公式(3)可以简化成:
Figure BDA0002648204370000042
对公式(4)进行线性拟合即可得到不同组织的衰减系数。
(二)反射系数提取;
忽略组织的吸收并且假设反射系数和散射系数成线性关系,即:
μt=α*μb (5)
公式(1)可以写成:
Figure BDA0002648204370000043
两边积分可得:
Figure BDA0002648204370000044
I(∞)=0,则:
Figure BDA0002648204370000045
通过公式(8)即可得到组织的反射系数;
公式(8)忽略了组织的吸收并且假设衰减和反射系数呈线性关系;考虑到组织的吸收也是和组织特性相关,即和反射或者散射系数相关,并且在实际情况中,反射系数与散射系数不一定成线性关系;可以进一步的将公式(8)扩展成:
Figure BDA0002648204370000051
通过调整a、b和c得到更加符合实际情况的反射系数;进一步的还可以增加到三次项和更高次项。
本发明未涉及部分与现有技术相同或可采用现有技术加以实现。

Claims (3)

1.一种OCT图像反射和衰减系数提取方法,其特征是它包括衰减系数提取方法和反射系数提取。
2.根据权利要求1所述的方法,其特征是所述的衰减系数提取方法是指:
Figure FDA0002648204360000011
考虑组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂吸收等因素,OCT信号可以用上式表达;其中I0是镜头出射光功率,T(z)是镜头的点扩散函数,η是造影剂吸收系数,S(z)表征光源的相干长度,μb是组织的反射系数,μt是组织衰减系数;公式1可以写成:
Figure FDA0002648204360000012
两边取对数:
Figure FDA0002648204360000013
组织的衰减系数主要和组织的散射和吸收相关,当光入射到组织上时会发生散射,在所有的散射光中只有在镜头的数值孔径(NA)之内的散射光才能被镜头收集,并形成OCT信号;因此,组织的反射系数μb远远小于衰减系数μt;再加上对数变换和组织的吸收,因此,公式(3)的第二项和第三项相比可以忽略不计;由于相同组织的衰减系数相同,因此公式(3)可以简化成:
Figure FDA0002648204360000014
对公式(4)进行线性拟合即可得到不同组织的衰减系数。
3.根据权利要求1所述的方法,其特征是所述的反射系数提取:
忽略组织的吸收并且假设反射系数和散射系数成线性关系,即:
μt=α*μb (5)
公式(1)可以写成:
Figure FDA0002648204360000021
两边积分可得:
Figure FDA0002648204360000022
I(∞)=0,则:
Figure FDA0002648204360000023
通过公式(8)即可得到组织的反射系数;
公式(8)忽略了组织的吸收并且假设衰减和反射系数呈线性关系;考虑到组织的吸收也是和组织特性相关,即和反射或者散射系数相关,并且在实际情况中,反射系数与散射系数不一定成线性关系;可以进一步的将公式(8)扩展成:
Figure FDA0002648204360000024
通过调整a、b和c得到更加符合实际情况的反射系数;进一步的还可以增加到三次项和更高次项。
CN202010861217.7A 2020-08-25 2020-08-25 一种oct图像反射和衰减系数提取方法 Pending CN112057044A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010861217.7A CN112057044A (zh) 2020-08-25 2020-08-25 一种oct图像反射和衰减系数提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010861217.7A CN112057044A (zh) 2020-08-25 2020-08-25 一种oct图像反射和衰减系数提取方法

Publications (1)

Publication Number Publication Date
CN112057044A true CN112057044A (zh) 2020-12-11

Family

ID=73659250

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010861217.7A Pending CN112057044A (zh) 2020-08-25 2020-08-25 一种oct图像反射和衰减系数提取方法

Country Status (1)

Country Link
CN (1) CN112057044A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090221920A1 (en) * 2008-01-18 2009-09-03 Boppart Stephen A Low-coherence interferometry and optical coherence tomography for image-guided surgical treatment of solid tumors
CN102046071A (zh) * 2008-06-02 2011-05-04 光学实验室成像公司 用于从光学相干断层扫描图像获得组织特性的定量方法
CN105996999A (zh) * 2016-05-19 2016-10-12 南京航空航天大学 基于oct测量样品深度分辨衰减系数的方法和系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090221920A1 (en) * 2008-01-18 2009-09-03 Boppart Stephen A Low-coherence interferometry and optical coherence tomography for image-guided surgical treatment of solid tumors
CN102046071A (zh) * 2008-06-02 2011-05-04 光学实验室成像公司 用于从光学相干断层扫描图像获得组织特性的定量方法
CN105996999A (zh) * 2016-05-19 2016-10-12 南京航空航天大学 基于oct测量样品深度分辨衰减系数的方法和系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHENYANG XU 等: "Characterization of atherosclerosis plaques by measuring both backscattering and attenuation coefficients in optical coherence tomography", 《JOURNAL OF BIOMEDICAL OPTICS》 *
K. A. VERMEER 等: "Depth-resolved model-based reconstruction of attenuation coefficients in optical coherence tomography", 《BIOMEDICAL OPTICS EXPRESS》 *
SHUANG CHANG 等: "Review of methods and applications of attenuation coefficient measurements with optical coherence tomography", 《J. BIOMED. OPT》 *

Similar Documents

Publication Publication Date Title
US10709419B2 (en) Dual modality imaging system for coregistered functional and anatomical mapping
AU2012332233B2 (en) Dual modality imaging system for coregistered functional and anatomical mapping
US9757092B2 (en) Method for dual modality optoacoustic imaging
US10433732B2 (en) Optoacoustic imaging system having handheld probe utilizing optically reflective material
WO2017133083A1 (zh) 基于全空间调制谱分割角度复合的微血管造影方法与系统
CN105996999B (zh) 基于oct测量样品深度分辨衰减系数的方法和系统
US11525666B2 (en) Spectral contrast optical coherence tomography angiography
Zhou et al. A noise reduction method for photoacoustic imaging in vivo based on EMD and conditional mutual information
Ng et al. Denoising of human speech using combined acoustic and EM sensor signal processing
Pickwell-MacPherson Practical considerations for in vivo THz imaging
CN112704470A (zh) 分光谱频域相干断层成像系统
Zhou et al. Wavelet de-noising method with adaptive threshold selection for photoacoustic tomography
CN106994006A (zh) 双模态成像系统
CN112057044A (zh) 一种oct图像反射和衰减系数提取方法
KR101424714B1 (ko) 레이저의 후산란 특성을 이용한 생체 조직의 분광학적 형태 추출 방법 및 시스템
CN115644811A (zh) 一种基于光学相干断层成像系统的快速的投影图像重建方法
CN105030202A (zh) 基于偏芯保偏光纤干涉增强的oct成像方法
Xue et al. Application of Wavelet Transform in the Denoising of Optical Weak Coherence Signal in Human Eye Examination
Lu et al. Low-sidelobe limited diffraction optical coherence tomography
CN113017569B (zh) 基于光谱子带时域自相关的皮肤伤口愈合情况检查系统
Hu et al. A hybrid denoising approach for PPG signals utilizing variational mode decomposition and improved wavelet thresholding
CN207590666U (zh) 双模态成像系统
Opilski et al. Optical diagnostic system for visualization subcutaneous blood vessels
Li et al. Enhancement of non-air conducted speech based on wavelet-packet adaptive threshold
Boroujeni et al. Microwave Sensor for Detection of Optical Transparent Foreign Body in Soft Tissue: Eye

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20201211