一种基于独立成分分析的全场光学相干显微成像去噪方法技术

本发明专利技术涉及一种基于独立成分分析的全场光学相干显微成像去噪方法，属于图像处理技术领域。该方法包括：采集由于系统低频振动产生运动伪影噪声的全场光学相干图像并进行预处理得到图像X

【技术实现步骤摘要】
一种基于独立成分分析的全场光学相干显微成像去噪方法
本专利技术属于图像处理
，涉及一种基于独立成分分析的全场光学相干显微成像去噪方法。
技术介绍
2002年法国Boccara团队基于时域OCT(OpticalCoherenceTomography)原理专利技术的全场光学相干层析成像，不需扫描，一次性凝视提取组织内部全视场图像，可原生环境下实现亚细胞分辨率细胞检测，穿透深度大、分辨率高、速度快，是一种非侵入性、非破坏性和内源性获取生物组织图像的方法。全场光学相干层析成像技术的原理是当样品被具有低相干长度的光源照射时由样品后向散射的光与光源干涉成像，该技术利用光源的低相干性来隔离由样品深度方向上的虚拟切片后向散射的光。干涉仪的使用使得可以通过干涉现象产生表示从样品的给定切片选择性地发出的光的干涉信号，并且消除源自样品的其余部分的光。全场光学相干层析成像可以获得具有微米量级分辨率的三维图像，使得可以观察到血管、血管壁、胶原蛋白、脂肪细胞等的大部分组织结构。同时，一次性凝视成像使其可以在几分钟内产生表面为几平方厘米，深度方向为几百微米的三维切片图像，可以潜在地用于诸如诊断癌症的医学应用。在手术期间的全场光学相干层析成像使得可以进行实时诊断，甚至在例如癌症的情况下优化手术操作，从而避免可能在切除后几天组织学分析揭示存在肿瘤细胞的情况下重复手术。然而，由于全场光学相干层析成像系统使用了迈克尔逊干涉仪原理，装置对振动十分敏感，周围环境微米量级的振动将使样品成像时产生轴向运动伪影噪声，严重影响图像质量。目前全场光学相干层析成像系统仍主要放置在实验室气浮隔振平台上使用，阻碍了其仪器化和设备化。由于周围环境振动噪声和样品光学相干信号具有相对独立性，且环境振动噪声主要集中在低频成分，因此，可以将样品光学相干信号和噪声分别看作互相独立的两个信号，利用独立成分分析原理将环境振动产生的轴向运动伪影噪声滤除。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种基于独立成分分析的全场光学相干显微成像去噪方法，利用独立成分分析原理将低频振动噪声和全场光学相干信号分离，去除系统振动导致样品轴向位移产生的运动伪影噪声，从而通过图像处理算法而不需要气浮隔振实验平台滤除系统噪声。为达到上述目的，本专利技术提供如下技术方案：一种基于独立成分分析的全场光学相干显微成像去噪方法，具有包括以下步骤：S1：获取由于系统低频振动产生运动伪影噪声的全场光学相干有噪图像xi，并对其预处理，即通过图像累加以增强有用信号并去除系统随机噪声；S2：对预处理后的图像X1进行二维小波变换，并通过置高频通道小波系数为零后重构，获得保留低频信息的全场光学相干有噪图像X2；S3：将图像X1和X2合并组成新的图像矩阵X3；S4：对图像矩阵X3进行去均值和白化预处理，获得白化后的图像矩阵Z；S5：将白化后的图像矩阵作为独立成分分析算法的输入进行计算，即求解混合矩阵W，使输出矩阵Y＝WZ的各行向量之间相互独立；S6：根据求解出的混合矩阵W，计算得到原始无噪干涉图像。进一步，步骤S1中，对图像xi预处理，即图像累加，其计算公式为：进一步，步骤S1中，通过全场光学相干显微成像装置获取全场光学相干图像，该成像装置由低相干光源(1)、成像透镜(3)、分光镜(4)、参考臂(5)、样品臂(6)、I/O驱动控制卡(7)和面阵CCD探测器(8)组成。进一步，步骤S2中，所述二维小波变换采用Daubechies8小波分解方法将原始有噪图像分解成不同频率范围的子图像集，包括LL低频系数部分、HL水平细节部分、LH垂直细节部分和HH高频系数部分。进一步，步骤S3中，将图像X1和X2合并组成新的图像矩阵X3，具体包括：设原始光学干涉无噪图像为s，系统振动时样品轴向位移产生的伪影噪声为n，则有：X1＝s+nX2＝a21s+a22n其中，{aij}为混合系数，i,j＝1,2，W为矩阵A的逆。进一步，所述步骤S4具体包括以下步骤：S41：去均值：以图像矩阵X3的某一行向量为对象，计算其元素的平均值，将该行向量的元素都减去此平均值，获得新向量，将矩阵其他行都按照相同方法处理，获得新矩阵X4；S42：白化预处理：寻找一个线性白化矩阵V，使得矩阵Z＝VX4的行向量相互正交，其中，X4为去均值后的矩阵，V＝D-1/2ET，D为矩阵X4的协方差矩阵E{X4X4T}的特征向量的对角矩阵，矩阵D-1/2通过简单的逐元素开方计算得到为E{X4X4T}的特征向量的正交矩阵。进一步，步骤S5中，求解混合矩阵W具体包括以下步骤：S51：随机初始化向量wi(0)，其元素个数为矩阵Z的行向量个数，i＝1,…,m，m为需要估计的独立成分的个数，且满足||wi(0)||＝1；S52：迭代：对每个i＝1,…,m，令k＝1，更新wi(k)＝E{z(wi(k-1)Tz)3}-3wi(k-1)，其中，z为矩阵Z的一列，E为矩阵z(wi(k-1)Tz)3的期望；S53：wi(k)归一化处理：wi(k)＝wi(k)/wi(k)，确保分离结果具有单位能量；S54：如果|wi(k)|趋近于1，则输出wi(k)，否则，令k＝k+1，并返回到步骤S52；S55：上述迭代过程结束后，得到的wi(k)再取转置为混合矩阵W的一个行向量，从而求得混合矩阵W。本专利技术的有益效果在于：本专利技术基于独立成分分析原理，利用小波变换产生一个保留低频伪影噪声信息的虚拟图像，和原图像一起作为独立成分分析算法的两个输入进行图像和噪声的分离，从而实现伪影噪声滤除。本专利技术仅依靠图像处理算法而不需要气浮隔振实验平台滤除系统噪声。本专利技术的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本专利技术的实践中得到教导。本专利技术的目标和其他优点可以通过下面的说明书来实现和获得。附图说明为了使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术作优选的详细描述，其中：图1为本专利技术基于独立成分分析的全场光学相干显微成像去噪方法的流程原理图；图2为全场光学相干层析成像装置原理示意图；图3为含运动伪影噪声的全场光学相干图像；图4为去除的运动伪影噪声图像；图5为去除噪声后的全场光学相干图像；附图标记：1-低相干光源，2-会聚透镜，3-成像透镜，4-分光镜，5-参考臂，6-样品臂，7-I/O驱动控制卡，8-面阵CCD探测器，9-补偿玻片，10-高数值孔径物镜。具体实施方式以下通过特定的具体实例说明本专利技术的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本专利技术的其他优点与功效。本专利技术还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本专利技术的精神下进行各种修饰或改本文档来自技高网...

【技术保护点】
1.一种基于独立成分分析的全场光学相干显微成像去噪方法，其特征在于，该方法具有包括以下步骤：/nS1：获取含有伪影噪声的全场光学相干图像x

【技术特征摘要】
1.一种基于独立成分分析的全场光学相干显微成像去噪方法，其特征在于，该方法具有包括以下步骤：
S1：获取含有伪影噪声的全场光学相干图像xi，并对其预处理；
S2：对预处理后的图像X1进行二维小波变换，并通过置高频通道小波系数为零后重构，获得保留低频信息的全场光学相干噪声图像X2；
S3：将图像X1和X2合并组成新的图像矩阵X3；
S4：对图像矩阵X3进行去均值和白化预处理，获得白化后的图像矩阵Z；
S5：将白化后的图像矩阵作为独立成分分析算法的输入进行计算，即求解混合矩阵W，使输出矩阵Y＝WZ的各行向量之间相互独立；
S6：根据求解出的混合矩阵W，计算得到原始去噪图像。


2.根据权利要求1所述的全场光学相干显微成像去噪方法，其特征在于，步骤S1中，对图像xi预处理，即图像累加，其计算公式为：


3.根据权利要求1所述的全场光学相干显微成像去噪方法，其特征在于，步骤S1中，通过全场光学相干显微成像装置获取全场干涉图像，该成像装置由低相干光源(1)、成像透镜(3)、分光镜(4)、参考臂(5)、样品臂(6)、I/O驱动控制卡(7)和面阵CCD探测器(8)组成。


4.根据权利要求1所述的全场光学相干显微成像去噪方法，其特征在于，步骤S2中，所述二维小波变换采用Daubechies8小波分解方法将原始有噪图像分解成不同频率范围的子图像集，包括LL低频系数部分、HL水平细节部分、LH垂直细节部分和HH高频系数部分。


5.根据权利要求1所述的全场光学相干显微成像去噪方法，其特征在于，步骤S3中，将图像X1和X2合并组成新的图像矩阵X3，具体包括：设原始光学干涉无噪图像为s，系统振动时样品轴向位移产生的伪影噪声...

【专利技术属性】
技术研发人员：熊亮，王金玉，尹韶云，杜凯，李俊，邹钱生，
申请(专利权)人：中国科学院重庆绿色智能技术研究院，
类型：发明
国别省市：重庆;50