While potteries are common artifacts excavated in archaeological sites, the restoration process relies on manual cleaning and reassembly of shattered pieces. Since the number of possible 3D configurations is considerably large, the exhaustive manual trial may result in abrasion on fractured surfaces and even failure to find the correct matches. As a result, many recent works suggest virtual reassembly from 3D scans of the fragments. The problem is challenging in the view of the conventional 3D geometric analysis, as it is hard to extract reliable shape features from the thin break lines. We propose to optimize for the global configuration by combining geometric constraints with information from noisy shape features. Specifically, we enforce bijection and continuity of sequence of correspondences given estimates of corners and pair-wise matching scores between multiple break lines. We demonstrate that our pipeline greatly increases the accuracy of correspondences, resulting in stable restoration of 3D configurations from irregular and noisy evidence.