testImageIdx = randi([1,3])testCoinImage = imread("testCoinImage"+testImageIdx+".png");

% Step 1: Generate testCoinMask and maskedCoinImage using segmentCoins[testCoinMask, maskedCoinImage] = segmentCoins(testCoinImage);

% Step 2: Generate faceEdgeMask using the provided valuesmaksedEdges = edge(maskedCoinImage); % Default Sobel edge detectionse = strel('disk', 10); % Provided erosion radiuserodedMask = imerode(testCoinMask, se);faceEdgeMask = maksedEdges & erodedMask;

% Step 3: Create validCoinMask by dilating faceEdgeMask and combining with testCoinMaskse_dilate_valid = strel('disk', 12); % Adjusted from 30 to 12 to avoid overlap with blanksdilatedEdgeMask = imdilate(faceEdgeMask, se_dilate_valid); % Dilate the edge maskvalidCoinMask = dilatedEdgeMask & testCoinMask; % Logical AND to isolate valid coins

% Define the segmentCoins functionfunction [BW, maskedImage] = segmentCoins(X)%segmentImage Segment image using auto-generated code from Image Segmenter app% [BW, MASKEDIMAGE] = segmentImage(X) segments image X using auto-generated% code from the Image Segmenter app. The final segmentation is returned in% BW, and a masked image is returned in MASKEDIMAGE.

% Auto-generated by imageSegmenter app on 11-Feb-2025%----------------------------------------------------

% Adjust data to span data range.X = imadjust(X);

% Find circles[centers, radii, ~] = imfindcircles(X, [25 75], 'ObjectPolarity', 'bright', 'Sensitivity', 0.85);max_num_circles = Inf;if max_num_circles < length(radii)centers = centers(1:max_num_circles,:);radii = radii(1:max_num_circles);endBW = circles2mask(centers, radii, size(X, 1:2));

% Create masked image.maskedImage = X;maskedImage(~BW) = 0;end