Flux point method

Implements #50 .
Use `scipy.stats.gmean` to compute log center positions