1. Feature detection
Localize extrema by fitting a quadratic
1) Sub-pixel/sub-scale interpolation using Taylor
expansion
2) Take derivative and set to zero