Allow `calculate_feature_density` to run arbitrary aggregate functions over grid cells

[atmorling]

Following on from #555, it would be nice to extend this functionality to arbitrary summary statistics