looking into Hogg formula for deVauc
switched from one-d to two-d Gaussian, Hogg profile works now
working on IDL code to grab pixels