Fix call to Map.ones in phasemask_fourquadrant.
Make create_fitsheader_for still handle scalar as 1-element 1D array. Update tests.