Regional frequency analysis (RFA) can reduce uncertainties in the estimations of return levels, provided that homogeneous regions can be delineated. In the framework of extreme marine events, a physically based method to form homogeneous regions by identifying typical storms footprints is proposed. First, a spatiotemporal declustering procedure is employed to detect storms generating marine extremes. Second, the identification of the most typical storms footprints relies on a clustering algorithm based on a criterion of storm propagation. The resulting regions are readily explicable: sites from a given region are likely to be impacted by the same storms, and any storm impacting a region is likely to remain enclosed in this region. This procedure is fairly simple to implement, as the only information required is the time of occurrence of the observed extremes. An application to the estimation of extreme significant wave heights from the numerical sea-state database ANEMOC-2 is given. Six regions, both physically and statistically homogeneous, are delineated in the North-East part of the Atlantic Ocean. It is shown that the identification of storms footprints allows the increase of the overall statistical homogeneity. Combined with RFA, the proposed method highlights regional differences in the spatial extent and intensity of storms.