一种基于多功能雷达慢扫增强的浒苔提取方法
技术领域
本发明涉及海洋遥感领域,特别是一种基于多功能雷达慢扫增强的浒苔提取方法。
背景技术
目前,国内外专家学者对浒苔遥感监测方法的研究主要有:光学遥感监测、辐射传输模型及微波遥感监测。光学遥感监测与辐射传输模型,是利用浒苔覆盖海水在可见光波段的波谱特性,提取反演信息。微波遥感监测方法,主要是利用浒苔与水体间的后向散射系数差异,通过图像分割算法,实现水面浒苔信息的提取。但由于卫星数据的时效性较差,不能实时监控海面浒苔信息。
申请人之前申报的申请号为201610822712.0的中国专利申请,其发明创造的名称为“基于相参雷达慢扫的自适应溢油信息提取方法”,利用相参雷达提取溢油信息,并且考虑了海面风场对溢油提取的影响,但是海面风场的影响,没有海浪波高和波长的影响更为直接、准确,并且该专利没有涉及幅度值对目标探测的影响。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足,而提供一种基于多功能雷达慢扫增强的浒苔提取方法,该基于多功能雷达慢扫增强的浒苔提取方法基于实际经验,并结合多功能雷达探测技术,经过信号处理、图像分割、特征提取,结合浪流信息、幅度特征建立浒苔样本库,采用BP神经网络模型训练样本。因此,相对于传统的算法,本发明具有时效性更强,识别精度更高,预测精度更准的特点。
为解决上述技术问题,本发明采用的技术方案是:
一种基于多功能雷达慢扫增强的浒苔提取方法,包括如下步骤。
步骤1,浪流信息提取:根据各散射单元积累的脉冲数,获得监测区域内浪流信息。
步骤2,图像分割:采用OSTU分割算法,将雷达视频数据进行分割,提取疑似浒苔区域。
步骤3,纹理特征与几何特征提取:对步骤2提取的疑似浒苔区域进行纹理特征和几何特征的提取;几何特征包括复杂度和形状因子。
步骤4,幅度特征提取:对疑似浒苔区域的雷达视频数据进行快速傅里叶变换,依据各径向上的频谱信息,提取该疑似浒苔区域的幅度特征。
步骤5,浒苔目标区域提取:根据提取的浪流信息、纹理特征、几何特征与幅度特征,组建浒苔识别样本库,并训练BP神经网络模型;利用训练好的BP神经网络模型,从疑似浒苔区域中识别并提取出浒苔目标区域。
步骤3中,采用灰度共生矩阵算法提取疑似浒苔区域的纹理特征和几何特征。
步骤4中,将疑似浒苔区域内的径向回波数据进行快速傅里叶变换,得到雷达回波的频谱信息,将频谱能量最大值对应的频率进行频域转换,得到该径向的径向幅度值C,继而提取疑似浒苔区域的幅度特性;幅度特性包括均值、标准差和偏度。
幅度特性中均值的计算公式为:
式中,μ为均值,C为径向幅度值,N为径向幅度值的个数。
幅度特性中标准差的计算公式为:
式中,σ为标准差,μ为均值,C为径向幅度值,N为径向幅度值的个数。
幅度特性中偏度的计算公式为:
式中,γ为偏度,σ为标准差,μ为均值,C为径向幅度值,N为径向幅度值的个数。
步骤5中,BP神经网络模型包括输入层、隐含层和输出层;其中,输入层节点数innode为11,11个输入层节点分别为波高、波长、角二阶矩、对比度、倒数差分矩、二阶熵、复杂度因子、形状因子、幅度均值、幅度标准差和幅度偏度;隐含层节点数hidenode的取值范围为10~16;输出层节点数为2;BP神经网络模型的信号正向传递,输入层、隐含层和输出层中各节点均选用S形函数为激活函数,S形函数的公式如下:
其中,
式中,y为每个节点的输出值,w为节点之间的连接权值,b为各节点之间的阈值,S为每个节点的输入值。
步骤5中,BP神经网络模型的误差信号反向传递,误差函数E(w,b)的计算公式为:
式中,yj为输出层的结果;dj为实际观察结果;E为误差,允许误差allowerror的取值范围为0.0001~0.001。
步骤5中,BP神经网络模型的训练方法,包括如下步骤。
步骤51,建立输入层矩阵:建立一个innode*N的矩阵,其中N为训练样本的数量,从第一列至第innode列依次为:波高、波长、角二阶矩、对比度、倒数差分矩、二阶熵,复杂度因子、形状因子、幅度均值、幅度标准差和幅度偏度。
步骤52,建立样本目标类型矩阵:建立一个1*N矩阵,为样本目标类型,与训练样本相对应。
步骤53,将每列数据进行归一化,取值范围在(0,1)之间。
步骤54,连接权值调整:依次将每一行数据导入到输入层,计算误差,并将计算得出的误差与允许误差进行比较,当超过允许误差时,调整连接权值。
步骤55,当计算误差满足训练要求中的允许误差时,确定连接权值。
步骤52中,样本目标类型t为0时,表示该样本为非浒苔区域;样本目标类型t为1时,表示该样本为浒苔区域。
本发明具有如下有益效果:能够全天候、不间断监测识别海面浒苔,并实时获取雷达扫描范围内的浒苔面积、移动方向与移动速度。本发明的优势在于:在建立浒苔识别模型时,加入了浪流的波长与波高影响因子,提高了浒苔目标的识别精度;在预测浒苔运动时,基于多普勒原理,获取浒苔目标区域的瞬时速度,提高了浒苔运动的预测精度。
附图说明
图1显示了本发明一种基于多功能雷达慢扫增强的浒苔提取方法的结构示意图。
图2显示了本发明中BP神经网络模型的模型示意图。
具体实施方式
下面结合附图和具体较佳实施方式对本发明作进一步详细的说明。
如图1所示,一种基于多功能雷达慢扫增强的浒苔提取方法,包括如下步骤。
步骤1,浪流信息提取:根据各散射单元积累的脉冲数,获得监测区域内浪流信息。监测区域内浪流信息的提取方法,具体参见申请号为201610766600.8的中国专利申请,其发明名称为“基于相参雷达慢扫模式下的浪流信息提取方法”。
步骤2,图像分割:采用OSTU分割算法,将雷达视频数据进行分割,提取疑似浒苔区域。疑似浒苔区域的提取方法,具体参见申请号201610822712.0的中国专利申请,其发明名称为“基于相参雷达慢扫的自适应溢油信息提取方法”。
步骤3,纹理特征与几何特征提取:优选采用灰度共生矩阵算法,对步骤2提取的疑似浒苔区域进行纹理特征和几何特征的提取;几何特征包括复杂度和形状因子等。
纹理特征与几何特征的具体提取方法,详细参见申请号201610822712.0的中国专利申请,其发明名称为“基于相参雷达慢扫的自适应溢油信息提取方法”。
步骤4,幅度特征提取:对疑似浒苔区域的雷达视频数据,优选为疑似浒苔区域内的径向回波数据,进行快速傅里叶变换,依据各径向上雷达回波的频谱信息,将频谱能量最大值对应的频率进行频域转换,得到该径向的径向幅度值C,继而提取疑似浒苔区域的幅度特性;幅度特性包括均值、标准差和偏度等。
其中,幅度特性中均值的计算公式为:
式中,μ为均值,C为径向幅度值,N为径向幅度值的个数。
幅度特性中标准差的计算公式为:
式中,σ为标准差,μ为均值,C为径向幅度值,N为径向幅度值的个数。
幅度特性中偏度的计算公式为:
式中,γ为偏度,σ为标准差,μ为均值,C为径向幅度值,N为径向幅度值的个数。
步骤5,浒苔目标区域提取:根据提取的浪流信息、纹理特征、几何特征与幅度特征,组建浒苔识别样本库,并训练BP神经网络模型;利用训练好的BP神经网络模型,从疑似浒苔区域中识别并提取出浒苔目标区域。
如图2所示,BP神经网络模型包括输入层、隐含层和输出层。其中,输入层节点数innode为11,11个输入层节点分别为波高、波长、角二阶矩(ASM)、对比度(CON)、倒数差分矩(IDF)、二阶熵(ESO),复杂度因子、形状因子、幅度均值、幅度标准差和幅度偏度。隐含层节点数hidenode的取值范围为10~16;输出层节点数为2。
BP神经网络模型的信号正向传递,输入层、隐含层和输出层中各节点均优选采用S形函数为激活函数,S形函数的公式如下:
其中,
式中,y为每个节点的输出值,w为节点之间的连接权值,b为各节点之间的阈值,S为每个节点的输入值。
BP神经网络模型的误差信号反向传递,误差函数E(w,b)的计算公式为:
式中,yj为输出层的结果;dj为实际观察结果;E为误差,允许误差allowerror的取值范围为0.0001~0.001。
BP神经网络模型训练的主要目的是反复修正权值和阈值,使得误差函数值达到最小。BP神经网络模型的训练方法,包括如下步骤。
步骤51,建立输入层矩阵:建立一个innode*N的矩阵,其中N为训练样本的数量,从第一列至第innode列依次为:波高、波长、角二阶矩、对比度、倒数差分矩、二阶熵,复杂度因子、形状因子、幅度均值、幅度标准差和幅度偏度。
建立的输入层矩阵如下表所示:
步骤52,建立样本目标类型矩阵:建立一个1*N矩阵,为样本目标类型,与训练样本相对应。
建立的样本目标类型矩阵,如下表所示:
[1 2 3 … N]
样本目标类型t为0时,表示该样本为非浒苔区域;样本目标类型t为1时,表示该样本为浒苔区域。
步骤53,将每列数据进行归一化,取值范围在(0,1)之间。
步骤54,连接权值调整:依次将每一行数据导入到输入层,计算误差,并将计算得出的误差与允许误差进行比较,当超过允许误差时,调整连接权值。
步骤55,当计算误差满足训练要求中的允许误差时,确定连接权值。
本发明相比于纹理特征与几个特征,浒苔区域与背景区域的幅度特征差异更加明显,故本发明采用FFT手段提取幅度信息,提升了浒苔的识别精度。在构建BP神经网络模型时,综合考虑了海浪对浒苔特征的影响,使浒苔监测更加精准,并且具有全天候监测、时效性强的特点
以上详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种等同变换,这些等同变换均属于本发明的保护范围。