We proposed an integration method based on two-dimensional Taylor expansion, in which the terms in the first two orders were used to approximate the height difference between two points. Combining the approximation relation between the first-order term and the second-order one, we derived a new difference operator, which was applied to the zone method. In addition, we proposed a non-iterative method to calculate non-matrix height. Numerical results show that when the sampling points are distributed in the barrel and pillow shapes, the reconstruction accuracy of the proposed method is 0.1-0.2 μm, which is better than that of other methods (0.5-0.7 μm). For the reconstruction of non-matrix data with 29236 sampling points, the calculation time of our method is about 0.1 s, significantly shorter than that of the iterative methods, and it is much shorter as the number of sampling points increases.