The function of enzymatic proteins is given by their ability to bind specific small molecules into their active sites. These sites can often be found in pockets on a hypothetical boundary between the protein and its environment. Detection, analysis, and visualization of pockets find its use in protein engineering and drug discovery. Many definitions of pockets and algorithms for their computation have been proposed. Kawabata and Go defined them as the regions of empty space into which a small spherical probe can enter but a large probe cannot and developed programs that can compute their approximate shape. In this article, this definition was slightly modified in order to capture the existence of large internal holes, and a Voronoi‐based method for the computation of the exact shape of these modified regions is introduced. The method first puts a finite number of large probes on the protein exterior surface and then, considering both large probes and atomic balls as obstacles for the small probe, the method computes the exact shape of the regions for the small probe. This is all achieved with Voronoi diagrams, which help with the safe navigation of spherical probes among spherical obstacles. Detected regions are internally represented as graphs of vertices and edges describing possible movements of the center of the small probe on Voronoi edges. The surface bounding each region is obtained from this representation and used for visualization, volume estimation, and comparison with other approaches. © 2019 Wiley Periodicals, Inc.