We introduce an approach for exploring eigenvector localization phenomena for a class of (unbounded) selfadjoint operators. More specifically, given a target region and a tolerance, the algorithm identifies candidate eigenpairs for which the eigenvector is expected to be localized in the target region to within that tolerance. Theoretical results, together with detailed numerical illustrations of them, are provided that support our algorithm. A partial realization of the algorithm is described and tested, providing a proof of concept for the approach.