具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

下面结合图1-图3描述本发明的水生植被的类型确定方法及装置。

图1是本发明实施例提供的水生植被的类型确定方法的流程示意图之一。

如图1所示，本发明实施例提供的一种水生植被的类型确定方法，包括以下步骤：

101、确定待检测下垫面是否为水体。

具体的，确定待检测下垫面是否为水体，包括：确定待检测下垫面的遥感影像的波段；根据波段，确定改进的归一化水体指数；根据改进的归一化水体指数，确定待检测下垫面是否为水体。其中，遥感影像的波段分为两种情况来进行分析，一种情况是为4个波段，另一种情况是超过4个波段。

当波段只有4个波段时，获取待检测下垫面的归一化植被指数NDVI和归一化水体指数NDWI；获取待检测下垫面的归一化植被指数包括获取待检测下垫面的近红外波段反射率R_nir和红波段反射率R_r；基于预设关系，根据近红外波段反射率R_nir和红波段反射率R_r，获取归一化植被指数NDVI。预设关系为式(1)：

式中，NDVI为归一化植被指数，R_nir为近红外波段反射率，R_r为红波段的反射率。

获取归一化水体指数便是根据红光波段反射率R_g和近红外波段反射率R_nir，获取公式如式(2)：

式中，NDWI为归一化植被指数，R_nir为近红外波段反射率，R_g为红光波段的反射率。

于是便计算得到了归一化植被指数NDVI和归一化水体指数NDWI，然后根据归一化植被指数NDVI和和归一化水体指数NDWI，确定待检测下垫面的水体指数WI作为改进的归一化水体指数，其中，WI的计算过程如式(3)：

WI＝NDVI-NDWI (3)

式中，WI为水体指数作为改进的归一化水体指数。

因此，便确定了当为4个波段的遥感影像时的改进的归一化水体指数，即WI。

当波段超过4个波段时，根据红光波段反射率R_g和短波红外波段反射率R_swir，确定改进的归一化水体指数MNDWI，如式(4)，或者是用WI模型计算改进的归一化水体指数。

式中，MNDWI为改进的归一化水体指数，R_swir为短波红外波段反射率，R_g为红光波段反射率。

于是，便分别计算得到了两种波段情况下对应的改进的归一化水体数值，然后再根据得到的改进的归一化水体指数，确定待检测下垫面是否为水体。具体的则是当MNDWI≥a₁或WI≤a₂时，则判定待检测下垫面为水体，否则为非水体，然后对非水体进行掩膜处理。其中a₁、a₂为设定阈值，可以根据具体数据分析进行人为确定，最终便确定了待检测下垫面是否为水体。

102、当待检测下垫面为水体时，基于归一化植被指数，确定水体中是否有水生植被。

当确定待检测下垫面为水体时，通过上述步骤便获取到了水体空间范围和水体NDVI的空间分布。然后确定NDVI与预设植被阈值a₃的大小，其中预设植被阈值a₃为人为设定值，根据实际情况进行判定，而确定水体中是否存在水生植被，则通过式(5)，如下：

NDVI≥a₃ (5)

当满足(5)中的条件时，便可以确定此时的水体中含有水生植被。

103、当水体中有水生植被时，通过微波遥感获取水生植被的后向散射系数。

当确定水体中有水生植被时，获取覆盖该水体区域的微波遥感波段及其极化方式，便得到了水生植被生长期的时间序列的后向散射系数δ_t。

104、基于归一化植被指数和后向散射系数，确定水生植被的类型，水生植被的类型包括：挺水植物、浮叶植物或潜水植物。

具体的，基于归一化植被指数和后向散射系数，确定水生植被的类型，包括：根据后向散射系数，确定水生植被的挺水植物识别指数BSI；若挺水植物识别指数BSI大于等于预设挺水植物指数阈值b₁，即BSI≥b₁，则确定水生植被为挺水植物，否则为浮叶植物或潜水植物。而确定水生植被的挺水植物识别指数BSI，则按照式(6)的方式进行确定，如下：

式中，BSI为挺水植物识别指数，δ_i为水生植物生长期时间序列的微波数据后向散射系数；i为参与计算的水生植被生长期微波遥感影像的第i景数据；n₁为遥感影像的景数；为时间序列的δ_i均值。

而在确定水生植被不是挺水植物时，也就是挺水植物识别指数BSI小于预设挺水植物指数阈值b₁时，则确定水生植被的反射系数RFI和水体植被差异指数VVI_i；若反射系数RFI大于等于预设反射系数阈值c₁，即式(7)，或水体植被差异指数VVI_i大于等于预设水体植被差异指数阈值c₂，即式(8)，则确定水生植被为浮叶植物，而执行式(7)、(8)的前提是NDVI≥a₃，且BSI＜b₁。

而当不满足上述式(7)和式(8)的条件时，则确定水生植被的温度差异指数TVI_i；根据温度差异指数TVI_i和水体植被差异指数VVI_i，确定温度植被差异指数DFI；若水体植被差异指数VVI_i小于预设水体植被差异指数阈值c₂，即式(9)；或温度植被差异指数DFI小于等于预设温度植被差异指数阈值c₃，即式(10)；则确定水生植被为潜水植物，其中，执行(9)(10)的前提是NDVI≥a₃，且BSI＜b₁。

DFI＝VVI_i-TVI_i≤c₃ (10)

式中，VVI_i为水体植被差异指数，DFI为温度-植被差异指数，TVI_i为温度差异指数。

而确定水生植被的温度差异指数TVI_i，包括：获取影像采集时的空气温度T_i；确定与遥感影像同步获取的时间序列的空气温度均值可以是通过地面气象站获取的空气温度均值，根据空气温度和空气温度均值，确定水生植被的温度差异指数TVI_i，如式(11)：

于是，便通过上述的公式计算之后，确定出水生植被的具体的植被类型，确定出具体水生植被的植被类型为挺水植物、浮叶植物还是具体的潜水植物，其中，潜水植物可以包括沉水植物和藻类等。

本发明实施例提供的一种水生植被的类型确定方法，通过确定待检测下垫面是否为水体；当所述待检测下垫面为水体时，基于归一化植被指数，确定所述水体中是否有水生植被；当所述水体中有水生植被时，通过微波遥感获取所述水生植被的后向散射系数；基于所述归一化植被指数和所述后向散射系数，确定所述水生植被的类型，所述水生植被的类型包括：挺水植物、浮叶植物或潜水植，通过在光学遥感的基础上加入微波遥感技术，既能考虑到水生植被的立体信息，又能考虑到水生植被的物候和光谱特征，能够实现准确、高效的对水生植被类型的识别确定。

进一步的，在上述实施例的基础上，本实施例中在确定水生植被的类型以后，还包括：对挺水植物、浮叶植物和潜水植物进行整合分类；根据整合分类的结果，生成水生植被类型图，使得可以更加直观地了解到具体的水生植被类型。

通过结合微波遥感和光学遥感技术的优势，既考虑挺水植被立体信息及其变化特征，又考虑水生植被的物候和光谱特征，使得在面向河湖水生植被和水环境管理时，充分考虑挺水植被、浮叶植被和潜水植被的后向散射系数和光学水体指数变化特征的基础上，能够更加准确、快速地完成对水生植被类型的识别。

基于同一总的发明构思，本申请还保护一种水生植被的类型确定装置。下面对本发明提供的水生植被的类型确定装置进行描述，下文描述的水生植被的类型确定装置与上文描述的水生植被的类型确定方法可相互对应参照。

图2是本发明实施例提供的水生植被的类型确定装置的结构示意图。

如图2所示，本实施例提供的一种水生植被的类型确定装置，包括：

水体确定模块10，用于确定待检测下垫面是否为水体；

植被确定模块20，用于当所述待检测下垫面为水体时，基于归一化植被指数，确定所述水体中是否有水生植被；

微波遥感模块30，用于当所述水体中有水生植被时，通过微波遥感获取所述水生植被的后向散射系数；

类型确定模块40，用于基于所述归一化植被指数和所述后向散射系数，确定所述水生植被的类型，所述水生植被的类型包括：挺水植物、浮叶植物或潜水植物。

本发明实施例提供的一种水生植被的类型确定装置，通过确定待检测下垫面是否为水体；当所述待检测下垫面为水体时，基于归一化植被指数，确定所述水体中是否有水生植被；当所述水体中有水生植被时，通过微波遥感获取所述水生植被的后向散射系数；基于所述归一化植被指数和所述后向散射系数，确定所述水生植被的类型，所述水生植被的类型包括：挺水植物、浮叶植物或潜水植，通过在光学遥感的基础上加入微波遥感技术，既能考虑到水生植被的立体信息，又能考虑到水生植被的物候和光谱特征，能够实现准确、高效的对水生植被类型的识别确定。

进一步的，本实施例中的水体确定模块10，具体用于：

确定待检测下垫面的遥感影像的波段；

根据波段，确定改进的归一化水体指数；

根据改进的归一化水体指数，确定待检测下垫面是否为水体。

进一步的，本实施例中水体确定模块10，具体用于：

当确定待检测下垫面的遥感影像的波段为4个时，所述根据所述波段，确定改进的归一化水体指数，包括：

获取待检测下垫面的归一化植被指数和归一化水体指数；

根据所述归一化植被指数和所述归一化水体指数，得到水体指数作为改进的归一化水体指数；

当确定待检测下垫面的遥感影像的波段超过4个时，所述根据所述波段，确定改进的归一化水体指数，包括：

根据短波红外波段反射率和红光波段反射率，确定改进的归一化水体指数；或者，

以所述水体指数作为改进的归一化水体指数。

进一步的，本实施例中的水体确定模块，具体还用于：

获取待检测下垫面的近红外波段反射率和红波段反射率；

基于预设关系，根据所述近红外波段反射率和所述红波段反射率，获取归一化植被指数；

根据所述红光波段反射率和所述近红外波段反射率，获取归一化水体指数。

进一步的，本实施例中的类型确定模块40，具体用于：

根据后向散射系数，确定水生植被的挺水植物识别指数；

若挺水植物识别指数大于等于预设挺水植物指数阈值，则确定水生植被为挺水植物。

进一步的，本实施例中的类型确定模块40，具体还用于：

若挺水植物识别指数小于预设挺水植物指数阈值，则确定水生植被的反射系数和水体植被差异指数；

若反射系数大于等于预设反射系数阈值，或水体植被差异指数大于等于预设水体植被差异指数阈值，则确定水生植被为浮叶植物。

进一步的，本实施例中的类型确定模块40，具体还用于：

确定水生植被的温度差异指数；

根据温度差异指数和水体植被差异指数，确定温度植被差异指数；

若水体植被差异指数小于预设水体植被差异指数阈值，或温度植被差异指数小于等于预设温度植被差异指数阈值，则确定水生植被为潜水植物。

进一步的，本实施例中的类型确定模块40，具体还用于：

获取影像采集时的空气温度；

确定与遥感影像同步获取的时间序列的空气温度均值；

根据空气温度和空气温度均值，确定水生植被的温度差异指数。

进一步的，本实施例中，还包括整合分类模块，用于：

对挺水植物、浮叶植物和潜水植物进行整合分类；

根据整合分类的结果，生成水生植被类型图。

图3是本发明实施例提供的电子设备的结构示意图。

如图3所示，该电子设备可以包括：处理器(processor)310、通信接口(Communications Interface)320、存储器(memory)330和通信总线340，其中，处理器310，通信接口320，存储器330通过通信总线340完成相互间的通信。处理器310可以调用存储器330中的逻辑指令，以执行水生植被的类型确定方法，该方法包括：确定待检测下垫面是否为水体；当所述待检测下垫面为水体时，基于归一化植被指数，确定所述水体中是否有水生植被；当所述水体中有水生植被时，通过微波遥感获取所述水生植被的后向散射系数；基于所述归一化植被指数和所述后向散射系数，确定所述水生植被的类型，所述水生植被的类型包括：挺水植物、浮叶植物或潜水植物。

此外，上述的存储器330中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法所提供的水生植被的类型确定方法，该方法包括：确定待检测下垫面是否为水体；当所述待检测下垫面为水体时，基于归一化植被指数，确定所述水体中是否有水生植被；当所述水体中有水生植被时，通过微波遥感获取所述水生植被的后向散射系数；基于所述归一化植被指数和所述后向散射系数，确定所述水生植被的类型，所述水生植被的类型包括：挺水植物、浮叶植物或潜水植物。

又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各提供的水生植被的类型确定方法，该方法包括：确定待检测下垫面是否为水体；当所述待检测下垫面为水体时，基于归一化植被指数，确定所述水体中是否有水生植被；当所述水体中有水生植被时，通过微波遥感获取所述水生植被的后向散射系数；基于所述归一化植被指数和所述后向散射系数，确定所述水生植被的类型，所述水生植被的类型包括：挺水植物、浮叶植物或潜水植物。

以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。