2017-05-26 19:10:29 coming_is_winter 阅读数 11329
  • 软件测试入门视频教程

    软件测试入门视频培训教程:该课程将带你走进“软件测试”的大门,具体内容包括软件测试环境搭建、软件开发模型、产品模型、CMM模型、测试用例、等价类划分、边界值划分、白盒测试、单元测试、bugfree搭建、系统测试、回归测试、验收测试等。本课程以接地气的语言来讲解,让你听的懂,学的会!本课程以全新的方式为你呈现教学内容,清新脱俗独具特色的授课方式将带给你新的体验。

    2158191 人正在学习 去看看 李晓鹏

                         边界跟踪:内边界跟踪算法解释

 

        边界跟踪是基于边缘的分割常用方法之一,用于区域已分出(二值或已标注),但边界未知的情况。分为内边界与外边界。内边界为区域的一个子集,外边界不是区域的一个子集。

        内边界跟踪算法具体的实施方法比较明确,但在学习的时候感觉稍微有点抽象,结合自己的理解拓展一下。
        具体算法:

                                                                                                       图1

 

                                                                                                         图2


       具体有四步,以上图片中包含四邻域与八邻域,下面主要以以八邻域为例,四邻域原理相似。

 

       第一步中图像处理是以整个图片的最左上角为坐标原点(0,0),所以最靠近左上角的P0为最小行数与最小列数之像素。定义的dir为边界移动的方向,具体为八邻域中的八个方向(0,1,2,3,4,5,6,7)中的一个(如图2中的(b)图),此时默认初始值为7(右下方向)。

       第二步中先判断dir的奇偶性(默认值为7),然后对dir进行计算更改(就是根据上下文调整邻域扫描的初始方向))进行下一轮边界扫描。具体操作看图,画的不是很好,有点乱。

  

 

      这里的奇偶判断及为何是+6,+7具体原理尚未能有效想明白,有待进一步学习。

      第三步比较好理解边界是一个闭合空间,搜索到首尾相接时即结束。

     第四步及将计算了两遍的Pn-1与Pn删除,最后边界为P0到Pn-2。

     水平有限,理解可能有误,加之文笔不好,凑合看哈。

     PS:后边发现有一个不确定的地方,搜索P0的时候dir是默认的7还是先判断dir的奇偶性再按第二步的步骤进行呢?这里先按初始7计算。

     参考文献:

     《图像处理、分析与机器视觉(第三版)》     作者 :Milan Sonka、艾海舟。

 

2017-05-14 22:19:10 PinkRiverside 阅读数 1010
  • 软件测试入门视频教程

    软件测试入门视频培训教程:该课程将带你走进“软件测试”的大门,具体内容包括软件测试环境搭建、软件开发模型、产品模型、CMM模型、测试用例、等价类划分、边界值划分、白盒测试、单元测试、bugfree搭建、系统测试、回归测试、验收测试等。本课程以接地气的语言来讲解,让你听的懂,学的会!本课程以全新的方式为你呈现教学内容,清新脱俗独具特色的授课方式将带给你新的体验。

    2158191 人正在学习 去看看 李晓鹏

边界跟踪

利用边缘检测算法可以检测出图像中的边界点。但是在很多情况下,仅仅检测出边界点是不够地,必须通过边界跟踪得到边界点序列等数据,为图像分析做准备。对二值图像的边界跟踪可基于八个方向码进行,如图所示。

这里写图片描述

P(x,y) 为物体的一个边界点,则P(x,y) 的下一个边界点必存在其八邻域内,因此可以根据八邻域信息进行外边界跟踪。在找到下一个边界点后,依次边界点为当前边界点继续搜索。这一搜索过程不断重复下去,直到搜索至起点。

图像分割

由边界跟踪结果可得到图像中目标区域边界的提取,将图像分为若干互不交迭的区域,并使这些特征在同一区域内呈现出相似性。

这里写图片描述

2017-03-09 10:40:26 zhangyibo123456789 阅读数 6794
  • 软件测试入门视频教程

    软件测试入门视频培训教程:该课程将带你走进“软件测试”的大门,具体内容包括软件测试环境搭建、软件开发模型、产品模型、CMM模型、测试用例、等价类划分、边界值划分、白盒测试、单元测试、bugfree搭建、系统测试、回归测试、验收测试等。本课程以接地气的语言来讲解,让你听的懂,学的会!本课程以全新的方式为你呈现教学内容,清新脱俗独具特色的授课方式将带给你新的体验。

    2158191 人正在学习 去看看 李晓鹏

 matlab图像处理之二值图像内外边界跟踪

注:原文链接:http://www.cnblogs.com/tiandsp/archive/2013/04/26/3045747.html
       目标内边界的像素全都在目标里面,目标外边界的像素全都不在目标上,是包围着目标的。
二值图像内外边界的计算都是有两种方法的,所以一共是4种算法,不过实际用到跟踪的只有一个而已。
 
       首先是内边界跟踪:
第一种方法不是跟踪方法。步骤是先对原图像腐蚀,然后用原图像减去腐蚀后的图像就得到边界了。
第二种方法是跟踪方法。步骤如下:
1.遍历图像。
2.标记第一个遇见像素块的前景像素(i,j)。
3.对这个像素周围八邻域逆时针搜索,如果搜索到周围有前景像素,那么更新坐标(i,j)为(i',j'),并标记。
4.不断执行第3步直到再次遇见此像素块第一次标记的像素。
5.继续执行第1步。
 
     然后是外边界跟踪:
第一种方法和求内边界第一种方法类似。先对原图像进行膨胀,然后用膨胀后的图像减去原图像即可。
     第二种也不算跟踪方法,只是标记算法而已。就是将图像中前景像素周围的非前景像素标记一下就行了。
 
效果如下:

原图:

内边界:

外边界:

matlab程序如下:

内边界:

clear all;
close all;
clc;

img=imread('rice.png');
img=img>128;
imshow(img);
[m n]=size(img);

imgn=zeros(m,n);        %边界标记图像
ed=[-1 -1;0 -1;1 -1;1 0;1 1;0 1;-1 1;-1 0]; %从左上角像素,逆时针搜索
for i=2:m-1
    for j=2:n-1
        if img(i,j)==1 && imgn(i,j)==0      %当前是没标记的白色像素
            if sum(sum(img(i-1:i+1,j-1:j+1)))~=9    %块内部的白像素不标记
                ii=i;         %像素块内部搜寻使用的坐标
                jj=j;
                imgn(i,j)=2;    %本像素块第一个标记的边界,第一个边界像素为2
                
                while imgn(ii,jj)~=2    %是否沿着像素块搜寻一圈了。
                    for k=1:8           %逆时针八邻域搜索
                        tmpi=ii+ed(k,1);        %八邻域临时坐标
                        tmpj=jj+ed(k,2);
                        if img(tmpi,tmpj)==1 && imgn(tmpi,tmpj)~=2  %搜索到新边界,并且没有搜索一圈
                            ii=tmpi;        %更新内部搜寻坐标,继续搜索
                            jj=tmpj;
                            imgn(ii,jj)=1;  %边界标记图像该像素标记,普通边界为1
                            break;
                        end
                    end
                end
                
            end
        end
    end
end

figure;
imgn=imgn>=1;
imshow(imgn,[]);

%不过要是真取二值图像内边界,通常是原图减去其腐蚀图就行了
se = strel('square',3); 
imgn=img-imerode(img,se);    
figure;
imshow(imgn)

外边界:

clear all;
close all;
clc;

img=imread('rice.png');
img=img>128;
imshow(img);
[m n]=size(img);

imgn=zeros(m,n);        %边界标记图像
ed=[-1 -1;0 -1;1 -1;1 0;1 1;0 1;-1 1;-1 0]; %从左上角像素判断
for i=2:m-1
    for j=2:n-1
        if img(i,j)==1      %如果当前像素是前景像素
            
            for k=1:8
                ii=i+ed(k,1);
                jj=j+ed(k,2);
                if img(ii,jj)==0    %当前像素周围如果是背景,边界标记图像相应像素标记
                    imgn(ii,jj)=1;
                end
            end
            
        end
    end
end
    
figure;
imshow(imgn,[]);

%不过要是真取二值图像外边界,通常是原图膨胀图减去原图就行了
se = strel('square',3); 
imgn=imdilate(img,se)-img;    
figure;
imshow(imgn)



2018-07-20 10:39:44 qq_38991255 阅读数 3200
  • 软件测试入门视频教程

    软件测试入门视频培训教程:该课程将带你走进“软件测试”的大门,具体内容包括软件测试环境搭建、软件开发模型、产品模型、CMM模型、测试用例、等价类划分、边界值划分、白盒测试、单元测试、bugfree搭建、系统测试、回归测试、验收测试等。本课程以接地气的语言来讲解,让你听的懂,学的会!本课程以全新的方式为你呈现教学内容,清新脱俗独具特色的授课方式将带给你新的体验。

    2158191 人正在学习 去看看 李晓鹏

前几篇文章中讲的都是对于离散点的操作,对于离散点的处理方法比较少,后来想到将离散点做出的曲线图保存成图像,对图像进行处理比较方便。在matlab图像处理工具箱中,有两个函数可以用来进行边界跟踪,其中一个是bwtraceboundary,另一个是bwboundaries。我们首先看一下bwtraceboundary 函数的常见调用方法:
B = bwtraceboundary(BW, p,fstep)
B = bwtraceboundary (BW, P,fstep, conn)
B = bwtraceboundary(BW, p, fstep, conn, N, dir)
其中,BW是二值图像,非0值构成了目标,0值构成了背景; p是一个二维向量,p中的元素是指开始搜寻边界点的行号和列号; fstep 是初始的搜索方向,其值可以为'N",即图像上面的元素,若值为'NE'则表示图像右上方的元素; conn刻画了跟踪边界的连续性,其值可以为4或者8; N是指边界的最大像素数目,默认值为inf; dir 是跟踪边界的方向,其值可以为'clockwise'或者counterclockwise';B是返回的边界像素的行号和列号。关于搜索方向的设定就不在这里加以赘述,在help中可以查到。

下面的例子首先读取一幅二值图像,从图像的左上方开始搜寻图像的第一个非0像素,边界的像素个数为50,然后在原图上显示边界。

利用bwtraceboundary函数跟踪边界,代码实现如下:

BW=imread('BW.jpg');                           %读取二值图像
imshow(BW,[]);                                 %显示二值图像
s=size(BW);                                    %图像矩阵的行和列
for row=2:55:s(1)
    for col=1:s(2)
        if BW(row,col),
            break;
        end
    end
    contour = bwtraceboundary(BW,[row,col],'W',8,50,...
        'counterclockwise');                                   %检测是否为边界
    if(~isempty(contour))
        hold on;
        plot(contour(:,2),contour(:,1),'g','LineWidth',2);     %画出边界
        hold on;
        plot(col,row,'gx','LineWidth',2);                      %画出边界起点
    else
        hold on;plot(col,row,'rx','LineWidth',2);
    end
end

另外一种方法将在下一篇博客中进行讲解。

2017-11-10 21:00:36 hustllx 阅读数 1106
  • 软件测试入门视频教程

    软件测试入门视频培训教程:该课程将带你走进“软件测试”的大门,具体内容包括软件测试环境搭建、软件开发模型、产品模型、CMM模型、测试用例、等价类划分、边界值划分、白盒测试、单元测试、bugfree搭建、系统测试、回归测试、验收测试等。本课程以接地气的语言来讲解,让你听的懂,学的会!本课程以全新的方式为你呈现教学内容,清新脱俗独具特色的授课方式将带给你新的体验。

    2158191 人正在学习 去看看 李晓鹏

matlab函数:B=bwboundaries(img);

使用栈实现代码:

close all;clear all;
img = imread('G:/bound.jpg');
img = rgb2gray(img);
img = im2bw(img);
img=im2bw(img);
B=bwboundaries(img);%库函数结果

[m,n]=size(img);
k=1;
flag=0;
for i=1:m %找到第一个不为0元素位置
    for j=1:n
        if img(i,j)>0
            stack(k,1)=i;
            stack(k,2)=j;
            starti=i;   %标记第一个点
            startj=j;
            dir=7;
            k=k+1;
            flag=1; %表示已经找到了第一个元素
            break;
        end
    end
    if flag==1
        break;
    end
end

while 1
    if mod(dir,2)==0
        dir=mod(dir+7,8); %偶数处理方式
    else
        dir=mod(dir+6,8); %奇数的处理方式
    end
    i=stack(k-1,1);
    j=stack(k-1,2);
    while 1
        flag=1;
        switch dir
            case 0
                if j+1<n+1&&img(i,j+1)>0%右方
                    stack(k,1)=i;
                    stack(k,2)=j+1;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                    dir=mod(dir+1,8);
                end
            case 1
                if i-1>0&&j+1<n+1&&img(i-1,j+1)>0%右上方
                    stack(k,1)=i-1;stack(k,2)=j+1;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                     dir=mod(dir+1,8);
                end
            case 2
                if i-1>0&&img(i-1,j)>0%上方
                    stack(k,1)=i-1;stack(k,2)=j;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                    dir=mod(dir+1,8);
                end
            case 3
                if i-1>0&&j-1>0&&img(i-1,j-1)>0%左上方
                    stack(k,1)=i-1;stack(k,2)=j-1;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                     dir=mod(dir+1,8);
                end
            case 4
                if j-1>0&&img(i,j-1)>0%左方
                    stack(k,1)=i;stack(k,2)=j-1;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                    dir=mod(dir+1,8);
                end
            case 5
                if i+1<m+1&&j-1>0&&img(i+1,j-1)>0%左下方
                    stack(k,1)=i+1;stack(k,2)=j-1;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                    dir=mod(dir+1,8);
                end
            case 6
                if i+1<m+1&&img(i+1,j)>0%下方
                    stack(k,1)=i+1;stack(k,2)=j;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                    dir=mod(dir+1,8);
                end
            case 7
                if i+1<m+1&&j+1<n+1&&img(i+1,j+1)>0%右下方
                    stack(k,1)=i+1;stack(k,2)=j+1;
                    k=k+1;
                    flag=0;%找到下一个元素位置
                else
                   dir=mod(dir+1,8);
                end
            otherwise
        end
        if flag==0%找到下一个点退出当前循环
            break;
        end
    end
    if  k-1>2&&stack(k-2,1)==starti&&stack(k-2,2)==startj&&stack(k-1,1)==stack(2,1)&&stack(k-1,2)==stack(2,2)%退出循环条件
        break;
    end
end
stack(k-1,:)=[];

figure;
subplot(1,2,1);
imshow(img);
subplot(1,2,2);
imshow(img);
hold on;
for i=1:length(stack)
    hold on; plot(stack(:,2),stack(:,1),'r');
end

使用递归实现:

没有更多推荐了,返回首页