["Geographical Analysis, Volume 58, Issue 2, April 2026. ", "\nABSTRACT\nSpatial interpolation is a crucial task in geography. As perhaps the most widely used interpolation methods, geostatistical models‐such as Ordinary Kriging (OK)‐assume spatial stationarity, which makes it difficult to capture the nonstationary characteristics of geographic variables. A common solution is trend surface modeling (e.g., Regression Kriging, RK), which relies on external explanatory variables to model the trend and then applies geostatistical interpolation to the residuals. However, this approach requires high‐quality and readily available explanatory variables, which are often lacking in many spatial interpolation scenarios‐such as estimating heavy metal concentrations underground. This study proposes a Focal Feature Regression Kriging (FFRK) method, which automatically extracts geospatial features to construct a regression‐based trend surface without requiring external explanatory variables. We conducted experiments on the spatial prediction of three heavy metals in a mining area in Australia. In comparison with 17 classical interpolation methods, the results indicate that FFRK, which relies solely on extracted geospatial features, consistently outperforms both conventional Kriging techniques and machine learning models that depend on explanatory variables. This approach effectively addresses spatial nonstationarity while reducing the cost of acquiring explanatory variables, improving both prediction accuracy and generalization ability.\n"]