一种SAR图像海洋内波参数提取方法与流程

本发明涉及图像特征提取方法，具体涉及一种sar图像海洋内波参数提取方法。
背景技术：
：海洋内波是一种较为普遍存在的海洋现象，其发生机制的复杂性及其时空上的随机性使得内波探测称为海洋领域的前沿性研究课题。合成孔径雷达作为一种高分辨率微波成像雷达，具有全天时、全天候、穿透性等优点，广泛应用于海洋内波遥感观测。借助于气象数据和水文数据，通过对sar图像中核心参数的提取可实现海洋内波动力参数的计算。受成像条件的影响，海洋sar图像的相干斑及背景杂波较为严重、图像辐射分辨率较低，造成海洋内波纹理特性的退化，甚至引起亮、暗条纹的丢失或损坏。因此参数提取是sar图像目标识别的重要环节。sar图像海洋内波参数通常指波长。目前sar图像海洋内波参数提取方法面临如下问题：(1)内波参数提取方法中内波流体动力学模型一般采用经典的kdv方程，该方程不是任何条件下都适用的；(2)同时现有方法中振幅的提取依赖于提取波长的精度，而实际中由于成像条件的影响，很多情况下sar图像中内波明暗相间的条带特征不明显、噪声干扰严重，甚至出现亮、暗带缺失或畸变等；(3)实测数据是验证提取精度的最好方法，但是很多sar图像不具有同步现场观测数据，由于现场水文同步测量数据的匮乏，现有方法必须借助于历史水文数据资料，而历史水文数据资料的精确度和时空一致性并不理想。技术实现要素：针对现有技术中的上述不足之处，本发明基于sar图像相干斑的乘性特性，以及海洋内波局部低秩的先验信息，构造sar图像海洋内波参数提取的优化模型，进而实现sar图像海洋内波参数提取。本发明“一种sar图像海洋内波参数提取方法”，包括如下步骤：s1、局部低秩优化建模基于局部低秩先验，sar图像海洋内波参数提取建模过程具体包括以下4个步骤：s1-1、检测海洋内波局部低秩先验对应图像片，故而必须检测含有海洋内波的图像片。首先利用高斯滤波抑制相干斑和背景杂波：其中，y0和g分别为原始观测sar图像和高斯卷积核，y为去噪后观测sar图像。将去噪后观测sar图像y分割为若干个无重叠的图像片yi(i＝1,2,…)，并通过计算图像片中梯度指数ηi确定是否存在海洋内波。其中梯度指数ηi定义为：其中，yi为第i个图像片，积分域ωi表示yi相应的图像域。设置阈值η0，当图像片的梯度指数ηi＞η0时则认为该图像片yi中存在海洋内波，否则认为该图像片yi中不存在海洋内波。挑出含有海洋内波的图像片后，对这些含有海洋内波的图像片进行特征提取。s1-2、构造局部低秩优化模型的逼近项sar图像相干斑为乘性噪声，其统计特性满足gamma分布：其中，y和x分别为去噪后观测sar图像和相应的雷达散射系数(也即无误差图像)，n为图像等效视数，γ(·)为gamma函数，由极大似然估计理论可知，局部低秩优化模型的逼近项可等价表示为：其中，积分域ω为y相应的图像域。假设经局部低秩优化后的sar图像被分割为k个无重叠图像片则每个图像片的局部低秩优化模型的逼近项可等价表示为：逼近项(4)和(5)根据相干斑的统计特性构造，可以抑制相干斑。s1-3、构造局部低秩优化模型的局部低秩约束项约束项用于保护海洋内波的结构信息。海洋图像满足局部低秩的先验条件，因此基于sar图像目标局部低秩的先验信息，考虑图像片xi的仿射秩，用xi仿射变换后的核范数进行度量：其中为仿射变换，||·||*表示核范数，即图像表示矩阵特征值的和。仅考虑旋转变换：其中，θ为旋转变换的旋转角。此时图像片xi的仿射秩可以表示为：(8)式即为构造的局部低秩约束项。s1-4、建立局部低秩优化模型由于我们只关注海洋内波，所以仅对含有海洋内波的图像片进行参数提取，此时基于局部低秩的sar图像海洋内波参数提取模型(也即局部低秩优化模型)为：其中，||·||*为核范数，r0为秩的阈值，r0的值由图像确定。s为含有海洋内波图像片对应的指标集。(9)式的含义是：模型逼近项最小化，使得局部低秩约束项小于阈值r0。由(9)式可以看出，对每个图像片都建立了局部低秩优化模型。s2、对图像片的局部低秩优化模型进行数值求解首先进行等价形式转换。由lagrange乘子法可将带约束问题(9)转化为无约束问题：其中，λi＞0为正则化参数。利用分离变量lagrange乘子法对(10)进行等价变换。仿射变换在图像域是线性的，满足且为恒等映射。但是关于旋转角θ并不是线性的，所以需要进行局部线性化。方便起见，记则由分离变量lagrange乘子法可得(10)式的等价形式：其中，λ与μ为非负正则化参数。||·||2表示2范数，即欧拉范数。对(11)式进行数值求解。一般而言，通过最小化所有变量获得最优解运算代价过大，因此利用交替最小化方法对算法进行加速。优选的，对(11)进行交替优化求解。交替优化求解具体包括三个步骤：s2-1、固定变量u和θ，优化变量x。关于(11)式中变量x无约束优化问题等价于：由euler-larange方程可得：可通过下述迭代求解：其中xk、θk、uk表示各参数第k次的迭代结果，初始条件为x0＝y，θ0的值由图像y决定，只要使得图像y变得低秩就可以；迭代次数k为大于等于0的整数。s2-2、固定变量x和θ，优化变量u。关于(11)式中变量u无约束优化问题等价于：上式可通过迭代求解，迭代过程如下式所示：其中，svd(·)表示矩阵的svd分解，即奇异值分解，(u,∑,v)为svd分解后的三个矩阵，其中∑ij表示∑的分量。γ为非负正则化参数，γ＝λ/μ。s2-3、固定变量x和u，优化变量θ。关于(11)式中变量θ无约束优化问题等价于：迭代过程如下式所示：s2-4、当||uk+1-uk||*≤ε时，迭代终止，其中ε为迭代终止常数；输出数值求解结果：最优旋转角低秩图像片和高辐射分辨率图像片进入步骤s3；否则继续进行s2-1。s3、参数提取与表示利用步骤s2的计算结果，对sar图像海洋内波参数(即波长)进行提取。最优旋转角反映了内波运动方向，低秩图像片中可以提取出内波波长ω。具体操作包括：首先对低秩图像片的矩阵按列求和等到一个一维向量g，g的表达式为：其中j表示图像片水平方向的像素点的标号。以j为横坐标，g(j)为纵坐标，可以作出一幅波形图像，该图像即为通过能量的累加获得的海洋内波波形图像。利用高斯滤波对g进行平滑，然后记海洋内波波形图像中局部最小值相应的指标为{i1,i2,…,id}，指标的个数由海洋内波波形图像决定，则内波波长ω可近似认为是相邻指标间的平均距离：其中，δ为sar图像中相邻像素点的实际距离。优选的，上述步骤中参数的取值满足：非负正则化参数μ的取值范围为0～20，λ的取值范围为0～20；时间步长δt的取值范围0～1。与传统sar图像海洋内波参数提取方法相比，本发明提供的方法具有以下优点：(1)基于局部低秩约束的sar图像海洋内波参数提取方法可以有效克服传统方法的不足，在海洋内波对比度较弱、波形损坏的情况下依然能够实现其特征提取。(2)由于该方法采用局部低秩优化模型，充分利用sar图像相干斑的统计特性、海洋内波的局部低秩特性，可以很好地抑制sar图像相干斑、进行目标结构信息保护与目标散射特性保持，因此可以有效的实现sar图像目标识别。(3)本发明中局部低秩优化模型中的逼近项计算了局部低秩优化前后的比值，因此本发明可以抑制背景杂波。(4)模型可稳健的数值求解，便于实现。附图说明为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本发明操作流程示意图；图2本发明中局部低秩优化建模流程图；图3为本发明中交替优化求解算法流程图；图4为本发明实施例中原始观测图像；图5为本发明实施例中去噪后观测图像；图6为本发明实施例中将去噪后观测图像划分为图像片的示意图；图7为本发明实施例中挑选出含有海洋内波的图像片的示意图；图8为本发明实施例中旋转后的高斯去噪图像；图9为本发明实施例中一个包含海洋内波的图像片；图10为本发明实施例中对图9所示图像片进行局部低秩优化后输出的低秩图像片；图11为本发明实施例中对图9所示图像片进行局部低秩优化后输出的高辐射分辨率图像片；图12为本发明实施例中从图9所示图像片中提取出的海洋内波波形图像。具体实施方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明“一种sar图像海洋内波参数提取方法”，其操作流程如图1所示。具体步骤如下：s1、局部低秩优化建模基于局部低秩先验，sar图像海洋内波参数提取建模(局部低秩优化建模)过程如图2所示，具体包括以下4个步骤：s1-1、检测海洋内波局部低秩先验对应图像片，故而必须检测含有海洋内波的图像片。如图4所示为原始观测图像y0。首先利用高斯滤波抑制相干斑和背景杂波：其中，y0和g分别为原始观测sar图像和高斯卷积核，y为去噪后观测sar图像(如图5所示)。本实施例中高斯卷积核中参数的取值为：高斯滤波大小为5×5，标准差为1。如图6所示，将去噪后观测sar图像y划分为100个50×50无重叠的图像片yi(不足50×50则忽略不计)，并通过计算图像片中梯度指数确定是否存在海洋内波。其中梯度指数ηi定义为：其中，yi为第i(i＝1,2,…,100)个图像片，积分域ωi表示yi相应的图像域。设置阈值η0＝5，当图像片的梯度指数ηi＞η0时则认为该图像片yi中存在海洋内波，否则认为该图像片yi中不存在海洋内波；挑选出含有海洋内波的图像片(如图7中白色网格标示的图像片)后，只针对存在海洋内波的图像片进行特征提取。s1-2、构造局部低秩优化模型的逼近项sar图像相干斑为乘性噪声，其统计特性满足gamma分布：其中，y和x分别为去噪后观测图像和相应的雷达散射系数(也即无误差图像)，n为图像等效视数，γ(·）为gamma函数，由极大似然估计理论可知，局部低秩优化模型的可等价表示为：其中，积分域ω为y相应的图像域。假设经局部低秩优化后的sar图像被分割为k个无重叠图像片每个图像片的局部低秩优化模型的逼近项可等价表示为：逼近项(4)和(5)根据相干斑的统计特性构造，可以抑制相干斑。s1-3、构造局部低秩优化模型的局部低秩约束项约束项用于保护海洋内波的结构信息。海洋图像满足局部低秩的先验条件，因此基于sar图像目标局部低秩的先验信息，考虑图像片xi的仿射秩，用xi仿射变换后的核范数进行度量：其中为仿射变换，||·||*表示核范数，即图像表示矩阵特征值的和。对高斯去噪后观测sar图像y进行旋转变换，使图像矩阵具有初步的低秩性：其中，θ为旋转变换的旋转角。图8所示为旋转后的高斯去噪图像。此时图像片xi的仿射秩可以表示为：(8)式即为构造的局部低秩约束项。s1-4、建立局部低秩优化模型由于我们只关注海洋内波，所以仅对含有海洋内波的图像片进行参数提取，此时基于局部低秩的sar图像海洋内波参数提取模型(也即局部低秩优化模型)为：其中，||·||*为核范数，r0为秩的阈值，r0的值由图像确定。s为含有海洋内波图像片对应的指标集。(9)式的含义是：模型逼近项最小化，使得局部低秩约束项小于阈值r0。由(9)式可以看出，对每个图像片都建立了局部低秩优化模型。s2、基于局部低秩优化模型数值求解首先进行等价形式转换。由lagrange乘子法可将带约束问题(9)转化为无约束问题：其中，λi＞0为正则化参数。利用分离变量lagrange乘子法对(10)进行等价变换。仿射变换在图像域是线性的，满足且为恒等映射。但是关于旋转角θ并不是线性的，所以需要进行局部线性化。方便起见，记则由分离变量lagrange乘子法可得(10)式的等价形式：其中，λ与μ为非负正则化参数。||·||2表示2范数，即欧拉范数。然后对(11)进行交替优化求解。一般而言，通过最小化所有变量获得最优解运算代价过大，因此利用交替最小化方法对算法进行加速。交替优化求解流程如图3所示，具体包括三个步骤：s2-1、固定变量u和θ，优化变量x。关于变量x无约束优化问题等价于：由euler-larange方程可得：可通过下述迭代求解：其中初始条件为x0＝y，θ0的值由图像y决定，只要使得图像y变得低秩就可以；迭代次数k为大于等于1的整数。s2-2、固定变量x和θ，优化变量u。关于(11)式中变量u无约束优化问题等价于：上式可通过迭代求解，迭代过程如下式所示：其中，svd(·)表示矩阵的svd分解，即奇异值分解，(u,∑,v)为svd分解后的三个矩阵，其中∑ij表示∑的分量。γ为非负正则化参数，γ＝λ/μ。s2-3、固定变量x和u，优化变量θ。关于(11)式中变量θ无约束优化问题等价于：迭代过程如下式所示：s2-4、当||uk+1-uk||*≤ε时，迭代终止，其中ε为迭代终止常数；输出数值求解结果：最优旋转角低秩图像片和高辐射分辨率图像片进入步骤s3；否则继续进行s2-1。本实施例中，对图7所示的含海洋内波的图像片中选取一个进行处理，所选取出来的图像片如图9所示。所设置参数如表1所示：表1数值求解参数设置旋转角θ非负正则化参数λ非负正则化参数μ非负正则化参数γ时间步长δt0°0.70.51.40.01通过代入上述参数数值计算得到该图像片的最优旋转角低秩图像片(如图10所示)，高辐射分辨率图像片(如图11所示)。s3、参数提取与表示利用步骤s2的计算结果，对sar图像海洋内波参数(即波长)进行提取。最优旋转角反映了内波运动方向，低秩图像片中可以提取出内波波长ω。具体操作包括：本实施例中首先对图10所示的低秩图像片的矩阵按列和等到一个一维向量g，g的表达式为：并利用高斯滤波对g进行平滑，可以得到海洋内波波形图像，如图12所示，其中横坐标j的最大值为50(即该图像片水平方向的像素点个数为50)，纵坐标表示g的数值。由图12可以判断出该图像片海洋内波波形中局部最小值相应的指标为{i1,i2}，其中i1＝19,i2＝35，则内波波长ω可近似认为是相邻指标间的平均距离：其中，δ＝25米为sar图像中相邻像素点的实际距离。当前第1页12