The two-point source longitudinal resolution of three-dimensional integral imaging depends on several factors including the number of sensors, sensor pixel size, pitch between sensors, and the lens point spread function. We assume the two-point sources to be resolved if their point spread functions can be resolved in any one of the sensors. Previous studies of integral imaging longitudinal resolution either rely on geometrical optics formulation or assume the point spread function to be of sub-pixel size, thus neglecting the effect of the lens. These studies also assume both point sources to be in focus in captured elemental images. More importantly, the previous analysis does not consider the effect of noise. In this manuscript, we use the Gaussian process-based two-point source resolution criterion to overcome these limitations. We compute the circle of confusion to model the out-of-focus blurring effect. The Gaussian process-based two-point source resolution criterion allows us to study the effect of noise on the longitudinal resolution. In the absence of noise, we also present a simple analytical expression for longitudinal resolution which approximately matches the Gaussian process-based formulation. Also, we investigate the dependence of the longitudinal resolution on the parallax of the integral imaging system. We present optical experiments to validate our results. The experiments demonstrate agreement with our Gaussian process-based two-point source resolution criteria.