适用于 v1.1.0+ 14 分钟阅读
Color Science 01 / 数学域解耦

负片反相的三大数学域

透射率域、光学密度域与正片输出域的严谨物理推导与非线性不可交换性定理。

# 透射率域:对数变换前的光电响应与设备基底

在胶片数字化工程中,相机传感器或扫描仪线阵电荷耦合器件接收到的物理量并非密度,而是通过胶片乳剂层后到达感光单元的光子通量。透射率域特指对数变换之前、与物理光强呈正比例关系的线性光域。此阶段的目标是获取具有明确物理意义、严格位于正值区间的相对透射率向量。

假定入射胶片平面的开场照明强度为 Iopen(λ)I_{\text{open}}(\lambda),穿透显影后底片的光强为 Itransmitted(λ)I_{\text{transmitted}}(\lambda),传感器内部暗电流与环境杂散光通量为 Idark(λ)I_{\text{dark}}(\lambda),则波长 λ\lambda 处的相对透射率定义为:

T(λ)=Itransmitted(λ)−Idark(λ)Iopen(λ)−Idark(λ)T(\lambda) = \frac{I_{\text{transmitted}}(\lambda) - I_{\text{dark}}(\lambda)}{I_{\text{open}}(\lambda) - I_{\text{dark}}(\lambda)}

在离散多通道采集设备中,传感器通道响应是透射光谱与滤光片光谱透过率及传感器量子效率的积分。令设备三通道原生输出为 Idevice=[R,G,B]TI_{\text{device}} = [R, G, B]^T。在透射率域内,所有物理运算必须严格保持光子通量的线性叠加性:

  1. 暗电平与暗场扣除:减去传感器读出热噪声与光机腔体内部杂散反射,确保完全无光透射时的物理读数为零。
  2. 平场增益归一化:通过无底片透光的背光场测量消除透镜边缘余弦四次方照度衰减及灯板各区域色温漂移。
  3. 设备通道分离:消除相机拜尔滤色片或扫描仪分光棱镜在三原色之间的光谱重叠,使得分离后的通道能够准确衡量底片青色、品红与黄色染料层的光吸收。

如果采集端未施加硬件级的通道分离校准,则直接进入后续计算会导致各染料层的吸收信息严重混叠。透射率数据必须钳位在微小正实数门限 ϵ\epsilon 以上,以防止后续对数运算出现奇异值:

Tvalid=max⁡(T,ϵ),ϵ=1.0×10−6T_{\text{valid}} = \max(T, \epsilon), \quad \epsilon = 1.0 \times 10^{-6}

重要

透射率域的数据本质是辐射度量学物理量。任何属于显示器特性的伽马编码、非线性对比度曲线或带有感知权重的色调映射,均严禁在此阶段介入。

# 光学密度域:对数变换后的吸收叠加与片基扣除

透射率信号经过常用对数变换后,进入光学密度域。根据比尔-朗伯吸收定律,平行单色光穿过均匀吸收介质时,光强呈指数级衰减。取对数运算将乘性衰减转换为加性吸收,从而为胶片三层染料吸收分析提供了线性代数基础。

未校正的原始采集密度向量 DrawD_{\text{raw}} 定义为相对透射率倒数的常用对数:

Draw=−log⁡10(Tvalid)=[−log⁡10(TR)−log⁡10(TG)−log⁡10(TB)]D_{\text{raw}} = -\log_{10}(T_{\text{valid}}) = \begin{bmatrix} -\log_{10}(T_R) \\ -\log_{10}(T_G) \\ -\log_{10}(T_B) \end{bmatrix}

在彩色负片工艺中,未经曝光但经过标准 C-41 工艺显影的胶卷边缘存在由未反应的彩色成色剂构成的橙黄色基底,即橙色色罩。该基底引入了显著的恒定偏置。设同卷同批次测得的未曝光片基密度为 DbaseD_{\text{base}},则净光学密度 DnetD_{\text{net}} 为原始采集密度与片基密度的逐通道代数差:

Dnet=Draw−Dbase=[Draw,R−Dbase,RDraw,G−Dbase,GDraw,B−Dbase,B]D_{\text{net}} = D_{\text{raw}} - D_{\text{base}} = \begin{bmatrix} D_{\text{raw}, R} - D_{\text{base}, R} \\ D_{\text{raw}, G} - D_{\text{base}, G} \\ D_{\text{raw}, B} - D_{\text{base}, B} \end{bmatrix}

在净光学密度域内,负片的显影产物与原始景物曝光量呈现出高度规律的反相响应。若胶片乳剂未进入饱和曝光肩部或严重欠曝趾部,净密度向量与感光卤化银晶体受光量 EE 的对数满足经典 Hurter-Driffield 关系:

Dnet,c≈γc⋅log⁡10(Ec)+Cc,c∈{R,G,B}D_{\text{net}, c} \approx \gamma_c \cdot \log_{10}(E_c) + C_c, \quad c \in \{R, G, B\}

其中 γc\gamma_c 为该染料层的显影反差系数。在光学密度域中,暗房工艺的核心操作得以在简单的代数加减法中完成:

密度域核心运算 对应物理现象与暗房操作 运算数学表达式
片基密度扣除 消除片基聚酯/三醋酸片基与未曝光橙色色罩吸收 Dnet=Draw−DbaseD_{\text{net}} = D_{\text{raw}} - D_{\text{base}}
印片光对数增益 模拟光学放大机色头红绿蓝加色/减色滤镜调节 Dgraded=Dnet+ΔPlightD_{\text{graded}} = D_{\text{net}} + \Delta P_{\text{light}}
密度动态范围缩放 将胶片记录的高反差密度映射至显示参考窗口 Dnorm=(Dhigh−D)/(Dhigh−Dlow)D_{\text{norm}} = (D_{\text{high}} - D) / (D_{\text{high}} - D_{\text{low}})
染料串扰去耦合 校正品红染料副吸收导致的高密度区域偏青发暗 Dcalibrated=Mcrosstalk⋅DnetD_{\text{calibrated}} = M_{\text{crosstalk}} \cdot D_{\text{net}}

# 正片与输出域:从光学密度到工作色域与显示呈现

经过片基扣除与印片光调整后的胶片净密度,依然属于物理吸收测量空间,而不是适合人眼观看或进行数字艺术创作的正片图像。正片与输出域的目标,是将经过反相映射后的相对曝光量转换到具有严格色度学定义的线性工作色彩空间,并最终通过色调映射输出给目标设备。

NexFilm 引擎采用线性 ProPhoto RGB 作为统一的内部正片工作空间。将归一化后的正片光强信号映射到该空间后,系统在此基底上展开后续的调色与交付处理:

  1. 工作空间线性化:将归一化对数密度转化为正片线性相对曝光 Epos=10−α⋅DgradedE_{\text{pos}} = 10^{-\alpha \cdot D_{\text{graded}}},使得正片信号具备与真实物理场景一致的辐射线性特性。
  2. 分通道与合成色调曲线:红、绿、蓝分通道曲线优先作用于各独立色阶,合成 RGB 曲线随后作用于整体亮度,模拟传统暗房相纸的感光反差特性曲线。
  3. 3D 相纸模拟查找表:导入经过严格光谱拟合的经典彩色印相纸(如 Kodak 2383、Kodak Vision Premier 2393)三维 LUT,赋予数字反相作品深沉的黑色层次与丰富的染料光泽。
  4. 输出色彩空间转换与 OETF 编码:经由矩阵变换适配目标色彩空间(如 sRGB、Display P3、Adobe RGB 或 ACEScg),并施加目标标准的光电传递函数(OETF)与色域量化裁切。
注意

在正片输出域中进行的饱和度微调、反差增强以及高光阴影压缩,属于感知色调调整范畴。这些操作绝不应回溯反噬光学密度的测量基底。

# 非线性不可交换性定理:矩阵作用域的严谨数学证明

在负片数字处理管线中,一个致命的架构错误是混淆透射率域矩阵与密度域矩阵。常有经验式软件直接在密度域施加普通显示器色彩空间转换矩阵,声称其等价于硬件分光校准。本节通过严格的泛函代数推导证明该做法在物理和数学上均属谬误。

设三通道输入透射率向量为 x=[T1,T2,T3]T∈R>03x = [T_1, T_2, T_3]^T \in \mathbb{R}_{>0}^3。设采集设备的光学通道分离校准矩阵为 A∈R3×3A \in \mathbb{R}^{3 \times 3},其元素为 aija_{ij}。

在透射率域,硬件校准矩阵的作用是混合光通量,校正滤镜串扰:

ytrans=Ax,yi=∑j=13aijxjy_{\text{trans}} = A x, \quad y_i = \sum_{j=1}^3 a_{ij} x_j

若在此之后施加对数变换进入密度域,所获得的密度响应向量为:

Dcorrect=−log⁡10(Ax)=[−log⁡10(∑j=13a1jxj)−log⁡10(∑j=13a2jxj)−log⁡10(∑j=13a3jxj)]D_{\text{correct}} = -\log_{10}(A x) = \begin{bmatrix} -\log_{10}\left(\sum_{j=1}^3 a_{1j} x_j\right) \\ -\log_{10}\left(\sum_{j=1}^3 a_{2j} x_j\right) \\ -\log_{10}\left(\sum_{j=1}^3 a_{3j} x_j\right) \end{bmatrix}

假定某种错误实现将矩阵操作挪移到对数之后,即先计算各通道独立密度 Draw=−log⁡10(x)D_{\text{raw}} = -\log_{10}(x),再应用矩阵 AA:

Dwrong=A[−log⁡10(x)]=[∑j=13a1j[−log⁡10(xj)]∑j=13a2j[−log⁡10(xj)]∑j=13a3j[−log⁡10(xj)]]=[−log⁡10(∏j=13xja1j)−log⁡10(∏j=13xja2j)−log⁡10(∏j=13xja3j)]D_{\text{wrong}} = A [-\log_{10}(x)] = \begin{bmatrix} \sum_{j=1}^3 a_{1j} [-\log_{10}(x_j)] \\ \sum_{j=1}^3 a_{2j} [-\log_{10}(x_j)] \\ \sum_{j=1}^3 a_{3j} [-\log_{10}(x_j)] \end{bmatrix} = \begin{bmatrix} -\log_{10}\left(\prod_{j=1}^3 x_j^{a_{1j}}\right) \\ -\log_{10}\left(\prod_{j=1}^3 x_j^{a_{2j}}\right) \\ -\log_{10}\left(\prod_{j=1}^3 x_j^{a_{3j}}\right) \end{bmatrix}

对比两式的真数部分,欲使 Dcorrect≡DwrongD_{\text{correct}} \equiv D_{\text{wrong}} 恒成立,必须对所有定义域内的 xx 满足:

∑j=13aijxj≡∏j=13xjaij,∀i∈{1,2,3}\sum_{j=1}^3 a_{ij} x_j \equiv \prod_{j=1}^3 x_j^{a_{ij}}, \quad \forall i \in \{1, 2, 3\}

对该恒等式两端求关于 xkx_k 的偏导数:

左端对 xkx_k 的一阶偏导数为常数:

∂∂xk(∑j=13aijxj)=aik\frac{\partial}{\partial x_k} \left(\sum_{j=1}^3 a_{ij} x_j\right) = a_{ik}

右端对 xkx_k 的一阶偏导数为:

∂∂xk(∏j=13xjaij)=aikxkaik−1∏j≠kxjaij\frac{\partial}{\partial x_k} \left(\prod_{j=1}^3 x_j^{a_{ij}}\right) = a_{ik} x_k^{a_{ik}-1} \prod_{j \neq k} x_j^{a_{ij}}

若使两端偏导数对所有 x∈R>03x \in \mathbb{R}_{>0}^3 恒等成立,矩阵 AA 必须且仅能退化为单位矩阵 II 或仅包含对角元素的缩放矩阵。一旦矩阵含有非对角互串扰元素(即存在 j≠kj \neq k 使得 aij≠0a_{ij} \neq 0),该恒等式必然破裂:

−log⁡10(Ax)≠A[−log⁡10(x)]-\log_{10}(A x) \neq A [-\log_{10}(x)]

结论:光通量的物理线性叠加发生在对数变换之前;染料吸光度的物理代数叠加发生对数变换之后。试图在密度域施加普通色度矩阵来代替透射率域分光校正,不仅无法校正硬件色偏,反而在高光与阴影过渡区域引入严重的伪色与非线性畸变。

# 信号流图与 NexFilm 真实物理能力边界

为了让开发者与专业调色师清晰掌握数据在引擎内的生命周期,下图展示了 NexFilm v1.1.0 的端到端信号流图。图中用不同颜色明确标识了三大数学域的分界:

透射率域 物理光子通量 · 线性辐射度 RAW / Scanner RGB 黑电平解马赛克 / 坏点消除 暗场扣除与平场归一化 T = (I - dark) / (open - dark) Capture Separation 3x3 相机x背光透射率通道解耦 对数变换门限钳位 D_raw = -log10(max(T, eps)) 光学密度域 比尔-朗伯吸光度 · 加性代数空间 片基采样与色罩扣除 D_net = D_raw - D_base 暗房印片光微调 Printer Lights 对数加减色平移 有界通道跨度响应 1.15x死区 / 1.35x上限 / 3x截断 密度反相与曝光归一化 E_pos = 10^(-alpha * D) 正片与输出域 (Print/Output) Linear ProPhoto · 色彩空间变换 Linear ProPhoto 统一基底 无损正片线性光强中间表示 四通道曲线与印相纸 LUT RGB分离曲线与 Kodak 2383 等 细节增强与色调压缩 高反差边缘锐化 / 纹理 / 阴影 显示与母带导出 sRGB / P3 / ACEScg / 16-bit TIFF

最后,本白皮书在此郑重界定 NexFilm v1.1.0 已落地能力与未来规划边界,拒绝任何技术浮夸:

  1. 已落地能力:三大数学域的彻底物理隔离;暗场与开场平场相对透射率归一化;透射色卡采样拟合与条件数安全回退;基于有界通道跨度的抗偏色统一管线;16-bit 线性母带与相机 RAW 原样重封装 DNG 导出。
  2. 当前技术边界:当前版本未内置物理 Status M 滤光片加权积分测定库,status_m_crosstalk_matrix() 属于历史经验参考;未实现基于具体胶片乳剂批次物理显影参数的非线性 H-D 曲线反演。所有反相计算均基于实测透射率与严密相对密度代数模型。