Fits画像から分散を求めるコード

TL; DR glasnsciです。astropyを使って天体写真等を撮影したfits画像から輝度分散を求めるコードのようなものです。Python3で動きます。