申请/专利权人：清华大学

申请日：2018-05-31

公开（公告）日：2020-11-20

公开（公告）号：CN108846834B

主分类号：G06T7/00(20170101)

分类号：G06T7/00(20170101);G06T7/11(20170101);G06T7/136(20170101);G06T7/40(20170101);G06T7/62(20170101)

优先权：

专利状态码：有效-授权

法律状态：2020.11.20#授权;2018.12.14#实质审查的生效;2018.11.20#公开

摘要：本发明提供一种医用图像处理装置及医用图像处理方法，包括：图像存储部，存储包含股骨头的第一图像；图像分割部，分割出仅含有股骨头部分的第二图像；特征检验部，检验第二图像边缘是否为圆，显示不为圆部分的图像；纹理检测部，包括：ROI区域选取单元，从第二图像选取ROI区域；灰度变换单元，所述区域作灰度梯度变换，大于第一阈值的点置为0；细分单元，将ROI区域细分成多个第二区域，ROI区域面积是第二区域面积的整数倍；计数单元，对第二区域内的0的个数进行计数；第一判断单元，判断第二区域0个数是否超过第二阈值；第一显示单元，显示超过第二阈值的第二区域的图像。上述装置和方法能够实现股骨头医用图像的自动化处理。

主权项：1.一种医用图像处理装置，其特征在于，包括：图像存储部，用于存储包含股骨头的第一图像;图像分割部，分割出仅含有股骨头部分的第二图像；特征检验部，检验第二图像的边缘是否为圆，显示不为圆的部分的图像；纹理检测部，包括：ROI区域选取单元，在第二图像上选取ROI区域；灰度变换单元，将ROI区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为0；细分单元，将ROI区域细分成多个第二区域，ROI区域面积是第二区域面积的整数倍；计数单元，对每个第二区域内的0的个数进行计数；第一判断单元，判断第二区域的0的个数是否超过第二阈值；第一显示单元，将异常部分在第二图像或第一图像内部呈现出来，所述异常部分是超过第二阈值的第二区域，其中，所述特征检验部包括：边缘提取单元，提取第二图像中股骨头边缘部分的点；拟合单元，采用拟合方法拟合上述各点对应的圆心和半径；确定单元，确定各点到拟合的圆心的距离以及所述距离与拟合的半径的差值的绝对值；第二判断单元，判断上述差值的绝对值是否大于第三阈值；第二显示单元，显示所述绝对值大于第三阈值的点，将大于第三阈值的绝对值对应的点作为股骨头异常的部分；加法器，将各点的所述绝对值相加；第三判断单元，判断上述相加的结果是否大于第四阈值，如果大于第四阈值，发出预警信号。

全文数据：医用图像处理装置及医用图像处理方法技术领域[0001]本发明涉及医疗技术领域，更具体地，涉及一种医用图像处理装置及医用图像处理方法。背景技术[0002]通常，关于股骨头的医用图像依靠医生的肉眼观测，并没有对股骨头的医用图像进行自动化处理的装置和方法。发明内容[0003]鉴于上述问题，本发明提供一种对股骨头的医用图像进行自动化处理的医用图像处理装置及医用图像处理方法。[0004]根据本发明的一个方面，提供一种医用图像处理装置，包括：[0005]图像存储部，用于存储包含股骨头的第一图像；[0006]图像分割部，分割出仅含有股骨头部分的第二图像；[0007]特征检验部，检验第二图像的边缘是否为圆，显示不为圆的部分的图像；[0008]纹理检测部，包括:R0I区域选取单元，在第二图像上选取R0I区域;灰度变换单元，将R0I区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为〇;细分单元，将R0I区域细分成多个第二区域，R0I区域面积是第二区域面积的整数倍;计数单元，对每个第二区域内的〇的个数进行计数;第一判断单元，判断第二区域的0的个数是否超过第二阈值;第一显示单元，显示超过第二阈值的第二区域的图像。[0009]根据本发明的另一个方面，提供一种医用图像处理方法，包括：[0010]步骤S1，存储包含股骨头的第一图像；[0011]步骤S2,分割出仅含有股骨头部分的第二图像；[0012]步骤S3,检验第二图像的边缘是否为圆，显示不为圆的部分的图像；[0013]步骤S4,检测第二图像的纹理特征，包括:在第二图像上选取R0I区域;将R0I区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为0;将R0I区域细分成多个第二区域，R〇I区域的面积是第二区域面积的整数倍;对每个第二区域内的0的个数进行计数;判断第二区域的〇的个数是否超过第二阈值，显示超过第二阈值的第二区域的图像。[0014]本发明所述医用图像处理装置及医用图像处理方法对股骨头的医用图像进行分割和不同区域划分显示自动化处理，方便观看，一目了然。[0015]上述医用图像处理装置实现R0I区域分割和疾病诊断的一体化，其分割准确、速度较快，可以为病人提供可靠的诊断结果;能够对医生诊断病变起参考作用;降低基层医院医生的误诊率，减少医生工作量;能够识别和诊断潜在病变。’附图说明[0016]通过参考以下具体实施方式及权利要求书的内容并且结合附图，本发明的其它目的及结果将更加明白且易于理解。在附图中：[0017]图1是本发明所述医用图像处理装置构成框图的示意图；[0018]图2是本发明所述医用图像处理方法流程图的示意图。具体实施方式[0019]在下面的描述中，出于说明的目的，为了提供对一个或多个实施例的全面理解，阐述了许多具体细节。然而，很明显，也可以在没有这些具体细节的情况下实现这些实施例。在其它例子中，为了便于描述一个或多个实施例，公知的结构和设备以方框图的形式示出。[0020]下面将参照附图来对根据本发明的各个实施例进行详细描述。[0021]图1是本发明所述医用图像处理装置构成框图的示意图，如图1所示，所述医用图像处理装置包括：[0022]图像存储部1，用于存储包含股骨头的第一图像；[0023]图像分割部2,分割出仅含有股骨头部分的第二图像，例如，用局部二元拟合LBF的方法，将人体髋关节医用图像中的股骨头部分分割出来；[0024]特征检验部3,检验第二图像的边缘是否为圆，显示不为圆的部分的图像；[0025]纹理检测部4,包括：[0026]R0I区域选取单元41，在第二图像上选取R0I区域感兴趣区域，预进行纹理检测的区域，也可以直接将整个第二图像作为R0I区域）；[0027]灰度变换单元42,将R0I区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为〇;[0028]细分单元43,将R0I区域细分成多个第二区域，R0I区域面积是第二区域面积的整数倍；[0029]计数单元44，对每个第二区域内的0的个数进行计数；[0030]第一判断单元45,判断第二区域的0的个数是否超过第二阈值，可以作为判断第二区域是否病变的依据；[0031]第一显示单元46,显示超过第二阈值的第二区域的图像，例如将异常部分在第二图像或第一图像内部呈现出来。[0032]优选地，所述细分单元43采取九宫格的方法，将R0I区域分为多个矩形第二区域，例如将R0I区域分为nXn个矩形第二区域，其中，n是可以根据判断精度自行设定的参数，为任何大于1的整数，优选地，将R0I区域分为4X4个矩形第二区域。[0033]在本发明的一个实施例中，图像分割部2包括：[G034]粗搜索模块M，对第一图像进行搜索，在股骨头边缘选取设定数量的轮廓点，形成初始轮廓；[0035]迭代模块22,通过全局线性回归函数和随机森林对初始轮廓进行迭代，直至得到最小化局部二值特征和最小化线性回归量，得到各轮廓点的最佳位置，形成最佳轮廓形状作为分割线，具体地，迭代模块22包括：[0036]随机森林单元，使用随机森林预测形状变化，在每一次迭代中，构建前一次迭代后的轮廓形状的各轮廓点的树，得到各轮廓点的局部二元特征；[0037]全局线性回归单元，在每一次迭代中，通过上一次迭代的轮廓形状后每个轮廓点对应的局部二值特征对全局线性回归器进行训练，确定形状增量，根据形状P量和卜一冷迭代后的轮廓形状，确定此次迭代的轮廓形状；’0[0038]最佳轮廓形状获得单兀，利用局部二值特征和全局线性回归器分别近似轮廊形状内、外部灰度，通过最小化能量函数进行迭代，得到最小化局部二值特征和最小化线性^归量，形成最佳轮廓形状，rqi2[0059]其中，i为像素点索引，r⑴为像素点i的灰度值的前缀，rq⑴为像素点i量化后的灰度级，Round是四舍五入函数，-表示采用rq⑴替代之前的i。[0060]上述各实施例中的医用图像处理装置将局部二元拟合LBF应用在股骨头医用图像处理，有足够的精确度和较快的分割速度，针对病变的两种分类方式:参数型和纹理型。对以后的病变分析有一定指导意义。[0061]图2是本发明所述医用图像处理方法流程图的示意图，如图2所示，所述医用图像处理方法包括：[0062]步骤S1，存储包含股骨头的第一图像；[0063]步骤S2，分割出仅含有股骨头部分的第二图像；[0064]步骤S3,检验第二图像的边缘是否为圆，显示不为圆的部分的图像；[0065]步骤S4，检测第二图像的纹理特征，包括:在第二图像上选取R0I区域;将R0I区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为0;将R0I区域细分成多个第二区域，R0I区域的面积是第二区域面积的整数倍;对每个第二区域内的0的个数进行计数;判断第二区域的0的个数是否超过第二阈值，显示超过第二阈值的第二区域的图像。[0066]优选地，在步骤S1和步骤S2之间还包括对第一图像进行滤波去噪和对比度增强的预处理步骤，所述预处理步骤包括：[0067]采用中值滤波器对第一图像进行滤波；[0068]统计滤波后的第一图像中各像素点的灰度值出现的次数；[0069]根据所有灰度值出现的总次数，得到各灰度值的出现概率；[0070]确定各灰度值出现概率的前缀，所述前缀是不大于灰度值的各灰度值的出现概率之和；[0071]根据下式⑴和⑵利用各像素点的灰度值的前缀对各像素点进行像素映射，[0072]rq⑴二Roundr⑴X255[0073]i-rqi[0074]其中，i为像素点索引，r⑴为像素点i的灰度值的前缀，rq⑴为像素点i量化后的灰度级，Round是四舍五入函数。[0075]优选地，所述步骤S2包括：[0076]对第一图像进行搜索，在股骨头边缘选取设定数量的轮廓点，形成初始轮廓；[0077]通过全局线性回归函数和随机森林队初始轮廓进行迭代，直至得到最小化局部特征映射量和最小化线性回归量，得到各轮廓点的最佳位置，形成最佳轮廓形状作为分割线。[0078]优选地，所述步骤S3包括：[0079]提取第二图像中股骨头边缘部分的点；[0080]采用拟合方法拟合上述各点对应的圆心和半径；[0081]确定各点到拟合的圆心的距离以及所述距离与拟合的半径的差值的绝对值；[0082]判断上述差值的绝对值是否大于第三阈值，显示所述绝对值大于第三阈值的点；[0083]将各点的所述绝对值相加；[0084]判断上述相加的结果是否大于第四阈值，如果大于第四阈值，发出预警信号。[0085]尽管前面公开的内容示出了本发明的示例性实施例，但是应当注意，在不背离权利要求限定的范围的前提下，可以进行多种改变和修改。根据这里描述的发明实施例的方法权利要求的功能、步骤和或动作不需以任何特定顺序执行。此外，尽管本发明的元素可以以个体形式描述或要求，但是也可以设想具有多个元素，除非明确限制为单个元素。

权利要求：1.一种医用图像处理装置，其特征在于，包括：图像存储部，用于存储包含股骨头的第一图像；图像分割部，分割出仅含有股骨头部分的第二图像；特征检验部，检验第二图像的边缘是否为圆，显示不为圆的部分的图像；纹理检测部，包括:R〇l区域选取单元，在第二图像上选取ROI区域；灰度变换单元，将ROI区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为0;细分单元，将ROI区域细分成多个第二区域，ROI区域面积是第二区域面积的整数倍;计数单元，对每个第二区域内的〇的个数进行计数;第一判断单元，判断第二区域的〇的个数是否超过第二阈值;第一显示单元，显示超过第二阈值的第二区域的图像。2.根据权利要求1所述的医用图像处理装置，其特征在于，所述特征检测部包括：边缘提取单元，提取第二图像中股骨头边缘部分的点；拟合单元，采用拟合方法拟合上述各点对应的圆心和半径；确定单元，确定各点到拟合的圆心的距离以及所述距离与拟合的半径的差值的绝对值；第二判断单元，判断上述差值的绝对值是否大于第三阈值；第二显示单元，显示所述绝对值大于第三阈值的点；加法器，将各点的所述绝对值相加；第三判断单元，判断上述相加的结果是否大于第四阈值，如果大于第四阈值，发出预警信号。3.根据权利要求1所述的医用图像处理装置，其特征在于，还包括预处理部，对第一图像进行滤波去噪和对比度增强。4.根据权利要求3所述的医用图像处理装置，其特征在于，所述预处理部包括中值滤波器和对比度增强模块，其中，所述对比度增强模块包括：读入单元，读入中值滤波器滤波后的第一图像；统计单元，统计第一图像中各像素点的灰度值出现的次数；出现概率确定单元，根据所有灰度值出现的总次数，得到各灰度值的出现概率；前缀确定单元，确定各灰度值出现概率的前缀，所述前缀是不大于灰度值的各灰度值的出现概率之和；均衡化单元，根据下式（1和2利用各像素点的灰度值的前缀对各像素点进行像素映射，rqi=Roundr⑴X2551i-rqi⑵其中，i为像素点索引，ri为像素点i的灰度值的前缀，rqi为像素点i量化后的灰度级，Round是四舍五入函数。5.根据权利要求1所述的医用图像处理装置，其特征在于，所述图像分割部包括：粗搜索模块，对第一图像进行搜索，在股骨头边缘选取设定数量的轮廓点，形成初始轮廓；_迭代模块，通过全局线性回归函数和随机森林对初始轮廓进行迭代，直至得到最小化局部二值特征和最小化线性回归量，得到各轮廓点的最佳位置，形成最佳轮廓形状作为分割线。6.根据权利要求1所述的医用图像处理装置，其特征在于，所述细分单元采取九宫格的方法，将ROI区域分为多个矩形第二区域。7.—种医用图像处理方法，其特征在于，包括：步骤S1，存储包含股骨头的第一图像；步骤S2,分割出仅含有股骨头部分的第二图像；步骤S3，检验第二图像的边缘是否为圆，显示不为圆的部分的图像；步骤S4,检测第二图像的纹理特征，包括:在第二图像上选取R〇l区域;将ROI区域作灰度梯度变换，将灰度梯度大于第一阈值的点置为0;将ROI区域细分成多个第二区域，ROIg域的面积是第二区域面积的整数倍;对每个第二区域内的0的个数进行计数;判断第二区域的〇的个数是否超过第二阈值，显示超过第二阈值的第二区域的图像。8.根据权利要求7所述的医用图像处理方法，其特征在于，所述步骤S3包括：提取第二图像中股骨头边缘部分的点；采用拟合方法拟合上述各点对应的圆心和半径；确定各点到拟合的圆心的距离以及所述距离与拟合的半径的差值的绝对值；判断上述差值的绝对值是否大于第三阈值，显示所述绝对值大于第三阈值的点；将各点的所述绝对值相加；判断上述相加的结果是否大于第四阈值，如果大于第四阈值，发出预警信号。9.根据权利要求7所述的医用图像处理方法，其特征在于，所述步骤S2包括：对第一图像进行搜索，在股骨头边缘选取设定数量的轮廓点，形成初始轮廓；通过全局线性回归函数和随机森林对初始轮廓进行迭代，直至得到最小化局部二值特征和最小化线性回归量，得到各轮廓点的最佳位置，形成最佳轮廓形状作为分割线。10.根据权利要求7所述的医用图像处理方法，其特征在于，在步骤S1和步骤S2之间还包括对第一图像进行滤波去噪和对比度增强的预处理步骤，所述预处理步骤包括：采用中值滤波器对第一图像进行滤波；统计滤波后的第一图像中各像素点的灰度值出现的次数；根据所有灰度值出现的总次数，得到各灰度值的出现概率；确定各灰度值出现概率的前缀，所述前缀是不大于灰度值的各灰度值的出现概率之和；根据下式⑴和⑵利用各像素点的灰度值的前缀对各像素点进行像素映射，rq⑴=RoundriX2551i-rqi2其中，i为像素点索引，r⑴为像素点i的灰度值的前缀，rqi为像素点i量化后的灰度级，Round是四舍五入函数。

百度查询： 清华大学 医用图像处理装置及医用图像处理方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD