辐射传输方程的分解

news2024/12/27 13:35:01

Decomposition of the Boundary Value Problem for Radiative Transfer Equation of MODIS and MISR instruments

在这里插入图片描述

0.Notions

Let L L L be the straming-collision operator, and S S S is scattering operator:
L I = Ω ⋅ ∇ I ( r , Ω ) + σ ( r , Ω ) I ( r , Ω ) S = ∫ 4 π σ s ( r , Ω ′ , Ω ) I ( r , Ω ′ ) d Ω ′ LI = \Omega\cdot \nabla I(r,\Omega)+\sigma(r,\Omega)I(r,\Omega)\\ S=\int_{4\pi}\sigma_s(r,\Omega',\Omega)I(r,\Omega')d\Omega' LI=ΩI(r,Ω)+σ(r,Ω)I(r,Ω)S=4πσs(r,Ω,Ω)I(r,Ω)dΩ
To descirbe reflective properties of boundary δ V \delta V δV, the scattering operator defined on the boundary δ V \delta V δV for radiation exiting the regine V V V is indroduced as
R I + = 1 π ∫ δ V d r B ′ ∫ n ( r B ′ ) ⋅ Ω ′ > 0 ρ B ( r B ′ , Ω ′ ; r B , Ω ) ∣ n ( r B ′ ) ⋅ Ω ′ ∣ I + ( r B ′ , Ω ′ ) d Ω ′ \mathcal{R}I^+=\frac{1}{\pi}\int_{\delta V} dr_B'\int_{n(r_B')\cdot \Omega'>0}\rho_B(r_B',\Omega';r_B,\Omega)|n(r_B')\cdot \Omega'|I^+(r_B',\Omega')d\Omega' RI+=π1δVdrBn(rB)Ω>0ρB(rB,Ω;rB,Ω)n(rB)ΩI+(rB,Ω)dΩ
Here, we use I + I^+ I+ to represent the radiance exiting and I − I^- I entering domina V. And then, the stationary radiative transfer equation could be expressed as
L I = S I + q I − = R I + q B LI = SI +q \quad I^-=\mathcal{R}I+q_B LI=SI+qI=RI+qB
And if R \mathcal{R} R is 0 and q B q_B qB is 0, the porblem is called standard porblem and the straming-collision oprator would be write as L 0 L_0 L0.

1.Black Soil Problem

The RTE is
L I = S I + q I − = R I + q B . LI = SI +q \quad I^-=\mathcal{R}I+q_B. LI=SI+qI=RI+qB.
As in Green function blog before, we have
I ( r , Ω ) = I b s ( r , Ω ) + I r e s t ( r , Ω ) I(r,\Omega) = I_{bs}(r,\Omega)+I_{rest}(r,\Omega) I(r,Ω)=Ibs(r,Ω)+Irest(r,Ω)
The first question is black soil problem, means that the bottom of canopy is non-reflecting surface (soil absorption is 1, black soil), so its the solution of
L I b s = S I b s I b s ( r t , Ω ) = c T δ ( Ω , Ω 0 ) + d ( r t , Ω ) I b s − ( r l , Ω ) = 0 I b s − ( r b , Ω ) = 0 \begin{aligned} &LI_{bs}=SI_{bs}\\ &I_{bs}^(r_t,\Omega)=c_T\delta(\Omega,\Omega_0)+d(r_t,\Omega)\\ &I_{bs}^-(r_l,\Omega)=0\\ &I_{bs}^-(r_b,\Omega)=0 \end{aligned} LIbs=SIbsIbs(rt,Ω)=cTδ(Ω,Ω0)+d(rt,Ω)Ibs(rl,Ω)=0Ibs(rb,Ω)=0
Here, we assumte the lateral effects is zero.

2.Canopy-Surface Interation

For notion, we use I r I_r Ir for the I r e s t I_{rest} Irest. We have that
L I r = S I r LI_r=SI_r LIr=SIr
and the boundary conditions is
I r − ( r t , Ω ) = 0 I r − ( r l , Ω ) = 0 I r − ( r b , Ω ) = R I + \begin{aligned} &I_r^-(r_t,\Omega)=0\\ &I_r^-(r_l,\Omega)=0\\ &I_r^-(r_b,\Omega)=\mathcal{R}I^+ \end{aligned} Ir(rt,Ω)=0Ir(rl,Ω)=0Ir(rb,Ω)=RI+
The solution I r I_r Ir of this boundary value problem depend on the solution of “complete tranport problem” I = I b s + I r I=I_{bs}+I_r I=Ibs+Ir (since R I + \mathcal{R}I^+ RI+) is also not known. It need further decompose.

The lower boundary conditions can be rewritten as
I r − ( r b , Ω ) = R I + T T I_r^-(r_b,\Omega)=\frac{\mathcal{R}I^+}{T}T Ir(rb,Ω)=TRI+T
Here, T T T​ is downward radiation flux density at canopy bottom, i.e.,
T ( r b ) = ∫ n ( r b ) ⋅ Ω ′ > 0 I ( r b , Ω ′ ) ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ ,    r b ∈ δ V T(r_b)=\int_{n(r_b)\cdot\Omega'>0}I(r_b,\Omega')|\Omega'\cdot n(r_b)|d\Omega', \ \ r_b\in\delta V T(rb)=n(rb)Ω>0I(rb,Ω)Ωn(rb)dΩ,  rbδV
The ratio R I + T \frac{\mathcal{R}I^+}{T} TRI+ is a factor of π \pi π smaller than HDRF (since there is no 1 / π 1/\pi 1/π in the dominator compared to HDRF).

And then ground BHR could be writen as
ρ e f f ( r b ) = ∫ n ( r b ) ⋅ Ω ′ > 0 ∫ n ( r b ) ⋅ Ω < 0 ∣ n ( r b ) ⋅ Ω ′ ∣ π − 1 ρ b ( r b , Ω ′ , Ω ) ∣ Ω ⋅ n ( r b ) ∣ I ( r b , Ω ′ ) d Ω ′ d Ω ∫ n ( r b ) ⋅ Ω ′ > 0 ∣ Ω ′ ⋅ n ( r b ) ∣ I ( r b , Ω ′ ) d Ω ′ \rho_{eff}(r_b)=\frac{\int_{n(r_b)\cdot\Omega'>0}\int_{n(r_b)\cdot\Omega<0}|n(r_b)\cdot\Omega'|\pi^{-1}\rho_b(r_b,\Omega',\Omega)|\Omega\cdot n(r_b)|I(r_b,\Omega')d\Omega'd\Omega}{\int_{n(r_b)\cdot\Omega'>0}|\Omega'\cdot n(r_b)|I(r_b,\Omega')d\Omega'} ρeff(rb)=n(rb)Ω>0Ωn(rb)I(rb,Ω)dΩn(rb)Ω>0n(rb)Ω<0n(rb)Ωπ1ρb(rb,Ω,Ω)∣Ωn(rb)I(rb,Ω)dΩdΩ
Here, ρ b \rho_b ρb is the BRF of canopy ground. So we could find that the denomator is actually the T T T, which means that ρ e f f = s o m e t h i n g T \rho_{eff}=\frac{something}{T} ρeff=Tsomething.

For horizontally inhomogeneous vegetation canopy, T T T vary a lot, but BHR would not so.

Now, we introduce the effective ground anisotropy:
d b ( r b , Ω ) = R I + T ρ e f f = 1 ρ e f f ⋅ ∫ 2 π + ρ b ( r b , Ω ′ , Ω ) ∣ n ( r b ) ⋅ Ω ′ ∣ I ( r b , Ω ′ ) d Ω ′ π ∫ 2 π − ∣ Ω ′ ⋅ n ( r b ) ∣ I ( r b , Ω ′ ) d Ω ′ \begin{aligned} d_b(r_b,\Omega)&=\frac{\mathcal{R}I^+}{T\rho_{eff}}\\ &=\frac{1}{\rho_{eff}}\cdot\frac{\int_{2\pi+}\rho_b(r_b,\Omega',\Omega)|n(r_b)\cdot\Omega'|I(r_b,\Omega')d\Omega'}{\pi\int_{2\pi-}|\Omega'\cdot n(r_b)|I(r_b,\Omega')d\Omega'} \end{aligned} db(rb,Ω)=TρeffRI+=ρeff1π2πΩn(rb)I(rb,Ω)dΩ2π+ρb(rb,Ω,Ω)n(rb)ΩI(rb,Ω)dΩ
This value seems strange. But if we use ρ e f f = s o m e t h i n g T \rho_{eff}=\frac{something}{T} ρeff=Tsomething​, we could find that this is actually tha ratio of reflected flux density and radiance. So its cosine-weighted intergral over upward directions is unity:
∫ 2 π + d b ( r b , Ω ) ∣ n ( r b ) ⋅ Ω ∣ d Ω = 1. \int_{2\pi+}d_b(r_b,\Omega)|n(r_b)\cdot\Omega|d\Omega=1. 2π+db(rb,Ω)n(rb)Ω∣dΩ=1.
So this term descibed the inhomogenous of the surface., its something with denonimator a flux density and numenator is reflected radiance with some direction.

By using the Equation (12), we obtain
I r − ( r b , Ω ) = R I + = ρ e f f T ( r b ) d b ( r b , Ω ) I_r^-(r_b,\Omega)=\mathcal{R}I^+=\rho_{eff}T(r_b)d_b(r_b,\Omega) Ir(rb,Ω)=RI+=ρeffT(rb)db(rb,Ω)
Using Green function method with q = I r − q=I_r^- q=Ir, we obtain
I r ( r , Ω ) = ∫ δ V d r b ′ ∫ n ( r b ′ ) ⋅ Ω ′ < 0 d Ω ′ G S ( r , Ω ; r b ′ , Ω ′ ) ρ e f f ( r b ′ ) d b ( r b ′ , Ω ′ ) T ( r b ′ ) I_r(r,\Omega)=\int_{\delta V}dr_b'\int_{n(r_b')\cdot\Omega'<0}d\Omega'G_S(r,\Omega;r_b',\Omega')\rho_{eff}(r_b')d_b(r_b',\Omega')T(r_b') Ir(r,Ω)=δVdrbn(rb)Ω<0dΩGS(r,Ω;rb,Ω)ρeff(rb)db(rb,Ω)T(rb)

Substituing this into Equation (5) we get that
I ( r , Ω ) = I b s ( r , Ω ) + ∫ δ V b d r b ′ ρ e f f ( r b ′ ) T ( r b ′ ) J ( r , Ω ; r b ′ ) I(r,\Omega)=I_{bs}(r,\Omega)+\int_{\delta V_b}dr_b'\rho_{eff}(r_b')T(r_b')J(r,\Omega;r_b') I(r,Ω)=Ibs(r,Ω)+δVbdrbρeff(rb)T(rb)J(r,Ω;rb)
where
J ( r , Ω ; r b ′ ) = ∫ n ( r b ′ ) ⋅ Ω ′ < 0 d Ω ′ G S ( r , Ω ; r ′ , Ω ′ ) d b ( r b ′ , Ω ′ ) J(r,\Omega;r_b')=\int_{n(r_b')\cdot\Omega'<0}d\Omega'G_S(r,\Omega;r',\Omega')d_b(r_b',\Omega') J(r,Ω;rb)=n(rb)Ω<0dΩGS(r,Ω;r,Ω)db(rb,Ω)
By the theory of Green function, J J J is the intensity of radiation field at r with derection Ω \Omega Ω generated by a point anisotropic source d b ( r b ′ , Ω ′ ) δ ( r b , r b ′ ) d_b(r_b',\Omega')\delta(r_b,r_b') db(rb,Ω)δ(rb,rb) located at r b ′ r_b' rb.

With T ( r b ) = ∫ n ( r b ) ⋅ Ω ′ > 0 I ( r b , Ω ′ ) ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ T(r_b)=\int_{n(r_b)\cdot\Omega'>0}I(r_b,\Omega')|\Omega'\cdot n(r_b)|d\Omega' T(rb)=n(rb)Ω>0I(rb,Ω)Ωn(rb)dΩ and (16), (17), we obtain
T ( r b ) = ∫ n ( r b ) ⋅ Ω ′ > 0 ( I b s ( r b , Ω ′ ) + ∫ δ V d r b ′ ρ e f f ( r b ′ ) T ( r b ′ ) J ( r b , Ω ′ ; r b ′ ) ) ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ \begin{aligned} T(r_b)=&\int_{n(r_b)\cdot\Omega'>0}(I_{bs}(r_b,\Omega') \\&+\int_{\delta V}dr_b'\rho_{eff}(r_b')T(r_b')J(r_b,\Omega';r_b'))|\Omega'\cdot n(r_b)|d\Omega' \end{aligned} T(rb)=n(rb)Ω>0(Ibs(rb,Ω)+δVdrbρeff(rb)T(rb)J(rb,Ω;rb))Ωn(rb)dΩ
By letting
T b s = ∫ n ( r b ) ⋅ Ω ′ > 0 I b s ( r b , Ω ) ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ T_{bs}=\int_{n(r_b)\cdot\Omega'>0}I_{bs}(r_b,\Omega)|\Omega'\cdot n(r_b)|d\Omega' Tbs=n(rb)Ω>0Ibs(rb,Ω)Ωn(rb)dΩ
and
G d ( r b , r b ′ ) = ∫ n ( r b ) ⋅ Ω ′ > 0 d Ω ′ J ( r b , Ω ′ ; r b ′ ) ∣ Ω ′ ⋅ n ( r b ′ ) ∣ = ∫ n ( r b ) ⋅ Ω ′ > 0 ∣ Ω ′ ⋅ n ( r b ′ ) ∣ d Ω ′ ∫ n ( r b ′ ) ⋅ Ω ′ < 0 d Ω ′ G S ( r , Ω ; r ′ , Ω ′ ) d b ( r b ′ , Ω ′ ) \begin{aligned} G_d(r_b,r_b')&= \int_{n(r_b)\cdot\Omega'>0}d\Omega'J(r_b,\Omega';r_b')|\Omega'\cdot n(r_b')|\\ &=\int_{n(r_b)\cdot\Omega'>0}|\Omega'\cdot n(r_b')|d\Omega'\int_{n(r_b')\cdot\Omega'<0}d\Omega'G_S(r,\Omega;r',\Omega')d_b(r_b',\Omega') \end{aligned} Gd(rb,rb)=n(rb)Ω>0dΩJ(rb,Ω;rb)Ωn(rb)=n(rb)Ω>0Ωn(rb)dΩn(rb)Ω<0dΩGS(r,Ω;r,Ω)db(rb,Ω)
Noting that G d G_d Gd is downward radiation flux density at r b r_b rb due to point anisotropic source d b ( r b , Ω ′ ) δ ( r b , r b ′ ) d_b(r_b,\Omega')\delta (r_b,r_b') db(rb,Ω)δ(rb,rb).

The (18) becomes
T ( r b ) = T b s + ∫ δ V G d ( r , r ′ ) ρ e f f ( r b ′ ) T ( r b ′ ) d r b ′ T(r_b)=T_{bs} + \int_{\delta V}G_d(r,r')\rho_{eff}(r_b')T(r_b') dr_b' T(rb)=Tbs+δVGd(r,r)ρeff(rb)T(rb)drb
Since the J J J is the radiance generated by d b ( r b ′ , Ω ′ ) δ ( r b , r b ′ ) d_b(r_b',\Omega')\delta(r_b,r_b') db(rb,Ω)δ(rb,rb), so G d G_d Gd​ is actually the flux density generated also by it.

Now, we note that
I S ( r , Ω ) = ∫ δ V d r b ′ J ( r , Ω ; r b ′ ) = ∫ δ V d r b ′ ∫ n ( r b ′ ) ⋅ Ω ′ < 0 d Ω ′ G S ( r , Ω ; r ′ , Ω ′ ) d b ( r b ′ , Ω ′ ) \begin{aligned} I_S(r,\Omega) &=\int_{\delta V}dr_b' J(r,\Omega;r_b')\\ &=\int_{\delta V}dr_b'\int_{n(r_b')\cdot\Omega'<0}d\Omega'G_S(r,\Omega;r',\Omega')d_b(r_b',\Omega') \end{aligned} IS(r,Ω)=δVdrbJ(r,Ω;rb)=δVdrbn(rb)Ω<0dΩGS(r,Ω;r,Ω)db(rb,Ω)
is the intensity of radiation field in canopy generated by anisotropic and heterogeneoussource source d b d_b db located at canopy bottom. It satiesfy the equation
L I S = S I S LI_S=SI_S LIS=SIS
and boundary condtions I S − ( r t , Ω ) = 0 I_S^-(r_t,\Omega)=0 IS(rt,Ω)=0, I S − ( r l , Ω ) = 0 I_S^-(r_l,\Omega)=0 IS(rl,Ω)=0 and I − ( r b , Ω ) = d b ( r b , Ω ) I^-(r_b,\Omega)=d_b(r_b,\Omega) I(rb,Ω)=db(rb,Ω).

This is the second “basic problem”, “S problem”. It follows that
∫ δ V G d ( r , r b ′ ) d r b ′ = ∫ n ( r b ) ⋅ Ω > 0 I S ( r b , Ω ) ∣ Ω ⋅ n ( r b ) ∣ d Ω \int_{\delta_V}G_d(r,r_b')dr_b'=\int_{n(r_b)\cdot\Omega>0}I_S(r_b,\Omega)|\Omega\cdot n(r_b)|d\Omega δVGd(r,rb)drb=n(rb)Ω>0IS(rb,Ω)∣Ωn(rb)dΩ
The integral of G d G_d Gd over boundary is downward flux at r b r_b rb which acount for contribution from all anisotropic sources at canopy bottom.

Rethinking the equation (16), the J is actually a weight to weight the all reflected flux density.

The ground BHR (Equation (11)) and the effective ground anisotropy is independent on vegetation canopy and are known.

Once G d G_d Gd is given, one can solve Equation (18) and substitude it into (15).

Solution of horizontally homogeneous vegetation canopy with Lambertian surface

If we assume that the gournd BRF is homegeneous Lambertian surface, then, T b s T_{bs} Tbs​ and T is independent to (x,y) . And gournd BHR become
ρ e f f = ∫ n ( r b ) ⋅ Ω ′ > 0 ∫ n ( r b ) ⋅ Ω < 0 ∣ n ( r b ) ⋅ Ω ′ ∣ π − 1 ρ b ( r b , Ω ′ , Ω ) ∣ Ω ⋅ n ( r b ) ∣ I ( r b , Ω ′ ) d Ω ′ d Ω ∫ n ( r b ) ⋅ Ω ′ > 0 ∣ Ω ′ ⋅ n ( r b ) ∣ I ( r b , Ω ′ ) d Ω ′ = ρ b ∫ n ( r b ) ⋅ Ω < 0 d Ω ∣ Ω ⋅ n ( r b ) ∣ π = ρ b \begin{aligned} \rho_{eff}&=\frac{\int_{n(r_b)\cdot\Omega'>0}\int_{n(r_b)\cdot\Omega<0}|n(r_b)\cdot\Omega'|\pi^{-1}\rho_b(r_b,\Omega',\Omega)|\Omega\cdot n(r_b)|I(r_b,\Omega')d\Omega'd\Omega}{\int_{n(r_b)\cdot\Omega'>0}|\Omega'\cdot n(r_b)|I(r_b,\Omega')d\Omega'}\\ &=\frac{\rho_b\int_{n(r_b)\cdot\Omega<0}d\Omega|\Omega\cdot n(r_b)|}{\pi}\\ &=\rho_b \end{aligned} ρeff=n(rb)Ω>0Ωn(rb)I(rb,Ω)dΩn(rb)Ω>0n(rb)Ω<0n(rb)Ωπ1ρb(rb,Ω,Ω)∣Ωn(rb)I(rb,Ω)dΩdΩ=πρbn(rb)Ω<0dΩ∣Ωn(rb)=ρb
and effective ground anisotropy is
d b ( r b , Ω ) = R I + T ρ e f f = 1 ρ e f f ⋅ ∫ 2 π + ρ b ( r b , Ω ′ , Ω ) ∣ n ( r b ) ⋅ Ω ′ ∣ I ( r b , Ω ′ ) d Ω ′ π ∫ 2 π − ∣ Ω ′ ⋅ n ( r b ) ∣ I ( r b , Ω ′ ) d Ω ′ = ρ b ρ e f f π = 1 π \begin{aligned} d_b(r_b,\Omega)&=\frac{\mathcal{R}I^+}{T\rho_{eff}}\\ &=\frac{1}{\rho_{eff}}\cdot\frac{\int_{2\pi+}\rho_b(r_b,\Omega',\Omega)|n(r_b)\cdot\Omega'|I(r_b,\Omega')d\Omega'}{\pi\int_{2\pi-}|\Omega'\cdot n(r_b)|I(r_b,\Omega')d\Omega'}\\ &=\frac{\rho_b}{\rho_{eff}\pi}=\frac{1}{\pi} \end{aligned} db(rb,Ω)=TρeffRI+=ρeff1π2πΩn(rb)I(rb,Ω)dΩ2π+ρb(rb,Ω,Ω)n(rb)ΩI(rb,Ω)dΩ=ρeffπρb=π1
Then the Equation (21) become
T = T b s + ∫ δ V G d ( r b , r b ′ ) ρ e f f T d r b ′ = T b s + ρ b T ∫ δ V G d ( r b , r b ′ ) d r b ′ → T ( 1 − ρ b T R S ) = T b s → T = T b s 1 − ρ b R s \begin{aligned} T &=T_{bs} + \int_{\delta V}G_d(r_b,r_b')\rho_{eff}T dr_b'\\ &= T_{bs} + \rho_bT\int_{\delta V}G_d(r_b,r_b')dr_b'\\ &\rightarrow T(1-\rho_bTR_S)=T_{bs}\\ &\rightarrow T=\frac{T_{bs}}{1-\rho_bR_s} \end{aligned} T=Tbs+δVGd(rb,rb)ρeffTdrb=Tbs+ρbTδVGd(rb,rb)drbT(1ρbTRS)=TbsT=1ρbRsTbs
where
R S = ∫ δ V G d ( r b , r b ′ ) d r b ′ R_S=\int_{\delta V}G_d(r_b,r_b')dr_b' RS=δVGd(rb,rb)drb
is the downward flux density at r b r_b rb generated by source 1 π \frac{1}{\pi} π1 (since d b d_b db has become 1 / π 1/\pi 1/π and combine the Equation (24) and (23)).

Substituting Eq. (22) into Eq. (16), we obtain:
I ( r , Ω ) = I b s ( r , Ω ) + I S ( r , Ω ) ρ e f f T I(r,\Omega)=I_{bs}(r,\Omega)+I_S(r,\Omega)\rho_{eff}T I(r,Ω)=Ibs(r,Ω)+IS(r,Ω)ρeffT
Using ρ e f f = ρ b \rho_{eff}=\rho_b ρeff=ρb, Then, by Eq.(27), we obtain:
I ( r , Ω ) = I b s ( r , Ω ) + T b s ρ b 1 − ρ b R S I S ( r , Ω ) I(r,\Omega)=I_{bs}(r,\Omega)+\frac{T_{bs}\rho_{b}}{1-\rho_bR_S}I_S(r,\Omega) I(r,Ω)=Ibs(r,Ω)+1ρbRSTbsρbIS(r,Ω)
So, the solution of original boundary problem is decomposed into the solution of black-soil problem and S problem.

Solution of horizontally inhomogeneous vegetation canopy

In this case, the solution of Eq. (21) is not directly Eq. (27), and approximation is performed. Consider the ratio
R S ( r b ) = ∫ δ V G d ( r b , r b ′ ) ρ e f f ( r b ′ ) T ( r b ′ ) d r b ′ ρ e f f ( r b ′ ) T ( r b ′ ) R_S(r_b)=\frac{\int_{\delta V}G_d(r_b,r_b')\rho_{eff}(r_b')T(r_b')dr_b'}{\rho_{eff}(r_b')T(r_b')} RS(rb)=ρeff(rb)T(rb)δVGd(rb,rb)ρeff(rb)T(rb)drb
which is BHR calculated for canopy illuminated by anisotropic sources ρ e f f ( r b ′ ) T ( r b ′ ) d b ( r b ′ , Ω ′ ) \rho_{eff}(r_b')T(r_b')d_b(r_b',\Omega') ρeff(rb)T(rb)db(rb,Ω) from below. To understand this, we could think that the numerator is the responsed flux density at r b r_b rb due to radiance ρ e f f ( r b ′ ) T ( r b ′ ) d b ( r b ′ , Ω ′ ) \rho_{eff}(r_b')T(r_b')d_b(r_b',\Omega') ρeff(rb)T(rb)db(rb,Ω) by Eq. (24), and the denomator is
∫ n ( r b ′ ) ⋅ Ω ′ < 0 ρ e f f ( r b ′ ) T ( r b ′ ) d ( r b ′ Ω ′ ) ∣ n ( r b ′ ) ⋅ Ω ′ ∣ d Ω ′ \int_{n(r_b')\cdot\Omega'<0}\rho_{eff}(r_b')T(r_b')d(r_b'\Omega')|n(r_b')\cdot\Omega'|d\Omega' n(rb)Ω<0ρeff(rb)T(rb)d(rbΩ)n(rb)ΩdΩ
which is equal to ρ e f f ( r b ′ ) T ( r b ′ ) \rho_{eff}(r_b')T(r_b') ρeff(rb)T(rb).

Since the T varies, so the R s R_s Rs varies. But linear operator theory has the conclusion that a continous positive linear operator B B B minimum m n m_n mn and maximum M n M_n Mn values of the function η n = ( B n u ) 1 / n \eta_n=(\sqrt{B^nu})^{1/n} ηn=(Bnu )1/n converge to mazimum eigenvalue ρ ( B ) \rho(B) ρ(B) for any chosen positive function u u u as n tend to infinity.

Here, the numerator in Eq. (35) can be treated as a positive intergral operator B B B with a kernel G d G_d Gd. We have n = 1 n=1 n=1, and we use the maximum eigenvalue R S ˉ \bar{R_S} RSˉ of B B B with kernel G d G_d Gd.

Then , we rewrite the Eq. (25) as
T ( r b ) = T b s ( r b ) + R S ˉ ρ e f f ( r b ) T ( r b ) . T(r_b)=T_{bs}(r_b)+\bar{R_S}\rho_{eff}(r_b)T(r_b). T(rb)=Tbs(rb)+RSˉρeff(rb)T(rb).
Then we have
T ( r b ) = T b s ( r b ) 1 − R S ˉ ρ e f f ( r b ) T(r_b)=\frac{T_{bs}(r_b)}{1-\bar{R_S}\rho_{eff}(r_b)} T(rb)=1RSˉρeff(rb)Tbs(rb)
Substitude this into Eq. (16), we obtain
I ( r , Ω ) = I b s ( r , Ω ) + ∫ δ V d r b ′ J ( r b ′ ; r b ′ , Ω ′ ) ρ e f f ( r b ′ ) T b s ( r b ′ ) 1 − R S ˉ ρ e f f ( r b ′ ) I(r,\Omega)=I_{bs}(r,\Omega)+\int_{\delta_V}dr_b'J(r_b';r_b',\Omega')\rho_{eff}(r_b')\frac{T_{bs}(r_b')}{1-\bar{R_S}\rho_{eff}(r_b')} I(r,Ω)=Ibs(r,Ω)+δVdrbJ(rb;rb,Ω)ρeff(rb)1RSˉρeff(rb)Tbs(rb)
where J has already mentioned below as
J ( r , Ω ; r b ′ ) = ∫ n ( r b ′ ) ⋅ Ω ′ < 0 d Ω ′ G S ( r , Ω ; r ′ , Ω ′ ) d b ( r b ′ , Ω ′ ) J(r,\Omega;r_b')=\int_{n(r_b')\cdot\Omega'<0}d\Omega'G_S(r,\Omega;r',\Omega')d_b(r_b',\Omega') J(r,Ω;rb)=n(rb)Ω<0dΩGS(r,Ω;r,Ω)db(rb,Ω)
For futher simlification, we replace the J J J​ with its mean value over the canopy bottom, then this becomes
I ( r , Ω ) = I b s ( r , Ω ) + ∫ δ V d r b ′ J ( r b ′ ; r b ′ , Ω ′ ) ρ e f f ˉ T b s ˉ 1 − R S ˉ ρ e f f ˉ I(r,\Omega)=I_{bs}(r,\Omega)+\int_{\delta_V}dr_b'J(r_b';r_b',\Omega')\frac{\bar{\rho_{eff}}\bar{T_{bs}}}{1-\bar{R_S}\bar{\rho_{eff}}} I(r,Ω)=Ibs(r,Ω)+δVdrbJ(rb;rb,Ω)1RSˉρeffˉρeffˉTbsˉ
Then the integral of J J J is I S I_S IS​ by Eq. (22), then we have
I ( r , Ω ) = I b s ( r , Ω ) + ρ e f f ˉ T b s ˉ 1 − R S ˉ ρ e f f ˉ I S ( r , Ω ) I(r,\Omega)=I_{bs}(r,\Omega)+\frac{\bar{\rho_{eff}}\bar{T_{bs}}}{1-\bar{R_S}\bar{\rho_{eff}}}I_S(r,\Omega) I(r,Ω)=Ibs(r,Ω)+1RSˉρeffˉρeffˉTbsˉIS(r,Ω)
Then, the results of original boundary problem become the solution of black soil problem and S problem.

Now, we could rewrite the HDRF of the surface :
H D R F = < I ( Ω ) > 0 1 π ∫ 2 π − < I ( Ω ′ ) > 0 ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ = < I b s ( r , Ω ) + ρ e f f ˉ T b s ˉ 1 − R S ˉ ρ e f f ˉ I S ( r , Ω ) > 1 π ∫ 2 π − < I ( Ω ′ ) > 0 ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ = R b s ( Ω ) + ρ e f f ˉ 1 − R S ˉ ρ e f f ˉ t b s T S ( Ω ) \begin{aligned} HDRF&=\frac{<I(\Omega)>_0}{\frac{1}{\pi}\int_{2\pi-}<I(\Omega')>_0 |\Omega'\cdot n(r_b)|d\Omega'}\\ &=\frac{<I_{bs}(r,\Omega)+\frac{\bar{\rho_{eff}}\bar{T_{bs}}}{1-\bar{R_S}\bar{\rho_{eff}}}I_S(r,\Omega)>}{\frac{1}{\pi}\int_{2\pi-}<I(\Omega')>_0 |\Omega'\cdot n(r_b)|d\Omega'}\\ &=R_{bs}(\Omega)+\frac{\bar{\rho_{eff}}}{1-\bar{R_S}\bar{\rho_{eff}}}t_{bs}T_S(\Omega) \end{aligned} HDRF=π12π<I(Ω)>0Ωn(rb)dΩ<I(Ω)>0=π12π<I(Ω)>0Ωn(rb)dΩ<Ibs(r,Ω)+1RSˉρeffˉρeffˉTbsˉIS(r,Ω)>=Rbs(Ω)+1RSˉρeffˉρeffˉtbsTS(Ω)
Here, transmitance is
t b s = T b s ˉ ∫ 2 π − < I ( Ω ′ ) > 0 ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ t_{bs}=\frac{\bar{T_{bs}}}{\int_{2\pi-}<I(\Omega')>_0 |\Omega'\cdot n(r_b)|d\Omega'} tbs=2π<I(Ω)>0Ωn(rb)dΩTbsˉ
and
T s = < I S ( Ω ) > 0 1 π ∫ 2 π − < d b ( Ω ′ ) > 0 ∣ Ω ′ ⋅ n ( r b ) ∣ d Ω ′ T_s=\frac{<I_S(\Omega)>_0}{\frac{1}{\pi}\int_{2\pi-}<d_b(\Omega')>_0 |\Omega'\cdot n(r_b)|d\Omega'} Ts=π12π<db(Ω)>0Ωn(rb)dΩ<IS(Ω)>0
which is the radiance generated by the isotropic source located at canopy bottom, noting that
∫ 2 π + d b ( r b , Ω ) ∣ n ( r b ) ⋅ Ω ∣ d Ω = 1. \int_{2\pi+}d_b(r_b,\Omega)|n(r_b)\cdot\Omega|d\Omega=1. 2π+db(rb,Ω)n(rb)Ω∣dΩ=1.
So its actually
T s = < I S ( Ω ) > 0 π T_s=<I_S(\Omega)>_0\pi Ts=<IS(Ω)>0π

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2234547.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

智会智展,活动必备

智会智展 APP 各大应用市场均可下载统一链接https://m.malink.cn/s/r6nQVf

Hive操作库、操作表及数据仓库的简单介绍

数据仓库和数据库 数据库和数仓区别 数据库与数据仓库的区别实际讲的是OLTP与OLAP的区别 操作型处理(数据库)&#xff0c;叫联机事务处理OLTP&#xff08;On-Line Transaction Processing&#xff09;&#xff0c;也可以称面向用户交易的处理系统&#xff0c;它是针对具体业务…

如何选择适合小团队的项目管理工具?免费与开源软件推荐

目录 一、小团队项目管理工具的重要性 二、热门项目管理工具介绍 &#xff08;一&#xff09;禅道 &#xff08;二&#xff09;Trello &#xff08;三&#xff09;Asana &#xff08;四&#xff09;JIRA 三、免费项目管理软件推荐 &#xff08;一&#xff09;ES 管理器 …

kafka如何获取 topic 主题的列表?

大家好&#xff0c;我是锋哥。今天分享关于【kafka如何获取 topic 主题的列表&#xff1f;】面试题&#xff1f;希望对大家有帮助&#xff1b; kafka如何获取 topic 主题的列表&#xff1f; 1000道 互联网大厂Java工程师 精选面试题-Java资源分享网 在Kafka中&#xff0c;可以…

Maven详解—(详解Maven,包括Maven依赖管理以及声明周期,Maven仓库、idea集成Maven)

文章目录 Maven详解一.初始Maven1.1 概述1.2 作用 二.Maven模型2.1 概述2.2 构建生命周期/阶段2.3 项目对象模型2.4 依赖管理模型 三.Maven仓库四.Maven安装4.1 下载4.2 安装步骤 五.Idea集成Maven Maven详解 一.初始Maven 1.1 概述 Maven是Apache旗下的一个开源项目&#x…

虚拟滚动 - 从基本实现到 Angular CDK

简介 在大数据列表的处理上&#xff0c;虚拟滚动是一种优化性能的有效方式。本篇文章将详细介绍两种常见的虚拟滚动实现方式&#xff1a;使用 transform 属性和 Intersection Observer。重点讲解如何通过 transform 属性实现高效的虚拟滚动&#xff0c;并对比Angular CDK中的实…

Spring Boot 配置文件启动加载顺序

前言 Spring Boot的启动加载顺序是一个涉及多个步骤和组件的过程。Spring Boot通过一系列默认设置简化了应用程序的配置&#xff0c;使得开发者能够快速地搭建和部署应用。为了实现这一目标&#xff0c;Spring Boot采用了一种分层和优先级机制来加载配置文件。 一、Spring Bo…

C# Modbus RTU通讯回顾

涉及技术&#xff1a; 1.使用NMdbus4 库 2.ushort[]转int 记得之前刚学习的时候&#xff0c;是ushort[] → Hex字符串→byte[] → 翻转byte[] →BitConverter.ToInt32()&#xff0c;饶了一大圈&#xff1b;实际上可以直接转&#xff1b;这里也有小细节&#xff1a;使用BitCo…

HFSS学习笔记(五)金属过孔、复制模型带激励等问题(持续更新...)

HFSS学习笔记&#xff08;五&#xff09;金属过孔、复制模型带激励等问题&#xff08;持续更新…&#xff09; 一、金属过孔设计 方法一&#xff1a;用介质减去金属圆柱体&#xff0c;然后再添加金属圆柱体 方法二&#xff1a;嵌入金属圆柱 圆柱过孔选择材料为“copper” HFS…

Late Chunking×Milvus:如何提高RAG准确率

01. 背景 在RAG应用开发中&#xff0c;第一步就是对于文档进行chunking&#xff08;分块&#xff09;&#xff0c;高效的文档分块&#xff0c;可以有效的提高后续的召回内容的准确性。而对于如何高效的分块是个讨论的热点&#xff0c;有诸如固定大小分块&#xff0c;随机大小分…

大屏可视化:舞动数据与美观的“设计秘籍”

大屏可视化鉴赏&#xff1a;踏入软件系统产品设计之旅&#xff0c;让我们一同鉴赏那些闪耀在智慧农业、智慧园区、智慧社区及智慧港口等领域的大屏可视化杰作。每一帧画面&#xff0c;都是科技与创新的完美融合&#xff0c;数据跃然屏上&#xff0c;智慧触手可及。 >> 数…

基于STM32的智能声音跟随小车设计

引言 本项目基于STM32微控制器设计了一个智能声音跟随小车&#xff0c;通过集成麦克风阵列实现声音源定位和跟随功能。该系统可以检测环境中的声音信号&#xff0c;如手掌拍击声或语音指令&#xff0c;驱动小车向声源方向移动。项目涉及硬件设计、声音信号处理算法以及电机控制…

Bruno解决SSL验证问题

在测试接口的时候&#xff0c;我使用的是Bruno这个软件&#xff0c;开源离线的API测试软件。 主页是这样子的 今天在测试一个HTTPS的接口时候&#xff0c;因为这个HTTPS接口是用的是自签证书&#xff0c;所以就报错误了。 Error invoking remote method send-http-request: …

【论文速读】| APOLLO:一种基于 GPT 的用于检测钓鱼邮件并生成警告用户的解释的工具

基本信息 原文标题&#xff1a;APOLLO: A GPT-based tool to detect phishing emails and generate explanations that warn users 原文作者&#xff1a;Giuseppe Desolda, Francesco Greco, Luca Vigan 作者单位&#xff1a;University of Bari “A. Moro”, Italy, King’…

jfrog artifactory oss社区版,不支持php composer私库

一、docker安装 安装环境&#xff1a;centos操作系统&#xff0c;root用户。 如果是mac或ubuntu等操作系统的话&#xff0c;会有许多安装的坑等着你。 一切都是徒劳&#xff0c;安装折腾那么久&#xff0c;最后还是不能使用。这就是写本文的初衷&#xff0c;切勿入坑就对了。 …

WindowsDocker安装到D盘,C盘太占用空间了。

Windows安装 Docker Desktop的时候,默认位置是安装在C盘,使用Docker下载的镜像文件也是保存在C盘,如果对Docker使用评率比较高的小伙伴,可能C盘空间,会被耗尽,有没有一种办法可以将Docker安装到其它磁盘,同时Docker的数据文件也保存在其他磁盘呢? 答案是有的,我们可以…

vue常见题型(1-10)

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 文章目录 2.2双向绑定的原理是什么vue框架采用的是数据双向绑定的方式&#xff0c;由三个重要部分构成2.2.1.ViewModel2.2.2 双向绑定2.2.3.1.编译Compile2.2.3.2.依赖收集 3…

python怎么将字符串转换为数字

python如何将列表中的字符串转为数字&#xff1f;具体方法如下&#xff1a; 有一个数字字符的列表&#xff1a; numbers [1, 5, 10, 8] 想要把每个元素转换为数字&#xff1a; numbers [1, 5, 10, 8] 用一个循环来解决&#xff1a; new_numbers []; for n in numbers:new_n…

大数据新视界 -- 大数据大厂之 Impala 性能优化:解锁大数据分析的速度密码(上)(1/30)

&#x1f496;&#x1f496;&#x1f496;亲爱的朋友们&#xff0c;热烈欢迎你们来到 青云交的博客&#xff01;能与你们在此邂逅&#xff0c;我满心欢喜&#xff0c;深感无比荣幸。在这个瞬息万变的时代&#xff0c;我们每个人都在苦苦追寻一处能让心灵安然栖息的港湾。而 我的…

Flutter 插件 sliding_up_panel 实现从底部滑出的面板

前言 sliding_up_panel 是一个 Flutter 插件&#xff0c;用于实现从底部滑出的面板。它在设计上非常灵活&#xff0c;能够适应多种 UI 场景&#xff0c;比如从底部滑出的菜单、可拖动的弹出面板等。以下是 sliding_up_panel 的详细用法&#xff0c;包括常用的参数说明和示例代…