CN113624720B - 基于傅里叶域光学相干层析成像的色散补偿方法 - Google Patents
基于傅里叶域光学相干层析成像的色散补偿方法 Download PDFInfo
- Publication number
- CN113624720B CN113624720B CN202110854092.XA CN202110854092A CN113624720B CN 113624720 B CN113624720 B CN 113624720B CN 202110854092 A CN202110854092 A CN 202110854092A CN 113624720 B CN113624720 B CN 113624720B
- Authority
- CN
- China
- Prior art keywords
- dispersion
- phase data
- interference signal
- depth
- peak
- 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.)
- Active
Links
Images
Classifications
-
- 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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于傅里叶域光学相干层析成像的色散补偿方法,包括:获取FD‑OCT系统对样品成像的A‑scan的干涉信号;对干涉信号进行预处理,以从干涉信号中获得仅含有样品光参考光相干信号的预处理后的干涉信号;提供仅含有二阶色散系数项和三阶色散系数项的色散相位数据组;根据预处理后的干涉信号和色散相位数据组,获得深度空间数据组;通过对深度空间数据组的寻峰操作,获得色散补偿相位数据;根据预处理后的干涉信号和色散补偿相位数据,获得色散补偿后的深度图像数据。本发明实现了对FD‑OCT系统成像的准确快捷的色散补偿,在确保数据处理实时性、满足实时成像要求的同时保证了FD‑OCT系统的高分辨率成像。
Description
技术领域
本发明涉及光学信号处理技术领域,特别涉及一种基于傅里叶域光学相干层析成像的色散补偿方法。
背景技术
光学相干层析成像(Optical Coherence Tomography,OCT)是上世纪末1991年由David Huang首次提出并发表在《自然》杂志上,是一种生物医学成像技术,具有无创、非侵入、高分辨率、实时成像等特点。自该技术提出至今的三十年内,OCT已成功应用在肿瘤科、牙科、心脏科、皮肤科等医学领域,最成功的是在眼科的应用,由于晶状体对光的透过率较高,使得OCT可以对眼底进行高分辨率成像,这有助于医生对眼底的病变进行无损检测并提出治疗方案,OCT技术已成为眼科诊断的“金标准”。
OCT技术分为时域OCT(TD-OCT)和傅里叶域OCT(FD-OCT),傅里叶域OCT较时域相比具有成像速度快,灵敏度高,信噪比低等优势,为OCT对生物样本的实时在体成像提供了可能。傅里叶域OCT按照系统组成可分为谱域OCT(SD-OCT)和扫频OCT(SS-OCT),二者的差别是分光方法不同,SD-OCT采用宽带光源和光栅进行空间分光,而SS-OCT采用波长随时间变化的扫频光源进行时间分光。SD-OCT和SS-OCT之间的系统结构不同,在各项成像性能指标上没有差别。
OCT的探测原理是部分相干光的干涉,系统的轴向分辨率是评价系统性能的重要指标,与单波长激光相比,宽带光源的相干长度大大降低,使得OCT具有了轴向上微米级的高分辨率。但是宽带光源的使用为系统带来了严重的色散失配,导致系统的点扩散函数(PSF)被展宽,从而降低轴向分辨率。具体来说,存在三个方面的因素导致色散失配:
首先,为了系统的灵活性,OCT系统多采用光纤结构,宽带光在光纤中的传播会为系统引入色散失配;其次,样品臂中为了消除色差,往往会采用由多个透镜组合而成的透镜组作为聚焦透镜,宽带光经过透镜组会为系统带来较大的色散失配;最后,由于系统探测的是样本内部的背向散射信号,光在样品中的传播同样会为系统引入色散。
色散的存在使系统的分辨率降低,所以色散补偿技术一直以来都是OCT的研究热点。色散补偿的方法分为物理(硬件)补偿和数值(软件)补偿。
物理补偿主要是针对样品臂中的消色差聚焦透镜组,在选购透镜组时,可以一同购买配套的色散补偿器,该补偿器的材料和厚度经过选择和设计,和透镜组的二阶、三阶色散系数相匹配,进而消除了样品臂和参考臂的色散失配。此种方法不需要额外的计算,方便快捷,但仅考虑了透镜的色散,光纤和样品带来的色散则未被考虑在内。此外也增加了系统成本和复杂性。
数值补偿是指获得了系统探测到的干涉信号后,对信号进行数据处理,利用算法求出系统的二阶、三阶色散系数。常用的方法有迭代补偿法、解卷积法、全深度补偿法等,这些算法性能上各有优势,但计算量较大,降低了系统的成像速度。
综上所述,色散失配会降低系统的分辨率,如何准确快捷地对OCT系统进行色散补偿一直是OCT的研究热点。
发明内容
有鉴于此,本发明提供一种基于傅里叶域光学相干层析成像的色散补偿方法,以实现准确快捷的色散补偿,在确保数据处理实时性、满足实时成像要求的同时保证FD-OCT系统的高分辨率成像。
本发明的技术方案是这样实现的:
一种基于傅里叶域光学相干层析成像的色散补偿方法,包括:
获取傅里叶域光学相干层析成像FD-OCT系统对样品成像的A-scan的干涉信号;
对所述干涉信号进行预处理,以获得预处理后的干涉信号,其中,所述预处理后的干涉信号中仅含有样品光参考光相干信号;
提供仅含有二阶色散系数项和三阶色散系数项的色散相位数据组;
将所述预处理后的干涉信号与所述色散相位数据组相乘后进行快速逆傅里叶变换IFFT,获得深度空间数据组;
通过对所述深度空间数据组的寻峰操作,获得色散补偿相位数据;
将所述色散补偿相位数据与所述预处理后的干涉信号相乘后进行IFFT,获得色散补偿后的深度图像数据。
进一步,所述的对所述干涉信号进行预处理,包括:
去除所述干涉信号中的参考光自相干信号和样品光自相干信号。
进一步,所述色散相位数据组中包括多个色散相位数据;
各个所述色散相位数据之间的二阶色散系数和/或三阶色散系数各不相同。
进一步,所述色散相位数据组中,各个所述二阶色散系数的取值均匀分布于其所属的数量级内所设定的第一阈值范围中,各个所述三阶色散系数的取值均匀分布于其所属的数量级内所设定的第二阈值范围中。
进一步,所述的将所述预处理后的干涉信号与所述色散相位数据组相乘后进行IFFT,获得深度空间数据组,包括:
将所述色散相位数据组中的每一个色散相位数据作为相位项分别与所述预处理后的干涉信号相乘并进行IFFT,以获得分别对应于每一个色散相位数据的每一个深度空间数据,其中,由所获得的所有深度空间数据共同组成所述深度空间数据组。
进一步,所述的通过对所述深度空间数据组的寻峰操作,获得色散补偿相位数据,包括:
针对所述深度空间数据组中的每一个深度空间数据:通过寻峰操作获得该深度空间数据中的各个峰的强度值和所述各个峰所在的深度位置;针对每个峰:获得该峰的强度值相对于该峰所在深度位置的前一个深度位置的强度值的变化率;将强度值的变化率大于第一设定阈值的峰的强度值进行记录;将所记录的所有强度值的总和作为该深度空间数据的总峰值强度;
选出所述深度空间数据组中的所有深度空间数据的总峰值强度最大的深度空间数据,并将获得该深度空间数据所使用的色散相位数据作为所述色散补偿相位数据。
一种非易失性计算机可读存储介质,所述非易失性计算机可读存储介质存储指令,所述指令在由处理器执行时使得所述处理器执行如上任一项所述的基于傅里叶域光学相干层析成像的色散补偿方法中的步骤。
一种电子设备,包括:
至少一个处理器;以及,
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行如上任一项所述的基于傅里叶域光学相干层析成像的色散补偿方法中的步骤。
从上述方案可以看出,采用本发明的基于傅里叶域光学相干层析成像的色散补偿方法,不需要引入额外的硬件设备,可基于已有的FD-COT系统,节约了系统成本;同时在进行色散补偿过程中仅需做多次快速逆傅里叶变换,算法的复杂度小,数据处理过程简单,耗时短,能够做到实时补偿;另外,在实际操作中对于同一种样品仅需进行一次确定色散系数的过程,即可用求得的色散系数对同一种样品进行精准的色散补偿,换言之,对样品的A-scan进行色散补偿后获得的色散系数可以直接用于B-scan(横向扫描)。采用本发明的基于傅里叶域光学相干层析成像的色散补偿方法能够快速对FD-OCT系系统进行色散补偿,提高FD-OCT系统的轴向分辨率。
附图说明
图1为本发明实施例的基于傅里叶域光学相干层析成像的色散补偿方法流程示意图;
图2为本发明实施例中所使用的FD-OCT系统光路结构示意图;
图3为本发明中一个具体实施的流程示意图;
图4为未进行色散补偿前的深度图像示意图;
图5为经过本发明实施例的色散补偿后的深度图像示意图;
图6为本发明实施例提供的电子设备结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明作进一步详细说明。
如图1所示,本发明实施例的基于傅里叶域光学相干层析成像的色散补偿方法,主要包括以下步骤:
步骤1、获取FD-OCT系统对样品成像的A-scan的干涉信号;
步骤2、对干涉信号进行预处理,以获得预处理后的干涉信号,其中,预处理后的干涉信号中仅含有样品光参考光相干信号;
步骤3、提供仅含有二阶色散系数项和三阶色散系数项的色散相位数据组;
步骤4、将预处理后的干涉信号与色散相位数据组相乘后进行IFFT(快速逆傅里叶变换),获得深度空间数据组;
步骤5、通过对深度空间数据组的寻峰操作,获得色散补偿相位数据;
步骤6、将色散补偿相位数据与预处理后的干涉信号相乘后进行IFFT,获得色散补偿后的深度图像数据。
其中,A-scan为轴向扫描。
在可选实施例中,干涉信号为波数空间的干涉信号。其中,若从FD-OCT系统直接获取的干涉信号若为波长空间的干涉信号,则需要将波长空间的干涉信号转换为波数空间的干涉信号。
在可选实施例中,步骤2的对干涉信号进行预处理,包括:
去除干涉信号中的参考光自相干信号和样品光自相干信号。
基于此,经过步骤2后所得到的预处理后的干涉信号,实际上是去除干涉信号中的参考光自相干信号和样品光自相干信号后,仅保留样品光参考光相干信号的干涉信号。
其中,在可选实施例中,将波长空间的干涉信号记为I0(λ),预处理后的样品光参考光相干信号记为I(λ),将I(λ)由波长空间转到波数空间后记为I(k)。在实际操作中,从FD-OCT系统直接得到的A-scan的干涉信号为波长空间的干涉信号。在波长空间去除干涉信号中的参考光自相干信号和样品光自相干信号,以获得波长空间的仅含有样品光参考光相干信号的干涉信号,再将其转换到波数空间,即为波数空间的仅含有样品光参考光相干信号的干涉信号,即预处理后波数空间的干涉信号;其中,λ为干涉信号的波长,k为干涉信号的波数,则可以依据波数波长关系公式
k=2π/λ
将波长空间的干涉信号转换为波数空间的干涉信号,在具体执行由波长空间向波数空间转换时,可采用等间距插值的方法,将去除参考光自相干信号和样品光自相干信号后的干涉信号等间距插值到波数空间,以得到波数空间的预处理后的干涉信号。可以利用诸如MATLAB等数学软件实现对干涉信号由波长空间向波数空间的转换。
在可选实施例中,色散相位数据组中包括多个色散相位数据;每个色散相位数据中均仅含有二阶色散系数项和三阶色散系数项,并且各个色散相位数据之间的二阶色散系数和/或三阶色散系数各不相同。
在可选实施例中,色散相位数据组中,各个二阶色散系数的取值均匀分布于其所属的数量级内所设定的第一阈值范围中,各个三阶色散系数的取值均匀分布于其所属的数量级内所设定的第二阈值范围中。
在本发明实施例中,色散相位公式表示为
ω=2kv
其中,k为波数,v为光速。
从本领域已有的公知技术中可以知晓,二阶色散系数a2的数量级为10-27,三阶色散系数a3的数量级为10-43。在可选实施例中,在此数量级范围内选定适当范围a2∈[b,c],a3∈[d,e],以在该范围内等间隔选取多个二阶色散系数的取值和多个三阶色散系数的取值,其中,[b,c]为关于二阶色散系数a2在10-27数量级中的第一阈值范围,[d,e]为关于三阶色散系数a3在10-43数量级中的第二阈值范围。其中,b、c、d、e的取值可以为[0,9]中的整数,例如,b=4、c=6、d=2、e=3。多个二阶色散系数和多个三阶色散系数的取值的数量可以为任意多个,取值数量越多则所得到的色散补偿相位数据越精确,取值数量越多则计算设备的进行计算所消耗的时间越长,从而使得成像实时性降低甚至丧失。所以,多个二阶色散系数和多个三阶色散系数的取值的数量需要同时兼顾尽量精确的图像数据和成像的实时性,因此,在优选实施例中,二阶色散系数的取值数量为100个,三阶色散系数的取值数量为100个。100个二阶色散系数在10-27的数量级中在[4,6]中等间隔选取(各个二阶色散系数的取值均匀分布于其所属的数量级10-27内所设定的第一阈值范围[4,6]中),100个三阶色散系数在10-43的数量级中在[2,3]中等间隔选取(各个三阶色散系数的取值均匀分布于其所属的数量级10-43内所设定的第二阈值范围[2,3]中)。
利用选取的100个二阶色散系数的取值,和100个三阶色散系数的取值,以及上述色散相位公式,则能够得到10000个色散相位数据,这10000个色散相位数据中,各个色散相位数据之间的二阶色散系数a2和/或三阶色散系数a3各不相同,进而这10000个色散相位数据之间没有完全相同的色散相位数据。由这10000个色散相位数据组成色散相位数据组。
在可选实施例中,步骤4的将预处理后的干涉信号与色散相位数据组相乘后进行IFFT,获得深度空间数据组,具体包括:
将色散相位数据组中的每一个色散相位数据作为相位项分别与预处理后的干涉信号相乘并进行IFFT,以获得分别对应于每一个色散相位数据的每一个深度空间数据,其中,由所获得的所有深度空间数据共同组成深度空间数据组。
在可选实施例中,步骤5的通过对深度空间数据组的寻峰操作,获得色散补偿相位数据,具体包括:
针对深度空间数据组中的每一个深度空间数据:通过寻峰操作获得该深度空间数据中的各个峰的强度值和各个峰所在的深度位置;针对每个峰:获得该峰的强度值相对于该峰所在深度位置的前一个深度位置的强度值的变化率;将强度值的变化率大于第一设定阈值的峰的强度值进行记录;将所记录的所有强度值的总和作为该深度空间数据的总峰值强度;
选出深度空间数据组中的所有深度空间数据的总峰值强度最大的深度空间数据,并将获得该深度空间数据所使用的色散相位数据作为色散补偿相位数据。
其中,总峰值强度最大代表着有多个深度位置成像清晰,是考虑全深度色散补偿的最优结果。
其中,将深度空间数据记为I(z),结合上述可选实施例,由于色散相位数据组中包括10000个色散相位数据,进而通过步骤4,利用预处理后的干涉信号I(k)数据分别与色散相位数据组中的每一个色散相位数据相乘并做IFFT,能够得到10000个分别对应于10000个色散相位数据的深度空间数据I(z),由这10000个深度空间数据I(z)组成了深度空间数据组。而这10000个色散相位数据分别对应于各个二阶色散系数a2和三阶色散系数a3的取值,进而这10000个深度空间数据I(z)分别对应了10000种二阶色散系数a2和三阶色散系数a3的取值。其中,z为样品的深度。
在步骤5中,针对这10000个深度空间数据I(z)分别进行寻峰操作,获得每一个深度空间数据I(z)中的所有的尖锐信号峰,并对峰值进行筛选,选出将强度值的变化率Q(表示峰的锐利程度)大于第一设定阈值的深度位置的强度值作为强度峰值进行记录,其中
Q=[a(i)-a(i-1)]/a(i)
其中,a(i)表示深度空间数据I(z)中的第i个位置(对应于样品的深度z)的信号强度值。上述公式表示第i个位置与相邻的前一个位置i-1的强度值变化差相对于第i个位置的强度值的程度。其中,关于寻峰操作中a(i)的具体位置的选取和相邻位置之间的间隔可以根据需要设定。具体实践中,可使用MATLAB中的“findpeaks”函数实现寻峰操作。
在可选实施例中,第一设定阈值可根据需要进行设定,例如第一设定阈值可以设置为0.4,满足Q值大于0.4的峰值作为高质量信号峰的强度值进行记录,并将记录的所有高质量信号峰的强度值进行求和记为H。
在步骤5中,针对这10000个深度空间数据I(z)分别进行寻峰操作得到10000个H,将这10000个H中最大的H所对应的深度空间数据I(z)中所使用的色散相位数据作为色散补偿相位数据,该色散补偿相位数据中的二阶色散系数a2的值和三阶色散系数a3的值分别作为补偿二阶色散系数值和补偿三阶色散系数值。
以下结合一个具体实施例,对本发明实施例的基于傅里叶域光学相干层析成像的色散补偿方法进行进一步说明。
如图2所示为本实施例所使用的FD-OCT系统光路结构示意图。其中,SLD光源1发出的光经过光纤环形器2传播到光纤准直器3。光纤准直器3将光纤光转换成空间准直光入射到分束器4,由于样品的背向散射信号远小于参考臂的反射信号,为了增大样品臂功率,选择的分束器4的分光比为90:10,经过分束器4后,光被分成了样品臂和参考臂两束。样品臂的光经过扫描振镜5后入射到聚焦透镜组6,待测样品7置于聚焦透镜组6的透镜焦点位置,升降台8用于调节样品高度,以便找到光束焦距位置。参考臂的光入射到参考镜9并反射回分束器4。样品臂和参考臂的探测光经光路原路返回至光纤环形器2,再由光纤环形器2进入光谱仪14。光谱仪14主要包括了光纤准直器10、光栅11、聚焦透镜12、线阵CCD 13。本实施例中,所用待测样品7为六层盖玻片,其中,单层盖玻片厚度为170微米,将六层盖玻片叠在一起,置于样品臂的聚焦透镜组6的透镜聚焦位置。调节参考镜9的位置,使样品臂和参考臂的零光程位置位于待测样品7的表面上方,光谱仪14采集到的干涉信号经由MATLAB进行数据处理得到色散补偿后的深度图像数据,具体参考图3所示的流程图,主要包括以下步骤。
步骤a、接收FD-OCT系统对待测样品成像的A-scan的波长空间干涉信号I0(λ),之后执行步骤b。
步骤b、将波长空间干涉信号I0(λ)中的参考光自相干信号和样品光自相干信号去除,并转换到波数空间,得到预处理后的干涉信号I(k),之后执行步骤c。
其中,依据波数波长关系公式
k=2π/λ
对干涉信号进行波长空间向波数空间的转换,可利用MATLAB采用等间距插值的方法实现。
其中,波数空间的A-scan干涉信号的函数表示为:
其中,ERei2kr为参考臂信号,为样品臂信号,ER为参考光振幅,k为波数,2r为参考臂光程长,a(z)为样品深度z处的振幅,z和z'为样品深度,n为样品折射率,2[r+nz]为样品臂光程长,S(k)为光源的光谱密度函数,表示的时间平均。i为虚数单位,R表示参考臂。
公式(1)中,ER 2为参考光自相干信号,为系统直流项,需要去除。去除直流项后的波数空间干涉信号(A-scan干涉信号)的函数表示为:
经过去除参考光自相干信号和样品光自相干信号并保留样品光参考光相干信号、波长空间向波数空间转换而得到的预处理后的干涉信号的函数表示为:
此时,预处理后的干涉信号变换到空域后的深度信号如图4所示,其中尖峰分布杂乱,代表着存在较大的色散适配,导致了轴向分辨率的严重降低。
步骤c、提供仅含有二阶色散系数项和三阶色散系数项的色散相位数据组,之后执行步骤d。
其中,关于色散的详细说明,路径长度为z的色散介质增加的相位为
其中,为零阶色散,为一阶色散,为二阶色散,为三阶色散,零阶色散是常数相位,一阶色散是群速度,表示波包的传播速度,将FD-OCT系统的相干长度更改为lC,Disp=lc/ng(lc为光源的相干长度,ng为群折射率),群折射率ng>1,在OCT系统中,深度分辨率定义为光源的相干长度,即深度分辨率等于光源的相干长度lc,所以系统相干长度lC,Disp缩短,深度分辨率能够得到改善。
其中,二阶色散和三阶色散会使空间信号展宽,降低FD-OCT系统所获得图像的轴向分辨率,所以本发明实施例中,需要对二阶色散系数和三阶色散系数进行补偿。
本实施例中,采用的色散相位数据的色散相位公式为
在本实施例中,色散相位数据组中包括10000个色散相位数据;每个色散相位数据中均仅含有二阶色散系数项和三阶色散系数项,如公式(5),并且各个色散相位数据之间的二阶色散系数a2和/或三阶色散系数a3各不相同。
在可选实施例中,色散相位数据组中,二阶色散系数a2的取值数量为100个,各个二阶色散系数a2的取值均匀分布于其所属的数量级10-27内的[4,6]范围中,三阶色散系数a3的取值数量为100个,各个三阶色散系数a3的取值均匀分布于其所属的数量级10-43内的[2,3]范围中。
步骤d、将色散相位数据组中的每一个色散相位数据作为相位项分别与预处理后的干涉信号相乘并做IFFT,得到分别对应于每一个色散相位数据的每一个深度空间数据,之后执行步骤e。
其中,色散相位数据组中包括10000个色散相位数据,经过步骤d后,得到10000个深度空间数据I(z),每一个深度空间数据I(z)分别对应于一个色散相位数据,也就是说,每一个深度空间数据I(z)分别对应一组二阶色散系数a2的值和三阶色散系数a3的值。
步骤e、对每一个深度空间数据执行寻峰操作,以得到总峰值强度H,之后执行步骤f。
其中,使用MATLAB中的“findpeaks”函数实现本步骤中的寻峰操作。
其中,通过如下公式得到各个峰的强度值的变化率
Q=[a(i)-a(i-1)]/a(i)
其中,Q为强度值的变化率,其反应了峰的锐利程度,a(i)表示深度空间数据I(z)中的第i个位置(对应于样品的深度z)的信号强度值。上述公式表示第i个位置与相邻的前一个位置i-1的强度值变化差相对于第i个位置的强度值的程度。,在本实施例中,a(i)表示寻峰操作得到的峰的强度值,a(i-1)表示a(i)峰所在深度位置的前一个深度位置的强度值。
其中,筛选出Q值大于0.4的峰值作为高质量信号峰的强度值进行记录,并将记录的所有高质量信号峰的强度值进行求和得到总峰值强度H。
步骤f、选出所有深度空间数据中的总峰值强度H最大的深度空间数据,之后执行步骤g。
本实施例中,共有10000个深度空间数据I(z),进而寻峰操作后得到10000个H,将这10000个H,在这10000个深度空间数据I(z)中,选出H最大的深度空间数据I(z)。
步骤g、将所选出的总峰值强度H最大的深度空间数据所对应的色散相位数据作为色散补偿相位数据,将色散补偿相位数据与预处理后的干涉信号相乘并做IFFT,得到色散补偿后的深度图像数据。
对峰值求和的解释说明如下,因为OCT对样品进行深度成像,每个峰可以看成是每一层的信号强度,如若信号强度较大,说明在本深度内的组织没有受到色散影响且被清晰成像,峰值求和取最大值即代表在多个深度内获得了最大值,即多深度清晰成像。所以,本发明实施例中,筛选出总峰值强度H最大时的二阶色散系数和三阶色散系数作为补偿能够得到清晰成像。
本实施例中,所选出的总峰值强度H最大的深度空间数据所对应的色散相位数据中的二阶色散系数和三阶色散系数分别为a2=4.71×10-27、a3=2.10×10-43。将以此结果为色散系数的色散相位数据与预处理后的干涉信号相乘,即为色散补偿后的干涉信号,表示为:
T(k)=I(k)eiφ(ω) (6)
公式(6)中,eiφ(ω)中的a2=4.71×10-27、a3=2.10×10-43。
对公式(6)得到的色散补偿后的干涉信号T(k)做IFFT,即可获得色散补偿后的深度图像数据,如图5所示。对该色散补偿后的深度图像数据进行测量得到信号峰的半高全宽(FWHM)为8μm。
本实施例中,所采用光源中心波长为λ0=840nm,带宽为Δλ=50nm,由深度分辨率公式:
可以得到,FD-OCT系统轴向分辨的理论值为6.2μm。
通过对比不难看出,经过色散补偿后的结果(8μm)接近于理论值(6.2μm),从而采用本实施例的色散补偿方法实现了FD-COT系统的高分辨率成像。
本发明实施例的基于傅里叶域光学相干层析成像的色散补偿方法中,即使采用遍历的方式进行寻峰操作,但执行该方法的过程中的复杂计算仅有快速逆傅里叶变换IFFT,所以整体上并不影响算法的实时性。
本发明实施例还同时提供一种非易失性计算机可读存储介质,该非易失性计算机可读存储介质存储指令,该指令在由处理器执行时使得所述处理器执行如上述说明中的基于傅里叶域光学相干层析成像的色散补偿方法中的各个步骤。
本发明实施例还同时提供一种的电子设备,如图6所示,该电子设备包括:至少一个处理器100以及存储器200。存储器200和至少一个处理器100通信连接,例如存储器200和至少一个处理器100通过总线连接。存储器200存储有可被至少一个处理器100执行的指令,所述指令被至少一个处理器100执行,以使至少一个处理器100执行如上述说明中的基于傅里叶域光学相干层析成像的色散补偿方法中的各个步骤。
采用本发明实施例的基于傅里叶域光学相干层析成像的色散补偿方法,不需要引入额外的硬件设备,可基于已有的FD-COT系统,节约了系统成本;同时在进行色散补偿过程中仅需做多次快速逆傅里叶变换,算法的复杂度小,数据处理过程简单,耗时短,能够做到实时补偿;另外,在实际操作中对于同一种样品仅需进行一次确定色散系数的过程,即可用求得的色散系数对同一种样品进行精准的色散补偿,换言之,对样品的A-scan进行色散补偿后获得的色散系数可以直接用于B-scan(横向扫描)。采用本发明实施例的基于傅里叶域光学相干层析成像的色散补偿方法能够快速对FD-OCT系系统进行色散补偿,提高FD-OCT系统的轴向分辨率。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (4)
1.一种基于傅里叶域光学相干层析成像的色散补偿方法,包括:
获取傅里叶域光学相干层析成像FD-OCT系统对样品成像的A-scan的干涉信号;
对所述干涉信号进行预处理,以获得预处理后的干涉信号,其中,所述预处理后的干涉信号中仅含有样品光参考光相干信号;
提供仅含有二阶色散系数项和三阶色散系数项的色散相位数据组;其中,所述色散相位数据组中包括多个色散相位数据,各个所述色散相位数据之间的二阶色散系数和/或三阶色散系数各不相同;所述色散相位数据组中,各个所述二阶色散系数的取值均匀分布于其所属的数量级内所设定的第一阈值范围中,各个所述三阶色散系数的取值均匀分布于其所属的数量级内所设定的第二阈值范围中;
将所述预处理后的干涉信号与所述色散相位数据组相乘后进行快速逆傅里叶变换IFFT,获得深度空间数据组;
通过对所述深度空间数据组的寻峰操作,获得色散补偿相位数据;
将所述色散补偿相位数据与所述预处理后的干涉信号相乘后进行IFFT,获得色散补偿后的深度图像数据;
其中,所述将所述预处理后的干涉信号与所述色散相位数据组相乘后进行IFFT,获得深度空间数据组,包括:
将所述色散相位数据组中的每一个色散相位数据作为相位项分别与所述预处理后的干涉信号相乘并进行IFFT,以获得分别对应于每一个色散相位数据的每一个深度空间数据,其中,由所获得的所有深度空间数据共同组成所述深度空间数据组;
其中,所述通过对所述深度空间数据组的寻峰操作,获得色散补偿相位数据,包括:
针对所述深度空间数据组中的每一个深度空间数据:通过寻峰操作获得该深度空间数据中的各个峰的强度值和所述各个峰所在的深度位置;针对每个峰:获得该峰的强度值相对于该峰所在深度位置的前一个深度位置的强度值的变化率;将强度值的变化率大于第一设定阈值的峰的强度值进行记录;将所记录的所有强度值的总和作为该深度空间数据的总峰值强度;
选出所述深度空间数据组中的所有深度空间数据的总峰值强度最大的深度空间数据,并将获得该深度空间数据所使用的色散相位数据作为所述色散补偿相位数据。
2.根据权利要求1所述的基于傅里叶域光学相干层析成像的色散补偿方法,其特征在于,所述的对所述干涉信号进行预处理,包括:
去除所述干涉信号中的参考光自相干信号和样品光自相干信号。
3.一种非易失性计算机可读存储介质,所述非易失性计算机可读存储介质存储指令,其特征在于,所述指令在由处理器执行时使得所述处理器执行如权利要求1或2所述的基于傅里叶域光学相干层析成像的色散补偿方法中的步骤。
4.一种电子设备,其特征在于,包括:
至少一个处理器;以及,
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行如权利要求1或2所述的基于傅里叶域光学相干层析成像的色散补偿方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110854092.XA CN113624720B (zh) | 2021-07-28 | 2021-07-28 | 基于傅里叶域光学相干层析成像的色散补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110854092.XA CN113624720B (zh) | 2021-07-28 | 2021-07-28 | 基于傅里叶域光学相干层析成像的色散补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113624720A CN113624720A (zh) | 2021-11-09 |
CN113624720B true CN113624720B (zh) | 2022-07-26 |
Family
ID=78381246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110854092.XA Active CN113624720B (zh) | 2021-07-28 | 2021-07-28 | 基于傅里叶域光学相干层析成像的色散补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113624720B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116091368B (zh) * | 2023-04-10 | 2023-08-22 | 佛山光微科技有限公司 | 一种色散补偿方法、装置、电子设备及存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN200957082Y (zh) * | 2006-07-14 | 2007-10-10 | 浙江大学 | 光学相干层析成像中的色散补偿系统 |
CN104771144B (zh) * | 2015-04-13 | 2017-07-14 | 中国科学院上海光学精密机械研究所 | 频域光学相干层析成像系统的色散补偿方法 |
CN106361279B (zh) * | 2016-08-30 | 2019-12-03 | 中国科学院上海光学精密机械研究所 | 光学相干层析成像系统全探测深度色散补偿方法 |
CN107661089A (zh) * | 2017-09-19 | 2018-02-06 | 北京工业大学 | 一种频域光学相干层析连续色散补偿成像方法和系统 |
CN110623644A (zh) * | 2019-10-09 | 2019-12-31 | 南京沃福曼医疗科技有限公司 | 一种导管偏振敏感光学相干层析成像解调方法用色散补偿方法 |
CN112630188A (zh) * | 2020-02-17 | 2021-04-09 | 北京理工大学 | 基于分数域参数检测的光学相干层析成像色散补偿方法 |
CN111568374A (zh) * | 2020-04-20 | 2020-08-25 | 南京师范大学 | 一种傅里叶域光学相干层析成像的自动数值色散补偿方法 |
CN112683848B (zh) * | 2020-12-21 | 2022-09-02 | 中国科学院上海光学精密机械研究所 | 光学相干层析成像系统色散补偿方法 |
CN112597947B (zh) * | 2020-12-29 | 2022-03-08 | 电子科技大学 | 一种基于傅里叶域光学相干层析成像技术的色散补偿方法 |
-
2021
- 2021-07-28 CN CN202110854092.XA patent/CN113624720B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113624720A (zh) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
De Boer et al. | Twenty-five years of optical coherence tomography: the paradigm shift in sensitivity and speed provided by Fourier domain OCT | |
Zhang et al. | Graphics processing unit accelerated non-uniform fast Fourier transform for ultrahigh-speed, real-time Fourier-domain OCT | |
Yun et al. | High-speed spectral-domain optical coherence tomography at 1.3 µm wavelength | |
JP4389032B2 (ja) | 光コヒーレンストモグラフィーの画像処理装置 | |
US8457440B1 (en) | Method and system for background subtraction in medical optical coherence tomography system | |
US8565499B2 (en) | Methods, systems, and computer readable media for synthetic wavelength-based phase unwrapping in optical coherence tomography and spectral domain phase microscopy | |
JP5371315B2 (ja) | 光干渉断層撮像方法および光干渉断層撮像装置 | |
EP1470410B1 (en) | Apparatus and method for rangings and noise reduction of low coherence interferometry (lci) and optical coherence tomography (oct) signals by parallel detection of spectral bands | |
US7872757B2 (en) | Apparatus and method for ranging and noise reduction of low coherence interferometry LCI and optical coherence tomography OCT signals by parallel detection of spectral bands | |
CN112597947B (zh) | 一种基于傅里叶域光学相干层析成像技术的色散补偿方法 | |
Ralston et al. | Real-time interferometric synthetic aperture microscopy | |
Kuo et al. | Balanced detection for spectral domain optical coherence tomography | |
CN112683848B (zh) | 光学相干层析成像系统色散补偿方法 | |
Yu et al. | High-resolution extended source optical coherence tomography | |
US20130301033A1 (en) | Submicron Resolution Spectral-Domain Optical Coherence Tomography | |
Wang et al. | Speckle reduction in optical coherence tomography using angular compounding by B-scan Doppler-shift encoding | |
JP2005283155A (ja) | 光干渉断層像撮像法における分散補正装置 | |
CN113624720B (zh) | 基于傅里叶域光学相干层析成像的色散补偿方法 | |
CN111289470B (zh) | 基于计算光学的oct测量成像方法 | |
Dhalla et al. | Complete complex conjugate resolved heterodyne swept-source optical coherence tomography using a dispersive optical delay line | |
Meleppat et al. | Optical frequency domain imaging with a rapidly swept laser in the 1300nm bio-imaging window | |
Cao et al. | Non-harmonic analysis applied to optical coherence tomography imaging | |
CN116626888A (zh) | 偏振敏感光学相干层析成像的全深度色散补偿方法、系统 | |
Zhang et al. | Self-spectral calibration for spectral domain optical coherence tomography | |
Ding et al. | A new spectral calibration method for Fourier domain optical coherence tomography |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20211109 Assignee: Weiwu Photon (Beijing) Technology Co.,Ltd. Assignor: TSINGHUA University Contract record no.: X2023980037583 Denomination of invention: Dispersion compensation method based on Fourier domain optical coherence tomography imaging Granted publication date: 20220726 License type: Common License Record date: 20230704 |