Background returns incorrect units

[cshanahan1]

The following should return flux in DN, not DN**2:

img = np.tile((np.arange(1., 13)), (13, 1)) * u.DN trace = FlatTrace(img, 2) background = Background(img, trace) background.bkg_image()