应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 分割图的邻接矩阵仿真
21/1返回列表
查看: 515  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiyue6911

新虫 (小有名气)

[交流] 分割图的邻接矩阵仿真

小弟最近研究基于对象的目标检测，第一步进行多尺度分割，
右图为待转化的分割图像，左图为邻接矩阵，行和列表示[1 2 3 4 5 6 7];
关键是，区域相连判断，思路：遍历区域边界，unique边界点的8邻域；
现存问题：如何将邻接矩阵进行可视化，要用labview吗？
自己写的将分割好的图像转化为邻接矩阵：
CODE:
x = repmat(uint8(0),[512 512]);
x(1:512,256)=255;x(170,1:256)=255;x(340,1:256)=255;x(85,256:356)=255;x(170,356:512)=255;x(340,356:512)=255;
x(85:340,356)=255;x(390,356:456)=255;x(450,356:456)=255;x(390:450,356)=255;x(390:450,456)=255;
x = ~im2bw(x);
[L num] = bwlabel(x,8);%邻接区域个数
regionmat = zeros(num);%存储邻接矩阵
imshow(L);
stats = regionprops(L);%构造区域个节点；
for i = 1:num
    node(i).id = i;
    node(i).data = stats(i);
    node(i).relaedge = [1 1 1];
    node(i).weight = [1 1 1];
end
[r c] = find(L==0);%遍历所有边界点
%遍历边界点的8邻域，构造区域邻接集合
[r1 c1] = size(L);
L_new = zeros(r1+2,c1+2);
L_new(2:(size(L_new,1)-1),2:(size(L_new,2)-1)) = L;
col = 8;%最多连通区域；
for i = 1:(size(r,1))
    temp = [L_new(r(i)-1+1,c(i)-1+1) L_new(r(i)-1+1,c(i)+1) L_new(r(i)-1+1,c(i)+1+1);
                    L_new(r(i)+1,c(i)-1+1) L_new(r(i)+1,c(i)+1) L_new(r(i)+1,c(i)+1+1);
                    L_new(r(i)+1+1,c(i)-1+1) L_new(r(i)+1+1,c(i)+1) L_new(r(i)+1+1,c(i)+1+1);];
    p = zeros(1,col);
    temp = unique(temp(:));
    temp = setdiff(temp,0);
    p(1:length(temp)) = temp;
    relas(i,:) = p;
end
relas = unique(relas,'rows');
%赋值邻接矩阵
for i = 1:size(relas,1)
    temp =setdiff(relas(i,:),0);
    cn = nchoosek(temp,2);
    for j =1:size(cn,1)
        regionmat(cn(j,1),cn(j,2)) = 1;
    end
end
disp('regionmat:');disp(regionmat);

分割图的邻接矩阵仿真
测试.JPG

[ Last edited by dbb627 on 2015-11-14 at 14:46 ]
回复此楼

» 猜你喜欢
1楼 2015-10-27 17:37:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiyue6911

新虫 (小有名气)
升级版本
直接运行可以生成可视化图形结构
%分割图的邻接矩阵仿真
CODE:
x = repmat(uint8(0),[512 512]);
x(1:512,256)=255;x(170,1:256)=255;x(340,1:256)=255;x(85,256:356)=255;x(170,356:512)=255;x(340,356:512)=255;
x(85:340,356)=255;x(390,356:456)=255;x(450,356:456)=255;x(390:450,356)=255;x(390:450,456)=255;
locations = [85,128;265,128;425,128;425,400;265,300;265,400;85,400];
x = ~im2bw(x);
[L num] = bwlabel(x,8);%邻接区域个数
regionmat = zeros(num);%存储邻接矩阵
figure(1),imshow(L);
for i = 1:size(locations,1)
    text(locations(i,2),locations(i,1),num2str(i));
end
stats = regionprops(L);%构造区域个节点；
for i = 1:num
    node(i).id = i;
    node(i).data = stats(i);
    node(i).relaedge = [1 1 1];
    node(i).weight = [1 1 1];
end
[r c] = find(L==0);%遍历所有边界点
%遍历边界点的8邻域，构造区域邻接集合
[r1 c1] = size(L);
L_new = zeros(r1+2,c1+2);
L_new(2:(size(L_new,1)-1),2:(size(L_new,2)-1)) = L;
col = 8;%最多连通区域；
for i = 1:(size(r,1))
    temp = [L_new(r(i)-1+1,c(i)-1+1) L_new(r(i)-1+1,c(i)+1) L_new(r(i)-1+1,c(i)+1+1);
                    L_new(r(i)+1,c(i)-1+1) L_new(r(i)+1,c(i)+1) L_new(r(i)+1,c(i)+1+1);
                    L_new(r(i)+1+1,c(i)-1+1) L_new(r(i)+1+1,c(i)+1) L_new(r(i)+1+1,c(i)+1+1);];
    p = zeros(1,col);
    temp = unique(temp(:));
    temp = setdiff(temp,0);
    p(1:length(temp)) = temp;
    relas(i,:) = p;
end
relas = unique(relas,'rows');
%赋值邻接矩阵
for i = 1:size(relas,1)
    temp =setdiff(relas(i,:),0);
    cn = nchoosek(temp,2);
    for j =1:size(cn,1)
        regionmat(cn(j,1),cn(j,2)) = 1;
    end
end
disp('regionmat:');disp(regionmat);
%邻接矩阵图形化；
vertexs = [1 1;1 3;2 .5;2 3.5;3 1;3 3;3.5 2];
[r c] = find(regionmat ~= 0);
%绘制节点
figure(2),
axis([0 5 0 5]);
for i = 1:size(vertexs,1)
    hold on,plot(vertexs(i,1),vertexs(i,2),'go','markersize',10);
    hold on,plot(vertexs(i,1),vertexs(i,2),'gx','markersize',10);
    hold on,text(vertexs(i,1),vertexs(i,2)-0.2,num2str(i));
end
%绘制边界
for i = 1:size(r,1)
    line([vertexs(c(i),1) vertexs(r(i),1)],[vertexs(c(i),2) vertexs(r(i),2)],'color','r','linestyle','-.');
end

分割图的邻接矩阵仿真-1
邻接图示例.JPG
赞一下 回复此楼
2楼2015-10-27 21:44:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiyue6911 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 8rmuugja8q 2026-02-22 6/300 2026-02-23 06:39 by w4l55oybr1
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 4/200 2026-02-23 06:36 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:34 by w4l55oybr1
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 8/400 2026-02-23 06:21 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 10/500 2026-02-23 04:58 by 5jlh3qtdvx
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 6/300 2026-02-23 02:08 by 5jlh3qtdvx
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +5 successhe 2026-02-16 6/300 2026-02-22 21:38 by 山西悬空寺空悬�
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[教师之家] 为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 5/250 2026-02-22 21:15 by 山西悬空寺空悬�
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
信息提示
关闭
请填处理意见
关闭