Grid coverage calculations are a commonly used analysis tool for Earth-observation and communication satellite missions. A dynamic coverage calculation by a maneuverable satellite instrument over a wide pointing domain requires a fast coverage calculator. Given an arbitrary target grid of points on the Earth’s surface, the basic coverage calculation problem is to compute the access intervals of the satellite instrument with each point on the grid over the mission period. The access is computed using either the line of sight or a particular instrument field of view (FOV) geometry. Grid coverage calculations can be computationally intensive for high-resolution grids due to the large number of point containment queries needed to establish access. In this article, we introduce a preprocessing algorithm to accelerate grid coverage calculations. The target grid is first sorted into the cells of a preprocessing grid, which has equal latitude spacing and equal longitude spacing along a given latitude, with the longitude spacing increasing towards the poles to ensure roughly equal cell area. Due to the equal spacing property, the grid cells can be indexed directly using the bounding box of the instrument footprint on the Earth’s surface in latitude and longitude. Our method considers rectangular, conical, and spherical-polygon instrument profiles, which together can account for most practical instrument geometries. The preprocessing algorithm ensures that the number of point containment operations is proportional only to the average number of points inside the satellite footprint, rather than the total number of points in the target grid. This leads to dramatic runtime improvements, especially when the access area is only a small fraction of the Earth’s surface, as is typically the case for Earth-observation satellites. Results show 1-2 orders of magnitude improvement in runtime over the commonly used brute-force approach. We analyze the runtime performance with respect to the resolution of the preprocessing grid for various footprint sizes, and suggest a way to obtain the optimal resolution. To support the algorithm development, we also provide an overview of several geometric formulas for representing instrument profiles and projecting them onto the Earth’s surface using a spherical Earth model.