[发明专利]视网膜OCT体数据三维层状边界的提取方法及系统

权利要求书

1.一种视网膜OCT体数据三维层状边界的提取方法，其特征在于，包括以下步骤：

S1：对视网膜OCT三维体数据进行竖直方向的差值滤波运算，以获得差值体数据；

S2：通过所述差值体数据、以及所述视网膜OCT三维体数据的各个像素在竖直方向的坐标值进行计算，以获得边界指标体数据；

S3：在所述边界指标体数据中进行查找，以获得当前边界的离散边界点；

S4：对所述离散边界点的集合进行平滑操作，以获得当前边界；

S5：判断是否已经获得了所述视网膜OCT三维体数据的所有边界，若是，则执行步骤S6，否则，根据当前边界更新所述视网膜OCT三维体数据，并更新所述差值滤波运算的卷积算子，返回步骤S1；

S6：输出所述视网膜OCT三维体数据的所有边界。

2.如权利要求1所述的方法，其特征在于，步骤S2中，所述边界指标体数据通过下式进行计算，

I_i，j，k＝w₁*D_i，j，k+w₂*Y_i，j，k

其中，Y_i，j，k为视网膜OCT三维体数据，D_i，j，k为差值体数据，k为所述视网膜OCT三维体数据的各个像素在竖直方向的坐标值，w₁，w₂均为与竖直坐标k成正比的非负实数，I_i，j，k为边界指标体数据。

3.如权利要求1所述的方法，其特征在于，步骤S2之前还包括以下步骤：

S0：对所述原视网膜OCT三维体数据进行平滑运算，以获得平滑体数据；

步骤S2中，所述边界指标体数据通过下式进行计算，

I_i，j，k＝w₁*D_i，j，k+w₂*S_i，j，k

其中，S_i，j，k为平滑体数据，D_i，j，k为差值体数据，k为所述视网膜OCT三维体数据的各个像素在竖直方向的坐标值，w₁，w₂均为与竖直坐标k成正比的非负实数，I_i，j，k为边界指标体数据。

4.如权利要求1～3中任一项所述的方法，其特征在于，所述视网膜OCT三维体数据的所有边界包括：RPE层边界、IS/OS层边界和ILM层边界。

5.如权利要求4所述的方法，其特征在于，步骤S1中，所述差值滤波运算具体包括以下步骤：

S11：以所述视网膜OCT三维体数据的某个像素为中心，与所述视网膜OCT三维体数据的其它相邻像素建立一个M×M×M个像素点组成的正方体，所述M为不等于1的正奇数；

S12：根据当前边界对应的卷积算子，使所述正方体的中心之上的像素的亮度值之和减去所述正方体的中心之下的像素的亮度值之和，并以相减值的平均值作为所述正方体的中心的亮度值；

S13：以其它像素作为所述正方体的中心，返回步骤S11，直至获得所述视网膜OCT三维体数据的所有像素的亮度值，所述所有像素的亮度值构成差值体数据。