Reconstructed maps of the opacity-corrected brightness temperatures, $T_{\text{b}*}^{\text{HI}}$ (top) and $T_{\text{b}*}^{\text{CO}}$ (bottom), integrated over velocity, toward the G139 region, after application of the Gaussian decomposition algorithm ROHSA (leftmost column), together with the contributions from the individual Gaussian kinematic components (columns 2–6). The maps can be rescaled to dust intensity at 353 GHz using the best-fit values of the conversion factors derived in Sect. 3.2 (see Fig. 3).