function watershedSeg(org, sob55);

subplot(3,3,1), imshow(org,[]), title('original');

%calculate gradient magnitude and use it as segmentaion function
sob55 = sobel55(org);
subplot(3,3,2), imshow(sob55,[]), title('sobel');

%After applying watershed transform on gradient magnitude -- oversegmented
W = watershed(sob55);
Lrgb = label2rgb(W);
subplot(3,3,3), imshow(Lrgb,[]), title('WT on gradient');

%marking the foregroung objects
%morphological operators transform the original image into another image through the iteration 
%with other image of a certain shape and size which is known as structuring element
structureElem = strel('disk', 7);

%opening by reconstruction 
imgErode = imerode(org, structureElem);
imgRecons = imreconstruct(imgErode, org);
%subplot(3,3,2), imshow(imgRecons), title('Open Reconstruction')

%closing by reconstr
imgDilate = imdilate(imgRecons, structureElem);
imgRecons2 = imreconstruct(imcomplement(imgDilate), imcomplement(imgRecons));
imgRecons2 = imcomplement(imgRecons2);
%subplot(3,3,3), imshow(imgRecons2), title('close reconstruction')

%calculate regional maxima
regMaxima = imregionalmax(imgRecons2);
subplot(3,3,4), imshow(regMaxima), title('Region maxima');
org2 = org;
org2(regMaxima) = 255;
subplot(3,3,5),imshow(org2,[]), title('maxima superimposed on original')

%Clean the edges of the markers using edge reconstruction
structureElem2 = strel(ones(5,5));
regMaxima2 = imclose(regMaxima, structureElem2);
regMaxima3 = imerode(regMaxima2, structureElem2);
regMaxima4 = bwareaopen(regMaxima3, 20);
org3 = org;
org3(regMaxima4) = 255;

%Compute background markers
bw = imbinarize(imgRecons2);
subplot(3,3,6), imshow(bw,[]), title('binarized')

%distance transform of binarized
distTrans = bwdist(bw);
distTransW = watershed(distTrans);
bgm = distTransW == 0;
subplot(3,3,7), imshow(bgm,[]), title('background maxima')

%Compute the Watershed Transform of the Segmentation Function
sob2 = imimposemin(sob55, bgm | regMaxima4);
W2 = watershed(sob2);

org4 = org;
org4(imdilate(W2 == 0, ones(3, 3)) | bgm | regMaxima4) = 255;
subplot(3,3,8),imshow(org4,[]), title('Markers/boundaries superimposed on original')
Lrgb2 = label2rgb(W2, 'jet', 'w', 'shuffle');
subplot(3,3,9), imshow(Lrgb2,[]),title('watershed label')