dd = np.diag(np.rot90(f_xy))