[发明专利]一种自动分离动静脉血管的内源光学成像方法

权利要求书

1.一种自动分离动静脉血管的内源光学成像方法，其特征在于包括以下步骤：

第一步，将窄带准单色光照射到被测对象上，用CCD或CMOS相机通过光学成像系统，以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的光学图像，每帧图像的曝光时间不高于100ms，帧间隔时间?T不高于500ms，采集的帧数N≥20；

第二步，按公式一计算平均时间序列                                               中第n个点的平均灰度值：

   公式一，

其中为第n幅图像的第i个像素的灰度值，为N帧图像所有像素的平均灰度值， R为图像矩阵的行数，C为图像矩阵的列数，单幅图像的尺寸为I_R×C，该图像共有R×C个像素，相应的图像时间序列中包含R×C个像素时间序列；

第三步，利用快速傅里叶变换算法将第二步所得的平均像素时间序列变换到频谱域，确定频率区间(0.05Hz, 1Hz)中具有最大幅值的频率值f_m，f_m位于0.1Hz附近，称为0.1-Hz低频振荡；

第四步，利用快速傅里叶变换算法将每个像素时间序列变换到频谱域，获得f_m频率处的特征频率幅值A_m(r , c)，其中r为像素所在位置的行，c为像素所在位置的列；以每个像素时间序列对应的特征频率幅值A_m(r , c)为灰度，组成二维的特征频率幅值图A_m；

第五步，利用快速傅里叶变换算法将每个像素时间序列变换到频谱域，计算频率区间[1Hz, 1/(2??T)]中所有离散频率点的谱功率总和W_g(r , c)，其中各频率点的谱功率为相应频率幅值的平方；计算特征频率比值R_a(r , c)，R_a(r , c)为各像素序列对应的特征频率f_m的谱值(A_m(r , c))²与W_g(r , c)的比值，R_a(r , c)=(A_m(r , c))²/W_g(r , c)；以各像素序列对应的特征频率比值R_a(r , c)为灰度，组成二维的特征频率比值图R_a；

第六步，利用局部对比度增强方法对由第四步和第五步得到的特征频率幅值图A_m和特征频率比值图R_a做局部对比度增强，具体操作如下：

以图像中第r行，第c列处的像素为中心，选取像素大小为49×49的空间滑动窗W，按公式二得到中心像素的新灰度值p(r, c)，

  公式二

其中，p_min为待操作图像所有像素灰度值的最小值，p_max为待操作图像所有像素灰度值的最大值，按公式三计算， 

 公式三

其中和分别为W滑动窗所覆盖的所有像素的灰度平均值和方差；

遍历图像A_m和R_a中所有的像素，将每个像素的灰度值按步骤6.1的方法进行更新，得到局部对比度增强后的特征频率幅值图I_Am和特征频率比值图I_Ra；

第七步，使用阈值分割方法将第六步中所得局部对比度增强后的特征频率幅值图I_Am作二值化处理，得到二维动脉血管结构图像I_artery，其中动脉血管区域像素值为1，其他区域像素值为0；

第八步，使用阈值分割方法将第六步中所得图像I_Ra作二值化处理，得到二维静脉血管结构图像I_vein，其中静脉血管区域像素值为1，其他区域像素值为0。

2.如权利要求1所述的一种自动分离动静脉血管的内源光学成像方法，其特征在于将图像I_artery和图像I_vein按如下方式构建新的RGB图像C_RGB，动脉区域和静脉区域分别用红色R、绿色G、蓝色B三基色中的一种来表示，设置动脉区域的颜色值为红色：（1，0，0），静脉区域的颜色值为蓝色（0，0，1），其他区域的颜色值为白色（1，1，1）；从而绘制出一幅以不同颜色标记动静脉血管的RGB图像。