Add function calculate_weights()

[KDolos]

... depending on sampling density as a simply way to consider differences in sampling design.
Consider additionally Germany nested design