silv_dominant_height function

Calculates the dominant height