Ghost imaging through turbulent atmospheres are theoretically studied. Based on the extended Huygens-Fresnel integral, we obtain an analytical imaging formula. The ghost image can be viewed as the convolution of the original object and a point-spread function (PSF). The imaging quality is determined by the size of the PSF. Increasing the turbulence strength and propagation distance, or decreasing the source size, will increase the size of the PSF, and lead to the degradation of the imaging quality.