A framework with a combination of the radial basis functions (RBFs) method and the least-squares integration method is proposed to improve the integration process from gradient to shape. The principle of the framework is described, and the performance of the proposed method is investigated by simulation. Improvement in accuracy is verified by comparing the result with the usual RBFs-based subset-by-subset stitching method. The proposed method is accurate, automatic, easily implemented, and robust and even works with incomplete data.