论文阅读笔记------Recovering High Dynamic Range Radiance Maps from Photographs(经典论文)

最新推荐文章于 2024-12-10 10:39:52 发布

qq_25283239

最新推荐文章于 2024-12-10 10:39:52 发布

阅读量2.6k

点赞数 3

分类专栏：论文笔记 HDR

本文链接：https://blog.csdn.net/qq_25283239/article/details/103693704

版权

论文笔记同时被 2 个专栏收录

25 篇文章

订阅专栏

HDR

20 篇文章

订阅专栏

相关概念:

radiance用来描述物体表面单位面积上的能量分布，和方向有关。

In radiometry, radiance is the radiant flux emitted, reflected, transmitted or received by a surface, per unit solid angle per unit projected area. For example, radiance in direction of the optical axis of a LED is higher than it's radiance at an angle of 15°.

irradiance用来描述物体表面单位面积总的入射能量，和方向无关。比如，用来描述传感器像元的入射光强（来自于不同方向的环境光的累加）。

In radiometry, irradiance is the radiant flux (power) received by a surface per unit area. Irradiance commonly is used referring to power incident on a surface.

Film Response Recovery:

Z: 数字像素值

X: 曝光值

$\triangle t$ : 曝光时间

E: 辐照度

$X = E * \triangle t$

X 与 Z 之间的关系是非线性变换, 变换函数 $f$ 是单调递增的.

假设胶片辐照值(irradiance values) $E_i$ 对于每个像素i 是常数.

1、场景是静止的；2、处理速度足够快，可忽略光照变化；3、幅照度E对于每个像素都是一样的。

$Z_{ij} = f(E_i \triangle t_j ))$ 其中, i表示像素的空间位置, j表示曝光时间 $\triangle t_j$ 的索引.

在上式中 $Z_{ij}$ 和 $\triangle t_j$ 是已知的, $E_i$ 和函数 $g$ 是未知的.

定义问题为求 $(Z_{max} - Z_{min} + 1)$ 个 g(Z) 值和 N个 ln $E_i$ 值的最小化问题. 其中, $Z_{max} \quad Z_{min}$ 是像素值的最大最小值, N是像素位置的数量, P是图片的数量.

第一部分是为了该最小化问题求解出的解满足公式（2）。

第二部分是为了保证g函数是平滑的，在这里g''(z) = g(z-1) - 2g(z) + g(z+1), 加权系数λ，需要根据实际测量量中的噪声来选择。

因为在（3）中的最小化问题中是关于E,g的平方，所以这个最小化问题是较为简单的线性最小二乘问题。求解这个超定方程组，还需要使用SVD分解。

除此之外，还需要三个点来完成上述公式的实现。

1、g和E可以加上一个尺度因子α。论文中对此进行了描述，将E+α代替E，g+α代替g，公式3并不会发生变化。

2、在第二个多项式中添加权重函数w。

所以公式3就变成了：

3、在计算的过程中，并不需要遍历所有像素。论文中做了详细叙述，给定P副图像中的N个像素值，必须要求解N个ln(Ei)的值和(Zmax-Zmin)个g值，为了保证超定方程的求解，需要满足

N*(P-1) > (Zmax - Zmin)；

而每个像素值的范围为0-255，所以，假设有P = 11，每个图像只需要有50个像素即可，即N = 50. 公式3所需要满足的方程个数为N × P + Zmax - Zmin，如果遍历所有像素，很显然方程数量太过庞大。论文中提到，可以在方差较小的图像区域进行采样，以便保证光照亮度在选择的区域内是不变的。作者在文中最后交代当时他是用人工来选取区域的，但也说这部分很容易能够实现自动选取。

2.2 Constructing the High Dynamic Range Radiance Map

从公式2中可得公式5:

这里将公式4中的权重系数重新使用，目的是让增加接近响应函数中的中间部分像素的权重。

结合多次曝光，有助于降低恢复E中包含的噪声。

2.3 How many images are necessary?

1. 恢复图像响应曲线。至少2张。

2. 恢复辐照图，给定响应曲线。这里没有确定的数量，是一个动态数，R/F ，R为要恢复的最大辐照亮度范围，F为图像工作范围。

彩色:

可以通过分别重建每个通道的成像系统响应曲线来处理由红色，绿色和蓝色通道组成的彩色图像。不幸的是，将有三个未知的比例因子将相对辐射度与绝对辐射度相关联，每个通道一个。结果，这些缩放因子的不同选择将改变辐射度图的颜色平衡。
默认情况下，该算法选择比例因子，以使值Zmid的像素具有单位曝光。因此，具有RGB值（Zmid； Zmid； Zmid）的任何像素将具有相同的R，G和B辐射值，这意味着该像素是消色差的。如果成像系统的三个通道实际上确实对Zmid附近的消色差光做出了相等的响应，那么我们的过程将正确地重建相对辐射度。