将上面的代码完整的写出来

下面是一个完整的MATLAB程序，用于将gprMax输出文件中的B-scan图像转换为去除直达波的图像：

% 读取.gprMax输出文件
data = load('output.mat');
scan = data.data;

% 定义均值过滤函数
function scan_new = remove_direct_wave(scan, kernel_size)
    % 进行水平方向上的平均滤波得到平均图像
    scan_hmean = imgaussfilt(scan, [kernel_size, 1]);

    % 进行垂直方向上的平均滤波得到全局平均值
    scan_vmean = imgaussfilt(scan_hmean, [1, kernel_size]);
    global_mean = mean(scan_vmean(:));

    % 对原始图像进行剪切
    scan_new = max(scan - global_mean, 0);
end

% 对B-scan图像进行去除直达波处理
kernel_size = 5;
scan_new = remove_direct_wave(scan, kernel_size);

% 绘制去除直达波后的B-scan图像
figure;
imagesc(scan_new);
colormap('gray');
axis image;
colorbar;
title('B-scan after removing direct wave');

在上面的代码中，`output.mat`表示gprMax的输出文件名，`kernel_size`为均值滤波器的大小，根据具体情况进行调整。函数`remove_direct_wave`将原始B-scan图像进行去除直达波处理，并返回去除直达波后的图像`scan_new`。程序最后会绘制出去除直达波后的B-scan图像，并添加颜色条和标题。