本研究使用影像處理技術發展雨滴譜儀。此系統使用攝影機擷取降落的雨滴，並且使用影像處理技術偵測分析雨滴特性諸如雨滴形狀與降落速度。藉由觀測出的資料可以用來偵測降雨型態以及增進小範圍降雨分析的準確度。在此系統中，使用一遠心鏡頭與攝影機搭配可以減少因失焦造成的模糊與變形現象進而增進量測雨滴的準確度。然而，使用遠心鏡頭必須提供相當足量的光照；因此，我們增設一個光源以補足亮度。由於攝影機的幀幅限制，無法於單次量測完整偵測並分析所有參數特性。正因如此，此系統使用兩種長短曝光時間來分別偵測雨滴形狀與雨滴速度。再者，此系統的影像處理步驟也會於本研究詳細介紹。
最後，我們採取三種不同實驗以驗證此系統︰玻璃珠實驗驗證了此系統架構與處理步驟之可行性；灑水實驗模擬了玻璃珠實驗所不及的較複雜與充滿不確定因素的降雨情形。最後的降雨實驗則是拍攝並分析實際降雨情形，而後將計算出的資料與中央氣象局新竹氣象站之統計資料進行比較。在降雨量的估測我們獲得了14.42%的誤差。

In this sturdy, an image-based disdrometer has been developed. The system uses a CCD camera to capture falling raindrops and applies proper image processing techniques to determine the characteristics of raindrops. The observed data can be used to analyze raining condition and increase the accuracy of analyzing precipitation in small area. In the system, a telecentric lens is applied with the camera to reduce the blur effect and distortion from out of focus, which is easier to measure the shape more accurately. However, the telecentric lens needs enough illuminance to capture raindrops, so a light source is also applied. Due to the limitation of camera frame rate, it is difficult to measure the raindrop shape and velocity in one measurement. That is, two different exposure time measurements, 100us and 2000us, are applied to measure shape and velocity respectively. Furthermore, the details of image processing procedures are also introduced.
Finally, three kinds of experiments were done to verify the system. Marble experiments verified the feasibility and performance of the system structure and image processing procedures. Water sprinkling experiments simulated the raining condition with more uncertainty than the marble experiment. Last, real raining condition pictures were taken and analyzed by the system. The analyzed data was compared with the data from Hsinchu weather station, Central Weather Bureau, Taiwan. An error rate of 14.42% in rain rate was obtained.

[1] K. Droegemeier, J. Smith, S. Businger, C. Doswell III, J. Doyle, C. Duffy, et al., "Hydrological aspects of weather prediction and flood warnings: Report of the Ninth Prospectus Development Team of the US Weather Research Program," Bulletin of the American Meteorological Society, vol. 81, pp. 2665-2680, 2000.
[2] H. Leijnse, R. Uijlenhoet, C. Van De Beek, A. Overeem, T. Otto, C. Unal, et al., "Precipitation measurement at CESAR, the Netherlands," Journal of Hydrometeorology, vol. 11, pp. 1322-1329, 2010.
[3] T. Kozu and K. Nakamura, "Rainfall parameter estimation from dual-radar measurements combining reflectivity profile and path-integrated attenuation," Journal of Atmospheric and Oceanic Technology, vol. 8, pp. 259-270, 1991.
[4] G. B. Foote and P. Du Toit, "Terminal velocity of raindrops aloft," Journal of Applied Meteorology, vol. 8, pp. 249-253, 1969.
[5] J. O. Laws and D. A. Parsons, "The relation of raindrop‐size to intensity," Eos, Transactions American Geophysical Union, vol. 24, pp. 452-460, 1943.
[6] A. Jameson and A. Kostinski, "What is a raindrop size distribution?," Bulletin of the American Meteorological Society, vol. 82, pp. 1169-1177, 2001.
[7] M. Fernandez-Raga, C. Palencia, C. Tomas, A. Calvo, A. Castro, and R. Fraile, "Rain research with disdrometers: a bibliometric review," Atmospheric Measurement Techniques Discussions, vol. 4, pp. 6041-6068, 2011.
[8] D. Atlas, R. Srivastava, and R. S. Sekhon, "Doppler radar characteristics of precipitation at vertical incidence," Reviews of Geophysics, vol. 11, pp. 1-35, 1973.
[9] H. Messer-Yaron, P. Alpert, A. Zinevich, and O. Goldshtein, "Monitoring and Mapping of Atmospheric Phenomena," ed: Google Patents, 2006.
[10] R. Doviak and D. Zrnic, "Doppler radar and weather observations," 1984.
[11] J. Joss and A. Waldvogel, "Raindrop size distribution and sampling size errors," Journal of the Atmospheric Sciences, vol. 26, pp. 566-569, 1969.
[12] X. Liu, T. Gao, and L. Liu, "A video precipitation sensor for imaging and velocimetry of hydrometeors," Atmospheric Measurement Techniques Discussions, vol. 6, pp. 10165-10189, 2013.
[13] N. A. Sivasubramanian and J. R. Saylor, "Application of a histogram modification algorithm to the processing of raindrop images," Optical Engineering, vol. 47, pp. 037011-037011-10, 2008.
[14] J. Saylor and N. Sivasubramanian, "Edge detection methods applied to the analysis of spherical raindrop images," Applied optics, vol. 46, pp. 5352-5367, 2007.
[15] D. Saxena and J. Saylor, "Use of thresholding algorithms in the processing of raindrop imagery," Applied optics, vol. 45, pp. 2672-2688, 2006.
[16] J. Zhidong, G. Taichang, Z. Shijun, L. Zhitian, and Z. Dongli, "High speed data acquisition and processing system for precipitation detection based on linear CCD," in 2013 the International Conference on Remote Sensing, Environment and Transportation
[17] J. S. Marshall and W. M. K. Palmer, "The distribution of raindrops with size," Journal of meteorology, vol. 5, pp. 165-166, 1948.
[18] C. W. Ulbrich, "Natural variations in the analytical form of the raindrop size distribution," Journal of Climate and Applied Meteorology, vol. 22, pp. 1764-1775, 1983.
[19] G. Feingold and Z. Levin, "The lognormal fit to raindrop spectra from frontal convective clouds in Israel," Journal of Climate and Applied Meteorology, vol. 25, pp. 1346-1363, 1986.
[20] R. Gunn and G. D. Kinzer, "The terminal velocity of fall for water droplets in stagnant air," Journal of Meteorology, vol. 6, pp. 243-248, 1949.
[21] L. J. Battan, "Some observations of vertical velocities and precipitation sizes in a thunderstorm," Journal of Applied Meteorology, vol. 3, pp. 415-420, 1964.
[22] K. Beard, "Terminal velocity and shape of cloud and precipitation drops aloft," Journal of the Atmospheric Sciences, vol. 33, pp. 851-864, 1976.
[23] D. L. Mitchell, "Use of mass-and area-dimensional power laws for determining precipitation particle terminal velocities," Journal of the atmospheric sciences, vol. 53, pp. 1710-1723, 1996.
[24] S. Niu, X. Jia, J. Sang, X. Liu, C. Lu, and Y. Liu, "Distributions of raindrop sizes and fall velocities in a semiarid plateau climate: Convective versus stratiform rains," Journal of Applied Meteorology and Climatology, vol. 49, pp. 632-645, 2010.
[25] M. Thurai, W. Petersen, A. Tokay, C. Schultz, and P. Gatlin, "Drop size distribution comparisons between Parsivel and 2-D video disdrometers," Advances in Geosciences, vol. 30, pp. 3-9, 2011.
[26] A. Tokay, A. Kruger, and W. F. Krajewski, "Comparison of drop size distribution measurements by impact and optical disdrometers," Journal of Applied Meteorology, vol. 40, pp. 2083-2097, 2001.
[27] M. Schönhuber, H. Urban, J. P. Baptista, W. Randeu, and W. Riedler, "Weather radar versus 2D-video disdrometer data," ed: Unesco Press, 1997, pp. 159-171.
[28] V. Bringi, T. Seliga, and E. Muelle, "First comparisons of rainfall rates derived from radar differential reflectivity and disdrometer measurements," Geoscience and Remote Sensing, IEEE Transactions on, pp. 201-204, 1982.
[29] E. A. Brandes, G. Zhang, and J. Vivekanandan, "An evaluation of a drop distribution-based polarimetric radar rainfall estimator," Journal of Applied Meteorology, vol. 42, pp. 652-660, 2003.
[30] J. Joss and A. Waldvogel, "Ein spektrograph für niederschlagstropfen mit automatischer auswertung," Pure and Applied Geophysics, vol. 68, pp. 240-246, 1967.
[31] P. Winder and K. Paulson, "The measurement of rain kinetic energy and rain intensity using an acoustic disdrometer," Measurement Science and Technology, vol. 23, p. 015801, 2012.
[32] M. Löffler-Mang and J. Joss, "An optical disdrometer for measuring size and velocity of hydrometeors," Journal of Atmospheric and Oceanic Technology, vol. 17, pp. 130-139, 2000.
[33] M. Schönhuber, G. Lammer, and W. Randeu, "One decade of imaging precipitation measurement by 2D-video-distrometer," Advances in Geosciences, vol. 10, pp. 85-90, 2007.
[34] J. N. Kapur, P. K. Sahoo, and A. K. Wong, "A new method for gray-level picture thresholding using the entropy of the histogram," Computer vision, graphics, and image processing, vol. 29, pp. 273-285, 1985.
[35] P.-Y. Yin, "Maximum entropy-based optimal threshold selection using deterministic reinforcement learning with controlled randomization," Signal Processing, vol. 82, pp. 993-1006, 2002. P.-Y. Yin, "Maximum entropy-based optimal threshold selection using deterministic reinforcement learning with controlled randomization," Signal Processing, vol. 82, pp. 993-1006, 2002.
[36] C.-K. Leung and F. Lam, "Performance analysis for a class of iterative image thresholding algorithms," Pattern Recognition, vol. 29, pp. 1523-1530, 1996.
[37] K. V. Beard and R. J. Kubesh, "Laboratory measurements of small raindrop distortion. Part 2: Oscillation frequencies and modes," Journal of the atmospheric sciences, vol. 48, pp. 2245-2264, 1991.
[38] K. Andsager, K. V. Beard, and N. F. Laird, "Laboratory measurements of axis ratios for large raindrops," Journal of the Atmospheric Sciences, vol. 56, pp. 2673-2683, 1999.
[39] H. Pruppacher and K. Beard, "A wind tunnel investigation of the internal circulation and shape of water drops falling at terminal velocity in air," Quarterly Journal of the Royal