We present results from the 2D anisotropic baryon acoustic oscillation (BAO) signal present in the final data set from the WiggleZ Dark Energy Survey. We analyse the WiggleZ data in two ways: first using the full shape of the 2D correlation function and secondly focusing only on the position of the BAO peak in the reconstructed data set. When fitting for the full shape of the 2D correlation function we use a multipole expansion to compare with theory. When we use the reconstructed data we marginalize over the shape and just measure the position of the BAO peak, analysing the data in wedges separating the signal along the line of sight from that parallel to the line of sight. We verify our method with mock data and find the results to be free of bias or systematic offsets. We also redo the pre-reconstruction angle-averaged (1D) WiggleZ BAO analysis with an improved covariance and present an updated result. The final results are presented in the form of Ω   , ( ), and ( ) for three redshift bins with effective redshifts = 0.44, 0.60, and 0.73. Within these bins and methodologies, we recover constraints between 5 and 22 per cent error. Our cosmological constraints are consistent with flat ΛCDM cosmology and agree with results from the Baryon Oscillation Spectroscopic Survey.