CN114748032A - 一种基于oct血管成像技术的运动噪声补偿方法 - Google Patents
一种基于oct血管成像技术的运动噪声补偿方法 Download PDFInfo
- Publication number
- CN114748032A CN114748032A CN202210180495.5A CN202210180495A CN114748032A CN 114748032 A CN114748032 A CN 114748032A CN 202210180495 A CN202210180495 A CN 202210180495A CN 114748032 A CN114748032 A CN 114748032A
- Authority
- CN
- China
- Prior art keywords
- image
- phase
- calibration
- scan
- oct
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 19
- 238000001228 spectrum Methods 0.000 claims abstract description 6
- 238000007781 pre-processing Methods 0.000 claims abstract description 3
- FCKYPQBAHLOOJQ-UHFFFAOYSA-N Cyclohexane-1,2-diaminetetraacetic acid Chemical compound OC(=O)CN(CC(O)=O)C1CCCCC1N(CC(O)=O)CC(O)=O FCKYPQBAHLOOJQ-UHFFFAOYSA-N 0.000 claims abstract 2
- 238000006073 displacement reaction Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 15
- 238000010587 phase diagram Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000017531 blood circulation Effects 0.000 abstract description 8
- 210000005252 bulbus oculi Anatomy 0.000 abstract description 5
- 238000002583 angiography Methods 0.000 abstract description 4
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 238000012014 optical coherence tomography Methods 0.000 description 49
- 210000001519 tissue Anatomy 0.000 description 8
- 230000008569 process Effects 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 210000001508 eye Anatomy 0.000 description 3
- 206010012689 Diabetic retinopathy Diseases 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000004424 eye movement Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 210000001210 retinal vessel Anatomy 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 206010003210 Arteriosclerosis Diseases 0.000 description 1
- 208000018152 Cerebral disease Diseases 0.000 description 1
- 206010016654 Fibrosis Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 208000035719 Maculopathy Diseases 0.000 description 1
- 206010061323 Optic neuropathy Diseases 0.000 description 1
- 208000017442 Retinal disease Diseases 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 208000011775 arteriosclerosis disease Diseases 0.000 description 1
- 230000004397 blinking Effects 0.000 description 1
- 210000003161 choroid Anatomy 0.000 description 1
- 230000007882 cirrhosis Effects 0.000 description 1
- 208000019425 cirrhosis of liver Diseases 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 208000030533 eye disease Diseases 0.000 description 1
- 210000004220 fundus oculi Anatomy 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 230000001631 hypertensive effect Effects 0.000 description 1
- 208000017169 kidney disease Diseases 0.000 description 1
- 208000002780 macular degeneration Diseases 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 208000020911 optic nerve disease Diseases 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000009711 regulatory function Effects 0.000 description 1
- 210000001525 retina Anatomy 0.000 description 1
- 210000003583 retinal pigment epithelium Anatomy 0.000 description 1
- 208000004644 retinal vein occlusion Diseases 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
- 210000005166 vasculature Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
- A61B3/1241—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes specially adapted for observation of ocular blood flow, e.g. by fluorescein angiography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Hematology (AREA)
- Vascular Medicine (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明涉及一种基于OCT血管成像技术的运动噪声补偿方法,包括如下步骤:步骤1,OCT扫描***执行B‑M‑B式扫描,CCD相机记录原始光谱数据;步骤2,对步骤1中CCD相机采集到的原始光谱数据,进行预处理,得到复数信号图像;根据图像不同的采集位置,将***获取的三维数据分割形成B‑scan簇,将每一组B‑scan簇内的相邻图像数据进行依次配对,形成配对组;步骤3,完成对整体波动的校准;步骤4,继续对配对组内的相位图像φ(z,x)的相位误差做进一步校准;步骤5,当B‑scan图像全部完成校准后,校准后的B‑scan簇用于改善基于振幅、相位或者复数形式的血流成像算法。本方法能有效减少眼球抖动对OCTA血管造影图像造成的运动噪声,提高眼底血流信号的信噪比和成像灵敏度。
Description
技术领域
本发明属于眼科成像技术领域,涉及一种运动噪声补偿方法,具体涉及一种基于OCT血管成像技术的运动噪声补偿方法。
背景技术
眼底的视网膜血管是人体唯一可见的血管,是了解其他脏器血管情况的重要观察窗口。眼底血管形态结构可能会因为一些血管性疾病发生改变,因此通过眼底图像可以帮助诊断各种眼病以及内科病和脑系科的疾病。例如糖尿病视网膜病变、视神经病变、黄斑病变、高血压动脉硬化、肝硬化、肾病、颅脑疾病等等,都可将视网膜血管作为疾病诊断的一个重要指标。快速、准确并且可靠的眼底血管病理特征提取方法在医学诊断领域具有重要的应用价值,对特殊病人以及眼睛不适患者定期进行眼底检查十分必要。
光学相干层析成像(Optical Coherence Tomography,OCT)是一种基于低相干干涉原理的光学成像技术,能对生物活体样品进行非接触、非侵入的横断面实时成像。这项技术可以实现对眼底组织进行微米级别的高分辨率成像,能清晰显示视网膜不同层的结构。糖尿病视网膜病变和视网膜静脉阻塞等疾病,早期征兆表现为眼底血管的改变。对于表征循环调节功能的眼底血管***,OCT技术只能提供有限的结构细节。光相干断层扫描血管成像(Optical Coherence Tomography Angiography,OCTA)是一种基于OCT的新型无创血流影像检测技术,无需荧光标记或其他外源造影对比剂,便可获得眼底毛细血管形态、灌注以及血流动力学信息,目前开始逐步应用于眼科疾病的研究中,应用前景非常广阔。
然而,与眼底照相和荧光血管造影不同,OCT作为一种高分辨率成像技术,对环境变化更加敏感,对信号数据的获取过程要求精度高,不可避免的存在大量干扰,OCTA血管成像会极大的受到患者运动变化引起的伪影和噪声的影响。因为OCTA血管成像的扫描机制需要对同一位置的眼底进行多幅采样,采样过程中眼球运动造成的采样点的幅间差异通常会在最终的血管图像形成运动噪声。运动伪影和噪声会严重妨碍眼底组织和血管结构的可视化,进而影响临床诊断。造成这些伪影的原因有:眼球晃动或震颤、眨眼、眼组织的内在特征、视网膜色素上皮脱落等视网膜疾病和脉络膜基质层的去相关信号对血管信号和血管内血流信号的干扰等。伪影对运动极其敏感,即使是最小的眼球运动,也会产生十分明显的运动伪影,表现为水平或垂直的暗或亮的条纹、图像错位、拉伸或扭曲。在伪影存在的情况下,成像结果会存在争议。只有对OCTA成像过程中的伪影问题进行针对性的优化,才能得到合格的眼底血管图像结果。
针对OCT***图像采集时遭遇的运动噪声问题,目前主流方法主要从硬件或软件方面入手进行运动校正处理。在硬件方面,通过增加一套与OCT***使用分离的光学追踪***,比如激光扫描成像***(Scanning Laser Ophthalmoscopy,SLO),来执行实时跟踪,这种方法不仅增加了***复杂度,同时让患者承受来自SLO的附加光学辐射。软件方面,通过表面分割从结构图中提取组织特征,来进行几何变换校正横向扭曲,完成对OCT图像的运动校准,但是此方法难以满足微血管造影的精准配准需求。
发明内容
本发明的目的在于克服现有技术的不足,提出一种对OCTA的原始信号进行信号补偿的方法,以此来降低由于眼球运动造成的伪影对血流信号的影响。
本发明解决其技术问题是采取以下技术方案实现的:
一种基于OCT血管成像技术的运动噪声补偿方法,其特征在于:包括如下步骤:
步骤1,当执行OCTA采样时,OCT扫描***执行B-M-B式的扫描方式,除了在相互垂直的快轴和慢轴方向进行正常的B-scan扫描外,还需要对眼底的同一位置处进行N次重复的横向扫描,CCD相机记录原始光谱数据;
步骤2,对步骤1中CCD相机采集到的原始光谱数据,进行初步波数校正和傅里叶变换的预处理,得到复数信号图像,复数信号图像的振幅和相位信息将分别用于后续步骤的图像校准;根据图像不同的采集位置,将***获取的三维数据分割形成B-scan簇,将每一组B-scan簇内的相邻图像数据进行依次配对,形成配对组;
步骤3,根据相邻图像振幅信息的相对变化,对步骤2中相邻配对组的图像进行像素移位来校正像素误差,完成对整体波动的校准;
步骤4,在配对组B-scan图像的整体像素误差得到校准的基础上,继续对配对组内的相位图像φ(z,x)的相位误差做进一步校准;
步骤5,在完成前一配对组的B-scan图像像素和相位校准后,将校准后的B-scan图像作为后一相邻配对组的参考帧来重复上述步骤的校准处理,依次重复完成对整个OCT三维数据组的运动噪声校准;当B-scan图像全部完成校准后,校准后的B-scan簇用于改善基于振幅、相位或者复数形式的血流成像算法。
而且,所述步骤2中得到的复数信号图像被解析成用振幅信息和相位信息表示的复数信号:
其中表示经过波数校正傅里叶变换后得到的复数信号,z表示深度方向坐标,x表示沿振镜快速扫描方向坐标,Ioct(k,x)表示波数校准后的OCT光谱信号,k表示波数坐标,I(z,x)表示振幅图像信号,i表示虚数单位,φ(z,x)表示相位图像信号。
而且,所述步骤3中对步骤2中相邻配对组的图像进行像素移位,来完成整体波动的校准的方法为:
步骤3-1,将配对组内的前一幅图像作为参考帧,第二幅图像作为当前帧;
步骤3-2,对B-scan配对组的两个振幅图像I(z,x)进行平滑处理,然后利用下述公式对当前帧Is(z,x)和参考帧IR(z,x)之间进行行和列的互相干信息的计算:
其中,χ(Δz,Δx)表示两幅图像计算的互相关值,Δz表示z方向上的偏移量,Δx表示x方向上的偏移量;在估计的最大位移范围内,通过上述公式来计算两帧图像在列和行方向上的互相关信息χ,因此上述互相关信息χ是与行或列的相对位移Δz和Δx有关的函数;
步骤3-3,通过比较行和列的互相关信息矩阵,找到分别使得行和列的互相关信息χ最大的相对位移量(Δzoptimal,Δxoptimal),即为图像的整体像素位移;
而且,所述步骤4中对配对组内的相位图像φ(z,x)的相位误差做进一步校准的方法为:
步骤4-1,利用参考帧的振幅图像IR来设置阈值生成掩膜,利用掩膜对相位差分图像进行相乘处理,仅提取保留生物组织区域的相位,来参与对整体相位误差的计算;利用下述公式针对掩膜处理后的相位差分图像的每一列,计算此列的中位数相位,得到整体相位误差:
步骤4-2,利用上式计算得到的整体相位误差,再通过下述公式对当前帧相位图里的对应列进行补偿,并依次完成当前帧相位图所有列的相位补偿;在完成对当前帧的相位补偿后,对补偿后的当前帧相位图进行2π折叠处理,从而得到当前帧位于(-π,π)之间的主值相位图,再结合像素校准图像得到最终的校准复信号图像,
步骤4-3,完成对当前帧的相位校准后,利用当前帧校准后的振幅和相位信息,可以完整还原出当前帧的复数信号图,由此完成B-scan图像配对组的运动噪声补偿。
本发明的优点和积极效果是:
本基于OCT血管成像技术的运动噪声补偿方法,能有效减少眼球抖动对OCTA血管造影图像造成的运动噪声,提高眼底血流信号的信噪比和成像灵敏度,此方法适用于改善基于振幅相位,或者复信号形式的血流成像算法。
附图说明
图1是本发明的流程图;
图2是本发明OCT***执行OCTA算法时的扫描模式;
图3是本发明通过整体相位误差对OCT相位图像进行相位补偿演示图。
具体实施方式
以下结合附图对本发明的实施例做进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
一种对OCT血管成像技术的运动噪声补偿方法,其中眼球抖动对OCT信号带来的运动噪声影响可以分为整体波动和相位波动。整体波动会改变OCT的灰度结构图,导致OCT图像的整体像素错移;相位波动虽然不改变结构图,但是同样会对OCTA图像造成伪影。因此,需要对这两种波动情况分别进行校准。本发明的目的在于克服现有技术的不足,提出一种对OCTA的原始信号进行信号补偿的方法,以此来降低由于眼球运动造成的伪影对血流信号的影响。
对于运动噪声给OCT图像带来的不同影响,针对B-scan图像的振幅和相位信息分别进行补偿。
首先针对振幅信息,通过分析相邻B-scan图像的互相干程度,计算两幅图像在深度方向和扫描快轴方向的整体像素位移,然后分别在两个方向上对B-scan图像进行像素移位来完成整体波动的校准。然后在B-scan图像整体像素误差得到校准的基础上,继续对图像的相位误差做进一步校准。计算相邻B-scan图像的相位误差,然后利用得到的相位误差对B-scan图像进行逐列相位补偿,由此消除相位波动并完成对OCT图像的的全部运动噪声补偿。完成了运动补偿的OCT图像,可以直接用于基于不同OCTA血管成像算法的图像运算处理,从而得到更高质量的OCTA影像。
本发明实施流程图如图1所示:
步骤1,当执行OCTA采样时,OCT扫描***执行B-M-B式的扫描方式(在同一横向位置进行重复性扫描的三维扫描模式),除了在相互垂直的快轴和慢轴方向进行正常的B-scan扫描外,需要对眼底的同一位置处进行N次重复的横向扫描,相机记录原始光谱数据;
步骤2,对相机采集到的原始光谱数据,进行初步波数校正和傅里叶变换的预处理,可以被解析成用振幅信息和相位信息表示的复数信号:
其中表示经过波数校正傅里叶变换后得到的复数信号,z表示深度方向坐标,x表示沿振镜快速扫描方向坐标,Ioct(k,x)表示波数校准后的OCT光谱信号,k表示波数坐标,I(z,x)表示振幅图像信号,i表示虚数单位,φ(z,x)表示相位图像信号。得到的OCT复数信号图像的振幅和相位信息将分别用于后续步骤的图像校准。根据图像不同的采集位置,将***获取的三维数据分割形成B-scan簇,将每一组B-scan簇内的相邻图像数据进行依次配对,形成配对组;
步骤3,对相邻配对组的图像进行整体波动校准。配对组内的前一幅图像作为参考帧,第二幅图像作为当前帧。首先对B-scan配对组的两个灰度图像I(z,x)进行平滑处理,然后利用公式(2)对当前帧Is(z,x)和参考帧IR(z,x)之间进行行和列的互相干信息的计算:
其中,χ(Δz,Δx)表示两幅图像计算的互相关值,Δz表示z方向上的偏移量,Δx表示x方向上的偏移量。在估计的最大位移范围内,通过上述公式来计算两帧图像在列和行方向上的互相关信息χ,上述互相关信息χ是与行或列的相对位移有关的函数。通过比较行和列的互相关信息矩阵,找到分别使得行和列的互相关信息χ最大的相对位移量(Δzoptimal,Δxoptimal),即为图像的整体像素位移。通过最大互相关信息得到相邻灰度图像的整体像素位移量后,对当前帧的复信号图像分别进行两个方向上的位移补偿。通常来说,由于眼球运动导致的整体图像像素错位,在快轴x方向较小,在深度z方向较大,所以必要时根据具体硬件扫描情况,如果对最终的OCTA运动补偿的结果影响不大,可以省去对快轴x方向的互相干计算补偿,节省数据处理速度;
步骤4,在配对组B-scan图像的整体像素误差得到校准的基础上,继续对配对组内的相位图像φ(z,x)的相位误差做进一步校准。由于OCT图像内没有生物组织的像素位置的相位随时间随机变化,与组织运动导致的整体相位误差无关,需要在计算过程中予以排除。利用参考帧的灰度图像IR来设置阈值生成掩膜,利用掩膜对相位差分图像进行相乘处理,仅提取保留生物组织区域的相位,来参与对整体相位误差的计算。如公式(3)所示,针对掩膜处理后的相位差分图像的每一列,计算此列的中位数相位,得到整体相位误差:
其中,表示第j幅图像x列的整体相位误差,φi(z,x)表示第j幅图像x列z行的相位,maskstr表示根据灰度图像生成的掩膜。利用上式计算得到的整体相位误差,通过公式4对当前帧相位图里的对应列进行补偿,并依次完成当前帧相位图所有列的相位补偿,如图3所示。在完成对当前帧的相位补偿后,对补偿后的当前帧相位图进行2π折叠处理,从而得到当前帧位于(-π,π)之间的主值相位图。
完成对当前帧的相位校准后,利用当前帧校准后的振幅和相位信息,可以完整还原出当前帧的复数信号图,由此完成B-scan图像配对组的运动噪声补偿;
步骤5,在完成前一配对组的B-scan图像幅值和相位校准后,将校准后的B-scan图像作为后一相邻配对组的参考帧来重复上述步骤的校准处理,依次重复完成对整个OCT三维数据组的运动噪声校准。当B-scan图像全部完成校准后,校准后的B-scan组可以用于改善基于振幅,相位,或者复数形式的血流成像算法,可以有效提高各类OCTA算法的成像效果。
相比其他运动降噪方法,本发明从软件入手从振幅和相位信息两方面对OCT影像进行综合性的运动噪声抑制消除处理,处理后的OCT信号可以适用于不同OCTA算法的运算处理。本方法在软件方面对运动噪声的处理方式简单直接,不需要额外的硬件***来对***的运动噪声进行追踪,同时对硬件的计算要求也更低,应用范围更加广泛。
尽管为说明目的公开了本发明的实施例和附图,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内,各种替换、变化和修改都是可能的,因此,本发明的范围不局限于实施例和附图所公开的内容。
Claims (4)
1.一种基于OCT血管成像技术的运动噪声补偿方法,其特征在于:包括如下步骤:
步骤1,当执行OCTA采样时,OCT扫描***执行B-M-B式的扫描方式,除了在相互垂直的快轴和慢轴方向进行正常的B-scan扫描外,还需要对眼底的同一位置处进行N次重复的横向扫描,CCD相机记录原始光谱数据;
步骤2,对步骤1中CCD相机采集到的原始光谱数据,进行初步波数校正和傅里叶变换的预处理,得到复数信号图像,复数信号图像的振幅和相位信息将分别用于后续步骤的图像校准;根据图像不同的采集位置,将***获取的三维数据分割形成B-scan簇,将每一组B-scan簇内的相邻图像数据进行依次配对,形成配对组;
步骤3,根据相邻图像振幅信息的相对变化,对步骤2中相邻配对组的图像进行像素移位来校正像素误差,完成对整体波动的校准;
步骤4,在配对组B-scan图像的整体像素误差得到校准的基础上,继续对配对组内的相位图像φ(z,x)的相位误差做进一步校准;
步骤5,在完成前一配对组的B-scan图像像素和相位校准后,将校准后的B-scan图像作为后一相邻配对组的参考帧来重复上述步骤的校准处理,依次重复完成对整个OCT三维数据组的运动噪声校准;当B-scan图像全部完成校准后,校准后的B-scan簇用于改善基于振幅、相位或者复数形式的血流成像算法。
3.根据权利要求1所述的一种基于OCT血管成像技术的运动噪声补偿方法,其特征在于:所述步骤3中对步骤2中相邻配对组的图像进行像素移位,来完成整体波动的校准的方法为:
步骤3-1,将配对组内的前一幅图像作为参考帧,第二幅图像作为当前帧;
步骤3-2,对B-scan配对组的两个振幅图像I(z,x)进行平滑处理,然后利用下述公式对当前帧Is(z,x)和参考帧IR(z,x)之间进行行和列的互相干信息的计算:
其中,χ(Δz,Δx)表示两幅图像计算的互相关值,Δz表示z方向上的偏移量,Δx表示x方向上的偏移量;在估计的最大位移范围内,通过上述公式来计算两帧图像在列和行方向上的互相关信息χ,因此上述互相关信息χ是与行或列的相对位移Δz和Δx有关的函数;
步骤3-3,通过比较行和列的互相关信息矩阵,找到分别使得行和列的互相关信息χ最大的相对位移量(Δzoptimal,Δxoptimal),即为图像的整体像素位移;
4.根据权利要求1所述的一种基于OCT血管成像技术的运动噪声补偿方法,其特征在于:所述步骤4中对配对组内的相位图像φ(z,x)的相位误差做进一步校准的方法为:
步骤4-1,利用参考帧的振幅图像IR来设置阈值生成掩膜,利用掩膜对相位差分图像进行相乘处理,仅提取保留生物组织区域的相位,来参与对整体相位误差的计算;利用下述公式针对掩膜处理后的相位差分图像的每一列,计算此列的中位数相位,得到整体相位误差:
步骤4-2,利用上式计算得到的整体相位误差,再通过下述公式对当前帧相位图里的对应列进行补偿,并依次完成当前帧相位图所有列的相位补偿;在完成对当前帧的相位补偿后,对补偿后的当前帧相位图进行2π折叠处理,从而得到当前帧位于(-π,π)之间的主值相位图,再结合像素校准图像得到最终的校准复信号图像,
步骤4-3,完成对当前帧的相位校准后,利用当前帧校准后的振幅和相位信息,可以完整还原出当前帧的复数信号图,由此完成B-scan图像配对组的运动噪声补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210180495.5A CN114748032A (zh) | 2022-02-25 | 2022-02-25 | 一种基于oct血管成像技术的运动噪声补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210180495.5A CN114748032A (zh) | 2022-02-25 | 2022-02-25 | 一种基于oct血管成像技术的运动噪声补偿方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114748032A true CN114748032A (zh) | 2022-07-15 |
Family
ID=82325348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210180495.5A Pending CN114748032A (zh) | 2022-02-25 | 2022-02-25 | 一种基于oct血管成像技术的运动噪声补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114748032A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115953542A (zh) * | 2023-03-14 | 2023-04-11 | 北京心联光电科技有限公司 | 一种光学相干断层扫描微血管成像方法、装置及设备 |
CN116548911A (zh) * | 2023-05-19 | 2023-08-08 | 北京至真互联网技术有限公司 | 一种基于oct技术的眼部血管成像方法及*** |
CN116725492A (zh) * | 2023-07-11 | 2023-09-12 | 江苏金视传奇科技有限公司 | 一种基于光学相干层析成像的血管成像方法及其*** |
-
2022
- 2022-02-25 CN CN202210180495.5A patent/CN114748032A/zh active Pending
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115953542A (zh) * | 2023-03-14 | 2023-04-11 | 北京心联光电科技有限公司 | 一种光学相干断层扫描微血管成像方法、装置及设备 |
CN115953542B (zh) * | 2023-03-14 | 2023-05-12 | 北京心联光电科技有限公司 | 一种光学相干断层扫描微血管成像方法、装置及设备 |
CN116548911A (zh) * | 2023-05-19 | 2023-08-08 | 北京至真互联网技术有限公司 | 一种基于oct技术的眼部血管成像方法及*** |
CN116548911B (zh) * | 2023-05-19 | 2023-12-08 | 北京至真互联网技术有限公司 | 一种基于oct技术的眼部血管成像方法及*** |
CN116725492A (zh) * | 2023-07-11 | 2023-09-12 | 江苏金视传奇科技有限公司 | 一种基于光学相干层析成像的血管成像方法及其*** |
CN116725492B (zh) * | 2023-07-11 | 2023-12-12 | 江苏金视传奇科技有限公司 | 一种基于光学相干层析成像的血管成像方法及其*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114748032A (zh) | 一种基于oct血管成像技术的运动噪声补偿方法 | |
US9418423B2 (en) | Motion correction and normalization of features in optical coherence tomography | |
CN104958061B (zh) | 双目立体视觉三维成像的眼底oct成像方法及其*** | |
US10463247B2 (en) | Automatic three-dimensional segmentation method for OCT and doppler OCT angiography | |
EP1976424B1 (en) | A method of eye examination by optical coherence tomography | |
CN107862724B (zh) | 一种改进的微血管血流成像方法 | |
Cheng et al. | Robust three-dimensional registration on optical coherence tomography angiography for speckle reduction and visualization | |
Qian et al. | In vivo visualization of eye vasculature using super-resolution ultrasound microvessel imaging | |
CN111436905B (zh) | 一种光学相干层析血管造影方法及装置 | |
Qian et al. | Super-resolution ultrasound localization microscopy for visualization of the ocular blood flow | |
Pan et al. | OCTRexpert: a feature-based 3D registration method for retinal OCT images | |
CN109171670B (zh) | 一种基于逆向主成分分析法的3d血管成像算法 | |
CN105608675A (zh) | 一种眼底组织oct图像运动伪差的矫正方法 | |
Solano et al. | Mapping pulsatile optic nerve head deformation using OCT | |
Kim et al. | Strain by virtual extensometers and video-imaging optical coherence tomography as a repeatable metric for IOP-Induced optic nerve head deformations | |
WO2023178078A1 (en) | Oct speckle velocimetry | |
CN112155511B (zh) | 一种基于深度学习补偿oct采集过程中人眼抖动的方法 | |
CN114202512A (zh) | 动态随访肿瘤局部放疗化疗效果的超声造影定量分析方法、***及计算机可读存储介质 | |
CN111493832A (zh) | 一种基于En face-OCT的内窥成像方法 | |
Viedma et al. | A deep learning multi-capture segmentation modality for retinal OCT imaging | |
CN111951315B (zh) | 皮肤光学相干层析图像配准的方法 | |
Liao et al. | Deep-Learning-based Vasculature Extraction for Single-Scan Optical Coherence Tomography Angiography | |
Ma et al. | Deep‐learning segmentation method for optical coherence tomography angiography in ophthalmology | |
Kim | Dynamic Biomechanical Characteristics of the Optic Nerve Head by OCT Imaging | |
Richer et al. | Noise-free one-cardiac-cycle OCT videos for local assessment of retinal tissue deformation |
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 |