In this paper, we present a new method for estimating the L 2 -gain of systems governed by 2nd order linear Partial Differential Equations (PDEs) in two spatial variables, using semidefinite programming. It has previously been shown that, for any such PDE, an equivalent Partial Integral Equation (PIE) can be derived. These PIEs are expressed in terms of Partial Integral (PI) operators mapping states in L 2 [Ω], and are free of the boundary and continuity constraints appearing in PDEs. In this paper, we extend the 2D PIE representation to include input and output signals in R n , deriving a bijective map between solutions of the PDE and the PIE, along with the necessary formulae to convert between the two representations. Next, using the algebraic properties of PI operators, we prove that an upper bound on the L 2 -gain of PIEs can be verified by testing feasibility of a Linear PI Inequality (LPI), defined by a positivity constraint on a PI operator mapping R n × L 2 [Ω]. Finally, we use positive matrices to parameterize the cone of positive PI operators on R n × L 2 [Ω], allowing feasibility of the L 2 -gain LPI to be tested using semidefinite programming. We implement this test in the MATLAB toolbox PIETOOLS, and demonstrate that this approach allows an upper bound on the L 2 -gain of PDEs to be estimated with little conservatism.