Skip to content
Snippets Groups Projects
Commit b421cfd5 authored by fritsch's avatar fritsch
Browse files

Update inPhase_Segmentation_v11.m

parent 326c9d12
No related branches found
No related tags found
No related merge requests found
......@@ -301,18 +301,22 @@ for k = 1 : length(folderExp)
currImageCropBF2 = imgradient(currImageCropBF);
currImageCropBF2(~emulsionMask) = mean(currImageCropBF2(emulsionMask));
[centerBF, radiusBF, ~] = imfindcircles(imadjust(currImageCropBF2), [10, round(size(currImageCropBF2, 1)/5)], 'Sensitivity', 0.6, 'ObjectPolarity', 'dark');
% remove misdetection of bright spots
if ~isempty(centerBF)
currImageCropBF_Mask = createCirclesMask(size(currImageCropBF), centerBF, radiusBF);
currImageCropBF_Mask = regionprops(currImageCropBF_Mask, currImageCropBF, 'MeanIntensity', 'Centroid');
currImageCropBF_Int = vertcat(currImageCropBF_Mask.MeanIntensity);
currImageCropBF_Centroid = vertcat(currImageCropBF_Mask.Centroid);
currIndx = currImageCropBF_Int <= mean(currImageCropBF(emulsionMask));
if sum(currIndx) > 0
currIndx = sum(abs(centerBF - currImageCropBF_Centroid(currIndx, :)), 2) <= 2;
try
if ~isempty(centerBF)
currImageCropBF_Mask = createCirclesMask(size(currImageCropBF), centerBF, radiusBF);
currImageCropBF_Mask = regionprops(currImageCropBF_Mask, currImageCropBF, 'MeanIntensity', 'Centroid');
currImageCropBF_Int = vertcat(currImageCropBF_Mask.MeanIntensity);
currImageCropBF_Centroid = vertcat(currImageCropBF_Mask.Centroid);
currIndx = currImageCropBF_Int <= mean(currImageCropBF(emulsionMask));
if sum(currIndx) > 0
currIndx = sum(abs(centerBF - currImageCropBF_Centroid(currIndx, :)), 2) <= 2;
end
centerBF = centerBF(currIndx, :);
radiusBF = radiusBF(currIndx);
end
centerBF = centerBF(currIndx, :);
radiusBF = radiusBF(currIndx);
catch
end
if saveImage
hf1 = figure('visible', 'off');
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment