MATLAB: How can i detect/localize the spine region and segment individual vertebra from the given image. I really need help to start with. With much thanks.

filterslocalizesegmentspinevertebra

I intend to use the segmented vertebra to further classify it as benign or malignant. Since many regions in the image have similar intensities i am finding it difficult to precisely segment vertebrae. could any shape features me used? plz suggest.

Best Answer

Find algorithms Here in this link 20.7.3.3 Medical Applications, Vertebra, Spine, Skeletons, Bone

Nothing is built in to MATLAB - no "Vertebrae Segmentation/isolation Toolbox" exists - so you'll have to build something up from published algorithms and lower level toolbox functions.

Of course if this is for publication or a university degree, you can't just use a published algorithm, you'll have to show how your new and totally different algorithm, or improved, adapted algorithm, is better than existing ones.

Related Solutions

MATLAB: How can i use region based segmentation to segment common carotid artery of ultrasound image and measure intima media thickness

Find an algorithm here: http://iris.usc.edu/Vision-Notes/bibliography/contentsmedical.html#Medical%20Applications,%20CAT,%20MRI,%20Ultrasound,%20Heart%20Models,%20Brain%20Models. Return here if you have any programming questions.

MATLAB: How to find the coordinates of the top and bottom end plate and calculate the height

The FEX submission pgonCorners (Download) can be used to find approximate corners for each of the vertebrae. Once you have the corner points, you can draw lines through the ones with similar vertical coordinates to demarcate the end plates. You can also average the vertical coordinates of the 4 vertices to define an overall height for each vertebra.

load spine
[m,n]=size(BW);
A=false(m,n);
R=regionprops(BW,'PixelIdxList');
Nv=numel(R); %number of vertebrae
corners=cell(1,Nv);
imshow(BW)
for i=1:Nv
    
   B=A;
   B(R(i).PixelIdxList)=1;
   
   corners{i} = pgonCorners(B,4,40);
   
   hold on
   
   plot( corners{i}(:,2),corners{i}(:,1),'ro','MarkerFaceColor','m');
   
   hold off
end
heights=cellfun(@(c) mean(c(:,1)),corners);
[heights,idx]=sort(heights);
corners=corners(idx);

Best Answer

Related Solutions

MATLAB: How can i use region based segmentation to segment common carotid artery of ultrasound image and measure intima media thickness

MATLAB: How to find the coordinates of the top and bottom end plate and calculate the height

Related Question