一种光相干层析成像系统图像运动伪差校正算法
【专利摘要】一种光相干层析成像系统图像运动伪差校正算法，包括：原始数据的获取：对OCT系统采集的原始数据进行预处理形成结构图像，对图像进行灰度化，将灰度化后的图像进行滤波预处理；对X方向的运动伪差进行校正；对Y方向的运动伪差进行校正；对小幅度旋转引起的运动伪差进行校正；再用归一化互相关算法对X方向和Y方向进行第二次对齐；在保持X方向和Y方向分别对齐的状态下，将图像恢复到原未拉平的形状；进行3D可视化处理。本发明可以对OCT系统在动态成像过程中产生的伪差进行校正，包括X方向、Y方向以及小幅度旋转伪差校正。能够更加真实的反应眼底视网膜中央凹附近的生理结构，降低临床诊断的误诊以及错诊，提高了临床诊断的准确性和科学性。
【专利说明】一种光相干层析成像系统图像运动伪差校正算法
【技术领域】
[0001] 本发明涉及一种光相干层析成像系统。特别是涉及一种光相干层析成像系统图像 运动伪差校正算法。
【背景技术】
[0002] 光相干层析成像（Optical Coherence Tomography, 0CT)技术可以对组织的内部 结构和形态进行二维和三维可视化。具有高分辨率、高灵敏度、活体、实时、非侵入、横向分 辨率和纵向分辨率相互独立等优点。自从0CT出现，它已经成了眼科疾病诊断和研究的一 种标准工具。频域0CT的出现极大的提高了成像的灵敏度和扫描速度，使0CT系统实时成 像成为可能。
[0003] 在获取眼底图像时，由于人眼的飘逸、眨眼、扫视以及头部等不自主运动，导致在 一些超高速的0CT系统中产生运动伪差。运动伪差产生图像畸变会使数据丧失连续性，三 维重建出的三维图像将会丢失一些组织的结构信息，使得视网膜的生理结构不能得到真实 的反应，实时成像过程中由于0CT系统的高度敏感性，甚至会导致成像的失败，从而导致一 些眼科疾病的误诊。所以图像运动伪差的校正在0CT技术中占据相当重要的地位，已经成 为了相关领域的研究热点。近年来有一些研究小组做了相关的研究。
[0004] 一般可以通过提高图像的扫描速度来降低运动伪差，而运动伪差在很大程度上取 决于受试者，在采集数据时，随着扫描时间的增加运动伪差也会增大。在理想状态下，即非 常快的体数据获取速度（小于looms)时会达到无运动伪差的模式。但是这种情况在目前 的眼科0CT设备是不可能实现，除非减少体数据中A扫描的数量，而减少A扫描的数量，则 会严重的限制眼底图像的细节，3D视网膜数据的获取和横向的抽样密度之间相互制约。在 扫描时可以像其他眼科成像设备一样固定下巴和前额减少头部运动，从而降低轴向的运动 伪差。使受试者凝视一个固定目标也可降低成像过程中的运动伪差。但是，这些方法都不 能完全的消除运动伪差，所以在0CT成像过程中运动伪差的校正是非常必须的。
[0005] 0CT工作原理，其主要由一个光纤式迈克尔逊干涉仪组成。光源一般选用低相干光 源或飞秒激光光源，指示光源作为初始的聚焦和准直光源。光源发出的低相干光先经过一 个2 X 1的光纤耦合器，再进入2 X 2的光纤耦合器后，被均匀地分成两束，分别照射在样品 臂的样品和参考臂上装的反光镜上。参考臂上的反光镜反射回来的参考光和样品臂上被样 品背景反射回来的信号光之间形成一定的光程，而参考臂上的参考镜可以轴向移动，以调 整两者之间的光程差，只有当其光程差处在光源的一个相干长度范围内时，这两束光才会 产生相干信号，而这种相干信号带有一定的组织背景信号。因此，只有在探测光束的焦点处 返回的光束才有最强的干涉信号。这种干涉信号被光电探测部分接收，再经过信号放大、滤 波、信号采集与处理，最后在计算机上呈现清晰的图像。
【发明内容】
[0006] 本发明所要解决的技术问题是，提供一种能够更加真实的反应眼底视网膜的生理 结构的光相干层析成像系统图像运动伪差校正算法。
[0007] 本发明所采用的技术方案是：一种光相干层析成像系统图像运动伪差校正算法， 包括如下步骤：
[0008] 1)原始数据的获取：
[0009] 2)对0CT系统采集的原始数据进行预处理形成结构图像，对图像进行灰度化，将 灰度化后的图像进行滤波预处理；
[0010] 3)对X方向的运动伪差进行校正；
[0011] 4)对Y方向的运动伪差进行校正；
[0012] 5)对小幅度旋转引起的运动伪差进行校正；
[0013] 6)再用归一化互相关算法对X方向和Y方向进行第二次对齐；
[0014] 7)在保持X方向和Y方向分别对齐的状态下，将图像恢复到原未拉平的形状；
[0015] 8)进行3D可视化处理。
[0016] 步骤3)所述的对X方向的运动伪差进行校正是，首先求取每一帧图像的质量重 心，以第二帧图像的质量重心为中心点画出一个像素为mXn大小的区域作为模板，以第一 帧图像为目标图像，用归一化互相关算法对第二帧图像进行模板匹配校正，依次类推，每帧 图像都以前一帧图像为目标图像进行匹配校正，将X方向对齐，从而使体数据集的整个X方 向对齐，其中，m、η分别是模板的列数和行数。
[0017] 步骤4)所述的对Υ方向的运动伪差进行校正是：对步骤3)X方向已经对齐的体数 据集提取与步骤1)中所述Κ个位置相对应的重建图像，然后寻找Κ帧重建图像和相应位置 直接扫描所获得的图像的设定层边界；再针对Κ个位置的Κ对相对应的图像的设定层边界 值分别做差值Ak[I]，得到的Κ组差值，对所得到的Κ组差值求平均值得到Α[Ι]，其中的平均 值A[i]为第i帧图像在Υ方向的位移量，依所得到的位移量对I帧图像进行Υ方向平移， 使原始数据在Y方向得到校正，其中A k[I]是在第k个位置相对应的两组图像的边界值的 差值，A[I] = 是体数据的第i帧图像的序号。
[0018] 步骤5)所述的对小幅度旋转引起的运动伪差进行校正是：寻找步骤4)校正后的 体数据集中的I帧图像的设定层边界，以所述的设定层边界最低点为基准，将图像的边界 拉平，完成小幅度旋转引起的运动伪差的校正。
[0019] 本发明的一种光相干层析成像系统图像运动伪差校正算法，在标准FD - 0CT系统 上，不需增加任何的追踪定位系统以及硬件的情况下，就可以对0CT系统在动态成像过程 中产生的伪差进行校正，包括X方向、Y方向以及小幅度旋转伪差校正。本发明能够更加真 实的反应眼底视网膜中央凹附近的生理结构。由于使眼底视网膜的生理结构得到了真实的 反应，从而可以降低临床诊断的误诊以及错诊，提高了临床诊断的准确性和科学性。同时也 为实现0CT动态成像中眼运动的实时追踪校正奠定了基础。
【专利附图】
【附图说明】
[0020] 图1是本发明方法的流程图；
[0021] 图2图像数据获取扫描示意图；
[0022] 图3沿着慢轴⑵快速获取K帧B扫描图像；
[0023] 图4是未采用本发明的方法进行校正的3D图像效果图；
[0024] 图5是采用本发明的方法进行校正后的3D图像效果图。
【具体实施方式】
[0025] 下面结合实施例和附图对本发明的一种光相干层析成像系统图像运动伪差校正 算法做出详细说明。
[0026] 在获取眼底图像时，原始数据（X方向）左右的位移使得眼底视网膜的成像的组织 错位，X方向的校正是对数据左右的位移进行补偿，校正X方向运动伪差产生的组织错位。 经过被测眼视轴方向（Y方向）运动伪差的校正，沿着慢轴（Z方向）重建出来的图像中很多 毛刺消失，表面变得光滑。对畸变的原始数据进行X和Y方向的校正之后，重建的3D图像 的表面不再是凹凸不平有褶皱，而是非常的光滑，几乎与真实的眼底视网膜图像没有差别， 视网膜en-face单层切片图像也变得非常清晰，更加真实的反应了眼底视网膜中央凹附近 的生理结构。
[0027] 如图1所示，本发明的一种光相干层析成像系统图像运动伪差校正算法，包括如 下步骤：
[0028] 1)原始数据的获取：
[0029] 首先利用0CT系统采集被测组织3D体数据块（如图2所示），获取3D体数据块的 大小为MXNXI。再沿着慢轴⑵扫描方向快速的获取K帧（如图3所示的）B扫描图像， 每帧图像的间隔为ΜΛΚ-1)个像素。其中M、N分别是沿着快轴（X)扫描获得的每帧图像像 素的列数和行数，I是体数据集的帧数，K是沿慢轴（Z)方向快速获取的B扫描图像帧数，为 大于或等于1小于或等于Μ的整数。
[0030] 2)对0CT系统采集的原始数据进行预处理形成结构图像，对图像进行灰度化，将 灰度化后的图像进行滤波预处理；
[0031] 3)对X方向的运动伪差进行校正；
[0032] 所述的对X方向的运动伪差进行校正是，首先求取每一帧图像的质量重心，以第 二帧图像的质量重心为中心点画出一个像素为mXn大小的区域作为模板，以第一帧图像 为目标图像，用归一化互相关算法对第二帧图像进行模板匹配校正，依次类推，每帧图像都 以前一帧图像为目标图像进行匹配校正，将X方向对齐，从而使体数据集的整个X方向对 齐，其中，m、η分别是模板的列数和行数。
[0033] 4)对Υ方向的运动伪差进行校正；
[0034] 所述的对Υ方向的运动伪差进行校正是：对步骤3)X方向已经对齐的体数据集提 取与步骤1)中所述Κ个位置相对应的重建图像，此时得到的Κ帧重建图像与在相应位置直 接扫描所获得的图像信息量是一致的，只是在每列像素 Υ方向上存在位移差。然后寻找Κ 帧重建图像和相应位置直接扫描所获得的图像的设定层边界（如ILM、RNFL、IPL、INL、OPL、 RPE_ChCap、IS/0S 等）；
[0035] 5)针对K个位置的K对相对应的图像的设定层边界值分别做差值Ak[I]，得到 的K组差值，对所得到的K组差值求平均值得到A[I]，其中的平均值A[i]为第i帧图像 在Y方向的位移量，依所得到的位移量对I帧图像进行Y方向平移，使原始数据在Y方 向得到校正。其中A k[I]是在第k个位置相对应的两组图像的边界值的差值，A[I]= (4[1]+· · · · +AK[I])/K，i是体数据的第i帧图像的序号；
[0036] 6)对小幅度旋转引起的运动伪差进行校正；
[0037] 所述的对小幅度旋转引起的运动伪差行校正是：寻找步骤5)校正后的体数据集 中的I帧图像的设定层边界（如ILM、RNFL、IPL、INL、OPL、RPE_ChCap、IS/OS等），以所述的 设定层边界最低点为基准，将图像的边界拉平，完成小幅度旋转引起的运动伪差的校正；
[0038] 7)再用归一化互相关算法对X方向和Y方向进行第二次对齐，
[0039] 8)在保持X方向和Y方向分别对齐的状态下，将图像恢复到原未拉平的形状；
[0040] 9)进行3D可视化处理。
【权利要求】
1. 一种光相干层析成像系统图像运动伪差校正算法，其特征在于，包括如下步骤： 1) 原始数据的获取： 2) 对OCT系统采集的原始数据进行预处理形成结构图像，对图像进行灰度化，将灰度 化后的图像进行滤波预处理； 3) 对X方向的运动伪差进行校正； 4) 对Y方向的运动伪差进行校正； 5) 对小幅度旋转引起的运动伪差进行校正； 6) 再用归一化互相关算法对X方向和Y方向进行第二次对齐； 7) 在保持X方向和Y方向分别对齐的状态下，将图像恢复到原未拉平的形状； 8) 进行3D可视化处理。
2. 根据权利要求1所述的一种光相干层析成像系统图像运动伪差校正算法，其特征在 于，步骤3)所述的对X方向的运动伪差进行校正是，首先求取每一帧图像的质量重心，以第 二帧图像的质量重心为中心点画出一个像素为mXn大小的区域作为模板，以第一帧图像 为目标图像，用归一化互相关算法对第二帧图像进行模板匹配校正，依次类推，每帧图像都 以前一帧图像为目标图像进行匹配校正，将X方向对齐，从而使体数据集的整个X方向对 齐，其中，m、η分别是模板的列数和行数。
3. 根据权利要求1所述的一种光相干层析成像系统图像运动伪差校正算法，其特征在 于，步骤4)所述的对Υ方向的运动伪差进行校正是：对步骤3) X方向已经对齐的体数据集 提取与步骤1)中所述Κ个位置相对应的重建图像，然后寻找Κ帧重建图像和相应位置直接 扫描所获得的图像的设定层边界；再针对Κ个位置的Κ对相对应的图像的设定层边界值分 别做差值A k[I]，得到的Κ组差值，对所得到的Κ组差值求平均值得到A[I]，其中的平均值 A[i]为第i帧图像在Y方向的位移量，依所得到的位移量对I帧图像进行Y方向平移，使原 始数据在Y方向得到校正，其中A k[I]是在第k个位置相对应的两组图像的边界值的差值， A[I] = %[1]+. . . .+AK[I])/K，i是体数据的第i帧图像的序号。
4. 根据权利要求1所述的一种光相干层析成像系统图像运动伪差校正算法，其特征在 于，步骤5)所述的对小幅度旋转引起的运动伪差进行校正是：寻找步骤4)校正后的体数据 集中的I帧图像的设定层边界，以所述的设定层边界最低点为基准，将图像的边界拉平，完 成小幅度旋转引起的运动伪差的校正。
【文档编号】A61B3/14GK104146682SQ201410455290
【公开日】2014年11月19日 申请日期:2014年9月9日 优先权日:2014年9月9日
【发明者】李跃杰, 张明蓉, 王立伟, 高振玉 申请人:中国医学科学院生物医学工程研究所