In this work, Cox processes and Convolutional Neural Net classifiers (CNNs) are used to estimate the number of instances of an object in an image. Poisson processes are well suited to model events that occur randomly in space, such as the location of objects in an image or the enumeration of objects in a scene. Doubly Stochastic Poisson (or Cox) processes offer increased flexibility, but their computational complexity and storage requirements do not easily scale with image size, typically requiring O(n 3) computation time and O(n 2) storage, where n is the number of pixels in an image. To mitigate this problem we employ the Kronecker algebra, which takes advantage of direct product structures. As the likelihood is non-Gaussian, the Laplace approximation is used for inference, employing the conjugate gradient and Newton's method. Our approach has then close to linear performance, requiring only O(n 3/2) computation time and O(n) memory. The proposed algorithm consists of selecting a subset of bounding boxes in the image domain, and querying them for the presence of the object of interest by running a pre-trained CNN classifier. The resulting observations are then aggregated and a posterior distribution