Skull striping without affecting tumor region

9 次查看（过去 30 天）
Sehrish
Sehrish 2012-7-5
评论： elif yoldas 2016-4-15
Can anyone help me to remove the skull part of the brain image without effecting the tumor region? the image is http://imageshack.us/photo/my-images/196/img6oe.jpg/ Actually the intensity level of tumor as well as the skull part is really high and whatever i tried, it effects the tumor region as well. kindly help me regarding this:

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

回答（3 个）

Matt Kindig
Matt Kindig 2012-7-5
编辑：Matt Kindig 2012-7-5
If you have the Image Processing Toolbox, you can use 'regionprops' to identify the various objects after your thresholding. For your particular image, I was able to identify the tumor using the following approach:
img = imread('/path/to/your/filename.jpg');
bw = (img > 0.5*255); %thresholded to remove background and white/gray matter
lbl = bwlabel(bw); %labelled regions
props = regionprops(lbl, 'Solidity', 'Area');
% solidity is the percentage "filled" of an area. For the skull,
%the solidity will be really low.
solidity = [props.Solidity];
area = [ props.Area];
hiSolid= solidity > 0.2; %get only high solidity objects
maxArea = max( area(hiSolid));
tumorLabel = find( area==maxArea); %label of tumor
tumor = ismember(lbl, tumorLabel); %b/w image of tumor
imshow(tumor); %this isolates tumor
  3 个评论
显示 1更早的评论
Matt Kindig
Matt Kindig 2012-7-5
编辑：Matt Kindig 2012-7-5
Sure, one way is to first extract everything inside the skull, and then remove everything that is the skull itself. Something like this seems to work. Again, not necessarily a general solution, but it should give you an overview of some common image processing functions that can help.
insideSkull = (img > 0.3*255); %isolate all dense material
insideSkull = imfill(insideSkull, 'holes');
bw = (img > 0.4*255);
lbl = bwlabel(skull);
skull = find(Solidity < 0.2);
skull = ismember(lbl, skull);
skullOrBorder = skull | ~insideSkull;
noSkull = img;
noSkull(skullOrBorder)=0;
imshow(noSkull); %this only contains gray matter and tumor
You might need to adjust the thresholds to get better performance. Also, I would recommend looking at the Demos for the Image Processing Toolbox. You can access them from the Help browser using Help->Image Processing Toolbox->Demos.
Sehrish
Sehrish 2012-7-5
its working but giving different results for different images. the problem is in threshhold values. How can i get automated threshholds?

请先登录，再进行评论。

Felix
Felix 2012-7-5
编辑：Felix 2012-7-5
Unless you have a huge number of images, I think the best solution is to mask it manually, because no algorithm is as good as the eye. For example I use the program MIPAV, which can read many different scanner formats.
Anton Semechko
Anton Semechko 2012-7-5
编辑：Anton Semechko 2012-7-5
The simplest way to segment the tumor in the sample image would be to use region growing algorithm. There are a few implementations of this method on FEX, but for the following demo I used this one:
http://www.mathworks.com/matlabcentral/fileexchange/19084-region-growing
1) Download the function titled 'regiongrowing' from the link given above into your working directory.
2) Segment the tumor using pixel position [380 480] as the seed point:
% get the sample image
im=imread('http://desmond.imageshack.us/Himg196/scaled.php?server=196&filename=img6oe.jpg&res=landing');
% segment the ROI
thr=15; % intensity similarity tolerance
bw=regiongrowing(double(im),380,480,thr);
% fill in the holes
bw=imfill(bw,'holes');
% visualize the segmentation
r=im; g=im; b=im;
r(bw)=255; g(bw)=255;
rgb=cat(3,r,g,b);
figure, imshow(rgb)
clear r g b
You will get different results depending on the value of 'thr' parameter.
Hope this helps.
  9 个评论
显示 7更早的评论
Image Analyst
Image Analyst 2013-1-30
You have to adapt it. That's what image processing is all about. Often people make an algorithm that works fine for just one image but is not robust enough for all similar images. That's where you skill as an image analyst comes in. You have to figure out what parameters need to be veried for the other images, for example, threshold values, minimum sizes, kernel width, etc. It's your project, and you have your images, so it's up to you, and perhaps the radiologist who's helping you that supplied you with the images.
elif yoldas
elif yoldas 2016-4-15
code is not done work and image is not found :(

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Image Data Workflows 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by