This work is devoted to deriving the Onsager-Machlup action functional for stochastic partial differential equations with (non-Gaussian) Levy process as well as Gaussian Brownian motion. This is achieved by applying the Girsanov transformation for probability measures and then by a path representation. This enables the investigation of the most probable transition path for infinite dimensional stochastic dynamical systems modeled by stochastic partial differential equations, by minimizing the Onsager-Machlup action functional.