def METHOD_NAME(self,psf):
    sz = psf.shape
    psfc = psf[:,:,sz[2]/2].squeeze()
    X, Y = np.meshgrid(np.arange(sz[0]),np.arange(sz[1]))
    xc = (X*psfc).sum()/psfc.sum()
    yc = (Y*psfc).sum()/psfc.sum()
    psfs = np.roll(np.roll(psf,int(np.round(sz[0]/2-xc)),0),int(np.round(sz[1]/2-yc)),1)
    return psfs