[发明专利]基于OCTA视网膜图像的黄斑中心无血管区体积计算方法及系统

权利要求书

1.一种基于OCTA视网膜图像的黄斑中心无血管区体积计算方法，其特征在于，所述方法包括以下步骤：

步骤1，采集OCTA视网膜图像；

步骤2，设置包含黄斑中心无血管区域即FAZ区域的限制区域；具体过程包括：

步骤2-1，分割OCTA视网膜图像的内界膜ILM和外丛状层OPL下边界；

步骤2-2，将ILM作为限制区域的上边界，将OPL下边界下方lμm处作为限制区域的下边界；

步骤3，拉平所述限制区域；具体过程包括：

步骤3-1，获取OPL下边界三维全局最低点，并将该最低点下方lμm处的行坐标作为下边界基准；

步骤3-2，将限制区域内的每一列移动至下边界基准，实现每一列下边界对齐，获得初步拉平限制区域；

步骤3-3，获取所述初步拉平限制区域的上边界三维全局最高点，并将该最高点作为上边界基准；

步骤3-4，将所述初步拉平限制区域内的每一列变换为固定长度，该固定长度为所述下边界基准、上边界基准之间的距离长度；

步骤4，在拉平后的限制区域内在轴向方向依次生成若干帧一个像素厚度的投影图像，并在每一帧投影图像上标注FAZ区域，记为标签图像；

步骤5，将所述标签图像还原至未拉平限制区域前的原始层位置；具体过程包括：

反向执行步骤3的反操作，将标签图像还原为符合OCTA视网膜图像原始层结构分布的标签图像；

步骤6，统计还原后的标签图像中标签像素的个数，由此求取FAZ区域的体积。

2.根据权利要求1所述的基于OCTA视网膜图像的黄斑中心无血管区体积计算方法，其特征在于，所述l的取值为OCTA图像采集设备划定的深层血管复合体DVC下边界值。

3.根据权利要求1所述的基于OCTA视网膜图像的黄斑中心无血管区体积计算方法，其特征在于，步骤3-4中将所述初步拉平限制区域内的每一列变换为固定长度，具体采用线性插值法。

4.根据权利要求1所述的基于OCTA视网膜图像的黄斑中心无血管区体积计算方法，其特征在于，步骤6中所述求取FAZ区域的体积，所用公式为：

V＝n×p

式中，V表示FAZ区域的体积，n表示标签像素的个数，p表示图像的分辨率参数。

5.一种基于OCTA视网膜图像的黄斑中心无血管区体积计算系统，其特征在于，所述系统包括：

采集模块，用于采集OCTA视网膜图像；

区域设置模块，用于设置包含黄斑中心无血管区域即FAZ区域的限制区域；具体过程包括：

步骤2-1，分割OCTA视网膜图像的内界膜ILM和外丛状层OPL下边界；

步骤2-2，将ILM作为限制区域的上边界，将OPL下边界下方lμm处作为限制区域的下边界；

区域变换模块，用于拉平所述限制区域；具体过程包括：

步骤3-1，获取OPL下边界三维全局最低点，并将该最低点下方lμm处的行坐标作为下边界基准；

步骤3-2，将限制区域内的每一列移动至下边界基准，实现每一列下边界对齐，获得初步拉平限制区域；

步骤3-3，获取所述初步拉平限制区域的上边界三维全局最高点，并将该最高点作为上边界基准；

步骤3-4，将所述初步拉平限制区域内的每一列变换为固定长度，该固定长度为所述下边界基准、上边界基准之间的距离长度；

标签图像生成模块，用于在拉平后的限制区域内在轴向方向依次生成若干帧一个像素厚度的投影图像，并在每一帧投影图像上标注FAZ区域，记为标签图像；

图像还原模块，用于将所述标签图像还原至未拉平限制区域前的原始层位置；具体过程包括：

反向执行区域变换模块的反操作，将标签图像还原为符合OCTA视网膜图像原始层结构分布的标签图像；

体积计算模块，用于统计还原后的标签图像中标签像素的个数，由此求取FAZ区域的体积。