Hi, i am trying to perform image segmentation in scilab. I am trying to segment the breast breast region in a digital mammogram. I have tried using the follwing code to implement the watershed transform but this does not work.

Image=imread('C:\users\Ciara\Pictures\mammograms\mammogram3.jpg'); //% load image

Image=RGB2Gray(Image);
global EDGE_SOBEL;
GradientImage=EdgeFilter(Image, EDGE_SOBEL);
Threshold=CalculateOtsuThreshold(Image);
EdgeImage=~SegmentByThreshold(GradientImage, Threshold);
DistanceImage=DistanceTransform(EdgeImage);
ThresholdImage=SegmentByThreshold(DistanceImage, 8);
MarkerImage=SearchBlobs(ThresholdImage);
SegmentedImage=uint8(Watershed(GradientImage, MarkerImage));

figure();ShowImage(SegmentedImage, 'SegementedImage');
endfunction


Could anyone please help me discover where i am going wrong or suggest a better way of doing this?

Thanks