How to measure the average thickness at the center region of binary image

Question

Timilehin Oyinloye 2023-3-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1932385-how-to-measure-the-average-thickness-at-the-center-region-of-binary-image

评论： Timilehin Oyinloye 2023-3-22

I have a binary image that has various objects (white pixels) in various positions. I want to determine the average thickness of each white object in the binary image at its center (as shown in the third image with red letters). I appreciate your support.

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

请先登录，再进行评论。

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

Answer 1

Antoni Garcia-Herreros 2023-3-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1932385-how-to-measure-the-average-thickness-at-the-center-region-of-binary-image#answer_1197665

在 MATLAB Online 中打开

Hello Timilehin,

You can refer to this thread for more info, but this should do the trick.

lear all 
close all
I=imread("image.jpeg");
if length(size(I))>2
    I=rgb2gray(I);
end
I=imbinarize(I);
I = imfill(I, 'holes');
BW2 = bwareaopen(I, 10); % Retain only large objects
r=regionprops(BW2,'Area','Centroid');
ApproxThick=zeros(size(r));
MeanCentralThick=ApproxThick;
Icopy=I;
for i=1:length(r) % Loop through all your objects
    bwi = bwareafilt(Icopy, 1);
    Icopy=logical(Icopy-bwi); % Subtract the objects from the image
    
    ri=regionprops(bwi,'Area','Centroid');
    BWSkel=bwskel(bwi);
    ApproxThick(i)=ri.Area/sum(BWSkel(:)); % This would give you a good idea of the mean thickness of all the object
    
    connec=bwconncomp(BWSkel); 
    [row,col] = ind2sub(connec.ImageSize,cell2mat(connec.PixelIdxList));
    Mask=false(size(BWSkel)); % Mask of the skeletonize image containing only the central region
    idxmin=floor(sum(BWSkel(:))*0.4); % Adjust these values to obtain the central region
    idxmax=floor(sum(BWSkel(:))*0.6);
    Mask(row(idxmin:idxmax),col(idxmin:idxmax))=BWSkel(row(idxmin:idxmax),col(idxmin:idxmax));
    
    EuclImage = bwdist(~bwi); % Euclidean distance
    Thickness=EuclImage(Mask);
    MeanCentralThick(i)=2*mean(Thickness); % Average thickness of the central region
end