Ocean Surface Non-cyclone Wind Block Ambiguity Removal Algorithm for Scatterometer
-
摘要: 利用最大似然法 (MLE) 对散射计数据反演得出的风矢量,一般存在多个模糊解,故需采用圆中数滤波法进行模糊去除。但传统圆中数滤波法难以解决风场反演中块状模糊问题。该文根据散射计非气旋性第1风场的空间分布特性,归纳出一套适合散射计的加强型圆中数滤波块状模糊去除方法,并探讨了其适应性。该方法定义简单,计算量小,且易收敛。利用欧洲EUMETSAT提供的部分ASCAT (Advanced SCATterometer) L2原始数据对该方法进行验证结果表明,对于非气旋风场分布情况,该方法在利用其他工具去除台风气旋覆盖区域后,能有效解决非气旋区域风场块状模糊问题。Abstract: The Maximum Likelihood Estimation (MLE) algorithm for scatterometer wind vector retrieval generates several wind vector ambiguities, so a circle median filter is needed to perform the ambiguity removal. But the traditional circle median filtering method can hardly solve the block ambiguity problem. According to the spatial distribution characteristic of the most likely ambiguities in each non-cyclone wind vector cell, a new enhanced circle median filtering method for block ambiguity removal is derived and discussed theoretically, with experiments carried out to check its adaptability. This method features simple definition, low computation and easiness to converge. Using some L2 raw data from EUMETSAT to validate the method, the results indicate that under non-cyclone wind distribution condition the new method is effective in resolving the problem of block ambiguity after eliminating the cyclone wind field with other reference data.The core of the enhanced circular median filter algorithm is to initialize the non-cyclone characteristics first, which can effectively solve the problem of massive fuzzy. Thus the data which may cause circular median filtering failure and the data that affect the neighborhood will be corrected. Then divide the two-dimensional space into M rows and N columns, then calculate the wind field in the open window. The wind vector in the center of a window is solved by selecting an alternative from corresponding fuzzy solutions. And then do this with the next location iteratively, until the wind field does not change or until the times of iterations reaches a preset maximum number. Finally, defective value in the wind field is smoothed.Compared with traditional circle median filtering method, this approach is better in several ways. First, by initializing the first wind field, the fuzziness of the second wind field is reduced; while with the traditional method the block fuzzy cannot be removed. Second, the calculation process is simple and need no statistical circular histogram, nor do they need to calculate the mean. Third, the definition of circular median is only one rather than getting multiple solutions. And last, the calculation will not be interrupted by the narrow wind element values in the boundary region.The enhanced circular median filtering method based on non-cyclone wind field vector distribution characteristics can overcome the harsh conditions of the traditional method (such as wind field must be randomly distributed, non-block fuzzy), extracting the true wind vector solutions to overcome the fuzzy block. In order to apply this method widely, a crucial issue is to determine and eliminate the coverage typhoon cloud by operational satellite equipment properly. The method provides a new idea for exacting data of non-cyclone wind field on the ocean surface.
-
Key words:
- ASCAT;
- non-cyclone wind;
- enhanced circle median filter;
- block ambiguity
-
表 1 本研究结果与EUMETSAT L2风场数据偏差和绝对偏差统计
Table 1 Statistics on biases between filtered wind vector field and L2 wind vector field from EUMETSAT
统计量 偏差 绝对偏差 风速/(m·s-1) 风向/(°) 风速/(m·s-1) 风向/(°) 最小值 -10.09 -48.4 0 0 最大值 3.31 57.1 10.09 57.1 平均值 -1.5858 1.0424 1.859 8.426 方差 3.3437 155.17 2.4035 57.1 表 2 本研究结果与EUMETSAT的L2风场数据偏差和绝对偏差统计
Table 2 Statistics on biases between filtered wind vector field and L2 wind vector field from EUMETSAT
统计量 偏差 绝对偏差 风速/(m·s-1) 风向/(°) 风速/(m·s-1) 风向/(°) 最小值 -8.09 -30.4 0 0 最大值 4.17 51.7 8.03 71.1 平均值 -1.932 2.1457 1.651 10.317 方差 5.1473 175.14 3.2156 78.42 -
[1] 魏应植, 汤达章, 许健民, 等.多普勒雷达探测"艾利"台风风场不对称结构.应用气象学报, 2007, 18(3):285-294. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20070350&flag=1 [2] 魏应植, 吴陈锋, 苏卫东.利用多普勒雷达径向速度提取台风环境风场信息.应用气象学报, 2010, 21(3):307-316. doi: 10.11898/1001-7313.20100306 [3] 钟刘军, 阮征, 葛润生, 等.风廓线雷达回波信号强度定标方法.应用气象学报, 2010, 21(5):598-605. doi: 10.11898/1001-7313.20100509 [4] 肖卫华, 符养, 高太长, 等.利用折射指数推算大气纬圈平均风场方法.应用气象学报, 2011, 22(3):346-355. doi: 10.11898/1001-7313.20110311 [5] 薛谌彬, 龚建东, 薛纪善, 等.FY-2E卫星云导风定高误差及在同化中的应用.应用气象学报, 2011, 22(6):681-690. doi: 10.11898/1001-7313.20110605 [6] [2011-06-09].http://www.esa.int/esaME/ascat.html. [7] Figa-Saldaña J, Wilson J J W, Attema E, et al. The advanced scatterometer (ASCAT) on the meteorological operational (MetOp) platform: A follow on for European wind scatterometers. Can J Remote Sensing, 2002, 28(3): 404-412. doi: 10.5589/m02-035 [8] 孙强, 孙瀛.我国首台天基微波散射计测量海面风场的应用研究.海洋技术, 2009, 28(1):70-74. http://www.cnki.com.cn/Article/CJFDTOTAL-HYJS200901018.htm [9] 李燕初, 孙瀛, 林明森, 等.用圆中数滤波器排除卫星散射计风场反演中的风向模糊.台湾海峡, 1999(1):42-48. http://www.cnki.com.cn/Article/CJFDTOTAL-TWHX199901006.htm [10] Schroeder L C, Grantham W L, Bracalente E M, et al. Removal of ambiguous wind directions for a Ku-band wind scatterometer using three different azimuth angles. IEEE Transactions on Geoscience and Remote Sensing, 1985(2):91-100. http://ieeexplore.ieee.org/document/4072256/ [11] Ebuchi N. Evaluation of Wind Vectors Observed by QuikSCAT/SeaWinds Using Ocean Buoy Data. Geoscience and Remote Sensing Symposium, 2001: 1082-1085. doi: 10.1175/1520-0426%282002%29019<2049%3AEOWVOB>2.0.CO%3B2 [12] Shaffer S J, Dunbar R S. A median-filter-based ambiguity removal algorithm for NSCAT. IEEE Transactions on Geoscience and Remote Sensing, 1991, 29(1):167-174. doi: 10.1109/36.103307 [13] Freilich M H, SeaWinds P I. Algorithm Theoretical Basis Document. 1999. [14] 解学通, 方裕, 陈晓翔, 等.基于最大似然估计的海面风场反演算法研究.地理与地理信息科学, 2005, 21 (1):30-33. http://www.cnki.com.cn/Article/CJFDTOTAL-DLGT200501009.htm [15] Chi C Y, Li F K. A comparative study of several wind estimation algorithms for spaceborne scatterometers. IEEE Transactions on Geoscience and Remote Sensing, 1988(2): 115-121. http://ieeexplore.ieee.org/document/3011/ [16] Howard S. A circlar median filter approach for resolving directional ambiguities in wind fields retrieved from spaceborne scatterometer data. J G R, 1990, 95(4):5291-5303. doi: 10.1029/JC095iC04p05291/abstract [17] Scott J, Shaffer R D, Vincent H S, et al. A median-filter-based ambiguity removal algorithm for NSCAT. IEEE Transaction on Geoscience and Remote Sensing, 1991, 29(1):167-174. doi: 10.1109/36.103307 [18] 解学通, 方裕, 陈克海, 等. SeaWind散射计海面风场模糊去除方法研究.北京大学学报:自然科学版, 2005, 41(6): 882-889. http://www.cnki.com.cn/Article/CJFDTOTAL-RDQX201005015.htm