Space-variant blur is encountered when objects extend beyond the isoplanatic patch associated with the intervening atmospheric turbulence. The method of phase diversity, used to estimate jointly the object and the aberrations, is generalized to accommodate turbulence-induced space-variant blur. This generalization utilizes a parametric model for the blur function that is constructed with multiple phase screens. Simulation results are presented that demonstrate the recovery of near-diffraction-limited imagery from phase-diversity imagery that has been degraded with rather severe anisoplanatism.