To overcome the problem that traditional ice sensors cannot detect multiple icing areas in a wide range of airfoils, a multiple icing area localization method based on the grayscale centroid method is proposed. The vibration spectrum of piezoelectric array before and after icing is collected, and the attenuation rate of spectrum amplitude is extracted as the positioning eigenvalue. Combined with the reconstruction algorithm based on multi-point icing probability detection (RAMIPD) and the grayscale centroid method, the center coordinates of the multiple icing areas are calculated to get the positions of the icing regions on a plane or a curved surface. When the sweep frequency range is the full band of 0. 5~ 5 kHz and the diameter of icing area is 70 mm, the average positioning error of single icing area on the aluminum plate is 27. 4 mm. The average positioning error of two icing areas is 29. 0 mm. The average lateral positioning error of icing areas on the airfoil is 22. 6 mm. To improve the positioning accuracy, a feature extraction method for sensitive frequency band selection based on wavelet packet decomposition (WPD) is proposed. Experimental results show that the single icing area positioning accuracy on the aluminum plate is increased by 34. 59% , the fluctuation of the positioning error is reduced by 45. 67% , the lateral positioning accuracy on the airfoil is increased by 54. 87% , and the fluctuation of the lateral positioning error is reduced by 46. 63% . The proposed icing area localization method can detect the multiple icing areas with high precision and provide a new idea for icing area detection.