The invention discloses rendering sound field signals, such as Higher-Order Ambisonics (HOA), for arbitrary loudspeaker setups, where the rendering results in highly improved localization properties and is energy preserving. This is obtained by a new type of decode matrix for sound field data, and a new way to obtain the decode matrix. In a method for rendering an audio sound field representation for arbitrary spatial loudspeaker setups, the decode matrix (D) for the rendering to a given arrangement of target loudspeakers is obtained by steps of obtaining a number (L) of target speakers, their positions (L), positions (S) of a spherical modeling grid and a HOA order (N), generating (141) a mix matrix (G) from the positions (S) of the modeling grid and the positions (L) of the speakers, generating (142) a mode matrix ({tilde over (Ψ)})from the positions (S) of the spherical modeling grid and the HOA order, calculating (143) a first decode matrix ({circumflex over (D)}) from the mix matrix (G) and the mode matrix ({tilde over (Ψ)}), and smoothing and scaling (144, 145) the first decode matrix ({circumflex over (D)}) with smoothing and scaling coefficients.