将下列代码调用的函数源代码具体实现补充进去，使代码量增多，并将扩写后的代码展示出来：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('提取结果');

% 定义imread函数
function img = imread(filename)
    % 打开文件
    fileID = fopen(filename,'r');

    % 读取文件头
    A = fread(fileID, 2, 'uint8');
    width = fread(fileID, 1, 'uint16');
    height = fread(fileID, 1, 'uint16');

    % 读取像素数据
    data = fread(fileID, width*height, 'uint8');

    % 关闭文件
    fclose(fileID);

    % 将像素数据转换为图像矩阵
    img = reshape(data, width, height)';
end

% 定义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

% 定义strel函数
function se = strel(shape, size)
    % 根据形状和大小生成结构元素矩阵
    if strcmp(shape, 'disk')
        radius = size;
        [x, y] = meshgrid(-radius:radius, -radius:radius);
        se = (x.^2 + y.^2 <= radius.^2);
    else
        se = ones(size);
    end
end

% 定义imopen函数
function binary_img = imopen(binary_img, se)
    % 获取结构元素的大小
    [m, n] = size(se);

    % 对二值图像进行膨胀操作
    dilated_img = zeros(size(binary_img));
    for i = 1:m
        for j = 1:n
            if se(i, j)
                dilated_img(i:end-m+i, j:end-n+j) = max(dilated_img(i:end-m+i, j:end-n+j), binary_img);
            end
        end
    end

    % 对膨胀后的图像进行腐蚀操作
    binary_img = zeros(size(binary_img));
    for i = 1:m
        for j = 1:n
            if se(i, j)
                binary_img(i:end-m+i, j:end-n+j) = min(binary_img(i:end-m+i, j:end-n+j), dilated_img);
            end
        end
    end
end

这里我们分别实现了`imread`、`rgb2gray`、`strel`和`imopen`的具体功能。`imread`函数实现了从文件中读取图像数据的功能，`rgb2gray`函数实现了将彩色图像转换为灰度图像的功能，`strel`函数实现了生成结构元素矩阵的功能，`imopen`函数实现了对二值图像进行开运算的功能。这样做的好处是将代码分成了多个函数，使得代码更加模块化，易于维护和修改。同时，这样做也增加了代码量，使得我们更容易理解和掌握这段代码的功能。