基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法技术

本发明专利技术提供一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其实现包括以下步骤：步骤1、采集SD-OCT视网膜图像；步骤2、分割步骤1所得图像的布鲁赫膜边界；步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。本发明专利技术的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法考虑和克服了脉络膜层血管对GA投影图像的影响，提高了GA投影图像的对比度和灰度一致性。

【技术实现步骤摘要】

本专利技术涉及投影图像生成和处理
，具体而言涉及一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法。
技术介绍
SD-OCT(频域光学相干断层扫描)视网膜图像是一种频域光学相干断层成像图像，它可以有效地呈现视网膜组织层的灰度和结构等变化，临床实验表明SD-OCT图像能够用于测量地图状萎缩(GA)的范围。由于SD-OCT图像是三维的，不利于医生直观而全面的观测GA的范围，因此需要将三维的SD-OCT图像投影成二维的眼底图像。现有的GA投影图像生成方法主要是基于GA区域脉络膜灰度值高于正常脉络膜区域的特性。现有技术中主要有如下两种GA投影图像生成方法：(1)基于像素求和的方法。该方法通过将轴方向的所有像素灰度求和GA投影图像，该方法没有排除视网膜神经纤维层(RNFL)等其他高亮视网膜结构，导致GA投影图像的对比度较低。(2)基于区域限制的像素求和方法。该方法是对上述第一种方法的不足加以改进的方法，将投影区域限制在视网膜色素上皮层(RPE)的下方，以排除RNFL等高亮区域的影响，但该方法没有考虑RPE层下方脉络膜层内血管的影响。由于脉络膜层含有大量的血管，血管在SD-OCT图像上表现为灰度值较低的暗区，这对GA投影图像的对比度和灰度一致性造成了影响。因此，现有的GA投影图像生成方法没有考虑脉络膜层血管的影响，导致生成的GA投影图像对比度不高。
技术实现思路
本专利技术目的在于提供一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，克服上述现有技术中的缺陷。本专利技术的上述目的通过独立权利要求的技术特征实现，从属权利要求以另选或有利的方式发展独立权利要求的技术特征。为达成上述目的，本专利技术所采用的技术方案如下：一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其实现包括以下步骤：步骤1、采集SD-OCT视网膜图像；步骤2、分割步骤1所得图像的布鲁赫膜边界；步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。进一步的实施例中，前述步骤2中，采用基于三维图搜索方法来分割布鲁赫膜边界。进一步的实施例中，前述步骤3中，如果当前像素点的灰度值大于该点后续所有像素点的灰度值，则该点即为至上而下灰度值递减的像素点。进一步的实施例中，前述步骤4的实现包括以下步骤：步骤4-1、直线连接由步骤3得到的相邻的灰度值递减的像素点，从而构成一曲面；步骤4-2、将曲面与横坐标之间的面积作为GA投影图像的值。进一步的实施例中，前述步骤5中，采用模板大小为3×3中值滤波来平滑前述步骤4生成的GA投影图像，得到最终的结果。由以上本专利技术的技术方案可知，本专利技术的有益效果在于：考虑和克服了脉络膜层血管对GA投影图像的影响，提高了GA投影图像的对比度和灰度一致性。附图说明图1是本专利技术一实施方式基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法的流程图。图2是由灰度值递减像素点构成曲面下方面积作为GA投影图像的值的流程图。图3是SD-OCT视网膜图像内部组织示意图。图4是BM边界下方0.5毫米厚的窄带区域示意图。图5是窄带区域的灰度值变化曲线示意图。图6是本专利技术方法得到的最终的GA投影图像。图7是基于像素求和方法得到的GA投影图像。图8是基于区域限制的像素求和方法得到的GA投影图像。图9是本专利技术方法与像素求和方法、基于区域限制的像素求和方法得到的GA投影图像的局部放大比较图。具体实施方式为了更了解本专利技术的
技术实现思路
，特举具体实施例并配合所附图式说明如下。如图1所示，根据本专利技术的较佳实施方式，一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其实现包括以下步骤：步骤1、采集SD-OCT视网膜图像；步骤2、分割步骤1所得图像的布鲁赫膜边界(即BM边界)；步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。本实施例中，以SD-OCT视网膜图像作为输入，采用图像处理手段由三维的输入图像生成二维的地图状萎缩投影图像。图3给出了SD-OCT视网膜图像内部组织的示意图。作为可选的实施方式，前述步骤2中，采用基于三维图搜索方法来分割布鲁赫膜边界。在一些实施例中，可以采用自动的分割方式，也可以采用半自动的方式来进行分割，例如计算机程序自动分割+修正的方式，当然并不以此为限制。结合图3、4所示，图3中示例性地给出了BM边界的示意图，图4为采用计算机程序自动分割+修正方式分割BM边界后BM边界下方0.5毫米厚的窄带区域示意图。在步骤3中，如果当前像素点的灰度值大于该点后续所有像素点的灰度值，则该点即为至上而下灰度值递减的像素点。结合图2所示，作为可选的实施方式，前述步骤4的实现包括以下步骤：步骤4-1、直线连接由步骤3得到的相邻的灰度值递减的像素点，从而构成一曲面；步骤4-2、将曲面与横坐标之间的面积作为GA投影图像的值从而生成投影图像。本实施例中，前述步骤5中，采用模板大小为3×3中值滤波来平滑前述步骤4生成的GA投影图像，得到最终的结果，该3×3大小的模板能够在去除噪声的同时保持目标的边缘。下面结合图3-图9所示，对本专利技术上述实施例及其具体应用实例进行进一步说明。结合图1所示的流程，通过OCT成像设备采集到SD-OCT视网膜图像，作为一个示例，其采集三维SD-OCT视网膜图像大小为1024×512×128，对应视网膜2mm×6mm×6mm的区域。图3给出了一帧SD-OCT视网膜图像的感兴趣区域，图中标注了视网膜的几个主要相关组织结构，即：RPE：视网膜色素上皮层，BM：布鲁赫膜，Choroidal vessel：脉络膜血管，Choroid：脉络膜，CSI：脉络膜巩膜边界，GA：地图状萎缩，Sclera：巩膜。这些组织机构及其名称和英文缩写是本领域所熟知的。图4为基于三维图搜索方法分割得到的BM边界下方0.5毫米(约250个像素)的窄带区域，然后定位每列中至上而下(即沿x轴方向)灰度值递减的像素点本文档来自技高网...

【技术保护点】
一种基于SD‑OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，其实现包括以下步骤：步骤1、采集SD‑OCT视网膜图像；步骤2、分割步骤1所得图像的布鲁赫膜边界；步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。

【技术特征摘要】
1.一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，其实现
包括以下步骤：
步骤1、采集SD-OCT视网膜图像；
步骤2、分割步骤1所得图像的布鲁赫膜边界；
步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减
的像素点；
步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积
作为GA投影图像的值生成投影图像；以及
步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。
2.根据权利要求1所述的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其
特征在于，前述步骤2中，采用基于三维图搜索方法来分割布鲁赫膜边界。
3.根据权利要求1所...

【专利技术属性】
技术研发人员：陈强，牛四杰，沈宏烈，陆圣陶，王玉萍，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏;32