大气散射模型
理论模型
在雾天拍照时，探测器接收到的光源会被雾干扰，此时的收集到的光源主要来自两部分
- 目标反射光经过粒子衰减到达探测系统的光
- 光源经过粒子散射形成的大气光
数学模型
I
(
x
,
λ
)
=
e
−
β
(
λ
)
d
(
x
)
R
(
x
,
λ
)
+
L
∞
(
1
−
e
−
β
(
λ
)
d
(
x
)
)
=
D
(
x
,
λ
)
+
A
(
x
,
λ
)
I(x, \lambda) = e^{-\beta(\lambda)d(x)}R(x, \lambda) + L_{\infty}(1-e^{-\beta(\lambda)d(x)}) = D(x, \lambda) + A(x, \lambda)
I(x,λ)=e−β(λ)d(x)R(x,λ)+L∞(1−e−β(λ)d(x))=D(x,λ)+A(x,λ)
- x x x表示像素点的位置， λ \lambda λ表示光的波长， L ∞ L_{\infty} L∞表示无穷远处的大气光值
- I ( x , λ ) I(x, \lambda) I(x,λ)：探测系统所获得的雾天图像
- R ( x , λ ) R(x, \lambda) R(x,λ)：要恢复的无雾图像
- e − β ( λ ) d ( x ) e^{-\beta(\lambda)d(x)} e−β(λ)d(x)：传输函数，也就是真实的目标图将要受到的粒子的影响，物理意义为经过粒子衰减能够达到探测系统的那部分光的比例
在图像处理中，大部分通过探测系统获得含雾图像并对其进行去雾图像处理均以上述公式为理论模型，思路是根据先验知识或者图像处理手段，从图像中估计传输函数 e − β ( λ ) d ( x ) e^{-\beta(\lambda)d(x)} e−β(λ)d(x)或者大气光 A ( x , λ ) A(x, \lambda) A(x,λ)，最后将其带入模型中求解目标图 R ( x , λ ) R(x, \lambda) R(x,λ)
引起探测系统成像质量下降的原因
- 目标反射光在传输过程中受到介质中悬浮粒子的吸收和散射作用，导致能量衰减：降低图像亮度，图像对比度下降
- 太阳光、天空光等环境光受到介质中粒子的散射作用形成杂散光：图像模糊，图像色彩不自然
公式推导
1. 平行光
反射光在传播的过程中，随着传输距离的增加光强逐渐衰减，假设光束具有单位横截面积，从
x
x
x位置，每传输一段距离
d
x
dx
dx强度变化为
d
E
(
x
,
λ
)
E
(
x
,
λ
)
=
−
β
(
λ
)
d
x
\frac{dE(x, \lambda)}{E(x, \lambda)} = -\beta(\lambda)dx
E(x,λ)dE(x,λ)=−β(λ)dx
- E ( x , λ ) E(x, \lambda) E(x,λ)表示衰减后的光强
- β ( λ ) \beta(\lambda) β(λ)是散射系数，描述介质对不同波长光的散射能力
当输入光为平行光，且未衰减光束光强为
E
0
(
λ
)
E_{0}(\lambda)
E0(λ)时，对上述微分方程从
x
=
0
→
x
=
d
x=0\rightarrow x=d
x=0→x=d积分，即可得到平行光束在
x
=
d
x=d
x=d处衰减后的光强
E
(
d
,
λ
)
=
E
0
(
λ
)
e
−
β
(
λ
)
d
E(d, \lambda) = E_{0}(\lambda)e^{-\beta(\lambda)d}
E(d,λ)=E0(λ)e−β(λ)d
具体计算过程：
d
E
(
x
,
λ
)
E
(
x
,
λ
)
=
−
β
(
λ
)
d
x
→
两
边
对
x
积
分
\frac{dE(x, \lambda)}{E(x, \lambda)} = -\beta(\lambda)dx~~~\rightarrow 两边对x积分
E(x,λ)dE(x,λ)=−β(λ)dx →两边对x积分
d
E
(
x
,
λ
)
E
(
x
,
λ
)
=
−
β
(
λ
)
d
x
∫
0
d
d
E
(
x
,
λ
)
E
(
x
,
λ
)
=
∫
0
d
−
β
(
λ
)
d
x
l
n
E
(
x
,
λ
)
∣
x
=
0
x
=
d
=
−
β
(
λ
)
x
∣
x
=
0
x
=
d
l
n
E
(
d
,
λ
)
E
0
(
λ
)
=
−
β
(
λ
)
d
E
(
d
,
λ
)
=
E
0
(
λ
)
e
−
β
(
λ
)
d
\begin{aligned} \frac{dE(x, \lambda)}{E(x, \lambda)} & = -\beta(\lambda)dx \\ \int_{0}^{d} \frac{dE(x, \lambda)}{E(x, \lambda)} & = \int_{0}^{d} -\beta(\lambda)dx \\ lnE(x, \lambda)|_{x=0}^{x=d} & = -\beta(\lambda)x|_{x=0}^{x=d} \\ ln\frac{E(d, \lambda)}{E_{0}(\lambda)} & = -\beta(\lambda)d \\ E(d, \lambda) &= E_{0}(\lambda)e^{-\beta(\lambda)d} \end{aligned}
E(x,λ)dE(x,λ)∫0dE(x,λ)dE(x,λ)lnE(x,λ)∣x=0x=dlnE0(λ)E(d,λ)E(d,λ)=−β(λ)dx=∫0d−β(λ)dx=−β(λ)x∣x=0x=d=−β(λ)d=E0(λ)e−β(λ)d
2. 点光源
当输入光为点光源，且未衰减光束光强为
I
0
I_{0}
I0时，对微分方程从
x
=
0
→
x
=
d
x=0\rightarrow x=d
x=0→x=d积分，即可得到点光源在
x
=
d
x=d
x=d处衰减后的光强
E
(
d
,
λ
)
=
I
0
(
λ
)
e
−
β
(
λ
)
d
d
2
=
L
∞
ρ
(
x
)
d
2
e
−
β
(
λ
)
d
E(d, \lambda) = \frac{I_{0}(\lambda)e^{-\beta(\lambda)d}}{d^{2}} = \frac{L_{\infty}\rho (x)}{d^{2}}e^{-\beta(\lambda)d}
E(d,λ)=d2I0(λ)e−β(λ)d=d2L∞ρ(x)e−β(λ)d
探测器接受到的大气光成分主要包括太阳直射光、大气漫反射光和地面反射光。如下图所示，体积微元
d
V
=
d
w
x
2
d
x
dV = dwx^{2}dx
dV=dwx2dx内的介质被看做成一个光源，强度为
d
I
(
x
,
λ
)
=
d
V
k
β
(
λ
)
dI(x, \lambda) = dVk\beta (\lambda)
dI(x,λ)=dVkβ(λ)，
k
k
k为光源常数。根据点光源的损耗公式，到达探测器后的光强为
d
L
(
x
,
λ
)
=
d
I
(
x
,
λ
)
e
−
β
(
λ
)
x
d
w
x
2
=
k
β
(
λ
)
e
−
β
(
λ
)
x
d
x
dL(x, \lambda) = \frac{dI(x, \lambda)e^{-\beta(\lambda)x}}{dwx^{2}} = k\beta(\lambda)e^{-\beta(\lambda)x}dx
dL(x,λ)=dwx2dI(x,λ)e−β(λ)x=kβ(λ)e−β(λ)xdx，在
x
=
0
→
x
=
d
x=0\rightarrow x=d
x=0→x=d进行积分可以得到总的大气光强值
L
(
d
,
λ
)
=
k
(
1
−
e
−
β
(
λ
)
d
)
L(d, \lambda) = k(1-e^{-\beta(\lambda)d})
L(d,λ)=k(1−e−β(λ)d)，由于光源来自无穷远处的太空，所以
k
k
k表示无穷远处的大气光强值，令
k
=
L
∞
(
λ
)
k = L_{\infty}(\lambda)
k=L∞(λ)，大气光强值表示为
L
(
d
,
λ
)
=
L
∞
(
1
−
e
−
β
(
λ
)
d
)
L(d, \lambda) = L_{\infty}(1-e^{-\beta(\lambda)d})
L(d,λ)=L∞(1−e−β(λ)d)
成像总光强为
I
(
x
)
=
L
∞
ρ
(
x
)
d
2
e
−
ρ
d
(
x
)
+
L
∞
(
1
−
e
−
β
d
(
x
)
)
I(x) = \frac{L_{\infty}\rho(x)}{d^{2}}e^{-\rho d(x)} + L_{\infty}(1-e^{-\beta d(x)})
I(x)=d2L∞ρ(x)e−ρd(x)+L∞(1−e−βd(x))
为方便计算，令大气透射率为
t
(
x
)
=
e
−
ρ
d
(
x
)
t(x) = e^{-\rho d(x)}
t(x)=e−ρd(x)，目标反射光为
J
(
x
)
=
L
∞
ρ
(
x
)
d
2
J(x) = \frac{L_{\infty}\rho(x)}{d^{2}}
J(x)=d2L∞ρ(x)，目标衰减反射光为
D
(
x
)
=
J
(
x
)
t
D(x) = J(x)t
D(x)=J(x)t，大气光为
A
=
L
∞
(
1
−
t
)
=
A
∞
(
1
−
t
)
A = L_{\infty}(1-t) = A_{\infty}(1-t)
A=L∞(1−t)=A∞(1−t)，最终大气散射模型的数学表达式为
I
(
x
)
=
D
+
A
=
J
(
x
)
t
(
x
)
+
A
∞
(
1
−
t
(
x
)
)
I(x) = D+A = J(x)t(x) + A_{\infty}(1-t(x))
I(x)=D+A=J(x)t(x)+A∞(1−t(x))
5568
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
