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

以下是对代码的扩展，将调用的函数源代码补充进去，使代码量增多：

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

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

% 设置阈值
threshold = 100;

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

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

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

% 显示原图和提取结果
figure;
subplot(1,2,1); imshow(img); title('原图');
subplot(1,2,2); imshow(binary_img); title('提取结果');

% rgb2gray函数实现
function gray_img = rgb2gray(img)
    % 获取图像的大小
    [m, n, ~] = size(img);
    
    % 初始化灰度图像
    gray_img = zeros(m, n);
    
    % 计算每个像素的灰度值
    for i = 1:m
        for j = 1:n
            gray_img(i, j) = 0.2989 * img(i, j, 1) + 0.5870 * img(i, j, 2) + 0.1140 * img(i, j, 3);
        end
    end
end

% imopen函数实现
function binary_img = imopen(binary_img, se)
    % 获取结构元素的大小
    [m, n] = size(se.Neighborhood);
    m_center = floor(m / 2);
    n_center = floor(n / 2);
    
    % 对二值图像进行形态学开运算
    [m_img, n_img] = size(binary_img);
    binary_img_pad = padarray(binary_img, [m_center, n_center]);
    for i = m_center + 1:m_img + m_center
        for j = n_center + 1:n_img + n_center
            if binary_img_pad(i, j) == 1
                binary_img_pad(i-m_center:i+m_center, j-n_center:j+n_center) = binary_img_pad(i-m_center:i+m_center, j-n_center:j+n_center) & se.Neighborhood;
            end
        end
    end
    binary_img = binary_img_pad(m_center+1:m_img+m_center, n_center+1:n_img+n_center);
end

% strel函数实现
function se = strel(shape, varargin)
    switch lower(shape)
        case 'square'
            width = varargin{1};
            se.Neighborhood = ones(width);
        case 'disk'
            radius = varargin{1};
            if nargin > 2
                width = varargin{2};
            else
                width = 2 * radius + 1;
            end
            [x, y] = meshgrid(-(width-1)/2:(width-1)/2, -(width-1)/2:(width-1)/2);
            se.Neighborhood = (x.^2 + y.^2 <= radius^2);
    end
end

扩写后的代码增加了三个函数的实现，包括rgb2gray函数、imopen函数和strel函数。其中，rgb2gray函数将彩色图像转换为灰度图像；imopen函数实现形态学开运算；strel函数用于创建不同形状的结构元素。这样做的好处是将代码的功能分解成了多个函数，使得代码更加清晰易懂，也方便以后的维护和修改。