On metric measure spaces with sub-Gaussian heat kernel behavior in small time, we obtain a sufficient condition to solve Wick renormalized stochastic quantization equations with polynomial interaction. Given the power of the nonlinearity, the local solution condition depends on the Hausdorff dimension $d_h$, the walk dimension $d_w$, and the maximal spatial Hölder regularity of the heat kernel $Θ$. A slightly more restrictive condition based on the same parameters is required for a global solution. For all global solutions, we construct an invariant measure for the Markov process defined by the solution. Our results apply to many rough spaces such as Barlow--Kigami type fractals as well as their Cartesian products and open up the possibility of making rigorous various structures in quantum field theory and statistical mechanics in non-integer dimensions. In the process, we build entirely from the short-time heat semigroup the necessary analytic framework that accommodates the issues which come with allowing rough local geometry.