code: def kcenter(self):
          
          return np.array([
              np.abs(np.fft.fftshift(np.fft.fftfreq(q))).argmin()
              for q in self.shape
          ]).astype('float')
docs: Return the frequency center of the tile (says fftshift)