import numpy as np
from spatialdm.stats import rbfweight
X_loc = np.vstack([np.repeat(range(10), 10), np.tile(range(10), 10)]).T
spatial_W, KNN_connect = rbfweight(X_loc, l=1.2, n_neighbors=16)
#
import matplotlib.pyplot as plt
plt.scatter(X_loc[:, 0], X_loc[:, 1], c=spatial_W.toarray()[35], s=100)
plt.show()