Generalize Sersic bounding box calc to n=1,4 and move into a function
Clean up calculation of threshold surface brightness
Add new --band option to control which AB magnitude to use from catalog