Expansion Method:
Use the trend surface function in place of a constant coefficient:
\[
y_i = f_{(lat_i, long_i)}^1 + f_{(lat_i, long_i)}^2 x_{i1} + \epsilon_i
\]
\[
\begin{array}{l}
f_{(lat_i, long_i)}^1 = \beta_{01} +\beta_{02}lat_i + \beta_{03}long_i\\
f_{(lat_i, long_i)}^2 = \beta_{11} +\beta_{12}lat_i + \beta_{13}long_i
\end{array}
\]
Geographically Weighted Regression:
Instead of fitting a regression model to all the observations, it fits to a subset of them each time based on the bandwidth defined.
It is like a spatial moving average; instead of calculating an average, we are fitting a regression model to each window.