Nghiên cứu hoàn thiện hệ thông lida quan trắc son khí

VIỆN NGHIÊN CỨU ĐIỆN TỬ, TIN HỌC, TỰ ĐỘNG HÓA U DUY TU N SOL Ĩ - 2021 VIỆN NGHIÊN CỨU ĐIỆN TỬ, TIN HỌC, TỰ ĐỘNG HÓA U DUY TU N - 2021 I CA Đ AN s u n nv t n tl n n v n usn xn m o n yl n trìn n n kết n n u trìn b y tron lu n n l trun t từn ể bảo v v n o n t u nt u– nn n u n nl m o uy n u n n usn xn lu n n ã m ơn, ỉ rõ n uồn ố m o n rằn m t n t n trí LIDAR nt , n , n v n n t m u n n u ,k v mìn , h quan s úp ỡ o v t n tron lu n n n y u , n y tháng năm 2021 ả lu n n D i T n IC r t n, n ỡ, usn xn u n Anh, n v nv k í N n l n l ảm ơn ến TS t y ã n n n uy n u n uy n m n, u s n tron suốt qu trìn t úp n lu n n t ến sĩ n ã usn xnb yt l n v n, k í usn t l n nn n n u ảm ơn t uy n u s n tron suốt t nt ,Tn ,T ế ruy n nt n y nn từ năm n ến tron n u sn nn n p v t n ã s n x n ảm ơn t u n nt ,Tn ,T on n usn n úp ỡ tốt n ất ể n n u sn n , n v to n t ể , n h n n tìn ảm n n b n ín t n nt n u, n t p tron n qu n nt n ảm ơn! , n y … tháng … năm 2021 ả lu n n D ii T C C Danh mụ k u, viết tắt v t ến n v n mụ b ểu xi n mở u C ơn n qu n v t ốn r qu n trắ sol k í t vấn 1.1 p ơn p p o sol k í tron k í 11 1.1.1 ơn p p o sol k í tr c tiếp 12 1.1.2 ơn p p o sol k í n 1.1.3 ụn n t ếp ( r tr n t ế v t ) 13 t m 16 r o sol k í 18 1.2 Các nghiên c u thiết kế, chế t o h 1.2.1 Trên gi i: 18 1.2.2 Ở t m 24 ơn p p x n cao l p biên khí số li u Lidar 27 1.3.1 Trên gi i 27 1.3.2 Ở Vi t Nam 29 ết lu n ơn 30 ơn ơn p p qu n trắ sol k í bằn r 32 t vấn 32 2.1 Nguyên lý ho t ng c r o c thơng số khí 39 2.2 ơn trìn Lidar 44 2.2.1 số K 48 2.2.2 G(R) hàm phụ thu c vào yếu tố hình h c c a h 48 2.1.3 β( ) 2.1.4 ( ) số tán x n tr sol k í 49 ng truy n 50 ết lu n ơn 51 ơn ơn p p t ết kế, chế t o tích h p h Lidar quan trắc sol khí tố u 53 3.1 Bài toán xây d ng h Lidar 53 iii 3.2 Thiết kế h Lidar 54 3.2.1 L a ch n m t số modul ph n c ng 55 3.2.2 Chế t o m t số modul ph n c ng 66 3.3 p n t n số r 70 3.4 u ỉn r 77 3.5 Xây d n ơn trìn p n m m o v l u tr số l u 84 ết lu n ơn 87 ơn y ng ph n m m x l kết o v p ơn p p tín o l p b n k í t k u v 88 t vấn 88 4.1 n m m x l số l u o 89 o t n mv n kết o rt 93 l u o sol k í t Lidar 93 n số k t u t v t số l u r 93 l u o sol k í t rn n u ế t o 94 ố l u o sol k í bằn b n t m k n 96 4.2.2 n tí kết o 97 ết o t t ết b 97 ết o t t ết b r 99 4.3 ín o l p b n k í bằn số l u r 103 ết lu n ơn 114 K T LU N 115 n mụ n trìn k o t ả l n qu n ến lu n n 117 l u t m k ảo 118 Tiếng Vi t 118 Tiếng Anh 118 126 iv D , N N ABL Atmospheric Boundary Layer pr n k í ADC Analog - to -Digital Converter uyển t ơn t - số AERONET Aerosol Robotic Network AFV Apparent Field of View AOD Aerosol Optical Depth APD Avalanche Photodiode AVHRR Advanced Very High Resolution Radiometer n l qu n kế o sol k í tron k í r n n ìn b ểu k ến y qu n ốt t EARLINET European Aerosol Lidar Network Field of View Laser Light Amplification by Stimulated Emission of Radiation LBKQ The planetary boundary layer Lidar Light Detection and Ranging LRS Lidar and Radiometric Sensing MODIS Moderate Resolution Imaging Spectroradiometer MPLNET The NASA Micro-Pulse Lidar Network NASA National Aeronautics and Space Administration v l ết b o b ả o x p n n r qu n trắ k í Châu Âu r FOV sol k í n n ìn k uế n s n bằn p t x kí t í p b n k í ết b p t n n s n v k oản rv ảm b ến p n x ết b o qu n p ìn ản n r qu n n k n v trụ o ỳ ol m ồn p OF Overlap Function PMT Photomultiplier Tube Modules QE Quantum Efficiency u suất l n t S Radiant Sensitivity n yb x SNR Signal-to-Noise Ratio UAV Unmanned Aerial Vehicle WCT Wavelet Covariance Transform B ến Wavelet h p p ơn s α Extinction Coefficient H số suy β Scatter Coefficient H số t n x n n n n qu n ỷ số tín vi n u tr n n tt ể b yk n n l o n s n u D Số hi u hình v Tên hình v Trang ìn Cấu trúc c a t ng khí ìn M ng tr m Lidar khu v ìn Minh h a hình th c nghiên c u khí 11 Ảnh minh h 12 ìn ìn p n m n n u khí ơn p p o sol k í tr c tiếp n b y m n t ết b o sol k í 10 13 ìn 1.6 ấu trú r t o k ểu -static 15 ìn 1.7 ấu trú r t o k ểu ồn trụ 16 ìn ìn Ảnh h Lidar t i Vi n V t lý,Vi n Hàn lâm KH&CNVN 17 ìn 10 Ảnh h Lidar LSA- 18 tr m 16 , r n i h c Bách Khoa Hà N i cấu trúc (a), ảnh h Lidar (b) t tr n ih v t lý Hành tinh t i Manoa, Honolulu, Hawai ìn ìn trí 12 un Cấu trúc h Lidar t , a 19 o n 21 ìn 1.13 Ảnh h Lidar NASA chế t o 23 ìn 14 Ảnh (a) (b) h Lidar LSA-2C 24 ìn Ản 25 ìn 2.1 idar c tài KC.01.21/06-10 nx n tr n s n n s n kí t í ìn í ìn ấu trú v n uy n l k uế ìn ìn M ìn t t kí t n uv b c m t số lo i sol khí ph biến n uy n l Ảnh m un ấp nt ế 33 34 ống n n qu n n o 36 36 n ấu nối bên cao áp C4900 36 s n n n s n k 37 t số lo ìn khối h Lidar ìn 2.8 biểu di n qu trìn ìn 2.9 di n tả trình tín hi u Lidar ìn 2.10 k n u 39 o tín u Lidar khí Minh ho hình h c h Lidar 42 43 44 vii ìn 2.11 ìn 2.12 ìn 3.1 khối mơ tả chất v t lý h Lidar Minh h a ản 45 ởng hàm chồng ch p tín hi u 49 khối h Lidar 54 Hình 3.2 Ảnh khối phát tia laser c a laser LOTIS 55 ìn 3.3 Ảnh khối phát tia laser c a laser Briliant 56 ìn 3.4 Ảnh khối phát tia laser c a laser Litron Nano 57 ìn 3.5 Ản 59 ìn 3.6 up tl s r l s r ấu trú k ối thu 60 ìn 3.7 Ản v sơ cấu trú kín t ìn 3.8 Bản vẽ hình h v -43.9-UV-355-532 ng truy n qua c a thấu kính LUD 50.0 ìn 3.9 Bản vẽ hình h 16 x 15 v ng truy n qua c a thấu kính DCX UV ìn 3.10 Ảnh kính l c F03-532.0-4-1.00 ìn 3.11 Ản v ng truy n t s n k PBSH- 450-1300-100 ìn 3.12 Cấu trúc hình h PMT v th n văn lo i Celestron 61 61 62 63 qu k ối tách chùm nh y b c x t o b c sóng 63 64 ìn 3.13 khối m ch cấp nguồn h Lidar 66 ìn 3.14 m ch nguyên lý cấp nguồn h Lidar 67 Hình 3.15 M ch in m ch cấp nguồn có linh ki n Hình 3.16 M ch nguyên lý, m ch in cấp nguồn v khuế i PMT Hình 3.17 Ảnh b cấp nguồn Hình 3.18 Ảnh chuẩn tr c ánh sáng dùng hai thấu kính 69 Hình 3.19 Bản vẽ k í b gá lắp hai thấu kính 69 Hình 3.20 Ảnh 3D c a b gá hai thấu kính 70 Hình 3.21 Ảnh b ph n quang h c thu 70 Hình 3.22 Mơ hình t ng qt mơ ph ng h Lidar 72 Hình 3.23 Hình 3.24 67 u khiển h số u khiển h số khuế i PMT u thu t toán m p n 68 68 75 (a) tín hi u Rayleigh mơ ph ng hi u t o viii cao 76 T T V uy n u n n , u n( ), " n u n ụn r tron n n u m y ",Tuyển tập cơng trình nghiên c u Vi n Vật ý a c u 2008, tr 1-10 Nguy n n ìn , n ăn run ( ), " t tr ển h thống Lidar s dụng nghiên c u sol khí t i Vi n v t lý – nh ng kết b n u, H i ngh khoa h c kỷ ni 35 ă VAS , pp.9–14 ăn Hả ( ), " ụn k t u t son k í tron t n k í quyển" l uy n ế ếu ( n ụn tron o ), " n u t ết kế, t n số k í quyển", Tr n Xuân Hoài, Âu Duy Tuấn ( l n – u k ển t n s ụn ni 35 ă VAS pp.133–138 m ồv rn o n ), " u n ơ( n u ết o tr n v t t ốn -10 r ết kế t ốn nt o n ", H i ngh khoa h c kỷ ), " sở k í t n " Lê Huy Minh, uy n u n n ( ), " ng dụng công ngh v trụ ể nghiên c u n ản ởng c a t n n ly t ng khí t i xác s dụng tín hi u v tinh khu v c Vi t Nam", Báo cáo ì c cơng ngh ũ ụ 2008-2011 m u n n v n s (2015), " ểm dày quang h c sol khí từ số li u tr m AERONET Vi t Nam so sánh chúng v i số li u MODIS",T p chí Các Khoa h c v t, ố , r –263 10 T uy n u n uấn ( ), " ụn k ảo s t p n bố n t n u v p t tr ển k t u t v m t k í quyển" run t m qu n trắ ( l n m tr ng khơng khí t i th ), " r n n t n v di n biến chất i" A 11 A.Chaikovsky, A Bril, O Dubovik et al (2004), "CIMEL and multiwavelength Lidar measurements for troposphere aerosol altitude distributions investigation, long-range transfer monitoring and regional ecological problems solution: field validation of retrieval techniques", Opt Pura Apl, Vol.37, pp.3241-3246 118 12 A.Jayaraman (2003), pp.1–26 "Lidar and Its Appilactions", Phys Res Lab 13 A P Van Ulden & J.Wieringa (1996), "Atmospheric boundary layer research at Cabauw", Boundary-Layer Meteorol, Vol.78, pp.39–69 14 A Stewart Walker (2018), Lidar Magazine, pp.1–52 15 Alessia Sannino, Antonella Boselli et al (2019), Optimization of the Lidar optical design for measurement of the aerosol extinction vertical profile, EPJ Web of Conferences, pp 1–4 16 Alexander Baklanov & Branko Grisogono (2007), "Atmospheric boundary layers: nature, theory and applications to environmental modelling and security", Springer Sci, pp.1–4 17 Ali H Omar, David M Winker, Mark A Vaughan et al (2009), "The CALIPSO automated aerosol classification and Lidar ratio selection algorithm", J Atmos Ocean Technol,Vol.20, pp.1994–2014, doi: 10.1017/CBO9781107415324.004 18 Anatoly P Chaikovsky, O Dubovik et al (2002), "Methodology to retrieve atmospheric aerosol parameters by combining ground-based measurements of multiwavelength Lidar and sun sky-scanning radiometer", Society of Photo-Optical Instrumentation Engineers (SPIE), Vol.4678, pp.257–268, doi: 10.1117/12.458450 19 Andrea Lammert and Jens Bosenberg (2006), "Determination of the convective boundary-layer height with laser remote sensing", Springer Sci, Vol.119, pp.159–160, doi: 10.1007/s10546-005-9020-x 20 Beyrich F (1997), "Mixing height estimation from sodar data", Atmos Environ, Vol.31, pp.3941–3953 21 C F Bohren and D R Huffman (1983), Absorption and Scattering of Light Small Particles 22 Caicedo V, Rappenglück B, Lefer B, Morris G, Toledo D, Delgado R (2017), "Comparison of aerosol Lidar retrieval methods for boundary layer height detection using ceilometer aerosol backscatter data", pp.1609–1622, doi: 10.5194/amt-10-1609-2017 119 23 Chi O Ao, Duane E Waliser et al (2012), "Planetary boundary layer heights from GPS radio occultation refractivity and humidity profiles", J Geophys Res Atmos, Vol.117, pp.1–18, doi: 10.1029/2012JD017598 24 Collis RTH, Russell PB (2012), "Lidar measurement of particles and gases by elastic backscattering and differential absorption", pp.71–151, doi: 10.1007/3-540-07743-x 25 David S Hall LAH (2011), "High Definition Lidar System", pp.1–26 26 Delbarre H, "Laboratory for Physico-Chemistry of the Atmosphere" 27 Dian J Seidel, Chi O Ao and Kun Li (2010), "Estimating climatological planetary boundary layer heights from radiosonde observations: Comparison of methods and uncertainty analysis", J Geophys Res Atmos, Vol.115, pp.1–15, doi: 10.1029/2009JD013680 28 Dongsong SUN, Sakae KAWATO and Takao KOBAYASHI (2000), "Ultraviolet Rayleigh Lidar for Wind and Temperature Measurements", Opt Rev, Vol.7, pp.555–560 29 E.D.Hinkley (1996), "Laser Monitoring of the Atmosphere", SpringerVerlag 30 Eduardo Quel et al (2007), "Using Lidar to Ozone Layer" 31 Ellsworth J.Welton et al (2014), The NASA Micro - Pulse Lidar Network (MPLNET)", pp.1–5 32 Francesc Rocadenbosch (2003), "Lidar - Aerosol Sensing", In: Encyclopedia of Optical Engineering, pp.1090 - 1102 33 Frank Raes, Timothy Bates, Frank McGovern and Marc Van Liedekerke (2000), "The 2nd Aerosol Characterization Experiment (ACE-2): general overview and main results" 34 Frederick G Fernald (1984), "Analysis of atmospheric Lidar observations: some comments", Appl Opt, Vol.23, pp.652–653 35 George C Holzworth (1964), "Estimates of Mean Maximum Mixing Depths in the Contiguous United States" Vol.92, pp.235–242 36 Gilles Larchevêque (2002), "Development of the Jungfraujoch MultiWavelength Lidar System for Continuous Observations of the Aerosol Optical Properties in the Free Troposphere" 120 37 Hamamatsu (1999), "Photomultiplier tubes Basics and Applications" 38 Hamamatsu (2006), "Photomultiplier Tubes Opening The Future with Photonics" 39 Hamamatsu Photonics (2001), "Metal Package Photomultiplier Tube R7400U Series" 40 Hamamatsu Photonics (2007), "PHOTOMULTIPLIER TUBES Basic and Applications" 41 Hamamatsu Photonics (2008), "Compact High Voltage Power Supply Units C4900", pp.40–41 42 Hamamatsu Photonics (2018), "Photosensor Modules H11526 Series" 43 HUANG Wei-Ping JK et al (2013), "Design of the photomultiplier tube bases for high dynamic range readout in WCDA", Chinese Phys C, Vol.37, doi: 10.1088/1674-1137/37/3/036001 44 Hyo Sang Lee, Silver Spring et al (2003), "Portable Digital Lidar System", No.6.593.582, pp.1–25 45 Hyo Sang Lee, Silver Spring et al (2010), "Enhanced Portable Digital Lidar System" No 7.741.618, pp.1–31 46 Ian M Brooks (2003), "Finding boundary layer top: Application of a Wavelet covariance transform to Lidar backscatter profiles', J Atmos Ocean Technol, Vol.20, pp.1092–1105, doi:10.1175/1520-0426 47 Ian M Brooks et al (2008), "Boundary Layer Structure and Aerosol Properties: Lidar Retrievals and In-Situ Measurements", pp.146–149 48 J N Porter, B R Lienert, S K Sharma, and H W Hubble (2002), " A small portable Mie-Rayleigh Lidar system to measure aerosol optical and spatial properties", J Atmos Ocean Technol, Vol.19, pp.1873–1877, doi: 10.1175/1520-0426 49 James D Klett (1981), "Stable analytical inversion solution for processing Lidar returns", Appl Opt, Vol.20, pp.211–220, doi: 10.1364/ao.20.000211 50 James D Klett (1985), "Lidar inversion with variable backscatter/extinction ratios", Appl Opt, Vol.24, pp.1638–1643, doi: 10.1364/ao.24.001638 121 51 John W Hair, Loren M Caldwell, David A Krueger, and Chiao-Yao She (2001), "High-Spectral-Resolution Lidar with Iodine-Vapor Filters: Measurement of Atmospheric-State and Aerosol Profiles", Appl Opt, Vol.40, pp.5280–5294, doi: 10.1364/ao.40.005280 52 Juan Carlos Fernandez Diaz, William E Carter, Ramesh L Shrestha, and Craig L Glennie (2016), "Lidar Remote Sensing" 53 K J Davis, N Gamage et al (2000), An objective method for deriving atmospheric structure from airborne Lidar observations, J Atmos Ocean Technol, Vol.17, pp.1455–1468, doi: 10.1175/1520-0426 54 Kaufman, Yoram J et al (2002), "A satellite view of aerosols in the climate system", Nature, Vol.419, pp.215–223 55 Kenneth Sassen and Gregory C Dodd (1982), "Lidar crossover function and misalignment effects", Appl Opt, Vol.21, pp.3162–3165, doi: 10.1364/ao.21.003162 56 Lammert A, Bösenberg J (2006), "Determination of the convective boundary-layer height with laser remote sensing" Boundary-Layer Meteorol, Vol.119, pp.159–170, doi: 10.1007/s10546-005-9020-x 57 Laser Lotis, "YAG LASER SYSTEM, LS-2137U" 58 Laurent Menut, Cyrille Flamant, Jacques Pelon and Pierre H Flamant (1999),"Urban boundary-layer height determination from Lidar measurements over the Paris area", Appl Opt, Vol.38, pp.945–954 doi: 10.1364/ao.38.000945 59 M J Granados-Muñoz, F Navas-Guzmán, J A Bravo-Aranda et al (2012), "Automatic determination of the planetary boundary layer height using Lidar: One-year analysis over southeastern Spain", J Geophys Res Atmos, Vol.117, pp.1–10, doi: 10.1029/2012 60 M Kulmala, A Asmi et al (2011), "European Integrated project on Aerosol Cloud Climate and Air Quality interactions (EUCAARI)", Atmos Chem Phys, Vol.11, pp.13061–130143, doi: 10.5194/acp-1113061-2011 61 Jiandong Mao (2012) "Noise reduction for lidar returns using local threshold wavelet analysis", Opt Quantum Electron, Vol.43, pp.59–68 doi: 10.1007/s11082-011-9503-6 122 62 Mego Pinandito et al (2000), "Mie Scattering Lidar Observation of Aerosol Vertical Profiles in Jakarta", Environ Sci, Vol.13, pp.205-216 63 Petra Seibert et al (2000), "Review and intercomparison of operational methods for the determination of the mixing height", Atmos Environ, Vol.34, pp.1001–1027, doi: 10.1016/0360-3016(85)90285-8 64 Photek (2016), "Photomultiplier & Photodiode User Guide" 65 QI YuLei, GE JinMing and HUANG JianPing (2013), "Spatial and temporal distribution of MODIS and MISR aerosol optical depth over northern China and comparison with AERONET", Chinese Sci Bull, Vol.58, pp.2497–2506, doi: 10.1007/s11434-013-5678 66 R J Charlson and J Heintzenberg (1995), "Aerosol forcing of climate", New-York, John Wiley Sons, Ltd 67 R M Measures (1984), "Laser remote sensing fundamentals and application" 68 R.G Barry and R J Chorley (1998), "Atmosphere, Weather & Climate" 69 R.M Banta, Y L Pichugina et al (2013), "Insight into Wind Properties in the Turbine - Rotor Layer of the Atmosphere from High -Resolution Doppler Lidar" 70 Rajitha Paleti , Y Bhavani Kumar and T Krishna Chaitanya (2013), "Wavelet transform method for deriving atmospheric boundary layer height from Lidar signals", Int J Eng Technol, Vol.5, pp.1465–1473 71 Shengzhen Zhou, Perry K Davy, Xuemei Wang et al (2016), "High time-resolved elemental components in fine and coarse particles in the Pearl River Delta region of Southern China: Dynamic variations and effects of meteorology", Sci Total Environ, Vol.572, pp.634–648 doi: 10.1016/j.scitotenv.2016.05.194 72 Shuyan Liu and Xin-Zhong Liang (2010), "Observed diurnal cycle climatology of planetary boundary layer height", J Clim, Vol.23, pp.5790–5809, doi: 10.1175/2010JCLI3552.1 73 T C SPANGLER and R A DIRKS (1974), "Meso-scale variations of the urban mixing height" 123 74 T.S Bates, B J Huebert, J.L Gras et al (1998), "International Global tmosp r m stry ( ) roj t’s rst rosol r t r z t on Experiment (ACE 1)", J Geophys Res Atmos, Vol 103, pp.16,297– 16,318 75 Timothy S Bates, Barry J Huebert et al (1998), "International Global tmosp r m stry ( ) roj t’s rst rosol r t r z t on Experiment (ACE 1) - Overview", J Geophys Res, Vol.103, pp.16297– 16318 76 Vishnu R, Bhavani Kumar Y, Rao T, et al (2016), "Development of Lidar sensor for cloud-based measurements during convective conditions", Remote Sens Atmos, Clouds, Precip VI, doi: 10.1117/12.2223674 77 Vladimir A Kovalev and William E Eichinger (2004), "ELASTIC LIDAR: Theory, Practice, and Analysis Methods" 78 Volker Freudenthaler (2008), "The telecover test: A quality assurance tool for the optical part of a Lidar system", 24th Int Laser Radar Conf Boulder, pp.1–2 79 Wandinger U (2005), "Introduction to Lidar" 80 Weitkamp Claus (2005), "Lidar Range-Resolved Optical Remote Sensing of the Atmosphere" 81 Y Bhavani Kumar and S Purusotham (2010), "Mathematical algorithms for determination of mixed layer height from laser radar signals", Int J Comput Sci Eng, Vol.2, pp.2059–2063 82 Z Wang, X Cao, L Zhang et al (2012), "Lidar measurement of planetary boundary layer height and comparison with microwave profiling radiometer observation", Atmos Meas Tech, Vol.5, pp.1965– 1972, doi: 10.5194/amt-5-1965 83 http://mplnet.gsfc.nasa.gov(2011) Micro - Pulse Lidar Network 84 http://weather.uwyo.edu/upperair/sounding.html weather.uwyo.edu 85 http://www litronlasers.com http://www litronlasers.com 86 http://www.lotis-tii.com (2006) http://www.lotis-tii.com 87 http://www.quantel-laser.com http://www.quantel-laser.com 124 88 http//www.celestron.com (2005) http//www.celestron.com 89 https://aeronet.gsfc.nasa.gov/ https://aeronet.gsfc.nasa.gov/ 90 https://sg.tek.com/oscilloscope/mdo4000c-mixed-domain-oscilloscope (2012) Mixed Domain Oscilloscopes 91 https://www-calipso.larc.nasa.gov/ The Cloud-Aerosol Infrared Pathfinder Satellite Observation (CALIPSO) Lidar 92 https://www.cvilaseroptics.com/ (2014) Optics & Lasers Catalog 93 https://www.earlinet.org zzhttps://www.earlinet.org 94 SigmaSpace (2016) MiniMPL-532 BUDGETARY QUOTATION 125 and H 1.1 ơn trìn m p p ơn trìn r n C t n số k t u t p n 126 n r theo 1.2 ố n n t u t p số l u 1.3 ố t ết l p ế ồn b n o 127 1.4 ố t ết l p ấu ìn o k n v o 128 o n k ố t t n số o Lidar C WCT function pblh = wct_pbl(f,a,dz,z,r,lb,lt) step = round(a/dz); for i=lb:lt s = a/dz/2; c = i+s; if c