【遥感学习整理 04】生态环境遥感监测:从植被覆盖度到城市绿地与气溶胶反演

一次实时数据链路丢数排查:从 MQTT 队列积压到 SQL 统计开销 本文通过一次实时数据链路压测,排查了从MQTT到ClickHouse的端到端数据丢失问题。核心方法是为每条消息添加序号,建立独立观察点(MQTT订阅端、Kafka消费者),结合转发服务阶段统计与线程栈分析,定位到下游处理速度不足源于Druid SQL合并统计的高开销。优化措施包括关闭SQL合并、扩容代理队列、增强生产侧失败处理,并验证15分钟内实现全链路对账一致,显著缓解积压与丢包。关键收获:数据可对账、明确各阶段“成功”定义、关注客户端队列状态。 阅读详情

生态环境遥感监测:从植被覆盖度到城市绿地与气溶胶反演

文章目录

第一章 生态环境遥感监测的统一框架

1.1 生态环境遥感监测解决什么问题

生态环境遥感监测的核心,不是单纯生成一幅分类图,而是回答以下问题:

  • 地表生态状态如何;
  • 植被覆盖是否增加或退化;
  • 草原和沙化区域如何变化;
  • 城市绿地分布是否合理;
  • 大气气溶胶浓度在空间上如何分布;
  • 多种生态因子如何综合形成区域评价结果。

不同专题看似差异明显,但都可以统一为“对象—信号—参数—评价—验证”的分析链条:

生态对象

遥感观测信号

预处理

指数、分类或反演

生态参数

时空分析

综合评价

地图、统计与管理成果

其中:

  • 对象包括植被、裸土、城市绿地、气溶胶和地形;
  • 信号包括反射率、热辐射、角度信息和纹理;
  • 参数包括植被覆盖度、气溶胶光学厚度、坡度和土壤指数;
  • 评价包括变化检测、等级划分和综合指数计算。

1.2 五类专题之间的关系

本篇整理的五个应用专题分别对应生态监测中的不同层次:

专题主要对象核心参数主要方法
自然生态环境监测区域生态系统植被、土壤、坡度综合指数多因子综合评价
城市绿地提取城市植被对象绿地范围和面积高分影像、面向对象提取
植被覆盖度反演像元内植被比例VFC/FVCNDVI像元二分模型
草原与沙漠化监测草地植被变化植被覆盖度与变化区多时相变化检测
MODIS气溶胶监测大气颗粒物AOD暗像元法与查找表

1.3 生态监测不能只依赖单景影像

生态过程具有明显的时间性。单景影像往往只能反映某一日期的状态,而无法区分:

  • 短期天气扰动;
  • 季节物候变化;
  • 长期退化;
  • 灾害影响;
  • 管理措施效果;
  • 传感器或大气差异。

因此,生态环境监测通常需要同时重视:

  1. 多时相一致性;
  2. 数据预处理质量;
  3. 参数物理意义;
  4. 阈值的区域适用性;
  5. 外部验证与误差分析。

第二章 生态环境监测的数据预处理基础

2.1 为什么预处理决定生态参数可信度

植被覆盖度、气溶胶光学厚度和综合生态指数都依赖像元值。如果输入影像受到大气、几何错位、云影或缩放因子影响,最终生态参数会产生系统偏差。

预处理通常包括:

原始遥感数据

数据读取

辐射定标

几何校正或配准

大气校正

镶嵌与裁剪

云、阴影和异常值处理

统一空间分辨率和投影

2.2 辐射定标

原始DN需要转换为具有物理意义的辐亮度或反射率。

线性定标可表示为:

L λ = G λ × D N + B λ L_{\lambda} = G_{\lambda}\times DN + B_{\lambda} Lλ​=Gλ​×DN+Bλ​

其中:

  • L λ L_{\lambda} Lλ​为谱辐亮度;
  • G λ G_{\lambda} Gλ​为增益;
  • B λ B_{\lambda} Bλ​为偏移量;
  • D N DN DN为原始数字量化值。

若进行植被指数、气溶胶反演或多时相比较,应优先使用表观反射率或地表反射率,而不是直接使用DN。

2.3 大气校正

生态监测中,大气校正的作用是削弱散射、吸收和路径辐射影响。

传感器观测辐亮度可以概念化表示为:

L λ s e n s o r = L λ p a t h + T λ u p ⋅ L λ s u r f a c e L_{\lambda}^{sensor} = L_{\lambda}^{path} + T_{\lambda}^{up} \cdot L_{\lambda}^{surface} Lλsensor​=Lλpath​+Tλup​⋅Lλsurface​

大气校正尤其影响:

  • 红光和近红外波段;
  • NDVI;
  • 植被覆盖度;
  • 高分辨率绿地提取;
  • MODIS气溶胶反演;
  • 多时相变化检测。

2.4 图像配准

多时相草原变化监测要求不同年份影像严格对齐。否则,边缘错位会形成伪变化。

配准误差可通过均方根误差评价:

R M S E = 1 n ∑ i = 1 n [ ( Δ x i ) 2 + ( Δ y i ) 2 ] RMSE = \sqrt{ \frac{1}{n} \sum_{i=1}^{n} \left[ (\Delta x_i)^2 + (\Delta y_i)^2 \right] } RMSE=n1​i=1∑n​[(Δxi​)2+(Δyi​)2] ​

变化检测中,RMSE不仅要“数值小”,还要检查道路、河流、地块和边界是否真正重合。

2.5 镶嵌与裁剪

研究区超过单景影像覆盖范围时,需要镶嵌。镶嵌前应统一:

  • 投影;
    -像元大小;
  • 波段顺序;
  • 辐射尺度;
  • NoData;
  • 成像季节;
  • 大气处理级别。

裁剪的作用不仅是减少数据量,还能:

  • 限定统计范围;
  • 避免研究区外极端值影响;
  • 提高运算效率;
  • 便于按行政区或生态区汇总。

第三章 植被覆盖度的概念与意义

3.1 植被覆盖度是什么

植被覆盖度是植被叶、茎、枝在地面的垂直投影面积占统计区总面积的比例。

设像元面积为 A A A,其中植被投影面积为 A v A_v Av​,则植被覆盖度为:

V F C = A v A VFC = \frac{A_v}{A} VFC=AAv​​

其值通常位于 [ 0 , 1 ] [0,1] [0,1]之间,也可表示为百分比。

3.2 植被覆盖度与NDVI的区别

NDVI是光谱指数,植被覆盖度是地表比例参数。二者相关,但并不等价。

NDVI计算为:

N D V I = N I R − R e d N I R + R e d NDVI = \frac{NIR-Red} {NIR+Red} NDVI=NIR+RedNIR−Red​

NDVI主要受以下因素影响:

  • 叶绿素吸收;
  • 冠层结构;
  • 土壤背景;
  • 大气条件;
  • 植被含水量;
  • 观测角度。

植被覆盖度则试图估计像元内植被面积比例。

3.3 植被覆盖度的应用

植被覆盖度可用于:

  • 草原退化;
  • 沙漠化监测;
  • 水土保持;
  • 森林和农田覆盖;
  • 生态质量评价;
  • 植被恢复;
  • 土地退化;
  • 城市绿地统计;
  • 气候与生态研究。

第四章 像元二分模型

4.1 模型假设

像元二分模型假设每个像元由两部分组成:

  1. 纯植被组分;
  2. 无植被组分。

传感器观测到的信号是两类组分的线性混合:

S = V F C ⋅ S v e g + ( 1 − V F C ) ⋅ S s o i l S = VFC\cdot S_{veg} + (1-VFC)\cdot S_{soil} S=VFC⋅Sveg​+(1−VFC)⋅Ssoil​

整理后得到:

V F C = S − S s o i l S v e g − S s o i l VFC = \frac{S-S_{soil}} {S_{veg}-S_{soil}} VFC=Sveg​−Ssoil​S−Ssoil​​

其中:

  • S S S为像元遥感信号;
  • S v e g S_{veg} Sveg​为纯植被组分信号;
  • S s o i l S_{soil} Ssoil​为纯土壤或无植被组分信号。

4.2 基于NDVI的改进模型

实际应用中常用NDVI替代一般遥感信号:

V F C = N D V I − N D V I s o i l N D V I v e g − N D V I s o i l VFC = \frac{NDVI-NDVI_{soil}} {NDVI_{veg}-NDVI_{soil}} VFC=NDVIveg​−NDVIsoil​NDVI−NDVIsoil​​

其中:

  • N D V I s o i l NDVI_{soil} NDVIsoil​表示无植被端元;
  • N D V I v e g NDVI_{veg} NDVIveg​表示纯植被端元。

若研究区可近似认为存在纯裸地和纯植被像元,则可写为:

V F C = N D V I − N D V I m i n N D V I m a x − N D V I m i n VFC = \frac{NDVI-NDVI_{min}} {NDVI_{max}-NDVI_{min}} VFC=NDVImax​−NDVImin​NDVI−NDVImin​​

4.3 模型的物理意义

当某像元的NDVI接近 N D V I s o i l NDVI_{soil} NDVIsoil​时,模型认为植被覆盖较低;当接近 N D V I v e g NDVI_{veg} NDVIveg​时,模型认为植被覆盖较高。

该模型的优势是:

  • 结构简单;
  • 参数少;
  • 计算速度快;
  • 便于大区域处理;
  • 适用于多种分辨率数据。

其限制包括:

  • 假设线性混合;
  • 端元依赖区域;
  • 高覆盖区易饱和;
  • 土壤背景差异影响明显;
  • 阴影和水体可能产生异常值。

第五章 NDVI端元参数的确定

5.1 为什么不能直接使用绝对最大值和最小值

影像中的绝对最大和最小NDVI可能来自:

  • 云影;
  • 水体;
  • 噪声;
  • 饱和像元;
  • 错误值;
  • 极少量异常目标。

因此,直接使用全影像极值会导致覆盖度偏差。

5.2 百分位法

常见方法是使用NDVI分布中的低百分位和高百分位作为端元。

例如:

  • 低端使用2%或5%分位数;
  • 高端使用95%或98%分位数。

这种方法比绝对极值稳定,但阈值仍需结合研究区地表组成确定。

5.3 分土地利用类型确定端元

若研究区包含林地、农田、城市、水体和裸地,不同类型的土壤背景和植被结构差异显著。

可以为不同土地利用类型分别统计:

  • N D V I s o i l NDVI_{soil} NDVIsoil​;
  • N D V I v e g NDVI_{veg} NDVIveg​;
  • 有效像元数量;
  • 置信范围。

然后构建空间变化的端元参数图。

这种做法比全区使用固定端元更合理。

5.4 基于实测覆盖度反推端元

若存在实测覆盖度,可利用以下关系估计端元:

N D V I s o i l = V F C m a x ⋅ N D V I m i n − V F C m i n ⋅ N D V I m a x V F C m a x − V F C m i n NDVI_{soil} = \frac{ VFC_{max}\cdot NDVI_{min} - VFC_{min}\cdot NDVI_{max} }{ VFC_{max}-VFC_{min} } NDVIsoil​=VFCmax​−VFCmin​VFCmax​⋅NDVImin​−VFCmin​⋅NDVImax​​

N D V I v e g = ( 1 − V F C m i n ) ⋅ N D V I m a x − ( 1 − V F C m a x ) ⋅ N D V I m i n V F C m a x − V F C m i n NDVI_{veg} = \frac{ (1-VFC_{min})\cdot NDVI_{max} - (1-VFC_{max})\cdot NDVI_{min} }{ VFC_{max}-VFC_{min} } NDVIveg​=VFCmax​−VFCmin​(1−VFCmin​)⋅NDVImax​−(1−VFCmax​)⋅NDVImin​​

这种方法利用已知覆盖度约束参数,比经验端元更有物理依据。

5.5 端元参数图

当不同土地利用类型采用不同端元时,可以分别生成:

  • N D V I s o i l NDVI_{soil} NDVIsoil​参数图;
  • N D V I v e g NDVI_{veg} NDVIveg​参数图。

像元级覆盖度计算为:

V F C ( x , y ) = N D V I ( x , y ) − N D V I s o i l ( x , y ) N D V I v e g ( x , y ) − N D V I s o i l ( x , y ) VFC(x,y) = \frac{ NDVI(x,y)-NDVI_{soil}(x,y) }{ NDVI_{veg}(x,y)-NDVI_{soil}(x,y) } VFC(x,y)=NDVIveg​(x,y)−NDVIsoil​(x,y)NDVI(x,y)−NDVIsoil​(x,y)​

这种空间化参数方法可减少全区统一阈值带来的偏差。


第六章 植被覆盖度计算与异常值处理

6.1 基本计算

若使用ENVI波段运算,表达式可概念化为:

(NDVI - NDVIsoil) / (NDVIveg - NDVIsoil)

6.2 值域约束

模型计算结果可能小于0或大于1。通常需要约束:

V F C f i n a l = { 0 , V F C < 0 V F C , 0 ≤ V F C ≤ 1 1 , V F C > 1 VFC_{final} = \begin{cases} 0, & VFC<0 \\ VFC, & 0\leq VFC\leq 1 \\ 1, & VFC>1 \end{cases} VFCfinal​=⎩ ⎨ ⎧​0,VFC,1,​VFC<00≤VFC≤1VFC>1​

6.3 水体、阴影与非植被区域

低NDVI不一定都是裸土,也可能是:

  • 水体;
  • 阴影;
  • 建筑;
  • 云影;
  • 深色地物。

因此,应结合土地覆盖掩膜排除不适合参与植被覆盖度解释的区域。

6.4 覆盖度分级

覆盖度可按实际需求划分等级,例如:

覆盖度等级
0–0.1极低覆盖
0.1–0.3低覆盖
0.3–0.5中低覆盖
0.5–0.7中高覆盖
0.7–1.0高覆盖

分级阈值应结合草原、森林、城市绿地或荒漠区具体背景调整。


第七章 植被覆盖度验证

7.1 为什么必须验证

像元二分模型简单,但端元选择和影像质量会显著影响结果。

验证数据可来源于:

  • 样方调查;
  • 无人机影像;
  • 高分辨率影像目视判读;
  • 地面照片;
  • 统计资料;
  • 同期生态调查。

7.2 常见评价指标

若存在实测值,可使用:

R M S E = 1 n ∑ i = 1 n ( y ^ i − y i ) 2 RMSE = \sqrt{ \frac{1}{n} \sum_{i=1}^{n} (\hat{y}_i-y_i)^2 } RMSE=n1​i=1∑n​(y^​i​−yi​)2 ​

M A E = 1 n ∑ i = 1 n ∣ y ^ i − y i ∣ MAE = \frac{1}{n} \sum_{i=1}^{n} |\hat{y}_i-y_i| MAE=n1​i=1∑n​∣y^​i​−yi​∣

还可以报告:

  • R 2 R^2 R2;
  • Bias;
  • 相对误差;
  • 不同覆盖等级精度。

7.3 尺度匹配

地面样方通常小于遥感像元,因此需要考虑:

  • 样方是否足够大;
  • 是否代表整个像元;
  • 是否与卫星成像时间一致;
  • 是否存在GPS误差;
  • 是否需要多样方聚合。

第八章 草原与沙漠化监测

8.1 草原退化的遥感表现

草原退化可能表现为:

  • 植被覆盖度下降;
  • 裸土比例增加;
  • NDVI下降;
  • 土壤亮度增强;
  • 地表纹理变化;
  • 沙地扩张;
  • 斑块破碎化增加。

8.2 多时相监测的基本条件

两个时相进行比较前,应保证:

  • 同季节或相近物候期;
  • 相似大气条件;
  • 统一辐射尺度;
  • 精确配准;
  • 相同投影;
  • 相同像元大小;
  • 一致分类规则;
  • 一致阈值或经过区域校准。

8.3 植被覆盖区提取

可根据覆盖度阈值提取植被区。

例如:

V e g e t a t i o n = { 1 , V F C > 0.3 0 , V F C ≤ 0.3 Vegetation= \begin{cases} 1, & VFC>0.3 \\ 0, & VFC\leq 0.3 \end{cases} Vegetation={1,0,​VFC>0.3VFC≤0.3​

阈值0.3属于经验值,应根据草原类型、年景和调查目标调整。

8.4 变化计算

设后期植被区为 V 2 V_2 V2​,前期植被区为 V 1 V_1 V1​,可计算:

Δ V = V 2 − V 1 \Delta V = V_2-V_1 ΔV=V2​−V1​

结果解释为:

值含义
-1植被减少或退化
0无变化
1植被增加

8.5 变化结果不能直接等同于沙漠化

植被减少可能来自:

  • 干旱;
  • 收割;
  • 放牧;
  • 火灾;
  • 土地利用变化;
  • 季节差异;
  • 影像质量问题。

只有结合裸土、沙地、土壤质地和长期趋势,才能判断是否属于沙漠化。

8.6 草原监测中的时序思维

两期差值只能描述端点差异。更稳定的方法应使用多年时间序列,分析:

  • NDVI趋势;
  • 覆盖度趋势;
  • 年际波动;
  • 持续退化;
  • 恢复时间;
  • 干旱响应;
  • 管理措施效果。

第九章 自然生态环境综合评价

9.1 单因子无法完整描述生态质量

生态环境质量不仅取决于植被,还与土壤、地形、水分和人类活动有关。

材料中的综合评价采用:

  • 植被覆盖度;
  • 土壤指数;
  • 坡度。

这种方法属于多因子综合评价。

9.2 植被因子

植被覆盖度通常作为正向指标:

  • 覆盖度越高;
  • 生态结构越稳定;
  • 水土保持能力通常越强;
  • 地表裸露越少。

但在城市或湿地环境中,高覆盖并不一定代表生态质量绝对更高,仍需结合类型解释。

9.3 土壤因子

材料使用裸土植被指数GRABS:

G R A B S = V I − 0.09178 × B I + 5.58959 GRABS = VI - 0.09178\times BI + 5.58959 GRABS=VI−0.09178×BI+5.58959

其中:

  • V I VI VI为穗帽变换绿度分量;
  • B I BI BI为穗帽变换亮度分量。

该指数通过组合植被与土壤亮度信息描述裸土状况。

9.4 坡度因子

坡度影响:

  • 地表径流;
  • 水土流失;
  • 土地利用适宜性;
  • 土壤稳定性;
  • 工程建设难度。

坡度通常作为负向指标,即坡度越大,综合生态适宜性评分越低。

9.5 指标归一化

不同指标量纲不同,必须统一尺度。

正向指标可采用:

X i ′ = X i − X m i n X m a x − X m i n X_i' = \frac{X_i-X_{min}} {X_{max}-X_{min}} Xi′​=Xmax​−Xmin​Xi​−Xmin​​

负向指标可采用:

X i ′ = X m a x − X i X m a x − X m i n X_i' = \frac{X_{max}-X_i} {X_{max}-X_{min}} Xi′​=Xmax​−Xmin​Xmax​−Xi​​

材料中使用1至10级编码,本质上也是归一化和等级化。

9.6 综合指数

材料中的综合模型为:

E = W 1 S v + W 2 S s + W 3 S t E = W_1S_v + W_2S_s + W_3S_t E=W1​Sv​+W2​Ss​+W3​St​

其中:

  • S v S_v Sv​为植被因子;
  • S s S_s Ss​为土壤因子;
  • S t S_t St​为地形因子;
  • W 1 W_1 W1​、 W 2 W_2 W2​、 W 3 W_3 W3​为权重。

示例权重为:

W 1 = 0.7 , W 2 = 0.2 , W 3 = 0.1 W_1=0.7,\quad W_2=0.2,\quad W_3=0.1 W1​=0.7,W2​=0.2,W3​=0.1

9.7 权重的含义

权重不是纯数学参数,而是对生态作用机制的判断。

不同地区权重应根据:

  • 生态主导因子;
  • 专家知识;
  • 统计相关性;
  • 主成分分析;
  • 层次分析;
  • 熵权法;
  • 实测生态质量;

进行确定或校准。

9.8 综合评价等级

综合指数可划分为:

等级含义
优生态结构稳定,自我恢复能力强
良生态功能较好,局部存在扰动
中生态系统受到一定破坏
差生态结构不合理,恢复能力弱

等级阈值必须结合研究区分布和管理标准解释。


第十章 坐标系与生态监测成果精度

10.1 地理坐标系与投影坐标系

地理坐标系使用经纬度描述位置,投影坐标系将地球曲面转换到平面。

生态面积、长度和坡度计算通常应在适合区域的投影坐标系中完成。

10.2 常见国内坐标基准

常见基准包括:

  • 北京54;
  • 西安80;
  • CGCS2000;
  • WGS84。

不同基准之间不能简单改标签,通常需要坐标转换参数。

10.3 面积统计的坐标要求

若直接在经纬度坐标系中计算面积,结果可能产生明显误差。

应使用:

  • 合适投影;
  • 米制单位;
  • 统一坐标;
    -明确分带;
  • 可靠行政边界。

第十一章 高分辨率城市绿地提取

11.1 高分辨率影像的优势

高分辨率影像可以识别:

  • 街旁绿地;
  • 公园;
  • 居住区绿化;
  • 道路绿化带;
  • 屋顶绿化;
  • 小型林地;
  • 河岸植被。

但高分辨率也带来:

  • 阴影多;
  • 类内差异大;
  • 光谱混乱;
  • 数据量大;
  • 单像元噪声显著。

11.2 WorldView-2波段优势

除蓝、绿、红和近红外外,WorldView-2还包含:

  • 海岸波段;
  • 黄色波段;
  • 红边波段;
  • 近红外2波段。

其中红边和近红外波段对植被识别和健康状况分析具有重要作用。

11.3 为什么采用面向对象方法

城市绿地不是单纯的绿色像元,而是具有空间连通性、形状和上下文的对象。

面向对象方法流程为:

高分辨率影像

影像分割

对象合并

计算NDVI、光谱和形状特征

规则分类

绿地对象

矢量后处理

11.4 分割与合并

分割尺度过小会造成绿地破碎;尺度过大会把绿地与道路、建筑或阴影合并。

分割参数应通过样区试验确定。

常见判断依据:

  • 树冠是否被过度分割;
  • 草坪是否保持完整;
  • 绿化带是否连续;
  • 建筑边界是否清晰;
  • 阴影是否独立成对象。

11.5 基于规则的绿地提取

材料中使用NDVI阈值和可见光波段阈值构建规则。

概念上可写为:

G r e e n = ( N D V I > T n d v i ) ∧ ( B a n d 2 < T b ) Green= (NDVI>T_{ndvi}) \land (Band_2<T_b) Green=(NDVI>Tndvi​)∧(Band2​<Tb​)

其中:

  • T n d v i T_{ndvi} Tndvi​为植被阈值;
  • T b T_b Tb​为排除高亮人工地物的阈值。

具体阈值不能跨影像直接复用,因为它受:

  • 大气;
    -太阳高度;
  • 地物材料;
  • 影像增益;
  • 季节;
  • 城市背景;

影响。

11.6 融合与正射校正

高分辨率绿地提取常先进行:

  • 全色—多光谱融合;
  • RPC正射校正;
  • 快速大气校正;
  • 面向对象提取。

融合提高空间细节,但可能改变光谱。若NDVI依赖融合影像,应检查红光和近红外是否保持。

11.7 矢量后处理

输出绿地矢量后通常需要:

  • 删除小碎片;
  • 填补孔洞;
  • 融合相邻图斑;
  • 裁剪行政边界;
  • 赋予绿地类型;
  • 检查拓扑;
  • 统计面积;
  • 关联街道或社区。

第十二章 城市绿地质量评价

12.1 仅提取面积还不够

绿地监测还可以进一步分析:

  • 人均绿地;
  • 绿地率;
  • 公园服务范围;
  • 绿地连通性;
  • 道路绿化连续性;
  • 热岛缓解;
  • 居住区可达性;
  • 生态廊道。

12.2 绿地率

绿地率可以表示为:

R g = A g A u r b a n × 100 % R_g = \frac{A_g} {A_{urban}} \times100\% Rg​=Aurban​Ag​​×100%

其中:

  • A g A_g Ag​为绿地面积;
  • A u r b a n A_{urban} Aurban​为城市统计范围面积。

12.3 绿地连通性

仅有较高绿地率,不代表生态结构一定合理。碎片化绿地的生态功能可能低于连续绿地。

可结合:

  • 斑块数量;
  • 平均斑块面积;
  • 边界密度;
  • 最近邻距离;
  • 连通度;
  • 景观聚集度;

评价绿地格局。


第十三章 MODIS气溶胶遥感监测

13.1 气溶胶是什么

气溶胶是液态或固态微粒在空气中的悬浮体系。

其来源包括:

  • 沙尘;
  • 海盐;
  • 硫酸盐;
  • 硝酸盐;
  • 有机碳;
  • 黑碳;
  • 烟尘;
  • 工业排放;
  • 生物质燃烧。

13.2 气溶胶光学厚度

气溶胶光学厚度AOD表示气溶胶对太阳辐射衰减的综合程度。

根据Beer-Lambert定律,可概念化表示为:

I = I 0 e − τ I = I_0e^{-\tau} I=I0​e−τ

其中:

  • I 0 I_0 I0​为入射辐射;
  • I I I为透过辐射;
  • τ \tau τ为总光学厚度。

AOD越大,说明气溶胶对辐射的消光作用越强。

13.3 MODIS数据结构

MODIS L1B数据通常包含:

  • 反射率数据集;
  • 发射率数据集;
  • 太阳角度;
  • 卫星角度;
  • 地理定位信息;
  • 质量信息。

不同数据集可能具有不同空间分辨率和缩放因子,因此不能直接叠加。

13.4 MODIS几何校正

处理流程中需要分别校正:

  • 反射率数据;
  • 发射率数据;
  • 角度数据。

角度数据还需要重采样到与反射率一致的空间分辨率。

13.5 角度数据缩放

MODIS角度数据常以整数存储,并使用缩放因子恢复真实角度。

若原始值扩大100倍,则:

θ r e a l = θ s t o r e d × 0.01 \theta_{real} = \theta_{stored}\times0.01 θreal​=θstored​×0.01

忽略缩放因子会使太阳和传感器几何完全错误。

13.6 云检测

云会显著影响气溶胶反演。

云检测需要排除:

  • 高亮厚云;
  • 卷云;
  • 云边缘;
  • 云影;
  • 雪和冰;
  • 高反射地表。

云掩膜质量直接影响AOD空间分布。


第十四章 暗像元法反演气溶胶

14.1 暗像元法的基本思想

暗像元法利用浓密植被在某些可见光波段反射率较低、在短波红外中相对稳定的特征,估算地表反射率,再将观测信号中多余部分归因于大气气溶胶。

传感器表观反射率可概念化表示为:

ρ T O A = ρ p a t h + T ⋅ ρ s u r f a c e \rho_{TOA} = \rho_{path} + T\cdot\rho_{surface} ρTOA​=ρpath​+T⋅ρsurface​

当 ρ s u r f a c e \rho_{surface} ρsurface​可以估计时,就可以反演气溶胶影响。

14.2 暗像元选择

暗像元通常来自:

  • 浓密植被;
  • 深色地表;
  • 非云区域;
  • 非水体;
  • 非高亮裸地。

若研究区植被稀少,暗像元法适用性会下降。

14.3 为什么冬季效果较差

冬季可能存在:

  • 植被覆盖低;
  • 裸土和建筑比例高;
  • 雪和冰;
  • 地表反射率不满足暗像元假设;
  • 湿度和气溶胶类型变化。

因此,暗像元法在冬季和干旱区可能产生较大误差。

14.4 查找表方法

查找表通过辐射传输模型预先模拟:

  • 不同AOD;
  • 太阳天顶角;
  • 观测天顶角;
  • 相对方位角;
  • 气溶胶模型;
  • 地表反射率;
  • 波段响应;

对应的表观反射率。

实际反演时,在查找表中寻找与观测条件最接近的AOD。

14.5 查找表的适用范围

一套查找表不能无条件用于所有时空场景。

需要明确:

  • 季节;
  • 气溶胶类型;
  • 高程;
  • 波段响应;
    -太阳角度范围;
  • 观测角度范围;
  • 水汽和臭氧参数。

第十五章 气溶胶反演结果验证

15.1 地面观测

AOD结果可使用地面太阳光度计验证。

验证应考虑:

  • 卫星过境时间;
  • 地面站时间窗口;
  • 空间匹配半径;
  • 云筛选;
  • 波长转换;
  • 地表类型。

15.2 波长差异

若地面观测与卫星AOD波长不同,可利用Ångström指数转换:

τ ( λ 2 ) = τ ( λ 1 ) ( λ 2 λ 1 ) − α \tau(\lambda_2) = \tau(\lambda_1) \left( \frac{\lambda_2}{\lambda_1} \right)^{-\alpha} τ(λ2​)=τ(λ1​)(λ1​λ2​​)−α

其中:

  • α \alpha α为Ångström指数;
  • τ ( λ ) \tau(\lambda) τ(λ)为对应波长AOD。

15.3 误差来源

AOD误差可能来自:

  • 云污染;
  • 地表反射率估计错误;
  • 气溶胶模型不匹配;
  • 几何角度错误;
  • 查找表范围不足;
  • 像元混合;
  • 高反射地表;
  • 数据缩放错误。

第十六章 五类专题的统一技术路线

16.1 植被覆盖度反演

多光谱影像

辐射与大气处理

NDVI

确定NDVIsoil和NDVIveg

像元二分模型

覆盖度约束与分级

验证与统计

16.2 草原与沙漠化监测

多时相影像

统一预处理

植被覆盖度

提取植被覆盖区

变化检测

退化与恢复分析

长期趋势判断

16.3 自然生态环境评价

遥感影像与DEM

植被覆盖度

土壤指数

坡度

归一化

加权综合指数

生态等级

16.4 城市绿地提取

高分辨率影像

融合与正射校正

大气校正

影像分割

NDVI、光谱和形状规则

绿地对象

矢量后处理

16.5 MODIS气溶胶反演

MODIS L1B

几何校正

反射率、发射率和角度合成

角度缩放

云检测

暗像元识别

查找表反演

AOD结果与验证


第十七章 生态监测中的常见误区

17.1 NDVI可以直接代表植被覆盖度

错误。NDVI是光谱指数,覆盖度是比例参数,两者需要通过模型建立关系。

17.2 端元取全影像最大值和最小值即可

错误。极值可能来自噪声、云影和异常地物。

17.3 覆盖度下降就表示沙漠化

错误。还可能是干旱、收割、放牧、火灾和季节差异。

17.4 综合生态指数越复杂越科学

错误。因子和权重必须具有生态意义,并经过验证。

17.5 城市绿地只需要NDVI阈值

错误。阴影、屋顶、人工草坪和水体会造成混淆,高分影像还需要对象和空间规则。

17.6 融合后影像更清晰,NDVI一定更准确

错误。融合可能改变红光和近红外光谱值。

17.7 MODIS角度数据可以直接使用

错误。角度数据通常包含缩放因子。

17.8 暗像元法适合所有季节和地区

错误。冬季、干旱区、高反射地表和低植被区适用性较差。

17.9 地图看起来合理就不需要验证

错误。生态参数和反演结果必须通过独立数据验证。


第十八章 术语表

术语说明
VFC/FVC植被覆盖度
NDVI归一化植被指数
像元二分模型将像元视为植被与非植被组分线性混合的模型
NDVIsoil无植被端元NDVI
NDVIveg纯植被端元NDVI
GRABS裸土植被指数
BI穗帽变换亮度分量
VI穗帽变换绿度分量
AOD气溶胶光学厚度
暗像元法利用低地表反射率区域反演气溶胶的方法
查找表预先模拟观测条件与参数关系的数据表
QUAC快速大气校正方法
FLAASH基于辐射传输模型的大气校正方法
面向对象提取以分割对象为基本单位的信息提取方法
分割尺度决定对象大小和边界的参数
变化检测识别多时相地表状态差异
综合指数将多个生态因子加权组合的评价指标
归一化将不同量纲指标转换到统一尺度
Ångström指数描述AOD随波长变化关系的参数
MODIS中分辨率成像光谱仪

mall 电商项目完整部署---实战 mall 后端 SpringBoot,mall‑admin‑web 后台管理 Vue 前端,mall‑app‑sc 移动端 uni‑app 前端;修改数据库、redis、minio、logstash 地址,容器之间 link 别名:mysql 容器别名。🔔排坑:登录报跨域 / 接口报错,检查 prod.env.js 地址,确认 8080 后端服务正常运行。能打开页面,接口列表展示出来,后端部署完成。,所以打包的时候会自动构建镜像,很慢而且你暂时不需要。,容器数据持久化,删除容器数据不会丢失。 阅读详情

相关推荐

ffmpeg 4.4 cenc-aes-ctr 加解密 MP4 工程性质分析

ffmpeg 4.4 cenc-aes-ctr 加解密 MP4 工程性质分析

Dream_Dynasty的博客 3045

Flink 学习日志 03:自定义一个SourceFunction,完成对Redis中数据的读取

本文介绍了在Flink中自定义SourceFunction读取Redis数据的完整流程:通过Jedis客户端连接WSL2中的Redis服务,编写RedisSource实现循环读取所有Key-value并逐条输出,结合Thread.sleep()实现每秒一条的流式效果。关键步骤包括环境准备、添加Jedis依赖、配置IDEA运行时加载provided依赖,并注意生产环境中应避免使用阻塞的keys *,改用scan命令。最终实现动态读取、实时打印,验证了Flink与Redis的数据集成能力。

木木的学习笔记 423

JeecgBoot漏洞利用工具解析:一键getshell原理与安全加固实践

在应用安全领域,远程代码执行(RCE)和文件上传漏洞是Web系统面临的高危威胁,其核心原理在于对用户输入的控制不足或安全校验缺失。攻击者通过构造恶意请求,可绕过系统边界,在服务器上执行任意命令或上传恶意文件,从而获取系统控制权。这类漏洞的技术价值在于揭示了应用层安全防护的薄弱环节,其利用方式自动化、集成化趋势明显,极大地提升了攻击效率。在低代码开发平台等快速构建的业务系统中,此类风险尤为突出。应用场景广泛覆盖安全测试、渗透评估与应急响应。本文聚焦于JeecgBoot平台的历史高危漏洞,深入解析其SQL注入、

weixin_30787531的博客 552

Day3 - Flink 窗口与时间语义:从滚动窗口到迟到数据处理

水位线是一个特殊的时间戳,代表"小于这个时间戳的数据都已经到齐了比如水位线推进到 10:00:05,就告诉 Flink:10:00:05 之前的数据应该都到了,如果还有晚于这个时间到达但时间戳小于 10:00:05 的数据,就是"迟到数据"。窗口大小固定,窗口之间不重叠,数据只属于一个窗口。时间线: 0 10 20 30 40 50 60窗口1: [--------)窗口2: [--------)窗口3: [--------)窗口4: [--------)

YEEMOLY的博客 251

YonSuite如何新增用户并分配角色角色授权图文教程

用友YonSuite产品基础教学

ct341688的博客 619

告别抖动与滞后:深入博能传动伺服DSC功能,用105报文提升多轴同步性能

本文深入解析博能传动伺服DSC功能如何通过105报文实现多轴同步性能提升,有效解决工业自动化中的抖动与滞后问题。详细介绍了DSC技术原理、105报文实现机制及多轴同步工程实践,帮助工程师优化伺服系统性能,实现绝对齿轮同步。

weixin_30394981的博客 184

智能问诊的流程

本文介绍基于FastAPI构建的医疗多Agent系统MediAgent,通过症状分析、诊断推理与治疗建议等模块协同,实现从问诊到报告生成的完整工作流。系统采用前后端分离架构,结合MySQL、RabbitMQ与Docker实现工程化部署,支持异步化验结果接入。虽当前状态暂存于内存,但已具备可扩展性,未来可通过Redis持久化、消息队列联动及医疗安全机制完善,应用于慢病管理、辅助决策等场景,推动AI在医疗流程中的落地应用。

Joey 的工作空间~ 774

Cortex-M TrustZone实战:从硬件隔离到安全交互的嵌入式开发指南

在嵌入式系统开发中,硬件安全隔离是构建可信执行环境(TEE)的核心技术。其原理是通过处理器硬件机制,将系统资源划分为安全与非安全域,实现内存、外设和代码的强制访问控制。这项技术的核心价值在于,它能在单一芯片上为资源受限的物联网(IoT)设备、智能传感器等提供芯片级的安全基础,无需依赖外置安全元件,从而降低BOM成本并简化设计。典型的应用场景包括保护设备密钥、实现安全启动、构建安全服务(如加密、密钥管理)以供应用程序调用。本文聚焦于Arm Cortex-M系列微控制器上的TrustZone技术实现,针对开发中

weixin_42639246的博客 300

行业研究报告采购避坑指南:从需求拆解到机构选型的实操逻辑

先说结论:买行业研究报告这件事,坑不在“哪家便宜”,而在“你拿到的数据能不能回答你的问题”。作为一个经常帮企业做采购决策的顾问,我见过太多花了几万块买回一堆PDF却用不上的案例。这篇文章把采购流程拆成五步,每一步告诉你该看什么、该避什么,尽量帮你少交学费。

2601_96585642的博客 225

Chrome红页误报?手把手教你通过Search Console快速解除安全警告

本文详细解析了Chrome红页误报的常见原因及解决方案,重点介绍如何通过Google Search Console快速解除安全警告。内容涵盖误报诊断、技术修复、申诉技巧及长效预防措施,帮助网站管理员有效应对Google Safe Browsing的误判问题,确保网站安全状态快速恢复。

a1b2c的博客 995

我的世界模组开发——前言

避坑:在开发过程中一定不要过分依赖AI,虽说现在的人工智能很厉害,但是也有很大的弊端,很多AI生成的代码错误百出,还得自己修改,而且像deepseek这种AI的训练数据还停留在一年前,给出的maven仓库地址都是过时的甚至是错误的,开发模组就是要多学习,多看源码,这样开发模组的水平才能越来越高。总而言之,build.gradle文件在我的世界Forge模组开发中扮演了重要的角色,它定义了项目的构建过程、依赖关系和编译选项,为开发者提供了灵活而强大的构建工具。开发工具推荐使用idea,因为可以使用通灵义码。

qq_24992377的博客 2594

大数据与审计专业哪些证书企业比较认可?

大数据与审计专业选证书,首先要明确就业方向:是进入会计师事务所做财务审计,还是进入企业从事内部审计、风险控制,或向数据审计、信息系统审计发展。不同岗位看重的专业能力不同,证书的价值也应放在具体招聘要求中判断。面向2027年的求职与备考规划,值得重点考虑的有六类:初级会计专业技术资格、注册会计师、审计专业技术资格、CDA数据分析师、国际注册内部审计师和国际注册信息系统审计师。选择时,宜采用“财审基础+数据能力+岗位专长”的组合,而不是单纯追求持证数量。

YangYang9YangYan的博客 296

uniapp - H5人脸识别认证与活体检测功能(纯前端免费方案+微信sdk人脸识别)完整源码,微信公众号网页/h5端人脸识别功能人脸核身(微信方案提供后端spring boot+thinkphp源码

uniapp实现人脸认证,uniAPP h5手机拍照+人脸识别,免费人脸识别插件,前端人脸识别解决方案,uniapp H5人脸识别带前后摄像头,人脸识别在H5页面进行采集,微信公众号H5人脸核身,H5端人脸识别功能实现,免费h5人脸识别插件方案,微信公众号H5网页实现人脸识别,uniapp人脸监测,uniapp人脸生物监测,使用uniapp开发H5人脸采集功能/人脸识别功能,uniapp项目—h5端人脸采集功能实现,uniapp vue3人脸核身,H5活体检测认证,人脸认证张嘴摇摇头指引动作,人脸核身 微信

🏆 前端领域优质创作者 3809

深入浅出聊聊 Linux 的 PageCache 与系统极致性能

本文深入浅出解析了 Linux PageCache 的核心机制,揭示其如何通过“延迟写、预读读”提升系统性能。以监控告警为切入点,阐明 free -h 中 buff/cache 并非真正占用,真正可用内存由 available 决定。结合 Kafka、Elasticsearch、MySQL 等实战场景,说明中间件如何依赖 PageCache 实现极致性能。强调不可手动清理缓存,并建议合理调优脏页参数,避免卡顿。最终指出:内存满≠危机,而是系统高效运行的体现。

小小默:进无止境 483

STM32F407中文手册

STM32f407中文手册 本参考手册面向应用开发人员,提供有关使用 STM32F405xx/07xx、STM32F415xx/17xx、 STM32F42xxx 和 STM32F43xxx 微控制器存储器与外设的完整信息。STM32F405xx/07xx、STM32F415xx/17xx、STM32F42xxx 和 STM32F43xxx 构成一个微 控制器系列,各产品具有不同的存储器大小、封装和外设。有关订购信息以及器件的机械与电气特性,请参见数据手册。有关 ARM Cortex™-M4F 内核的信息,请参见《Cortex™-M4F 技术参考手册》。

假人有限元模型

假人模型,车辆碰撞假人有限元模型,适合做乘员安全的研究者利用

ubuntu安装hadoop系统

本文详细介绍了基于VMware克隆虚拟机搭建三节点Hadoop分布式集群的全过程:使用同一Ubuntu镜像克隆出master、worker1、worker2三台虚拟机,依次配置主机名、固定IP(如192.168.82.10/11/12)、SSH服务及免密登录,确保节点间网络互通。随后安装Java,下载配置Hadoop并分发至各节点,最终初始化启动集群,验证其存储与计算功能,实现高效的大数据环境部署。

qq_43592674的博客 338

文档:亲测可用解决vnc灰屏,黑屏,鼠标大黑叉情况

主要解决vnc-viewer灰屏,黑屏,鼠标大黑叉等现象,详细见文档。

Elasticsearch:稠密向量字段类型

dense_vector字段用于存储稠密向量,支持k近邻搜索(kNN),默认启用索引并使用余弦相似度。可接受浮点数数组或Base64编码字符串输入,支持自动量化(如int8、int4、bbq)以降低内存占用。通过index_options配置索引类型(如hnsw、bbq_disk),并可启用过采样与重新评分提升精度。向量默认不存储于_source中以节省空间,可通过fields或_source.exclude_vectors显式返回。支持bit向量及GPU加速索引,适用于高维向量场景。

Elastic 中国社区官方博客 288

TMS_WebGMaps_1.9.8.0_Full_Source

TMS_WebGMaps_1.9.8.0_Full_Source 全源码 delphi cbuilder

Day2 - Flink DataStream API 与转换算子:从 Source 到 Sink

实现或@Override// 发射数据@Override// 使用用自定义逻辑聚合,比 sum/min/max 更灵活。// 求每个用户的总消费金额reduce 的输入和输出类型必须相同。Flink 会把满足条件的相邻算子合并成一个 Task,在同一个线程里执行,减少线程切换和网络开销。合并条件上下游并行度相同;分区策略是 Forward;没有被用户手动断开。

YEEMOLY的博客 162

【消息队列】ACK确认了什么:消息队列可靠投递的完整链路

本文剖析消息系统中ACK机制在生产者、Broker、消费者三段链路中的作用,指出即使三段确认均正确配置,系统仍只能提供“至少一次”可靠性。原因在于重试、崩溃、重平衡等场景会导致重复投递,而“恰好一次”需依赖业务侧幂等性保障。关键结论:ACK仅保证不丢,不重复需靠唯一键、去重表或状态机等业务设计实现。

2503_94545876的博客 699

kafka指定分区消费

【代码】kafka指定分区消费。

hy_208的博客 64

数据科学与大数据技术专业有哪些比较好考的证书?

如果目前还没有明确方向,建议先完成Python与SQL基础学习,再选择一张与目标岗位最接近的证书。不要为了增加简历条目,在同一阶段同时准备多个不同技术方向。以数据分析方向为例,可以参考2025版CDA证书大纲安排基础学习:先掌握业务数据结构、数据库应用和统计分析,再进入指标体系、看板和用户分析。具备一定项目经验后,再学习二级涉及的统计检验、模型评估及数据治理。2025版教材介绍强调了思维、技术和实践的结合,适合用于组织学习;2027年的具体考试范围、报名条件和安排,则应依据届时发布的信息制定计划。

YangYang9YangYan的博客 331

如何在用户认证系统中使用 ElasticSearch 提升

本文介绍了如何在一个典型的企业级应用——用户认证与权限管理系统内引入并运用 ElasticSearch 技术来提高查找效率并降低数据库负载压力。读者可以通过了解其基础架构原理、正确配置方式以及典型应用场景从而掌握该工具的基本用法。如果你正在寻找一种更灵活高效的检索方案用于自己的项目当中,请考虑从以下几个方面入手实践: 1. 确认是否有必要引入新的全文搜索引擎;2. 设计好适用于该项目的数据模型以及索引规则;3. 持续优化你的 DSL 查询语句并监控运行状态;

wuchangjian的专栏 176

Flink SQL 窗口表值函数实战:滚动、滑动、累计与会话窗口

在 Flink SQL 中,流数据是无限的,要统计“最近一段时间”的数据,就必须把无界流切分成一个个有界窗口。Flink 1.13 之后引入了窗口表值函数(Window TVF),把窗口本身当作一张表来使用,通过TABLE(...)语法把原表转换成带窗口信息的表,再配合GROUP BY完成聚合。相比早期的语义统一:窗口是一张表,可以像普通表一样查询、过滤、关联。结果可读:窗口表会自动带上window_end三个字段,方便直接输出。性能更好:底层基于窗口聚合算子优化,支持切片共享等机制。

2401_83857914的博客 362

Flink 知识大总结:从核心概念到生产实践

Flink 的知识体系可以概括为「一个核心、两大抽象、三层保障」:以流为核心,以 DataStream 和 Table/SQL 为两大编程抽象,以状态管理、Checkpoint 容错和端到端一致性为三层保障。掌握这些主干知识后,再结合具体业务场景深入调优,就能在生产环境中稳定发挥 Flink 的实时计算能力。

2401_83857914的博客 328

数据库教程FGMT79:国产数据库-TiDB分布式数据库工程师培训实战教程(PB级大数据平台、大规模分布式集群架构)

本文为国产分布式数据库TiDB实战教程,系统讲解其架构原理与生产级部署。涵盖从环境准备、单机测试到十节点OLTP/HTAP集群部署、扩缩容、两地三中心容灾等全流程操作。基于64G内存8CPU配置,以TiUP工具实现自动化管理,深入解析存储、计算、调度三大核心机制。强调兼容MySQL协议,支持金融级高可用与弹性扩展,适用于PB级大数据平台与HTAP场景。适合DBA、运维及架构师快速掌握TiDB生产实践,无需重写业务代码即可实现平滑迁移。

风哥oracle/mysql/pg技术教程 160

最新大数据毕业设计选题推荐-基于大数据的蛋黄酱零售价格区域与时间多维数据可视化分析-大数据-Spark-Hadoop-Bigdata

本文介绍《基于大数据的蛋黄酱零售价格区域与时间多维数据可视化分析》系统,适用于大四计算机专业毕业设计。系统采用Hadoop+Spark进行分布式数据处理,结合Django/Spring Boot后端与Vue+ECharts前端,实现从数据采集、清洗、存储到多维分析与可视化展示的全流程。功能涵盖时间走势、规格对比、区域分析等,支持按时间、区域、规格多维度下钻,界面清晰,技术路线完整,适合作为大数据方向实践课题,提供源码及文档支持。

IT研究室的博客 380

大数据Spark(八十五):Stage及Spark计算模式

本文详解Spark中Stage的划分机制与Pipeline计算模式。Stage依据RDD宽依赖(Shuffle)划分为ShuffleMapStage和ResultStage,同一Stage内通过迭代器串联多个RDD,实现数据内存流转,减少落盘开销,提升效率。并行度由最后一步RDD的分区数决定,可通过Shuffle算子调优。结合SparkUI可直观分析性能瓶颈,是优化作业的关键基础。

Lansonli(蓝深李)的博客 764
上一篇: 【遥感学习整理 03】高分一号与Hyperion数据预处理:从辐射定标到融合与高光谱校正
下一篇: 【遥感学习整理 05】水质火点森林农业与矿物蚀变遥感监测
ThsPool
博客等级 码龄3年 2万+粉丝 93原创
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付:10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值