一种基于热红外成像技术估算土壤蒸发强度的方法技术

本发明专利技术属于农田水文与水环境领域，涉及土壤蒸发强度的估算，具体来说是一种简便、快速且精度较高的基于热红外成像技术估算土壤蒸发强度的方法，包括如下步骤：步骤1、获取地表温度图像；步骤2、数据预处理；步骤3、去除温度矩阵干扰项；步骤4、估算蒸发强度；步骤5、绘制土壤蒸发强度等值线图；本发明专利技术应用广泛，可运用于大部分地区土壤蒸发强度的估算。地表温度可采用便携式热红外成像仪获取，对实验器材的限制较小，可在野外实现现场观测，省时省力。本发明专利技术克服了微型蒸发器测量土壤蒸发的局限性，避免了土壤异质性对测量值的干扰，克服了频繁更换土壤蒸发器所带来的测量不确定性。更换土壤蒸发器所带来的测量不确定性。更换土壤蒸发器所带来的测量不确定性。

【技术实现步骤摘要】
一种基于热红外成像技术估算土壤蒸发强度的方法


[0001]本专利技术属于农田水文与水环境领域，涉及土壤蒸发强度的估算，具体来说是一种简便、快速且精度较高的基于热红外成像技术估算土壤蒸发强度的方法。

技术介绍

[0002]土壤蒸发是指水分从土壤表面以气态的形式向大气逸散的过程。土壤蒸发过程受气象条件、土壤水分状态、地下水位、土壤理化特性等多种因素的影响，涉及陆地和大气之间的质量和热量交换，并直接影响地下水资源和地上水资源管理，是全球水文循环的重要过程之一。
[0003]在我国西北干旱灌区，由于大量引黄灌溉，导致地下水埋深普遍较浅，而强烈的蒸发过程导致地下水大量流失，造成了水资源的大量浪费。盐随水去，大量的盐分逐渐积累在地表，使土壤盐碱化更加严重，导致灌区大面积减产，并对当地生态环境造成了不可估量的影响。因此，土壤蒸发的精确监测及估算对灌区水资源高效利用和土壤盐渍化控制有着十分重要的作用。
[0004]当前的土壤蒸发测量方法主要有直接测量法、遥感反演法和水量平衡法三种。直接测量法多使用微型蒸发器进行测量，该方法具有简便、灵活等特点，但频繁更换土体难以确保测量的代表性，且存在以点代面的局限性，不适用于空间异质性较高的地区。遥感反演法是一种较为成熟的区域尺度的蒸散发估算方法，但不同模型之间反演的土壤蒸发量存在较大差异，在土壤蒸发的连续性监测和精确量化分析方面还存在一定的差距。水量平衡法则由于各平衡要素的精确获取存在一定的难度、所需参数多、计算过程相对较为复杂，因此也难以实现大范围的推广应用。
[0005]因此，有必要提出一种快速准确估算农田土壤蒸发强度的方法，以便指导农田灌溉，提高水分利用效率及降低土壤盐渍化风险。

技术实现思路

[0006]针对灌区土壤蒸发的测量存在精度低，代表性差以及计算过程复杂等问题，本专利技术的目的是提供一种基于热红外成像技术估算土壤蒸发强度的方法，该方法选用高分辨率的热红外成像仪作为主要工具，并依据该仪器所拍摄数据进行计算和定量分析。该方法可根据实际情况快速识别并计算出所测区域的土壤蒸发强度，具有测量精度高，范围广等特点，可为灌区水分管理及作物灌溉时间提供一定的技术支持。
[0007]为了达到上述目的，本专利技术提供了如下技术方案：
[0008]一种基于热红外成像技术估算土壤蒸发强度的方法，所述方法包括如下步骤：
[0009]步骤1、获取地表温度图像；
[0010]使用便携式热红外成像仪对所测区域进行拍摄，得到热红外图像；
[0011]步骤2、数据预处理；
[0012]将热红外图像导入软件FLUKE中；软件FLUKE依据实际气温矫正背景温度，并将拍
摄时的气温记为T0；将图像温度数据导出为txt格式，并从中提取温度矩阵P；对温度矩阵P进行预处理，使用小波去噪原理去除温度矩阵P中的高频部分，得到新的温度矩阵P1；
[0013]步骤3、去除温度矩阵干扰项；
[0014]土壤蒸发干扰项主要为植物及地面石块，根据干扰项的分布特点，使用两种聚类法依次进行聚类；
[0015]3a.温度矩阵聚类
[0016]首先采用基于密度的聚类方法进行聚类，得出若干个聚类，并分别计算各个小类的温度平均值；
[0017]其次，根据各个温度平均值进行基于距离的二次聚类，将聚类结果分为三类，温度最低的一类为植物温度矩阵，温度最高的一类为地面石块温度矩阵；
[0018]3b.将植物温度矩阵从温度矩阵P1原始矩阵中删除，得到温度矩阵P2；
[0019]3c.将地面石块温度矩阵从温度矩阵P2中删除，得到温度矩阵P3；
[0020]3d.计算温度矩阵P3的霍普金斯常量，若温度矩阵P3的霍普金斯常量为0.75～1.0之间，则判定仍然具有明显的聚类趋势，重复步骤3a到3d，直到温度矩阵P3的霍普金斯常量为0.5～0.75之间；若温度矩阵P3的霍普金斯常量为0.5～0.75之间，则判定分布较为均匀，进行下一步；
[0021]步骤4、估算蒸发强度；
[0022]4a.使用软件SPSS遍历出温度矩阵P3中的最大温度PX，将最大温度PX与空气温度T0之差记作M；
[0023]4b.估算土壤蒸发强度，包括如下步骤：
[0024]4b.1：将温度矩阵P3各项减去空气温度T0，得到温度矩阵P4；
[0025]4b.2：将所得温度矩阵P4各项均除以M，得到温度矩阵P5；
[0026]4b.3：根据当地气象数据获取净辐射值Rn，将净辐射值Rn除以2.45
×
106，记作R；
[0027]4b.4：将温度矩阵P5各项均除以R，所得矩阵即为蒸发强度矩阵P6；
[0028]步骤5、绘制土壤蒸发强度等值线图；
[0029]首先，将蒸发强度矩阵P6导入软件matlab中，绘制土壤蒸发强度热力图；其次，利用工具包contour绘制等值线图，调整图像的大小和显示效果。
[0030]步骤2中，对温度矩阵P进行预处理，使用小波去噪原理去除温度矩阵P中的高频部分，得到新的温度矩阵P1，包括如下步骤：
[0031]步骤2b.1：将温度矩阵P导入软件matlab中；
[0032]步骤2b.2：对温度矩阵P进行小波分解；设置小波分解层数为3层，并获得低频分量WT1和高频分量WT2；
[0033]步骤2b.3：分尺度去噪；对高频分量WT2进行阈值处理，当小波系数小于阈值时，将其置为零；当小波系数大于阈值时，将其减去阈值；阈值为：其中，σ为温度矩阵P的方差，M为温度矩阵P的行数，N为温度矩阵的列数；
[0034]步骤2b.4：小波逆变换；对处理后的分量进行小波重构，得到新的温度矩阵P1。
[0035]步骤3a中，包括如下步骤：
[0036]步骤3a.1：首先进行基于密度的聚类，即只要临近区域的密度超过某一阈值，就进
行聚类；根据聚类结果，将矩阵分别记为N1、N2、N3
…
，并计算每一分类的温度平均值，分别记为A1、A2、A3
…
；
[0037]步骤3a.2：将计算所得的A1、A2、A3
…
进行基于距离的k
‑
means聚类，将聚类结果设置为3类，记为T1、T2、T3；
[0038]具体聚类过程为：
[0039]1、随机从A1、A2、A3
…
选择其中三个，其值记作S1、S2、S3；
[0040]2、分别计算A1、A2、A3
…
到S1、S2、S3之间的距离；
[0041]3、分别判断A1、A2、A3
…
到S1、S2、S3的距离，若距离S1近，则归为T1；若距离S2近，则归为T2；若距离S3近，则归为T3；
[0042]4、计算T1、T2、T3的平均值TA1、TA2、TA3，计算与S1、S2、S3与TA1、TA2、TA3的距离，若距离小于1，则结束本步骤；若大于本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于热红外成像技术估算土壤蒸发强度的方法，其特征在于：所述方法包括如下步骤：步骤1、获取地表温度图像；使用便携式热红外成像仪对所测区域进行拍摄，得到热红外图像；步骤2、数据预处理；将热红外图像导入软件FLUKE中；软件FLUKE依据实际气温矫正背景温度，并将拍摄时的气温记为T0；将图像温度数据导出为txt格式，并从中提取温度矩阵P；对温度矩阵P进行预处理，使用小波去噪原理去除温度矩阵P中的高频部分，得到新的温度矩阵P1；步骤3、去除温度矩阵干扰项；土壤蒸发干扰项主要为植物及地面石块，根据干扰项的分布特点，使用两种聚类法依次进行聚类；3a.温度矩阵聚类首先采用基于密度的聚类方法进行聚类，得出若干个聚类，并分别计算各个小类的温度平均值；其次，根据各个温度平均值进行基于距离的二次聚类，将聚类结果分为三类，温度最低的一类为植物温度矩阵，温度最高的一类为地面石块温度矩阵；3b.将植物温度矩阵从温度矩阵P1原始矩阵中删除，得到温度矩阵P2；3c.将地面石块温度矩阵从温度矩阵P2中删除，得到温度矩阵P3；3d.计算温度矩阵P3的霍普金斯常量，若温度矩阵P3的霍普金斯常量为0.75～1.0之间，则判定仍然具有明显的聚类趋势，重复步骤3a 到3d，直到温度矩阵P3的霍普金斯常量为0.5～0.75之间；若温度矩阵P3的霍普金斯常量为0.5～0.75之间，则判定分布较为均匀，进行下一步；步骤4、估算蒸发强度；4a.使用软件SPSS遍历出温度矩阵P3中的最大温度PX，将最大温度PX与空气温度T0之差记作M；4b.估算土壤蒸发强度，包括如下步骤：4b.1：将温度矩阵P3各项减去空气温度T0，得到温度矩阵P4；4b.2：将所得温度矩阵P4各项均除以M，得到温度矩阵P5；4b.3：根据当地气象数据获取净辐射值Rn，将净辐射值Rn除以2.45
×
106，记作R；4b.4：将温度矩阵P5各项均除以R，所得矩阵即为蒸发强度矩阵P6；步骤5、绘制土壤蒸发强度等值线图；首先，将蒸发强度矩阵P6导入软件matlab中，绘制土壤蒸发强度热力图；其次，利用工具包contour绘制等值线图，调整图像的大小和显示效果。...

【专利技术属性】
技术研发人员：黄权中，李金帅，黄冠华，刘佳伟，李泽，刘娜，徐旭，霍再林，熊云武，
申请(专利权)人：中国农业大学，
类型：发明
国别省市：