Pith. sign in

REVIEW 2 major objections 4 minor 99 references

Harnessing the Potential of Spatial Statistics for Spatial Omics Data with pasta

T0 review · 2 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read This paper argues that classical spatial statistics, matched to whether spatial omics data form point patterns or lattices, can quantify biological structures, recapitulating known breast cancer receptor zones and invasion patterns.

desk verdict Solid, honest resource paper; the k=6 neighbourhood choice needs a sensitivity check but doesn't sink the paper. read the letter →

arxiv 2412.01561 v3 pith:HQMT23IP submitted 2024-12-02 q-bio.QM q-bio.GN

classification q-bio.QMq-bio.GN
keywords spatialomicsstatisticspointpatternslatticedataautocorrelationMoran'sIBesag'sLco-localization
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Spatial omics data come in two statistical shapes: point patterns, where each cell or transcript is a stochastic point, and lattice data, where measurements sit on fixed spots or segmented cells. The paper argues that the choice of statistical lens should follow this data representation rather than the instrument, and that classical spatial statistics can then quantify biology from local gene expression to tissue organization. Re-analyzing a breast cancer imaging dataset, local Moran's I and the Moran scatter plot identify single-, double- and triple-positive zones for three hormone-receptor genes, and Besag's L function quantifies how differently two ductal carcinoma in situ subtypes associate with invasive tumour cells. Both results recapitulate the original publication's findings and add quantitative, significance-aware detail. The argument is accompanied by pasta, a vignette that demonstrates the analyses in R and Python.

What carries the argument

The machinery is the neighbourhood relation. In lattice analysis it is the weight matrix — a matrix of connection strengths between cells or spots — and the local autocorrelation statistics computed on it, especially local Moran's I with the Moran scatter plot for classifying 'high-high', 'low-low' and heterogeneous zones. In point pattern analysis it is the r-neighbourhood, a disk of radius r around each point, with complete spatial randomness (a uniform, independent scatter of points) as the reference; Besag's L compares observed neighbour counts at each radius to that random baseline. Observation windows, homogeneity assumptions, and edge corrections determine what the functions see, and the paper shows that changing the window or the neighbourhood definition can alter interpretations.

What would settle it

Recompute the breast-cancer receptor analysis with $k=4$, $k=8$, a distance-based weight matrix, and a contiguity-based matrix; if the union of significant high-high clusters no longer delineates the same triple-positive DCIS region, then the identification depends on the chosen neighbourhood rather than on the biology.

Watch

Extended reading notes

Core claim

The central claim is that the two streams of spatial statistics — point pattern analysis, which models the stochastic process generating point locations, and lattice analysis, which treats locations as fixed and models dependence among features through neighbourhood weights — are directly applicable to spatial omics and should be selected by data modality and mark type. The paper demonstrates that imaging-based data can be represented either way: cell centroids as a point pattern, or segmented cells as an irregular lattice of expression measurements. Using a six-nearest-neighbour weight matrix, local Moran's I plus the Moran scatter plot recover triple-positive ERBB2/ESR1/PGR regions in breast cancer, and bivariate Lee's L and multivariate Geary's c add pairwise and multivariate views of spatial correlation. On the point-pattern side, Besag's L in a restricted observation window shows DCIS1 cells are spaced apart from invasive tumour cells while DCIS2 cells mix with them at chance levels, recapitulating the original findings.

Load-bearing premise

The load-bearing premise is that the six-nearest-neighbour neighbourhood defines the right scale for identifying receptor-positive regions; a different $k$ or a distance-based weight matrix could shift the triple-positive classification.

Editorial extensions

If this is right

  • Analysts can pick point-pattern or lattice methods from how the data are represented, not from the brand of technology, so imaging data can be analyzed with both streams and spot-based data can be segmented into lattices or points.
  • Local spatial autocorrelation with a Moran scatter plot turns diffuse gene-expression maps into discrete regions, such as hormone-receptor-positive tumour zones, with per-cell significance attached.
  • Besag's cross L gives a scale-resolved, quantitative readout of whether cell types cluster, mix, or repel, replacing visual inspection of co-localization.
  • The choice of weight matrix, whether contiguous, k nearest neighbours, or distance-based, changes local Moran's I values for some cells, so reporting the choice and its rationale is part of the analysis.
  • Observation scale controls biological interpretation: islet cells that look clustered in a whole field of view look randomly distributed within an islet, so scale must follow the question.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The sensitivity of local Moran's I to neighbourhood choice shown on the lung cancer example implies that the triple-positive breast-cancer regions should be checked for the same sensitivity; a sweep over $k$ or a distance-based matrix would tell whether those regions are an artifact of $k=6$.
  • Because imaging data can be represented as both a point pattern and a lattice, the same biological question could be cross-checked across the two streams, for example comparing cell-type co-localization from Besag's L with spatial autocorrelation of categorical cell-type marks.
  • Cells with no contiguous neighbours produce zero local Moran's I values, which suggests spatial statistics could double as a segmentation-quality diagnostic in imaging-based assays.
  • The paper's scale-dependent results point to a testable extension: normalization choices upstream, spatially aware versus global, may change which cells are classified as high-high, so combining normalization and neighbourhood choice in one sensitivity analysis would clarify how much of the biological readout is preprocessing-driven.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. This paper argues that classical spatial statistics, partitioned into point pattern analysis and lattice data analysis, provides a powerful and underused toolkit for spatially resolved omics data. The authors illustrate this claim through several public datasets: Xenium breast cancer data (receptor expression and tumour cell co-localization), IMC pancreas islets (scale and homogeneity effects), Visium mouse brain (local Moran's I), and CosMx lung cancer (weight matrix sensitivity). They also introduce pasta, a collection of R and Python vignettes that demonstrates the methods on real data, and discuss technical caveats such as window sampling, homogeneity assumptions, confounding between intensity and interaction, and weight matrix construction.

Significance. If the illustrative analyses are robust, the paper makes a valuable educational contribution by connecting two mature statistical literatures (point processes and lattice data) to the rapidly growing spatial omics field. The analyses are transparent, use standard estimators, and are accompanied by reproducible code, versioned source code, and a public vignette; the paper also candidly discusses scale, homogeneity, and confounding. The central risk is that the main biological demonstration, the triple-positive receptor map in breast cancer, depends on an analyst-chosen neighbourhood size that is not sensitivity-tested, and the recapitulation claim is not quantitatively anchored to the original publication. These issues are fixable and do not undermine the overall educational value of the resource.

major comments (2)
  1. [Biological applications, 'Spatial autocorrelation reveals triple positive regions in breast tumours'] The triple-positive receptor map in Figure 2B is built from univariate local Moran's I values with the neighbourhood fixed to the six nearest neighbours (k=6) of each cell, but the paper provides no sensitivity analysis for this choice on the Xenium data. This is load-bearing because the paper itself demonstrates in Figure 5D-F and Supplementary Figure S4C-D that local Moran's I values and even the presence of neighbours change when the weight matrix is constructed differently, and in the 'Definition of the neighbourhood' section it acknowledges that 'it remains to be investigated how much the construction of the weight matrix influences downstream analyses in spatial omics data.' Given that cell density varies across tissue compartments, a fixed k=6 mixes widely different physical scales, and small perturbations of the local Moran's I near the significance or high-high quadrant thresholds could shift the cells labelled triple-positive, especially at region boundaries. Please add a sensitivity analysis (e.g., k=4, 8, 10; distance-based and contiguity-based weights) and report the stability of the triple-positive regions and whether the qualitative recapitulation of Janesick et al. persists.
  2. [Biological applications, 'Spatial autocorrelation reveals triple positive regions in breast tumours'] The statement that spatial statistics 'was able to recapitulate the original findings in Janesick et al. [44]' is not accompanied by any quantitative comparison to the original annotations (e.g., overlap of the inferred triple-positive cell set with the pathologist-annotated regions), and the cross-reference to '(Figure 5)' is evidently incorrect because Figure 5 displays the Visium and CosMx analyses rather than the breast cancer receptor map. Please either provide a quantitative evaluation of the agreement with the original findings or temper the claim to a qualitative demonstration, and correct the figure citation.
minor comments (4)
  1. [Figure 3 legend and Supplementary Figure S3] The legend of Figure 3 contains the typo 'DICS 2' for 'DCIS 2', and Supplementary Figure S3 contains 'ductual carcinoma in situ' for 'ductal carcinoma in situ'; these should be corrected.
  2. [Definition of the neighbourhood is critical in lattice data analysis] The sentence beginning 'Figure 5C-D show the difference' refers to panels D, E, and F for the contiguity, 10-nearest-neighbour, and 1000-pixel-distance weight matrices; the panel reference should read 'Figure 5D-F' to match the figure.
  3. [Methods and Figure 3] The paper does not report the exact intensity threshold used to restrict the observation window in Figure 3; please provide this value in the Methods or figure caption, since the comparison in Supplementary Figure S3 is appreciated but the specific threshold remains an analyst choice.
  4. [Figure 4] In Figure 4, panels B and C have different y-axis ranges, which can make the visual comparison of clustering strength between the homogeneous and inhomogeneous K-functions misleading; consider using a common y-axis range.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation: the spatial-statistics demonstrations use literature methods on external datasets, and the authors' own packages serve only as implementations, not as load-bearing evidence.

full rationale

The paper's central demonstrations—recapitulating triple-positive receptor regions and DCIS invasion patterns in the Xenium breast cancer data of Janesick et al.—are applications of established spatial statistics (local Moran's I with Moran scatter plots, Besag's L functions) to externally published data. No quantity is fitted to a subset of data and then reported as a prediction; the k=6 nearest-neighbour choice is an analyst decision, not a fitted parameter, and the paper explicitly acknowledges that weight-matrix construction is an open question and that 'it remains to be investigated how much the construction of the weight matrix influences downstream analyses' (Section 'Definition of the neighbourhood is critical in lattice data analysis'). This is an acknowledged limitation, not a circular step. Citations to the authors' own packages (sosta, spatialFDA) are implementation references for window determination and metric computation; the biological conclusions are validated against the independent external findings of Janesick et al., so those self-citations are not load-bearing evidence. The point-pattern/lattice framing is explicitly imported from the spatial-statistics literature (e.g., Baddeley et al., Cressie), not derived from the paper's own outputs. No equation defines its target result by construction, and no fitted value is renamed as a prediction. The paper is therefore self-contained as a demonstration of existing methods, and no significant circularity is present.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

This is a perspective and resource paper, not a derivation; it introduces no new entities or fitted parameters for a predictive model. The listed free parameters are analyst-chosen neighbourhood and window settings that materially affect the illustrative results. The axioms are the standard assumptions of spatial statistics that the paper imports from the cited literature (Baddeley et al., Cressie, and others) and discusses explicitly.

free parameters (4)
  • k-nearest neighbours for Xenium local Moran's I = k = 6
    Chosen by hand for the breast cancer analysis; no sensitivity analysis reported for this dataset.
  • k-nearest neighbours for CosMx local Moran's I = k = 10
    One of three neighbourhood definitions compared in Figure 5E.
  • Distance threshold for CosMx local Moran's I = 1000 pixels (about 180 micrometers)
    Chosen as an alternative neighbourhood definition in Figure 5F.
  • Intensity threshold for restricted observation window (Xenium and IMC) = Not specified numerically
    Used to define analysis windows via the sosta package; threshold value is an analyst choice affecting L-function results.
assumptions (5)
  • domain assumption Points are realizations of a stochastic point process (event-based view).
    Underlies all point pattern analyses (Section 2.4); the paper acknowledges this may not hold in biological development.
  • domain assumption Homogeneity of the point process for homogeneous K and L functions.
    Used in Figure 3B and Figure 4B/E; the paper discusses confounding of intensity and interaction (Box Point Processes).
  • domain assumption The observation window is a representative sample (small world vs window sampling).
    Section 3.2; results depend on window choice.
  • domain assumption Tobler's first law holds for the spatial relationship (near things more related).
    Basis for spatial autocorrelation; the paper notes it is not generally applicable in compartmentalized anatomy (Step 3).
  • standard math Edge-correction methods (isotropic, translation) provide unbiased estimates.
    Used in K and Besag L calculations; not derived in the paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Harnessing the Potential of Spatial Statistics for Spatial Omics Data with pasta." pith.science (2026). https://pith.science/paper/HQMT23IP

@misc{pith2026241201561,
  author       = {Pith},
  title        = {Pith review of: Harnessing the Potential of Spatial Statistics for Spatial Omics Data with pasta},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HQMT23IP}},
  note         = {Machine review of arXiv:2412.01561}
}
read the original abstract

Spatial omics assays allow for the molecular characterisation of cells in their spatial context. Notably, the two main technological streams, imaging-based and high-throughput sequencing-based, can give rise to very different data modalities. The characteristics of the two data types are well known in adjacent fields such as spatial statistics as point patterns and lattice data, and there is a wide range of tools available. This paper discusses the application of spatial statistics to spatially-resolved omics data and in particular, discusses various advantages, challenges, and nuances. This work is accompanied by a vignette, pasta, that showcases the usefulness of spatial statistics in biology using several R packages.

Figures

Figures reproduced from arXiv: 2412.01561 by the authors.

Figure 1
Figure 1. Overview of the two streams of spatial omics data and corresponding analysis strategies. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Analysis of three receptor genes (ERBB2, ESR1 and PGR) in a Xenium human breast cancer data set [44]. A) H&E image of a breast cancer section. The black dashed line indicates the region profiled by Xenium (hand annotated). B) Moran’s scatter plot [7] highlighting the status (+/-) of each cell for the receptors ERBB2, ESR1 and PGR. Double and triple positive results were derived from the union of the univariate “high… view at source ↗
Figure 3
Figure 3. Point pattern analysis of different cell types in a Xenium human breast cancer data set [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: A) Single image of an IMC proteomics dataset showing islets in the human pancreas [27]. The points are the centroids of cells coloured by cell type category. The coloured dashed line indicates that the analysis window is set to be the entire FOV. Axes on the µm scale. …
Figure 5
Figure 5. Figure 5: Spatial autocorrelation analysis. Panel A) shows log-transformed counts of Nrgn ex￾pression in the Visium mouse coronal brain section data [1]. We note regions with no expression (white), regions with low expression (blue) and regions with high expression (violet). B) …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

99 extracted references · 46 canonical work pages

  1. [44]

    High Resolution Mapping of the Tumor Microenvironment Using Integrated Single-Cell, Spatial and in Situ Analysis

    Amanda Janesick et al. “High Resolution Mapping of the Tumor Microenvironment Using Integrated Single-Cell, Spatial and in Situ Analysis”. In: Nature Communications 14.1 (Dec. 2023), p. 8353. issn: 2041-1723. doi: 10.1038/s41467-023-43458-x . (Visited on 04/30/2024)

  2. [1]

    Mouse Brain Section (Coronal), Spatial Gene Expression Dataset Analyzed Using Space Ranger 1.0.0

    10x Genomics. Mouse Brain Section (Coronal), Spatial Gene Expression Dataset Analyzed Using Space Ranger 1.0.0 . Dec. 2019. (Visited on 11/20/2024)

  3. [2]

    GraphCompass: Spatial Metrics for Differential Analyses of Cell Organization across Conditions

    Mayar Ali et al. “GraphCompass: Spatial Metrics for Differential Analyses of Cell Organization across Conditions”. In: Bioinformatics 40.Supplement 1 (July 2024), pp. i548–i557. issn: 1367-

  4. [3]

    Matrix-Assisted Laser Desorption Ionization Imag- ing Mass Spectrometry: In Situ Molecular Mapping

    Peggi M. Angel and Richard M. Caprioli. “Matrix-Assisted Laser Desorption Ionization Imag- ing Mass Spectrometry: In Situ Molecular Mapping”. In: Biochemistry 52.22 (June 2013), pp. 3818–3828. issn: 0006-2960. doi: 10.1021/bi301519p. (Visited on 04/24/2025)

  5. [4]

    A Local Indicator of Multivariate Spatial Association: Extending Geary’s c

    Luc Anselin. “A Local Indicator of Multivariate Spatial Association: Extending Geary’s c”. In: Geographical Analysis 51.2 (2019), pp. 133–150. issn: 1538-4632. doi: 10.1111/gean.12164 . (Visited on 02/05/2024)

  6. [5]

    An Introduction to Spatial Data Science with GeoDa: Volume 1: Exploring Spatial Data

    Luc Anselin. An Introduction to Spatial Data Science with GeoDa: Volume 1: Exploring Spatial Data. CRC Press, 2024

  7. [6]

    Local Indicators of Spatial Association—LISA

    Luc Anselin. “Local Indicators of Spatial Association—LISA”. In: Geographical Analysis 27.2 (1995), pp. 93–115. issn: 1538-4632. doi: 10.1111/j.1538-4632.1995.tb00338.x . (Visited on 11/11/2023)

  8. [7]

    The Moran Scatterplot as an ESDA Tool to Assess Local Instability in Spatial Association

    Luc Anselin. “The Moran Scatterplot as an ESDA Tool to Assess Local Instability in Spatial Association”. In: Spatial Analytical Perspectives on GIS . Routledge, 2019, pp. 111–126

Show all 99 references
  1. [8]

    What Is Special About Spatial Data? Alternative Perspectives on Spatial Data Analysis

    Luc Anselin. “What Is Special About Spatial Data? Alternative Perspectives on Spatial Data Analysis”. In: Spring 1989 Symposium onSpatial Statistics, Past, Present and Future . 1989

  2. [9]

    Segment Anything for Microscopy

    Anwai Archit et al. Segment Anything for Microscopy . Aug. 2023. doi: 10.1101/2023.08. 21.554208. (Visited on 07/25/2024)

  3. [10]

    Spatially Resolved Transcriptomes— Next Generation Tools for Tissue Exploration

    Michaela Asp, Joseph Bergenstr ˚ ahle, and Joakim Lundeberg. “Spatially Resolved Transcriptomes— Next Generation Tools for Tissue Exploration”. In: BioEssays 42.10 (2020), p. 1900221. issn: 1521-1878. doi: 10.1002/bies.201900221. (Visited on 11/09/2022)

  4. [11]

    Spatial Point Patterns

    Adrian Baddeley, Ege Rubak, and Rolf Turner. Spatial Point Patterns. 1st ed. CRC Interdisci- plinary Statistics Series. CRC Press, Taylor & Francis Group, Dec. 2015.isbn: 9781482210200. 24

  5. [12]

    Spatstat: An R Package for Analyzing Spatial Point Pat- terns

    Adrian Baddeley and Rolf Turner. “Spatstat: An R Package for Analyzing Spatial Point Pat- terns”. In: Journal of Statistical Software 12 (Jan. 2005), pp. 1–42. issn: 1548-7660. doi: 10.18637/jss.v012.i06. (Visited on 02/20/2024)

  6. [13]

    Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing

    Yoav Benjamini and Yosef Hochberg. “Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing”. In: Journal of the Royal Statistical Society: Series B (Methodological) 57.1 (Jan. 1995), pp. 289–300. issn: 2517-6161. doi: 10.1111/j.2517- 6161.1...

  7. [14]

    Contribution to the Discussion on Dr Ripley’s Paper

    Julian Besag. “Contribution to the Discussion on Dr Ripley’s Paper”. In: JR Stat Soc B 39 (1977), pp. 193–195. doi: 10.1111/j.2517-6161.1977.tb01616.x

  8. [15]

    R Packages for Analyzing Spatial Data: A Comparative Case Study with Areal Data

    Roger Bivand. “R Packages for Analyzing Spatial Data: A Comparative Case Study with Areal Data”. In: Geographical Analysis 54.3 (2022), pp. 488–518. issn: 1538-4632. doi: 10.1111/ gean.12319. (Visited on 01/24/2024)

  9. [16]

    Bivand, Edzer Pebesma, and Virgilio G´ omez-Rubio

    Roger S. Bivand, Edzer Pebesma, and Virgilio G´ omez-Rubio. Applied Spatial Data Analysis with R, Second Edition . Springer, NY, 2013

  10. [17]

    CODEX Multiplexed Tissue Imaging with DNA-conjugated Antibodies

    Sarah Black et al. “CODEX Multiplexed Tissue Imaging with DNA-conjugated Antibodies”. In: Nature Protocols 16.8 (Aug. 2021), pp. 3802–3835. issn: 1750-2799. doi: 10.1038/s41596- 021-00556-8. (Visited on 04/30/2024)

  11. [18]

    The Dawn of Spatial Omics

    Dario Bressan, Giorgia Battistoni, and Gregory J. Hannon. “The Dawn of Spatial Omics”. In: Science 381.6657 (Aug. 2023), eabq4964. doi: 10.1126/science.abq4964 . (Visited on 03/06/2024)

  12. [19]

    Bull et al

    Joshua A. Bull et al. MuSpAn: A Toolbox for Multiscale Spatial Analysis . Dec. 2024. doi: 10.1101/2024.12.06.627195. (Visited on 05/14/2025)

  13. [21]

    scFeatures: Multi-View Representations of Single-Cell and Spatial Data for Dis- ease Outcome Prediction

    Yue Cao et al. “scFeatures: Multi-View Representations of Single-Cell and Spatial Data for Dis- ease Outcome Prediction”. In: Bioinformatics 38.20 (Oct. 2022). Ed. by Olga Vitek, pp. 4745–

  14. [22]

    Spatiotemporal Transcriptomic Atlas of Mouse Organogenesis Using DNA Nanoball-Patterned Arrays

    Ao Chen et al. “Spatiotemporal Transcriptomic Atlas of Mouse Organogenesis Using DNA Nanoball-Patterned Arrays”. In: Cell 185.10 (May 2022), 1777–1792.e21. issn: 0092-8674. doi: 10.1016/j.cell.2022.04.003. (Visited on 04/23/2025)

  15. [23]

    Spatially Resolved, Highly Multiplexed RNA Profiling in Single Cells

    Kok Hao Chen et al. “Spatially Resolved, Highly Multiplexed RNA Profiling in Single Cells”. In: Science 348.6233 (Apr. 2015), aaa6090. doi: 10 . 1126 / science . aaa6090. (Visited on 04/05/2024)

  16. [24]

    Noel A. C. Cressie. Statistics for Spatial Data . Rev. ed. Wiley Series in Probability and Math- ematical Statistics. New York, NY: Wiley, 1993. isbn: 978-0-471-00255-0. 25

  17. [25]

    Mark R. T. Dale and Marie-Jos´ ee Fortin. Spatial Analysis: A Guide for Ecologists . Second Edition. Cambridge ; New York: Cambridge University Press, 2014. isbn: 978-0-521-14350-9

  18. [26]

    Imcdatasets: Collection of Publicly Available Imaging Mass Cytometry (IMC) Datasets

    Nicolas Damond. Imcdatasets: Collection of Publicly Available Imaging Mass Cytometry (IMC) Datasets. Manual. 2024. doi: 10.18129/B9.bioc.imcdatasets

  19. [27]

    A Map of Human Type 1 Diabetes Progression by Imaging Mass Cytometry

    Nicolas Damond et al. “A Map of Human Type 1 Diabetes Progression by Imaging Mass Cytometry”. In: Cell Metabolism 29.3 (Mar. 2019), 755–768.e5. issn: 1550-4131. doi: 10 . 1016/j.cmet.2018.11.014. (Visited on 10/02/2024)

  20. [28]

    Systematic Assessment of Tissue Dissociation and Storage Biases in Single-Cell and Single-Nucleus RNA-seq Workflows

    Elena Denisenko et al. “Systematic Assessment of Tissue Dissociation and Storage Biases in Single-Cell and Single-Nucleus RNA-seq Workflows”. In: Genome Biology 21.1 (June 2020), p. 130. issn: 1474-760X. doi: 10.1186/s13059-020-02048-6 . (Visited on 04/23/2025)

  21. [29]

    Systematic Comparison of Single-Cell and Single-Nucleus RNA-sequencing Methods

    Jiarui Ding et al. “Systematic Comparison of Single-Cell and Single-Nucleus RNA-sequencing Methods”. In: Nature Biotechnology 38.6 (June 2020), pp. 737–746.issn: 1087-0156, 1546-1696. doi: 10.1038/s41587-020-0465-8 . (Visited on 04/23/2025)

  22. [30]

    Robinson

    Martin Emons, Samuel Gunz, and Mark D. Robinson. spatialFDA: A Tool for Spatial Multi- Sample Comparisons. Manual. 2025. doi: 10.18129/B9.bioc.spatialFDA

  23. [31]

    Spatial Analysis with SPIAT and spaSim to Characterize and Simulate Tissue Microenvironments

    Yuzhou Feng et al. “Spatial Analysis with SPIAT and spaSim to Characterize and Simulate Tissue Microenvironments”. In: Nature Communications 14.1 (May 2023), p. 2697. issn: 2041-

  24. [32]

    Macrophage and Neutrophil Heterogeneity at Single-Cell Spatial Resolution in Human Inflammatory Bowel Disease

    Alba Garrido-Trigo et al. “Macrophage and Neutrophil Heterogeneity at Single-Cell Spatial Resolution in Human Inflammatory Bowel Disease”. In: Nature Communications 14.1 (July 2023), p. 4506. issn: 2041-1723. doi: 10.1038/s41467-023-40156-6 . (Visited on 04/05/2024)

  25. [33]

    Bioconductor: Open Software Development for Computational Biology and Bioinformatics

    Robert C Gentleman et al. “Bioconductor: Open Software Development for Computational Biology and Bioinformatics”. In: Genome biology 5 (2004), pp. 1–16

  26. [34]

    Spatial Weights Matrices

    Arthur Getis. “Spatial Weights Matrices”. In: Geographical Analysis 41.4 (2009), pp. 404–410. issn: 1538-4632. doi: 10.1111/j.1538-4632.2009.00768.x. (Visited on 01/18/2024)

  27. [35]

    The Analysis of Spatial Association by Use of Distance Statis- tics

    Arthur Getis and J. K. Ord. “The Analysis of Spatial Association by Use of Distance Statis- tics”. In: Geographical Analysis 24.3 (1992), pp. 189–206. issn: 1538-4632. doi: 10.1111/j. 1538-4632.1992.tb00261.x. (Visited on 01/03/2024)

  28. [36]

    Highly Multiplexed Imaging of Tumor Tissues with Subcellular Resolu- tion by Mass Cytometry

    Charlotte Giesen et al. “Highly Multiplexed Imaging of Tumor Tissues with Subcellular Resolu- tion by Mass Cytometry”. In: Nature Methods 11.4 (Apr. 2014), pp. 417–422. issn: 1548-7105. doi: 10.1038/nmeth.2869. (Visited on 04/03/2024)

  29. [37]

    Analyzing Spatial Point Patterns in Digital Pathology: Immune Cells in High-Grade Serous Ovarian Carcinomas

    Jonatan A Gonz´ alez et al. “Analyzing Spatial Point Patterns in Digital Pathology: Immune Cells in High-Grade Serous Ovarian Carcinomas”. In:The American Statistician (2025), pp. 1– 26

  30. [38]

    Combining Incompatible Spatial Data

    Carol A Gotway and Linda J Young. “Combining Incompatible Spatial Data”. In: Journal of the American Statistical Association 97.458 (June 2002), pp. 632–648. issn: 0162-1459. doi: 10.1198/016214502760047140. (Visited on 11/27/2024). 26

  31. [39]

    Whole-Cell Segmentation of Tissue Images with Human-Level Per- formance Using Large-Scale Data Annotation and Deep Learning

    Noah F. Greenwald et al. “Whole-Cell Segmentation of Tissue Images with Human-Level Per- formance Using Large-Scale Data Annotation and Deep Learning”. In: Nature Biotechnology 40.4 (Apr. 2022), pp. 555–565. issn: 1546-1696. doi: 10.1038/s41587-021-01094-0 . (Visited on 07/25/2024)

  32. [40]

    Robinson

    Samuel Gunz and Mark D. Robinson. Sosta: A Package for the Analysis of Anatomical Tissue Structures in Spatial Omics Data . Manual. 2025. doi: 10.18129/B9.bioc.sosta

  33. [41]

    Multiplexed Protein Maps Link Subcellular Organization to Cellular States

    Gabriele Gut, Markus D. Herrmann, and Lucas Pelkmans. “Multiplexed Protein Maps Link Subcellular Organization to Cellular States”. In: Science 361.6401 (Aug. 2018), eaar7042. doi: 10.1126/science.aar7042. (Visited on 04/30/2024)

  34. [42]

    Dictionary Learning for Integrative, Multimodal and Scalable Single-Cell Analysis

    Yuhan Hao et al. “Dictionary Learning for Integrative, Multimodal and Scalable Single-Cell Analysis”. In: Nature Biotechnology (2023). doi: 10.1038/s41587-023-01767-y

  35. [43]

    High-Plex Imaging of RNA and Proteins at Subcellular Resolution in Fixed Tissue by Spatial Molecular Imaging

    Shanshan He et al. “High-Plex Imaging of RNA and Proteins at Subcellular Resolution in Fixed Tissue by Spatial Molecular Imaging”. In: Nature Biotechnology 40.12 (Dec. 2022), pp. 1794–

  36. [45]

    MIBI-TOF: A Multiplexed Imaging Platform Relates Cellular Phenotypes and Tissue Structure

    Leeat Keren et al. “MIBI-TOF: A Multiplexed Imaging Platform Relates Cellular Phenotypes and Tissue Structure”. In: Science Advances 5.10 (Oct. 2019), eaax5851. doi: 10 . 1126 / sciadv.aax5851. (Visited on 04/30/2024)

  37. [46]

    Spatiomic - Scalable Spatial Proteomics Analyses for Pathology

    Malte Kuehl et al. Spatiomic - Scalable Spatial Proteomics Analyses for Pathology . June 2025. (Visited on 06/13/2025)

  38. [47]

    Immunological Method for Mapping Genes on Drosophila Polytene Chromosomes

    Pennina R. Langer-Safer, Michael Levine, and David C. Ward. “Immunological Method for Mapping Genes on Drosophila Polytene Chromosomes.” In:Proceedings of the National Academy of Sciences of the United States of America 79.14 (July 1982), pp. 4381–4385. issn: 0027-8424. doi: 1...

  39. [48]

    Developing a Bivariate Spatial Association Measure: An Integration of Pearson’s r and Moran’s I

    Sang-Il Lee. “Developing a Bivariate Spatial Association Measure: An Integration of Pearson’s r and Moran’s I”. In: Journal of Geographical Systems 3.4 (Dec. 2001), pp. 369–385. issn: 1435-5930. doi: 10.1007/s101090100064. (Visited on 01/24/2024)

  40. [49]

    SpatialDM for Rapid Identification of Spatially Co-Expressed Ligand– Receptor and Revealing Cell–Cell Communication Patterns

    Zhuoxuan Li et al. “SpatialDM for Rapid Identification of Spatially Co-Expressed Ligand– Receptor and Revealing Cell–Cell Communication Patterns”. In: Nature Communications 14.1 (July 2023), p. 3995. issn: 2041-1723. doi: 10 . 1038 / s41467 - 023 - 39608 - w. (Visited on 04/03/2024)

  41. [50]

    Spatially Resolved Epigenomic Profiling of Single Cells in Complex Tissues

    Tian Lu, Cheen Euong Ang, and Xiaowei Zhuang. “Spatially Resolved Epigenomic Profiling of Single Cells in Complex Tissues”. In: Cell 185.23 (Nov. 2022), 4448–4464.e17. issn: 0092-8674. doi: 10.1016/j.cell.2022.09.035. (Visited on 04/04/2024). 27

  42. [52]

    SpatialData: An Open and Universal Data Framework for Spatial Omics

    Luca Marconato et al. “SpatialData: An Open and Universal Data Framework for Spatial Omics”. In: Nature Methods 22.1 (2025), pp. 58–62

  43. [53]

    Scater: Pre-Processing, Quality Control, Normalisation and Visual- isation of Single-Cell RNA-seq Data in R

    Davis J. McCarthy et al. “Scater: Pre-Processing, Quality Control, Normalisation and Visual- isation of Single-Cell RNA-seq Data in R”. In: Bioinformatics (Oxford, England) 33.8 (2017), pp. 1179–1186. doi: 10.1093/bioinformatics/btw777

  44. [54]

    Characterizing Spatial Gene Expression Heterogeneity in Spatially Resolved Single-Cell Transcriptomics Data with Nonuniform Cellular Densities

    Brendan F. Miller et al. “Characterizing Spatial Gene Expression Heterogeneity in Spatially Resolved Single-Cell Transcriptomics Data with Nonuniform Cellular Densities”. In: Genome Research (May 2021), gr.271288.120. issn: 1088-9051, 1549-5469. doi: 10.1101/gr.271288

  45. [55]

    The Emerging Landscape of Spatial Profiling Technologies

    Jeffrey R. Moffitt, Emma Lundberg, and Holger Heyn. “The Emerging Landscape of Spatial Profiling Technologies”. In: Nature Reviews Genetics (July 2022). issn: 1471-0056, 1471-0064. doi: 10.1038/s41576-022-00515-3 . (Visited on 11/07/2022)

  46. [56]

    Fast and Robust Feature- Based Stitching Algorithm for Microscopic Images

    Fatemeh Sadat Mohammadi, Hasti Shabani, and Mojtaba Zarei. “Fast and Robust Feature- Based Stitching Algorithm for Microscopic Images”. In: Scientific Reports 14.1 (June 2024), p. 13304. issn: 2045-2322. doi: 10.1038/s41598-024-61970-y . (Visited on 07/25/2024)

  47. [57]

    Notes on Continuous Stochastic Phenomena

    P. A. P. Moran. “Notes on Continuous Stochastic Phenomena”. In: Biometrika 37.1/2 (1950), pp. 17–23. issn: 0006-3444. doi: 10.2307/2332142. JSTOR: 2332142. (Visited on 01/18/2024)

  48. [58]

    Museum of Spatial Transcriptomics

    Lambda Moses and Lior Pachter. “Museum of Spatial Transcriptomics”. In: Nature Methods 19.5 (May 2022), pp. 534–546. issn: 1548-7105. doi: 10.1038/s41592-022-01409-2 . (Visited on 11/08/2022)

  49. [59]

    Voyager: Exploratory Single-Cell Genomics Data Analysis with Geospatial Statistics

    Lambda Moses et al. Voyager: Exploratory Single-Cell Genomics Data Analysis with Geospatial Statistics. Aug. 2023. doi: 10.1101/2023.07.20.549945. (Visited on 08/29/2023)

  50. [60]

    Oliveira et al

    Michelli F. Oliveira et al. Characterization of Immune Cell Populations in the Tumor Mi- croenvironment of Colorectal Cancer Using High Definition Spatial Profiling . June 2024. doi: 10.1101/2024.06.04.597233. (Visited on 10/09/2024)

  51. [61]

    The Modifiable Areal Unit Problem

    Stan Openshaw. The Modifiable Areal Unit Problem . Concepts and Techniques in Modern Geography 38. Norwich: Geo, 1984. isbn: 978-0-86094-134-7

  52. [62]

    Spatial Components of Molecular Tissue Biology

    Giovanni Palla et al. “Spatial Components of Molecular Tissue Biology”. In: Nature Biotech- nology 40.3 (Mar. 2022), pp. 308–318. issn: 1546-1696. doi: 10.1038/s41587-021-01182-1 . (Visited on 11/09/2022)

  53. [63]

    Squidpy: A Scalable Framework for Spatial Omics Analysis

    Giovanni Palla et al. “Squidpy: A Scalable Framework for Spatial Omics Analysis”. In: Nature Methods 19.2 (Feb. 2022), pp. 171–178. issn: 1548-7105. doi: 10.1038/s41592-021-01358-2 . (Visited on 04/08/2024). 28

  54. [64]

    Spatial Omics Technologies at Multimodal and Single Cell/Subcellular Level

    Jiwoon Park et al. “Spatial Omics Technologies at Multimodal and Single Cell/Subcellular Level”. In: Genome Biology 23.1 (Dec. 2022), p. 256. issn: 1474-760X. doi: 10.1186/s13059- 022-02824-6. (Visited on 07/25/2024)

  55. [65]

    Spatial Analysis for Highly Multiplexed Imaging Data to Identify Tissue Microenvironments

    Ellis Patrick et al. “Spatial Analysis for Highly Multiplexed Imaging Data to Identify Tissue Microenvironments”. In: Cytometry Part A 103.7 (May 2023), pp. 593–599. issn: 1552-4930. doi: 10.1002/cyto.a.24729. (Visited on 04/18/2023)

  56. [66]

    Spatial Data Science: With Applications in R

    Edzer Pebesma and Roger Bivand. Spatial Data Science: With Applications in R . 1st ed. New York: Chapman and Hall/CRC, May 2023. isbn: 978-0-429-45901-6. doi: 10 . 1201 / 9780429459016. (Visited on 01/18/2024)

  57. [67]

    Cell Segmentation in Imaging-Based Spatial Transcriptomics

    Viktor Petukhov et al. “Cell Segmentation in Imaging-Based Spatial Transcriptomics”. In: Nature Biotechnology 40.3 (Mar. 2022), pp. 345–354. issn: 1546-1696. doi: 10.1038/s41587- 021-01044-w. (Visited on 01/18/2024)

  58. [68]

    Bin2cell Reconstructs Cells from High Resolution Visium HD Data

    Krzysztof Pola´ nski et al. “Bin2cell Reconstructs Cells from High Resolution Visium HD Data”. In: Bioinformatics 40.9 (Sept. 2024), btae546.issn: 1367-4811. doi: 10.1093/bioinformatics/ btae546. (Visited on 11/28/2024)

  59. [69]

    Comparison of Spatial Transcriptomics Technologies Using Tumor Cryosections

    Anne Rademacher et al. Comparison of Spatial Transcriptomics Technologies Using Tumor Cryosections. Apr. 2024. doi: 10.1101/2024.04.03.586404. (Visited on 10/09/2024)

  60. [71]

    PySAL: A Python Library of Spatial Analytical Methods

    Sergio J. Rey and Luc Anselin. “PySAL: A Python Library of Spatial Analytical Methods”. In: Handbook of Applied Spatial Analysis . Springer, Berlin, Heidelberg, 2010, pp. 175–193. isbn: 978-3-642-03647-7. doi: 10.1007/978-3-642-03647-7_11 . (Visited on 12/11/2024)

  61. [72]

    SpatialExperiment: Infrastructure for Spatially-Resolved Transcriptomics Data in R Using Bioconductor

    Dario Righelli et al. “SpatialExperiment: Infrastructure for Spatially-Resolved Transcriptomics Data in R Using Bioconductor”. In: Bioinformatics 38.11 (May 2022), pp. 3128–3131. issn: 1367-4803. doi: 10.1093/bioinformatics/btac299. (Visited on 05/23/2024)

  62. [73]

    The Second-Order Analysis of Stationary Point Processes

    Brian D. Ripley. “The Second-Order Analysis of Stationary Point Processes”. In: Journal of Applied Probability 13.2 (1976), pp. 255–266. issn: 0021-9002. doi: 10.2307/3212829. JSTOR: 3212829. (Visited on 03/11/2024)

  63. [74]

    Finding the Edge of a Poisson Forest

    Brian D. Ripley and Jean-Paul Rasson. “Finding the Edge of a Poisson Forest”. In: Journal of Applied Probability 14.3 (Sept. 1977), pp. 483–491. issn: 0021-9002, 1475-6072. doi: 10.2307/ 3213451. (Visited on 03/11/2024)

  64. [75]

    Slide-Seq: A Scalable Technology for Measuring Genome-Wide Expression at High Spatial Resolution

    Samuel G. Rodriques et al. “Slide-Seq: A Scalable Technology for Measuring Genome-Wide Expression at High Spatial Resolution”. In: Science 363.6434 (Mar. 2019), pp. 1463–1467. doi: 10.1126/science.aaw1219. (Visited on 04/30/2024)

  65. [76]

    SpaNorm: Spatially-Aware Normalisation for Spatial Transcriptomics Data

    Agus Salim et al. SpaNorm: Spatially-Aware Normalisation for Spatial Transcriptomics Data . June 2024. doi: 10.1101/2024.05.31.596908. (Visited on 05/30/2025). 29

  66. [77]

    Spatial Atlas of the Mouse Central Nervous System at Molecular Resolu- tion

    Hailing Shi et al. “Spatial Atlas of the Mouse Central Nervous System at Molecular Resolu- tion”. In: Nature 622.7983 (Oct. 2023), pp. 552–561. issn: 1476-4687. doi: 10.1038/s41586- 023-06569-5. (Visited on 08/15/2024)

  67. [78]

    BANKSY Unifies Cell Typing and Tissue Domain Segmentation for Scalable Spatial Omics Data Analysis

    Vipul Singhal et al. “BANKSY Unifies Cell Typing and Tissue Domain Segmentation for Scalable Spatial Omics Data Analysis”. In: Nature Genetics (2024). doi: 10.1038/s41588- 024-01664-3

  68. [79]

    Visualization and Analysis of Gene Expression in Tissue Sections by Spatial Transcriptomics

    Patrik L. St ˚ ahl et al. “Visualization and Analysis of Gene Expression in Tissue Sections by Spatial Transcriptomics”. In:Science 353.6294 (July 2016), pp. 78–82.doi: 10.1126/science. aaf2403. (Visited on 04/08/2024)

  69. [80]

    Highly Sensitive Spatial Transcriptomics at Near-Cellular Resolution with Slide-seqV2

    Robert R. Stickels et al. “Highly Sensitive Spatial Transcriptomics at Near-Cellular Resolution with Slide-seqV2”. In: Nature Biotechnology 39.3 (Mar. 2021), pp. 313–319. issn: 1546-1696. doi: 10.1038/s41587-020-0739-1 . (Visited on 04/30/2024)

  70. [81]

    Cellpose: A Generalist Algorithm for Cellular Segmentation

    Carsen Stringer et al. “Cellpose: A Generalist Algorithm for Cellular Segmentation”. In: Nature Methods 18.1 (Jan. 2021), pp. 100–106. issn: 1548-7105. doi: 10.1038/s41592-020-01018-x . (Visited on 07/25/2024)

  71. [82]

    A Computer Movie Simulating Urban Growth in the Detroit Region

    W. R. Tobler. “A Computer Movie Simulating Urban Growth in the Detroit Region”. In: Economic Geography 46 (1970), pp. 234–240. issn: 0013-0095. doi: 10.2307/143141. JSTOR: 143141. (Visited on 03/06/2024)

  72. [83]

    Quantitative Multiplex Immunohistochemistry Reveals Myeloid- Inflamed Tumor-Immune Complexity Associated with Poor Prognosis

    Takahiro Tsujikawa et al. “Quantitative Multiplex Immunohistochemistry Reveals Myeloid- Inflamed Tumor-Immune Complexity Associated with Poor Prognosis”. In: Cell Reports 19.1 (Apr. 2017), pp. 203–217. issn: 2211-1247. doi: 10.1016/j.celrep.2017.03.037 . (Visited on 07/25/2024)

  73. [84]

    Effects of Changing Spatial Scale on the Analysis of Landscape Pattern

    Monica G. Turner et al. “Effects of Changing Spatial Scale on the Analysis of Landscape Pattern”. In: Landscape Ecology 3.3 (Dec. 1989), pp. 153–162. issn: 1572-9761. doi: 10.1007/ BF00131534. (Visited on 10/02/2024)

  74. [85]

    Methods and Applications for Single-Cell and Spatial Multi-Omics

    Katy Vandereyken et al. “Methods and Applications for Single-Cell and Spatial Multi-Omics”. In: Nature Reviews Genetics 24.8 (Aug. 2023), pp. 494–515. issn: 1471-0064. doi: 10.1038/ s41576-023-00580-2 . (Visited on 11/20/2024)

  75. [86]

    Anndata: Access and Store Annotated Data Matrices

    Isaac Virshup et al. “Anndata: Access and Store Annotated Data Matrices”. In: Journal of Open Source Software 9.101 (2024), p. 4371. doi: 10.21105/joss.04371

  76. [87]

    Cell Segmentation for Image Cytometry: Advances, Insufficiencies, and Challenges

    Zhenzhou Wang. “Cell Segmentation for Image Cytometry: Advances, Insufficiencies, and Challenges”. In: Cytometry Part A 95.7 (2019), pp. 708–711. issn: 1552-4930. doi: 10.1002/ cyto.a.23686. (Visited on 01/18/2024)

  77. [88]

    Weber and Yixing E

    Lukas M. Weber and Yixing E. Dong. STexampleData: Collection of Spatial Transcriptomics Datasets in SpatialExperiment Bioconductor Format. Manual. 2024. doi: 10.18129/B9.bioc. STexampleData. 30

  78. [89]

    Ggplot2: Elegant Graphics for Data Analysis

    Hadley Wickham. Ggplot2: Elegant Graphics for Data Analysis . Springer-Verlag New York,

  79. [90]

    An End-to-End Workflow for Multiplexed Image Processing and Analysis

    Jonas Windhager et al. “An End-to-End Workflow for Multiplexed Image Processing and Analysis”. In: Nature Protocols 18.11 (Nov. 2023), pp. 3565–3613. issn: 1750-2799. doi: 10. 1038/s41596-023-00881-0 . (Visited on 10/02/2024)

  80. [91]

    SpottedPy Quantifies Relationships between Spatial Tran- scriptomic Hotspots and Uncovers Environmental Cues of Epithelial-Mesenchymal Plastic- ity in Breast Cancer

    Eloise Withnell and Maria Secrier. “SpottedPy Quantifies Relationships between Spatial Tran- scriptomic Hotspots and Uncovers Environmental Cues of Epithelial-Mesenchymal Plastic- ity in Breast Cancer”. In: Genome Biology 25.1 (Nov. 2024), p. 289. issn: 1474-760X. doi: 10.1186...

  81. [92]

    SCANPY: Large-Scale Single-Cell Gene Expression Data Analysis

    F Alexander Wolf, Philipp Angerer, and Fabian J Theis. “SCANPY: Large-Scale Single-Cell Gene Expression Data Analysis”. In: Genome biology 19 (2018), pp. 1–5

  82. [93]

    Mxfda: A Comprehensive Toolkit for Functional Data Analysis of Single- Cell Spatial Data

    Julia Wrobel et al. “Mxfda: A Comprehensive Toolkit for Functional Data Analysis of Single- Cell Spatial Data”. In: Bioinformatics Advances 4.1 (Jan. 2024), vbae155. issn: 2635-0041. doi: 10.1093/bioadv/vbae155. (Visited on 11/27/2024)

  83. [94]

    Mass Spectrometry Imaging: The Rise of Spatially Resolved Single-Cell Omics

    Hua Zhang, Daniel G. Delafield, and Lingjun Li. “Mass Spectrometry Imaging: The Rise of Spatially Resolved Single-Cell Omics”. In: Nature Methods 20.3 (Mar. 2023), pp. 327–330. issn: 1548-7105. doi: 10.1038/s41592-023-01774-6 . (Visited on 04/23/2025)

  84. [95]

    Mapping the Transcriptome: Realizing the Full Potential of Spatial Data Analysis

    Eleftherios Zormpas et al. “Mapping the Transcriptome: Realizing the Full Potential of Spatial Data Analysis”. In: Cell 186.26 (Dec. 2023), pp. 5677–5689. issn: 00928674. doi: 10.1016/j. cell.2023.11.003. (Visited on 12/11/2023)

  85. [96]

    Zuur, Elena N

    Alain F. Zuur, Elena N. Ieno, and Graham M. Smith. Analysing Ecological Data. Statistics for Biology and Health. New York: Springer, 2007. isbn: 978-0-387-45967-7. 31 Supplementary data Table S1: Non-exhaustive table of technologies and the data modality they can represent. Fo...

  86. [120]

    (Visited on 04/03/2024)

  87. [1723]

    (Visited on 04/08/2024)

    doi: 10.1038/s41467-023-37822-0 . (Visited on 04/08/2024)

  88. [1806]

    doi: 10.1038/s41587-022-01483-z

    issn: 1546-1696. doi: 10.1038/s41587-022-01483-z . (Visited on 09/27/2023)

  89. [2016]

    isbn: 978-3-319-24277-4

  90. [4753]

    issn: 1367-4803, 1367-4811. doi: 10 . 1093 / bioinformatics / btac590. (Visited on 12/01/2022)

  91. [4811]

    (Visited on 10/02/2024)

    doi: 10.1093/bioinformatics/btae242. (Visited on 10/02/2024)

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.