Image segementation of cement paste complex structure!
13 次查看(过去 30 天)
显示 更早的评论
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.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/153576/image.jpeg)
Joanna
2 个评论
Maria Huguett
2019-4-5
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!
采纳的回答
Image Analyst
2013-12-9
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 个评论
Sean de Wolski
2013-12-13
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...
更多回答(1 个)
Sean de Wolski
2013-12-9
编辑:Sean de Wolski
2013-12-9
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.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/174056/image.png)
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'});
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/174061/image.png)
7 个评论
Sean de Wolski
2015-3-5
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
2015-3-5
To get rid of holes in blobs smaller than a certain amount, say 500 pixels, I use bwareaopen
binaryImage = ~bwareaopen(~binaryImage, 500);
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!