Characterizing gas bubbles in liquids is crucial to many biomedical, environmental and industrial applications. In this paper a novel method is proposed for the classification of bubble sizes using ultrasound analysis, which is widely acknowledged for being non-invasive, non-contact and inexpensive. This classification is based on 2D templates, i.e. the average spectrum of events representing the trace of bubbles when they cross an ultrasound field. The 2D patterns are obtained by capturing ultrasound signals reflected by bubbles. Frequency-domain based features are analyzed that provide discrimination between bubble sizes. These features are then fed to an artificial neural network, which is designed and trained to classify bubble sizes. The benefits of the proposed method are that it facilitates the processing of multiple bubbles simultaneously, the issues concerning masking interference among bubbles are potentially reduced and using a single sinusoidal component makes the transmitter-receiver electronics relatively simpler. Results from three bubble sizes indicate that the proposed scheme can achieve an accuracy in their classification that is as high as 99%.