Change default spectral averaging method to 'median'
Created by: alurban
This PR fixes #1228 (closed) by changing the default spectral averaging method to 'median'
across the codebase. This includes changes to reference values in the unit tests for gwpy.astro
to reflect median averaging.