these maps is $\sim$12$''$. Using the Bayesian [*PPMAP*]{} procedure operated on the [*Herschel*]{} data at wavelengths of 70, 160, 250, 350 and 500 $\mu$m [@marsh15; @marsh17], the [*Herschel*]{} temperature and column density maps were produced for the [*EU-funded ViaLactea project*]{} [@molinari10b].

In Figures \[fig2\]a and \[fig2\]b, we display the [*Herschel*]{} temperature and column density ($N(\mathrm H_2)$) maps of our selected target area, respectively. In both the [*Herschel*]{} maps, a broken box shows the area studied by @kohno18. The [*Herschel*]{} temperature map shows the existence of embedded structures with T$_\mathrm{d}$ $\sim$17–19 K in the system, and is overlaid with a temperature contour of 18.6 K, allowing us to depict filamentary structures (see labels fs1 and fs2 in Figure \[fig2\]a). The [*Herschel*]{} temperature map also exhibits extended features with T$_\mathrm{d}$ $\sim$21–29 K in the direction of the previously known H[ii]{} regions (see W33 A (G012.90$-$0.24, G012.90$-$0.26), W33 B (G012.68$-$0.18), W33 Main (G012.81$-$0.19), MIR bubbles, and G013.210$-$0.144).

In Figure \[fig2\]b, the column density contour (in white) is also overlaid on the $N(\mathrm H_2)$ map with a level of 3.35 $\times$ 10$^{22}$ cm$^{-2}$, indicating the presence of materials with high column densities. Using this $N(\mathrm H_2)$ contour level, we trace the boundary of the 