How to fill interior as different color

Question

M.S. Khan 2019-9-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/480725-how-to-fill-interior-as-different-color

评论： M.S. Khan 2021-2-28

Hi Mathwork community,

Using below code, i can get boundary as blue color and gaps filled as yellow color but i want to fill the interior as green color.

Thanks in advance for all help.

s = load('Result_matrix.mat');
Aorg = s.Aorg;
gapmax = 10; 
[m,n] = size(Aorg);
A = double(Aorg>0); 
P=[1 2 3;
   8 0 4;
   7 6 5];
Apad = zeros(m+2,n+2);
maxattempt = 10;
for niter = 1:maxattempt  
    Apad(2:end-1,2:end-1) = A>0;
    B = zeros(m,n,8);
    for k=1:8
        [i,j] = find(P==k);
        B(:,:,k) = Apad(i-2+(2:end-1),j-2+(2:end-1));
    end
    As = A>0 & sum(diff(B(:,:,[1:8 1]),1,3)==1,3)<=1;
    [y,x] = find(As);
    if isempty(x)
        break
    end
    p = length(x);
    xy = [x,y];
    xy1 = reshape(xy,1,[],2);
    xy2 = reshape(xy,[],1,2);
    d = sum(abs(xy1-xy2),3);
    d(d==0) = NaN;
    i = (1:size(d,1))';
    [dmin,j] = min(d,[],2);
    keep = dmin <= gapmax;
    i = i(keep);
    j = j(keep);
    for k=1:length(i)
        xyi = xy(i(k),:);
        xyj = xy(j(k),:);
        dxy = xyi-xyj;
        if abs(dxy(1)) > abs(dxy(2))
            xr = linspace(xyi(1),xyj(1),abs(dxy(1))+1);
            yr = round(interp1([xyi(1),xyj(1)],[xyi(2),xyj(2)],xr));
        else
            yr = linspace(xyi(2),xyj(2),abs(dxy(2))+1);
            xr = round(interp1([xyi(2),xyj(2)],[xyi(1),xyj(1)],yr));
        end
        A(sub2ind(size(A),yr(2:end-1),xr(2:end-1))) = 3; 
    end
end

4 个评论
显示 2更早的评论隐藏 2更早的评论

KALYAN ACHARJYA 2019-9-17

Thanks @Bruno

@M.S. Khan

Are you looking for following one? Bruno provided the answer, please do confirm.

M.S. Khan 2019-9-17

Yes Mr. Kalyan. i want background as blue, boundary as yellow, filles holes as red and interior as green.

Thanks for cooperation.

请先登录，再进行评论。

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

Answer 1

Bruno Luong 2019-9-17

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/480725-how-to-fill-interior-as-different-color#answer_392214

编辑：Bruno Luong 2019-9-17

在 MATLAB Online 中打开

Fill the hole with imfill

Next you can color your integer image by defining the appropriate colormap

% Generate test image contains values in 0:4
[~,~,I]=histcounts(peaks,5);
I = I-1;
imagesc(I)
colormap([1 0 0; % color of 0 value
          0 1 0; % color of 1 value
          0 0 1; % color of 2 value
          0 0 0; % color of 3 value
          1 1 1])% color of 4 value
      
set(gca,'Clim',[-0.5 max(I(:))+0.5]);
colorbar

16 个评论
显示 14更早的评论隐藏 14更早的评论

Bruno Luong 2019-9-17

编辑：Bruno Luong 2019-9-17

在 MATLAB Online 中打开

I already gave you a clear guidance about coloring. You must use IMFILL and then apply my code after filling.

I have done with my own function ROSENFELD that does more or less IMFILL

s = load('Result_matrix.mat');
% Code to fill the gap, just recopy here but IRRELEVANT to coloring
Aorg = s.Aorg;
gapmax = 10; 
[m,n] = size(Aorg);
A = double(Aorg>0); 
P=[1 2 3;
   8 0 4;
   7 6 5];
Apad = zeros(m+2,n+2);
maxattempt = 10;
for niter = 1:maxattempt  
    Apad(2:end-1,2:end-1) = A>0;
    B = zeros(m,n,8);
    for k=1:8
        [i,j] = find(P==k);
        B(:,:,k) = Apad(i-2+(2:end-1),j-2+(2:end-1));
    end
    As = A>0 & sum(diff(B(:,:,[1:8 1]),1,3)==1,3)<=1;
    [y,x] = find(As);
    if isempty(x)
        break
    end
    p = length(x);
    xy = [x,y];
    xy1 = reshape(xy,1,[],2);
    xy2 = reshape(xy,[],1,2);
    d = sum(abs(xy1-xy2),3);
    d(d==0) = NaN;
    i = (1:size(d,1))';
    [dmin,j] = min(d,[],2);
    keep = dmin <= gapmax;
    i = i(keep);
    j = j(keep);
    for k=1:length(i)
        xyi = xy(i(k),:);
        xyj = xy(j(k),:);
        dxy = xyi-xyj;
        if abs(dxy(1)) > abs(dxy(2))
            xr = linspace(xyi(1),xyj(1),abs(dxy(1))+1);
            yr = round(interp1([xyi(1),xyj(1)],[xyi(2),xyj(2)],xr));
        else
            yr = linspace(xyi(2),xyj(2),abs(dxy(2))+1);
            xr = round(interp1([xyi(2),xyj(2)],[xyi(1),xyj(1)],yr));
        end
        A(sub2ind(size(A),yr(2:end-1),xr(2:end-1))) = 3; 
    end
end
I = rosenfeld(~(A>0)); % replace this with IMFILL if you own IMP toolbox, I don't
% Coloring part
imagesc(I)
colormap([0 1 0; % yellow, color of 1 value
          0 0 1; % blue, color of 2 value
          1 0 0])% red, color of 3 value
      
set(gca,'Clim',[min(I(:))-0.5 max(I(:))+0.5]);
h = colorbar;
set(h, 'Ytick', 1:3, 'YtickLabel', {'Boundary','Exterior','Interior'});

M.S. Khan 2019-9-18

在 MATLAB Online 中打开

Hi Bruno, the same file i using the same method but i could not succeed. can you guide me where i am making mistake. i used fill function but not working. The last part of the coding is making problem. Thanks

gapmax = 10; 
[m,n] = size(Aorg);
A = double(Aorg>0); 
P=[1 2 3;
   8 0 4;
   7 6 5];
Apad = zeros(m+2,n+2);
maxattempt = 10;
for niter = 1:maxattempt  
    Apad(2:end-1,2:end-1) = A>0;
    B = zeros(m,n,8);
    for k=1:8
        [i,j] = find(P==k);
        B(:,:,k) = Apad(i-2+(2:end-1),j-2+(2:end-1));
    end
    As = A>0 & sum(diff(B(:,:,[1:8 1]),1,3)==1,3)<=1;
    [y,x] = find(As);
    if isempty(x)
        break
    end
    p = length(x);
    xy = [x,y];
    xy1 = reshape(xy,1,[],2);
    xy2 = reshape(xy,[],1,2);
    d = sum(abs(xy1-xy2),3);
    d(d==0) = NaN;
    i = (1:size(d,1))';
    [dmin,j] = min(d,[],2);
    keep = dmin <= gapmax;
    i = i(keep);
    j = j(keep);
    for k=1:length(i)
        xyi = xy(i(k),:);
        xyj = xy(j(k),:);
        dxy = xyi-xyj;
        if abs(dxy(1)) > abs(dxy(2))
            xr = linspace(xyi(1),xyj(1),abs(dxy(1))+1);
            yr = round(interp1([xyi(1),xyj(1)],[xyi(2),xyj(2)],xr));
        else
            yr = linspace(xyi(2),xyj(2),abs(dxy(2))+1);
            xr = round(interp1([xyi(2),xyj(2)],[xyi(1),xyj(1)],yr));
        end
        A(sub2ind(size(A),yr(2:end-1),xr(2:end-1))) = 3; 
    end
end
if isempty(x)
    fprintf('Success\n')
else
    fprintf('Fail\n')
end
    
figure(1);
clf;
ax=subplot(1,3,1);
imagesc([-20 280],[-140 140],Aorg>0); 
%----
%----
set(ax,'xtick',-20:10:280,'ytick',-140:10:140,'clim',[0 3]); % change here for Q1
axis('on','image') 
axis equal
% title('Apple 10,000 point cloud @ width 10');
ax=subplot(1,3,2);
imagesc([-20 280],[-140 140],A); % change here for Q1
set(ax,'xtick',-20:10:280,'ytick',-140:10:140,'clim',[0 3]); % change here for Q1
axis('on','image') 
axis equal
title('Apple 10,000 point cloud after filling gaps');
dlmwrite('Boundary_closed_1s_3s.txt',A)
% new attempt to pun 3rd figure
%% interior filling
A = load('Boundary_closed_1s_3s.txt')
Perimeter=bwperim(A)
% dlmwrite('Closed Boundary filled with 1s.txt',Perimeter)
interior_filled = imfill(A,'holes')
A = 3*double(A)
interior_filled = 2*double(interior_filled)
% interior_filled(Perimeter)=1 
interior_filled(Perimeter)=A(Perimeter)
interior_filled(interior_filled==3)=1
interior_filled(interior_filled==9)=3
% imshow(interior_filled,[])
dlmwrite('Boundary_closed_1s_3s_interior_2s.txt',interior_filled)
%% to count number of filled voxels on bounary and interior 
m = load('Boundary_closed_1s_3s_interior_2s.txt')
%----- filling colors
ax=subplot(1,3,3); %  i want to fill 
I = imfill(~(m>0),'holes'); % replace this with IMFILL if you own IMP toolbox, I don't
% Coloring part
imagesc(I)
colormap([0 1 0; % yellow, color of 1 value
          0 0 1; % blue, color of 2 value
          1 0 0])% red, color of 3 value
      
set(gca,'Clim',[min(I(:))-0.5 max(I(:))+0.5]);
h = colorbar;
set(h, 'Ytick', 1:3, 'YtickLabel', {'Boundary','Exterior','Interior'});