We explore the QCD equation of state at finite baryon density through an expansion in baryon number fugacity, making use of the recent lattice data on the four leading Fourier coefficients of the expansion. A state-of-the-art description of the lattice data at both imaginary $\mu_B$ and $\mu_B = 0$ is provided by the cluster expansion model (CEM). The CEM pressure function has three temperature dependent input parameters, which we fix by parameterizing the available lattice QCD data. This results in a crossover model of the QCD equation of state at finite baryon density, which can be used in fluid dynamical simulations of heavy-ion collisions.