Allow overscan region to be along either axis

As reported in #70, ``subtract_overscan`` assumed that the axis along which the overscan should be aggregated (with either median or mean) was the second of the two axes. This allows the overscan to be in either axis (though it will fail in the unlikely case that the overscan is square).