Image segementation of cement paste complex structure!
7 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Dear all,
I have a cement paste gray image of the CT, now I want to do some quantative analysis of the porosity and the microstructure characterization. From the image we at least can get 4 main composition, how can I use the histogram or threshold to finely distinguish different parts in this image? Can anybody help me with the Matlab programming? I tried the standard segmentation but it doesn´t work so well since the boundary of cement material is almost invisible. Thank you so much for your help.
Joanna
2 commentaires
Maria Huguett
le 5 Avr 2019
Could you tell me please, where can I find the "Sean's intensity thresholding method"? I would like to try it, since I need a method to segmentate pores, grains and paraffin of my sediment sample, to get porosity values. Thanks!
Réponse acceptée
Image Analyst
le 9 Déc 2013
I'm not sure what you meant by "the boundary of cement material is almost invisible" - perhaps you could explain.
If Sean's intensity thresholding method doesn't work, like say you had some structures that were of the same brightness but just different textures, you could combine that with some texture segmentation. Try stdfilt() or entropyfilt(). Or you could try simulated annealing (I have a demo fro that if you want). Again, only for cases where intensity thresholding doesn't do the job well enough.
If you want an intensity thresholding segmentation demo, see my Image Segmentation Tutorial: http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862
15 commentaires
Sean de Wolski
le 13 Déc 2013
I guess I'm having a hard time understanding why there isn't more contrast and why the background (air) isn't the darkest part...
Plus de réponses (1)
Sean de Wolski
le 9 Déc 2013
Modifié(e) : Sean de Wolski
le 9 Déc 2013
There are a few things a little weird about your image:
- First, for a CT image, I'd expect the background to have a lower intensity than everything else since this is void space (unless is was scanned in a liquid of some kind). What do you know about the setup of the system and do you have control over it?
- Second, it appears the histogram has already been modified in some way. I would expect this to be fairly bimodal for an object in the forground.
As far as identifying the cement matrix, here's a rough first pass to get you started:
I = imread('cct.jpg');
mask of the concrete part
disk = getnhood(strel('disk',4)); %disk
Istd = stdfilt(I,disk); %std filter tor emove backgoryund
Mcyl = Istd>2;
C = (conv2(double(I),double(disk),'same'));
levels = multithresh(C,2); %otsu thresholding
L = imquantize(C,levels); %apply thresholds
L = L+1; %increment
L(~Mcyl) = 0; %remove background
%%visualize
cmap = [1 1 1;lines(3)]; % colormap
Lrgb = label2rgb(L,cmap); % to view
% view it
imshow(Lrgb);
colormap(cmap);
hCb = colorbar;
set(hCb,'YTick',0:3);
set(hCb,'YTickLabel',{'Background','Porous','Cement Matrix','Aggregate'});
7 commentaires
Sean de Wolski
le 5 Mar 2015
A morphological opening ( imopen) with a small structuring element should do it. Alternatively, you could create a mask of the pores and then use imfill(pores,'holes') to fill in all of the holes inside a pore.
Image Analyst
le 5 Mar 2015
To get rid of holes in blobs smaller than a certain amount, say 500 pixels, I use bwareaopen
binaryImage = ~bwareaopen(~binaryImage, 500);
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!