Close issue 3

Rewrite objective function to update using incremental changes, only
recomputing the pixels that have been modified in the current increment.
This was done to speed up simulated annealing as the objective function,
especially calculating the regional energy, was extremely slow.