该方法由Debevec在1997年的论文《Recovering High Dynamic Range Radiance Maps from Photographs》中提出。原文及翻译可以参考博客《Recovering High Dynamic Range Radiance Maps from Photographs》。
变量约定
X X X:曝光量(the exposure),单位: J m − 2 Jm^{-2} Jm−2;
E E E:传感器或胶片上的辐照度(the irradiance);
Δ t \Delta t Δt:曝光时间(the exposure time);
X X X, E E E, Δ t \Delta t Δt的关系: X = E Δ t X=E\Delta t X=EΔt。
Z Z Z:数字图像的像素值;
f f f:相机的响应函数(Camera response function);
Z Z Z、 X X X、 f f f的关系: Z = f ( X ) Z=f(X) Z=f(X)。
算法
对同一场景使用不同的曝光时间进行拍摄,得到 P P P张图像,每一张的曝光时间记为 Δ t j \Delta t_{j} Δtj, 1 ≤ j ≤ P 1\leq j\leq P 1≤j≤P。每张图片上都取相同位置的 N N N个点, i i i为像素点的索引, 1 ≤ i ≤ N 1\leq i\leq N 1≤i≤N。
假设场景是静态的,并且曝光过程完成得足够快,可以安全地忽略灯光的变化。然后可以假设像素 i i i的辐照度值 E i E_{i} Ei是常数。
令 Z i j Z_{ij} Zij表示第 j j j张图片的第 i i i个像素的值,则有
Z i j = f ( E i Δ t j ) Z_{ij}=f(E_{i}\Delta t_{j}) Zij=f(EiΔtj)
假设函数 f f f单调可逆,则有
f − 1 ( Z i j ) = E i Δ t j f^{-1}(Z_{ij})=E_{i}\Delta t_{j} f−1(Zij)=EiΔtj
对上式两边取对数,则有
ln f − 1 ( Z i j ) = ln E i + ln Δ t j \ln f^{-1}(Z_{ij})=\ln E_{i}+\ln \Delta t_{j} lnf−1(Zij)=ln