Color-magnitude diagrams of Gaia (upper left), SkyMapper (upper right), NSC (bottom left), and M2002 (bottom right) datasets. In each diagram, left panel: CMD overlapped with MIST evolutionary tracks of 7, 9, 12, 15, 20, 25, 32 and 40 M_⊙ and color-coded as BSG (blue), YSG (yellow), and RSG (red) phases. The regions of each type of evolved massive star are outlined by the dashed lines. Targets without errors are not shown in the CMDs, while the average photometric uncertainties are indicated when available. Right panel: selected targets for each type of massive stars with the same color convention. The RSG region is empirically extended as shown by the horizontal dashed lines. Green contours represent the number density. See text for details.