基于sd-oct视网膜图像的地图状萎缩投影图像生成方法
【专利摘要】本发明提供一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其实现包括以下步骤：步骤1、采集SD-OCT视网膜图像；步骤2、分割步骤1所得图像的布鲁赫膜边界；步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。本发明的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法考虑和克服了脉络膜层血管对GA投影图像的影响，提高了GA投影图像的对比度和灰度一致性。
【专利说明】基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法
【技术领域】
[0001]本发明涉及投影图像生成和处理【技术领域】，具体而言涉及一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法。
【背景技术】
[0002]SD-OCT(频域光学相干断层扫描)视网膜图像是一种频域光学相干断层成像图像，它可以有效地呈现视网膜组织层的灰度和结构等变化，临床实验表明SD-OCT图像能够用于测量地图状萎缩(GA)的范围。由于SD-OCT图像是三维的，不利于医生直观而全面的观测GA的范围，因此需要将三维的SD-OCT图像投影成二维的眼底图像。现有的GA投影图像生成方法主要是基于GA区域脉络膜灰度值高于正常脉络膜区域的特性。
[0003]现有技术中主要有如下两种GA投影图像生成方法:
[0004](I)基于像素求和的方法。该方法通过将轴方向的所有像素灰度求和GA投影图像，该方法没有排除视网膜神经纤维层(RNFL)等其他高亮视网膜结构，导致GA投影图像的对比度较低。
[0005](2)基于区域限制的像素求和方法。该方法是对上述第一种方法的不足加以改进的方法，将投影区域限制在视网膜色素上皮层(RPE)的下方，以排除RNFL等高亮区域的影响，但该方法没有考虑RPE层下方脉络膜层内血管的影响。由于脉络膜层含有大量的血管，血管在SD-OCT图像上表现为灰度值较低的暗区，这对GA投影图像的对比度和灰度一致性造成了影响。
[0006]因此，现有的GA投影图像生成方法没有考虑脉络膜层血管的影响，导致生成的GA投影图像对比度不高。
【发明内容】
[0007]本发明目的在于提供一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，克服上述现有技术中的缺陷。
[0008]本发明的上述目的通过独立权利要求的技术特征实现，从属权利要求以另选或有利的方式发展独立权利要求的技术特征。
[0009]为达成上述目的，本发明所采用的技术方案如下:
[0010]一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其实现包括以下步骤:
[0011]步骤1、采集SD-OCT视网膜图像；
[0012]步骤2、分割步骤I所得图像的布鲁赫膜边界；
[0013]步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；
[0014]步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及
[0015]步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。
[0016]进一步的实施例中，前述步骤2中，采用基于三维图搜索方法来分割布鲁赫膜边界。
[0017]进一步的实施例中，前述步骤3中，如果当前像素点的灰度值大于该点后续所有像素点的灰度值，则该点即为至上而下灰度值递减的像素点。
[0018]进一步的实施例中，前述步骤4的实现包括以下步骤:
[0019]步骤4-1、直线连接由步骤3得到的相邻的灰度值递减的像素点，从而构成一曲面；
[0020]步骤4-2、将曲面与横坐标之间的面积作为GA投影图像的值。
[0021]进一步的实施例中，前述步骤5中，采用模板大小为3X3中值滤波来平滑前述步骤4生成的GA投影图像，得到最终的结果。
[0022]由以上本发明的技术方案可知，本发明的有益效果在于:考虑和克服了脉络膜层血管对GA投影图像的影响，提高了 GA投影图像的对比度和灰度一致性。
【专利附图】
【附图说明】
[0023]图1是本发明一实施方式基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法的流程图。
[0024]图2是由灰度值递减像素点构成曲面下方面积作为GA投影图像的值的流程图。
[0025]图3是SD-OCT视网膜图像内部组织示意图。
[0026]图4是BM边界下方0.5毫米厚的窄带区域示意图。
[0027]图5是窄带区域的灰度值变化曲线示意图。
[0028]图6是本发明方法得到的最终的GA投影图像。
[0029]图7是基于像素求和方法得到的GA投影图像。
[0030]图8是基于区域限制的像素求和方法得到的GA投影图像。
[0031]图9是本发明方法与像素求和方法、基于区域限制的像素求和方法得到的GA投影图像的局部放大比较图。
【具体实施方式】
[0032]为了更了解本发明的技术内容，特举具体实施例并配合所附图式说明如下。
[0033]如图1所示，根据本发明的较佳实施方式，一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其实现包括以下步骤:
[0034]步骤1、采集SD-OCT视网膜图像；
[0035]步骤2、分割步骤I所得图像的布鲁赫膜边界(即BM边界)；
[0036]步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点；
[0037]步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及
[0038]步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。
[0039]本实施例中，以SD-OCT视网膜图像作为输入，采用图像处理手段由三维的输入图像生成二维的地图状萎缩投影图像。图3给出了 SD-OCT视网膜图像内部组织的示意图。
[0040]作为可选的实施方式，前述步骤2中，采用基于三维图搜索方法来分割布鲁赫膜边界。
[0041]在一些实施例中，可以采用自动的分割方式，也可以采用半自动的方式来进行分害I]，例如计算机程序自动分割+修正的方式，当然并不以此为限制。
[0042]结合图3、4所示，图3中示例性地给出了 BM边界的示意图，图4为采用计算机程序自动分割+修正方式分割BM边界后BM边界下方0.5毫米厚的窄带区域示意图。
[0043]在步骤3中，如果当前像素点的灰度值大于该点后续所有像素点的灰度值，则该点即为至上而下灰度值递减的像素点。
[0044]结合图2所示，作为可选的实施方式，前述步骤4的实现包括以下步骤:
[0045]步骤4-1、直线连接由步骤3得到的相邻的灰度值递减的像素点，从而构成一曲面；
[0046]步骤4-2、将曲面与横坐标之间的面积作为GA投影图像的值从而生成投影图像。
[0047]本实施例中，前述步骤5中，采用模板大小为3X3中值滤波来平滑前述步骤4生成的GA投影图像，得到最终的结果，该3X3大小的模板能够在去除噪声的同时保持目标的边缘。
[0048]下面结合图3-图9所示，对本发明上述实施例及其具体应用实例进行进一步说明。
[0049]结合图1所示的流程，通过OCT成像设备采集到SD-OCT视网膜图像，作为一个示例，其采集三维SD-OCT视网膜图像大小为1024X512X128，对应视网膜2mmX6mmX6mm的区域。
[0050]图3给出了一帧SD-OCT视网膜图像的感兴趣区域，图中标注了视网膜的几个主要相关组织结构，即=RPE:视网膜色素上皮层，BM:布鲁赫膜，Choroidal vessel:脉络膜血管，Choroid:脉络膜，CS1:脉络膜巩膜边界，GA:地图状萎缩，Sclera:巩膜。这些组织机构及其名称和英文缩写是本领域所熟知的。
[0051]图4为基于三维图搜索方法分割得到的BM边界下方0.5毫米(约250个像素)的窄带区域，然后定位每列中至上而下(即沿X轴方向)灰度值递减的像素点。将灰度值递减像素点构成曲面下方的多边形面积作为GA投影图像的值。
[0052]图5示例性地给出了第100列窄带区域的灰度值变化曲线，曲线上的带黑色星号的圆圈即为定位得到的至上而下灰度值递减的像素点，这些像素点用直线连接而成的曲线下方面积即作为GA投影图像的值，以此生成GA投影图像。
[0053]最后利用模板大小为3X3的中值滤波平滑所生成的GA投影图像即可得到最终的结果，如图6所示。
[0054]图7和图8分别为现有技术中基于像素求和方法和基于区域限制的像素求和方法得到的GA投影图像，图9为三种方法得到GA投影图像的局部放大(如图8中的黑色虚线框)比较。从图6-图9可知:本发明得到的GA投影图像能够在二维眼底图像中清晰地显示GA区域，且克服了脉络膜层血管对GA投影的影响(如图9中的左列区域所示)，GA区域的内部灰度一致性和对比度都优于现有的方法，这为临床医生的疾病诊断和后续的GA分割提供了方便。
[0055]虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属【技术领域】中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。
【权利要求】
1.一种基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，其实现包括以下步骤: 步骤1、采集SD-OCT视网膜图像； 步骤2、分割步骤I所得图像的布鲁赫膜边界； 步骤3、在布鲁赫膜边界下方的0.5毫米厚的窄带区域，定位每列自上而下中灰度值递减的像素点； 步骤4、将步骤3得到的灰度值递减像素点构成一曲面，并将曲面与横坐标之间的面积作为GA投影图像的值生成投影图像；以及 步骤5、用中值滤波平滑前述步骤4生成的GA投影图像得到最终的结果。
2.根据权利要求1所述的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，前述步骤2中，采用基于三维图搜索方法来分割布鲁赫膜边界。
3.根据权利要求1所述的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，前述步骤3中，如果当前像素点的灰度值大于该点后续所有像素点的灰度值，则该点即为至上而下灰度值递减的像素点。
4.根据权利要求1所述的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，前述步骤4的实现包括以下步骤: 步骤4-1、直线连接由步骤3得到的相邻的灰度值递减的像素点，从而构成一曲面； 步骤4-2、将曲面与横坐标之间的面积作为GA投影图像的值。
5.根据权利要求1所述的基于SD-OCT视网膜图像的地图状萎缩投影图像生成方法，其特征在于，前述步骤5中，采用模板大小为3 X 3中值滤波来平滑前述步骤4生成的GA投影图像，得到最终的结果。
【文档编号】G06T7/00GK104050672SQ201410295666
【公开日】2014年9月17日 申请日期:2014年6月26日 优先权日:2014年6月26日
【发明者】陈强, 牛四杰, 沈宏烈, 陆圣陶, 王玉萍 申请人:南京理工大学