In this paper, we introduce a feature preserving surface reconstruction algorithm to produce a high fidelity triangulated mesh from an input point set. The concept of local Delaunay triangulation is applied to speed up the reconstruction procedure and to preserve features. The proposed algorithm has running time complexity of O(nlogn), where n is the number of points. Additionally, the local Delaunay triangulation improves the memory efficiency of the proposed method compared to global Delaunay-based methods. We introduced the concept of minimum circum-radius triangle to select the prime-triangle in the local Delaunay mesh. Based on the local projection and minimum circum-radius, our triangle selection criteria ensures the output mesh will be free from non-manifold edges and fold-over triangles. On top of that, we have provided the theoretical correctness of the proposed algorithm with the assumption of ϵ-sampling model. The experimental results show that the proposed method is capable of producing high fidelity meshes from large real-world non-uniform data. We also show the effectiveness of the proposed method compared to the state-of-the-art methods in terms of visual and quantitative analysis. © 2021 Elsevier B.V.