It has just occurred to me that this area averaged bin weighting is not correct for the higher frequency bins in which the lower `ell` edge is greater than `pi / self.dx`, due to the annular `ell` region being cut off by the square grid edges.  This seems like something that is a pain to calculate and correct for, but maybe we should mention this in the dox?  I'm happy to add a few words if people agree.