The calculation of earthwork plays a major role in the plan or design phase of many civil engineering projects, such as seashore reclamation; and thus, it has become very important to improve upon its accuracy. There have been common drawbacks to earlier methods of ground profiling, such as dialing with sharp corners or the grid points of any tow straight lines. In this paper, we prepose an algorithm for finding a terrain surface using the natural boundary conditions and the both direction spline method, which interpolates the given three-dimensional data by using spline. As a result of this study, the algorithm of the proposed two methods to estimate pit excavation volume should provide a better accuracy than Spot height, Chambers, Chen, or Lin method. Also, the mathematical model mentioned offers maximum accuracy in estimating the volume of a pit excavation.