CN112057044A - 一种oct图像反射和衰减系数提取方法 - Google Patents
一种oct图像反射和衰减系数提取方法 Download PDFInfo
- 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
Links
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
- 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/0073—Measuring 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex 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采用宽带光源,通过参考臂调整z轴延迟,对待测物的不同深度成像,速度低,在临床实用度较低;频域OCT分为两种:一种是基于宽带光源加上分光光谱仪结构的频域OCT系统;另一种是基于扫频激光器和点探测器结构的频域OCT(OFDI)。由于OFDI的成像速度快,结构简单,在临床得到了广泛的应用和研究。基于OFDI的OCT信号与组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂的吸收等因素相关,干扰因素较多。其中组织的反射和衰减系数是真正有临床价值的信号。这里提出一种反射和衰减系数提取方法。
发明内容
本发明的目的是针对现有OCT图像提取中使用的OFDI方法因人体组织干扰较多,影响处理质量和速度的问题,发明一种OCT图像反射和衰减系数提取方法。
本发明的技术方案是:
一种OCT图像反射和衰减系数提取方法,其特征是它包括:
(一)衰减系数提取:
考虑组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂吸收等因素,OCT信号可以用上式表达;其中I0是镜头出射光功率,T(z)是镜头的点扩散函数,η是造影剂吸收系数,S(z)表征光源的相干长度,μb是组织的反射系数,μt是组织衰减系数;公式1可以写成:
两边取对数:
组织的衰减系数主要和组织的散射和吸收相关,当光入射到组织上时会发生散射,在所有的散射光中只有在镜头的数值孔径(NA)之内的散射光才能被镜头收集,并形成OCT信号;因此,组织的反射系数μb远远小于衰减系数μt;再加上对数变换和组织的吸收,因此,公式(3)的第二项和第三项相比可以忽略不计;由于相同组织的衰减系数相同,因此公式(3)可以简化成:
对公式(4)进行线性拟合即可得到不同组织的衰减系数。
(二)反射系数提取;
忽略组织的吸收并且假设反射系数和散射系数成线性关系,即:
μt=α*μb (5)
公式(1)可以写成:
两边积分可得:
I(∞)=0,则:
通过公式(8)即可得到组织的反射系数;
公式(8)忽略了组织的吸收并且假设衰减和反射系数呈线性关系;考虑到组织的吸收也是和组织特性相关,即和反射或者散射系数相关,并且在实际情况中,反射系数与散射系数不一定成线性关系;可以进一步的将公式(8)扩展成:
通过调整a、b和c得到更加符合实际情况的反射系数;进一步的还可以增加到三次项和更高次项。
本发明的有益效果是:
本发明能明显提高图像处理速度和质量,降低计算难度,提高处理速度。
附图说明
图1是本发明的组织的反射与散射示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示。
一种OCT图像反射和衰减系数提取方法,其特征是它包括:
(一)衰减系数提取:
考虑组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂吸收等因素,OCT信号可以用上式表达;其中I0是镜头出射光功率,T(z)是镜头的点扩散函数,η是造影剂吸收系数,S(z)表征光源的相干长度,μb是组织的反射系数,μt是组织衰减系数;公式1可以写成:
两边取对数:
组织的衰减系数主要和组织的散射和吸收相关。从图1可以看出,当光入射到组织上时会发生散射,在所有的散射光中只有在镜头的数值孔径(NA)之内的散射光才能被镜头收集,并形成OCT信号;因此,组织的反射系数μb远远小于衰减系数μt;再加上对数变换和组织的吸收,因此,公式(3)的第二项和第三项相比可以忽略不计;由于相同组织的衰减系数相同,因此公式(3)可以简化成:
对公式(4)进行线性拟合即可得到不同组织的衰减系数。
(二)反射系数提取;
忽略组织的吸收并且假设反射系数和散射系数成线性关系,即:
μt=α*μb (5)
公式(1)可以写成:
两边积分可得:
I(∞)=0,则:
通过公式(8)即可得到组织的反射系数;
公式(8)忽略了组织的吸收并且假设衰减和反射系数呈线性关系;考虑到组织的吸收也是和组织特性相关,即和反射或者散射系数相关,并且在实际情况中,反射系数与散射系数不一定成线性关系;可以进一步的将公式(8)扩展成:
通过调整a、b和c得到更加符合实际情况的反射系数;进一步的还可以增加到三次项和更高次项。
本发明未涉及部分与现有技术相同或可采用现有技术加以实现。
Claims (3)
1.一种OCT图像反射和衰减系数提取方法,其特征是它包括衰减系数提取方法和反射系数提取。
2.根据权利要求1所述的方法,其特征是所述的衰减系数提取方法是指:
考虑组织的反射、衰减、光源相干长度、镜头的点扩散函数和造影剂吸收等因素,OCT信号可以用上式表达;其中I0是镜头出射光功率,T(z)是镜头的点扩散函数,η是造影剂吸收系数,S(z)表征光源的相干长度,μb是组织的反射系数,μt是组织衰减系数;公式1可以写成:
两边取对数:
组织的衰减系数主要和组织的散射和吸收相关,当光入射到组织上时会发生散射,在所有的散射光中只有在镜头的数值孔径(NA)之内的散射光才能被镜头收集,并形成OCT信号;因此,组织的反射系数μb远远小于衰减系数μt;再加上对数变换和组织的吸收,因此,公式(3)的第二项和第三项相比可以忽略不计;由于相同组织的衰减系数相同,因此公式(3)可以简化成:
对公式(4)进行线性拟合即可得到不同组织的衰减系数。
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)
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测量样品深度分辨衰减系数的方法和系统 |
-
2020
- 2020-08-25 CN CN202010861217.7A patent/CN112057044A/zh active Pending
Patent Citations (3)
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)
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 |