Iterative kernel background estimator

Made a new PR as couldn't push to the old one: #173

Still in progress. To do:
 - [x] Add unit tests
 - [x] Finish tutorial
 - [x] Include initial image and mask image additionally in the tutorial
 - [x] Clean up code formatting
 - [x] Add docstrings
(Anything else?)

Please comment if you see anything to be done or changed.
You can get the `EnergyDependentTablePSF` object directly ... I don't know if it's always been like this or if I changed it, but see [here](https://gammapy.readthedocs.org/en/latest/api/gammapy.datasets.FermiGalacticCenter.html#gammapy.datasets.FermiGalacticCenter.psf).