A cluster-based method for hydrometeor classification using polarimetric variables. Part I: Interpretation and analysis

Hydrometeor classification methods using polarimetric radar variables rely on probability density functions (PDFs) or membership functions derived empirically or by using electromagnetic scattering calculations. This paper describes an objective approach based on cluster analysis to deriving the PDFs. An iterative procedure with K-means clustering and expectation-maximization clustering based on Gaussian mixture models is developed to generate a series of prototypes for each hydrometeor type from several radar scans. The prototypes are then grouped together to produce a PDF for each hydrometeor type, which is modeled as a Gaussian mixture. The cluster-based method is applied to polarimetric radar data collected with the CP-2 S-band radar near Brisbane, Queensland, Australia. The results are illustrated and compared with theoretical classification boundaries in the literature. Some notable differences are found. Automated hydrometeor classification algorithms can be built using the PDFs of polarimetric variables associated with each hydrometeor type presented in this paper.