Use of Representative Pixels of Motion for Wind Vector Height Assignment of Semi-transparent Clouds
-
摘要: 半透明云风矢量高度指定是卫星风矢量算法的重要部分,需要使用来自半透明云体的辐射和云下背景辐射两个变量。为了更精确地获得来自半透明云体的辐射,采用像元对追踪相关系数的贡献和红外亮温两个变量进行统计分析,将追踪图像块中的像元分为3个部分:小贡献像元、暖区段大贡献像元和冷区段大贡献像元;对暖区段、冷区段内大贡献像元的特征进行考察表明,冷区段的大贡献像元更能代表追踪图像块运动,称为代表运动像元。选择代表运动像元参与计算来自半透明云体的辐射,对半透明云风矢量高度指定算法进行改进,在该算法中分别使用和不使用代表运动像元,计算FY-2气象卫星风矢量,将结果与欧洲中期天气预报中心(ECMWF)分析场进行对比表明,在半透明云高度指定算法中使用代表运动像元,FY-2气象卫星风矢量误差明显降低。
-
关键词:
- 风矢量;
- 半透明云风矢量高度指定;
- 代表运动像元
Abstract: Height assignment of semi-transparent clouds needs two variables:Radiation from semitransparent clouds themselves and radiation from the background underneath semitransparent clouds. Generally, the former radiation is calculated using algorithms based on all pixels in tracer boxes. A new algorithm, however, is proposed considering only pixels with great contribution to correlation coefficients during wind vector tracing stage, and only allow these pixels to join the calculation. Scatter diagrams are made with infrared brightness temperature against their contribution to the feature tracking correlation coefficients. Scatter diagrams show a "C" shape pattern and can be used to recognize representative pixels of motion. Pixels with relative less contribution to correlation coefficients are at the left side of the "C" shape with infrared brightness temperature in the middle range. The other pixels with relative more contributions are naturally separated into two branches with warmer and colder infrared brightness temperature respectively. When tracer boxes are partly covered by clouds, the warm branch pixel brightness temperature is with typical ground feature at both the scope and the variability. And when tracer boxes are totally covered by clouds, warm branch pixels are in the slot of clouds with typical cloud features among cloud cells. Thus, no matter tracer boxes are partly covered or totally covered by clouds, cold branch pixels with greater contribution to correlation coefficients are better at representing the cloud motion and can be named as pixels with contribution to the motion. The algorithm using different pixel's contribution to the feature tracking improves height assignment quality by effectively removing rapid changing pixels in tracer boxes when estimating the cloud height. Using representative pixels of motion at the estimation of radiation from semi-transparent clouds themselves, FY-2 wind vector height assignment quality is improved.-
Key words:
- wind vector;
- height assignment;
- pixels representative of motion
-
图 2 追踪区部分被云覆盖案例 (a) 追踪图像块位于30.94°N, 120.44°E, (b) 追踪图像块位于15.86°N, 118.81°E, (c) 追踪图像块位于33.41°N, 123.83°E
Fig. 2 Typical examples for the tracer partly covered by clouds (a) tracer box located at 30.94°N, 120.44°E, (b) tracer box located at 15.86°N, 118.81°E, (c) tracer box located at 33.41°N, 123.83°E
图 3 追踪区全部被云覆盖典型案例 (a) 追踪图像块位于34.56°N, 112.82°E, (b) 追踪图像块位于6.27°N, 139.42°E, (c) 追踪图像块位于5.33°N, 93.30°E,(d) 追踪图像块位于2.42°N, 90.35°E, (e) 追踪图像块位于29.27°S, 97.89°E
Fig. 3 Typical examples for the tracer totally covered by cloud (a) tracer box located at 34.56°N, 112.82°E, (b) tracer box located at 6.27°N, 139.42°E, (c) tracer box located at 5.33°N, 93.30°E, (d) tracer box located at 2.42°N, 90.35°E, (e) tracer box located at 29.27°S, 97.89°E
图 6 暖区段和冷区段像元数之比直方图
(a) 暖区段和冷区段总贡献之比小于0.3,(b) 暖区段和冷区段总贡献之比为0.8~1.2,(c) 暖区段和冷区段总贡献之比大于2
Fig. 6 Histogram of pixel number ratio between warm and cold segments
(a) total contribution ratio between warm and cold segments less than 0.3, (b) total contribution ratio between warm and cold segments in 0.8-1.2, (c) total contribution ratio between warm and cold segments greater than 2
图 7 追踪图像块中最高红外亮温的直方图
(a) 暖区段和冷区段总贡献之比小于0.3,(b) 暖区段和冷区段总贡献之比为0.8~1.2,(c) 暖区段和冷区段总贡献之比大于2
Fig. 7 Histogram of IR maximum temperature in tracer boxes
(a) total contribution ratio between warm and cold segments less than 0.3, (b) total contribution ratio between warm and cold segments in 0.8-1.2, (c) total contribution ratio between warm and cold segments greater than 2
表 1 半透明云高度指定算法中是否使用代表运动像元FY-2E卫星红外高层风风速与ECMWF分析场比较
Table 1 Bias and standard deviations between FY-2E wind vectors and ECMWF analyses for height assignment algorithms whether using contribution to the motion
统计量 未使用代表运动像元 使用代表运动像元 偏差/(m·s-1) -4.239 -2.103 标准差/(m·s-1) 7.787 7.436 风矢量数 13895 13536 -
[1] Xu Jianmin, Zhang Qisong.Calculation of Cloud Motion Wind with GMS-5 Images in China//Proceedings of the Third International Winds Workshop, EUMETSAT Publication.EUM P18, 1996:45-52. [2] Xu Jianmin, Zhang Qisong, Fang Xiang, et al.Cloud Motion Winds from FY-2 and GMS-5 Meteorological Satellites//Proceedings of the 4th International Winds Workshop, EUMETSAT Publication.EUM P24, 1998:41-48. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.475.9694&rep=rep1&type=pdf [3] Xu Jianmin, Zhang Qisong.Recent Works Aimed at Operational FY-2C AMVS//Proceedings of the 7th International Winds, 2004. [4] Xu Jianmin, Zhang Qisong, Zhang Xiaohu, et al.Status of operational AMVs from FY-2C//Proceedings of the 8th International Winds Workshop, 2006. [5] Zhang Xiaohu, Xu Jianmin.Two Applications of Improvements for AMVs of NSMC/CMA-Re-navigation Based on Full Earth Disk Image & Calibration of Radiation Transfer Using NWP Data. The 10th International Winds Workshop, 2010. [6] 许健民, 张其松, 王大昌.云迹风计算中的两个几何问题.应用气象学报, 1997, 8(1):11-18. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=19970102&flag=1 [7] 许健民, 张其松.卫星风推导和应用综述.应用气象学报, 2006, 17(5):574-582. doi: 10.11898/1001-7313.20060515 [8] 薛谌彬, 龚建东, 薛纪善, 等.FY-2E卫星云导风定高误差及在同化中的应用.应用气象学报, 2011, 22(6):681-690. doi: 10.11898/1001-7313.20110605 [9] 王金成, 陆慧娟, 韩威, 等.GRAPES全球三维变分同化业务系统性能.应用气象学报, 2017, 28(1):11-24. doi: 10.11898/1001-7313.20170102 [10] 黄丽萍, 陈德辉, 邓莲堂, 等.GRAPES_Meso V4.0主要技术改进和预报效果检验.应用气象学报, 2017, 28(1):15-37. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20170103&flag=1 [11] 朱立娟, 龚建东, 黄丽萍, 等.GRAPES三维云初始场形成及在短临预报中的应用.应用气象学报, 2017, 28(1):38-51. doi: 10.11898/1001-7313.20170104 [12] Menzel W P, Smith W L, Stewart T.Improved cloud motion wind vector and altitude assignment using VAS.J Climate Appl Meteor, 1983, 22:377-384. doi: 10.1175/1520-0450(1983)022<0377:ICMWVA>2.0.CO;2 [13] Nieman N J, Schmetz J, Menzel W P.A comparison of several techniques to assign heights to cloud tracers.J Appl Meteor, 1993, 32:1559-1568. doi: 10.1175/1520-0450(1993)032<1559:ACOSTT>2.0.CO;2 [14] 许健民, 张其松, 方翔.用红外和水汽两个通道的卫星测值指定云迹风的高度.气象学报, 1997, 55(4):408-417. doi: 10.11676/qxxb1997.041 [15] Szejwach G.Determination of semitransparent cirrus cloud temperature from infrared radiances:Application to Meteosat.J Appl Meteor, 1982, 21(3):384-393. doi: 10.1175/1520-0450(1982)021<0384:DOSTCC>2.0.CO;2 [16] Menzel W P.Cloud tracking with satellite imagery:From the pioneering work of Ted Fujita to the present.Bull Amer Meteor Soc, 2001, 82:33-47. doi: 10.1175/1520-0477(2001)082<0033:CTWSIF>2.3.CO;2 [17] Velden C.Recent innovations in deriving tropospheric winds from meteorological satellites.Bull Amer Meteor Soc, 2005, 86:205-223. doi: 10.1175/BAMS-86-2-205 [18] Xu J, Holmlund K, Zhang Q, et al.Comparison of two schemes for derivation of atmospheric motion vectors.J Geophys Res, 2002, 107(D14):ACL 4-1-ACL 4-15. [19] Büche G.Water vapor structure displacements from cloud-free Meteosat scenes and their interpretation for the wind field.J Appl Meteor, 2006, 45:556-575. doi: 10.1175/JAM2343.1 [20] Borde R.A Direct Link Between Feature Tracking and Height Assignment of Operational Atmospheric Motion Vectors.Proc the Ninth Int Winds Workshop, 2008. [21] Oyama R, Borde R, Schmetz J, et al.Development of Height Assignment Directly Linked to Feature Tracking at JMA//Proceedings of 9IWW, 2008. [22] 张晓虎, 张其松, 许健民.半透明云风矢量高度算法中云下背景辐射的估计.应用气象学报, 2017, 28(3):283-291. doi: 10.11898/1001-7313.20170303