We need to produce figures showing the rms maps

• rms maps
• rms histograms
• for all regions in DR1

• v_LSR, sigma_v, Tk, Tex, para-N(NH3). Others?
• generate for all regions

@rfriesen We need to include the script that is used to generate the images with Herschel data and the footprint of the GAS data, plus the YSO catalogues.

• add script that create figures
• include YSO catalogues (edited if needed)
• check that YSO catalogues are fully plotted
• Add Herschel data to dataverse
• Add figures to repo

• integrated intensity for all transitions
• plot for all regions

For aesthetics, it may be a good idea to clean up the images used in contour maps to kill off the noise spikes. One way to do this might be binary opening on the parameter maps or on a masked map that is masked at the threshold of the lowest contour. This just eliminates small regions in the map.

from scipy.ndimage import binary_opening
selem = np.array([0,1,0],[1,1,1],[0,1,0])
mask = binary_opening(map > LowestContour,selem)
MaskedMap = mask*map