基于X波段雷达图像反演近岸水深技术的创新论文.docx

基于X波段雷达图像反演近岸水深技术的创新论文0 引 言对近岸浅水区域的海底地形的掌握对于近岸活动比如捕鱼，码头建设，铺设石油管道以及形态动力学的研究等都有着至关重要的作用传统获取水深的方法包括船载声呐探测、机载激光测探、潜水器测量、超光谱图像测量等但是它们普遍具有成本高、效率低的缺点，并且测量精度会受海水清澈度的影响为了克服传统测量方法所存在的不足，基于X波段雷达图像序列反演水深的方法得到发展在浅水区域，局部的海底地形对表面波的传播有着重要的影响当波移动至浅水区域，波的周期不发生变化，但是波的传播速度会发生变化，继而波长减小，波数增加表面流的存在也影响表面波的传播，因此水深场以及表面流速场的反演方法都是基于这种传播变化——在物理学中被描述为表面波的色散关系1998年 Paul Bell运用连续的雷达图像序列之间的互相关性推导出了空间变化的表面波速【1】，同时利用从浮标中获取的频率信息，通过运用线性重力波的色散关系计算出了空间变化的水深，但是没有考虑海流的存在Hessner等人运用一维FFT变换实现了对图像序列的频率分解【2】，某一固定频率的波所对应的波长通过确定局部空间的相位梯度计算得到。

但是此方法的局限性在于它不能运用到包含同一频率但不同传播方向的波的波场，同时也没有考虑到海流的存在对时间序列的雷达图像进行3?D FFT变化，并取模的平方得到三维图像谱，由于波数和频率被色散关系联系在一起，因此线性表面波的信号应该很好地分布于其确定的三维形状上色散关系的形状取决于水深和表面流速，因此通过拟合理论的色散关系和三维图像谱的坐标分布便可反演出大的空间范围内的平均水深以及流速但是此方法中的3?D FFT是针对全局范围的算子，因此假定了波场的均匀性以及稳定性如果在深水区存在变化的流速或者浅水区存在变化的水深，波的折射将会产生，波场变成了非均匀场，以上方法不再适用，因此需要在局部空间范围内对波参数进行分析自1999年以来，Seemann等人针对非均匀波场做了一系列研究，推导出了局部三维图像谱，同时反演出了局部范围内的水深以及流速本文将利用模拟的X波段雷达图像展开近岸浅水区域的水深的反演工作，该工作考虑到了表面波场的非均匀性，因此采用了局部反演算法，反演出了局部的水深值1 色散关系与水深以及流速的关系色散关系描述了波数和角频率之间的动力学关系，正常的色散关系适用于海表面重力波，线性色散关系可表示为： 式中：表示为重力加速度；为水深；为近表面流速。

在式其中三维谱的谱分辨率为： 图4 一组方向滤波器中的前四个 单波空间场包含了幅值及相位模式信息： 与单波空间场对应的梯度图像： 其中代表复数值的波数向量，其实部代表局部的波数值局部区域的大小选为88个像素点，因此要得到局部区域的波数，需要分析局部点所包含的所有像素点 位于色散关系滤波器带宽内的背景噪声重新分布在了单波波数场中，因此为了消除噪声的影响，运用方差最小拟合法得到复数值的波数向量 其中向量通过行扫描局部区域内的单波空间场及其梯度图像获得，向量是向量的共轭向量 式中：表示谱坐标集中元素的个数，通过设置阈值从局部能量谱中选取出谱坐标集：表示能量阈值加权方差函数是一个非线性函数，含有三个变量，求该函数最小值属于优化问题，本文采用拟牛顿法搜索最小值，并得到局部的水深及流速4 数据处理结果模拟雷达图像时输入的非等值水深场如图5（a）所示，每88个像素点设置一个水深值，为减少模拟时的计算量，水深值只沿一维变化反演得到的水深场如图5（b）所示，反演时选择的局部区域的大小为88个像素点反演的水深值与输入的水深值吻合较好，平均误差约为2%，相比于过去的均匀场水深反演方法，该反演方法可将水深值的分辨率缩小到88个像素点。

5 结 语利用X波段雷达图像可提取出重要的海态信息，比如水深、流速等等均匀场的水深及流速的反演方法已相对成熟，本文的工作是针对非均匀场反演浅水水深值由于实际的海况比较复杂，并且还没有得到可以用于比测的实际水深值，本文采用数值模拟的方法，通过输入非等值的水深仿真出非均匀波场及其雷达图像利用仿真的雷达图像反演出局部水深值，并与输入的水深值进行对比，结果吻合较好，对利用实际的雷达图像反演非均匀场的水深具有重要的指导意义本文的工作是基于岸基X波段雷达，对于船基X波段雷达来说，还要考虑运动补偿等因素，并且实际海况复杂多变，水深的反演过程有待进一步分析研究图5 输入的水深场与反演得到的水深场对比图参考文献. Coast. Eng.， 1999， 37（3）： 513?527.， 1999： 16?20.【3】 SENET C， SEEMANN J， ZIEMER F. An iterative technique to determine the near surface current velocity from time series of sea surface images // Proceedings of Oceans MTS/IEEE Conference?500 Years of Ocean Exploration. ： IEEE，1997： 66?72.IEEE Transactions on Geoscience and Remote Sensing， 2002， 40（12）： 2606?2618.. Procee?dings of SPIE 1999， 3808： 536?546.： IEEE， 2000： 1898?1900.【7】 SEEMANN J， SENET C， ZIEMER F. Local analysis of inhomogeneous sea surfaces in coastal waters using nautical radar image sequences. Berlin， Germany： Springer?Verlag， 2000： 179?186.. Proceedings of SPIE， 1999， 3808： 148?158.： IEEE， 2000： 1769?1774.。