用matlab实现.编写一个局域子窗滤波函数(边缘位置数值不变，其它每个位置，计算周围9个点高程的中位数，作为该位置的输出值)

以下是使用Matlab实现的局域子窗滤波函数，使用中值滤波器，对输入矩阵进行滤波，输出滤波后的矩阵：

function output = local_median_filter(input, window_size)
    % 计算窗口的半径
    r = floor(window_size / 2);
    % 获取输入矩阵的尺寸
    [height, width] = size(input);
    % 创建输出矩阵
    output = zeros(height, width);
    % 对每个像素进行滤波
    for i = 1:height
        for j = 1:width
            % 计算窗口的位置
            x_min = max(j - r, 1);
            x_max = min(j + r, width);
            y_min = max(i - r, 1);
            y_max = min(i + r, height);
            % 将窗口内的像素值复制到窗口数组中
            window = input(y_min:y_max, x_min:x_max);
            % 对窗口数组进行排序
            window = sort(window(:));
            % 计算中位数
            mid = floor(numel(window) / 2) + 1;
            if mod(numel(window), 2) == 0
                output(i, j) = (window(mid - 1) + window(mid)) / 2;
            else
                output(i, j) = window(mid);
            end
        end
    end
end

该函数使用了两重循环对每个像素进行滤波，对于每个像素，计算出窗口的位置，并将窗口内的像素值复制到窗口矩阵中，然后对窗口矩阵进行排序，计算中位数并将其作为该像素的输出值。注意，在边缘位置，由于窗口会超出矩阵边界，因此需要对窗口进行裁剪，只考虑矩阵内部的像素。