We introduce a novel approach for solving the problem of identifying regions in the framework of Method of Regions by considering singularities and the associated Landau equations given a multi-scale Feynman diagram. These equations are then analyzed by an expansion in a small threshold parameter via the Power Geometry technique. This effectively leads to the analysis of Newton Polytopes which are evaluated using a Mathematica based convex hull program. Furthermore, the elements of the Gröbner Basis of the Landau Equations give a family of transformations, which when applied, reveal regions like potential and Glauber. Several one-loop and two-loop examples are studied and benchmarked using our algorithm which we call ASPIRE. © 2019, The Author(s).