This paper develops and implements a new algorithm for calculating wave trace invariants of a bounded plane domain around a periodic billiard orbit. The algorithm is based on a new expression for the localized wave trace as a special multiple oscillatory integral over the boundary, and on a Feynman diagrammatic analysis of the stationary phase expansion of the oscillatory integral. The algorithm is particularly effective for Euclidean plane domains possessing a ℤ2 symmetry which reverses the orientation of a bouncing ball orbit. It is also very effective for domains with dihedral symmetries. For simply connected analytic Euclidean plane domains in either symmetry class, we prove that the domain is determined within the class by either its Dirichlet or Neumann spectrum. This improves and generalizes the best prior inverse result that simply connected analytic plane domains with two symmetries are spectrally determined within that class.