Fatal编程技术网
  • matlab/
  • Matlab 如何使用邻接矩阵创建图形

Matlab 如何使用邻接矩阵创建图形

matlab image-processing graph

Matlab 如何使用邻接矩阵创建图形,matlab,image-processing,graph,Matlab,Image Processing,Graph,我从图像中找到邻接，然后尝试使用graph函数绘制图形： g = graph(adj); plot(g); 但我没有得到预期的图表 输入图像 这是我的密码。我将邻接矩阵存储在变量adj中 clc; clear all; close all; X=imread('one.bmp'); imfinfo('one.bmp') figure,imshow(X) b = imresize(X,[100,100]); si = size(b,1); sj = size(b,2); %figure;

我从图像中找到邻接，然后尝试使用graph函数绘制图形：

g = graph(adj); 
plot(g);
但我没有得到预期的图表

输入图像

这是我的密码。我将邻接矩阵存储在变量
adj
中

clc;
clear all;
close all;
X=imread('one.bmp');
imfinfo('one.bmp')
figure,imshow(X)

b = imresize(X,[100,100]);
si = size(b,1);
sj = size(b,2);
%figure;imshow(b);

% Binarization
% th = graythresh(b);
% I = im2bw(b,th);
%thinning
kl=bwmorph(~b,'thin',inf);
figure,imshow(kl)
R(:,:)=kl(:,:);
%grid size 
t1=25;
D=100;
 I=1;
U1=t1;
J=1;
U2=t1;
E=1;
t2=D/t1;
%Z=1;
for iir=1:t2
    for jjr=1:t2
 B(I:U1,J:U2)=R(I:U1,J:U2);
%  vc=sum(B(I:U1,J:U2));
%  Fd=sum(vc);
 [x,y]=find(B==1);
 CX=mean(x);
 CY=mean(y);
 CXXX(E)=CX;
 CYYY(E)=CY;
 CXX(iir,jjr)=CX;
 CYY(iir,jjr)=CY;

 T(I:U1,J:U2)=B(I:U1,J:U2);
    J=J+t1;
  U2=U2+t1;
E=E+1;
 clear B x y 

    end

I=I+t1;
U1=U1+t1;
J=1;
  U2=t1;

end
%plot and grid
figure,imshow(R)
hold on
M10 = size(R,1);
N10 = size(R,2);

a=t1; 
b=t1;
for k = 1:a:M10
    x = [1 N10]; 
   y = [k k]; 
   plot(x,y,'Color','white');
    set(findobj('Tag','MyGrid'),'Visible','on')
end
for k = 1:b:N10 
    x = [k k]; 
    y = [1 M10];
    plot(x,y,'Color','white');
    set(findobj('Tag','MyGrid'),'Visible','on')
end

plot(CYY,CXX,'g*')
%line(CYY,CXX)
%CC=bwconncomp(CXX,4)

hold off
%node neighbourhoood analyssis
N1=t2;
for I2=1:N1
    for J2=1:N1
        %last row
       if(I2>=N1)
             W1=CXX(I2,J2);
           W2=CXX(I2-1,J2);
           W3=CYY(I2,J2);
           W4=CYY(I2-1,J2);    
           W6=[W1,W2];
           W7=[W3,W4];
           line(W7,W6);
           if(J2>=N1)
               Z=CXX(I2,J2);
           else
           if (CXX(I2,J2+1)>1)&& ((CYY(I2,J2+1)>1))
            TXX=CXX(I2,J2);
            TYY=CXX(I2,J2+1);
               TTX=CYY(I2,J2);
               TTY=CYY(I2,J2+1);

            IY=[TXX,TYY];
            IIY=[TTX,TTY];
            line(IIY,IY);
           end
           end


       else
           if(J2>=N1);
             W1=CXX(I2,J2);
           W2=CXX(I2+1,J2);
           W3=CYY(I2,J2);
           W4=CYY(I2+1,J2);    
           W6=[W1,W2];
           W7=[W3,W4];
           line(W7,W6);
           else   
               if (CXX(I2,J2+1)>1)&& ((CYY(I2,J2+1)>1))
            TXX=CXX(I2,J2);
            TYY=CXX(I2,J2+1);
               TTX=CYY(I2,J2);
               TTY=CYY(I2,J2+1);

            IY=[TXX,TYY];
            IIY=[TTX,TTY];
            line(IIY,IY);


               end
         if (CXX(I2+1,J2)>1)&& ((CYY(I2+1,J2)>1))
           W1=CXX(I2,J2);
           W2=CXX(I2+1,J2);
           W3=CYY(I2,J2);
           W4=CYY(I2+1,J2);    
           W6=[W1,W2];
           W7=[W3,W4];
           line(W7,W6);
           J2=J2+1

         end
        end
    end
    end
end

    A=zeros(t2,t2);   
    ttt=1;
   for rt=1:t2
      for  rt1=1:t2
          if(CXX(rt,rt1)>1)
              A(rt,rt1)=ttt
          end
          ttt=ttt+1;
      end
   end



g=1;
jk=1;
um=t2-1;
um1=t2;
for iir=1:t2
    for jjr=1:t2
        if(A(iir,jjr)>=0)
            BB(jk)=0;
                       DD(g)=0;
                      FF(g)=0;
                      HH(g)=0;
        end
        if(A(iir,jjr)>=1)
          if(iir==um1)&&(jjr==1)
             GG(g)=A(iir,jjr);
            HH(g)=A(iir-1,jjr); 
            BB(jk)=A(iir,jjr+1);
            DD(g)=0;
          else
              if(iir==um1)&&(jjr>1)&&(jjr<=um)
               FF(g)=A(iir,jjr-1); 
                BB(jk)=A(iir,jjr+1);
                HH(g)=A(iir-1,jjr);
              else

              if(iir==um1)&&(jjr==um1)
                  HH(g)=A(iir-1,jjr);
                  FF(g)=A(iir,jjr-1);
                  DD(g)=0;
                  BB(jk)=0;
              else

        if(iir==1)&&(jjr==um1)
            FF(g)=A(iir,jjr-1);
            DD(g)=A(iir+1,jjr);
            BB(jk)=0;

        else
            if(iir>=1)&&(iir<=um)&&(jjr==um1)
             HH(g)=A(iir-1,jjr);
             DD(g)=A(iir+1,jjr);
            FF(g)=A(iir,jjr-1);
            BB(jk)=0;
            else
              if(iir==1)&&(jjr==1)
                BB(jk)=A(iir,jjr+1);  
                DD(g)=A(iir+1,jjr);
              else
              if(iir==1)&&(jjr>=1)&&(jjr<=um)  
              FF(g)=A(iir,jjr-1);    
              DD(g)=A(iir+1,jjr); 
                BB(jk)=A(iir,jjr+1);
              else
               if(iir>1)&&(iir<=um)&&(jjr==1)   
                HH(g)=A(iir-1,jjr);
                DD(g)=A(iir+1,jjr);
                BB(jk)=A(iir,jjr+1);
               else
                   if(iir>1)&&(iir<=um)&&(jjr>1)&&(jjr<=um) 
                       BB(jk)=A(iir,jjr+1);
                   DD(g)=A(iir+1,jjr);
                   HH(g)=A(iir-1,jjr);
                    FF(g)=A(iir,jjr-1);
                   end
               end
              end
              end
            end
        end
              end
          end
          end
        end
    g=g+1;
    jk=jk+1;
    end

end
adj=zeros(t2*t2);
H9=size(adj);
Y=1;
for ll=1:1
        for ii=1:H9(1,1)
            for jj=1:H9(1,1)
                if (ii>=4)
                    if (jj==DD(1,Y))
                    adj(ii,jj)=1;
                end
                if (jj==FF(1,Y))
                    adj(ii,jj)=1;
                end
                if (jj==BB(1,Y))
                    adj(ii,jj)=1;
                end
                if (jj==HH(1,Y))
                    adj(ii,jj)=1;
                end
                else
                if (jj==BB(1,Y))
                    adj(ii,jj)=1;
                end
                if (jj==DD(1,Y))
                    adj(ii,jj)=1;
                end
                if (jj==FF(1,Y))
                    adj(ii,jj)=1;
                end
                if (jj==HH(1,Y))
                    adj(ii,jj)=1;
                end
                end
            end
            Y=Y+1;
        end

end
问题是graph对象不知道要绘制的点的位置。需要使用质心坐标来确定从邻接矩阵绘制边的位置

假设邻接矩阵中点的质心由
CX
和
CY
给出，则需要将此信息传递给
绘图
函数

g = graph(adj); 
plot(g, 'XData', CX, 'YData', CY);
请参见

如何获得邻接矩阵？它应该代表什么？像素邻接？没有足够的条目来表示所有像素。连通分量邻接？您的输入图像看起来已完全连接，因此条目太多。你能给我这个程序的完整代码吗？这是完整代码。我唯一遗漏的是将邻接矩阵赋给变量
A
。您只需将
CX
和
CY
变量替换为每个图形节点对应的坐标即可。请务必代我向shivu和raghu问好。：）我不知道如何用对应于每个图形节点的坐标替换CX和CY变量你已经计算出了质心，对吗？然后在
adj
中确定这些质心之间的连通性。因此，您应该知道
adj
中的哪个节点与每个质心关联。 
g = graph(adj); 
plot(g, 'XData', CX, 'YData', CY);