一种基于像素信号能量的图像边缘roi的确定算法
【专利摘要】本发明涉及边缘检测技术领域，公开了一种基于像素信号能量的图像边缘ROI的确定算法。具体包括以下步骤：图像矩阵记为I，选取图像矩阵上第x行第y列为中心像素点I(x,y)；以像素点I(x,y)为中心的像素块的块内灰度方差用下式表示：V(x，y)＝var{G(x+i，y+j)},i,j∈{?1,0,1}，其中G(x,y)为像素点I(x,y)的灰度值；记像素点I(x,y)的信号能量En(x,y)＝V(x,y)2；引入信号能量的百分位数函数Percentile(u,v)，其中u,v∈[0,100]，表示矩阵I中像素点I(x,y)的信号能量En(x,y)小于总能量百分之v的像素点个数占到了矩阵I的像素点总个数的百分之u；根据信号能量的百分位数函数Percentile(u,v)，选取显著性水平α的值，计算1?α取值，剔除信号能量的百分位数函数中像素能量在总能量中占比低于1?α的像素点，所述选出的像素点有1?α的可能性属于图像边缘，从而可以获得图像边缘ROI。
【专利说明】
-种基于像素信号能量的图像边缘ROI的确定算法
技术领域
[0001] 本发明设及边缘检测技术领域，特别是一种基于像素信号能量的图像边缘R0I的 确定算法。
【背景技术】
[0002] 目前传统图像边缘检测技术采用数字图像边缘检测技术，主要目的是为了检测或 者提取出一张数字图像中有明显变化的边缘或者不联系区域包括一阶微分边缘算子、Log 算子、Canny算子、Sobel算子等，虽然数字图像边缘检测技术可W较好的得到图像的边缘， 但是检测算法多根据像素的局部情况进行判定，而且检测到的边缘轮廓常常是不连续的， 甚至包含一些人们不感兴趣的区域边缘，因此检测出的图像边缘往往存在个别的误检，比 如局部脱离边缘主轮廓的点明显是非边缘区域被判断成边缘，或者图像中目标内部出现了 小幅灰度变化的区域被误判为边缘区域。如图1所示，该图为CT扫描得到的牙齿图像，可W 看出图1中的背景区域有大量区域因其附近像素的灰度出现了波动而被误判为边缘，如图2 所示为利用传统化nny算子处理图1后所获得的边缘图像且在牙齿区域内部也出现了很多 脱离边缘主轮廓的非边缘区域。同样在其他物体的图像边缘区域同样也会出现上述的误判 情况。
[0003] 故有必要设计算法去除非边缘区域，在边缘检测前设定大致的图像边缘R0I，在确 定边缘可能出现的区域后，只在该区域进行边缘检测，W提高图像边缘检测的准确性。
【发明内容】
[0004] 本发明的目的是确定图像边缘可能出现的区域，从而设置图像边缘R0I，减少传统 的数字图像边缘检测算法带来的边缘误检，本发明提出了一种基于像素信号能量的图像边 缘R0I的确定算法。
[0005] 本发明采用的技术方案如下：
[0006] -种基于像素信号能量的图像边缘R0I的确定算法，其特征在于包括W下步骤：
[0007] 步骤一、图像矩阵记为I，选取图像矩阵上第X行第y列为中屯、像素点I(x，y);
[000引步骤二、W像素点I(x，y)为中屯、的像素块的块内灰度方差用下式表示:V(x，y) = va;r{G(x+i，y+j)}，i, je {-1,0,1}，其中G(x，y)为像素点I(x，y)的灰度值；
[0009] 步骤S、记像素点I(x，y)的信号能量6。山7)=^《，7)2，像素点是边缘的可能性大 小建模为像素点I(x，y)的信号能量；
[0010] 步骤四、引入信号能量的百分位数函数Pe;rcentile(u，v)，其中U，ve [0,100]，表 示矩阵I中像素点I(x，y)的信号能量En(x，y)小于总能量百分之V的像素点个数占到了矩阵 I的像素点总个数的百分之U;
[0011] 步骤五、用二维坐标图表示百分位数函数化rcentile(u，V)，选取显著性水平α，可 求得1-α的值；
[0012] 步骤六、在百分位数函数化;rcentile(u，v)中将像素点I(x，y)信号能量En(x，y)占 总能量的比低于1-α的像素点全部剔除，即将其灰度置零，可W获得图像边缘ROI。
[0013] 进一步地，上述信号能量有W下关系式
其中 I二{I (义，5〇|6。（义，5〇<化1；虹6311。1(1},化1；虹6311。1(1为像素{曰亏能里阔值。
[0014] 与现有技术相比，采用上述技术达到W下有益效果:根据图像像素块的整体情况 进行判定，W像素块的灰度方差作为判断标准，并且引入信号能量和能量百分位数函数，可 有效剔除属于图像边缘的可能性较小的区域。
【附图说明】
[0015] 图1为通过CT机扫描牙齿得到的原始图像。
[0016] 图2为利用传统化nny算子处理图1后所获得的边缘图像。
[0017] 图3为本发明图像边缘的像素信号能量图。
[001引图4为本发明信号能量的百分位数函数化rcenti le (U，V)的二维坐标图。
[0019] 图5为本发明基于像素信号能量的图像边缘R0I图。
【具体实施方式】
[0020] 现结合附图对本发明做进一步地描述。
[0021] -种基于像素信号能量的图像边缘R0I的确定算法，具体包括W下步骤：
[0022] 本实施例W牙齿为例，除了牙齿本发明还适用于机械部件、精密零件等其他需要 计算图像边缘R0I的物体；
[0023] 步骤一、如图1所示为通过CT机扫描物体得到的原始图像，得到需要进行边缘检测 的初始图像，图像中可W看到被检测物体的整体轮廓，图像矩阵记为I，选取图像矩阵上第X 行第y列为中屯、像素点I(x，y);
[0024] 步骤二、W像素点I(x，y)为中屯、的像素块的块内灰度方差用下式表示:V(x，y) = va;r{G(x+i，y+j)}，i，je {-1,0,1}，其中G(x，y)为像素点I(x，y)的灰度值，图像的边缘区域 的各个像素点之间会呈现出较大的灰度差，而像素块的灰度方差可W用来描述运一灰度变 化的程度，且用到了像素块内的全部像素点的信息，因此可W将像素块的灰度方差作为某 像素点属于图像边缘可能性大小的评判标准，V(x，y)的数值越小代表像素点I(x，y)越有可 能不是图像边缘，应该删去；
[0025] 步骤Ξ、为确定需要删去的部分，需要确定删除阔值，运里引入像素点的信号能量 的概念，记像素点I(x，y)的信号能量6。(^，7)=八^，7)2，像素点1(^，7)是边缘的可能性大小 建模为像素点I(x，y)的信号能量，矩阵I中能量积累小的区域应当考虑被剔除，可得到像素 信号能量图如图3所示；
[00%]步骤四、为了进一步确定图像边缘R0I，引入信号能量的百分位数函数化rcentile
(u，v)，其中！ Γ = {Γ (x，y) |En(x, ,， y)<E：nthresh〇id} Jnthreshoid为像素信号能量阔值，百分位数函数化rcentile(u,v)表示矩阵I 中像素点I(x，y)的信号能量En(x，y)小于总能量百分之V的像素点个数占到了矩阵I的像素 点总个数的百分之U，可得到像素信号能量的百分位数函数化rcenti 1 e (U，V);
[0027] 步骤五、如图4中所示，用二维坐标图表示百分位数函数化;rcentile(u,v)，横坐标 为数量百分比U，纵坐标为能量百分比V，选取显著性水平α，能量百分比V在能量百分比为 1%时候突然急速增大，所W此实施例中显著性水平α取1%，可求得1-α的值为99%。在百分 位数函数中将像素点的信号能量占总能量的比低于1-0 = 99%对应的像素点全部剔除，即 将其灰度值置零，即可W得到图像边缘R0I，如图4中对应能量百分比V为1%的数量百分比 为95.5%，说明整个图像中95.5%的像素的能量只占到了总信号能量的1%，运些点属于图 像边缘的可能性低，可W剔除，剩余4.5%的像素占到总能量的99%可W保留，由此可W得 到图像边缘R0I，如图5所示。
[0028] 虽然本发明的技术内容已经W较佳的实例公开如上，然而并非用W限定本发明， 如果本领域技术人员，在不脱离本发明的精神所做的非实质性改变或改进，都应该属于本 发明权利要求保护的范围。
【主权项】
1. 一种基于像素信号能量的图像边缘ROI的确定算法，其特征在于包括以下步骤： 步骤一、图像矩阵记为I，选取图像矩阵上第X行第y列为中心像素点I(x，y); 步骤二、以像素点I(x，y)为中心的像素块的块内灰度方差用下式表示:V(X，y)=var{G (x+i，y+j)}，i, je {-1,0,1}，其中G(x，y)为像素点I(x，y)的灰度值； 步骤三、记像素点I(x，y)的信号能量En(x，y)=V(x，y)2,像素点是边缘的可能性大小建 模为像素点I(x，y)的信号能量； 步骤四、引入信号能量的百分位数函数Percentile(u，v)，其中u，ve [〇, 1〇〇]，表示矩 阵I中像素点I(x，y)的信号能量En(x，y)小于总能量百分之v的像素点个数占到了矩阵I的 像素点总个数的百分之u; 步骤五、用二维坐标图表示百分位数函数PerCentile(U，V)，选取显著性水平α，可求得 l-α的值； 步骤六、在百分位数函数Percentile(u，v)中将像素点I(x，y)信号能量En(x，y)占总能 量的比低于l-α的像素点全部剔除，即将其灰度置零，可以获得图像边缘R0I。2. 如权利要求1所示的基于像素信号能量的图像边缘R0I的确定算法，其特征在于所述 信号能量有以下关系式，其中r = {I'（X，y)|En(X，y) <Enthreshold}，Enthreshold为像素信号能量阈值。
【文档编号】G06T7/00GK106067176SQ201610409897
【公开日】2016年11月2日
【申请日】2016年6月12日 公开号201610409897.2, CN 106067176 A, CN 106067176A, CN 201610409897, CN-A-106067176, CN106067176 A, CN106067176A, CN201610409897, CN201610409897.2
【发明人】罗敏, 余翊森, 蒋千军, 王浩正
【申请人】成都甄识科技有限公司