how to segment trus prostate image?

Hi. I am trying to segment the prostate from TRUS image. This is what I have so far (see attachment). Does anyone know how to :

1. remove the unrelated object (ie. the white curve on top of prostate in last figure)
2. outline/connect the white pixel to form the prostate boundary?

This is the code:

clc;
clear;
close all;
i=imread('trus1.jpg');
i=im2uint8(rgb2gray(i));
[m,n]=size(i);
I=imresize(i,[2*m,2*n]);

figure, imshow(I), title('original image');

[~, threshold] = edge(I, 'sobel');
fudgeFactor = .5;
BWs = edge(I,'sobel', threshold * fudgeFactor);
figure, imshow(BWs), title('binary gradient mask');

se90 = strel('line', 3, 90);
se0 = strel('line', 3, 0);

BWsdil = imdilate(BWs, [se90 se0]);
figure, imshow(BWsdil), title('dilated gradient mask');
BWdfill = imfill(BWsdil,'holes');
figure, imshow(BWdfill);
title('binary image with filled holes');

BWnobord = imclearborder(BWdfill,6);
figure, imshow(BWnobord), title('cleared border image');

Answer 1

[Image Analyst]

THere is nothing built in to do that so you'll have to search the literature to see if anyone has done that before. Try PubMed, or VisBib.

http://www.visionbib.com/bibliography/contentsmedical.html#Medical%20Applications,%20CAT,%20MRI,%20Ultrasound,%20Heart%20Models,%20Brain%20Models