Can anyone help me to segment the stone from the kidney image

Question

MINO GEORGE 2021-6-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/862085-can-anyone-help-me-to-segment-the-stone-from-the-kidney-image

评论： Image Analyst 2022-4-28

st1.PNG

在 MATLAB Online 中打开

Hi

I am using k means clustering algorithm and this is code i have used

grayLevels = double(grayImage(:)); 
[clusterIndexes, clusterCenters] = kmeans(grayLevels, numberOfClusters,...
	'distance', 'sqEuclidean', ...
	'Replicates', 2);
labeledImage = reshape(clusterIndexes, rows, columns);
subplot(2, 2, 2);
 figure,imshow(labeledImage,[]);
 caption = sprintf('k-means with %d clusters', numberOfClusters);
 title(caption,'FontSize', fontSize);
 axis on image;
 colorbar
hp = impixelinfo();

And the input and output image are attached here. Can anyone help me to segment the stone region(the white dot) from the image.

Any help is aapreciated.

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

请先登录，再进行评论。

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

Answer 1

Saurabh Gund 2022-4-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/862085-can-anyone-help-me-to-segment-the-stone-from-the-kidney-image#answer_951870

在 MATLAB Online 中打开

clc
clear all
close all
warning off
[filename,pathname]=uigetfile('*.*','Pick a MATLAB code file');
filename=strcat(pathname,filename);
a=imread(filename);
%subplot(3,4,1),
imshow(a);
%Convert to grayscale image
%b=rgb2gray(a);
b=im2gray(a);
%subplot(3,4,2),
imshow(b);
impixelinfo;
c=imbinarize(b,20/255);
%c=b>50;
%subplot(3,4,3),
imshow(c);
%Fill the holes in the image
d=imfill(c,'holes');
%subplot(3,4,4),
imshow(d);
%delete writings
e=bwareaopen(d,1000);
%subplot(3,4,5),
imshow(e);
PreprocessedImage=uint8(double(a).*repmat(e,[1 1 3]));
%subplot(3,4,6),
imshow(PreprocessedImage);

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

Image Analyst 2022-4-28

在 MATLAB Online 中打开

st1.PNG

Here is what your code produces. Where is the kidney?

clc
clear all
close all
warning off
% [filename,pathname]=uigetfile('*.*','Pick a MATLAB code file');
filename = '\st1.png';
pathname = pwd;
filename=strcat(pathname,filename);
a=imread(filename);
%subplot(3,4,1),
imshow(a);
%Convert to grayscale image
%b=rgb2gray(a);
b=im2gray(a);
%subplot(3,4,2),
imshow(b);
impixelinfo;
c=imbinarize(b,20/255);
%c=b>50;
%subplot(3,4,3),
imshow(c);
%Fill the holes in the image
d=imfill(c,'holes');
%subplot(3,4,4),
imshow(d);
%delete writings
e=bwareaopen(d,1000);
%subplot(3,4,5),
imshow(e);
PreprocessedImage=uint8(double(a).*repmat(e,[1 1 3]));
%subplot(3,4,6),
imshow(PreprocessedImage);