土壤是陆地上万千生物赖以生存的“家”。土壤的分层信息是反映其质量的重要指标之一。探地雷达具有良好的穿透性，无需剖面就可以掀开土壤分层信息的盖头，让我们清晰了解土壤结构。
土壤分层示意
基于包络检波和STFT谱分析的探地雷达土壤分层
电磁波在土壤介质中传播将发生衰减，并且由于地面坡度、粗糙度等几何特性的不均一性，以及地表植被散射影响，接收机接收到的回波信号往往具有复杂的散射过程，包含各种各样的杂波干扰，使得土层识别十分困难。
因此，本研究首先对雷达回波信号进行增益处理来弥补信号在传播过程中的能量损失，并进行滤波与降噪处理以降低噪声，从而提高回波信号信噪比。具体滤波处理过程为：通过均值法去除背景的影响；通过带通滤波去除外界射频信号的干扰；通过中值法去除来自地面反射的直达波。然后分别在时域和频域对回波信号进行处理以获取时域分层信息和估计各层电磁波传播速度：一方面采用包络检波的方法进行土层分界面的时域确定，进而确定土层分界面对应传播时间信息；另一方面采用短时傅里叶变换对信号进行时频域分析，计算基于时域分层的各土壤层的频域振幅峰值加权值，利用建立的频域振幅峰值加权值与土壤介电常数的线性回归关系计算各层的土壤介电常数，从而确定电磁波在各土层的传播速度，最后结合土层界面时间信息确定各土层厚度信息。
为研究探地雷达对不同土层深度和不同土层界面的敏感性，在北京通州中国农业大学实验站搭建模拟实验环境。设计10 m×10 m的区域，等分成5 m×10 m的两块（左边以石子为障碍层，右边以瓷砖为障碍层），挖掘出10条长度为10 m、宽度为1 m、高度以10 cm为步长逐渐递增的阶梯状平面，分别以2.0系列石子和瓷砖为界面在长度10 m的阶梯上左右两边均匀铺设5 m×1 m（长×宽）的障碍层，最后覆土整平。覆土后最深处障碍层位于土壤下1 m位置，障碍层厚度约1 cm，其上下土壤为实验田常规砂质壤土，如下图所示。
北京通州实验站模拟环境搭建及测线分布
考虑到模拟环境的阶梯式土层厚度最深处为1 m，因此选用探测深度为0.8~1.2 m的CAS-S800（工作频率800 MHz）探地雷达在其上方以轮测法沿4条测线进行数据采集，每条测线长度为10 m左右，采样时窗设置为40 ns，采样频率设置为26 324 MHz进行数据采集。电磁波在介质中的平均速率根据第一障碍层位置估算。
可以看出，三维土壤分层识别图中准确识别出了前6层台阶的位置，信号随着障碍层深度的增加会逐渐减弱，障碍层深度不大于65 cm时都能提取出障碍层位置，障碍层深度大于65 cm时回波信号极为微弱。从二维结果中可以看出，每个天线位置获取的数据中包含了当地对应台阶的位置，通过包络检波得到的上下包络线已经比回波曲线本身更能直观地显示出回波跳变的位置，但只通过看上下包络线距离的变化有时候还不够直观，特别是当障碍层深度增加，信号变化较小时。因此，本文进一步通过识别其Hilbert谱瞬时相位跳变顺利找到了造成跳变的障碍层位置，结果如上图的三维结果所示。
在复杂的实际农田中本方法有效性如何？一起来看看吧。本研究选取位于114°50′22″E —115°13′27″E，36°35′43″N —36°57′N之间的河北省曲周县部分农田区域作为实验采样区域作为实验采样区域。研究区面积约为667 km 2，土壤类型以潮土为主，质地类型主要为轻壤土。土壤在垂直方向上砂壤黏交替沉积呈层状，土壤呈现出一定的分层特征。
位于河北省曲周县的研究区农田环境采样图
在2016年4月开展了外业实验，根据项目区常见的4种土体结构类型（通体粘、夹粘型、倒蒙金型、底漏型）分别选取了8个采样点进行了探地雷达数据采集、剖面数据采集和土壤样品采样。由于调研时间为2016年4月，项目区大部分农田生长着小麦、玉米等作物，有的样点数据不是很理想，因此剔除干扰、噪声较大的样点，最终选择了裸露地表的样点进行测算和验证，每种土体构型选择一个样本，共计4个采样点。
实际农田环境中选用探测深度更深、工作频率为400 MHz的GR-IV型便携式探地雷达进行数据采集。采样点位置如图所示（底图来自于Google Earth高清影像），每个采样点根据田块大小设置多条10~25 m长度的测线进行往返数据采集，采样时窗设置为25 ns，采样频率设置为100 000 MHz。
基于包络检波土壤分层信息识别研究区4个采样点的土壤层位置。实验表明，算法能在复杂的农田环境中检测出不同土壤层的存在。以样点1为代表给出了分层界面三维显示。为以更直观的二维图形式展现回波数据及其上下包络线情况，图中也给出了天线在随机选取的几个位置上样点1的二维数据及包络线。
河北省曲周县研究区典型样点土壤分层情况
实际农田环境中，存在着一些不确定性因素，从而为土层厚度测量带来一定误差。这些不确定因素包括：
1
农田环境的不确定性
由于农田环境中地表的复杂性和非均一性，例如地表植被、地上覆膜、土壤空洞以及土层盐分或者水分的不均匀等，使得探地雷达回波图像中出现异常的信号突变。数据采集时间为2016年4月，虽然我们选择了裸露地表或者植被稀疏地表，但杂草等地表覆盖物依然不可避免，有的地方有前一年的地表残膜存在。另外曲周地区以前存在大量的盐碱地，经过多年的治理，虽然有了很大改善，但个别样点可能依然存在盐度较高的土块。这些都将带来探地雷达回波信号的不稳定。预处理中我们通过信号滤波降低了一部分由这些突变带来的噪声，但由土壤中较大空洞等造成的回波变化依然会引起其下部土层的层厚测定误差。
2
土壤层的渐变性
在农田环境中土壤质地的变化通常是渐变的，介电常数也因此是渐变的，为了测算方便，本文将每层介电常数设为一常数进行计算，且实际通过介电常数对速度的估计也并非十分准确，而土层厚度的计算是根据回波时延和所估算出的速度来确定的，因此速度估算的误差使得土层厚度的测算也存在一定误差。然而，这种局限性和约束性对土层厚度测算的影响在实际应用中通常处于可接受的范围。
3
土层深度的影响
然而，应用的道路并不是坦途一片。探地雷达信号在不同土壤环境中穿透性、反射性都有比较大的差异，本研究在位于华北平原的河北曲周进行了实际农田环境实验，验证了本方法在华北平原地区的适用性。对于我国其它区域，由于土壤质地差异较大，雷达信号的穿透性和受干扰情况不尽相同，具体扩展到其它地区的实验将在下一步研究。此外，探地雷达采样时应注意避开作物覆盖的时间或样点，根据研究区土壤黏度和采样时土壤水分情况合理设置探地雷达参数。因此，在以后的工作中我们将考虑采用更高功率或更低频率的探地雷达进行数据获取，以提高信号的穿透性。同时，还将考虑尝试其它从弱信号中检测土层信息的方法。
原文可在中国知网或期刊官网下载
作者简介
中国农业大学土地科学与技术学院副教授、博士生导师。从事微波遥感及其在土地和农业监测中的应用研究。农业农村部农业灾害遥感重点实验室和自然资源部农用地质量与监控实验室骨干成员。先后获得2019年中国测绘学会测绘科技进步奖一等奖、2016-2017年中国农学会神农中华农业科技奖一等奖、2011年与2012年国土资源科学技术奖二等奖等奖项。主持和参加国家自然科学基金、国土资源部和农业部业务项目等多项项目。近年，在Remote Sensing 、International Journal of Remote Sensing、《地球信息科学学报》《资源科学》《农业工程学报》《农业机械学报》等国内外期刊发表论文30余篇，申请专利6项，软件著作权20余项。
《地球信息科学学报》是由中国科学院主管，由中国科学院地理科学与资源研究所、中国地理学会联合主办，由资源与环境信息系统国家重点实验室（中国科学院地理科学与资源研究所）、虚拟地理环境教育部重点实验室（南京师范大学）、三维信息获取与应用教育部重点实验室（首都师范大学）联合协办的地球信息科学领域综合性学术期刊，月刊。
学报创刊于1996年，创始人为陈述彭院士。学报现任主编徐冠华院士，执行主编陆锋研究员。
学报主要刊登以地球系统信息流为主要研究对象、以地球时空信息机理、时空信息认知与分析方法、时空信息可视化与信息图谱等为研究主题的原创性学术论文及相关评论与简讯。特别重视地球信息科学及相关交叉学科的最新研究成果的报道。学报设置了地理信息科学理论与方法、地理空间分析综合应用、遥感科学与应用技术等栏目，并面向国家战略需求，依托重大科技计划项目成果组织专刊，全面反映地球信息科学的最新研究进展。
学报是中国科学引文数据库（CSCD）核心期刊、中国科技核心期刊、全国中文核心期刊。2020年《中国学术期刊影响因子年报》显示，学报影响因子位列测绘科学技术期刊第2名。在2020年版中国科学院科学出版基金中文科技期刊排行榜中，学报排名第25位。
学报坚持贯彻百花齐放、百家争鸣的精神，秉持公平公正客观、严谨认真的原则，不断提高办刊质量。学报坚持追踪本学科领域的前沿与热点，突出地球信息科学理论与方法及地球观测技术的办刊特色，兼顾学科交叉的示范性应用，并不断扩大国内和国际影响力。
转载自地球信息科学学报
经作者授权转载
文章仅代表作者观点，与本公众号无关，版权归原作者所有
原文标题：借探地雷达看土壤分层
排版：倪宁淇返回搜狐，查看更多