Geospatial fire behaviour and fire hazard simulators, fire effects models and smoke emission software commonly use standard fuel models in order to simplify data collection and the inclusion of complex fuel scenarios. These fuel models are often mapped using remotely sensed data. However, given the great complexity of fuelbeds, with properties that vary widely in both time and space, the use of these standard fuel models can greatly limit accurate fuel mapping. This affects fuel hazard assessment, fuel reduction treatment plans, fire management decision-making and evaluation of the environmental impact of wildfire. In this study, we developed unique customized fire behaviour fuel models for shrub and bracken communities, by using k-medoids clustering analysis based on both fuel structural characteristics and potential fire behaviour. We used an original database of 722 destructive sample plots in nine different shrub and bracken communities covering the entire distribution area in Galicia (NW Spain), one of the regions in Europe most affected by forest fires. Measurements of cover, height and fuel fractions loads differentiated by size and vegetative state (live or dead) were used to estimate the potential rate of fire spread with five different models including fireline intensity, heat per unit area and the flame length for each sampling site and considering extreme environmental conditions. The optimal number of clusters was established by combining practical knowledge about the shrubland communities under study and their associated fire behaviour, with maximization of the mean value of the silhouette variable and minimization of the within-cluster sum of squares. The structural characteristics of the medoids derived from the analysis were associated with each of the proposed customized fuel models. Finally, a simple dichotomous classification based only on shrub height was developed to enable construction of spatially explicit fuel model maps based on remotely sensed data. Thus, the methodology applied allows generation of a more realistic representation of fuel distribution in the landscape, based on fuel structure measurements of natural regional ecosystems rather than on the use of standard models. We believe that the proposed methodology is generally applicable to communities composed of other shrub and fern species in different biogeographical regions.