基于双抛物线型NDVI‑Ts特征空间的土壤湿度遥感监测方法技术

本发明专利技术公开了一种基于双抛物线型NDVI‑Ts特征空间的土壤湿度遥感监测方法，涉及土壤检测技术领域，本发明专利技术顾及NDVI＜0.15，发现NDVI‑Ts特征空间呈双抛物线型，提出基于双抛物线型NDVI‑Ts特征空间的TVDI土壤湿度遥感监测方法，通过计算双抛物线型NDVI‑Ts特征空间中干、湿边方程，获得研究区TVDI影像，并根据TVDI值的大小将土壤湿度进行划分，得到待监测地区的地表土壤湿度。该方法在监测地表0‑10cm深土壤湿度状况方面优于传统的三角形NDVI‑Ts特征空间的TVDI土壤湿度监测方法，能更好地反映地表0‑5cm深土壤湿度，为土壤湿度遥感监测提供依据。

【技术实现步骤摘要】

本专利技术涉及土壤检测
，尤其涉及一种基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法。
技术介绍
干旱是人类面临的主要自然灾害之一，而土壤湿度作为陆面水资源形成、转化、消耗过程的基本参数，是反映陆地干旱的重要指标。多时相、多光谱、高光谱遥感数据反映了大面积的地表信息，使得快速、及时、动态监测旱情成为可能。国内外学者基于多源、多光谱、多时相遥感数据，提出了多种土壤湿度遥感监测的方法和模型。总体上，这些方法主要可分为可见光——红外法和微波遥感法。利用可见光和红外波段遥感监测土壤湿度的具体方法有距平植被指数法、条件植被指数法、条件温度指数法、归一化温度指数法、条件植被温度指数法、植被供水指数法、作物缺水指数法、热惯量和表观热惯量法、光谱特征空间法等。温度植被干旱指数(TemperatureVegetationDrynessIndex，TVDI)是一种应用较为广泛的方法。Price和Carlson等发现，当研究区域的植被覆盖和土壤湿度变化范围较大时，根据遥感资料得到的以归一化植被指数(NormalizedDifferenceVegetationIndex，NDVI)为横轴和以地表辐射温度为纵轴的散点图呈三角形特征；Moran等认为NDVI-Ts(地表温度)之间为梯形特征空间；Sandholt等提出了TVDI的计算方法；Patel等基于MODIS产品数据的NDVI-Ts特征空间计算TVDI，并估算了印度亚湿润地区的土壤湿度，结果表明TVDI与实测土壤水分之间存在强负相关关系；Schnur等采用250m分辨率MODIS产品NDVI和增强型植被指数(EnhancedVegetationIndex，EVI)数据评估美国西南区域的土壤水分，发现NDVI与土壤湿度的相关性略高于EVI和土壤湿度的相关性。国内外相关研究大多采用1km或5km分辨率MODIS产品数据来构建NDVI-Ts特征空间，并认为此特征空间呈三角形或梯形，但是，采用NDVI-Ts三角形或梯形特征空间计算TVDI的原理为Tsmax与NDVI呈线性关系，随着NDVI的增加，Tsmax线性减小，且一般认为，当NDVI＜0.15时陆地表面是裸地，线性拟合时不予考虑。
技术实现思路
针对上述缺陷或不足，本专利技术的目的在于提供一种基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法。本专利技术的技术方案为：一种基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法，包括：1)获取研究地区遥感影像数据；2)对遥感影像数据进行预处理，所述预处理包括对遥感影像数据进行辐射校正、大气校正、以及几何校正，得到预处理后的遥感影像图；3)根据遥感影像图，获取研究区归一化植被指数(NDVI)与地表温度(Ts)数据，利用温度植被干旱指数(TVDI)计算公式获取每个遥感影像图中像元的TVDI值，TVDI值表达式为：TVDI=Ts-TsminTsmax-Tsmin---(1)]]>其中，Ts表示地表温度；Tsmin表示相同NDVI值对应的最小地表温度，是NDVI-Ts特征空间中的湿边；Tsmax为研究区相同NDVI值对应的最大地表温度，代表特征空间中的干边；4)构建双抛物线型NDVI-Ts特征空间散点图，并根据双抛物线型NDVI-Ts特征空间获取TVDI中Tsmax、Tsmin算法表达式：Tsmax＝a1×NDVI2+b1×NDVI+c1Tsmin＝a2×NDVI2+b2×NDVI+c2(2)其中，a1、b1、c1、a2、b2、c2为方程拟合系数；5)根据公式(1)和公式(2)得到研究区每个像元的TVDI值；6)根据TVDI值的大小，依据标准将土壤湿度进行划分，得到待监测地区的土壤湿度分布状况。所述步骤6)中根据获得的TVDI，结合前人研究成果，将土壤湿度划分为5类，具体为：当0＜TVDI≤0.2时，则为极湿润；当0.2＜TVDI≤0.4时，则为湿润；当0.4＜TVDI≤0.6时，则为正常；当0.6＜TVDI≤0.8时，则为干旱；当0.8＜TVDI≤I时，则为极干旱。与现有技术比较，本专利技术的有益效果为：本专利技术提供了一种基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法，分析了基于1km、500m、250m分辨率MODIS数据的NDVI-Ts特征空间，发现NDVI-Ts特征空间均具有双抛物线型特征。与基于1km、500m、250m分辨率MODIS数据的NDVI-Ts三角形特征空间进行对比分析，并将得到的TVDI数据与实测土壤湿度进行相关性分析，结果表明：基于双抛物线型NDVI-Ts特征空间的TVDIc在反映地表浅层土壤湿度方面要优于基于三角形NDVI-Ts特征空间的TVDIt；且250m双抛物线型NDVI-Ts特征空间在监测地表0-5cm土壤湿度方面是最好的，而500m双抛物线型NDVI-Ts特征空间在监测10cm深土壤湿度时更有优势。附图说明图1是本专利技术基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法流程图；图2是TVDI在三角形NDVI-Ts特征空间中的定义图；图3是TVDI在双抛物线型NDVI-Ts特征空间中的定义图；图4是本专利技术实施案例中10cm深度土壤湿度与TVDI的关系图。具体实施方式下面结合附图对本专利技术做详细描述。如图1所示，本专利技术提供了一种基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法，包括：1)获取研究地区遥感影像数据；2)对遥感影像数据进行预处理，所述预处理包括对遥感影像数据进行辐射校正、大气校正、以及几何校正，得到预处理后的遥感影像图；3)根据遥感影像图，获取研究区归一化植被指数(NDVI)与地表温度(Ts)数据，利用温度植被干旱指数(TVDI)计算公式获取每个遥感影像图中像元的TVDI值，TVDI值表达式为：TVDI=Ts-TsminTsmax-Tsmin---(1)]]>其中，Ts表示地表温度；Tsmin表示相同NDVI值对应的最小地表温度，是NDVI-Ts特征空间中的湿边；Tsmax为研究区相同NDVI值对应的最大地表温度，代表特征空间中的干边；TVDI值的取值范围是0到1。TVDI值越大，Ts越接近干边，越干旱；反之TVDI值越小，Ts越接近湿边，越湿润。在现有技术中，三角形或梯形NDVI-Ts特征空间计算TVDI原理是：随着NDVI增加，Tsmax呈线性减小趋势；当NDVI＜0.15时，认为陆地表面是裸地，无植被覆盖，在线性拟合Tsmax时不予考虑。本专利技术发现，NDVI-Ts特征空间呈双抛物线型，随着NDVI增加，Tsmax呈非线性减小趋势，如图2所示，与NDVI-Ts三角形或梯形特征空间不一致。但若将双抛物线型特征空间中的干湿边在拟合时不考虑NDVI＜0.15部分，则NDVI-Ts散点图呈三角形。利用上述两种特征空间分别计算TVDI，经对比发现两种特征空间得到的结果相似，但与实测土壤湿度数据(SoilMoisture，SM)进行相关分析时发现，由双抛物线型特征空间计算的TVDI与0-5cm和10cm深SM相关性要优于由三角形特征空间计算的TVDI与SM的相关性。因而，认为NDVI＜0.15部分应被包含在双抛物线型NDVI-Ts特征空间中，不应舍弃。本专利技术将NDVI扩展到0.15以内，构建双抛物线型本文档来自技高网...

【技术保护点】
一种基于双抛物线型NDVI‑Ts特征空间的土壤湿度遥感监测方法，其特征在于，包括：1)获取研究地区遥感影像数据；2)对遥感影像数据进行预处理，所述预处理包括对遥感影像数据进行辐射校正、大气校正、以及几何校正，得到预处理后的遥感影像图；3)根据遥感影像图，获取研究区归一化植被指数(NDVI)与地表温度(Ts)数据，利用温度植被干旱指数(TVDI)计算公式获取每个遥感影像图中像元的TVDI值，TVDI值表达式为：TVDI=Ts-TsminTsmax-Tsmin---(1)]]>其中，Ts表示地表温度；Tsmin表示相同NDVI值对应的最小地表温度，是NDVI‑Ts特征空间中的湿边；Tsmax为研究区相同NDVI值对应的最大地表温度，代表特征空间中的干边；4)构建双抛物线型NDVI‑Ts特征空间散点图，并根据双抛物线型NDVI‑Ts特征空间散点图，获取TVDI中Tsmax、Tsmin算法表达式：Tsmax＝a1×NDVI2+b1×NDVI+c1Tsmin＝a2×NDVI2+b2×NDVI+c2        (2)其中，a1、b1、c1、a2、b2、c2为方程拟合系数；5)根据公式(1)和公式(2)得到研究区每个像元的TVDI值；6)根据TVDI值的大小，依据标准将土壤湿度进行划分，得到待监测地区的土壤湿度分布状况。...

【技术特征摘要】
1.一种基于双抛物线型NDVI-Ts特征空间的土壤湿度遥感监测方法，其特征在于，包括：1)获取研究地区遥感影像数据；2)对遥感影像数据进行预处理，所述预处理包括对遥感影像数据进行辐射校正、大气校正、以及几何校正，得到预处理后的遥感影像图；3)根据遥感影像图，获取研究区归一化植被指数(NDVI)与地表温度(Ts)数据，利用温度植被干旱指数(TVDI)计算公式获取每个遥感影像图中像元的TVDI值，TVDI值表达式为：TVDI=Ts-TsminTsmax-Tsmin---(1)]]>其中，Ts表示地表温度；Tsmin表示相同NDVI值对应的最小地表温度，是NDVI-Ts特征空间中的湿边；Tsmax为研究区相同NDVI值对应的最大地表温度，代表特征空间中的干边；4)构建双抛物线型NDVI-Ts特征空间散点图，并根据双抛...

【专利技术属性】
技术研发人员：刘英，岳辉，侯恩科，
申请(专利权)人：西安科技大学，
类型：发明
国别省市：陕西;61