Modal cross coupling frequently occurs in modal approaches from wavefront gradient data such as lateral shearing measurement through Zernike circle polynomials, since the gradients of Zernike circle polynomials are not orthogonal. We use a modal approaches incorporating the Gram matrix, using the orthogonality of angular derivative of m≠0 modes with respect to weight function w(ρ) = ρ (polar coordinates), and the orthogonality of radial derivative of m = 0 modes with respect to weight function w(ρ) = ρ(1-ρ2) (polar coordinates). The Gram matrix method needs no auxiliary vector functions. The Zernike coefficients can be obtained with no modal cross coupling. The simulation results are given, which indicate that the modal cross coupling is avoided by using Gram matrix method. This method can be easily extended to annulus, and the coefficients of Zernike annular polynomials with no modal cross coupling can be obtained.