《基于SDOCT的挫伤性视网膜内外节缺失检测方法.pdf》由会员分享，可在线阅读，更多相关《基于SDOCT的挫伤性视网膜内外节缺失检测方法.pdf（21页完整版）》请在专利查询网上搜索。
1、10申请公布号CN104143087A43申请公布日20141112CN104143087A21申请号201410356040X22申请日20140724G06K9/00200601G06K9/36200601G06K9/46200601A61B3/1020060171申请人苏州大学地址215000江苏省苏州市工业园区仁爱路199号72发明人陈新建朱伟芳张立春陈浩宇石霏向德辉王莉芸张莉高恩婷74专利代理机构南京纵横知识产权代理有限公司32224代理人董建林54发明名称基于SDOCT的挫伤性视网膜内外节缺失检测方法57摘要本发明公开了一种基于SDOCT图像的挫伤性视网膜内节/外节缺失三维自动检测。
2、方法，包括以下步骤（1）图像预处理，采用多尺度三维图割方法将视网膜内部自动分割成11个表面，提取位于第7表面和第8表面之间的内节/外节区域作为感兴趣区域，进行平坦化处理和双线性滤波增强；（2）对感兴趣区域中每个体素提取5类（共计57个）特征；（3）采用主成份分析方法对特征进行优化选择；（4）将特征样本分成训练集和测试集，对训练集样本采用ADABOOST算法训练集成分类器；（5）对测试集样本进行缺失/非缺失识别；（6）对识别结果进行血管轮廓影响排除和孤立点排除等后处理，计算相应的缺失体积，缺失体积识别误差小，准确性好。51INTCL权利要求书3页说明书13页附图4页19中华人民共和国国家知识产权。
3、局12发明专利申请权利要求书3页说明书13页附图4页10申请公布号CN104143087ACN104143087A1/3页21一种基于SDOCT的挫伤性视网膜内外节缺失检测方法，其特征在于包括以下步骤1图像预处理采用多尺度三维图割方法进行视网膜内部分层，将视网膜SDOCT图像自动分割成10个分层，产生11个表面，将包含内节/外节区域的第7层和第8层之间的数据提取出来，参考第11层进行平坦化，将平坦化后第7层和第8层之间的内节/外节区域的SDOCT图像作为感兴趣区域，进行双线性滤波增强；2体素特征提取对步骤1所述平坦化后的内节/外节区域的所有体素的灰度值进行归一化处理，对内节/外节区域内的缺失体。
4、素和非缺失体素分别提取5类特征，所述5类特征包括归一化的灰度值、分块均值、分块标准差、13个方向的灰度差绝对值和基于灰度共生矩阵的特征；3特征优化选择对步骤2中所提取的特征采用主成份分析法进行特征选择；4集成分类器的训练将经步骤3选择的特征分成训练集和测试集，在训练集中，按缺失样本数量与非缺失样本数量的比例为11的比例从非缺失样本中随机抽取数量与缺失样本数量相同的特征样本，采用ABABOOST算法进行集成分类器的训练；5体素缺失/非缺失的分类识别对测试集样本进行缺失/非缺失识别，采用步骤4训练得到的集成分类器对测试集样本进行分类，对每个体素的缺失与否进行分类识别；6后处理与缺失体积计算对步骤5。
5、分类识别得到的结果进行血管轮廓影响排除和孤立点排除，并计算视网膜内节/外节的缺失体积。2根据权利要求1所述的基于SDOCT的挫伤性视网膜内外节缺失检测方法，其特征在于所述步骤1具体包括以下步骤11视网膜内部分层采用多尺度三维图搜索方法将视网膜SDOCT图像自动分割成10个分层，产生11个表面；12内节/外节区域提取和平坦化将包含内节/外节区域的第7层和第8层之间的数据提取出来，参考第11层进行平坦化，将平坦化后的第7层与第8层之间的SDOCT图像作为感兴趣区域，进行内节/外节的缺失检测。3根据权利要求1所述的基于SDOCT的挫伤性视网膜内外节缺失检测方法，其特征在于所述步骤2所述5类特征的提取。
6、具体包括以下步骤21归一化灰度值表示每个体素的灰度级，对所述感兴趣区域内的所有体素的灰度值采用如式1所示的方法进行线性归一化至0255；其中，I,J,K表示体素的坐标，INORMALIZEDI,J,K表示体素I,J,K归一化后的灰度值，IORIGINALI,J,K表示体素I,J,K原始的灰度值，IMIN和IMAX分别表示感兴趣区域内体素灰度值的最小值和最大值；所述内节/外节缺失区域的灰度值低于非缺失区域的灰度值；22分块均值用于表示以体素I,J,K为中心的分块的灰度均值，所述分块为555分块，所述555分块的灰度均值MBLOCKI,J,K，计算公式如式2权利要求书CN104143087A2/3。
7、页323分块标准差STDBLOCKI,J,K用于表示以体素I,J,K为中心的555分块的灰度值变化情况，表征灰度值在555分块中的分散程度；根据分块标准差STDBLOCKI,J,K识别分块体素全部位于缺失区域或者分块体素全部位于非缺失区域或者分块中同时包含缺失区域和非缺失区域的体素，计算公式如下2413个方向上的灰度差绝对值之和定义为中心点体素I,J,K的灰度值与13个方向周围体素的灰度差绝对值之和，反映中心点体素在缺失区域或非缺失区域的位置，依据所述13个方向上的灰度差绝对值之和AIDSTEPI,J,K判断中心点体素I,J,K位于缺失与非缺失交界处或者缺失区域或者非缺失区域，选择两种距离步长。
8、计算13个方向上的灰度差绝对值之和AIDSTEPI,J,K，两种距离步长为STEP1，2，具体计算公式如公式4表示13个方向上的灰度差绝对值之和；所述13个方向的选取方法为假设1表示X轴与XY平面投影图的夹角，2表示Z轴与XY平面投影图的夹角，那么所述13个方向分别为1,20,90、45,90、90,90、135,90、0,45、180,45、90,45、90,45、0,0、45,45、135,45、45,45、135,45；25基于灰度共生矩阵特征描述若干立体切面之间灰度级的空间相关性，所述灰度共生矩阵为三维体数据的灰度共生矩阵，采用三维方法在13个方向上跨切面搜索不同灰度级，构造13个灰度。
9、共生矩阵，13个灰度共生矩阵构造方法具体包括以下步骤对每个555分块构造所述13个不同方向上的灰度共生矩阵；分别提取4个基于所述灰度共生矩阵的特征，所述4个基于灰度共生矩阵的特征为A对比度、B相关性、C能量和D均匀性；对比度CONTRASTDIRECTIONMI,J,K、相关性CORRELATIONDIRECTIONMI,J,K、能量ENERGYDIRECTIONMI,J,K和均匀性HOMOGENEITYDIRECTIONMI,J,K分别采用式5、式6、式7和式8进行计算权利要求书CN104143087A3/3页4其中，PX,Y表示灰度共生矩阵中第X,Y个元素；NG表示量化后的灰度级，取NG8。
10、；X和X分别表示第X行的均值和标准差；Y和Y分别表示第Y列的均值和标准差。4根据权利要求1所述的基于SDOCT的挫伤性视网膜内外节缺失检测方法，其特征在于，所述步骤3具体包括以下步骤对于步骤2所述提取的5类特征采用主成份分析方法进行特征选择，将经过特征选择后的结果作为新的特征，根据新特征所代表的数据总方差的比例大于或者等于90的原则，选择前12个主成分作为新特征；完成特征优化选择后，感兴趣区域内的每个体素用一个12维的特征向量表征，在特征空间中进行分类。5根据权利要求1所述的基于SDOCT的挫伤性视网膜内外节缺失检测方法，其特征在于，所述步骤4具体包括以下步骤将经步骤3选择的特征分成训练集和测。
11、试集，在训练集中，通过下采样方式以缺失样本数量与非缺失样本数量的比例为11构造新的训练集，采用分类器集成的方法进行分类，以决策树为弱分类器，采用ADABOOST算法进行集成分类器的训练。6根据权利要求1所述的基于SDOCT的挫伤性视网膜内外节缺失检测方法，其特征在于，所述步骤6具体包括以下步骤61血管轮廓影响的排除将第1至第10层的全部三维SDOCT数据投影至XY平面，对投影图进行血管轮廓检测，将步骤5的缺失识别结果的X坐标值等于血管轮廓的X坐标值，则该X坐标值处的体素被认为是缺失误识别，将对应的识别结果修改为非缺失区域；62将步骤61的缺失识别结果采用形态学方法进行孤立点的排除；63将步骤6。
12、2的孤立点的排除后的缺失识别结果乘以体素的分辨率，转化为检测到的内节/外节缺失体积。权利要求书CN104143087A1/13页5基于SDOCT的挫伤性视网膜内外节缺失检测方法技术领域0001本发明属于视网膜图像处理与分析方法，尤其是对SDOCT频域光学相干断层成像视网膜图像中的内节/外节缺失的三维自动定量检测方法。背景技术0002眼外伤是导致视力损伤甚至致盲的主要因素之一，是我国单眼盲的首位致盲原因。当眼球表面的损伤传递到后端的视网膜时，视网膜震荡伤的特点是灰白变色，或是当外伤闭合后，视网膜仍然浑浊。人眼的组织病理学研究和动物研究表明，光感受器的内节/外节连接区损伤是视网膜震荡伤的发病原。而。
13、光感受器的内节/外节与光传导紧密相关，因此，内节/外节的完整性与视力密切相关。0003SDOCT频域光学相干断层成像是一种无损评估视网膜生理和病理指标的强有力技术，能快速获得高分辨的断层三维图像。在SDOCT图像中，光感受器内节/外节是介于外界膜和色素上皮层之间的一个高亮反射层。内节/外节完整性的遭破坏程度是光感受器受损程度的一个代表性指标，通常与视网膜震荡伤导致的视力低下以及其他视网膜疾病有关。0004在SDOCT图像中，通过定量检测内节/外节区域缺失范围来定量评估光感受器的损伤程度是一个非常有意义的研究方向。有文献报道在一个B扫描面中进行内节/外节缺失长度的手动检测。但是，这是基于单个二维。
14、界面图像。此外，有文献提出采用部分OCT投影图像或者冠状面投影图像等技术更好地可视化光感受器的整体性和缺失区域。但是，这些检测内节/外节缺失面积的技术仍然是基于二维图像的，而且是手动的，这在选择缺失区域边界时可能会引入主观因素。关于内节/外节缺失区域边界的自动检测和内节/外籍缺失面积计算方法也有文献报道，但是，这种方法是基于缺失区域是圆形的假设，而实际上，缺失区域的形状可能是任意形状。0005更重要的是，内节/外节是一层有厚度的薄膜，应该以内节/外节缺失的体积来衡量光感受器的损伤程度。这方面的研究还处于起步阶段，现有技术对视网膜内外节缺失检测也仅仅只针对局部较小的范围进行检测，特征提取不全面，。
15、并且分类不平衡，例如，CN103679198A公开了一种基于K近邻的视网膜内节外节缺失自动检测，提取以黄斑中心为中心、直径1MM的内节外节区域为感兴趣区域，进行5个特征提取，提取特征少，特征不够典型，且由于K近邻分类器的约束，采用的K近邻分类存在分类不平衡问题，分类性能差；其次，在后处理与缺失体积计算中，对识别分类得到的结果没有进行血管轮廓影响排除和孤立点排除，缺失体积误差大，准确性差。发明内容0006本发明克服了现有技术的不足，提供了一种具有可行性和有效性的视网膜内节/外节缺失三维自动定量检测方法，能够将内节/外节全部区域作为感兴趣区域，全面提取体素特征，充分考虑不平衡分类问题的特点，精确计。
16、算视网膜内外节缺失区域。说明书CN104143087A2/13页60007为解决上述技术问题，本发明采用的技术方案为0008一种基于SDOCT的挫伤性视网膜内外节缺失检测方法，该方法主要包括6个步骤图像预处理，体素特征提取，特征选择，集成分类器的训练，体素缺失识别，后处理与缺失体积计算。0009一种基于SDOCT的挫伤性视网膜内外节缺失检测方法，即包括以下步骤00101图像预处理采用多尺度三维图割方法进行视网膜内部分层，将视网膜SDOCT图像自动分割成10个分层，产生11个表面，将包含内节/外节区域的第7层和第8层之间的数据提取出来，参考第11层进行平坦化，将平坦化后第7层和第8层之间的内节/。
17、外节区域的SDOCT图像作为感兴趣区域，进行双线性滤波增强；00112体素特征提取对步骤1所述平坦化后的内节/外节区域的所有体素的灰度值进行归一化处理，根据手动标记的金标准，对内节/外节区域内的缺失体素和非缺失体素分别提取5类特征，所述5类特征包括归一化的灰度值、分块均值、分块标准差、13个方向的灰度差绝对值、基于灰度共生矩阵的特征；00123特征优化选择对步骤2中所提取的特征采用主成份分析法进行特征选择；00134集成分类器的训练将经步骤3选择的特征分成训练集和测试集；在训练集中，按11比例，从非缺失样本中随机抽取数量与缺失样本数量相同的特征样本，采用ABABOOST算法进行集成分类器的训练。
18、；00145体素缺失/非缺失的分类识别对测试集样本进行缺失/非缺失识别，采用步骤4训练得到的集成分类器对测试样本进行分类，对每个体素的缺失与否进行分类识别；00156后处理与缺失体积计算对步骤5分类识别得到的缺失结果进行血管轮廓影响排除和孤立点排除，并计算相应的缺失体积。0016本发明技术方案具体描述如下。00171图像预处理0018图像预处理主要包括以下三个步骤视网膜内部分层，包含内节/外节区域的感兴趣区域提取和平坦化，以及图像的双线性滤波增强。001911视网膜内部分层0020人类视网膜是一个相当薄的组织，厚度不到05MM，是眼睛最重要的组成部分。视网膜本身具有相当复杂的结构，基本结构类似。
19、于一个“三层蛋糕”三层神经细胞被两层类神经连接层隔开，可以更精细地分成10个层。视网膜内部分层对于分析视网膜病变如眼外伤的严重程度、黄斑水肿的形成等具有重要作用。采用多尺度三维图搜索方法将视网膜SDOCT图像自动分割成10个分层，产生11个表面。所述10个分层自上至下依次为神经纤维层，神经节细胞层，内丛状层，内核层，外丛状层，外核层，内节层，外节层，维尔赫夫氏膜和视网膜色素上皮层；0021三维图搜索方法是基于图论的方法，采用从粗糙到精细的不同分辨率，来检测视网膜的内部各表面。图搜索采用基于边界的代价函数，当代价函数最小时，即找到各个表面；002212内节/外节区域提取和平坦化说明书CN1041。
20、43087A3/13页70023完成了各分层的分割后，可以将包含内节/外节区域的第7层和第8层之间的数据提取出来，并参考视网膜色素上皮层的下界即第11层进行平坦化，因为第11层的分层效果相对鲁棒。本发明主要平坦化后的第7层与第8层之间的SDOCT图像作为感兴趣区域，研究内节/外节的缺失检测。即将包含内节/外节区域的第7层和第8层之间的数据提取出来，参考第11层进行平坦化，将平坦化后第7层和第8层之间的内节/外节区域的SDOCT图像作为感兴趣区域，进行双线性滤波增强。00242特征提取0025本发明是一种基于纹理特征分类的内节/外节缺失检测方法，首先，对感兴趣区域内的所有体素的灰度值进行归一化处。
21、理，然后提取了5类特征，5类特征包括归一化的灰度值、分块均值、分块标准差、13个方向的灰度差绝对值、基于灰度共生矩阵的特征；5类特征共计57个特征归一化的灰度值特征1，分块均值特征2，分块标准差特征3，13个方向的灰度差绝对值步长分别取1和2，特征45，基于灰度共生矩阵的特征，基于灰度共生矩阵的特征包括对比度13个方向，特征618、相关性13个方向，特征1931、能量13个方向，特征3244、均匀性13个方向，特征4557。002621归一化灰度值表示每个体素的灰度级，对感兴趣区域内的所有体素的灰度值采用如式1所示的方法进行线性归一化至0255。00270028其中，I,J,K表示体素的坐标，。
22、INORMALIZEDI,J,K表示体素I,J,K归一化后的灰度值，IORIGINALI,J,K表示体素I,J,K原始的灰度值，IMIN和IMAX分别表示感兴趣区域内体素灰度值的最小值和最大值。内节/外节缺失区域的一个重要特点是灰度值低于非缺失区域，即内节/外节缺失区域的灰度值低于非缺失区域的灰度值。002922分块均值用于表示以体素I,J,K为中心的分块的灰度均值，555分块的灰度均值MBLOCKI,J,K，计算公式如式200300031由于内节/外节缺失区域的灰度值整体比非缺失区域的灰度值小，所以内节/外节缺失区域的分块均值MBLOCKI,J,K小于非缺失区域的分块均值。003223分块标。
23、准差STDBLOCKI,J,K用于表示以体素I,J,K为中心的555分块的灰度值变化情况，表征灰度值在555分块中的分散程度；可根据分块标准差STDBLOCKI,J,K识别分块体素全部位于缺失区域、分块体素全部位于非缺失区域、或者分块中同时包含缺失区域和非缺失区域的体素，计算公式如公式200330034依据分块标准差STDBLOCKI,J,K的值识别分块体素位于缺失区域或者非缺失区域或者同时包含缺失区域和非缺失区域的体素；如果某分块中的体素灰度值都很接近，说明书CN104143087A4/13页8即该分块全部位于缺失区域或者全部位于非缺失区域，即STDBLOCKI,J,K的值属于一个经验值偏小。
24、的范围，STDBLOCKI,J,K值的范围与图像本身体素灰度值有关，例如小于所有体素灰度值STDBLOCKI,J,K平均值；如果分块中同时包含缺失区域和非缺失区域的体素，则对应的STDBLOCKI,J,K的值属于一个经验值偏大的范围，例如，大于所有体素灰度值STDBLOCKI,J,K平均值。通过STDBLOCKI,J,K的值识别分块体素位于缺失区域或者非缺失区域或者同时包含缺失区域和非缺失区域的体素，STDBLOCKI,J,K的值识别范围没有一个明确的范围，基于操作者的经验值识别，一般如果某个分块的STDBLOCKI,J,K小于所有体素灰度值STDBLOCKI,J,K平均值，分块全部位于缺失区。
25、域或者全部位于非缺失区域，否则，分块中同时包含缺失区域和非缺失区域的体素；00352413个方向上的灰度差绝对值之和定义为中心点体素I,J,K的灰度值与周围体素13个方向的灰度差绝对值之和，13个方向上的灰度差绝对值反映中心点体素在缺失区域或非缺失区域的位置，依据所述13个方向上的灰度差绝对值之和AIDSTEPI,J,K判断中心点体素I,J,K位于缺失与非缺失交界处或者缺失区域或者非缺失区域；若该特征值属于一个比较大的范围，一般大于或者等于05，说明该中心点体素位于缺失与非缺失交界处；若该特征值接近于零，一般小于05，说明中心点体素位于缺失区域或者非缺失区域，这里选择了两种距离步长，STEP1。
26、，2，即STEP为1或者2，具体计算公式如下00360037其中所述13个方向的选取方法为假设1表示X轴与XY平面投影图的夹角，2表示Z轴与XY平面投影图的夹角，那么所述13个方向分别为1,20,90、45,90、90,90、135,90、0,45、180,45、90,45、90,45、0,0、45,45、135,45、45,45、135,45；13个方向对应的L、M、N取值分别为ISTEP,J,K、ISTEP,JSTEP,K、I,JSTEP,K、ISTEP,JSTEP,K、ISTEP,J,KSTEP、ISTEP,J,KSTEP、I,JSTEP,KSTEP、I,JSTEP,KSTEP、I,J,。
27、KSTEP、ISTEP,JSTEP,KSTEP、ISTEP,JSTEP,KSTEP、ISTEP,JSTEP,KSTEP、ISTEP,JSTEP,KSTEP；表示13个方向上灰度差绝对值之和。003825基于灰度共生矩阵的特征描述多个立体切面之间灰度级的空间相关性，基于灰度共生矩阵的特征为三维体数据的灰度共生矩阵，采用了一种三维方法在13个方向上跨切面搜索不同灰度级，可以构造13个灰度共生矩阵；13个灰度共生矩阵构造方法具体包括以下步骤0039对每个555分块构造所述13个不同方向步骤24中选取的13个方向上的灰度共生矩阵；然后，分别提取4个基于该灰度共生矩阵的特征A对比度，用于度量立体图像的局。
28、部对比度，当大灰度差情况发生概率大时，对比度值较大；B相关性，用于表征体素对中两个体素的相关性，当体素对的灰度值越相关时，该相关性特征值越大；C能量，用于度量体素对的重复出现次数，当体素对重复次数越多时，能量值越大；D均匀性，用于表征体素对的局部均匀性，当每个体素对的灰度级越相近时，该均匀性特征值越大。说明书CN104143087A5/13页90040对比度CONTRASTDIRECTIONMI,J,K、相关性CORRELATIONDIRECTIONMI,J,K、能量ENERGYDIRECTIONMI,J,K和均匀性HOMOGENEITYDIRECTIONMI,J,K分别采用式5、式6、式7和。
29、式8进行计算00410042004300440045其中，PX,Y表示灰度共生矩阵中第X,Y个元素；NG表示量化后的灰度级一般基于经验值选取量化后的主要灰度级，取NG8；X和X分别表示第X行的均值和标准差；Y和Y分别表示第Y列的均值和标准差。00463特征优化选择0047特征选择目的是降低特征空间的维数，从大的特征集中选择主要的特征，对于训练样本，本发明采用主成份分析方法进行特征选择。对于步骤2所述提取的5类特征采用主成份分析方法进行特征选择，将经过特征选择后的结果作为新的特征，根据新特征所代表的数据总方差的比例大于等于90的原则，选择前12个主成分作为新特征；0048完成特征优化选择后，感兴。
30、趣区域内的每个体素用一个12维的特征向量来表征，步骤4的集成分类器将在特征空间中进行。00494集成分类器的训练0050由于缺失体素的数量远少于非缺失体素的数量，两类样本的类间分布不平衡，属于典型的非平衡分类问题。传统的分类器，大都建立在训练样本集分布基本均衡的前提下，且以总体分类精度作为目标，在本问题中，由于缺失类样本相对稀缺以及不可避免的噪声问题，传统的分类器如支持向量机、贝叶斯分类器、分类与回归树方法、决策树、线性判别分析、双线性判别分析等，都会出现分类器对数量上占优势的多数类非缺失类存在明显的偏好，使得多数类的分类精度提高，而少数类缺失类的分类精度降低。本发明通过下采样方式舍弃大部分多。
31、数类样本，以缺失样本数量与非缺失样本数量的比例为11比例构造新的训练集，并采用分类器集成的方法来进一步克服不平衡分类问题，以决策树为弱分类器，采用ADABOOST算法进行集成分类器的训练。00515体素缺失/非缺失的分类识别0052采用步骤4训练得到的集成分类器，对内节/外节区域内的体素进行分类识别，分成缺失/非缺失两类。说明书CN104143087A6/13页1000536后处理与缺失体积计算0054由于视网膜血管轮廓的影响，SDOCT图像的B扫描图的相应位置上会出现上下贯穿的阴影，在IS/OS连接区内也会出现血管阴影，这些阴影往往会被误判为缺失，因此需要进行排除。本设计首先将三维OCT数据。
32、投影至XY平面，并对投影图进行血管轮廓检测，若步骤5得到的缺失识别结果的X坐标值与血管轮廓的X坐标值相等，则认为是缺失误识别，将对应的识别结果改为非缺失，可有效地排除血管轮廓对内节/外节缺失检测的影响。0055由于内节/外节缺失的生理特点，孤立的缺失点或者孤立的非缺失点应该不存在，因此，对检测结果采用形态学方法进行了孤立点的排除。0056经后处理的检测结果是以体素个数为单位的，最后检测结果乘以体素的分辨率，转化为检测到的内节/外节缺失体积。0057步骤6具体包括以下步骤005861血管轮廓影响的排除将第1至第10层的全部三维SDOCT数据投影至XY平面，对投影图进行血管轮廓检测，将步骤5的缺失。
33、识别结果的X坐标值等于血管轮廓的X坐标值，则认为是缺失误识别，将对应的识别结果修改为非缺失区域；005962将步骤61的缺失识别结果采用形态学方法进行孤立点的排除；006063将步骤62的孤立点的排除后缺失识别结果乘以体素的分辨率，转化为检测到的内节/外节缺失体积。0061本发明有益效果包括0062本发明提供了一种基于SDOCT图像、具有可行性和有效性的视网膜内节/外节缺失三维自动检测方法，将感兴趣区域扩展至整个内节/外节连接区；提取57个体素特征，采用了主成份分析方法进行了特征选择与优化，充分考虑到了不平衡分类问题的存在，采用了ADABOOST算法训练集成分类器，从测试结果来看，分类性能比之。
34、前申请有较大提高，视网膜内节外节检测准确，同时，本发明对识别分类得到的结果进行血管投影和孤立点排除，缺失体积识别误差小，准确性好；0063进一步地，本申请扩大了训练和测试的数据集，包括15个外伤眼和15个正常眼；对于眼外伤等临床常见眼科疾病的诊断与治疗起到了重要的辅助作用。附图说明0064图1为本发明的内节/外节缺失三维自动检测的整体框图；0065图2为图2视网膜的分层效果图SDOCT图像的XZ平面图，分成了10个层11个表面，自上至下依次为神经纤维层，神经节细胞层，内丛状层，内核层，外丛状层，外核层，内节层，外节层，维尔赫夫氏膜和视网膜色素上皮层；0066图3为平坦化后的SDOCT图像和感兴。
35、趣区域XZ平面图；0067图4为双线性滤波后的SDOCT图像和感兴趣区域XZ平面图；0068图5为血管轮廓的影响，A箭头所示为血管轮廓在OCT图像XZ平面图中的位置；B视网膜OCT图像在XY平面上的投影图以及血管；A和B图中箭头处的血管位置一一对应；0069图6内节/外节缺失检测结果与金标准。图A表示由观察者1手动标记的一个外伤眼数据的内节/外节缺失三维金标准，图B表示对该外伤眼数据采用本发明方法检说明书CN104143087A107/13页11测到的内节/外节缺失区域；图C表示采用本发明方法对正常眼数据的检测结果；0070图7为外伤数据和正常数据缺失体积检测结果的比较。具体实施方式0071下。
36、面结合附图对本发明作更进一步的说明。0072本方法的基本框图如附图1所示，主要包括6个步骤图像预处理，体素特征提取，特征优化选择，集成分类器训练，体素缺失/非缺失分类识别，后处理与缺失体积计算。具体描述如下。00731图像预处理0074图像预处理主要包括以下三个步骤视网膜内部分层，包含内节/外节区域的感兴趣区域的提取和平坦化，图像的双线性滤波增强。0075A视网膜内部分层0076人类视网膜是一个相当薄的组织，厚度不到05MM，是眼睛最重要的组成部分。视网膜本身具有相当复杂的结构，基本结构类似于一个“三层蛋糕”三层神经细胞被两层类神经连接层隔开，可以更精细地分成10个层。视网膜内部分层对于分析视。
37、网膜病变如眼外伤的严重程度、黄斑水肿的形成等具有重要作用。采用多尺度三维图搜索方法将视网膜SDOCT图像自动分割成10个分层，产生11个表面，如附图2所示。三维图搜索方法是基于图论的方法，采用从粗糙到精细的不同分辨率，来检测视网膜的内部各表面。图搜索采用基于边界的代价函数，当代价函数最小时，即找到各个表面。0077B内节/外节区域提取和平坦化0078完成了各分层的分割后，可以将包含内节/外节区域的第7层和第8层之间的数据提取出来，并参考视网膜色素上皮层第10层的下界即第11层进行平坦化，因为第11层的分层效果相对鲁棒。本发明主要平坦化后的第7层与第8层之间的SDOCT图像作为感兴趣区域，研究内。
38、节/外节的缺失检测。平坦化后的SDOCT图像和感兴趣区域XZ平面图如附图3所示，双线性滤波后的SDOCT图像和感兴趣区域XZ平面图如图4所示，滤除图像中的散斑噪音。00792体素特征提取0080本发明是一种基于纹理特征分类的内节/外节缺失检测方法。首先，对感兴趣区域内的所有体素的灰度值进行归一化处理，然后提取了5类特征，共计57个特征归一化的灰度值特征1，分块均值特征2，分块标准差特征3，13个方向的灰度差绝对值步长分别取1和2，特征45，基于灰度共生矩阵的特征，包括对比度13个方向，特征618，相关性13个方向，特征1931，能量13个方向，特征3244，均匀性13个方向，特征4557。00。
39、81归一化灰度值表示每个体素的灰度级，对感兴趣区域内的所有体素的灰度值采用如式1所示的方法进行线性归一化至0255。00820083其中，I,J,K表示体素的坐标，INORMALIZEDI,J,K表示体素I,J,K归一化后的灰度值，IORIGINALI,J,K表示体素I,J,K原始的灰度值，IMIN和IMAX分别表示感兴趣区域内说明书CN104143087A118/13页12体素灰度值的最小值和最大值。内节/外节缺失区域的一个重要特点是灰度值低于非缺失区域。0084分块均值用于表示以体素I,J,K为中心的小分块的灰度均值，这里计算了555小分块的灰度均值MBLOCKI,J,K，计算公式如下00。
40、850086由于内节/外节缺失区域的灰度值整体比非缺失区域的灰度值小，所以内节/外节缺失区域的分块均值MBLOCKI,J,K小于非缺失区域的分块均值。0087分块标准差STDBLOCKI,J,K用于表示以体素I,J,K为中心的555小分块的灰度值变化情况，可以表征灰度值在555分块中的分散程度。可根据分块标准差识别分块体素全部位于缺失区域、分块体素全部位于非缺失区域、或者分块中同时包含缺失区域和非缺失区域的体素，计算公式如下0088008913个方向上的灰度差绝对值之和AIDSTEPI,J,K定义为中心点体素I,J,K的灰度值与周围体素13个方向的灰度差绝对值之和，该特征可以很好地反映中心点体。
41、素在缺失区域或非缺失区域的大概位置，若该特征值大于05，说明该中心点体素位于缺失与非缺失交界处；若该特征值小于05，说明中心点体素位于较大块缺失区域或者非缺失区域，这里选择了两种距离步长，STEP1，2，具体计算公式如下00900091其中所述13个方向的选取方法为假设1表示X轴与XY平面投影图的夹角，2表示Z轴与XY平面投影图的夹角，那么所述13个方向分别为1,20,90、45,90、90,90、135,90、0,45、180,45、90,45、90,45、0,0、45,45、135,45、45,45、135,45；13个方向对应的L、M、N取值分别为ISTEP,J,K、ISTEP,JSTE。
42、P,K、I,JSTEP,K、ISTEP,JSTEP,K、ISTEP,J,KSTEP、ISTEP,J,KSTEP、I,JSTEP,KSTEP、I,JSTEP,KSTEP、I,J,KSTEP、ISTEP,JSTEP,KSTEP、ISTEP,JSTEP,KSTEP、ISTEP,JSTEP,KSTEP、ISTEP,JSTEP,KSTEP；表示13个方向上灰度差绝对值之和。0092三维体数据的灰度共生矩阵可以描述多个立体切面之间灰度级的空间相关性。本设计采用了一种三维方法在13个方向上跨切面搜索不同灰度级，可以构造13个灰度共生矩阵。假设1表示X轴与XY平面投影图的夹角，2表示Z轴与XY平面投影图的夹角。
43、，那么13个方向分别为1,20,90、45,90、90,90、135,90、0,45、180,45、90,45、90,45、0,0、45,45、135,45、45,45、135,45。这里，对每个555小分块构造了以上13个不同方向上的灰度共生矩阵。然后，分别提取4个基于该灰度共生矩阵的特征A对说明书CN104143087A129/13页13比度，用于度量立体图像的局部对比度，当大灰度差情况发生概率大时，对比度值较大；B相关性，用于表征体素对中两个体素的相关性，当体素对的灰度值越相关时，该相关性特征值越大；C能量，用于度量体素对的重复出现次数，当体素对重复次数越多时，能量值越大；D均匀性，用于。
44、表征体素对的局部均匀性，当每个体素对的灰度级越相近时，该均匀性特征值越大。对比度CONTRASTDIRECTIONMI,J,K、相关性CORRELATIONDIRECTIONMI,J,K、能量ENERGYDIRECTIONMI,J,K和均匀性HOMOGENEITYDIRECTIONMI,J,K分别采用式5、式6、式7和式8进行计算。00930094009500960097其中，PX,Y表示灰度共生矩阵中第X,Y个元素；NG表示量化后的灰度级，取NG8；X和X分别表示第X行的均值和标准差；Y和Y分别表示第Y列的均值和标准差。00983特征优化选择0099特征选择是模式识别中常用的一种技术，目的是。
45、降低特征空间的维数，从大的特征集中选择主要的特征。本发明采用了主成份分析方法进行特征选择。根据新特征所代表的数据总方差的比例大于等于90，选择前12个主成分作为新特征。0100完成特征优化选择后，感兴趣区域内的每个体素可用一个12维的特征向量来表征，接下来的分类将在特征空间中进行。01014集成分类器的训练0102由于缺失体素的数量远少于非缺失体素的数量，两类样本的类间分布不平衡，属于典型的非平衡分类问题。传统的分类器，大都建立在训练样本集分布基本均衡的前提下，且以总体分类精度作为目标，在本问题中，由于缺失类样本相对稀缺以及不可避免的噪声问题，传统的分类器如支持向量机、贝叶斯分类器、分类与回归。
46、树方法、决策树、线性判别分析、双线性判别分析等，都会出现分类器对数量上占优势的多数类非缺失类存在明显的偏好，使得多数类的分类精度提高，而少数类缺失类的分类精度降低。本发明通过下采样方式舍弃大部分多数类样本，以缺失样本数量与非缺失样本数量11的比例构造训练集，并采用分类器集成的方法来进一步克服不平衡分类问题，以决策树为弱分类器，采用ADABOOST算法进行集成分类器的训练。01035体素缺失/非缺失的分类识别0104采用上述训练得到的集成分类器，对内节/外节区域内的体素进行分类识别，分说明书CN104143087A1310/13页14成缺失/非缺失两类。01056后处理与缺失体积计算0106如附。
47、图5A所示，箭头所示为血管轮廓在OCT图像XZ平面图中的位置，由于视网膜血管轮廓的影响，SDOCT图像的B扫描图的相应位置上会出现上下贯穿的阴影，在IS/OS连接区内也会出现血管阴影，这些阴影往往会被误判为缺失，因此需要进行排除。首先将第1至第10层的全部三维OCT数据投影至XY平面，如图5B所示，视网膜OCT图像在XY平面上的投影图以及血管，图5A和图5B在图中箭头处的血管位置一一对应，本申请对投影图进行血管轮廓检测，如步骤5中的缺失识别结果的X坐标值与血管轮廓的X坐标值相等，则认为是缺失误检测，将相应的识别结果改为为非缺失，可有效地排除血管轮廓对内节/外节缺失检测的影响。0107由于内节/。
48、外节缺失的生理特点，孤立的缺失点或者孤立的非缺失点应该不存在，因此，对缺失识别结果进一步采用形态学方法进行了孤立点的排除。0108经后处理的缺失识别结果是以体素个数为单位的，最后缺失识别结果乘以体素的分辨率，转化为检测到的内节/外节缺失体积。01097实验方法与结果0110内节/外节缺失的金标准由2个观察者在眼科医生指导下独立手动标注。对15个存在内节/外节缺失的眼外伤病人和15个正常人数据进行了本发明方法的测试。采用了留一交叉验证方法来检验本发明方法的可行性和有效性。由于在训练阶段，对非缺失样本进行了随机降采样，为了减少随机误差的影响，采用ADABOOST算法训练了5个集成分类器模型，并对1。
49、5个外伤数据和15个正常数据的内节/外节区域分别进行了5次缺失检测。0111对外伤病人和正常人的内节/外节缺失检测体积的均值及其95置信区间进行了比较，对采用T检验方法评估两组数据缺失体积的统计差异性。0112为了进一步评估本发明方法的性能，将检测结果与第1个观察者手动标记的金标准进行比较。采用灵敏度SENSITIVITY，SEN、特异度SPECICITY，SPE和均衡正确检测率作为评估的客观指标，分别采用式9、式10和式11计算0113011401150116其中，TP、TN、FP和FN分别表示真阳、真阴、假阳和假阴。0117附图6A给出了由观察者1手动标记的一个外伤眼数据的内节/外节缺失三维金标准结果附图中心的阴影部分即为内节/外节缺失区域，附图6B给出了对该外伤眼数据采用本发明方法进行内节/外节缺失检测的结果，附图6C给出了本发明方法对一个正常眼数据进行检测的结果。附图6A理解为挫伤性视网膜内节/外节缺失区域期望的检测结果，附图6B采用本发明技术方案检测的结构，对比附图6A的期望的检测结果，本发明技术方案检测的结。