摘要(英) |
The distance between the Chunge-Li VHF radar station and Taiwan Taoyuan International Airport is about 12km. In the meanwhile, the Chunge-Li VHF radar station locates in the lane of the taking-off or landing airplane. When the missions of observing the atmosphere are performed by Chunge-Li VHF radar, the echo signal is always polluted by the echo power of the crossing-radar-beam airplane, so the quality of the atmospheric data would be influenced and low quality. However, in the paper, [洪萱芸, 2020], an algorithm was proposed firstly. The algorithm can filter out the signal of airplane by Hilbert-Huang Transform (HHT). The result of the algorithm is better than the result by using Wavelet Transform. The algorithm greatly improves the quality of the observation data. However, the condition is only for the vertical beam and the result of the algorithm isn’t applied to be calculated the 3-D wind field by Velocity-Azimuth Display (VAD) technology、raindrop size distribution …etc. Thus, this study attempts to find the adapting parameter in the algorithm to apply in the observation of the oblique beam. After the aircraft echo is filtered out, the result will try to calculate the 3-D wind field by Velocity-Azimuth Display (VAD) technology. Finally, the wind field would be compared with Central Weather Bureau wind profiler(449MHz) result, and the application would be further discussed.
Through the algorithm, the high frequency component of aircraft echo would be filtered out. However, if the aircraft echo contained the low-frequency component or the problem of mixing wave from doing Empirical Mode Decomposition (EMD), there would be anomalous strong signal at low-frequency of the spectrum. With these problems, the algorithm should be modified. The way is that, with the DC terms of the real part (In-Phase, I) and imaginary part (Quadrature, Q), the amplitude at the time of aircraft echo on the I and Q is replaced. The time of aircraft echo is still detected by the HHT. This algorithm is named “DC term replacement”. Finally, the root-mean-square error(RMSE) and correlation coefficient between the 3-D wind field by “filtering out by HHT” or “DC term replacing” and the wind field of the Central Weather Bureau (CWB) wind profiler can be calculated. In the meanwhile, two indexes between the wind field before filtered out and the wind field of CWB wind profiler can also be calculated. A comparison with two parameters of three different conditions would be made. Not only “filtering out” but also “DC term replacing” result can be more greatly improved than the result before filtered out. To the other condition, the comparison between “no aircraft”、“filtering out” and “DC term replacing” should be made. The results of “filtering out” and “DC term replacing” still need to be improved by comparing with the result of “no aircraft”. To the results of “filtering out” and “DC term replacing”, it is comparable each other. But, visually, the spectrum of “DC term replacement” is more reasonable. Eventually, DC term replacement would be recommended.
The distance between the Chunge-Li VHF radar station and Taiwan Taoyuan International Airport is about 12km. In the meanwhile, the Chunge-Li VHF radar station locates in the lane of the taking-off or landing airplane. When the missions of observing the atmosphere are performed by Chunge-Li VHF radar, the echo signal is always polluted by the echo power of the crossing-radar-beam airplane, so the quality of the atmospheric data would be influenced and low quality. However, in the paper, [洪萱芸, 2020], an algorithm was proposed firstly. The algorithm can filter out the signal of airplane by Hilbert-Huang Transform (HHT). The result of the algorithm is better than the result by using Wavelet Transform. The algorithm greatly improves the quality of the observation data. However, the condition is only for the vertical beam and the result of the algorithm isn’t applied to be calculated the 3-D wind field by Velocity-Azimuth Display (VAD) technology、raindrop size distribution …etc. Thus, this study attempts to find the adapting parameter in the algorithm to apply in the observation of the oblique beam. After the aircraft echo is filtered out, the result will try to calculate the 3-D wind field by Velocity-Azimuth Display (VAD) technology. Finally, the wind field would be compared with Central Weather Bureau wind profiler(449MHz) result, and the application would be further discussed.
Through the algorithm, the high frequency component of aircraft echo would be filtered out. However, if the aircraft echo contained the low-frequency component or the problem of mixing wave from doing Empirical Mode Decomposition (EMD), there would be anomalous strong signal at low-frequency of the spectrum. With these problems, the algorithm should be modified. The way is that, with the DC terms of the real part (In-Phase, I) and imaginary part (Quadrature, Q), the amplitude at the time of aircraft echo on the I and Q is replaced. The time of aircraft echo is still detected by the HHT. This algorithm is named “DC term replacement”. Finally, the root-mean-square error(RMSE) and correlation coefficient between the 3-D wind field by “filtering out by HHT” or “DC term replacing” and the wind field of the Central Weather Bureau (CWB) wind profiler can be calculated. In the meanwhile, two indexes between the wind field before filtered out and the wind field of CWB wind profiler can also be calculated. A comparison with two parameters of three different conditions would be made. Not only “filtering out” but also “DC term replacing” result can be more greatly improved than the result before filtered out. To the other condition, the comparison between “no aircraft”、“filtering out” and “DC term replacing” should be made. The results of “filtering out” and “DC term replacing” still need to be improved by comparing with the result of “no aircraft”. To the results of “filtering out” and “DC term replacing”, it is comparable each other. But, visually, the spectrum of “DC term replacement” is more reasonable. Eventually, DC term replacement would be recommended. |
參考文獻 |
1. Balsley, B. B. and K. S. Gage, 1980: “The MST radar technique: Potential for middle atmospheric studies”, Pure and Applied Geophys., Vol. 118, P. 452-493.
2. Gage, K. S., W. L. Ecklund, and B. B. Balsley, 1985: “A modified Fresnel scattering model for the parameterization of Fresnel returns”, Radio Science, 20, 1493-1501.
3. Wait, J. R., 1962: “Electromagnetic waves in stratified media”, Pergamon Press, New York, pp 305.
4. Yeh, K. C., and C. H. Liu, 1974: “Theory of ionosphere waves, Academic Press, New York”, pp 464.
5. Y. H. Chu, T. S. Hsu, L. H. Chen, J. K. Chao, C. H. Liu, and J. Rottger, 1990:”A Study of The Characteristics of VHF Radar Echo Power in The Taiwan Area”, Radio Sci., Vol. 25, pp. 527-538
6. Ottersten, H.,1969: “Radar backscattering from the turbulent clear atmosphere”, Radio Sci., 4, 1251-1255.
7. Tatarskii, V. I., 1961: “Wave Propagation in a Turbulent Medium”, McGraw-Hill, New York, 285pp.
8. Browning, K. A. and Wexler, R., 1968: “The determination of kinematic properties of a wind field using Doppler radar”, J. Appl. Meteorol., 7, 105–113.
9. Woodman, R.F. and W.E. Guillen, 1974: “Radar Observations of Winds and Turbulence in the Stratosphere and Mesosphere”, J. Atmos. Sci., Vol.31, p.495-505
10.Koichiro Wakasugi, Shoichiro Fukao, Susumu Kato, Akiyoshi Mizutani, and Masaru Mastuo, 1985:”Air and Precipitation Particle Motions within a Cold Front Measured by the MU VHF Radar”, Radio Sci. Vol. 20, Num. 6, pp. 1233-1240
11.Koichiro Wakasugi, Akiyoshi Mizutani, and Masaru Mastuo, 1986:”A Direct Method for Deriving Drop-Size Distribution and Vertical Air Velocities from VHF Doppler Radar Spectra”, Journal of Atmospheric and Oceanic Technology, Vol.3, pp. 623-629
12. Y.H. Chu and J. K. Chao, 1990:”Aspect Sensitivity at Tropospheric Heights Measured with Vertically Pointed Beam of the Chung-Li VHF Radar”, Radio Sci., Vol. 25, Num. 4, pp. 539-550
13. Norden E. Huang et al, 1998:“The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis”, Royal Society, Volume 454, Issue 1971
14. Norden E. Huang et al, 2003:” A confidence limit for the empirical mode decomposition and Hilbert spectral analysis”, Royal Society, Volume 459, Issue 2037
15. Jidong Gao, Kelvin K. Droegemeier, Jiandong Gong, and Qin Xu, 2004:”A Method for Retrieving Mean Horizontal Wind Profiles from Single-Doppler Radar Observations Contaminated by Aliasing”, GAO ET AL, pp. 1399-1409
16. Norden E. Huang etc, 2009:” On Instantaneous Frequency “, World Scientific, Vol. 1, No. 2, pp. 177-229
17. Zhaohua Wu and Norden E. Huang, 2009:” Ensemble Empirical Mode Decomposition: A Noise-Assisted Data Analysis Method”, World Scientific, Vol. 1, No. 1, pp. 1-41
18. J.-C. BOISSE, V. KLAUS, AND J.-P. AUBAGNAC, 1998:” A Wavelet Transform Technique for Removing Airplane Echoes from ST Radar Signals”, Journal of Atmospheric and Oceanic Technology, Vol.16, pp. 334-346
19. B. L. Weber, D. B. Wuertz, D. C. Welsh, 1992:” Quality Controls for Profiler Measurements of Winds and RASS Temperatures”, Journal of Atmospheric and Oceanic Technology, Vol.10, pp. 452-464
20. 陳孟遠,利用中壢特高頻雷達研究對流降水系統雨滴粒徑與速度之關係,碩士論文,國立中央大學,1990。
21. 陳昭宇,多頻段剖風儀雷達觀測結果之比對與分析,碩士論文,國立中央大學,2019。
22. 林廷翰,中壢特高頻雷達系統初始相位偏差估計與應用,博士論文,國立中央大學,2019。
23. 洪萱芸,以希爾伯特-黃轉換辨識並濾除特高頻雷達飛機雜波,碩士論文,國立中央大學,2020。
ICAO. “Technical Provisions for Mode S Services and Extended Squitter.” Doc/ International Civil Aviation Organization ; 9871-AN/464. ICAO, Montréal, QC,Canada, 1 edition, 2008. |