cs.CVOct 20, 2026
SaveSAIL: A Learnable Generalisation of Spatial Autocorrelation over Graph-Structured Data
Organizations: Predictive Systems in Biomedicine (PRISM) Lab, Department of Computer Science, University of Warwick, Coventry CV4 7AL, United Kingdom
Abstract
In vision and scientific domains involving spatial data, a key question is whether node feature values vary independently across space or form coherent spatial structures. Classical spatial autocorrelation statistics such as Moran's and Geary's quantify this phenomenon but operate on a single scalar variable with fixed neighbourhood weights, limiting their use for high-dimensional graph-like data. We introduce SAIL (Spatial Autocorrelation Index Learner), a differentiable framework that generalises spatial autocorrelation into a learnable operator by jointly learning feature projections defining candidate scalar fields and inter-node neighbourhood affinities. SAIL produces interpretable local and global spatial coherence descriptors, enabling supervised prediction and unsupervised discovery of spatial autocorrelation modes and hotspots. We show that SAIL exactly recovers classical spatial statistics under appropriate parametrisations and, for tasks driven by latent spatial autocorrelation, can approximate the Bayes-optimal predictor with an inductive bias matched to such tasks. Experiments on synthetic data and large-scale biomedical graphs demonstrate statistical correctness and discovery of spatial coherence with applications across vision, geospatial analysis, and biomedicine. The code repository is available at: https://github.com/pkeller00/SAIL.
Figures & tables
Fig. 1: Overview of the SAIL framework. (i) SAIL takes graph-structured data whose nodes carry -dimensional features. Inside the SAIL module, a learnable feature projection maps node features to scalar fields (modes), and a learnable neighbourhood affinity assigns strong or weak weights to node pairs. Each projected field is combined with these affinities to give node-wise local spatial autocorrelation (SAC) scores for every mode, which can also be aggregated into an -dimensional global SAC descriptor per graph. (ii) SAIL modes can be learned through supervised or unsupervised optimisation. In the unsupervised setting, SAIL identifies feature projections that maximise spatial autocorrelation (SAC), revealing which combinations of features exhibit the strongest spatial coherence across the dataset. Conversely, in supervised learning, global SAC descriptors are used to train a downstream predictor, enabling SAIL to identify spatial autocorrelation patterns associated with the target. In both settings, the learned feature projections provide interpretable combinations of features that define each mode. Finally, a formal statistical framework enables the detection of significant hot and cold spots on the learned modes.
Fig. 2: SAIL recovers classical SAC on synthetic graphs with known structure. Left: a graph with unstructured, independently sampled node values (colour bar shows feature value), for which SAIL correctly reports near-zero global SAC (Moran’s , Geary’s ) and classifies most nodes as non-significant under hotspot analysis. Right: a spatially homogeneous graph with a smooth imposed gradient, for which SAIL correctly detects strong positive autocorrelation (Moran’s , Geary’s ) and hotspot analysis correctly localises High–High and Low–Low clusters matching the imposed gradient.
Fig. 3: Results demonstrating distinction of SAIL from traditional message-passing architectures. In a regression benchmark where models predict Moran’s of heterogeneous test graphs, SAIL achieves significantly higher Pearson correlation between predicted and true values than classical message-passing networks with comparable or larger parameter counts ( , box plots over 1000 train–validation–test runs).
| Region/DICE Score | SAIL | Univariate | Mean | MULTISPATI | Hotspot |
| Nucleolus | |||||
| Endosomes/Golgi 1 | |||||
| Cell periphery 1 | |||||
| Cell periphery 2 |
TABLE I: SAIL segmentation performance on the 4i dataset. We compare SAIL against classical Moran’s computed on the best-performing individual protein among the proteins with the highest global Moran’s (Univariate), the multivariate extension of Moran’s obtained by averaging the univariate Moran’s statistic across all proteins (Mean), MULTISPATI-PCA [ 15 ] and Hotspot [ 21 ] . For each compartment, we report the Dice score between LISA-style hotspot mask and region annotation.
Fig. 4: Unsupervised SAIL Modes recover known subcellular compartments in the 4i spatial proteomics dataset. (a) Two data-driven SAIL Modes are visualised over the field of view. SAIL Mode 1 (showing the highest Dice score with Cell periphery 2) is defined by positive weights on p-MEK, TSG101 and Climp63 and a negative weight on ISG15. SAIL Mode 8 (showing the highest Dice score with the Nucleolus) is defined by positive weights on GSK3A/B, CAV1 and SARA and negative weights on PCNA, FBL and LC3B. (b) Visualisations of LISA High–High and Low–Low scoring locations for both SAIL Modes, which were used for defining spatial regions. We also show the expression of the most important proteins for each SAIL Mode and the reference annotation mask.
| Model | AUC (95% CI) | Mann–Whitney ( -value) |
| SAIL (Mode 18) | 0.86 (0.64–1.00) | 0.005 |
| Favourable Niches | 0.78 (0.55–0.95) | 0.021 |
| PD-L1+ (untreated slide) | 0.78 (0.49–1.00) | 0.072 |
| PD-L1+ (all slides) | 0.69 (0.42–0.92) | 0.19 |
TABLE II: Prediction of response to PD-1/PD-L1 blockade in the cSCC cohort. We report discriminative performance as the area under the ROC curve (AUC) with 95% confidence interval, together with the two-sided Mann–Whitney U p-value for separating responders from non-responders. SAIL is compared against the Favourable Niches classifier and pre-treatment PD-L1 status [ 22 ] .
Fig. 5: Association of an unsupervised SAIL Mode with response to PD-1/PD-L1 blockade. (a) Per-patient global Geary-style dissimilarity SAIL score for SAIL Mode 18, compared between responders and non-responders. Each dot represents a patient; responders show a higher score ( ), suggesting greater local cell heterogeneity. (b) Gene weights defining SAIL Mode 18. (c) Visualisation of responder and non-responder patients. We show the local and Geary-style dissimilarity score per cell.
| Metric | SAIL | Univariate | Mean | MULTISPATI | Hotspot |
| C-Index | |||||
| RMST (days) | |||||
| HR (95% CI) | (1.45–18.89) | (0.15–1.31) | (0.98–10.08) | (1.03–10.56) | (1.36–17.59) |
| Log-rank -value |
TABLE III: Prognostic performance of SAIL on the ORION cohort. We compare SAIL against classical Geary’s computed on individual protein markers (Univariate) and the multivariate extension of Geary’s obtained by averaging the univariate Geary’s (Mean), MULTISPATI-PCA [ 15 ] , and Hotspot [ 21 ] . Patients are stratified into groups via a median split on each score; C-Index, hazard ratio (HR) with 95% CI, restricted mean survival time difference ( RMST), and log-rank -value are reported.
Fig. 6: Unsupervised SAIL Modes link spatial coherence to outcome. (a) A SAIL Mode discovered in ORION-CRC (mIF). Bar chart shows weights. The corresponding projected scalar field is overlaid on a tissue section to visualise its spatial pattern, and the associated global SAIL scores stratify patients into distinct survival groups. (b) SAIL learns a scalar field on TCGA-BRCA WSI deep-feature patch graphs (Mode 1); slides are visualised by alongside hotspot clusters. The corresponding global SAIL score (median split) defines risk groups with significant Kaplan–Meier separation in disease-specific survival. (c) Similarly to (b), we analyse the TCGA-CRC cohort.
| Model | AUC | AP |
| All features ( ) | ||
| SAIL | ||
| MesoGraph | ||
| PINS | ||
| Max-MIL | ||
| Interpretable features only ( ) | ||
TABLE IV: Performance of supervised SAIL for mesothelioma subtype prediction [ 37 ] , using all cell features (top) or only the interpretable features (bottom). We report AUC and average precision (AP) with their standard deviation ( ) across cross-validation folds. Comparisons include graph-based state of the art (MesoGraph) and standard multiple-instance learning pipelines (PINS, Max-MIL).
Fig. 7: Supervised SAIL discriminates mesothelioma subtype from the spatial coherence of interpretable cell features (model trained on the interpretable features only). (a) loadings of the two SAIL Modes with the largest weights in the prediction head (Mode 11, weight ; Mode 12, weight ; both towards non-epithelioid), shown for the features with the largest absolute loading. (b) Local predicted value and local Geary-style SAIL scores for these modes, shown for a representative epithelioid and non-epithelioid sample.
Appendix figures & tables29 assets
Supplementary material from the paper’s appendix.
Appendix
| Compartment | Mode | LISA cluster | Dice |
| Nucleolus | Low–Low | ||
| Cell periphery 2 | Low–Low | ||
| Endosomes/Golgi 1 | High–High | ||
| Cell periphery 1 | Low–Low | ||
| Nucleus | High–High | ||
| ER/mitochondria 1 | High–High |
TABLE V: Best-matching unsupervised SAIL Mode per annotated subcellular compartment on the 4i dataset. For each compartment, we report the mode whose LISA hotspot mask (High–High or Low–Low) attains the highest Dice overlap with the expert annotation.
| Compartment | SAIL | Univ. | Mean | MULTISPATI | Hotspot |
| Nucleolus | (DAPI) | ||||
| Cell periphery 2 | (CRT) | ||||
| Endosomes/Golgi 1 | (SARA) | ||||
| Cell periphery 1 | (CRT) | ||||
| Nucleus | (PCNA) | ||||
| ER/mitochondria 1 | (ABCD3) |
TABLE VI: Full per-compartment Dice comparison on the 4i dataset, extending the four-compartment summary of the main text to all annotated compartments. SAIL is the best-matching unsupervised mode (Supplementary Table V ); Univariate (Univ.) is classical Moran’s with LISA computed on the single best-performing individual protein (named in parentheses) among the proteins with the highest global Moran’s ; Mean is the traditional multivariate extension of Moran’s , i.e. the mean of the univariate Moran’s statistics across all proteins; MULTISPATI is MULTISPATI-PCA [ 1 ] with local Moran’s LISA computed on the axes and Hotspot is Hotspot [ 2 ] with local Moran’s LISA computed on the scores of its protein modules, all using the same log-distance kNN graph. For each compartment we report the Dice overlap between the derived High–High/Low–Low hotspot mask and the expert annotation.
Fig. 8: All unsupervised SAIL Modes learned on the 4i dataset. For each mode we show the learned projection weights over the protein markers (bar chart, as in Figure 4 ), alongside the corresponding High–High and Low–Low LISA cluster maps overlaid on the tissue. Together these define, for every mode, both the molecular program it emphasises and the spatial regions where that program is coherently high (High–High) or low (Low–Low). The best-matching mode per annotated compartment is quantified in Supplementary Table V .
Fig. 9: All proteins in the 4i dataset. For each protein we show its localised expression, to help understand its expected spatial distribution and how it contributes to SAIL Modes that overlap with known subcellular compartments. We also show the ground-truth annotation masks for all compartments.
| Mode | Med. | Med. | AUC | ||
| 18 | 0.692 | 0.556 | 0.864 | 95 | 0.005 |
| 11 | 0.589 | 0.554 | 0.809 | 89 | 0.018 |
| 4 | 0.588 | 0.513 | 0.791 | 87 | 0.027 |
| 8 | 0.634 | 0.541 | 0.782 | 86 | 0.032 |
| 14 | 0.638 | 0.549 | 0.773 | 85 | 0.038 |
| 7 | 0.611 | 0.519 | 0.764 | 84 | 0.045 |
TABLE VII: SAIL niche discrimination of PD-1/PD-L1 blockade response (see Supplementary Section H ). Each row is one SAIL Mode, scored at the subject level across 11 responder and 10 non-responder patients and ranked by discriminative AUC. Med. and Med. are the median per-subject global Geary-style SAIL scores among responder ( ) and non-responder ( ) patients, respectively. AUC is the area under the ROC curve for separating the two groups by SAIL score; and are the Mann–Whitney–Wilcoxon rank-sum statistic and its two-sided -value.
Fig. 10: Gene weights defining all unsupervised SAIL Modes learned on the cSCC PD-1/PD-L1 CosMx cohort. For each mode, bars show the top genes by absolute loading in the learned projection, with positively and negatively weighted genes coloured separately; these weights specify the gene combination whose spatial Geary-style dissimilarity each mode minimises.
| Mode | C-Index | RMST | |
| 1 | 0.68 | 213 | |
| 2 | 0.66 | 213 | |
| 3 | 0.52 | 288 | |
| 4 | 0.72 | 502 | |
| 5 | 0.68 | 482 | |
| 6 | 0.51 | -153 |
TABLE VIII: Survival association (overall survival) of all unsupervised SAIL Modes in ORION-CRC. Patients are split at the median of each mode’s global SAIL score. C-Index, log-rank -value and RMST (restricted mean survival time of the low- minus the high-score group, in days) is reported.
| TCGA-BRCA | TCGA-CRC | |||||||
| Mode | C-Index | HR (95% CI) | RMST | C-Index | HR (95% CI) | RMST | ||
| 1 | 0.64 | 2.39 (1.28–4.46) | 386 | 0.59 | 1.45 (0.91–2.30) | 102 | ||
| 2 | 0.55 | 1.26 (0.69–2.28) | 195 | 0.60 | 1.56 (0.98–2.48) | 653 | ||
| 3 | 0.58 | 1.37 (0.76–2.47) | 159 | 0.54 | 1.14 (0.72–1.80) | 107 | ||
| 4 | 0.57 | 1.80 (0.98–3.28) | 302 | 0.56 | 1.55 (0.98–2.47) | 365 | ||
| 5 | 0.54 | 1.13 (0.62–2.07) | 104 | 0.55 | 1.74 (1.09–2.79) | 689 | ||
TABLE IX: Survival association of all unsupervised SAIL Modes in TCGA-BRCA (disease-specific survival) and TCGA-CRC (overall survival). Patients are split at the median of each mode’s global SAIL score. C-Index is reported as ; is the log-rank -value; HR (95% CI) compares the high- with the low-score group; RMST is the restricted mean survival time of the low- minus the high-score group, in days.
| Covariate | C-Index | Log-rank -value | HR (high vs low) | RMST (low–high, days) |
| Mean aggregate | 0.68 | 0.04 | 3.15 | 362 |
| CD68 | 0.61 | 0.05 | 3.00 | 643 |
| CD4 | 0.69 | 0.13 | 2.28 | 385 |
| Grade | 0.61 | 0.13 | 2.20 | 393 |
| E-cadherin | 0.59 | 0.14 | 2.22 | 440 |
| CD3e | 0.63 | 0.15 | 2.18 | 460 |
TABLE X: Baseline survival associations in the ORION cohort using inverse-distance spatial weighting to compute univariate global spatial autocorrelation using Geary’s . For the multivariate case, we report the mean of the univariate spatial autocorrelations across proteins (Mean aggregate). Results are also compared with immune groups (IFM1–4) defined in the original ORION study.
Fig. 11: Unsupervised SAIL Modes link spatial coherence to clinical outcome. A SAIL Mode was identified in the BRAF-CRC cohort by minimising Geary’s over one-hot encoded cell types (cancer epithelial cells, normal epithelial cells, fibroblasts, lymphocytes, neutrophils, macrophages, endothelial cells, and myocytes). The bar chart shows the learned weights assigned to each cell type. Global SAIL scores derived from this mode stratify patients into distinct survival groups.
| Covariate | C-Index | Log-rank -value | HR (high vs low) | RMST (low–high, months) |
| Mean aggregate | 0.59 | 0.01 | 1.58 | 10.4 |
| Endothelial cell | 0.58 | 0.02 | 1.57 | 11.8 |
| Lymphocyte | 0.59 | 0.06 | 1.42 | 7.2 |
| Cancer epithelial cell | 0.57 | 0.06 | 1.42 | 9.2 |
| Fibroblast | 0.57 | 0.07 | 1.41 | 7.5 |
| Normal epithelial cell | 0.55 | 0.26 | 1.23 | 0.9 |
TABLE XI: Baseline survival associations in the BRAF-CRC cohort using inverse-distance spatial weighting to compute univariate global spatial autocorrelation using Geary’s on each one-hot cell-type channel. The multivariate covariate (Mean aggregate) corresponds to the mean of the univariate cell-type spatial autocorrelations.
Fig. 12: Zoomed-in regions of SAIL outputs. (a) Patch graphs for TCGA-BRCA and TCGA-CRC. We show two zoomed-in regions per patient, including the score and LISA clusters for each graph node. (b) Cell graphs for BRAF-CRC and ORION-CRC. We show two zoomed-in regions per patient, including the score and LISA clusters for each graph node.
Fig. 13: 16 learned SAC descriptors (trained using interpretable features) on the mesothelioma dataset used in Section J . For each descriptor we report the three features with the most positive and negative loadings. Here “nbhd” denotes a neighbourhood average.
| Train step (ms) | Inference (ms) | Peak mem. (GB) | Max. | ||
| Linear | none | ||||
| Orthogonal | none | ||||
| Linear | shared-map | ||||
| Orthogonal | shared-map | ||||
| Linear | query–key | ||||
| Orthogonal | query–key |
TABLE XII: Computational cost of each SAIL variant at the reference configuration ( , , , ) on one NVIDIA A100-80GB GPU; mean standard deviation over repeats. Max. is the largest graph size in the sweep that fitted in GPU memory.
Fig. 14: Computational scaling of SAIL in its default configuration (orthogonal , no learned feature affinity, Moran-style objective with decorrelation and sparsity regularisation) on synthetic spatial data (Supplementary Section K-B ). Each panel varies one factor with the others held at nodes, features, modes and neighbours. Train step : one full-graph forward pass, backward pass and Adam update (one epoch); Inference : one forward pass of the trained model. (a) Runtime versus number of nodes : a fixed, -independent cost dominates up to . (b) Runtime versus feature dimension : training time grows approximately as , owing to the matrix exponential of the Stiefel parametrisation, while inference is essentially independent of . (c) Runtime versus number of modes and (d) versus neighbourhood size : both have only a minor effect. (e) Peak GPU memory versus grows linearly. The dashed line marks GPU capacity (NVIDIA A100, GB); marks the first size that did not fit, joined by a dotted line to the largest size that did. Points are means over repeats with standard deviation bands; legend exponents are log–log least-squares slopes. All axes are logarithmic.
Fig. 15: Computational cost of SAIL’s optional components. Top row (a–c): linear ; bottom row (d–f): orthogonal (Stiefel) , the default. Colours denote the feature affinity : none (blue; attention from spatial distance only), shared-map (orange) and query–key (green) (Supplementary Section K-B ). Other factors are held at , , , . (a, d) Time per training step versus feature dimension : nearly flat with the linear map, approximately cubic with the orthogonal map. (b, e) Peak GPU memory versus number of nodes : linear for all variants, but both learned affinities require approximately more memory, because they gather the -dimensional features of every neighbour ( ), and therefore exceed GPU memory at . (c, f) Peak GPU memory versus neighbourhood size : linear in for the learned affinities. The shared-map and query–key affinities have near-identical memory (and near-identical runtime with the linear map), since both are dominated by the same neighbour-feature gather; the shared-map curve is drawn wider underneath so that both remain visible. marks the first configuration that did not fit in GPU memory. The dashed line marks GPU capacity. Means over repeats with standard deviation bands; all axes are logarithmic.
Fig. 16: Additional synthetic tests of SAIL. (a) Toy graph examples, coloured by both LISA clusters and the projected scalar field . Top-left: a smooth, organically varying spatial field generated with Perlin noise, designed to produce complex but strongly positive local spatial autocorrelation. Top-right: a checkerboard graph with a sharply alternating high-contrast pattern that induces strong negative spatial autocorrelation through repeated adjacent sign flipping. Bottom-left: a plateau graph with a two-level step surface with an abrupt spatial boundary, producing clustered positive autocorrelation within each region and a sharp discontinuity at the divide. Bottom-right: concentric rings with an oscillating radial wave pattern that alternates between positive and negative spatial dependence as distance increases from the centre. (b) A synthetic 3D point set, coloured by both LISA clusters and the projected scalar field . It is generated by sampling a grid and assigning heights using 2D Perlin noise, giving a smooth synthetic terrain; the resulting points are used to build a -nearest-neighbour graph for SAIL analysis.
| Category | Setting | Mean Dice |
| 0 | 0.42 | |
| 1 † | 0.44 | |
| 5 | 0.42 | |
| 10 | 0.36 | |
| 100 | 0.32 | |
| 0 | 0.44 |
TABLE XIII: Ablation study of SAIL components on the 4i dataset (Supplementary Section M ). For each setting we report the best Dice per subcellular region (maximum over the 20 modes and High–High/Low–Low LISA clusters), averaged over the 10 regions. Each block varies one component while holding all others at their default setting (indicated via † ).
| Category | Setting | AUC (95% CI) | -value |
| 0 | 0.84 (0.62–1.00) | 0.010 | |
| 1 † | 0.86 (0.64–1.00) | 0.005 | |
| 5 | 0.82 (0.60–0.98) | 0.015 | |
| 10 | 0.75 (0.49–0.95) | 0.062 | |
| 100 | 0.78 (0.55–0.96) | 0.032 | |
| 0 | 0.84 (0.63–0.99) | 0.010 |
TABLE XIV: Ablation study of SAIL components on the cSCC PD-1/PD-L1 blockade cohort. For each setting, we report the mode that best separates responders from non-responders (highest AUC over the modes): its AUC with a bootstrap confidence interval (patients resampled with replacement, resamples) and the Mann–Whitney -value. Each block varies one component while holding all others at their default setting (indicated via † ).
| Category | Setting | C-Index | -value | RMST |
| 0 | 0.67 | 0.001 | 365 | |
| 1 † | 0.64 | 0.005 | 386 | |
| 5 | 0.60 | 0.046 | 243 | |
| 10 | 0.60 | 0.464 | 96 | |
| 100 | 0.60 | 0.123 | 188 | |
| 0 | 0.62 | 0.020 | 347 |
TABLE XV: Ablation study of SAIL components on the TCGA-BRCA cohort ( patients; Supplementary Section M ). For each setting we report the SAIL Mode with the highest C-Index; patients are split into high/low groups via a median split on that mode’s score, and we report the log-rank -value and RMST (days, 10-year horizon) between groups. Each block varies one component while holding all others at their default setting (indicated via † ).
| Category | Setting | C-Index | -value | RMST |
| 0 | 0.61 | 0.015 | 474 | |
| 1 † | 0.60 | 0.057 | 653 | |
| 5 | 0.60 | 0.091 | 178 | |
| 10 | 0.57 | 0.320 | 112 | |
| 100 | 0.57 | 0.021 | 373 | |
| 0 | 0.60 | 0.002 | 807 |
TABLE XVI: Ablation study of SAIL components on the TCGA-CRC cohort ( patients; Supplementary Section M ). For each setting we report the SAIL Mode with the highest C-Index; patients are split into high/low groups via a median split on that mode’s score, and we report the log-rank -value and RMST (days, 10-year horizon) between groups. Each block varies one component while holding all others at their default setting (indicated via † ).
Fig. 17: All SAIL Modes learned on the Global Macro Database. For each mode, bars show the top macroeconomic indicators by absolute loading in the learned projection, with positively and negatively weighted indicators coloured separately; these weights specify the indicator combination whose spatial Geary-style dissimilarity each mode minimises.
Fig. 18: Temporal trends of all SAIL Modes learned on the Global Macro Database. For each mode, the line shows the global Geary-style SAIL score of the mode’s projected field on each yearly country graph from 2000 to 2023; lower values indicate that neighbouring countries are more similar along that mode.
Fig. 19: Unsupervised SAIL for COVID-19 trend analysis. Here we show a line graph that reports the global Moran’s of new COVID-19 cases (y-axis) against a given week (x-axis).
Fig. 20: A choropleth map for the world coloured by LISA clustering. Interestingly, we notice that across time, COVID-19 hotspots tend to form in Western and Central European countries, reflecting a consistent pattern of spatially concentrated High–High clusters and regional synchronisation of outbreaks.
Fig. 21: SAIL tracks the mixing of two gases. (i) Loadings of the single SAIL Mode learned on a training run, over species identity, particle speed and four nuisance noise features. (ii) SAIL score over time (seconds) for the three scenarios (partition closed, small hole, partition removed); lines show the mean and shaded bands standard deviation over seeds. (iii) Gas states (top) and the corresponding LISA clusters of the learned mode (bottom) at selected time points.
| Model | #Params | Pearson | MSE |
| SAIL (dense) | |||
| GCN (X, FC) | |||
| GCN (X+coords, FC) | |||
| GAT (X, FC) | |||
| GAT (X+coords, FC) | |||
| MLP (X) |
TABLE XVII: Per-model results on the synthetic Moran’s regression benchmark. Each model is evaluated at its Stage-1–selected configuration over seeds; we report the mean standard deviation of test Pearson and MSE. Rows are ordered by descending Pearson correlation.
| SAIL | GCN | GAT | MLP | |
TABLE XVIII: Spatial-versus-content blend ablation. Mean test Pearson correlation (subscript: standard deviation over seeds) as the target shifts from pure spatial autocorrelation ( ) to pure feature content ( ). GCN, GAT and MLP use node coordinates (X+coords), and the GNNs use the fully connected graph (FC). SAIL is strongest at and degrades monotonically as the target loses its spatial component, while the content-recoverable baselines improve. The best model per row is in bold, with a clear crossover near – .
Fig. 22: Spatial-versus-content blend ablation. Mean Pearson correlation (shaded bands standard deviation over seeds) as a function of the blend weight , interpolating the regression target between pure spatial autocorrelation ( ) and pure feature content ( ). SAIL is strongest at and degrades monotonically as the spatial component is removed, while GCN, GAT, and MLP improve, crossing SAIL near – .