Single Image Haze Removal Using Dark Channel Prior

Kaiming HeJian SunX. Tang

article2011TPAMI1,770 citationsBest Paper Award (CVPR 2009)

Introduces the dark channel prior, a simple statistical discovery that enables accurate estimation of haze thickness to recover clear outdoor scenes and generate depth maps from a single degraded photograph.

Listen

Outdoor images frequently suffer from atmospheric haze, fog, and smoke, which scatter light, reduce contrast, and shift natural colors. This degradation impairs both human visual quality and downstream computer vision algorithms that assume clear scene radiance. Because haze density varies with unknown depth, restoring a clear image from a single picture has traditionally been an ill-posed problem requiring multiple exposures, polarization filters, or user-provided three-dimensional models.

The article demonstrates a novel, physically grounded method for haze removal from a single image using a newly identified statistical property called the dark channel prior. The core objective is to directly estimate haze thickness and restore clear scene radiance and depth information without requiring supplementary hardware or external 3D data.

The approach relies on an empirical observation across natural outdoor scenes: in local non-sky patches, at least one color channel typically exhibits very low intensity due to shadows, colorful surfaces, or dark objects. The authors validated this statistical prior on a dataset of 5,000 haze-free outdoor landscape and cityscape images. By integrating this prior into standard atmospheric scattering models, the method calculates haze transmission, refines edge boundaries using soft matting, and estimates atmospheric light from the most haze-opaque regions of the image.

The key findings confirm the effectiveness and robustness of this method. First, statistical analysis verified that approximately 75% of dark channel pixels in haze-free images have zero intensity, and 90% have intensities below 25 out of 255. Second, the algorithm effectively restores vivid colors and sharp structures even in heavily hazy scenes where previous methods fail due to faint color variation. Third, the process simultaneously extracts an accurate relative depth map as a direct by-product of dehazing. Fourth, the method prevents over-saturation and significantly suppresses halo artifacts around depth discontinuities, matching or exceeding the performance of techniques that rely on pre-existing 3D terrain models.

These findings provide immediate practical value for computer vision systems, aerial imaging, and photography by enhancing visibility and depth perception without expensive multi-sensor setups. However, the dark channel prior becomes invalid when scene objects are inherently bright, uniform, and match the atmospheric light without casting shadows (such as white marble buildings), which can lead to localized underestimations of transmission. Future work should integrate more complex atmospheric scattering models to handle non-uniform atmospheric effects, such as direct sunlight and horizon-specific color shifts.

He et al (2011).pdf
  • Paper: Single image dehazing, Raanan Fattal (2008). It formulates the core single-image dehazing problem using atmospheric scattering and surface shading statistics, establishing foundational concepts that the dark channel prior directly refines and improves upon.
  • Paper: A Closed-Form Solution to Natural Image Matting, Anat Levin et al. (2006). It introduces the matting Laplacian formulation used in the source paper for soft matting and refining coarse transmission maps.
  • Paper: Guided Image Filtering, Kaiming He et al. (2010). It introduces the guided filter, an edge-preserving smoothing filter designed as a fast alternative to the matting Laplacian for refining transmission maps in single image dehazing.
Cover for Single Image Haze Removal Using Dark Channel Prior

Abstract

In this paper, we propose a simple but effective image prior - dark channel prior to remove haze from a single input image. The dark channel prior is a kind of statistics of the haze-free outdoor images. It is based on a key observation - most local patches in haze-free outdoor images contain some pixels which have very low intensities in at least one color channel. Using this prior with the haze imaging model, we can directly estimate the thickness of the haze and recover a high quality haze-free image. Results on a variety of outdoor haze images demonstrate the power of the proposed prior. Moreover, a high quality depth map can also be obtained as a by-product of haze removal.

Table of Contents

  • 1. Introduction
  • 2. Background
  • 3. Dark Channel Prior
  • 4. Haze Removal Using Dark Channel Prior
  • 4.1. Estimating the Transmission
  • 4.2. Soft Matting
  • 4.3. Recovering the Scene Radiance
  • 4.4. Estimating the Atmospheric Light
  • 5. Experimental Results
  • 6. Discussions and Conclusions
  • References

Knowls

  1. Knowl 1 — Dark Channel and Dark Channel Prior Definition

    definition

    For an arbitrary color image JJ, the dark channel Jdark(x)J^{dark}(x) at pixel xx is defined as:

    Jdark(x)=min⁡c∈{r,g,b}(min⁡y∈Ω(x)(Jc(y)))J^{dark}(x) = \min_{c \in \{r, g, b\}} \left( \min_{y \in \Omega(x)} \left( J^c(y) \right) \right)

    where JcJ^c is the color channel c∈{r,g,b}c \in \{r, g, b\} of JJ, and Ω(x)\Omega(x) is a local square patch centered at pixel xx.

    The dark channel prior states that for natural, haze-free outdoor images, the intensity of Jdark(x)J^{dark}(x) is low and tends toward zero in non-sky regions:

    Jdark(x)≈0J^{dark}(x) \approx 0

    This low intensity stems from three primary image characteristics in outdoor scenes:

    1. Shadows (such as cast shadows from buildings, vehicles, trees, and rocks).
    2. Colorful objects or surfaces that have low reflectance in at least one color channel (such as green plants, red flowers, or blue water).
    3. Inherently dark objects or surfaces with low reflectance across all channels (such as dark tree trunks and soil).
  2. Knowl 2 — Empirical Statistics of Dark Channels in Outdoor Images

    empirical result

    The statistical validity of the dark channel prior was verified on a dataset of 5,000 haze-free outdoor cityscape and landscape images collected from image search engines using 150 popular user tags. The images were resized so that their maximum dimension was 500 pixels, sky regions were manually masked out, and dark channels were computed using a local patch size of 15×1515 \times 15 pixels.

    The statistical findings are:

    • Approximately 75% of the pixels in the dark channels have an intensity value equal to 0 (on an 8-bit scale [0,255][0, 255]).
    • Approximately 90% of the pixels have an intensity value below 25.
    • The histogram of the average dark channel intensity per image shows that the vast majority of haze-free outdoor images have near-zero average intensity, confirming that deviation from the prior is rare.
  3. Knowl 3 — Coarse Transmission Estimation via Dark Channel Prior

    equation

    In the standard haze image formation model I(x)=J(x)t(x)+A(1−t(x))I(x) = J(x)t(x) + A(1 - t(x)), II is the observed hazy image, JJ is the haze-free scene radiance, A=(Ar,Ag,Ab)TA = (A^r, A^g, A^b)^T is the global atmospheric light vector, and t(x)t(x) is the medium transmission. Under the assumption that transmission is locally constant t~(x)\tilde{t}(x) within patch Ω(x)\Omega(x), normalizing by AcA^c and taking the spatial and channel minimum yields the coarse transmission estimate:

    t~(x)=1−ωmin⁡c∈{r,g,b}(min⁡y∈Ω(x)(Ic(y)Ac))\tilde{t}(x) = 1 - \omega \min_{c \in \{r, g, b\}} \left( \min_{y \in \Omega(x)} \left( \frac{I^c(y)}{A^c} \right) \right)

    where min⁡c(min⁡y∈Ω(x)(Ic(y)/Ac))\min_c (\min_{y \in \Omega(x)} (I^c(y)/A^c)) is the dark channel of the normalized haze image.

    The parameter ω∈(0,1]\omega \in (0, 1] is an aerial perspective constant (fixed to ω=0.95\omega = 0.95) that adaptively retains a small fraction of haze for distant objects to preserve depth cues.

    In sky regions, Ic(x)≈AcI^c(x) \approx A^c, so the normalized dark channel approaches 1 and t~(x)→0\tilde{t}(x) \to 0, naturally representing infinite scene depth without requiring prior sky segmentation.

  4. Knowl 4 — Transmission Map Refinement Using Soft Matting Laplacian

    model/method

    Because the coarse transmission map t~(x)\tilde{t}(x) computed patch-wise contains block artifacts, it is refined to an edge-preserving continuous transmission map t(x)t(x) by minimizing the soft matting cost function:

    E(t)=tTLt+λ(t−t~)T(t−t~)E(t) = t^T L t + \lambda (t - \tilde{t})^T (t - \tilde{t})

    where tt and t~\tilde{t} are the vectorized forms of the refined and coarse transmission maps, λ=10−4\lambda = 10^{-4} is a regularization parameter, and LL is the Matting Laplacian matrix. The (i,j)(i, j)-th element of LL is defined as:

    Li,j=∑k∣(i,j)∈wk(δi,j−1∣wk∣(1+(Ii−μk)T(Σk+ε∣wk∣U3)−1(Ij−μk)))L_{i,j} = \sum_{k \mid (i,j) \in w_k} \left( \delta_{i,j} - \frac{1}{|w_k|} \left( 1 + (I_i - \mu_k)^T \left( \Sigma_k + \frac{\varepsilon}{|w_k|} U_3 \right)^{-1} (I_j - \mu_k) \right) \right)

    where Ii,IjI_i, I_j are the RGB colors at pixels ii and jj, δi,j\delta_{i,j} is the Kronecker delta, μk\mu_k is the 3×13 \times 1 mean color vector in window wkw_k, Σk\Sigma_k is the 3×33 \times 3 color covariance matrix in window wkw_k, U3U_3 is the 3×33 \times 3 identity matrix, ∣wk∣|w_k| is the number of pixels in window wkw_k, and ε\varepsilon is a regularizing parameter.

    The optimal transmission vector tt is computed by solving the linear system:

    (L+λU)t=λt~(L + \lambda U)t = \lambda \tilde{t}

    where UU is an identity matrix of the same size as LL.

  5. Knowl 5 — Atmospheric Light Estimation Using the Dark Channel

    algorithm

    Instead of selecting the globally brightest pixel in an image (which can be a white car, wall, or light source), the atmospheric light AA is estimated by finding the brightest pixel within the most haze-opaque regions identified by the dark channel:

    Input: Haze image II, dark channel IdarkI^{dark}
    Output: Atmospheric light vector A=(Ar,Ag,Ab)TA = (A^r, A^g, A^b)^T
    Identify the candidate set PP containing the top 0.1%0.1\% brightest pixels in IdarkI^{dark}
    Find pixel x∗=arg⁡max⁡x∈P(max⁡c∈{r,g,b}Ic(x))x^* = \arg\max_{x \in P} (\max_{c \in \{r,g,b\}} I^c(x))
    A←I(x∗)A \leftarrow I(x^*)
    return AA

    Restricting the candidate set to the top 0.1% pixels in IdarkI^{dark} ensures the search is localized to dense haze regions where airlight dominates scene radiance.

  6. Knowl 6 — Scene Radiance Recovery and Lower Bounding

    equation

    Given the input hazy image I(x)I(x), the atmospheric light AA, and the refined transmission t(x)t(x), the scene radiance J(x)J(x) is recovered by inverting the haze imaging equation:

    J(x)=I(x)−Amax⁡(t(x),t0)+AJ(x) = \frac{I(x) - A}{\max(t(x), t_0)} + A

    where t0t_0 is a fixed transmission lower bound set to t0=0.1t_0 = 0.1.

    The threshold t0t_0 prevents division by values near zero in dense haze regions, protecting the recovered radiance from noise amplification while preserving a small amount of haze in extremely dense areas. Because haze-free scene radiance is typically dimmer than atmospheric light, the exposure of J(x)J(x) is adjusted upward for display.

  7. Knowl 7 — Single Image Haze Removal Procedure

    algorithm

    The complete end-to-end dehazing algorithm recovers the haze-free radiance JJ and refined transmission map tt from a single RGB image II:

    Input: RGB haze image II, patch size 15×1515 \times 15, ω=0.95\omega = 0.95, λ=10−4\lambda = 10^{-4}, t0=0.1t_0 = 0.1
    Output: Dehazed image JJ, transmission map tt
    for each pixel xx in II do
        Idark(x)←min⁡c∈{r,g,b}min⁡y∈Ω(x)Ic(y)I^{dark}(x) \leftarrow \min_{c \in \{r,g,b\}} \min_{y \in \Omega(x)} I^c(y)
    P←top 0.1% brightest pixels in IdarkP \leftarrow \text{top } 0.1\% \text{ brightest pixels in } I^{dark}
    x∗←arg⁡max⁡x∈P(max⁡c∈{r,g,b}Ic(x))x^* \leftarrow \arg\max_{x \in P} (\max_{c \in \{r,g,b\}} I^c(x))
    A←I(x∗)A \leftarrow I(x^*)
    for each pixel xx in II do
        t~(x)←1−ωmin⁡c∈{r,g,b}min⁡y∈Ω(x)(Ic(y)/Ac)\tilde{t}(x) \leftarrow 1 - \omega \min_{c \in \{r,g,b\}} \min_{y \in \Omega(x)} (I^c(y) / A^c)
    Construct Matting Laplacian matrix LL from II
    Solve (L+λU)t=λt~(L + \lambda U)t = \lambda \tilde{t} using Preconditioned Conjugate Gradient (PCG)
    for each pixel xx and color channel cc do
        Jc(x)←(Ic(x)−Ac)/max⁡(t(x),t0)+AcJ^c(x) \leftarrow (I^c(x) - A^c) / \max(t(x), t_0) + A^c
    return J,tJ, t

    The local minimum operation is executed in O(N)O(N) linear time with respect to image size NN using van Herk's fast algorithm. Solving the linear system with PCG takes 10–20 seconds for a 600×400600 \times 400 image on a 3.0 GHz processor.

  8. Knowl 8 — Depth Map Reconstruction from Medium Transmission

    model/method

    In a homogeneous atmosphere with scattering coefficient β\beta, the medium transmission t(x)t(x) decays exponentially with scene depth d(x)d(x):

    t(x)=e−βd(x)t(x) = e^{-\beta d(x)}

    Using the refined transmission map t(x)t(x), the relative scene depth d(x)d(x) is reconstructed directly up to an unknown global scaling parameter β\beta:

    d(x)=−1βln⁡t(x)d(x) = -\frac{1}{\beta} \ln t(x)

    This produces a dense, depth-edge-preserving depth map of the scene as a direct byproduct of single-image haze removal.

  9. Knowl 9 — Limitations and Failure Cases of the Dark Channel Prior

    limitation

    The dark channel prior fails under two principal scenarios:

    1. Large, inherently bright objects without shadows: When a scene region contains objects whose intrinsic radiance is bright and uniform across all color channels (such as a large white marble surface or white car) with no cast shadows, Jdark(x)J^{dark}(x) is non-zero. The dark channel method interprets this high intensity as haze, underestimating transmission t~(x)\tilde{t}(x) and producing unnaturally dark, over-saturated output in those regions.
    2. Invalidation of the atmospheric scattering model: The standard model I(x)=J(x)t(x)+A(1−t(x))I(x) = J(x)t(x) + A(1 - t(x)) assumes uniform atmospheric light and fails to capture non-homogeneous illumination effects, such as direct solar glare in the sky or wavelength-dependent Rayleigh scattering that introduces a bluish hue near the horizon.
  10. Knowl 10 — Grayscale Image Dehazing Using Dark Channel Prior

    model/method

    The dark channel dehazing method applies to single-channel grayscale images I(x)I(x) when sufficient local shadows or dark objects exist. Because color channels are absent, the cross-channel minimum is dropped, and the coarse transmission is given by:

    t~(x)=1−ωmin⁡y∈Ω(x)(I(y)A)\tilde{t}(x) = 1 - \omega \min_{y \in \Omega(x)} \left( \frac{I(y)}{A} \right)

    where AA is the scalar atmospheric light estimated from the top 0.1% brightest pixels of the local minimum filtered image. Refinement uses the grayscale Matting Laplacian, and the grayscale scene radiance J(x)J(x) is recovered via J(x)=(I(x)−A)/max⁡(t(x),t0)+AJ(x) = (I(x) - A)/\max(t(x), t_0) + A.

Coverage note — No substantial contributed material was omitted; all core components including prior definition, statistical verification, transmission estimation, matting refinement, atmospheric light estimation, radiance recovery, depth extraction, grayscale extension, and failure cases are included.

References

  1. 1.P. Chavez. An improved dark-object substraction technique for atmospheric scattering correction of multispectral data. Remote Sensing of Environment, 24:450–479, 1988.
  2. 2.R. Fattal. Single image dehazing. In SIGGRAPH, pages 1–9, 2008.
  3. 3.E. B. Goldstein. Sensation and perception. 1980.
  4. 4.E. Hsu, T. Mertens, S. Paris, S. Avidan, and F. Durand. Light mixture estimation for spatially varying white balance. In SIGGRAPH, pages 1–7, 2008.
  5. 5.J. Kopf, B. Neubert, B. Chen, M. Cohen, D. Cohen-Or, O. Deussen, M. Uyttendaele, and D. Lischinski. Deep photo: Model-based photograph enhancement and viewing. SIGGRAPH Asia, 2008.
  6. 6.H. Koschmieder. Theorie der horizontalen sichtweite. Beitr. Phys. Freien Atm., 12:171–181, 1924.
  7. 7.A. Levin, D. Lischinski, and Y. Weiss. A closed form solution to natural image matting. CVPR, 1:61–68, 2006.
  8. 8.S. G. Narasimhan and S. K. Nayar. Chromatic framework for vision in bad weather. CVPR, pages 598–605, 2000.
  9. 9.S. G. Narasimhan and S. K. Nayar. Vision and the atmosphere. IJCV, 48:233–254, 2002.
  10. 10.S. G. Narasimhan and S. K. Nayar. Contrast restoration of weather degraded images. PAMI, 25:713–724, 2003.
  11. 11.S. G. Narasimhan and S. K. Nayar. Interactive deweathering of an image using physical models. In Workshop on Color and Photometric Methods in Computer Vision, 2003.
  12. 12.S. K. Nayar and S. G. Narasimhan. Vision in bad weather. ICCV, page 820, 1999.
  13. 13.A. J. Preetham, P. Shirley, and B. Smits. A practical analytic model for daylight. In SIGGRAPH, pages 91–100, 1999.
  14. 14.Y. Y. Schechner, S. G. Narasimhan, and S. K. Nayar. Instant dehazing of images using polarization. CVPR, 1:325, 2001.
  15. 15.S. Shwartz, E. Namer, and Y. Y. Schechner. Blind haze separation. CVPR, 2:1984–1991, 2006.
  16. 16.R. Tan. Visibility in bad weather from a single image. CVPR, 2008.
  17. 17.M. van Herk. A fast algorithm for local minimum and maximum filters on rectangular and octagonal kernels. Pattern Recogn. Lett., 13:517–521, 1992.

Citation

MLA
Kaiming He, et al. “Single Image Haze Removal Using Dark Channel Prior”. IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 33, no. 12, 2011, pp. 2341–53, https://doi.org/10.1109/TPAMI.2010.168.
APA
Kaiming He, Jian Sun, & Xiaoou Tang. (2011). Single Image Haze Removal Using Dark Channel Prior. IEEE Transactions on Pattern Analysis and Machine Intelligence, 33(12), 2341–2353. https://doi.org/10.1109/TPAMI.2010.168
Chicago
Kaiming He, Jian Sun, and Xiaoou Tang. 2011. “Single Image Haze Removal Using Dark Channel Prior”. IEEE Transactions on Pattern Analysis and Machine Intelligence 33 (12): 2341–53. https://doi.org/10.1109/TPAMI.2010.168.
Harvard
Kaiming He, Jian Sun and Xiaoou Tang (2011) “Single Image Haze Removal Using Dark Channel Prior”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 33(12), pp. 2341–2353. Available at: https://doi.org/10.1109/TPAMI.2010.168.
Vancouver
1. Kaiming He, Jian Sun, Xiaoou Tang (2011) Single Image Haze Removal Using Dark Channel Prior. IEEE Transactions on Pattern Analysis and Machine Intelligence 33:2341–2353

BibTeX

@article{Kaiming_He_2011, title={Single Image Haze Removal Using Dark Channel Prior}, volume={33}, ISSN={2160-9292}, url={http://dx.doi.org/10.1109/TPAMI.2010.168}, DOI={10.1109/tpami.2010.168}, number={12}, journal={IEEE Transactions on Pattern Analysis and Machine Intelligence}, publisher={Institute of Electrical and Electronics Engineers (IEEE)}, author={Kaiming He and Jian Sun and Xiaoou Tang}, year={2011}, month=Dec, pages={2341–2353} }
Metadata:Crossref

Access the Paper

This paper is available from its original source. Click below to access the PDF.

Open PDF