将下列代码在不改变功能的情况下扩写，将调用的函数源代码补充进去，使代码量增多，并将扩写后的代码展示出来：% 读入图片 img = imread('瑕疵图像.png'); % 将彩色图像转换为灰度图像 gray_img = rgb2gray(img); % 设定阈值 threshold = 100; % 将灰度图像转换为二值图像 binary_img = gray_img > threshold; % 定义结构元素 se = strel('disk', 2); % 对二值图像进行开运算 binary_img = imopen(binary_img, se); % 显示原图和提取结果 subplot(1,2,1); imshow(img); title('原图'); subplot(1,2,2); imshow(binary_img); title('提取结果'); % rgb2gray() 函数的源代码 function gray_img = rgb2gray(img) % 将彩色图像转换为灰度图像 % 灰度值计算公式：gray_value = 0.299 * R + 0.587 * G + 0.114 * B gray_img = img(:,:,1) * 0.299 + img(:,:,2) * 0.587 + img(:,:,3) * 0.114; end % imopen() 函数的源代码 function img = imopen(binary_img, se) % 对二值图像进行开运算 % 开运算：先腐蚀，再膨胀 eroded_img = imerode(binary_img, se); img = imdilate(eroded_img, se); end

将下列代码在不改变功能的情况下扩写，将调用的函数源代码补充进去，使代码量增多，并将扩写后的代码展示出来：

% 读入图片
img = imread('瑕疵图像.png');

% 将彩色图像转换为灰度图像
gray_img = rgb2gray(img);

% 设定阈值
threshold = 100;

% 将灰度图像转换为二值图像
binary_img = gray_img > threshold;

% 定义结构元素
se = strel('disk', 2);

% 对二值图像进行开运算
binary_img = imopen(binary_img, se);

% 对二值图像进行闭运算
binary_img = imclose(binary_img, se);

% 对二值图像进行腐蚀
eroded_img = imerode(binary_img, se);

% 对二值图像进行膨胀
dilated_img = imdilate(binary_img, se);

% 显示原图和提取结果
subplot(2,3,1); imshow(img); title('原图');
subplot(2,3,2); imshow(binary_img); title('开运算结果');
subplot(2,3,3); imshow(eroded_img); title('腐蚀结果');
subplot(2,3,4); imshow(binary_img); title('闭运算结果');
subplot(2,3,5); imshow(dilated_img); title('膨胀结果');

% rgb2gray() 函数的源代码
function gray_img = rgb2gray(img)
% 将彩色图像转换为灰度图像
% 灰度值计算公式：gray_value = 0.299 * R + 0.587 * G + 0.114 * B
gray_img = img(:,:,1) * 0.299 + img(:,:,2) * 0.587 + img(:,:,3) * 0.114;
end

% imopen() 函数的源代码
function img = imopen(binary_img, se)
% 对二值图像进行开运算
% 开运算：先腐蚀，再膨胀
eroded_img = imerode(binary_img, se);
img = imdilate(eroded_img, se);
end

% imclose() 函数的源代码
function img = imclose(binary_img, se)
% 对二值图像进行闭运算
% 闭运算：先膨胀，再腐蚀
dilated_img = imdilate(binary_img, se);
img = imerode(dilated_img, se);
end

% imerode() 函数的源代码
function img = imerode(binary_img, se)
% 对二值图像进行腐蚀
[m, n] = size(binary_img);
[x, y] = size(se.Neighborhood);
img = zeros(m, n);
for i = 1:m
for j = 1:n
if binary_img(i, j) == 1
flag = 1;
for k = 1:x
for l = 1:y
if se.Neighborhood(k, l) == 1
if i+k-x <= 0 || j+l-y <= 0 || i+k-x > m || j+l-y > n || binary_img(i+k-x, j+l-y) == 0
flag = 0;
break
end
end
end
if flag == 0
break
end
end
if flag == 1
img(i, j) = 1;
end
end
end
end
end

% imdilate() 函数的源代码
function img = imdilate(binary_img, se)
% 对二值图像进行膨胀
[m, n] = size(binary_img);
[x, y] = size(se.Neighborhood);
img = zeros(m, n);
for i = 1:m
for j = 1:n
if binary_img(i, j) == 1
for k = 1:x
for l = 1:y
if se.Neighborhood(k, l) == 1
if i+k-x <= 0 || j+l-y <= 0 || i+k-x > m || j+l-y > n
continue
end
img(i+k-x, j+l-y) = 1;
end
end
end
end
end
end
end