[发明专利]一种光学相干断层扫描微血管成像方法、装置及设备

权利要求书

1.一种光学相干断层扫描微血管成像方法，其特征在于，包括以下步骤：

获取待成像的三维Bscan图像信号；

对所述Bscan图像信号进行快速傅里叶变换得到第一图像信号；

计算所述第一图像信号的体素复数方差；

根据所述体素复数方差去除所述第一图像信号的噪声得到第二图像信号；

根据所述第二图像信号生成二维微血管断层图像；

将所述二维微血管断层图像进行三维重建得到三维En-face图像。

2.根据权利要求1所述的方法，其特征在于：所述计算所述第一图像信号的体素复数方差，体素复数方差计算公式为：

，其中，，；

VCND为体素复数方差，N为同一位置扫描的总次数，z_n为第n次扫描时图像的体素，μ为同一位置所有体素的平均值，n^th为在第n次扫描时图像中预设顺序的序号；为B扫描图像数据集，i为B扫描图像横坐标，j为B扫描图像纵坐标；

将同一位置的所述体素复数方差求和得到总体素复数方差，计算公式如下：

其中，SSCVD为总体素复数方差。

3.根据权利要求2所述的方法，其特征在于：所述根据所述体素复数方差去除所述第一图像信号的噪声得到第二图像信号，包括：

使用窗口较大的median_filter去除血管，使血管与周围组织融为一体，得到组织掩码，方法如下：

组织掩码 ← median_filter (SSCVD，kernal_size=C) ；

根据得到的组织掩码去除组织保留血管，方法如下：

Filtered_SSCVD←SSCVD-组织掩码*α；

采用percentile对处理后的图像信号进行优化，方法如下：

Filtered_SSCVD [Filtered_SSCVD  阈值] = 阈值；

式中， 阈值← percentile（Filtered_SSCVD, β）；

其中，← 表示其经过右边操作后得到左边，median_filter为算法，percentile为python库，SSCVD为总体素复数方差，kernal_size表示内核尺寸，C为预设的具体数值；Filtered_为算法，α是预设权重，α大于0且小于1，β为预设值。

4.根据权利要求1所述的方法，其特征在于：所述获取待成像的三维Bscan图像信号，包括：

获取维度为Y×Z×X×N的四维OCTA原始扫描信号，其中，Z 是Ascan 深度，X 是Bscan的横向长度，N 是相同位置 Bscan 的重复次数，Y是C scan 中的纵向宽度和生成的三维Bscan图像信号的总数；

将所述四维OCTA原始扫描信号分割成Y个维度为Z×X×N的三维Bscan图像信号；

将提取的每一个三维Bscan图像信号进行去直流信号干扰；

将去直流信号干扰后的三维Bscan图像信号进行窗函数计算；

将计算结果进行色散校正后得到待成像的三维Bscan图像信号。

5.根据权利要求1所述的方法，其特征在于：所述对所述Bscan图像信号进行快速傅里叶变换得到第一图像信号，包括：

将所述Bscan图像信号直接进行快速傅里叶变换得到第一图像信号。

6.根据权利要求1所述的方法，其特征在于：所述对所述Bscan图像信号进行快速傅里叶变换得到第一图像信号，包括：

将所述Bscan图像信号进行划分，生成至少两个分裂谱图像信号；

将每个分裂谱图像信号进行快速傅里叶变换得到一个第一图像信号。

7.根据权利要求1所述的方法，其特征在于：所述对所述Bscan图像信号进行快速傅里叶变换得到第一图像信号后，还包括：

对所述第一图像信号进行运动校准和相位补偿处理。