Documenting the heterogeneity of rainfall regimes is a prerequisite for water resources management, mitigation of risks associated to extremes weather events and for impact studies. In this paper, we present a method for regionalization of rainfall over the Peruvian Pacific slope and coast, which is the main economic zone of the country and concentrates almost 50% of the population. Our approach is based on a two-step process based on k-means clustering followed by the regional vector method (RVM) applied to a network of 145 rainfall stations covering the period 1964–2011. The advantage of combining cluster analysis and RVM is demonstrated compared with just applying each of these methods. Nine homogeneous regions are identified that depict the salient features of the rainfall variability over the study area. A detailed characterization of the rainfall regime in each of the identified regions is presented in response to climate variability at seasonal and interannual timescale. They are shown to grasp the main modes of influence of the El Niño Southern Oscillation (ENSO), that is, increased rainfall over downstream regions in northern Peru during extreme El Niño events and decreased rainfall over upstream regions along the Pacific slope during central Pacific El Niño events. Overall our study points to the value of our two-step regionalization procedure for climate impact studies.