读书笔记【顾行发，航天光学遥感器辐射定标原理与方法】

P4：像元级辐射定标可狭义上理解为图像的调制传递函数补偿MTFC (Modulation Transfer Function Compensation)

P9：敦煌辐射校正场→VIS+NIR；青海湖辐射校正场→TIR

P21：TERRA上午10:30过境，被称为地球观测第一颗上午星 (EOS-AM1)；AQUA下午01:30过境，被称为地球观测下午第一颗星 (EOS-PM1)

P22：为了监视太阳反射波段辐射定标的稳定度，MODIS定标小组还定期地进行遥感器对月球的测量。当月球相位角约为55.5°时，通过发射命令要求整个飞行器进行姿态变化，使得MODIS能通过对地观测区域或空间观测区域对月球进行测量，测量的结果直接用于遥感器的在轨辐射定标及监视遥感器的长期辐射响应的稳定性。

P23：MODIS的1，2，4波段的场地在规定标达到预期结果，两个遥感器都可用于监测在轨的其他性能。

P25：SPOT-1首次运用了线性阵列传感器和推扫技术。

P60：Ground Resolution: $GR=2\cdot H\cdot \tan \frac{IFOV}{2}$

P173：VIS+NIR的内置定标灯采用小型钨丝灯作为灯源，为了使光源稳定，可利用带滤光片的硅探测器监视光源变化，通过反馈电路控制灯的电流以维持灯的辐射稳定。
“灯+积分球”星上定标可以实现遥感器的全孔径，全视场定标，而“灯+漫射板”难以实现。
标准灯与太阳的光谱差异较大。

P175：“太阳+漫射板”→漫射板在外太空，经太阳直射后反射率严重衰减；“太阳+衰减板+漫射板”→衰减板位于漫射板前方，但衰减板也会由于外太空环境而引起反射率变化。

P178：到目前为止，月球辐射的光谱特性仍不清楚。

P179：由于月球盈缺，平衡摆动和日月距离变化引起的月球表面反射率的变化可以忽略，使得其有效反射率可认为基本恒定。

P186：MODIS-AQUA的性能较优于MODIS-TERRA。

P193：通过分析影响月球表观辐亮度变化的五个因子，可以得到星上定标时刻的月球表观辐亮度。

1. 第一个归一化因子为太阳的月球的距离校正因子
   $k_1=(\frac{D_{SM}}{U}{)}^2$, $D_{SM}=D_{SE}+R$
   $D_{SM}$：太阳和月球的距离
   $U$：日地平均距离
   $D_{SE}$：日地距离
   $R$：月球轨道半径
2. 第二个归一化因子是仪器和月球的距离校正
   $k_2=(\frac{D_{IM}}{R}{)}^2$, $D_{IM}=D_{EM}-A-H$
   $D_{IM}$：仪器到月球的距离
   $A$：地球半径
   $H$：仪器离地高度
3. 第三个正常化因子是月球表面的发光部分，可作为一个相角函数
   $k_3=\frac{0.9611}{f_1(\theta )}$, $f_1(\theta )=a_0+a_1\theta$
   $a_0=1$, $a_1=\frac{-1}{180{\mathrm{deg}}^{-1}}$
   $\theta$：0°相位角时月球表面被完全照亮，90°相位角时月球表面亮起一半，180°时月球表面是黑暗的。
4. 第四个归一化因子随月球相位角变化而改变其亮度，是一个月球相位角与反射率变化的函数
   $k_4=\frac{f_2(7)}{f_2(\theta )}$, $f_2(\theta )=b_0+b_1\theta +b_2{\theta }^2$
   $b_0=0.1287$, $b_1=-6.702\times 10^{-3}{\mathrm{deg}}^{-1}$
   $b_2=2.163\times 10^{-4}{\mathrm{deg}}^{-2}$, $f_2(7)=0.09238$
5. 第五个归一化因子是月球观测期间的仪器俯仰速率
   $k_5=\frac{25}{L_M\cdot \frac{D_{IM}}{R}}=(\frac{25}{L_M})\cdot (\frac{R}{D_{IM}})$
   $L_M$：图像最长垂直部分1％相应点的间隔

P199：场地定标的方法有反射率基法，辐照度基法，辐亮度基法。前两者为地基准同步，需要高朗伯体，后者为空基准同步，也可以定标TIR，更适用于低空间分辨率。

P290：撒哈拉沙漠和沙特阿拉伯沙漠作为稳定场景用于监测遥感器的衰减变化。

P291：稳定场景定标分为沙漠场景(气象卫星等低空间分辨率的时序定标)，极地场景，海洋场景(海洋，气象卫星)和云场景。

• 沙漠场景法时间序列定标的两个步骤：①根据图像的成像日期，观测几何和场地DN值，建立不同时间遥感器衰减关系，得到遥感器的时间序列相对定标系数。②以实地同步测量数据为标准，计算出当天的绝对辐射定标系数，根据不同时间定标系数相对的衰减关系，转换得到卫星运行期间时间序列的绝对辐射定标系数。
• 极地场景法：最理想的定标场是南极的 Dome Concordia(DomeC)站点，因为①DomeC比南极中心点更容易获取影像；②DomeC有实验站，可获取当地气象资料。注意①因极昼极夜现象，南极只能在11月-2月进行，格陵兰为5月-8月；②极地太阳天顶角大，BRDF影响严重；③极地臭氧含量变化显著，可以用臭氧产品去除影响；④云的影响，可以采用红外通道区分云和雪。
• 海洋场景法：主要对蓝绿波段进行绝对辐射定标。
• 云场景法：主要实现两个通道的相对辐射定标。

P313：一个成像系列可有效地看作一个空间频率滤波器，而它的成像特征和像质评价，能以物象之间的频率之比来表示，这个频率对比特性就是光学传递函数OTF(中国大百科全书定义的OTF：以空间频率为变量，表征成像过程中调制度与横向相移的相对变化的函数)。光学传递函数的幅值就是调制传递函数。
调制度：$M=\frac{L_{max}-L_{min}}{L_{max}+L_{min}}$
MTF：$MTF=\frac{M_i}{M_0}$
$M_i$：经过系统的图像调制度
$M_0$：未经系统的目标调制度
MTF的值越大，表明M_像与M_物越接近，光能的分配改变程度越小，光学系统的成像性能越好。

P316：线阵推扫CCD相机整个成像系统的MTF是静态MTF和动态MTF之积：
$MT{F}_{system}=MT{F}_{optics}\cdot MT{F}_{detector}\cdot MT{F}_{motion}\cdot MT{F}_{atmosphere}$
其system为成像系统，optics为光学系统，detector为探元，motion为像移，atmosphere为大气。大气MTF为大气气溶胶与湍流MTF和综合。

P361：干涉光谱成像仪(HSI)的辐射定标需要标定各个像元的光谱辐射响应函数，而色散型光谱成像仪只需要标定一行的像元。

P368：HJ1A-HSI的前20个通道有严重的纵向条纹。

1. 加性噪声去除
   $D{N}_{new}=D{N}_{org}+P_K-p_{iK}$
   $P_K$：第K波段均值
   $p_{i}K$：第K波段第i列均值
2. 乘性噪声去除
   $D{N}_{new}=D{N}_{org}\cdot \frac{p_{iK}}{P_K}$
3. 单波段图像匹配法
   $N{G}_i=\frac{{\sigma }_{i}K}{{\sigma }_K}$, ${\sigma }_K=\frac{1}{N}\sum _{j=1}^N{\sigma }_j$
   $B_i=p_{iK}-\frac{{\sigma }_{iK}\cdot P_K}{{\sigma }_K}$
   $$\begin{gathered} D{N}_{new}=\frac{D{N}_{org}-(p_{iK}-\frac{{\sigma }_{iK}\cdot P_K}{{\sigma }_K})}{\frac{{\sigma }_{iK}}{{\sigma }_K}} \\ =\frac{{\sigma }_K}{{\sigma }_{iK}}\cdot D{N}_{org}-\frac{{\sigma }_K}{{\sigma }_{iK}}\cdot p_{i}K+P_K \end{gathered}$$
4. 多波段图像匹配法
   $P_{iK}^{\prime }=P_K\cdot \frac{p_{iR}}{P_R}$
   ${\sigma }_{iK}^{\prime }={\sigma }_K\cdot \frac{{\sigma }_{iR}}{{\sigma }_R}$
   其R表示匹配的第R波段，并将新的$P_{iK}^{\prime }$和${\sigma }_{iK}^{\prime }$带入单波段匹配法中
5. 多波段窗口匹配法(最为理想)
   $D{N}_{new}(i,j,k)=D{N}_{org}(i,j,k)\cdot \frac{\mu (K)}{\mu (R)}$

P391：国外高光谱往往采用高斯型函数构造光谱响应函数。

函数名称	函数公式	宽度	备注
脉冲函数	$\{\begin{array}{ll}1& \lambda ={\lambda }_i\\ 0& \lambda \ne {\lambda }_i\end{array}$	独立点	${\lambda }_i$为中心波长
矩形函数	$\{\begin{array}{ll}1& \lambda \in ({\lambda }_1,{\lambda }_2)\\ 0& \lambda \notin ({\lambda }_1,{\lambda }_2)\end{array}$	$1\times FWHM$	${\lambda }_1$和${\lambda }_2$为起始和终止波长
三角型函数	$\{\begin{array}{ll}1-|\frac{2\lambda -{\lambda }_2-{\lambda }_1}{{\lambda }_2-{\lambda }_1}|& \lambda \in ({\lambda }_1,{\lambda }_2)\\ 0& \lambda \notin ({\lambda }_1,{\lambda }_2)\end{array}$	$2\times FWHM$	${\lambda }_1$和${\lambda }_2$为起始和终止波长
高斯函数	$\frac{1}{\sigma \lambda \sqrt{2\pi }}\cdot e^{-\frac{(\lambda -{\lambda }_i{)}^2}{2{\sigma }^2}}$	$6\times FWHM$	$\sigma =\frac{FWHM}{2\sqrt{2\ln 2}}$
$\sin c$函数	$\sin c=\frac{\sin x}{x}$	$2\times FWHM$	$x=\frac{\lambda -{\lambda }_i}{FWHM}$, ${\lambda }_i$为中心波长
$\sin c^2$函数	$\sin c^2=(\frac{\sin x}{x}{)}^2$	$2\times FWHM$	$x=\frac{\lambda -{\lambda }_i}{FWHM}$, ${\lambda }_i$为中心波长

P394：光谱响应函数对定标结果的影响主要和大气吸收密切相关，如果通道位于大气窗口，则不论采用何种光谱响应函数，得到的定标结果都具有一致性；如果通道位于大气吸收波段，光谱响应函数的差异对定标结果有显著的影响。

P474：偏振遥感不仅利用了电磁波具有能量大小的特点，还利用偏振信息来反演地球物化参数。
定量遥感的第一步就是遥感器的定标。

P483：通常认为偏振光主要是由气溶胶的散射特性造成的。

P484：植被，土壤和太阳耀光是偏振遥感研究中最广泛的三种。
BPDF：二向偏振分布函数

P492：偏振定标必须建立观测值与矢量化的偏振态之间的计算关系。

自己的感悟

• 学习不只是实践自学和上网自学，书本上的理论也应该在新媒体时代引起重视。
• 改变了对于MODIS的看法，其在定标上确实很优秀，同时对我国的航天遥感和对地观测寄予厚望，不要只顾着把卫星送上天就草草了事了。
• 码了一晚上博客，回顾了笔记本上的内容，更多地，学会了好多$\text{KATEX}$的书写格式，公式渲染真的很强大。