In this paper, we first prove that the local time associated with symmetric $α$-stable processes is of bounded $p$-variation for any $p>\frac{2}{α-1}$ partly based on Barlow's estimation of the modulus of the local time of such processes.\,\,The fact that the local time is of bounded $p$-variation for any $p>\frac{2}{α-1}$ enables us to define the integral of the local time $\int_{-\infty}^{\infty}\triangledown_-^{α-1}f(x)d_x L_t^x$ as a Young integral for less smooth functions being of bounded $q$-varition with $1\leq q<\frac{2}{3-α}$. When $q\geq \frac{2}{3-α}$, Young's integration theory is no longer applicable. However, rough path theory is useful in this case. The main purpose of this paper is to establish a rough path theory for the integration with respect to the local times of symmetric $α$-stable processes for $\frac{2}{3-α}\leq q< 4$.