Regularity for subelliptic PDE through uniform estimates in multi-scale geometries