一种基于时间序列遥感影像的茶园打顶时间检测方法技术

本发明专利技术公开了一种基于时间序列遥感影像的茶园打顶时间检测方法，包括：获取茶园区域内Landsat

【技术实现步骤摘要】
一种基于时间序列遥感影像的茶园打顶时间检测方法


[0001]本专利技术涉及遥感目标识别
，特别是涉及一种基于时间序列遥感影像的茶园打顶时间检测方法。

技术介绍

[0002]茶树生长过程中，顶芽生长旺盛，侧芽生长较弱或缓慢，呈现明显的顶端优势。茶园打顶是调节顶端优势，促进树冠培养最重要的技术措施，对于茶园优质高产和持续健康发展具有十分重要的作用。茶园打顶会造成短期内显著抑制茶树的生殖生长，茶树茶叶脱落，导致独特的“人工物候”现象。传统茶园打顶时间确认方法依靠人工经验判断，费时费力。
[0003]卫星遥感对地观测具有重访性特点，时间序列遥感影像数据能较为真实地反映地表在一个长时间范围内的动态变化情况。基于时间序列遥感影像可以检测到土地覆盖变化的区域和程度、变化类型的分布、变化的类型。因此，加强卫星遥感技术在茶园打顶时间检测上的应用具有重要意义。

技术实现思路

[0004]本专利技术的目的是提供一种基于时间序列遥感影像的茶园打顶时间检测方法，可有效利用茶园在时间序列遥感影像上的时序特征，快速检测茶园的打顶时间，弥补了传统本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于时间序列遥感影像的茶园打顶时间检测方法，其特征在于，包括以下步骤：S1，获取茶园研究区域中研究期内所有Landsat
‑
8 OLI与Sentinel
‑
2 MSI卫星传感器的地表反射率数据；S2，根据获取的地表反射率数据构建时间分辨率为5天的时间序列遥感影像集，并对时间序列遥感影像集进行预处理；S3，基于时间序列遥感影像集进行每个像元的NDVI、LSWI、EVI、RGRI指数的计算以及平滑处理；S4，计算每个像元在时间序列遥感影像集上相邻时间的指数变化强度系数及NDVI、LSWI、EVI、RGRI指数的变化方向；S5，基于时间序列遥感影像集，获取像元上NDVI、LSWI、EVI指数变化方向为下降同时RGRI指数变化方向为上升且指数变化强度系数最大处的相应时间；S6，步骤S5获得的时间即为茶园打顶时间，完成对茶园打顶时间的检测。2.根据权利要求1所述的基于时间序列遥感影像的茶园打顶时间检测方法，其特征在于，所述步骤S1中，地表反射率数据的来源是：基于谷歌云平台获取Landsat
‑
8 OLI与Sentinel
‑
2 MSI卫星传感器的地表反射率数据，且所述地表反射率数据经过大气校正、辐射校正、地形校正的预处理。3.根据权利要求1所述的基于时间序列遥感影像的茶园打顶时间检测方法，其特征在于，所述步骤S2，根据获取的地表反射率数据构建时间分辨率为5天的时间序列遥感影像集，并对时间序列遥感影像集进行预处理，具体包括：S201，对时间序列遥感影像集进行云量筛选；S202，对时间序列遥感影像集进行去云处理；S203，对时间序列遥感影像集进行裁剪拼接；S204，协调Landsat
‑
8 OLI与Sentinel
‑
2 MSI之间的卫星传感器误差；S205，对时间序列遥感影像集进行影像合成；S206，利用时序线性插值方法填补时间序列遥感影像集由于影像去云造成的空洞；S207，在时间序列遥感影像上增加时间属性，时间属性用DOY表示，DOY定义为影像时间为一年中的第几天。4.根据权利要求1所述的基于时间序列遥感影像的茶园打顶时间检测方法，其特征在于，所述步骤S3，基于时间序列遥感影像集进行每个像元的NDVI、LSWI、EVI、RGRI指数的计算以及平滑处理，具体包括：S301，对归一化植被指数NDVI、陆表水分指数LSWI、增强型植被指数EVI、红绿比值指数RGRI进行计算，计算公式如下：RGRI进行计算，计算公式如下：EVI＝2.5*(B
nir
‑
B
red
)/(B
nir
+6*B
red
‑
7.5*B
blue
...

【专利技术属性】
技术研发人员：徐伟恒，康倩文，王雷光，黄邵东，刘运，洪泽湖，钱瑞，
申请(专利权)人：西南林业大学，
类型：发明
国别省市：