Various lunar digital topographic models (DTMs) have been generated from the data collected from earlier and recent lunar missions. There are usually inconsistencies among them due to differences in sensor configurations, data acquisition periods, and production techniques. To obtain maximum value for science and exploration, the multi-source lunar topographic datasets must be co-registered in a common reference frame. Only such an effort will ensure the proper calibration, registration, and analysis of the datasets, which in turn will permit the full comparative and synergistic use of them. This study presents a multi-feature-based surface matching method for the co-registration of multiple lunar DTMs that incorporates feature points, lines, and surface patches in surface matching to guarantee robust surface correspondence. A combined adjustment model is developed for the determination of seven transformation parameters (one scale factor, three rotations, and three translations), from which the multiple DTMs could be co-registered. The lunar DTMs derived from the Chang’E-1, SELENE, and LRO laser altimeter data in the Apollo 15 landing area and the Sinus Iridum area are examined in this study. Small offsets were found among the Chang’E-1, SELENE, and LRO DTMs. Through experimental analysis, the developed multi-feature-based method was proven able to effectively co-register multiple lunar DTMs. The performances of the multi-feature-based surface matching method were compared with the point-based method, and the former was proven to be superior to the latter.