In this study, we modelled the incidence of COVID-19 cases and hospitalisations by basic health areas (ABS) in Catalonia. Spatial, temporal and spatio-temporal incidence trends were described using estimation methods that allow to borrow strength from neighbouring areas and time points. Specifically, we used Bayesian hierarchical spatio-temporal models estimated with Integrated Nested Laplace Approximation (INLA). An exploratory analysis was conducted to identify potential ABS factors associated with the incidence of cases and hospitalisations. High heterogeneity in cases and hospitalisation incidence was found between ABS and along the waves of the pandemic. Urban areas were found to have a higher incidence of COVID-19 cases and hospitalisations than rural areas, while socio-economic deprivation of the area was associated with a higher incidence of hospitalisations. In addition, full vaccination coverage in each ABS showed a protective effect on the risk of COVID-19 cases and hospitalisations.