[发明专利]盐渍土光谱曲线可见光波段遥感反演构建方法

说明书

技术领域

本发明涉及遥感与光谱学分析领域，尤其涉及一种盐渍土光谱曲线可见光波段遥感反演构建方法。

背景技术

21世纪以来，各国的土地资源越来越匮乏，盐渍化土地作为潜在的可利用土地越来越得到重视，而借助遥感技术监测盐渍化土壤资源成为了研究的热点。盐渍土的遥感监测是基于盐渍土光谱曲线来实现的。目前，用于高光谱测量的仪器ASDPro FR地物光谱仪可以获取350-2500nm波长范围的盐渍化土壤光谱曲线，其光谱分辨率在可见光波段为1.4nm，在近红外波段为2nm，AvaField-2地物光谱仪，其光谱响应范围300-1700nm，在300-1040nm范围内采样间隔为0.59nm，在1040-1700nm范围内采样间隔为3.8nm等。

基于光谱特征的盐渍土遥感监测精度与传感器的光谱分辨率的紧密相关。为提高识别精度通常采用提高传感器精度的方法，这往往造成光谱资源的浪费和信息处理量的无谓增加，给盐渍土的遥感监测带来不必要的麻烦。光谱曲线遥感反演构建的最佳分辨率是指既能满足监测需求又不会造成数据冗余的采样间隔，最佳分辨率可以保证实现遥感信息的无失真监测。无失真盐渍土光谱曲线的获得是盐渍土遥感信息获取的重要内容之一，获得无失真盐渍土光谱曲线的方法是遥感技术研究前沿和热点，因此研发确定光谱曲线遥感反演构建最佳分辨率的方法具有重要的意义。

从信息论的角度看：确定光谱曲线遥感反演构建最佳光谱分辨率的本质是无失真重建信息（信号）的采样间隔确定。

针对目前基于光谱特征进行盐渍土遥感监测存在的问题是如何选择确定遥感监测最佳光谱分辨率的需求，开发一种既能够满足遥感监测需求，又不会造成数据冗余，还能无失真保存完整信息的光谱曲线遥感反演构建最佳分辨率的确定方法就极其重要。

发明内容

有鉴于现有技术的不足，本发明所要解决的技术问题是提供一种盐渍土光谱曲线可见光波段遥感反演构建方法，其能够满足盐渍土监测需求且能够大幅降低盐渍土光谱冗余资源占用、无失真保存完整信息。

为实现上述目的，本发明提供了一种盐渍土光谱曲线可见光波段遥感反演构建方法，具体地，本发明提供的技术方案如下：

一种盐渍土光谱曲线可见光波段遥感反演构建方法，包括如下步骤：

a)通过具有特定采样间隔的光谱仪器检测盐渍土样本，从而得到波长λ与盐渍土样本的反射率R变化的光谱曲线，光谱曲线具有多个数据点；

b)在多个数据点处分别计算一阶导数和二阶导数；

c)根据一阶导数和二阶导数分别在多个数据点处求光谱曲线的曲率半径，比较求得的多个曲率半径，以得到光谱曲线的最小曲率半径r_min和最小曲率半径的对应波长λ_rmin；

d)计算λ_rmin处的变化频率，得到光谱曲线的最大变化频率f_max；

e)根据香农采样定理，确定盐渍土样本光谱曲线可见光波段最佳光谱分辨率利用所述最佳光谱分辨率Δλ_OR得到重构光谱曲线。

进一步地，光谱曲线的波长选取范围大于等于400nm，小于1650nm。

进一步地，对光谱曲线各个数据点进行一阶导数、二阶导数运算采用如下方法：

R′(λ_i)=[R(λ_i+1)-R(λ_i-1)]/(2Δλ)；

R″(λ_i)=[R(λ_i+1)-2R(λ_i)+R(λ_i-1)]/(Δλ)²；

λ_i为数据点i处的波长，R(λ_i)为所述数据点i处的反射率，R′(λ_i)和R″(λ_i)分别为光谱曲线在波长λ_i处的一阶和二阶导数，Δλ为λ_i-1到λ_i的波长间隔。

进一步地，曲率半径的计算方法为：

r=|[(1+R′²)(^3/2)]/R″|

进一步地，最大频率变化按照如下方法计算：

①当r_min≥1时，f_max在数值上等于r_min，最佳光谱分辨率