We report on a new method to extract thermodynamic properties from the density of partition function zeroes on finite lattices. This allows direct determination of the order and strength of phase transitions numerically. Furthermore, it enables efficient distinguishing between first- and second-order transitions, elucidates crossover between them and illuminates the origins of finite-size scaling. The power of the method is illustrated in typical applications for both Fisher and Lee-Yang zeroes.