计算灰度图像的暗通道函数matlab程序

以下是计算灰度图像的暗通道函数的matlab程序：

function darkChannel = computeDarkChannel(img, patchSize)

% Convert the image to grayscale
if size(img, 3) == 3
    img = rgb2gray(img);
end

% Pad the input image to avoid boundary issues
padSize = floor(patchSize/2);
img = padarray(img, [padSize padSize], Inf);

% Compute the dark channel for each pixel
darkChannel = zeros(size(img));
for i = 1:size(img, 1)
    for j = 1:size(img, 2)
        patch = img(i-padSize:i+padSize, j-padSize:j+padSize);
        darkChannel(i,j) = min(patch(:));
    end
end

% Remove the padding from the dark channel
darkChannel = darkChannel(padSize+1:end-padSize, padSize+1:end-padSize);

end

该函数接受两个输入参数：原始图像和用于计算暗通道函数的局部块大小。该程序首先将彩色图像转换为灰度图像，然后对输入图像进行填充以避免边界问题。接下来，程序遍历图像中的每个像素，并在该像素周围的局部块中计算最小值。最后，程序从计算出的暗通道函数中删除填充像素，并将其作为输出返回。