The volume of tissue activated (VTA) is a well-established approach to model the direct effects of deep brain stimulation (DBS) on neural tissue and previous studies have pointed to its potential clinical applications. However, the elevated computational time required to estimate the VTA with standard techniques used in biological neural modeling limits its suitability for practical use. The goal of this project was to develop a novel methodology to reduce the computation time of VTA estimation. To that end, we built a Gaussian process emulator. It combines a field of multi-compartment axon models coupled to the stimulating electric field with a Gaussian process classifier (GPC); following the premise that computing the VTA from a field of axons is in essence a binary classification problem. We achieved a considerable reduction in the average time required to estimate the VTA, under both ideal isotropic and realistic anisotropic brain tissue conductive conditions, limiting the loss of accuracy and overcoming other drawbacks entailed by alternative methods.