In this paper, we propose a spatially constrained clustering problem belonging to the family of p-regions problems. Our formulation is motivated by the recent developments of economic complexity on the evolution of the economic output through key interactions among industries within economic regions. The objective of this model consists in aggregating a set of geographic areas into a prescribed number of regions (so-called innovation ecosystems) such that the resulting regions preserve the most relevant interactions among industries. We formulate the p-Innovation Ecosystems model as a mixed-integer programming (MIP) problem and propose a heuristic solution approach. We explore a case involving the municipalities of Colombia to illustrate how such a model can be applied and used for policy and regional development.