Method, apparatus, and computer program product are provided for estimating crop pest risk and/or crop disease risk at sub-farm level. In some embodiments, a farm region is determined based on farm definition data, and input data associated with the farm region are retrieved from a plurality of data sources. The input data may include a plurality of pixel sets. In some embodiments, for each pixel set, crop risk data are determined based on the input data using one or more spatiotemporal regression models. The crop risk data may include an estimate of crop pest risk and/or crop disease risk for each pixel set. In some embodiments, the farm region is categorized into a plurality of sub-farms each defining a risk level category for that sub-farm based on the crop risk data. In some embodiments, the sub-farms are displayed as a visual heat-map, along with recommended antidote options.