0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

基于 Landsat 与地物波谱的稀土矿尾砂地植被

莱森光学 来源:莱森光学 作者:莱森光学 2024-12-10 15:09 120次阅读

稀土开采产生大量堆积尾砂,出现强烈水土流失过程,伴生严重水环境和地质灾害。采用工程和非工程措施结合手段开展稀土尾砂地植被修复是尾砂地灾害防治的常见手段。评价矿区植被恢复效果及成因可为尾砂治理措施的规划设计提供理论依据。

NDVI时程变化

一、引言

稀土开采产生大量堆积尾砂,出现强烈水土流失过程,伴生严重水环境和地质灾害。采用工程和非工程措施结合手段开展稀土尾砂地植被修复是尾砂地灾害防治的常见手段。评价矿区植被恢复效果及成因可为尾砂治理措施的规划设计提供理论依据。本文以江西省寻乌县离子型轻稀土尾砂区为研究区,收集研究区1989-2015年Landsat数据和寻乌站气象数据,实地调查研究区典型尾砂适生植被和堆积尾砂的地物波谱数据,解算典型植被的NDVI值,分析矿区NDVI的年际和年内时程变化规律。

二、基于Landsat数据的NDVI时程变化

在利用遥感技术对植被信息的研究中,NDVI是植被指数中运用最为广泛的一种,它能很好地反映了植被生长状态以及植被空间分布情况,与植被分布密度呈线性相关。许多研究表明NDVI与植被覆盖度、叶面积指数(LAI)和光合作用等植被参数有关,NDVI的时间变化曲线能有效地表现出植被的季节变化、年际变化和人为活动影响变化。因此,在监测废弃矿区植被恢复情况的过程中NDVI起到有效地作用。本文选取归一化植被指数(NDVI)来反映稀土矿尾砂地植被变化情况,其计算形式为近红外波段(NIR)与红光波段(R)的差值比上这两个波段数值之和,公式为:

wKgZO2dX6TGAcuSnAAAGnLpIa0s653.png

表1Landsat TM/OLI NDVI 计算公式

wKgZPGdX6TGAKnqOAAAej8SLHJ4698.png

研究区NDVI时程变化分析分为两个部分进行,第一部分为整体矿区NDVI时程变化分析,第二部分选取植被修复时间较早的“稀土尾砂试验区”双茶亭矿区分析NDVI时程变化规律。

2.1整体矿区NDVI时程变化

基于ENVI软件中BandMath 计算NDVI结果,将整体矿区NDVI值提取出来并导出成ASCII文件,方便之后分析使用。计算NDVI平均值,诊断NDVI时序趋势,揭示整体矿区NDVI时程变化规律,见图1。图1表明:1989-2015年整体矿区NDVI值在0.12-0.5间变化,平均值为0.31;NDVI最大值和最小值所在年份分别为1995年和2013年,分别为0.5和0.12,变化幅度为0.38;NDVI值相对时间的线性回归决定系数R2为0.79,取显著性水平为0.05,相应临界值R205.0 为0.08,决定系数大于相应临界值,线性回归趋势显著,NDVI值的倾向率为-0.08(10a)-1,表明NDVI值显著减小;NDVI值的Mann-kendall趋势分析检验值Z为-6.68,取显著性水平为0.05,相应临界值Z05.0 为1.96,检验值的绝对值Z大于相应临界值,Mann-kendall趋势检验显著,表明整体年NDVI值呈显著性下降趋势。

wKgZO2dX6TKAeTvzAAA67NyGACM376.png

图1 整体矿区 NDVI 时程变化

NDVI时程变化可分为三个阶段,第一阶段为稳定期,时间为1989-2005年,此阶段NDVI值相对时间的线性回归决定系数R2为0.052,取显著性水平为0.05,相应临界值R205.0 为0.157,决定系数小于相应临界值,NDVI变化趋势不显著,平均NDVI为0.43,整体稀土矿区存在一定的开采,但是整体NDVI值未大幅度减小;第二阶段为剧烈开采期,时间为2006-2008年,此阶段NDVI值相对时间的线性回归决定系数R2为0.69,取显著性水平为0.05,相应临界值R205.0 为0.658,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.3(10a)-1,平均NDVI为0.29,较稳定期下降了0.14,此阶段稀土矿区进行剧烈的采矿活动,人类活动对稀土矿区影响较大,严重破坏了矿区植被,导致NDVI值迅速下降;第三阶段为植被修复期,时间为2009-2015年,此阶段NDVI值相对时间的线性回归决定系数R2为0.395,取显著性水平为0.05,相应临界值R205.0 为0.197,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.046(10a)-1,整个时间段NDVI呈现较低值,平均NDVI为0.19,较剧烈开采期下降了0.1,寻乌县2009年对稀土矿区进行植被修复治理,使研究区内的2009年和2010年的植被数量有所增加,对于整体矿区恢复效果一般。

2.2双茶亭矿区NDVI时程变化

在整体矿区NDVI影像数据中提取以种植桉树为主的双茶亭矿区NDVI,计算NDVI平均值,诊断NDVI时序趋势,揭示双茶亭矿区NDVI时程变化规律,见图2。按照整体矿区NDVI变化阶段对双茶亭矿区NDVI时程变化进行分析。

wKgZPGdX6TKAKddMAAA5JwTQF8s360.png

图2 双茶亭矿区 NDVI 时程变化

结果表明:1989-2015年双茶亭矿区NDVI值在0.17-0.65间变化,平均值为0.4;NDVI最大值和最小值所在年份分别为2015年和2008年,分别为0.65和0.17,变化幅度为0.48;稳定期NDVI值相对时间的线性回归决定系数R2为0.321,取显著性水平为0.05,相应临界值R205.0 为0.156,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.049(10a)-1,平均NDVI为0.44,稀土开采对双茶亭矿区内植被存在一定影响,植被呈现减少趋势;剧烈开采期NDVI值相对时间的线性回归决定系数R2为0.732,取显著性水平为0.05,相应临界值R205.0 为0.57,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.63(10a)-1,平均NDVI为0.28,较稳定期减少0.16,和整体矿区情况类似,剧烈稀土开采活动严重破坏双茶亭矿区内植被,导致植被迅速减少;植被修复期NDVI相对时间的线性回归决定系数R2为0.884,取显著性水平为0.05,相应临界值R205.0 为0.187,决定系数大于相应临界值,NDVI呈显著上升趋势,倾向率为0.25(10a)-1,平均NDVI为0.39,较剧烈开采期增加0.11,双茶亭植被修复效果较好,对比整体矿区可得,油茶配置模式植被修复治理效果较差,桉树配置模式效果较好。

三、基于地物波谱的典型地物NDVI特征值分析

基于研究区各典型地物的地物波谱特性,利用红光波段和近红外波段的平均反射率计算相应地物的NDVI,见图3。图3表明:在研究区典型地物中,尾砂相应NDVI值为0.03,芒萁为0.69,芭茅为0.53,芒草为0.52,油茶为0.69,木荷为0.8,竹子为0.79,马尾松为0.75,桉树为0.6,樟树为0.52,其中木荷NDVI值最大,尾砂NDVI值最小,尾砂存在一定的背景值。由于研究区地物波谱数采集时间为2019年1月与7月,分别属于冬季和夏季,则获取各地物波谱特性计算出的NDVI会有一定的波动范围,其中波动范围较大的为芒萁、芭茅和芒草,这三种植被属草本类植被,说明对于NDVI,草本类植被受季节变化影响较大,乔木类植被对相对较小,即草本NDVI在年内变化更为剧烈,乔木类植被NDVI年内较为稳定。

wKgZO2dX6TKAXfVUAAAqHcZGfKg894.png

图3典型地物NDVI

四、基于小波分解的NDVI时程变化

4.1小波分析理论概述

傅里叶变换是最常用的分析信号工具,用于推导不同幅度和频率下特定正弦波函数的组合。但是傅里叶变换存在不足之处,在分析信号特征时会忽略信号的时间信息,通过变换结果无法判断某个信号发生的时间。以傅里叶变换为数学基础,理论物理学家A·Grossman建立了小波变换的理论体系。小波变换分析信号时可通过小波基函数的变换分析信号的局部特征,信号分析结果具有比傅里叶变换更良好的局部特性。分解后的信号呈现低频和高频特性两部分内容,低频信号具有较高的频率分辨率和较低的时间分辨率,高频信号具有较高的时间分辨率和较低的频率分辨率,所以有着“数学显微镜”之称。正是小波变换拥有着这些特性,使小波变换在信号处理、数值计算和流体力学等领域得到了广泛的应用。

4.2基于小波分解的NDVI多尺度时程变化

为了更好研究NDVI的时程变化规律,基于Matlab软件,本文将NDVI变化曲线作为一种信号曲线,对其进行小波分解,分析其低频信号与高频信号的特点。根据NDVI数据特点,采用db3小波函数对整体矿区NDVI和双茶亭NDVI进行分解,结果如图4和图5所示。

wKgZPGdX6TOAJIQEAABpKVNrq80718.png

图4 整体矿区 NDVI 信号和小波分解部分注:(a)表示原始信号,(b)表示低频信号,(c)表示高频信号。

(1)低频信号按照NDVI变化阶段对低频信号进行分析。图4表明:稳定期整体矿区低频NDVI值处于较高水平,进入剧烈开采期后低频NDVI值迅速下降,在植被修复期保持稳定;低频信号整体趋势与整体矿区NDVI时程变化趋势一致,说明低频信号表示NDVI的年际变化特征,将年内波动信息滤除,呈现与原始信号一致的三个变化阶段。图5表明:双茶亭矿区低频NDVI值呈现微下降-剧烈下降-上升的三个阶段变化,与双茶亭矿区NDVI时程变化趋势一致。

(2)高频信号高频信号表示为NDVI的年内变化,按照NDVI变化阶段对高频信号进行分析。图4表明:稳定期整体矿区高频NDVI值呈波动变化,波动幅度较大,说明在此阶段整体矿区年内NDVI变化较为剧烈;剧烈开采期整体矿区高频NDVI值也呈波动变化,但是波动幅度较小,说明在此阶段矿区年内NDVI呈波动变化,但是由于人类活动较为剧烈,使植被数量迅速较少,导致NDVI变化幅度变小;植被修复期整体矿区高频NDVI值在此阶段初期呈水平变化,后期呈波动变化,说明研究区年内NDVI在植被修复期初期年内变化较小,这是由于剧烈开采期过后植被破坏过于严重,导致矿区内植被数量稀少,年内植被变化较小,后期由于植被修复措施的实施,植被得到一定的恢复,然后呈现出了年内波动变化。图5表明:双茶亭矿区高频NDVI值三个时期都呈波动变化,这是由于在剧烈开采后,植被措施的迅速实施,使双茶亭矿区内植被得到良好的补充,导致在剧烈开采期之后、植被恢复期初植被表现出良好的年内波动变化。

wKgZO2dX6TSAfrfpAABmzq6r3Xw438.png

图5双茶亭矿区NDVI信号和小波分解部分注:(a)表示原始信号,(b)表示低频信号,(c)表示高频信号。

结合草本植被NDVI在年内变化更为剧烈,乔木类植被NDVI年内较为稳定的结论可得:稳定期与剧烈开采期NDVI波动较大是由于草本植被变化剧烈引起的,说明矿区植被变化的主要原因为草本植被的变化;在植被修复期初期高频NDVI值呈水平变化特征,说明此时矿区内草本类植被数量稀少,原因为剧烈稀土开采将草本植被破坏严重,几乎不存在草本类植被,在经过人工植被修复后,矿区植被得到一定的补充,草本植被开始增长,最后呈现波动特征。总体表明:矿区内对NDVI影响最大的为草本类植被,基于此情况,矿区植被恢复措施的植被配置模式可针对草本类植被进行选择,并针对草本植被进行维养。

四、讨论

本章基于LandsatTM/OLI 数据和研究区典型地物波谱数据,计算了研究区1989-2015年的NDVI值以及典型地物的NDVI值,总体分为三个部分:

(1)整体矿区NDVI时程变化和双茶亭矿区NDVI时程变化分析。整体矿区NDVI呈下降趋势,可分为三个阶段,即稳定期、剧烈开采期和植被修复期,其中整体矿区NDVI在稳定期呈波动变化,在剧烈开采期和植被修复期呈显著下降趋势,双茶亭矿区NDVI在稳定期和剧烈开采期呈显著下降趋势,在植被修复期呈显著上升趋势。油茶配置模式尾砂修复治理效果较差,桉树配置模式效果较好。

(2)典型地物NDVI值的分析。典型地物可分为三种稀土尾砂、草本植被和乔木类植被,其中尾砂和乔木类植被年内NDVI值变化幅度较小,草本植被年内NDVI值变化较大。

(3)对NDVI时程变化曲线进行小波分解。结果分解为低频信号和高频信号,且分别表示NDVI的年际变化和年内变化特征。通过分析整体矿区NDVI高频信号发现了对研究区NDVI产生影响的主要因素为草本类植被,此结论可指导后期继续进行植被修复,针对草本植被,选择更良好的植被配置模式

推荐

地物光谱仪iSpecField-HH/NIR/WNIR

地物光谱仪是莱森光学专门用于野外遥感测量、土壤环境、矿物地质勘探等领域的最新明星产品,独有的光路设计,噪声校准技术、可以实时自动校准暗电流,采用了固定全息光栅一次性分光,测试速度快,最短积分时间最短可达20μs,操作灵活、便携方便、光谱测试速度快、光谱数据准确,广泛应用于遥感测量、农作物监测、森林研究、海洋学研究和矿物勘察等各领域。

wKgZPGdX6TqAQA-jAAL-2GIR5Jc752.png


审核编辑 黄宇

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • 地物光谱仪
    +关注

    关注

    0

    文章

    75

    浏览量

    3241
收藏 人收藏

    相关推荐

    电磁波波谱图简介

    电磁波波谱图简介    电磁波波谱图   [此贴子已经被作者于2010-3-13 15:36:29编辑过]
    发表于 03-13 15:36

    上海本诺银胶;江南五荧光粉

    `深圳市恒晶电子材料有限公司是日本本诺电子材料、中国五集团下属子公司江南稀土材料研发的LED荧光粉在华南片区的指定代理商,主营本诺银胶/绝缘胶、江南YAG荧光粉。恒晶与原厂保持非常紧密的合作关系
    发表于 08-02 14:29

    回收蚂蚁机,回收比特大陆机,收购比特币

    回收蚂蚁机,回收比特大陆机,收购比特币机全国上门回收蚂蚁机比特大陆机包括(蚂蚁机T9
    发表于 06-01 10:32

    gis专业词汇汇总

    Band 光或热反射的电磁波谱( 如紫外线、兰、绿、红、近外线、远红外线、热、雷达等) 范围的多波谱图像表示数据的一层 。也指其他的用户指定的由原始图像运算得到的数值。多波谱彩色图想显示有三个波段
    发表于 05-24 07:53

    地物光谱仪功能及所需的测量条件是什么

    地物光谱仪的电磁波响应特性随电磁波长改变而变化的规律,称为地物波谱地物光谱仪是电磁辐射与地物相互作用的结果。不同的物质反射、透射、吸收、散
    发表于 12-13 15:55 1259次阅读

    由时事谈稀土与音圈电机模组的关系

    12月22日,五稀土对外发布公告称,公司接到中国五集团有限公司通知
    发表于 12-27 14:12 429次阅读

    关于地物光谱仪的设计优点

    的重要仪器,也是提取地物光谱特征信息的重要工具。 地物光谱仪可以实时测量原始数据,反射、透射辐射和辐照度波谱曲线。地物光谱仪全部使用固定全息光栅和全线阵列光学探测器。这带来了高灵敏度和
    发表于 09-23 10:28 386次阅读
    关于<b class='flag-5'>地物</b>光谱仪的设计优点

    高光谱遥感在植被监测中的研究综述

    高光谱遥感已成为地表植被地学过程对地观测的强有力的工具,其特点是在特定光谱域以高光谱分辨率同时获取连续的地物光谱图像,使得遥感应用着重于在光谱维上进行空间信息展开,定量分析地球表层生物物理化学
    的头像 发表于 10-24 11:18 1460次阅读
    高光谱遥感在<b class='flag-5'>植被</b>监测中的研究综述

    金属矿产遥感找原理与方法-莱森光学

    1、金属矿产遥感找原理 遥感的理论基础是电磁波辐射,地物的光谱辐射特性则是遥感技术赖以鉴别和区分的主要基础。航空航天遥感传感器接收的是地表各种地物的反射光谱信息,这些地物反射光谱实际
    的头像 发表于 12-12 14:02 2115次阅读
    金属矿产遥感找<b class='flag-5'>矿</b>原理与方法-莱森光学

    地物光谱仪应用于监测植被的变化

    植被的变化是地理环境变化的重要指标,可以反映地表土壤、气候、植物群落等环境因子的变化情况。为了更准确地监测植被变化,地物光谱仪可以用来监测植被变化。
    发表于 03-10 10:21 566次阅读
    <b class='flag-5'>地物</b>光谱仪应用于监测<b class='flag-5'>植被</b>的变化

    地面成像和非成像地物光谱仪在不同水分环境下湿地植被光谱特征对比研究

    地面成像和非成像地物光谱仪在不同水分环境下湿地植被光谱特征对比研究——以东洞庭湖湿地区域苔草植被光谱特征为例   引言 在较为复杂的湿地生态系统中,高光谱遥感凭借其成百上千波段数量、极高的光谱分辨率
    的头像 发表于 04-18 18:18 1256次阅读
    地面成像和非成像<b class='flag-5'>地物</b>光谱仪在不同水分环境下湿地<b class='flag-5'>植被</b>光谱特征对比研究

    手持式地物光谱仪对水体叶绿素的光谱特性测试研究

    为什么用地物光谱仪对地物进行测试? 对于地物的测试,传统的方法有化学法等。随着光谱学的发展,现在可以运用光谱法去对地物进行测试,获取地物自身
    的头像 发表于 07-03 15:51 1096次阅读
    手持式<b class='flag-5'>地物</b>光谱仪对水体叶绿素的光谱特性测试研究

    地物光谱匹配模型研究

    都与对波谱特性的分析密切相关。光谱匹配模型通过对地物光谱与参考光谱的匹配或地物光谱与数据库的比较,求算它们之间的相似性或差异性,突出特征谱段,有效地提取光谱维信息,以便对地物特性进行详
    的头像 发表于 07-07 14:38 526次阅读
    <b class='flag-5'>地物</b>光谱匹配模型研究

    地物光谱仪常见植被指数总结

    植被指数:利用卫星不同波段探测数据组合而成的,能反映植物生长状况的指数。植物叶面在可见光红光波段有很强的吸收特性,在近红外波段有很强的反射特性,这是植被遥感监测的物理基础,通过这两个波段测值的不同组
    的头像 发表于 07-19 16:55 2136次阅读
    <b class='flag-5'>地物</b>光谱仪常见<b class='flag-5'>植被</b>指数总结

    地物光谱仪:地物的反射光谱与地物波谱特性

    遥感波谱范围不断拓宽,空间、波谱和时间分辨率不断提高,遥感影像处理,地物识别和信息提取技术不断完善,遥感数据获取由多光谱,高光谱至超高光谱,信息挖掘由“粗糙”到“精细”,遥感地物识别由
    的头像 发表于 08-22 15:49 918次阅读