CN112116242B - 一种结合多种遥感指标的裸土变化识别方法 - Google Patents
一种结合多种遥感指标的裸土变化识别方法 Download PDFInfo
- Publication number
- CN112116242B CN112116242B CN202010977898.3A CN202010977898A CN112116242B CN 112116242 B CN112116242 B CN 112116242B CN 202010977898 A CN202010977898 A CN 202010977898A CN 112116242 B CN112116242 B CN 112116242B
- Authority
- CN
- China
- Prior art keywords
- bare soil
- index
- trend
- indexes
- vegetation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002689 soil Substances 0.000 title claims abstract description 162
- 230000008859 change Effects 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000011160 research Methods 0.000 claims abstract description 46
- 238000000926 separation method Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 230000007423 decrease Effects 0.000 claims description 7
- 230000001174 ascending effect Effects 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 230000000630 rising effect Effects 0.000 claims description 6
- 238000002310 reflectometry Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 10
- 238000011156 evaluation Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 238000004162 soil erosion Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Multimedia (AREA)
- Operations Research (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种结合多种遥感指标的裸土变化识别方法，一种结合多种遥感指标的裸土变化识别方法，包括以下步骤：步骤S1：建立植被近红外反射指数时序数据集；步骤S2：建立裸土指数时序数据集；步骤S3：建立特征指标时序数据集；步骤S4：特征指标评估和选取；步骤S5：计算特征指标的变化趋势；步骤S6：建立裸土变化识别方案；步骤S7：对研究区裸土变化进行识别，获得研究区裸土变化空间分布图：根据步骤S6所建立的裸土变化识别方案逐像元判断裸土是否发生变化，最终生成研究区裸土变化识别结果图。本发明适用于大范围裸土变化监测，时空连续性强，不需要事先提取裸土，且能够快速获取长时间序列裸土变化空间分布图。
Description
技术领域
本发明涉及时序遥感变化监测技术领域，特别是一种结合多种遥感指标的裸土变化识别方法。
背景技术
水土流失是当今世界面临的一个严峻的生态环境问题，由其诱发的各种灾害也给人们的生产生活带来了威胁和损害，而其中裸土是水土流失区的主要地表景观。对此，快速并准确地定位裸土区域，进行地表裸土的时空变化分析，对于水土流失的治理工作尤其重要。近年来，遥感技术已成为水土流失治理和评价不可或缺的高新技术。对此，众多学者采用遥感技术对裸土变化进行研究，其中主要方法包括三类：(1)通过人工目视解译方法进行对研究区遥感影像进行裸土识别；(2)通过采取数据样本，采用监督分类方法对研究区进行裸土识别，结合多个时相分类结果对比分析裸土时空变化特征。(3)为了解决建筑用地和裸土混杂的情况，有学者提出一种双重指数法，该方法首先利用裸土指数并阈值法提取裸土，接受利用归一化建筑指数来去除裸土中混杂的建设用地。
当前众多学者对裸土变化识别方法进行研究，但依然存在以下不足：1，基于人工目视解译的方法存在工作量大，获取大范围裸土信息比较困难等问题。2，基于监督分类的方法，在训练样区选取上难度大，而非监督分类方法获取的分类结果精度低。3，由于存在类内异质性，该方法不适于大范围裸土识别；4，进行裸土时空演变分析时，采用分类后比较方法会带来较大误差。
发明内容
有鉴于此，本发明的目的是提供一种结合多种遥感指标的裸土变化识别方法，结合植被丰度和裸土指数变化趋势，最终确定裸土增加或减少的区域；简单易操作且结果准确度高，不需要设置阈值，能很好地应用于大范围长时间序列裸土变化区域快速监测。
本发明采用以下方案实现：一种结合多种遥感指标的裸土变化识别方法，包括以下步骤：
步骤S1：建立植被近红外反射指数时序数据集；
步骤S2：建立裸土指数时序数据集；
步骤S3：建立特征指标时序数据集；
步骤S4：特征指标评估和选取；
步骤S5：计算特征指标的变化趋势；
步骤S6：建立裸土变化识别方案；
步骤S7：进行裸土变化识别，获得研究区裸土变化空间分布图：根据步骤S6所建立的裸土变化识别方案逐像元判断裸土是否发生变化，最终生成研究区裸土变化识别结果图。
进一步地，步骤S1的具体内容为：在研究区范围内，对采集的图像逐像元逐年建立研究时段内植被近红外反射指数的8天最大化合成多年时序数据集；利用WhittakerSmoother平滑方法，对原始8天最大化合成植被近红外反射指数多年时序数据平滑处理，从而逐像元获得多年平滑的逐日植被近红外反射指数时序数据集。
进一步地，所述步骤S2的具体内容为：
裸土指数(简称RIBS)，其计算公式为：
其中，NNDSI指归一化差分积雪指数，NBrightness是通过缨帽变换获取的亮度指标归一化得到，其计算公式分别为：
逐像元逐年逐期计算裸土指数，生成研究时段内研究区裸土指数时序数据集。
进一步地，所述步骤S3的具体内容为：对研究区，逐年逐像元提取包括植被近红外反射指数和裸土指数的年最大值、年均值和50th分位值即植被丰度特征指标，并建立研究时段内研究区各特征指标时序数据集。
进一步地，所述步骤S4的具体内容为：选取若干不同地物类型的稳定点位，具体地物类型包括：植被、不透水面和裸土，分别评估基于植被近红外反射率指数和裸土指数的年最大值、年均值和50th分位值即植被丰度特征指标在不同地物类型中的分离度，分离度计算公式：
其中，μi、μs分别代表地物类型A和地物类型B的均值，σi、σs分别代表地物类型A和地物类型B的标准差，SDI用于比较两种不同地物类型的分离度；SDI越大，分离度越大；如果值小于1，表明分离度差；最终选择拥有高分离度的特征指标，经过计算选取植被丰度和裸土指数年均值作为裸土变化识别的指标。
进一步地，所述步骤S5的具体内容为：利用一元线性拟合的方法，分别基于研究时段内特征指标时序数据集，逐像元依次计算研究时段内特征指标的变化趋势K；计算公式如下：
y＝Kx+b
当K>0时，表示该时序曲线呈上升趋势；当K<0时，表示该时序曲线具有下降趋势；基于Mann-Kendall方法，进一步判断该时序曲线的变化趋势是否显著；依据显著性检验的结果，将变化趋势分为三种情况：显著的上升趋势、没有趋势和显著的下降趋势；
基于Mann-Kendall方法，进一步判断该时序曲线的变化趋势是否显著；假设有n个样本量x1……xn的时间序列，对于所有k,j≤n,且k≠j，xk和xj的分布是不同的，计算检验统计量s，公式如下:
其中，S为正态分布，均值为0，方差var(s)＝n(n-1)(2n+5)/18；当n>10时，标准的正态统计变量通过下式计算:
对于统计值Z来说，大于0时，表示增加趋势；小于0时，表示减少趋势；Z的绝对值在大于1.96时，表示通过了置信度95％的显著性检验；当Z>1.96时，表示该时序曲线呈显著上升趋势；当z<-1.96时，表示该时序曲线具有显著下降趋势；当Z在正负1.96之间时为无趋势。
进一步地，步骤S6中所述的建立裸土变化识别方案的具体内容为：根据特征指标的变化趋势进行组合，将裸土指数和植被丰度均呈显著下降趋势的像元，判断为裸土增加；将裸土指数和植被丰度均呈显著上升趋势的像元，判断为裸土减少。
与现有技术相比，本发明具有以下有益效果：
(1)本发明适用于大范围裸土变化监测，时空连续性强，不需要事先提取裸土，且能够快速获取长时间序列裸土变化空间分布图。
(2)本发明充分考虑到不同地物类型在不同遥感指标中呈现的差异性，即植被丰度和裸土指数均最低，而其他地物通常在植被丰度较低的情况下，具有较高的裸土指数，进而依据植被丰度、裸土指数判断是否发生裸土变化。
(3)本发明充分利用多指标多年变化趋势，判断是否发生裸土增加或减少情况，相较于现有技术，能够提高识别结果的精度。
附图说明
图1为本发明实施例的流程图。
图2为本发明实施例的特征指标评估图；其中图2(a)为不同地物类型点位的植被指数年最大值箱线图；图2(b)为不同地物类型点位的植被指数年均值箱线图，图2(c)为不同地物类型点位的植被丰度(植被指数计算的)箱线图，图2(d)为不同地物类型点位的裸土指数年最大值箱线图，图2(e)为不同地物类型点位的裸土指数年均值箱线图，图2(f)为不同地物类型点位的植被丰度(裸土指数计算的)箱线图。
图3为本发明实施例的裸土减少点位各遥感指标时序信号图,其中，图3(a)为裸土减少点位裸土指数原始信号图，图3(b)为裸土指数年均值年际趋势图，图3(c)为植被指数原始信号图，图3(d)为植被丰度指标年际趋势图。
图4为本发明实施例的裸土增加点位各遥感指标时序信号图,其中，图4(a)为裸土增加点位裸土指数原始信号图，图4(b)为裸土指数年均值年际趋势图，图4(c)为植被指数原始信号图，图4(d)为植被丰度指标年际趋势图。
图5为本发明实施例的特征指标趋势组合图。
图6为本发明实施例的裸土变化识别流程图。
图7为本发明实施例的研究区裸土变化空间分布图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示，本实施例提供一种结合多种遥感指标的裸土变化识别方法，
首先基于植被近红外反射指数和裸土指数构建了多种特征指标，并进一步评估和选取最能体现不同地物类型差异的特征指标，根据不同地物类型之间的分离度结果表明，基于植被近红外反射指数构建的植被丰度和裸土指数的年均值最能体现不同地物类型的差异。通过判断植被丰度、裸土指数是否存在显著的变化趋势，最终依据植被丰度、裸土指数的变化趋势，建立一种结合多种遥感指标的裸土变化识别方法。
包括以下步骤：
步骤S1：建立植被近红外反射指数时序数据集；
步骤S2：建立裸土指数时序数据集；
步骤S3：建立特征指标时序数据集；
步骤S4：特征指标评估和选取；
步骤S5：计算特征指标的变化趋势；
步骤S6：建立裸土变化识别方案；
步骤S7：对研究区裸土变化进行识别，获得研究区裸土变化空间
分布图：根据步骤S6建立的裸土变化识别方案逐像元判断裸土是否发生变化，最终生成研究区裸土变化识别结果图。
在本实施例中，步骤S1的具体内容为：在研究区范围内，对采集的图像逐像元逐年建立研究时段内植被近红外反射指数的8天最大化合成多年时序数据集；利用WhittakerSmoother平滑方法，对原始8天最大化合成植被近红外反射指数多年时序数据平滑处理，从而逐像元获得多年平滑的逐日植被近红外反射指数时序数据集。
在本实施例中，在遥感领域，大多时序变化监测方法一般时间跨度都会在10年以上，但这个也没有具体的规定，时间跨度越长，方法监测结果精度会更高。在本实施例中，研究时段在10年以上。
在本实施例中，所述步骤S2的具体内容为：
裸土指数(简称RIBS)，其计算公式为：
其中，NNDSI指归一化差分积雪指数，NBrightness是通过缨帽变换获取的亮度指标归一化得到，其计算公式分别为：
逐像元逐年逐期计算裸土指数，生成研究时段内研究区裸土指数时序数据集。
在本实施例中，所述步骤S3的具体内容为：对研究区，逐年逐像元提取包括植被近红外反射指数和裸土指数的年最大值、年均值和50th分位值即植被丰度特征指标，并建立研究时段内研究区各特征指标时序数据集。
在本实施例中，所述步骤S4的具体内容为：选取若干不同地物类型的稳定点位，具体地物类型包括：植被、不透水面和裸土。探索植被、不透水面和裸土在裸土指数和植被指数不同分位值上表现出的差异性；分别评估基于植被近红外反射率指数和裸土指数的年最大值、年均值和50th分位值即植被丰度特征指标在不同地物类型中的分离度如图2，并选择拥有最佳分离的特征指标；根据评估结果表明：基于植被近红外反射指数构建植被丰度指标和裸土指数年均值能更好区分三种地物类型，因此，最终选取植被丰度和裸土指数年均值作为裸土变化识别的指标。
分离度计算公式：
其中μi、μs分别代表地物类型A和地物类型B的均值，σi、σs分别代表地物类型A和地物类型B的标准差，SDI用于比较两种不同地物类型的分离度；SDI越大，分离度越大；如果值小于1，表明分离度差；
根据对比结果表明，裸土指数年均值计算的不同地物类型两两之间的SDI数值要高于年最大值和50th分位值；而植被指数50th分位值计算的不同地物类型两两之间的SDI数值要高于年最大值和年均值。因此。最终选取植被丰度和裸土指数年均值作为裸土变化识别的指标。
在本实施例中，所述步骤S5的具体内容为：利用一元线性拟合的方法，分别基于研究时段内特征指标时序数据集，逐像元依次计算研究时段内特征指标的变化趋势K；计算公式如下：
y＝Kx+b
当K>0时，表示该时序曲线呈上升趋势；当K<0时，表示该时序曲线具有下降趋势；基于Mann-Kendall方法，进一步判断该时序曲线的变化趋势是否显著；依据显著性检验的结果，将变化趋势分为三种情况：显著的上升趋势、没有趋势和显著的下降趋势；
基于Mann-Kendall方法，进一步判断该时序曲线的变化趋势是否显著；假设有n个样本量x1……xn的时间序列，对于所有k,j≤n,且k≠j，xk和xj的分布是不同的，计算检验统计量s，公式如下:
其中，S为正态分布，均值为0，方差var(s)＝n(n-1)(2n+5)/18；当n>10时，标准的正态统计变量通过下式计算:
对于统计值Z来说，大于0时，表示增加趋势；小于0时，表示减少趋势；Z的绝对值在大于1.96时，表示通过了置信度95％的显著性检验；在本实施例选择了Z的绝对值大于1.96，即置信度为95％。当Z>1.96时，表示该时序曲线呈显著上升趋势；当z<-1.96时，表示该时序曲线具有显著下降趋势；当Z在正负1.96之间时为无趋势。
在本实施例中，步骤S6中所述的建立裸土变化识别方案的具体内容为：根据特征指标的变化趋势进行组合，将裸土指数和植被丰度均呈显著下降趋势的像元，判断为裸土增加；将裸土指数和植被丰度均呈显著上升趋势的像元，判断为裸土减少。
对于植被、不透水面、裸地和水体四种不同地物类型而言，裸土的裸土指数最低，其次为植被，再次为不透水面，而水体的裸土指数最高；植被丰度也存在明显差异。具体来说，植被的植被丰度最高，其次为不透水面、裸地，水体的植被丰度最低；裸土变化识别方案主要依据特征指标的数值在不同地物类型中存在明显的差异，因此，假设当地物类型发生变化时比如裸土变成植被，特征指标在年际水平上存在明显上升或者下降的变化趋势；因此，选取若干已知变化类型的点位，验证上述假设的合理性，并随机选取两个点位进行说明，第一个为裸土减少点位，第二个为裸土增加点位。对于第一个点位，当裸土减少时，基于植被近红外反射指数构建的植被丰度指标和裸土指数年均值呈现显著上升趋势如图3；对于第二个点位，当裸土增加时，基于植被近红外反射指数构建的植被丰度指标和裸土指数年均值呈现显著下降趋势如图4。依据实验结果，假设成立。
由于水体不稳定，易受包括降雨或冰川融化因素影响，所以不考虑水体，因此需对其进行掩膜处理；根据特征指标趋势组合的结果(图5)，采用研究时段开始年份和结束年份均值大于-0.05掩模掉水体发生变化的区域和变成水体的区域；最后得到裸土变化识别的方案：将裸土指数和植被丰度均呈显著下降趋势的像元，判断为裸土增加；将裸土指数和植被丰度均呈显著上升趋势的像元，判断为裸土减少。裸土变化识别方法技术流程图见图6。
以上上述流程，逐像元判断裸土是否发生变化，最终生成研究区裸土变化识别结果图。依据本实施例中提供的方法，以我国云南省为例，所获得的研究区2001-2018年裸土变化结果分布图见图7。
在本实施例中，在特征指标评估和选取时本实施例没有随意使用某一个特征指标，而是经过评估后才挑选的。特征指标评估和选取的步骤如下：
(a)选取若干不同地物类型的稳定点位，具体地物类型包括：植被、不透水面和裸土。
(b)根据不同地物类型点位分布位置提取不同特征指标的数值。
(c)分别评估基于植被近红外反射率指数和裸土指数的年最大值、年均值、50th分位值(植被丰度)等特征指标在不同地物类型中的分离度(图2)，并选择拥有最佳分离的特征指标。
在本实施例中，在建立裸土变化识别方法时。主要步骤如下：(a)根据特征指标评估时所表现出来的特点，特征指标在不同地物类型中存在明显的差异，因此，提出假设，当地物类型发生变化时比如裸土变成植被，特征指标在年际水平时序信号上存在明显上升或者下降的变化趋势。
(b)选取若干已知变化类型的点位，验证上述假设的合理性。
(c)在假设成立的基础上，继续统计裸土增加和裸土减少点位所满足的特征指标趋势组合的情况，并把90％以上点位均满足的规则作为本实施例裸土变化识别的方法。
依据本实施例提供的方法，以我国云南省为例，所获得的研究区2001-2018年裸土变化结果分布图见图7。结果显示，2001-2018裸土增加的区域相对较为集中，主要分布在大理、曲靖和九乡市；裸土减少的区域则比较分散，分布在云南省各个地方，没有相对较集中的区域。为了验证该结果的精度，随机挑选几处裸土发生变化的区域，与谷歌地图高分影像进行目视解译对比验证。通过谷歌历史影像与挑选的出来的区域进行
对比发现，该区域确实发生明显变化，与变化识别的结果一致。以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。
Claims (3)
1.一种结合多种遥感指标的裸土变化识别方法，其特征在于：包括以下步骤：
步骤S1：建立植被近红外反射指数时序数据集；
步骤S2：建立裸土指数时序数据集；
步骤S3：建立特征指标时序数据集；
步骤S4：特征指标评估和选取；
步骤S5：计算特征指标的变化趋势；
步骤S6：建立裸土变化识别方案；
步骤S7：对研究区裸土变化进行识别，获得研究区裸土变化空间分布图：根据步骤S6所建立的裸土变化识别方案逐像元判断裸土是否发生变化，最终生成研究区裸土变化识别结果图；
所述步骤S4的具体内容为：选取若干不同地物类型的稳定点位，具体地物类型包括：植被、不透水面和裸土；分别评估基于植被近红外反射率指数和裸土指数的年最大值、年均值和50th分位值即植被丰度特征指标在不同地物类型中的分离度，分离度计算公式：
其中，μi、μs分别代表地物类型A和地物类型B的均值，σi、σs分别代表地物类型A和地物类型B的标准差，SDI用于比较两种不同地物类型的分离度；SDI越大，分离度越大；如果值小于1，表明分离度差；最终选择拥有最高分离度的特征指标，经过计算选取植被丰度和裸土指数年均值作为裸土变化识别的指标；
步骤S1的具体内容为：在研究区范围内，对研究时间范围内采集的图像逐像元逐年建立研究时段内植被近红外反射指数的8天最大化合成多年时序数据集；利用WhittakerSmoother平滑方法，对原始8天最大化合成植被近红外反射指数多年时序数据平滑处理，从而逐像元获得多年平滑的逐日植被近红外反射指数时序数据集；
所述步骤S2的具体内容为：
裸土指数，其计算公式为：
其中，NNDSI指归一化差分积雪指数，NBrightness是通过缨帽变换获取的亮度指标归一化得到，其计算公式分别为：
逐像元逐年逐期计算裸土指数，生成研究时段内研究区裸土指数时序数据集；
所述步骤S3的具体内容为：对研究区，逐年逐像元提取包括植被近红外反射指数和裸土指数的年最大值、年均值和50th分位值即植被丰度特征指标，并建立研究时段内研究区各特征指标时序数据集。
2.根据权利要求1所述的一种结合多种遥感指标的裸土变化识别方法，其特征在于：所述步骤S5的具体内容为：利用一元线性拟合的方法，分别基于研究时段内特征指标时序数据集，逐像元依次计算研究时段内特征指标的变化趋势K；计算公式如下：
y＝Kx+b
当K>0时，表示时序曲线呈上升趋势；当K<0时，表示时序曲线具有下降趋势；
依据Mann-Kendall显著性检验的结果，将变化趋势分为三种情况：显著的上升趋势、没有趋势和显著的下降趋势；
基于Mann-Kendall方法，进一步判断该时序曲线的变化趋势是否显著；假设有n个样本量x1……xn的时间序列，对于所有k,j≤n,且k≠j，xk和xj的分布是不同的，计算检验统计量s，公式如下:
其中，S为正态分布，均值为0，方差var(s)＝n(n-1)(2n+5)/18；当n>10时，标准的正态统计变量通过下式计算:
对于统计值Z来说，大于0时，表示增加趋势；小于0时，表示减少趋势；Z的绝对值在大于1.96时，表示通过了置信度95％的显著性检验；当Z>1.96时，表示该时序曲线呈显著上升趋势；当z<-1.96时，表示该时序曲线具有显著下降趋势。
3.根据权利要求1所述的一种结合多种遥感指标的裸土变化识别方法，其特征在于：步骤S6中所述的建立裸土变化识别方案的具体内容为：根据特征指标的变化趋势进行组合，将裸土指数和植被丰度均呈显著下降趋势的像元，判断为裸土增加；将裸土指数和植被丰度均呈显著上升趋势的像元，判断为裸土减少。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010977898.3A CN112116242B (zh) | 2020-09-17 | 2020-09-17 | 一种结合多种遥感指标的裸土变化识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010977898.3A CN112116242B (zh) | 2020-09-17 | 2020-09-17 | 一种结合多种遥感指标的裸土变化识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112116242A CN112116242A (zh) | 2020-12-22 |
| CN112116242B true CN112116242B (zh) | 2022-08-16 |
Family
ID=73803208
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010977898.3A Active CN112116242B (zh) | 2020-09-17 | 2020-09-17 | 一种结合多种遥感指标的裸土变化识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112116242B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113052922B (zh) * | 2021-03-26 | 2023-03-14 | 重庆紫光华山智安科技有限公司 | 裸土识别方法、系统、设备及介质 |
| CN113486713A (zh) * | 2021-06-02 | 2021-10-08 | 广州林猫自然科技有限公司 | 一种红外相机动物图像识别训练集最佳组合的构建方法 |
| CN113378754B (zh) * | 2021-06-24 | 2023-06-20 | 中国计量大学 | 一种工地裸土监测方法 |
| CN116052017B (zh) * | 2023-01-17 | 2023-11-10 | 二十一世纪空间技术应用股份有限公司 | 基于双重指数法和自适应策略绿网苫盖信息遥感监测方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005005926A1 (en) * | 2003-07-09 | 2005-01-20 | Fabrizio Ferrucci | Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation |
| CN106198434A (zh) * | 2016-07-26 | 2016-12-07 | 新疆农垦科学院 | 一种北疆灰漠土土壤有机质含量的遥感监测方法 |
| CN107463775A (zh) * | 2017-07-24 | 2017-12-12 | 福州大学 | 基于多遥感指数变化趋势的植被流失去向识别方法 |
| CN109840516A (zh) * | 2019-03-06 | 2019-06-04 | 福州大学 | 一种基于时序遥感影像的水体变化自动识别方法 |
| CN110378576A (zh) * | 2019-07-01 | 2019-10-25 | 中国环境科学研究院 | 城市化植被效应有效距离的定量化检测方法 |
| CN110852381A (zh) * | 2019-11-11 | 2020-02-28 | 四川航天神坤科技有限公司 | 一种林火火烧迹地的提取方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10234439B2 (en) * | 2012-11-07 | 2019-03-19 | Airscout Inc. | Methods and systems for analyzing a field |
| US9390331B2 (en) * | 2014-04-15 | 2016-07-12 | Open Range Consulting | System and method for assessing riparian habitats |
| CN108388832A (zh) * | 2018-01-11 | 2018-08-10 | 福州大学 | 一种基于多时序指标变化趋势的退耕还林自动识别方法 |
| CN110580474B (zh) * | 2019-09-30 | 2020-08-14 | 生态环境部卫星环境应用中心 | 一种基于多源数据的农田重金属高风险区遥感快速识别方法 |
-
2020
- 2020-09-17 CN CN202010977898.3A patent/CN112116242B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005005926A1 (en) * | 2003-07-09 | 2005-01-20 | Fabrizio Ferrucci | Method and apparatus for automatically detecting and mapping, particularly for burnt areas without vegetation |
| CN106198434A (zh) * | 2016-07-26 | 2016-12-07 | 新疆农垦科学院 | 一种北疆灰漠土土壤有机质含量的遥感监测方法 |
| CN107463775A (zh) * | 2017-07-24 | 2017-12-12 | 福州大学 | 基于多遥感指数变化趋势的植被流失去向识别方法 |
| CN109840516A (zh) * | 2019-03-06 | 2019-06-04 | 福州大学 | 一种基于时序遥感影像的水体变化自动识别方法 |
| CN110378576A (zh) * | 2019-07-01 | 2019-10-25 | 中国环境科学研究院 | 城市化植被效应有效距离的定量化检测方法 |
| CN110852381A (zh) * | 2019-11-11 | 2020-02-28 | 四川航天神坤科技有限公司 | 一种林火火烧迹地的提取方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 福建省永泰县地表裸土动态变化分析1概述;林冬凤等;《海峡科学》;20170228(第02期);第3-6页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112116242A (zh) | 2020-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112116242B (zh) | 一种结合多种遥感指标的裸土变化识别方法 | |
| CN106023185B (zh) | 一种输电设备故障诊断方法 | |
| CN102662949B (zh) | 一种基于多特征融合的指定目标检索方法及系统 | |
| CN108052980B (zh) | 基于图像的空气质量等级检测方法 | |
| Etzelmüller et al. | Automatic regional classification of topography in Norway. | |
| CN111860351B (zh) | 一种基于行列自注意力全卷积神经网络的遥感图像鱼塘提取方法 | |
| CN101789005A (zh) | 一种基于感兴趣区域的图像检索方法 | |
| CN103714349A (zh) | 一种基于颜色和纹理特征的图像识别方法 | |
| CN114331031B (zh) | 一种建筑传统特征识别评价方法和系统 | |
| CN106844739A (zh) | 一种基于神经网络协同训练的遥感图像变化信息检索方法 | |
| CN108549961B (zh) | 一种基于cmip5预估海浪有效波高的方法 | |
| CN111461067B (zh) | 基于先验知识映射及修正的零样本遥感影像场景识别方法 | |
| CN113378971A (zh) | 近红外光谱的分类模型训练方法、系统及分类方法、系统 | |
| CN109447153A (zh) | 用于非均衡数据分类的散度-激励自编码器及其分类方法 | |
| CN111667473A (zh) | 基于改进Canny算法的绝缘子憎水性等级判断方法 | |
| CN114776304A (zh) | 一种深海矿产异常区识别的方法、装置、设备及介质 | |
| CN101876993B (zh) | 一种地基数字云图的纹理特征抽取和检索方法 | |
| CN113671599A (zh) | 一种基于全球气候模式的登陆气旋识别方法 | |
| CN107608938B (zh) | 基于增强回归树算法的面向二值分类的因子筛选方法 | |
| CN117078077A (zh) | 一种高速公路路域的生态脆弱性评价方法 | |
| CN116519710A (zh) | 一种用于复合绝缘子表面污秽状态的检测方法及系统 | |
| CN115497006A (zh) | 基于动态混合策略的城区遥感图像变化深度监测方法及系统 | |
| CN110251076A (zh) | 融合视觉注意力基于对比度的显著性检测方法及装置 | |
| CN112256766A (zh) | 用于能源采集终端的用电行为分析方法 | |
| CN111898579A (zh) | 基于极限梯度提升的高分遥感影像无偏半监督分类模型 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |