For a wavefront tested by Shack-hartmann wavefront sensor, the zonal integration method is often chosen by researchers to solve the reconstruction problem. But it has shown an unacceptable result when the phase derivative data is distributed on an unconnected domain, the obtained wavefront will contain different piston error on each subdomain. Therefore, a new zonal wavefront estimation algorithm is proposed to deal with this drawback, which uses a simultaneous fit method to correct piston error of each subdomain. The validity of the algorithm is verified by a numerical simulation and experimental results. ? 2019 SPIE.