HDR技术研究之Radiometric Self Calibration

一种基于论文Radiometric Self Calibration的相机响应函数计算方法

本文是Jesse Chen的原创文章。
在研究HDR技术时，会有一个比较成熟的HDR分支领域—相机响应函数的计算。作者在通读这个领域的经典论文时，注意到了经典论文Radiometric Self Calibration中给出的计算方法存在的计算不单调和不合理问题。针对这些问题，本博文的作者提出了新的计算方案，并编写程序，进行的仿真。

---

论文Radiometric Self Calibration 的回顾

我们都知道可通过一个响应函数$M=g(I)$描述成像系统最终产生的亮度观测值$M$和（经缩放的）环境光辐射$I$之间的关系。该响应函数一般都是单调递增的，因此我们可以用函数$g$的逆函数$f=g^{-1}$将所有的测量值M映射到环境光辐射值$I$，$I=f(M)$。

论文Radiometric Self Calibration描述了如何从一组同一场景的不同曝光图像中找到函数$f$。

值得注意的是，一个图像系统和另一个图像系统的响应曲线会有明显的差异，但是响应曲线不具有奇异的形式。首要原因是，图像传感器的输出是单调或半单调的。除非在图像系统中剩余模块中有不同寻常的映射，随着环境光辐射（曝光）的增加测量到的亮度要么也增加或者维持不变。因此，论文作者认为几乎所有的响应函数都可以用高阶多项式建模。

$$I=f(M)=\sum_{n=0}^{N}c_{n}M^{n}\tag{1}$$

因此只要确定了上述多项式的最高阶数$N$以及对应的多项式系数$c_n$就可以确定响应函数。

假定同一场景，分别用不同的曝光$e_q$和$e_{q+1}$得到了两张照片，设曝光比$R_{q,q+1}$为

$$R_{q,q+1}=\frac{e_q}{e_{q+1}}\tag{2}$$

因此图像系统的响应函数和曝光比之间满足：

$$\frac{f\left(M_{p,q}\right)}{f\left(M_{p,q+1}\right)}=R_{q,q+1}\tag{3}$$

对图像进行排序以满足$e_q<e_{q+1}$，因此$0<R_{q,q+1}<1$。将多项式模式带入上式，我们得到

∑Nn=0cnM