异平均数值而取得。4. 如权利要求1所述的量化方法，其中该组像素点的范围为一kXk方阵，其中k为自 然数。5.-种回音纹理特征的可视化方法，包含： 取得一超声波影像，所述超声波影像包含数个像素点，每一该像素点具有一像素点数 值； 接收所述超声波影像中的一选取区域，并计算该选取区域内所有该像素点数值以得到 一区域标准差值； 排除该选取区域中像素点数值小于一基准值的像素点，其中所述基准值为该区域标准 差值的倍率值； 针对剩余的像素点，分别统计以一第η像素点为中心的一组像素点，计算出对应的一 组像素点数值的一第η局部统计参数，η为1至剩余像素点的数目，其中η为自然数，所述 局部统计参数包含数个局部平均值； 根据所述局部统计参数取得至少一回音纹理指标； 计算剩余像素点数值，以取得一区域平均值，且分别计算所述局部平均值与该区域平 均值的绝对差值以取得每一该像素点的一回音纹理强度可视化值；以及 结合回音纹理强度可视化值与回音纹理指标，绘制一回音纹理影像化色阶图。6. 如权利要求5所述的可视化方法，进一步包含： 比较每一回音纹理强度可视化值与一上界值与一下界值，以得到一可视化呈现值； 当任一像素点的回音纹理强度可视化值大于所述上界值时，可视化呈现值为〇 ; 当任一像素点的回音纹理强度可视化值小于所述下界值时，可视化呈现值为一色阶 值；以及 当任一像素点的回音纹理强度可视化值小于所述上界值且大于所述下界值时，可视化 呈现值为该上界值减去回音纹理强度可视化值除以所述上界值减去所述下界值，再乘以所 述色阶值。7. 如权利要求6所述的可视化方法，其中所述局部统计参数进一步包含数个局部变异 数值，所述方法进一步包含：计算剩余像素点数值，以取得一区域标准差参数；以及计算所 述局部平均值以及所述局部变异数值以得到一局部变异平均数值及一局部变异的变异值， 所述上界值为局部变异平均数值除以局部变异的变异值，再乘以区域标准差参数除以一上 界参数，且所述下界值为该局部变异平均数值除以该局部变异的变异值，再乘以所述区域 标准差参数除以一下界参数。8. 如权利要求6所述的可视化方法，其中所述局部统计参数进一步包含数个局部变异 系数，所述方法进一步包含：计算剩余像素点数值，以取得一区域标准差参数；以及计算所 述局部变异系数以得到一局部变异平均系数，所述上界值为该局部变异平均系数的倒数， 再乘以所述区域标准差参数除以一上界参数，且该下界值为所述局部变异平均系数的倒 数，再乘以所述区域标准差参数除以一下界参数。9. 如权利要求7或8所述的量化方法，其中所述上界参数及下界参数为一对比强度参 数，所述上界参数及该下界参数介于0. 5~10. 0之间，且该下界参数大于该上界参数。10. 如权利要求5所述的可视化方法，其中该组像素点的范围为一kXk方阵，其中k为 自然数。11. 一种回音纹理特征的量化系统，包含： 一超声波影像擷取装置，用于取得一超声波影像，该超声波影像包含数个像素点，每一 该像素点具有一像素点数值； 一影像处理模块，连接该超声波影像擷取装置，用于接收及处理该超声波影像，所述影 像处理模块包含一分析单元及一运算单元，该分析单元与该运算单元相互连接； 一显示装置，连接影像处理模块，用于显示该超声波影像；以及 一输入装置，连接影像处理模块，用于选取该超声波影像中的一选取区域； 其中，分析单元接收该选取区域，运算单元计算该选取区域内所有该像素点数值以得 到一区域标准差值，该分析单元排除该选取区域中像素点数值小于一基准值的该像素点， 其中所述基准值为该区域标准差值的倍率值，该运算单元针对剩余像素点，分别统计以一 第η像素点为中心的一组像素点，计算出对应的一组像素点数值的一第η局部统计参数，η 为1至剩余像素点的数目，其中η为自然数，该运算单元根据该些局部统计参数取得至少一 回音纹理指标。12. 如权利要求11所述的量化系统，其中所述局部统计参数包含数个局部变异系数， 所述运算单元计算所述局部变异系数以得到一局部变异平均系数而取得该回音纹理指标。13. 如权利要求11所述的量化系统，其中所述局部统计参数包含数个局部平均值以及 数个局部变异数值，所述运算单元计算所述局部平均值以及该些局部变异数值以得到一局 部变异平均数值及一局部变异的变异值，再将所述局部变异的变异值的平方根除以所述局 部变异平均数值而取得该回音纹理指标。14.一种回音纹理特征的可视化系统，包含： 一超声波影像擷取装置，用于取得一超声波影像，所述超声波影像包含数个像素点，每 一像素点具有一像素点数值； 一影像处理模块，连接该超声波影像擷取装置，用于接收及处理所述超声波影像，该影 像处理模块包含一分析单元及一运算单元，分析单元与运算单元相互连接； 一显示装置，连接所述影像处理模块，用以显示该超声波影像；以及 一输入装置，连接所述影像处理模块，用以选取该超声波影像中的一选取区域； 其中该分析单元接收该选取区域，该运算单元计算该选取区域内所有该像素点数值 以得到一区域标准差值，该分析单元排除该选取区域中像素点数值小于一基准值的该像素 点，其中该基准值为该区域标准差值的倍率值，该运算单元针对剩余像素点，分别统计以一 第η像素点为中心的一组像素点，计算出对应的一组像素点数值的一第η局部统计参数， η为1至剩余像素点的数目，其中η为自然数，所述局部统计参数包含数个局部平均值，该 运算单元根据所述局部统计参数取得至少一回音纹理指针，该运算单元计算剩余像素点数 值，以取得一区域平均值，且分别计算所述局部平均值与该区域平均值的绝对差值以取得 每一该像素点的一回音纹理强度可视化值，该分析单元结合该回音纹理强度可视化值与该 回音纹理指标用以绘制一回音纹理影像化色阶图。15. 如权利要求15所述的可视化系统，其中所述分析单元进一步比较每一该回音纹理 强度可视化值与一上界值与一下界值以得到一可视化呈现值，当任一像素点的回音纹理强 度可视化值大于上界值时，该可视化呈现值为〇 ;当任一像素点的该回音纹理强度可视化 值小于下界值时，该可视化呈现值为一色阶值；以及当任一像素点的回音纹理强度可视化 值小于上界值且大于下界值时，该可视化呈现值为上界值减去该回音纹理强度可视化值除 以上界值减去下界值，再乘以该色阶值。16. 如权利要求15的可视化系统，其中所述局部统计参数进一步包含数个局部变异数 值，所述运算单元计算剩余像素点数值，以取得一区域标准差参数，并计算所述局部平均值 以及所述局部变异数值以得到一局部变异平均数值及一局部变异的变异值，所述运算单元 将所述局部变异平均数值除以所述局部变异的变异值，再乘以所述区域标准差参数除以一 上界参数以取得该上界值，且将该局部变异平均数值除以该局部变异上网变异值，再乘以 所述区域标准差参数除以一下界参数以取得该下界值。17. 如权利要求15所述的可视化系统，其中所述局部统计参数进一步包含数个局部变 异系数，所述运算单元计算剩余像素点数值，以取得一区域标准差参数，并计算所述局部变 异系数以得到一局部变异平均系数，该运算单元将所述局部变异平均系数的倒数，再乘以 该区域标准差参数除以一上界参数以取得该上界值，且将该局部变异平均系数的倒数，再 乘以该区域标准差参数除以一下界参数以得到该下界值。18. 如权利要求16或17所述的可视化系统，其中所述上界参数及所述下界参数为一对 比强度参数，所述上界参数及所述下界参数介于〇. 5~10. 0间，且下界参数大于上界参数。
【专利摘要】本发明公开了一种回音纹理特征的量化及可视化方法，所述方法包含：取得超声波影像；计算超声波影像中选取区域内所有像素点数值以得到区域标准差值；排除选取区域中，像素点数值小于区域标准差值的倍率值的像素点；针对剩余的像素点，分别统计以第n像素点为中心的一组像素点，计算出对应的一组像素点数值的第n局部平均值、第n局部变异数值及第n局部变异系数，n为1至剩余像素点的数目；以及根据这些局部平均值、局部变异数值、局部变异系数或其组合以取得回音纹理指标。可视化则以连续色阶的呈现来显示纹理分布情形，并透过可视化对比强度参数来调整不同的可视化呈现的对比度。
【IPC分类】A61B8/00
【公开号】CN105286912
【申请号】CN201410365753
【发明人】张金坚, 陈正刚, 陈炯年, 陈坤源, 王昱欣, 黄国祯
【申请人】安克生医股份有限公司
【公开日】2016年2月3日
【申请日】2014年7月29日