YOLA# 朗伯假设<根据双色反射模型的体反射项>根据双色反射模型的体反射项>, C p i C_{p_i} C p i 的值可用离散形式表示为:
C p i = m ( n ⃗ p i , l ⃗ p i ) e C p i ( λ ) ρ C p i ( λ ) C_{p_i}=m(\vec{n}_{p_i},\vec{l}_{p_i})e^{C_{p_i}}(\lambda)\rho^{C_{p_i}}(\lambda) C p i = m ( n p i , l p i ) e C p i ( λ ) ρ C p i ( λ ) 交叉颜色比率<考虑两个相邻像素>考虑两个相邻像素>,分别表示为p 1 p_{1} p 1 和p 2 p_{2} p 2 ,以及红色( R ) (R) ( R ) 和蓝色( B ) (B) ( B ) 通道,我 们可以通过以下计算步骤确定红蓝通道之间的比率M r b M_{rb} M r b :
M r b = R p 1 B p 2 R p 2 B p 1 M_{rb}=\frac{R_{p_1}B_{p_2}}{R_{p_2}B_{p_1}} M r b = R p 2 B p 1 R p 1 B p 2 对 M r b M_{rb} M r b 取对数并用公式C p i = m ( n ⃗ p i , l ⃗ p i ) e C p i ( λ ) ρ C p i ( λ ) C_{p_i}=m(\vec{n}_{p_i},\vec{l}_{p_i})e^{C_{p_i}}(\lambda)\rho^{C_{p_i}}(\lambda) C p i = m ( n p i , l p i ) e C p i ( λ ) ρ C p i ( λ ) 替换像素值后,我们得到:
l o g ( M r b ) = l o g ( m ( n p 1 ⃗ , l p 1 ⃗ ) ) − l o g ( m ( n p 1 ⃗ , l p 1 ⃗ ) ) + l o g ( e R p 1 ( λ ) ) − l o g ( e R p 2 ( λ ) ) + l o g ( ρ R p 1 ( λ ) ) − l o g ( ρ R p 2 ( λ ) ) + l o g ( m ( n p 2 ⃗ , l p 2 ⃗ ) ) − l o g ( m ( n p 2 ⃗ , l p 2 ⃗ ) ) + l o g ( e B p 2 ( λ ) ) − l o g ( e B p 1 ( λ ) ) + l o g ( ρ B p 2 ( λ ) ) − l o g ( ρ B p 1 ( λ ) ) \begin{aligned}\begin{aligned}log(M_{rb})=log(m(\vec{n_{p_1}},\vec{l_{p_1}}))-log(m(\vec{n_{p_1}},\vec{l_{p_1}}))\end{aligned}\\+log(e^{R_{\boldsymbol{p}_1}}(\lambda))-log(e^{R_{\boldsymbol{p}_2}}(\lambda))\\+log(\rho^{R_{\boldsymbol{\underset{1}{p}}}}(\lambda))-log(\rho^{R_{\boldsymbol{\underset{2}{p}}}}(\lambda))\\+log(m(\vec{n_{p_2}},\vec{l_{p_2}}))-log(m(\vec{n_{p_2}},\vec{l_{p_2}}))\\+log(e^{B_{\boldsymbol{p}_2}}(\lambda))-log(e^{B_{\boldsymbol{p}_1}}(\lambda))\\+log(\rho^{B_{p_2}}(\lambda))-log(\rho^{B_{p_1}}(\lambda))\end{aligned} l o g ( M r b ) = l o g ( m ( n p 1 , l p 1 )) − l o g ( m ( n p 1 , l p 1 )) + l o g ( e R p 1 ( λ )) − l o g ( e R p 2 ( λ )) + l o g ( ρ R 1 p ( λ )) − l o g ( ρ R 2 p ( λ )) + l o g ( m ( n p 2 , l p 2 )) − l o g ( m ( n p 2 , l p 2 )) + l o g ( e B p 2 ( λ )) − l o g ( e B p 1 ( λ )) + l o g ( ρ B p 2 ( λ )) − l o g ( ρ B p 1 ( λ )) 基于e C p 1 ≈ e C p 2 e^{C_{p_1}}\approx e^{C_{p_2}} e C p 1 ≈ e C p 2 的光照假设,上述方程可进一步简化为光照不变形式:
l o g ( M r b ) = l o g ( ρ R p 1 ( λ ) ) − l o g ( ρ R p 2 ( λ ) ) + l o g ( ρ B p 2 ( λ ) ) − l o g ( ρ B p 1 ( λ ) ) \begin{aligned}log(M_{rb})&=log(\rho^{R_{p_1}}(\lambda))-log(\rho^{R_{p_2}}(\lambda))\\&+log(\rho^{B_{p_2}}(\lambda))-log(\rho^{B_{p_1}}(\lambda))\end{aligned} l o g ( M r b ) = l o g ( ρ R p 1 ( λ )) − l o g ( ρ R p 2 ( λ )) + l o g ( ρ B p 2 ( λ )) − l o g ( ρ B p 1 ( λ )) 可学习核。目标是将固定的光照不变特征转化为可学习形式。具体而言,我们旨在学习一 组卷积核W 1 , W 2 , ⋯ W n ∈ k × k \mathcal{W}_1,\mathcal{W}_2,\cdots\mathcal{W}_n^{\in k\times k} W 1 , W 2 , ⋯ W n ∈ k × k ,其中n n n 表示核的数量, k k k 表示核大小。此处,我们将固定特 征扩展为更具通用性和泛化性的形式。设p i p_{i} p i 和w i w_{i} w i 表示核W n \mathcal{W}_{n} W n 内的一组像素位置及其对应权 重,其中i = 0 , 1 , ⋯ k 2 i=0,1,\cdots k^2 i = 0 , 1 , ⋯ k 2 。这些参数使我们能够将交叉颜色比率 (CCR) 演变为可适应形式, 从而提升其有效处理不同光照条件的能力。请注意w i w_{i} w i 是可训练的,这使得正负极性变得无关紧要
M r b = ∏ ( i , j = 1 i ≠ j ) k 2 ( R p i B p i ) w i ( B p j R p j ) w j M_{rb}=\prod_{\binom{i,j=1}{i\neq j}}^{k^2}\left(\frac{R_{p_i}}{B_{p_i}}\right)^{w_i}\left(\frac{B_{p_j}}{R_{p_j}}\right)^{w_j} M r b = ( i = j i , j = 1 ) ∏ k 2 ( B p i R p i ) w i ( R p j B p j ) w j 为使扩展形式仍满足光照不变性, M r b M_{rb} M r b 的对数需满足以下约束条件:
{ ∑ i k 2 w i l o g ( e R p i ( λ ) ) = 0 ∑ i k 2 w i l o g ( e B p i ( λ ) ) = 0 \begin{cases}\sum_i^{k^2}w_ilog(e^{R_{p_i}}(\lambda))=0\\\sum_i^{k^2}w_ilog(e^{B_{p_i}}(\lambda))=0&&\end{cases} { ∑ i k 2 w i l o g ( e R p i ( λ )) = 0 ∑ i k 2 w i l o g ( e B p i ( λ )) = 0 如果上述等式成立,e 项和 m 项将被消除。最终特征可以用以下广义形式表示:
log ( M r b ) = ∑ i k 2 w i log ( ρ R p i ( λ ) ) − ∑ i k 2 w i log ( ρ B p i ( λ ) ) \log(M_{rb})=\sum_i^{k^2}w_i\log(\rho^{R_{p_i}}\left(\lambda\right))-\sum_i^{k^2}w_i\log(\rho^{B_{p_i}}\left(\lambda\right)) log ( M r b ) = i ∑ k 2 w i log ( ρ R p i ( λ ) ) − i ∑ k 2 w i log ( ρ B p i ( λ ) ) 将核 W i \mathcal{W}_i W i 应用于图像 I I I 所得到的特征,记为 f W i ( I ) f_{\mathcal{W}_i}(I) f W i ( I ) ,可以表示为:
f W i ( I ) = [ W i ⊛ log ( R ) + ( − W i ) ⊛ log ( B ) W i ⊛ log ( R ) + ( − W i ) ⊛ log ( G ) W i ⊛ log ( G ) + ( − W i ) ⊛ log ( B ) ] f_{\mathcal{W}_i}(I)=\left[\begin{array}{c}\mathcal{W}_i\circledast \log(R)+(-\mathcal{W}_i)\circledast \log(B)\\\mathcal{W}_i\circledast \log(R)+(-\mathcal{W}_i)\circledast \log(G)\\\mathcal{W}_i\circledast \log(G)+(-\mathcal{W}_i)\circledast \log(B)\end{array}\right] f W i ( I ) = W i ⊛ log ( R ) + ( − W i ) ⊛ log ( B ) W i ⊛ log ( R ) + ( − W i ) ⊛ log ( G ) W i ⊛ log ( G ) + ( − W i ) ⊛ log ( B ) 零均值约束(Zero mean constraint) :根据公式 { ∑ i k 2 w i log ( e R p i ( λ ) ) = 0 ∑ i k 2 w i log ( e B p i ( λ ) ) = 0 \begin{cases}\sum_i^{k^2}w_i\log(e^{R_{p_i}}(\lambda))=0\\\sum_i^{k^2}w_i\log(e^{B_{p_i}}(\lambda))=0&&\end{cases} { ∑ i k 2 w i log ( e R p i ( λ )) = 0 ∑ i k 2 w i log ( e B p i ( λ )) = 0 以及近似 e R p i ≈ e B p i e^{R_{p_{i}}}\approx e^{B_{p_{i}}} e R p i ≈ e B p i ,在卷积核的语境下,我们只需确保 W n ∈ k × k \mathcal{W}_n^{\in k\times k} W n ∈ k × k 的均值为 0,如下所示:
W n ‾ = 1 k 2 ∑ i = 1 k 2 w i = 0 \overline{\mathcal{W}_n}=\frac{1}{k^2}\sum_{i=1}^{k^2}w_i=0 W n = k 2 1 i = 1 ∑ k 2 w i = 0
FRBNet# 受 Phong 光照模型中加性分解(additive decomposition)的启发,我们引入了一种适应真实低光场景的朗伯模型(Lambertian model)扩展版本。我们将局部光源重新解释为非均匀高光,其表示如下:
I C ( x , y ) = m [ n ⃗ ( x , y ) , l ⃗ ( x , y ) ] ⋅ φ C ( x , y ) ⋅ ρ C ( x , y ) + S C ( x , y ) , (2) I_C(x, y) = m[\vec{n}(x, y), \vec{l}(x, y)] \cdot \varphi_C(x, y) \cdot \rho_C(x, y) + S_C(x, y), \tag{2} I C ( x , y ) = m [ n ( x , y ) , l ( x , y )] ⋅ φ C ( x , y ) ⋅ ρ C ( x , y ) + S C ( x , y ) , ( 2 ) 其中 S C S_C S C 代表空间不规则的高光分量,可进一步定义为:
S C ( x , y ) = H C ( x , y ) ⋅ m [ n ⃗ ( x , y ) , l ⃗ ( x , y ) ] ⋅ φ C ( x , y ) ⋅ ρ C ( x , y ) , (3) S_C(x, y) = H_C(x, y) \cdot m[\vec{n}(x, y), \vec{l}(x, y)] \cdot \varphi_C(x, y) \cdot \rho_C(x, y), \tag{3} S C ( x , y ) = H C ( x , y ) ⋅ m [ n ( x , y ) , l ( x , y )] ⋅ φ C ( x , y ) ⋅ ρ C ( x , y ) , ( 3 ) 这里 H C H_C H C 表示高光干扰的相对强度。为简化符号,我们定义 D C ( x , y ) = m [ n ⃗ ( x , y ) , l ⃗ ( x , y ) ] ⋅ φ C ( x , y ) ⋅ ρ C ( x , y ) D_C(x, y) = m[\vec{n}(x, y), \vec{l}(x, y)] \cdot \varphi_C(x, y) \cdot \rho_C(x, y) D C ( x , y ) = m [ n ( x , y ) , l ( x , y )] ⋅ φ C ( x , y ) ⋅ ρ C ( x , y ) 为标准漫反射分量。将其代入公式 (2) 并重新排列各项,我们得到一个更简洁的表达式:
I C ( x , y ) = D C ( x , y ) + S C ( x , y ) = D C ( x , y ) ⋅ ( 1 + H C ( x , y ) ) . (4) I_C(x, y) = D_C(x, y) + S_C(x, y) = D_C(x, y) \cdot (1 + H_C(x, y)). \tag{4} I C ( x , y ) = D C ( x , y ) + S C ( x , y ) = D C ( x , y ) ⋅ ( 1 + H C ( x , y )) . ( 4 ) 利用通道比(Channel Ratios, CR)来分离光照不变特征已被证明对低光视觉任务有效 [44, 17, 5]。以红通道 R R R 和绿通道 G G G 之间的通道比为例,根据我们的扩展广义低光模型,其对数变换公式可表示为:
CR R G = log ( I R I G ) = log ( φ R ⋅ ρ R ⋅ ( 1 + H R ) φ G ⋅ ρ G ⋅ ( 1 + H G ) ) = log φ R − log φ G + log ρ R − log ρ G + log ( 1 + H R ) − log ( 1 + H G ) . (5) \begin{aligned} \text{CR}_{RG} &= \log \left(\frac{I_R}{I_G}\right) = \log \left(\frac{\varphi_R \cdot \rho_R \cdot (1 + H_R)}{\varphi_G \cdot \rho_G \cdot (1 + H_G)}\right) \\ &= \log \varphi_R - \log \varphi_G + \log \rho_R - \log \rho_G + \log(1 + H_R) - \log(1 + H_G). \end{aligned} \tag{5} CR RG = log ( I G I R ) = log ( φ G ⋅ ρ G ⋅ ( 1 + H G ) φ R ⋅ ρ R ⋅ ( 1 + H R ) ) = log φ R − log φ G + log ρ R − log ρ G + log ( 1 + H R ) − log ( 1 + H G ) . ( 5 ) 如公式 (5) 所示,来自高光项的非线性残差破坏了光照和反射率的清晰分离,限制了空间域通道比方法的有效性。为了克服这些限制,我们将分析转移到频域。在频域中,光照和反射分量自然地占据不同的频带 [60],从而能够更有效地分离光照不变特征。受先前空间域通道比工作 [44, 17, 5] 的启发,我们创新性地提出了频域通道比(Frequency-domain Channel Ratio, FCR) :
FCR R G = F [ log ( I R I G ) ] = F [ log φ R − log φ G ] + F [ log ρ R − log ρ G ] + F [ log ( 1 + H R ) − log ( 1 + H G ) ] , (6) \begin{aligned} \text{FCR}_{RG} &= \mathcal{F}[\log(\frac{I_R}{I_G})] \\ &= \mathcal{F}[\log \varphi_R - \log \varphi_G] + \mathcal{F}[\log \rho_R - \log \rho_G] + \mathcal{F}[\log(1 + H_R) - \log(1 + H_G)], \end{aligned} \tag{6} FCR RG = F [ log ( I G I R )] = F [ log φ R − log φ G ] + F [ log ρ R − log ρ G ] + F [ log ( 1 + H R ) − log ( 1 + H G )] , ( 6 ) 其中 F [ ⋅ ] \mathcal{F}[\cdot] F [ ⋅ ] 代表傅里叶变换算子。为了处理非线性残差项 Δ = F [ log ( 1 + H R ) − log ( 1 + H G ) ] \Delta = \mathcal{F}[\log(1 + H_R) - \log(1 + H_G)] Δ = F [ log ( 1 + H R ) − log ( 1 + H G )] ,我们应用了一阶泰勒展开(first-order Taylor expansion)。鉴于数据中的显著贡献通常是稀疏且局部的,我们假设 H C ∈ [ 0 , 1 ) H_C \in [0, 1) H C ∈ [ 0 , 1 ) 具有相对较小的幅度,允许我们将 log ( 1 + H C ) \log(1 + H_C) log ( 1 + H C ) 近似为 H C + O ( H C 2 ) H_C + \mathcal{O}(H_C^2) H C + O ( H C 2 ) 。
在上述假设下,通过忽略高阶项,我们可以得到 Δ \Delta Δ 的线性化近似如下:
Δ = F [ H R − H G ] = H R − H G , (7) \Delta = \mathcal{F}[H_R - H_G] = \mathcal{H}_R - \mathcal{H}_G, \tag{7} Δ = F [ H R − H G ] = H R − H G , ( 7 ) 其中 a R , a G a_R, a_G a R , a G 代表幅度项,θ R , θ G \theta_R, \theta_G θ R , θ G 表示相位分量。为了表征通道间的相位关系,我们引入了频域相关系数 C o r R G = e i ( θ G − θ R ) Cor_{RG} = e^{i(\theta_G - \theta_R)} C o r RG = e i ( θ G − θ R ) (推导自 A.2,见 [56]),它量化了频域中通道响应之间的角位移。这使我们能够将 Δ \Delta Δ 重写为:
Δ = e i θ R ⋅ ( a R − a G ⋅ e i ( θ G − θ R ) ) = e i θ R ⋅ ( a R − a G ⋅ C o r R G ) , (9) \Delta = e^{i\theta_R} \cdot \left(a_R - a_G \cdot e^{i(\theta_G - \theta_R)}\right) = e^{i\theta_R} \cdot (a_R - a_G \cdot Cor_{RG}), \tag{9} Δ = e i θ R ⋅ ( a R − a G ⋅ e i ( θ G − θ R ) ) = e i θ R ⋅ ( a R − a G ⋅ C o r RG ) , ( 9 ) 这种因式分解揭示了残差项被构建为一个相位调制(phase-modulated)分量,其中 e i θ R e^{i\theta_R} e i θ R 作为载波相位,而 ( a R − a G ⋅ C o r R G ) (a_R - a_G \cdot Cor_{RG}) ( a R − a G ⋅ C o r RG ) 编码了由通道间相位相关性调制的幅度差异。
最后,频域通道比的最终公式可以总结为:
FCR R G = F [ log φ R − log φ G ] ⏟ illumination (光照) + F [ log ρ R − log ρ G ] ⏟ reflectance (反射率) + e i θ R ( a R − a G ⋅ C o r R G ) ⏟ high-lit residual (高光残差) . (10) \text{FCR}_{RG} = \underbrace{\mathcal{F}[\log \varphi_R - \log \varphi_G]}_{\text{illumination (光照)}} + \underbrace{\mathcal{F}[\log \rho_R - \log \rho_G]}_{\text{reflectance (反射率)}} + \underbrace{e^{i\theta_R}(a_R - a_G \cdot Cor_{RG})}_{\text{high-lit residual (高光残差)}}. \tag{10} FCR RG = illumination ( 光照 ) F [ log φ R − log φ G ] + reflectance ( 反射率 ) F [ log ρ R − log ρ G ] + high-lit residual ( 高光残差 ) e i θ R ( a R − a G ⋅ C o r RG ) . ( 10 ) 利用谱分离的固有特性和残差项的相位调制结构,我们设计了专门的滤波策略,旨在鲁棒地提取光照不变特征,从而提高在不同光照条件下特征提取的可靠性和有效性。
频域中的光照不变特征增强过程# 为了增强光照不变特征,所提出的 FRBNet 首先将通道比的操作转换到频域。根据第 3.2 节中提出的 FCR 函数,在频域中利用通道间的关系。定义空间域中的输入图像为 I ( x , y ) \mathbf{I}(x, y) I ( x , y ) ,对于每一对通道,FCR 通过带有可学习频率参数 ( u , v ) (u, v) ( u , v ) 的频域对数差分来实现:
{ dif R G ( u , v ) = F [ log I R ( x , y ) ] − F [ log I G ( x , y ) ] dif G B ( u , v ) = F [ log I G ( x , y ) ] − F [ log I B ( x , y ) ] dif B R ( u , v ) = F [ log I B ( x , y ) ] − F [ log I R ( x , y ) ] . (11) \begin{cases} \text{dif}^{RG}(u, v) = \mathcal{F}[\log I_R(x, y)] - \mathcal{F}[\log I_G(x, y)] \\ \text{dif}^{GB}(u, v) = \mathcal{F}[\log I_G(x, y)] - \mathcal{F}[\log I_B(x, y)] \\ \text{dif}^{BR}(u, v) = \mathcal{F}[\log I_B(x, y)] - \mathcal{F}[\log I_R(x, y)]. \end{cases} \tag{11} ⎩ ⎨ ⎧ dif RG ( u , v ) = F [ log I R ( x , y )] − F [ log I G ( x , y )] dif GB ( u , v ) = F [ log I G ( x , y )] − F [ log I B ( x , y )] dif BR ( u , v ) = F [ log I B ( x , y )] − F [ log I R ( x , y )] . ( 11 ) 接下来,我们设计了一个可学习频域滤波器(Learnable Frequency-domain Filter, LFF) ,用于减少低光图像中光照和高光残差项对每一对通道鲁棒特征提取的影响。它由一个零直流频率窗口(zero-DC frequency window)和一个改进的径向基滤波器组成。频率响应特征 F inv ( u , v ) \mathbf{F}_{\text{inv}}(u, v) F inv ( u , v ) 可以表示为:
{ F inv R G ( u , v ) = L F F R G ( u , v ) ⋅ dif R G ( u , v ) F inv G B ( u , v ) = L F F G B ( u , v ) ⋅ dif G B ( u , v ) F inv B R ( u , v ) = L F F B R ( u , v ) ⋅ dif B R ( u , v ) . (12) \begin{cases} F_{\text{inv}}^{RG}(u, v) = LFF^{RG}(u, v) \cdot \text{dif}^{RG}(u, v) \\ F_{\text{inv}}^{GB}(u, v) = LFF^{GB}(u, v) \cdot \text{dif}^{GB}(u, v) \\ F_{\text{inv}}^{BR}(u, v) = LFF^{BR}(u, v) \cdot \text{dif}^{BR}(u, v). \end{cases} \tag{12} ⎩ ⎨ ⎧ F inv RG ( u , v ) = L F F RG ( u , v ) ⋅ dif RG ( u , v ) F inv GB ( u , v ) = L F F GB ( u , v ) ⋅ dif GB ( u , v ) F inv BR ( u , v ) = L F F BR ( u , v ) ⋅ dif BR ( u , v ) . ( 12 ) 然后,滤波后的频谱特征被变换回空间域。所有通道对(R & G, G & B, B & R)的结果特征被拼接在一起:
F inv ( x , y ) = Cat ( F − 1 [ F inv R G ( u , v ) ] ; F − 1 [ F inv G B ( u , v ) ] ; F − 1 [ F inv B R ( u , v ) ] ) , (13) \mathbf{F}_{\text{inv}}(x, y) = \text{Cat} \left(\mathcal{F}^{-1} \left[F_{\text{inv}}^{RG}(u, v)\right] ; \mathcal{F}^{-1} \left[F_{\text{inv}}^{GB}(u, v)\right] ; \mathcal{F}^{-1} \left[F_{\text{inv}}^{BR}(u, v)\right]\right), \tag{13} F inv ( x , y ) = Cat ( F − 1 [ F inv RG ( u , v ) ] ; F − 1 [ F inv GB ( u , v ) ] ; F − 1 [ F inv BR ( u , v ) ] ) , ( 13 ) 其中 F − 1 \mathcal{F}^{-1} F − 1 代表傅里叶逆变换,Cat 代表拼接操作。为了进一步将来自频域的增强光照不变特征与来自原始图像的空间域特征相结合,我们采用了一个参考 [5] 的通用融合模块进行整合:
F out = Conv { CB [ Cat ( CB [ F inv ( x , y ) ] ; CB [ I ( x , y ) ] ) ] } , (14) \mathbf{F}_{\text{out}} = \text{Conv} \{\text{CB} [\text{Cat} (\text{CB}[\mathbf{F}_{\text{inv}}(x, y)]; \text{CB}[\mathbf{I}(x, y)])]\}, \tag{14} F out = Conv { CB [ Cat ( CB [ F inv ( x , y )] ; CB [ I ( x , y )])]} , ( 14 ) 其中 Conv 是卷积,而 CB 是卷积后接批归一化(Batch Normalization, BN)。最后,输出特征 F out \mathbf{F}_{\text{out}} F out 被送入下游任务网络。
可学习频域滤波器 (Learnable Frequency-domain Filter)# 我们方法的核心是可学习频域滤波器(LFF),它自适应地处理频谱分量。该滤波器由两个互补的元素组成:用于衰减低频光照的零直流频率窗口(Zero-DC Frequency Window) W g \mathbf{W_g} W g ,以及用于编码谱距离和方向信息的改进径向基滤波器(Improved Radial Basis Filter) H ( u , v ) \mathbf{H}(u, v) H ( u , v ) ,其公式如下:
L F F ( u , v ) = W g ⋅ H ( u , v ) . (15) \mathbf{LFF}(u, v) = \mathbf{W_g} \cdot \mathbf{H}(u, v). \tag{15} LFF ( u , v ) = W g ⋅ H ( u , v ) . ( 15 ) 零直流频率窗口 (Zero-DC Frequency Window)。 为了在保留结构信息的同时抑制不需要的光照,我们采用了一个以频率平面原点为中心的高斯窗口:
W g ( u , v ) = exp ( − r ( u , v ) 2 σ w 2 ) , r ( u , v ) = u 2 + v 2 , (16) \mathbf{W_g}(u, v) = \exp \left( -\frac{\mathbf{r}(u, v)^2}{\sigma_w^2} \right), \quad \mathbf{r}(u, v) = \sqrt{u^2 + v^2}, \tag{16} W g ( u , v ) = exp ( − σ w 2 r ( u , v ) 2 ) , r ( u , v ) = u 2 + v 2 , ( 16 ) 其中 σ w \sigma_w σ w 是可学习的带宽参数,r ( u , v ) \mathbf{r}(u, v) r ( u , v ) 表示归一化的径向频率坐标。为了消除直流(DC)分量,显式地设定 W g ( 0 , 0 ) = 0 \mathbf{W_g}(0, 0) = 0 W g ( 0 , 0 ) = 0 ,这确保了滤波器在去除全局亮度偏差的同时,保留用于局部结构线索的中高频信息。
改进径向基滤波器 (Improved Radial Basis Filter)。 为了构建一个具有光谱自适应性和方向选择性的滤波器,我们采用了一组可学习的径向基函数(RBFs)并结合角度调制。RBF 可以捕捉频率幅度选择性,而角度项可以引入方向敏感性,从而在傅里叶域实现各向异性滤波。定义一组以预定义频率半径 μ k ∈ [ 0 , 1 ] \mu_k \in [0, 1] μ k ∈ [ 0 , 1 ] 为中心的 K K K 个径向基函数 ϕ ( u , v ) \phi(u, v) ϕ ( u , v ) :
ϕ k ( u , v ) = exp ( − ( r ( u , v ) − μ k ) 2 2 σ h 2 ) , k = [ 1 , 2 , ⋯ , K ] (17) \phi_k(u, v) = \exp \left( -\frac{(r(u, v) - \mu_k)^2}{2\sigma_h^2} \right), k = [1, 2, \cdots, K] \tag{17} ϕ k ( u , v ) = exp ( − 2 σ h 2 ( r ( u , v ) − μ k ) 2 ) , k = [ 1 , 2 , ⋯ , K ] ( 17 ) 其中 r ( u , v ) r(u, v) r ( u , v ) 是如前定义的归一化径向频率,σ h \sigma_h σ h 是所有基函数共享的可学习带宽参数。通过加权线性组合的可学习系数 a k a_k a k ,最终的径向响应为:
Φ ( u , v ) = ∑ k = 1 K a k ⋅ ϕ k ( u , v ) , k = [ 1 , 2 , ⋯ , K ] (18) \Phi(u, v) = \sum_{k=1}^{K} a_k \cdot \phi_k(u, v), k = [1, 2, \cdots, K] \tag{18} Φ ( u , v ) = k = 1 ∑ K a k ⋅ ϕ k ( u , v ) , k = [ 1 , 2 , ⋯ , K ] ( 18 ) 此外,参考第 3.2 节中的相位导向残差结构,干扰项表现出主导的方向分量。径向响应进一步通过由方向角的正弦谐波构建的角度项进行调制,以捕捉方向选择性:
M ( u , v ) = 1 + λ ⋅ ∑ n = 1 N [ cos ( n θ ( u , v ) ) + sin ( n θ ( u , v ) ) ] , θ ( u , v ) = arctan ( v u + ϵ ) , (19) M(u, v) = 1 + \lambda \cdot \sum_{n=1}^{N} [\cos(n\theta(u, v)) + \sin(n\theta(u, v))], \quad \theta(u, v) = \arctan \left(\frac{v}{u + \epsilon}\right), \tag{19} M ( u , v ) = 1 + λ ⋅ n = 1 ∑ N [ cos ( n θ ( u , v )) + sin ( n θ ( u , v ))] , θ ( u , v ) = arctan ( u + ϵ v ) , ( 19 ) 其中 N N N 是角频率的数量,λ \lambda λ 控制调制强度。最终的频域径向基滤波器响应由下式给出:
H ( u , v ) = Φ ( u , v ) ⋅ M ( u , v ) . (20) \mathbf{H}(u, v) = \Phi(u, v) \cdot M(u, v). \tag{20} H ( u , v ) = Φ ( u , v ) ⋅ M ( u , v ) . ( 20 ) 通过整合角度谐波,改进后的径向基滤波器既具有光谱局部性又具有方向响应性,能够以数据驱动的方式对齐或抑制这些定向残差,这对于在衰减结构化干扰的同时隔离光照不变特征至关重要。