In this work, we have introduced an orbital pair selection that offers a fully-automatized way to partition orbital pairs into sets which are treated at different accuracy levels within multi-level local coupled-cluster approaches. The selection scheme is based on orbital maps between structures on the reaction coordinate that are constructed by comparing the localization of orbitals. As this comparison itself requires atom maps, we have developed an atom mapping algorithm minimizing the number of formed and broken bonds while considering the orientation of atoms through their dihedral angles. The performance of this algorithm was demonstrated for 233 reactions of five subsets from the BH9 set. For one case, an unexpectedly high number of formed and broken bonds was found and could be traced back to an incorrect geometry. Therefore, the atom mapping algorithm is valuable in its own right as a tool for the verification of molecular geometries in given reaction path.