CN110742584A - 一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法 - Google Patents
一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法 Download PDFInfo
- Publication number
- CN110742584A CN110742584A CN201910955656.1A CN201910955656A CN110742584A CN 110742584 A CN110742584 A CN 110742584A CN 201910955656 A CN201910955656 A CN 201910955656A CN 110742584 A CN110742584 A CN 110742584A
- Authority
- CN
- China
- Prior art keywords
- matrix
- mueller
- polarization
- phase delay
- depolarization
- 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
-
- 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/21—Polarisation-affecting properties
- G01N21/23—Bi-refringence
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法,其特征是将导管PS‑OCT成像系统得到的测量琼斯矩阵转换成穆勒矩阵,再对穆勒矩阵进行乘法分解,得到相位延迟矩阵,利用相位延迟矩阵的迹计算出组织双折射的相位延迟,实现组织双折射成像。本发明相比于现有的琼斯矩阵法具有速度快,干扰小,准确率高的优点。
Description
技术领域
本发明涉及导管光学相干层析成像(Optical Coherence Tomography,OCT)覆盖心血管、脑血管等导管成像领域,尤其涉及一种导管偏振敏感OCT(Polarization-sensitive OCT)即PS-OCT图像中解调样品的双折射信息,具体地说是一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法。
背景技术
导管OCT成像技术是目前图像分辨率最高的血管成像方法,尤其是导管PS-OCT成像技术,能够解决动脉粥样硬化斑块的稳定性难以在体、实时、快速判断的医学难题,能够提高动脉粥样硬化类疾病的防治效果。然而现有OCT系统在分辨率上已经达到了可能判断组织斑块性质的水平,但在组织穿透能力、图像清晰度和组织斑块类型判断的准确性上仍有不足,采用PS-OCT技术,改善相关技术性能是OCT系统发展的关键方向,也是解决前述关键科学问题的必由之路。
在导管OCT方面,导管PS-OCT是导管OCT技术的延伸,它提供了一种组织双折射性质的定量测量方法。光的双折射改变了光的偏振态且能够与具有定向结构的蛋白质和生物大分子如胶原蛋白、肌动蛋白等相关。导管PS-OCT增强的双折射现象与大量的厚胶原蛋白纤维或内膜平滑肌细胞的存在密切有关,因此导管PS-OCT成像的高分辨率探测可应用于增强后的斑块稳定性测量。此外,导管PS-OCT系统具有评估斑块胶原蛋白和区分正常内膜、纤维斑块、脂质斑块和钙化斑块等方面的潜质。
为此一种导管偏振敏感光学相干层析成像系统及解调方法很好地解决了上述问题,在系统中光源采用快速扫描光源,系统中采用保偏光纤产生正交偏振态的延迟,通过偏振分束器进行偏振分集采集,在一幅图像中同时呈现正交两种输入偏振态的偏振分集成像。解调方法上通过偏振调平、背景信号消除、光谱整形、双态色散消除、插值傅里叶变换、参考面选取、偏振解算、极坐标转换为笛卡尔坐标等一系列步骤最终实现血管双折射成像。
发明内容
本发明的目的是针对现有的OCT导管光学相干层析成像系统只能获取血管强度信息的问题,发明一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法。
本发明的技术方案是:
一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法,其特征是将导管PS-OCT成像系统得到的测量琼斯矩阵转换成穆勒矩阵,再对穆勒矩阵进行乘法分解,得到相位延迟矩阵,利用相位延迟矩阵的迹计算出组织双折射的相位延迟,实现组织双折射成像。
具体包括以下步骤:
首先,将系统输出的四个电场分量定义为琼斯矩阵形式:
其次,进行琼斯矩阵和穆勒矩阵的转换,转换关系为:
其中表示克罗内克积,U表示变换矩阵:
根据琼斯矩阵和穆勒矩阵的转换关系,将系统输出的测量琼斯矩阵J变换成测量穆勒矩阵Q:
第三,对测量穆勒矩阵QZ的元素做欧拉变换,可得:
任取单个A-Scan,假设该条A-Line上的参考点在样品中深度为z,则在z处的测量所得的参考穆勒矩阵QZref可以表示为:
数据采集卡在偏振分集处测量到电信号,是由H通道和V通道的平衡光电探测器经过光电转换输出的,将采集到的电信号分别经过背景噪声去除、双态数值色散补偿、光谱整形滤波和傅里叶插值变换和得到H通道和V通道在z域上的两组信号H1+H2和V1+V2,经过图像分割得到H1,H2,V1,V2个电场信号,构成四个偏振态;以单个A-Line为例,构建成像系统在样品深度为z位置处得到的测量琼斯矩阵J(z)为:
假设该条A-Line上的参考点在样品深度为zref位置处,其对应的参考琼斯矩阵J(zref)为:
其中,对于导管PS-OCT成像系统而言,在解算生物组织的偏振信息时,参考面选择导管的外表面或样品前表面,利用公式(2-7)和公式(2-8)将测量琼斯矩阵J(z)和参考琼斯矩阵J(zref)转换成测量穆勒矩阵S(z)和参考穆勒矩阵S(zref);
设MST是生物组织或样品的往返传输矩阵,Min和Mout分别表示系统光路的穆勒矩阵,在偏振分集处测量到的样品参考面位置处的参考穆勒矩阵S(zref)表示为:
在偏振分集处测量到的样品z位置处的测量穆勒矩阵S(z)表示为:
将S(zref)和S(z)做乘法运算,消除导管在高速旋转时对组织双折射的相位影响,重新整理公式(2-14)和公式(2-15),得到生物组织或样品在深度为z的测量穆勒矩阵M(zref,z):
为了构建M(zref,z)与MS,T为一对相似矩阵,则需要满足QZrefMout必须为可逆矩阵,但如果QZrefMout包含退偏振效应和二向衰减效应,则相似矩阵的条件不成立;实际上,参考穆勒矩阵QZref为酉矩阵,但在导管PS-OCT成像系统中导管中光纤高速旋转必然带来较强的退偏振效应和二向衰减效应,所以Mout中也必然包含退偏振效应和二向衰减效应,则需要通过矩阵分解的方法来消除系统光路中的退偏振效应和二向衰减效应;
对于SMM偏振解算方法的另一个重要环节是通过矩阵分解消除退偏振效应和二向衰减效应,进而得到测量相位延迟穆勒矩阵,计算生物组织或样品的相位延迟,实现组织双折射成像。
所述的穆勒矩阵的分解方法包括乘法分解和求和分解两种,而乘法分解适合于处理先验序列的偏振单元和系统,虽然生物组织的偏振单元是离散的,是复杂多样的,且明显不具备先验序列的特点,但使用L-C穆勒矩阵分解法得到分解后的穆勒矩阵,其结构基本上可以正确地描述生物组织的基本偏振特性;基于导管PS-OCT成像系统,也可以使用对称分解法分解穆勒矩阵,在已知穆勒已退偏的情况下,对称分解法同样能够准确的表达生物组织的基本偏振特性;本发明采用的是利用L-C穆勒矩阵分解法得到测量穆勒矩阵的相位延迟穆勒矩阵MR。
任意描述传光介质的穆勒矩阵,包含目标传光介质调制入射光的三个偏振效应,包括二向衰减、相位延迟或退偏;任意穆勒矩阵都可以分解成下面的形式:
M=MΔMRMD (2-17)
其中,M为任意的穆勒矩阵,MΔ表示目标传光介质调制入射光的退偏矩阵,它描述了目标传光介质对于入射光的退偏振的能力;MR表示目标传光介质调制入射光的相位延迟矩阵,它描述了入射光在经过偏振器件后,产生的两个本征偏振态在快慢轴中传播速率不同,导致了这两个本征偏振态的时域信号存在时延,转换到频域即为相位延迟;MD表示目标传光介质调制入射光的二向衰减矩阵,它描述了入射光经过传光介质的反射光强随入射光偏振态的变化关系;
相位延迟矩阵MR可以表示为:
其中mR是相位延迟矩阵MR的一个3×3的子矩阵,它可以通过公式(2-19)计算得到:
其中m'和mΔ分别是退偏矩阵MΔ和退偏矩阵与相位延迟矩阵的乘积M'矩阵,M'可以表示为:
在未知退偏矩阵MΔ和相位延迟矩阵MR的情况下,可以通过二向衰减矩阵MD和测量穆勒矩阵M的乘积计算得到,二向衰减矩阵MD可以表示为:
二向衰减矩阵的子矩阵mD可以表示为:
其中||·||表示欧几里得向量范数;以上,可以计算出m'矩阵;
退偏矩阵MΔ可以表示为:
在相位延迟矩阵MR的计算过程中,并不需要求解退偏矩阵MΔ,只需要求解其子矩阵mΔ即可,退偏矩阵的子矩阵mΔ它可以通过公式(2-28)计算得到:
其中λi(i=1,2,3)表示子矩阵m'的特征值。
至此,求得退偏矩阵的子矩阵mΔ,将mΔ和m'代入到公式(2-19),联立公式(2-18)可以求解出相位延迟矩阵MR。在导管PS-OCT成像系统中,将测量穆勒矩阵M(zref,z)进行矩阵乘法分解,得到已消除退偏振效应和二向衰减效应的只包含组织双折射的穆勒矩阵,相位延迟矩阵MR(zref,z)的表达式为:
其中测量相位延迟矩阵MR(zref,z),生物组织或样品传输矩阵和系统传输输出矩阵都是只包含双折射成分的穆勒矩阵,且此时为酉矩阵,满足相似矩阵条件,MR(zref,z)与即为一对相似矩阵。总的相位延迟量与分解得到的相位延迟矩阵MR的关系可以表示为:
本发明的有益效果:
本发明相比于琼斯矩阵法具有速度快,干扰小,准确率高的优点。
附图说明
图1是本发明涉及的导管偏振敏感光学相干层析成像解调方法流程图。
图2是本发明的偏振解算方法流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1-2所示。
一种导管偏振敏感光学相干层析成像解调方法,如图1所示,它包括偏振调平、背景信号消除、光谱整形、双态色散消除、插值傅里叶变换、参考面选取、偏振解算、极坐标转换为笛卡尔坐标等一系列步骤最终实现血管双折射成像。
其中的偏振解算方法如图2所示。PS-OCT常用的偏振信息提取方法是相似琼斯矩阵(SJM)方法,该方法是在台式PS-OCT的基础上,针对导管PS-OCT的特点进行改进,但该方法的实际偏振解算结果并不理想,无法正确解算出生物组织或样品的双折射信息。为此本发明提出了相似穆勒矩阵(SMM)方法,SMM方法的核心思想是将导管PS-OCT成像系统得到的测量琼斯矩阵转换成穆勒矩阵,再对穆勒矩阵进行乘法分解,得到相位延迟矩阵,利用相位延迟矩阵的迹计算出组织双折射的相位延迟,实现组织双折射成像。
推算过程如下:
将系统输出的四个电场分量定义为琼斯矩阵形式:
琼斯矩阵和穆勒矩阵的转换关系为:
根据琼斯矩阵和穆勒矩阵的转换关系,将系统输出的测量琼斯矩阵J变换成测量穆勒矩阵Q:
对测量穆勒矩阵QZ的元素做欧拉变换,可得:
任取单个A-Scan,假设该条A-Line上的参考点在样品中深度为z,则在z处的测量所得的参考穆勒矩阵QZref可以表示为:
数据采集卡在偏振分集处测量到电信号,是由H通道和V通道的平衡光电探测器经过光电转换输出的,将采集到的电信号分别经过背景噪声去除、双态数值色散补偿、光谱整形滤波和傅里叶插值变换和得到H通道和V通道在z域上的两组信号H1+H2和V1+V2,经过图像分割得到H1,H2,V1,V2个电场信号,构成四个偏振态。以单个A-Line为例,构建成像系统在样品深度为z位置处得到的测量琼斯矩阵J(z)为:
假设该条A-Line上的参考点在样品深度为zref位置处,其对应的参考琼斯矩阵J(zref)为:
其中,对于导管PS-OCT成像系统而言,在解算生物组织的偏振信息时,参考面可以选择导管的外表面或样品前表面,利用公式(2-7)和公式(2-8)将测量琼斯矩阵J(z)和参考琼斯矩阵J(zref)转换成测量穆勒矩阵S(z)和参考穆勒矩阵S(zref)。
设MST是生物组织或样品的往返传输矩阵,Min和Mout分别表示系统光路的穆勒矩阵,在偏振分集处测量到的样品参考面位置处的参考穆勒矩阵S(zref)表示为:
在偏振分集处测量到的样品z位置处的测量穆勒矩阵S(z)表示为:
将S(zref)和S(z)做乘法运算,消除导管在高速旋转时对组织双折射的相位影响,可以重新整理公式(2-14)和公式(2-15),得到生物组织或样品在深度为z的测量穆勒矩阵M(zref,z):
为了构建M(zref,z)与MS,T为一对相似矩阵,则需要满足QZrefMout必须为可逆矩阵,但如果QZrefMout包含退偏振效应和二向衰减效应,则相似矩阵的条件不成立。实际上,参考穆勒矩阵QZref为酉矩阵,但在导管PS-OCT成像系统中导管中光纤高速旋转必然带来较强的退偏振效应和二向衰减效应,所以Mout中也必然包含退偏振效应和二向衰减效应,则需要通过矩阵分解的方法来消除系统光路中的退偏振效应和二向衰减效应。
对于SMM偏振解算方法的另一个重要环节是通过矩阵分解消除退偏振效应和二向衰减效应,进而得到测量相位延迟穆勒矩阵,计算生物组织或样品的相位延迟,实现组织双折射成像。
穆勒矩阵的常用分解方法包括乘法分解和求和分解两种,而乘法分解适合于处理先验序列的偏振单元和系统,虽然生物组织的偏振单元是离散的,是复杂多样的,且明显不具备先验序列的特点,但使用L-C穆勒矩阵分解法得到分解后的穆勒矩阵,其结构基本上可以正确地描述生物组织的基本偏振特性。基于导管PS-OCT成像系统,也可以使用对称分解法分解穆勒矩阵,在已知穆勒已退偏的情况下,对称分解法同样能够准确的表达生物组织的基本偏振特性。本方法采用的是利用L-C穆勒矩阵分解法得到测量穆勒矩阵的相位延迟穆勒矩阵MR。
任意描述传光介质的穆勒矩阵,包含目标传光介质调制入射光的三个偏振效应,包括二向衰减、相位延迟或退偏。任意穆勒矩阵都可以分解成下面的形式:
M=MΔMRMD (2-17)
其中,M为任意的穆勒矩阵,MΔ表示目标传光介质调制入射光的退偏矩阵,它描述了目标传光介质对于入射光的退偏振的能力。MR表示目标传光介质调制入射光的相位延迟矩阵,它描述了入射光在经过偏振器件后,产生的两个本征偏振态在快慢轴中传播速率不同,导致了这两个本征偏振态的时域信号存在时延,转换到频域即为相位延迟。MD表示目标传光介质调制入射光的二向衰减矩阵,它描述了入射光经过传光介质的反射光强随入射光偏振态的变化关系。
相位延迟矩阵MR可以表示为:
其中mR是相位延迟矩阵MR的一个3×3的子矩阵,它可以通过公式(2-19)计算得到:
其中m'和mΔ分别是退偏矩阵MΔ和退偏矩阵与相位延迟矩阵的乘积M'矩阵,M'可以表示为:
在未知退偏矩阵MΔ和相位延迟矩阵MR的情况下,可以通过二向衰减矩阵MD和测量穆勒矩阵M的乘积计算得到,二向衰减矩阵MD可以表示为:
二向衰减矩阵的子矩阵mD可以表示为:
其中||·||表示欧几里得向量范数。以上,可以计算出m'矩阵。
退偏矩阵MΔ可以表示为:
在相位延迟矩阵MR的计算过程中,并不需要求解退偏矩阵MΔ,只需要求解其子矩阵mΔ即可,退偏矩阵的子矩阵mΔ它可以通过公式(2-28)计算得到:
其中λi(i=1,2,3)表示子矩阵m'的特征值。
至此,求得退偏矩阵的子矩阵mΔ,将mΔ和m'带入到公式(2-19),联立公式(2-18)可以求解出相位延迟矩阵MR。在导管PS-OCT成像系统中,将测量穆勒矩阵M(zref,z)进行矩阵乘法分解,得到已消除退偏振效应和二向衰减效应的只包含组织双折射的穆勒矩阵,相位延迟矩阵MR(zref,z)的表达式为:
其中测量相位延迟矩阵MR(zref,z),生物组织或样品传输矩阵和系统传输输出矩阵都是只包含双折射成分的穆勒矩阵,且此时为酉矩阵,满足相似矩阵条件,MR(zref,z)与即为一对相似矩阵。总的相位延迟量与分解得到的相位延迟矩阵MR的关系可以表示为:
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
Claims (3)
1.一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法,其特征是将导管PS-OCT成像系统得到的测量琼斯矩阵转换成穆勒矩阵,再对穆勒矩阵进行乘法分解,得到相位延迟矩阵,利用相位延迟矩阵的迹计算出组织双折射的相位延迟,实现组织双折射成像。
2.根据权利要求1所述的方法,其特征是它包括以下步骤:
首先,将系统输出的四个电场分量定义为琼斯矩阵形式:
其次,进行琼斯矩阵和穆勒矩阵的转换,转换关系为:
根据琼斯矩阵和穆勒矩阵的转换关系,将系统输出的测量琼斯矩阵J变换成测量穆勒矩阵Q:
第三,对测量穆勒矩阵QZ的元素做欧拉变换,可得:
任取单个A-Scan,假设该条A-Line上的参考点在样品中深度为z,则在z处的测量所得的参考穆勒矩阵QZref可以表示为:
数据采集卡在偏振分集处测量到电信号,是由H通道和V通道的平衡光电探测器经过光电转换输出的,将采集到的电信号分别经过背景噪声去除、双态数值色散补偿、光谱整形滤波和傅里叶插值变换和得到H通道和V通道在z域上的两组信号H1+H2和V1+V2,经过图像分割得到H1,H2,V1,V2个电场信号,构成四个偏振态;以单个A-Line为例,构建成像系统在样品深度为z位置处得到的测量琼斯矩阵J(z)为:
假设该条A-Line上的参考点在样品深度为zref位置处,其对应的参考琼斯矩阵J(zref)为:
其中,对于导管PS-OCT成像系统而言,在解算生物组织的偏振信息时,参考面选择导管的外表面或样品前表面,利用公式(2-7)和公式(2-8)将测量琼斯矩阵J(z)和参考琼斯矩阵J(zref)转换成测量穆勒矩阵S(z)和参考穆勒矩阵S(zref);
设MST是生物组织或样品的往返传输矩阵,Min和Mout分别表示系统光路的穆勒矩阵,在偏振分集处测量到的样品参考面位置处的参考穆勒矩阵S(zref)表示为:
在偏振分集处测量到的样品z位置处的测量穆勒矩阵S(z)表示为:
将S(zref)和S(z)做乘法运算,消除导管在高速旋转时对组织双折射的相位影响,重新整理公式(2-14)和公式(2-15),得到生物组织或样品在深度为z的测量穆勒矩阵M(zref,z):
为了构建M(zref,z)与MS,T为一对相似矩阵,则需要满足QZrefMout必须为可逆矩阵,但如果QZrefMout包含退偏振效应和二向衰减效应,则相似矩阵的条件不成立;实际上,参考穆勒矩阵QZref为酉矩阵,但在导管PS-OCT成像系统中导管中光纤高速旋转必然带来较强的退偏振效应和二向衰减效应,所以Mout中也必然包含退偏振效应和二向衰减效应,则需要通过矩阵分解的方法来消除系统光路中的退偏振效应和二向衰减效应;
对于SMM偏振解算方法的另一个重要环节是通过矩阵分解消除退偏振效应和二向衰减效应,进而得到测量相位延迟穆勒矩阵,计算生物组织或样品的相位延迟,实现组织双折射成像。
3.根据权利要求2所述的方法,其特征是所述的穆勒矩阵的分解方法是利用L-C穆勒矩阵分解法得到测量穆勒矩阵的相位延迟穆勒矩阵MR。
任意描述传光介质的穆勒矩阵,包含目标传光介质调制入射光的三个偏振效应,包括二向衰减、相位延迟或退偏;任意穆勒矩阵都可以分解成下面的形式:
M=MΔMRMD (2-17)
其中,M为任意的穆勒矩阵,MΔ表示目标传光介质调制入射光的退偏矩阵,它描述了目标传光介质对于入射光的退偏振的能力;MR表示目标传光介质调制入射光的相位延迟矩阵,它描述了入射光在经过偏振器件后,产生的两个本征偏振态在快慢轴中传播速率不同,导致了这两个本征偏振态的时域信号存在时延,转换到频域即为相位延迟;MD表示目标传光介质调制入射光的二向衰减矩阵,它描述了入射光经过传光介质的反射光强随入射光偏振态的变化关系;
相位延迟矩阵MR可以表示为:
其中mR是相位延迟矩阵MR的一个3×3的子矩阵,它可以通过公式(2-19)计算得到:
其中m'和mΔ分别是退偏矩阵MΔ和退偏矩阵与相位延迟矩阵的乘积M'矩阵,M'可以表示为:
在未知退偏矩阵MΔ和相位延迟矩阵MR的情况下,可以通过二向衰减矩阵MD和测量穆勒矩阵M的乘积计算得到,二向衰减矩阵MD可以表示为:
二向衰减矩阵的子矩阵mD可以表示为:
退偏矩阵MΔ可以表示为:
在相位延迟矩阵MR的计算过程中,并不需要求解退偏矩阵MΔ,只需要求解其子矩阵mΔ即可,退偏矩阵的子矩阵mΔ它可以通过公式(2-28)计算得到:
其中λi(i=1,2,3)表示子矩阵m'的特征值。
至此,求得退偏矩阵的子矩阵mΔ,将mΔ和m'代入到公式(2-19),联立公式(2-18)可以求解出相位延迟矩阵MR。在导管PS-OCT成像系统中,将测量穆勒矩阵M(zref,z)进行矩阵乘法分解,得到已消除退偏振效应和二向衰减效应的只包含组织双折射的穆勒矩阵,相位延迟矩阵MR(zref,z)的表达式为:
其中测量相位延迟矩阵MR(zref,z),生物组织或样品传输矩阵和系统传输输出矩阵都是只包含双折射成分的穆勒矩阵,且此时为酉矩阵,满足相似矩阵条件,MR(zref,z)与即为一对相似矩阵。总的相位延迟量与分解得到的相位延迟矩阵MR的关系可以表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910955656.1A CN110742584A (zh) | 2019-10-09 | 2019-10-09 | 一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910955656.1A CN110742584A (zh) | 2019-10-09 | 2019-10-09 | 一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110742584A true CN110742584A (zh) | 2020-02-04 |
Family
ID=69277865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910955656.1A Pending CN110742584A (zh) | 2019-10-09 | 2019-10-09 | 一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110742584A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111965114A (zh) * | 2020-08-15 | 2020-11-20 | 天津大学 | 一种导管偏振敏感光学相干层析成像本地双折射解调方法 |
CN113888478A (zh) * | 2021-09-15 | 2022-01-04 | 天津大学 | 一种对血管内导管偏振敏感相干层析成像优化去消偏方法 |
CN115035210A (zh) * | 2022-08-10 | 2022-09-09 | 天津恒宇医疗科技有限公司 | 基于偏振多参量融合的ps-oct可视度提升方法及系统 |
WO2023236406A1 (zh) * | 2022-06-09 | 2023-12-14 | 温州医科大学 | 背向散射Mueller矩阵的可逆极分解方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101513366A (zh) * | 2009-03-18 | 2009-08-26 | 天津大学 | 基于ps-oct的三维可视化口腔龋齿检测装置及检测方法 |
WO2013001032A1 (en) * | 2011-06-28 | 2013-01-03 | Medlumics S.L. | Integrated delay line for optical coherence tomography |
KR101461235B1 (ko) * | 2008-02-29 | 2014-11-19 | 고려대학교 산학협력단 | 생체조직 검사를 위한 편광감도-광간섭 영상획득용 프로브 및 편광 민감 광 간섭 영상검출용 프로브 및 생체조직 검사를 위한 편광감도-광간섭 영상시스템의 구동방법 |
CN104568951A (zh) * | 2015-01-12 | 2015-04-29 | 南京理工大学 | 全场光学相干层析三维医学成像装置及方法 |
EP3085302A1 (en) * | 2015-04-23 | 2016-10-26 | Tomey Corporation | Optical coherence tomography using polarization information |
US20160320171A1 (en) * | 2015-04-30 | 2016-11-03 | Nidek Co., Ltd. | Optical coherence tomography apparatus |
CN106383087A (zh) * | 2016-10-09 | 2017-02-08 | 福建师范大学 | 一种偏振光学相干层析与偏振多光子显微的双模成像系统 |
CN107064001A (zh) * | 2017-02-28 | 2017-08-18 | 福建师范大学 | 基于光开关的单光谱仪偏振频域光学相干层析成像系统 |
US20180120084A1 (en) * | 2016-10-27 | 2018-05-03 | Hitachi-Lg Data Storage, Inc. | Optical image measurement apparatus and optical image measurement method |
CN109164048A (zh) * | 2018-09-18 | 2019-01-08 | 天津大学 | 一种对导管偏振敏感光学相干层析成像偏振解调方法 |
CN109171659A (zh) * | 2018-09-28 | 2019-01-11 | 南京航空航天大学 | 基于琼斯矩阵的光纤型扫频偏振敏感oct成像方法及系统 |
US10191145B2 (en) * | 2013-06-23 | 2019-01-29 | Eric Swanson | Integrated optical system with optical phased array photonic integrated circuit |
-
2019
- 2019-10-09 CN CN201910955656.1A patent/CN110742584A/zh active Pending
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101461235B1 (ko) * | 2008-02-29 | 2014-11-19 | 고려대학교 산학협력단 | 생체조직 검사를 위한 편광감도-광간섭 영상획득용 프로브 및 편광 민감 광 간섭 영상검출용 프로브 및 생체조직 검사를 위한 편광감도-광간섭 영상시스템의 구동방법 |
CN101513366A (zh) * | 2009-03-18 | 2009-08-26 | 天津大学 | 基于ps-oct的三维可视化口腔龋齿检测装置及检测方法 |
WO2013001032A1 (en) * | 2011-06-28 | 2013-01-03 | Medlumics S.L. | Integrated delay line for optical coherence tomography |
US10191145B2 (en) * | 2013-06-23 | 2019-01-29 | Eric Swanson | Integrated optical system with optical phased array photonic integrated circuit |
CN104568951A (zh) * | 2015-01-12 | 2015-04-29 | 南京理工大学 | 全场光学相干层析三维医学成像装置及方法 |
EP3085302A1 (en) * | 2015-04-23 | 2016-10-26 | Tomey Corporation | Optical coherence tomography using polarization information |
US20160320171A1 (en) * | 2015-04-30 | 2016-11-03 | Nidek Co., Ltd. | Optical coherence tomography apparatus |
CN106383087A (zh) * | 2016-10-09 | 2017-02-08 | 福建师范大学 | 一种偏振光学相干层析与偏振多光子显微的双模成像系统 |
US20180120084A1 (en) * | 2016-10-27 | 2018-05-03 | Hitachi-Lg Data Storage, Inc. | Optical image measurement apparatus and optical image measurement method |
CN107064001A (zh) * | 2017-02-28 | 2017-08-18 | 福建师范大学 | 基于光开关的单光谱仪偏振频域光学相干层析成像系统 |
CN109164048A (zh) * | 2018-09-18 | 2019-01-08 | 天津大学 | 一种对导管偏振敏感光学相干层析成像偏振解调方法 |
CN109171659A (zh) * | 2018-09-28 | 2019-01-11 | 南京航空航天大学 | 基于琼斯矩阵的光纤型扫频偏振敏感oct成像方法及系统 |
Non-Patent Citations (3)
Title |
---|
FIRDOUS, S.,ET AL: "Imaging of biological tissues with optical coherence tomography system using jones-mueller calculus", 《PROCEEDINGS OF THE SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING》 * |
TAO KUIYUAN等: "Catheter-Based Polarization Sensitive Optical Coherence Tomography Using Similar Mueller Matrix Method", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
肖悦娱等: "基于斯托克斯空间的任意正交基矢下琼斯矩阵的测量", 《激光与光电子学进展》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111965114A (zh) * | 2020-08-15 | 2020-11-20 | 天津大学 | 一种导管偏振敏感光学相干层析成像本地双折射解调方法 |
CN111965114B (zh) * | 2020-08-15 | 2023-08-29 | 天津大学 | 一种导管偏振敏感光学相干层析成像本地双折射解调方法 |
CN113888478A (zh) * | 2021-09-15 | 2022-01-04 | 天津大学 | 一种对血管内导管偏振敏感相干层析成像优化去消偏方法 |
CN113888478B (zh) * | 2021-09-15 | 2022-11-11 | 天津大学 | 一种对血管内导管偏振敏感相干层析成像优化去消偏方法 |
WO2023236406A1 (zh) * | 2022-06-09 | 2023-12-14 | 温州医科大学 | 背向散射Mueller矩阵的可逆极分解方法 |
CN115035210A (zh) * | 2022-08-10 | 2022-09-09 | 天津恒宇医疗科技有限公司 | 基于偏振多参量融合的ps-oct可视度提升方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109164048B (zh) | 一种对导管偏振敏感光学相干层析成像偏振解调方法 | |
CN110584613A (zh) | 一种导管偏振敏感光学相干层析成像系统及解调方法 | |
CN110742584A (zh) | 一种导管偏振敏感光学相干层析成像解调方法用偏振解算方法 | |
JP5135324B2 (ja) | サンプルの偏光感応性光周波数領域画像形成のための方法、構成およびシステム | |
US7903256B2 (en) | Methods, systems, and computer program products for performing real-time quadrature projection based Fourier domain optical coherence tomography | |
JP2015166747A (ja) | 角度分解型のフーリエドメイン光干渉断層撮影法を遂行する方法及びシステム | |
CN106137134B (zh) | 多角度复合的血流成像方法及系统 | |
EP3153843A1 (en) | Polarization sensitive optical image measurement system, and program loaded into said system | |
US20150216415A1 (en) | Method and apparatus for performing multidimensional velocity measurements using amplitude and phase in optical interferometry | |
CN112097678A (zh) | 基于频率盲估计的多表面面形测量方法 | |
CN104995530B (zh) | 被检体信息获取装置和被检体信息获取方法 | |
CN110954481B (zh) | 导管偏振敏感光学相干层析成像优化平均偏振解调方法 | |
CN211934010U (zh) | 一种导管偏振敏感光学相干层析成像系统 | |
CN110731755A (zh) | 一种导管偏振敏感光学相干层析成像系统的偏振调平方法 | |
Sathyanarayana et al. | Recovery of blood flow from undersampled photoacoustic microscopy data using sparse modeling | |
CN110768661A (zh) | 一种基于神经网络的锁相放大器 | |
CN110623644A (zh) | 一种导管偏振敏感光学相干层析成像解调方法用色散补偿方法 | |
CN110742582A (zh) | 一种导管偏振敏感光学相干层析成像解调用双态图像分割方法 | |
CN113888478B (zh) | 一种对血管内导管偏振敏感相干层析成像优化去消偏方法 | |
JP2010223670A (ja) | 光断層画像表示システム | |
CN110742583A (zh) | 一种导管偏振敏感光学相干层析成像解调用光谱整形方法 | |
CN111965114B (zh) | 一种导管偏振敏感光学相干层析成像本地双折射解调方法 | |
Zhu et al. | Displacement measurement using digital speckle multi-frequency harmonic wave correlation method | |
Gianto et al. | Study of robustness of 2D fringe processing in coherence scanning interferometry for the characterization of a transparent polymer film | |
Sudhakar Reddy et al. | Universal discrete finite rate of innovation scheme for sparse signal reconstruction |
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: 20200204 |