Unify metrics computation for original and surrogate timeseries

So far, `estimate_indicator_changes` and `significant_transitions` have very similar code but don't call the same, easier to maintain, functions. We should unify this such that both of them call something like:

```
initialize_indicator_changes()
estimate_indicator_changes!()
```