Data generated by airborne laser scans enable the characterization of surface vegetation for any application that might need it, such as forest management, modeling for numerical weather prediction, or wind energy estimation. In this work we present a new algorithm for calculating the vegetation density using data from airborne laser scans. The new routine is more robust than earlier methods, and an implementation in popular programming languages accompanies the article to support new users.