function vessel = extractVessel(data,datMean,thr)
    vesselRegion = datMean<thr;
    vesselRegion = bwareaopen(vesselRegion,1000,4);
    vesselRegion = myGrow2D(1-datMean,vesselRegion,2,150);
    BW = boundarymask(vesselRegion);
    figure;imshow(imoverlay(datMean,BW,'blue'),'InitialMagnification',67);
    vessel = find(vesselRegion>0);

end