In this paper, we present a scheme for constructing density functions for systems that are almost globally asymptotically stable (i.e., systems for which all trajectories converge to an equilibrium except for a set of measure zero) using navigation functions (NFs). Although recently-proven converse theorems guarantee the existence of density functions for such systems, such results are only existential and the construction of a density function for almost globally asymptotically stable systems remains a challenging task. We show that for a specific class of dynamical systems that are defined based on an NF, a density function can be easily derived from the system's underlying NF.