使用提供的何凯明暗通道去雾MATLAB代码实现图像去雾

何凯明的暗通道先验（Dark Channel Prior，DCP）是一种用于单图像去雾的经典算法。在MATLAB中，你可以通过以下步骤来实现这个过程：

1. **导入所需库**:
首先，你需要安装`imregionalmax`函数，这在早期版本的MATLAB中可能已经包含，如果没有，可以使用`iptregionalmax`代替。如果需要下载，可以访问`https://www.mathworks.com/help/images/ref/imregionalmax.html`。

2. **读取图像并计算暗通道**:

   % 读取输入图像
   img = imread('your_image.jpg');
   
   % 转换为灰度图，如果原图是彩色的
   gray_img = rgb2gray(img);
   
   % 计算每个像素的暗通道值 (通常取R通道)
   dark_channel = im(:,:,1); % 或者其他通道，如dark_channel = max(gray_img,[],3);

3. **检测透明区域**:

   % 找出局部最小值作为潜在的清晰区域候选
   [D, idx] = imregionalmin(dark_channel, ones([5 5])); % 选择5x5大小的窗口
   
   % 提取清晰区域
   clear_mask = false(size(gray_img));
   clear_mask(idx) = true;

4. **恢复去雾后的图像**:

   % 使用暗通道信息估计大气光照
   atmospheric_light = mean(mean(img .* ~clear_mask), 3);
   
   % 去除大气影响
   dehazed_img = atmospheric_light ./ dark_channel;
   
   % 限制结果在0到1之间
   dehazed_img = min(max(dehazed_img, 0), 1);

5. **显示结果**:

   imshowpair(img, dehazed_img, 'montage', []);
   title('Original Image vs Dehazed Image');