Three-dimensional (3D) kinetic maps of the Milky Way interstellar medium are an essential tool in studies of its structure and of star formation. We aim to assign radial velocities to Galactic interstellar clouds now spatially localized based on starlight extinction and star distances from Gaia and stellar surveys. We developed an automated search for coherent projections on the sky of clouds isolated in 3D extinction density maps on the one hand, and regions responsible for CO radio emissions at specific Doppler shifts on the other hand. The discrete dust structures were obtained by application of the Fellwalker algorithm to a recent 3D extinction density map. For each extinction cloud, a technique using a narrow sliding spectral window moved along the contour-bounded CO spectrum and geometrical criteria was used to select the most likely velocity interval. We applied the new contour-based technique to the 3D extinction density distribution within the volume encompassing the Taurus, Auriga, Perseus, and California molecular complexes. From the 45 clouds issued from the decomposition, 42 were assigned a velocity. The remaining structures correspond to very weak CO emission or extinction. We used the non-automated assignments of radial velocities to clouds of the same region presented in Paper I (Ivanova et al., 2021A&A...652A..22I) and based on KI absorption spectra as a validation test. The new fully automated determinations were found to be in good agreement with these previous measurements. Our results show that an automated search based on cloud-contour morphology can be efficient and that this novel technique may be extended to wider regions of the Milky Way and at larger distance. We discuss its limitations and potential improvements after combination with other techniques.

Cone search capability for table J/A+A/677/A107/table2 (3D position and assigned LSR velocity of each channel barycenter of the clumps)