The discovery of the Fermi Bubbles has led to a lot of interesting speculation about high energy phenomena originating at the center of the Milky Way. The analysis tools applied both in the original discovery and for testing subsequent physical hypotheses leave room for improvement, particularly in terms of moving from qualitative to quantitative assessment of the results. Here we describe one potential approach leveraging Haar wavelets to capture the expected multi-scale nature of structure in all-sky gamma-ray maps, coupled with statistics to quantify the tradeoffs between detection of detail and "false detections" resulting from photon counting noise.
