将下列代码在不改变功能的情况下扩写，将调用的函数源代码补充进去，使代码量增多，并将扩写后的代码展示出来：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() 和 imopen() 函数的源代码：

% 读入图片
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