how to calculate the area binary if some images have different pattern

Question

mohd akmal masud 2021-11-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1577570-how-to-calculate-the-area-binary-if-some-images-have-different-pattern

评论： mohd akmal masud 2021-11-4

采纳的回答： yanqi liu

在 MATLAB Online 中打开

Hi all I have 4 images that all images have 4 area(like circle) as below.

Then, I used this coding below to get the area for all the 4 images that every images have 4 circle

%% TO GET THE VOLUME SEGMENTATION AFTER DEEP LEARNING PERFORM(blob all the sphere binary)
alldice=[]
acc=[]
Ts = [];
Data = [];
for ii=1:4
    subplot(1,4,ii)
    I = readimage(IMDS,ii);
    [C,scores] = semanticseg(I,net1);
    
    outt2=C=="foreground";
    st2=strel('disk',5);
    outt22=imerode(outt2,st2);
%       outt22(:,:,ii) = imerode(outt2,st2);
    title('output')
    imshow(outt22)
info = regionprops(outt22,'Boundingbox') ;
imshow(outt22)
hold on
allbb=[]
for k = 1 : length(info)
     BB = info(k).BoundingBox;
     rectangle('Position', [BB(1),BB(2),BB(3),BB(4)],'EdgeColor','r','LineWidth',2) ;
     allbb=[allbb; BB]
end
	vol = [];
	for i=1:4
    
		zz=imcrop(outt22,allbb(i,:));
		
		stats = regionprops('table',zz,'Centroid',...
			'MajorAxisLength','MinorAxisLength');
		T = regionprops('table', outt22,'Area','Centroid')
		m1=stats.MajorAxisLength;
		m2=stats.MinorAxisLength;
		v1=stats.MinorAxisLength;
		
		vol(i)=m1*m2*(v1/2);					
    
	end
	ss=sprintf('the volume %d is %6f,%d is %6f,%d is %6f,%d is %6f',1,vol(1),2,vol(2),3,vol(3),4,vol(4));
		
	msgbox(ss);
end

Its work. And the answer like below

T =
4×2 table
Area        Centroid    
____    ________________
15     105.73    132.53
21     126.52    145.38
64     147.12    108.86
31     147.42     133.9
T =
4×2 table
Area        Centroid    
____    ________________
15     105.73    132.53
21     126.52    145.38
64     147.12    108.86
31     147.42     133.9
T =
4×2 table
Area        Centroid    
____    ________________
15     105.73    132.53
21     126.52    145.38
64     147.12    108.86
31     147.42     133.9
T =
4×2 table
Area        Centroid    
____    ________________
15     105.73    132.53
21     126.52    145.38
64     147.12    108.86
31     147.42     133.9

But when I tried the same code into this images below. It have error.

**NOTE: the first and second image just have 3 circle, the 2 others have 4 circle.

T =
3×2 table
Area        Centroid    
____    ________________
6     125.67    146.33
92     147.77    109.37
52     147.15    134.06
T =
3×2 table
Area        Centroid    
____    ________________
6     125.67    146.33
92     147.77    109.37
52     147.15    134.06
T =
3×2 table
Area        Centroid    
____    ________________
6     125.67    146.33
92     147.77    109.37
52     147.15    134.06
Index in position 1 exceeds array bounds (must not exceed 3).

anyone can help me?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

yanqi liu 2021-11-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1577570-how-to-calculate-the-area-binary-if-some-images-have-different-pattern#answer_822780

在 MATLAB Online 中打开

%% TO GET THE VOLUME SEGMENTATION AFTER DEEP LEARNING PERFORM(blob all the sphere binary)
alldice=[]
acc=[]
Ts = [];
Data = [];
for ii=1:4
    subplot(1,4,ii)
    I = readimage(IMDS,ii);
    [C,scores] = semanticseg(I,net1);
    
    outt2=C=="foreground";
    st2=strel('disk',5);
    outt22=imerode(outt2,st2);
    %       outt22(:,:,ii) = imerode(outt2,st2);
    
    title('output')
    imshow(outt22)
    
    info = regionprops(outt22,'Boundingbox') ;
    imshow(outt22)
    hold on
    
    allbb=[]
    for k = 1 : length(info)
        BB = info(k).BoundingBox;
        rectangle('Position', [BB(1),BB(2),BB(3),BB(4)],'EdgeColor','r','LineWidth',2) ;
        allbb=[allbb; BB]
    end
    
    vol = [];
    for i=1:size(allbb,1)
        
        zz=imcrop(outt22,round(allbb(i,:)));
        
        stats = regionprops('table',zz,'Centroid',...
            'MajorAxisLength','MinorAxisLength');
        T = regionprops('table', outt22,'Area','Centroid')
        m1=stats.MajorAxisLength;
        m2=stats.MinorAxisLength;
        v1=stats.MinorAxisLength;
        
        vol(i)=m1*m2*(v1/2);
        
    end
    ss=sprintf('the volume %d is %6f,%d is %6f,%d is %6f,%d is %6f',1,vol(1),2,vol(2),3,vol(3),4,vol(4));
    
    msgbox(ss);
end

3 个评论
显示 1更早的评论隐藏 1更早的评论

yanqi liu 2021-11-4

在 MATLAB Online 中打开

%% TO GET THE VOLUME SEGMENTATION AFTER DEEP LEARNING PERFORM(blob all the sphere binary)
alldice=[]
acc=[]
Ts = [];
Data = [];
for ii=1:4
    subplot(1,4,ii)
    I = readimage(IMDS,ii);
    [C,scores] = semanticseg(I,net1);
    
    outt2=C=="foreground";
    st2=strel('disk',5);
    outt22=imerode(outt2,st2);
    %       outt22(:,:,ii) = imerode(outt2,st2);
    
    title('output')
    imshow(outt22)
    
    info = regionprops(outt22,'Boundingbox') ;
    imshow(outt22)
    hold on
    
    allbb=[]
    for k = 1 : length(info)
        BB = info(k).BoundingBox;
        rectangle('Position', [BB(1),BB(2),BB(3),BB(4)],'EdgeColor','r','LineWidth',2) ;
        allbb=[allbb; BB]
    end
    
    vol = [];
    for i=1:size(allbb,1)
        
        zz=imcrop(outt22,round(allbb(i,:)));
        
        stats = regionprops('table',zz,'Centroid',...
            'MajorAxisLength','MinorAxisLength');
        T = regionprops('table', outt22,'Area','Centroid')
        m1=stats.MajorAxisLength;
        m2=stats.MinorAxisLength;
        v1=stats.MinorAxisLength;
        
        vol(i)=m1*m2*(v1/2);
        
    end
    ss = 'the volume';
    for i = 1 : length(vol)
        ss=sprintf('%s %d is %6f ', ss, i, vol(i));
    end
    msgbox(ss);
end