[发明专利]一种基于OCT扫描系统的三维血管类型区分方法

权利要求书

1.一种基于OCT扫描系统的三维血管类型区分方法，其特征在于，包括：

步骤1、通过OCT扫描系统对血管所在的区域进行连续的C扫描，得到血管在时间方向的三维立体图组；

步骤2、从所述三维立体图组中抽取单位深度所对应的二维图，并将所述二维图转成二维矩阵，将所述二维矩阵和其所对应的时间组成三维矩阵T；

步骤3、对所述三维矩阵T进行降噪处理，得到矩阵G；

步骤4、对矩阵G进行傅里叶变换将矩阵G从时域空间转换为频域空间，得到矩阵K；

步骤5、对矩阵K通过参量IR进行量化处理，得到矩阵A，所述矩阵A所表示的图像即为动脉分布图；

步骤6、将矩阵T减去矩阵A，得到矩阵B；

步骤7、对矩阵B以高斯函数作为卷积核进行卷积，得到矩阵C；

步骤8、计算矩阵C的Hessian矩阵H(x,y)；

步骤9、计算Hessian矩阵H(x,y)的特征值和对应的特征向量，并把两个特征向量的夹角方向作为血管的方向，旋转90度后得到血管径向搜索方向；

步骤10、通过血管径向搜索方向搜索，得到血管直径；

步骤11、通过血管直径得到静脉分布图和毛细血管分布图；

其中，所述参量IR的表达式为：

I(υ_f0)为I(υ)在特定频率f₀处的强度值，I(υ₀)为I(υ)在零频处的强度值，I(υ)表示为矩阵K每一点的强度值，所述特定频率f₀为预先设置；

在步骤10中，血管直径的确定方法包括：进行血管分割，然后进行二值化，沿着血管径向搜索方向和反方向分别开始寻找血管边界，相邻血管边界之间的距离即为血管直径，该血管直径记为第一直径；采用梯度极值方法计算沿血管径向搜索方向和其的反方向的梯度值变化情况，确定血管边缘，血管边缘梯度值从零到梯度增高，将这个突变的起始端和终止端作为血管的两个边缘点，将这两个边缘点的欧式距离作为血管直径，该血管直径记为第二直径；将第一直径和第二直径求平均得到血管直径。

2.根据权利要求1所述的一种基于OCT扫描系统的三维血管类型区分方法，其特征在于，在步骤3中，所述降噪处理为八领域降噪处理。

3.根据权利要求2所述的一种基于OCT扫描系统的三维血管类型区分方法，其特征在于，所述八领域降噪处理具体为：采用了一个值全为1的3*3的窗口轮流对矩阵T做卷积，之后沿时间方向依次对单列八邻域使用Savitzky-Golay滤波算子平滑数据。

4.根据权利要求1所述的一种基于OCT扫描系统的三维血管类型区分方法，其特征在于，在步骤11中，静脉分布图和毛细血管分布图的获取方式包括：将血管直径＞15μm的血管标识出来，得到静脉分布图，将血管直径≤15μm的血管标记出来，得到毛细血管分布图。