Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

PySpatial: A High-Speed Whole Slide Image Pathomics Toolkit

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read PySpatial claims that skipping patch segmentation and processing annotated regions directly with R-tree indexing and matrix batch computation speeds up whole slide pathomics feature extraction by roughly 10-fold on small objects and…

desk verdict A sensible engineering contribution that likely delivers real speedups, but the benchmark reporting is currently too thin to trust the headline numbers. read the letter →

arxiv 2501.06151 v1 pith:GLKWHEIM submitted 2025-01-10 eess.IV cs.CV

Yuechen Yang , Yu Wang , Tianyuan Yao , Ruining Deng , Mengmeng Yin , Shilin Zhao , Haichun Yang , Yuankai Huo This is my paper · ORCID
classification eess.IVcs.CV
keywords PySpatialwholeslideimagepathomicsfeatureextractionR-treeindexmatrixbatchcomputationdigitalpathology
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

PySpatial is a toolkit that tries to prove that whole slide image feature extraction can skip the conventional patch-split-and-merge step entirely. It processes only the annotated computational regions, indexes them with an R-tree, and computes features in matrix batches, reporting a nearly 10-fold speedup on small dense objects (PEC) and about a 2-fold speedup on large sparse objects (KPMP) relative to a CellProfiler-based pipeline. The paper also claims that the features it produces are consistent with CellProfiler's, based on matching frequency distributions. If these claims hold, large-scale digital pathology studies could cut one of the most expensive computational steps in tissue analysis without losing spatial context.

What carries the argument

The load-bearing mechanism is the pairing of an R-tree spatial index with matrix-based batch computation. The R-tree maps every annotated object to its original coordinates on the whole slide, so the pipeline can bypass patch tiling and post-hoc coordinate reassembly. Matrix encoding lets geometric and intensity computations run simultaneously across all objects, and a customizable matrix-size parameter plus a direct per-object API prevents memory overflows when objects are very large. This combination is what the paper says produces the speedup while keeping features spatially linked to the slide.

What would settle it

Run both pipelines on the same whole slide images on identical hardware, fixing and reporting CellProfiler's patch parameters and run count; if PySpatial's total processing time is not consistently lower across object types, the central speedup claim fails. A second check is to compute per-object feature differences for exactly matching objects rather than comparing aggregate frequency distributions.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that WSI-level pathomics analysis does not require intermediate patch segmentation and coordinate remapping. PySpatial directly operates on computational regions defined by annotations, builds an R-tree spatial index to preserve each object's location, converts the regions into matrix form, and extracts a 247-dimensional feature vector per object across four categories (Size & Shape, Texture, Intensity, Intensity Distribution). The reported outcome is a consistent decrease in processing time across both datasets, with the largest gain on small, abundant objects, and feature frequency distributions that track CellProfiler's closely. The paper attributes the speed advantage to eliminating redundant processing steps and to batch matrix computation, with a fallback API for large objects that exceed memory.

Load-bearing premise

The speed comparison assumes the CellProfiler baseline is configured fairly, but the paper does not report patch size, overlap, hardware, or number of runs, so a poorly configured baseline would inflate PySpatial's reported speedups.

Editorial extensions

If this is right

  • Processing time scales with the total area of annotated regions rather than with the full slide area, so slides with sparse annotations become much cheaper to process.
  • The same R-tree plus matrix-batch design can be reused for any region-based WSI analysis, including deep-learning segmentation outputs, without reworking the feature extraction core.
  • For large-object datasets, memory limits force the fallback to direct per-object computation, reducing but not eliminating the speed gain.
  • Because features are indexed to original coordinates, downstream spatial statistics and visualization can consume PySpatial output directly, preserving the slide's spatial layout.

Reading between the lines

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

  • The paper leaves implicit that this design could shift the compute bottleneck in large cohort studies from feature extraction to the segmentation step that produces the annotations.
  • A natural next experiment would be to benchmark PySpatial against another patch-free pipeline to separate the benefit of skipping patches from the benefit of the specific matrix kernel.
  • The consistency check compares distributions, not per-object values; a per-object correlation on identical objects would strengthen the accuracy claim beyond what is reported.
  • For very large slides, the matrix-memory tradeoff suggests a hybrid strategy—batch small regions, compute large regions individually, then merge—that the paper mentions but does not benchmark.
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

3 major / 4 minor

Summary. The manuscript introduces PySpatial, a pathomics toolkit that extracts object-level features directly from whole-slide images (WSIs) using R-tree spatial indexing and matrix-based batch computation, bypassing the patch-level processing used in CellProfiler workflows. The paper describes a 247-dimensional feature vector spanning Size & Shape, Texture, Intensity, and Intensity Distribution, and reports experiments on two datasets: PEC (small, dense objects) and KPMP (large, sparse glomeruli and arteries). It claims approximately 10x and 2x speedups over CellProfiler on these datasets and supports the accuracy claim with side-by-side frequency distributions for four selected features.

Significance. The proposed architecture is plausible and addresses a genuine bottleneck: repeated patch I/O and sequential per-object loops are known to dominate pathomics pipelines. The two real datasets are complementary, and the idea of focusing computation on sparse annotated regions is practically valuable. However, the central quantitative claims are not yet supported by a reproducible benchmark: the CellProfiler baseline is underspecified, no run-level variance is reported, and the feature-consistency validation is qualitative and limited to four features. If the missing experimental details and quantitative agreement metrics are supplied, the toolkit could be a useful contribution to digital pathology; as written, the evidence is insufficient to verify the reported speedups or the claim of maintained accuracy.

major comments (3)
  1. [Consistency Between PySpatial and CellProfiler] The central speedup claims in Figure 6 are not reproducible because the CellProfiler baseline is underspecified. The Experimental Setup section states only that the CellProfiler workflow divides the WSI into patches and records coordinates; it omits the patch size, patch overlap, CellProfiler version and mode, hardware specifications (CPU model, RAM), number of WSI slides per category, number of repeated runs, and whether the timing includes patching, coordinate mapping, feature merging, or only feature extraction. It is also unclear whether both pipelines consumed the same precomputed object masks or whether CellProfiler performed segmentation inside the comparison. Without these details, the reported 10-fold (PEC) and 2-fold (KPMP) ratios cannot be distinguished from an arbitrarily slow baseline; if small patches with high overlap were used, or if segmentation rather than feature extraction was being timed, the comparison would overstate PySpatial's advantage. The authors should report the full baseline configuration, use identical precomputed masks in both pipelines, provide per-slide timings with standard deviations, and preferably include CellProfiler's whole-slide or parallel mode as an additional baseline.
  2. [Consistency Between PySpatial and CellProfiler] The accuracy/consistency validation is purely qualitative and covers only four of the 247 features. Figure 5 shows side-by-side frequency distributions for Max Feret Diameter, Eccentricity, Hu Moment, and Mean Intensity, but no numerical agreement metric is provided, and the Conclusion incorrectly refers to 'area feature frequency distributions' even though area is not among the displayed features. The paper should report quantitative agreement on the same objects (for example, correlation, mean relative error, or Bland-Altman limits) and should extend the check beyond four features in order to support the claim that the full 247-dimensional feature set maintains 'high accuracy' relative to CellProfiler.
  3. [Memory Constraints and Matrix Optimization] The timing comparison does not state which of PySpatial's two computation modes (matrix-based batch or per-object fallback) was used for each dataset and category. The Memory Constraints section explains that large KPMP objects can force a fallback to per-object computation, so the 2-fold speedup may represent a mixture of modes rather than a single algorithmic comparison. The authors should report the mode used for each category, the memory limits encountered, and the impact of the 'customizable matrix size' parameter on the measured times.
minor comments (4)
  1. [Running Performance Comparison] Figure 6 appears to show raw processing-time bars without visible numeric y-axis labels or error bars; the exact times and the number of slides per category should be stated so that the speedup ratios can be verified.
  2. [Methods] The 247 features are not enumerated and no formulas or parameters are given for the Size & Shape, Texture, Intensity, and Intensity Distribution modules; for a toolkit paper, the exact feature definitions and a link to the software repository are needed for reproducibility.
  3. [References] Reference [13] contains the placeholder text 'Accessed Month Day, Year' and should be completed with the actual access date.
  4. [General] There are several language issues that should be corrected, such as 'utilizes matrix-level operations' and 'PySpatial also provides API, which allows users to skip matrixization and directly compute object-level features'.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: PySpatial is an empirical performance benchmark, not a derivation whose outputs are equivalent to its inputs.

full rationale

PySpatial does not present an analytic derivation chain or a first-principles prediction; its central claims are empirical measurements of processing time and feature-consistency against CellProfiler. The speed comparison is an experimental benchmark, not a quantity obtained by fitting parameters to the same data. The feature-consistency check compares PySpatial's outputs with CellProfiler's as a reference implementation; although the paper could be strengthened by an external ground truth, agreement with a widely used tool is a standard implementation-fidelity check and does not make the timing result equivalent to an input. The only self-citation, reference [4] from the same Vanderbilt group, is used illustratively to describe a traditional patching workflow and is not load-bearing for any of the paper's conclusions. No equations are reused to define the result, no fitted parameters are renamed as predictions, and no uniqueness theorem is imported from prior work. The absent details about the CellProfiler baseline (patch size, overlap, hardware, runs) are reproducibility concerns, not circularity. Therefore no circularity pattern is present.

Assumptions & free parameters 1 free parameters · 3 assumptions · 0 invented entities

The paper's central speedup claim rests on an unspecified baseline configuration (patch size, hardware, repetitions) and on the assumption that CellProfiler features are the reference for accuracy. No new entities are introduced.

free parameters (1)
  • CellProfiler patch size
    The baseline pipeline's patch size is not reported; the speedup ratio depends on this hand-chosen value.
assumptions (3)
  • domain assumption CellProfiler features are the ground truth for pathomics feature extraction.
    PySpatial's accuracy is validated by comparing feature distributions to CellProfiler (Figure 5), treating CellProfiler as the reference without an external ground truth.
  • domain assumption The annotated computational regions (from manual or deep learning masks) correctly capture all objects of interest.
    The pipeline processes only annotated regions; if annotations miss objects, features are missing. The paper assumes annotations are accurate.
  • domain assumption The hardware and software environment is consistent for both pipelines.
    The experimental setup does not describe hardware, memory limits, or whether both pipelines run on the same machine; runtime comparisons require a controlled environment.

how reviews work

0 comments
Cite this review

Pith. "Pith review of PySpatial: A High-Speed Whole Slide Image Pathomics Toolkit." pith.science (2026). https://pith.science/paper/GLKWHEIM

@misc{pith2026250106151,
  author       = {Pith},
  title        = {Pith review of: PySpatial: A High-Speed Whole Slide Image Pathomics Toolkit},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GLKWHEIM}},
  note         = {Machine review of arXiv:2501.06151}
}
read the original abstract

Whole Slide Image (WSI) analysis plays a crucial role in modern digital pathology, enabling large-scale feature extraction from tissue samples. However, traditional feature extraction pipelines based on tools like CellProfiler often involve lengthy workflows, requiring WSI segmentation into patches, feature extraction at the patch level, and subsequent mapping back to the original WSI. To address these challenges, we present PySpatial, a high-speed pathomics toolkit specifically designed for WSI-level analysis. PySpatial streamlines the conventional pipeline by directly operating on computational regions of interest, reducing redundant processing steps. Utilizing rtree-based spatial indexing and matrix-based computation, PySpatial efficiently maps and processes computational regions, significantly accelerating feature extraction while maintaining high accuracy. Our experiments on two datasets-Perivascular Epithelioid Cell (PEC) and data from the Kidney Precision Medicine Project (KPMP)-demonstrate substantial performance improvements. For smaller and sparse objects in PEC datasets, PySpatial achieves nearly a 10-fold speedup compared to standard CellProfiler pipelines. For larger objects, such as glomeruli and arteries in KPMP datasets, PySpatial achieves a 2-fold speedup. These results highlight PySpatial's potential to handle large-scale WSI analysis with enhanced efficiency and accuracy, paving the way for broader applications in digital pathology.

Figures

Figures reproduced from arXiv: 2501.06151 by the authors.

Figure 1
Figure 1. Comparison between PySpatial and CellProfiler pipelines for WSI feature extraction. The top panel illustrates the PySpatial [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. PySpatial’s computational workflow for WSI feature extraction. The figure illustrates the computational workflow of PySpatial [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Overview of KPMP Dataset Object Categories.The [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Overview of Feature Categories and Subtypes Extracted [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]
Figure 5
Figure 5. Figure 5: Comparison of Feature Frequency Distributions between PySpatial and CellProfiler. The figure presents frequency distributions [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 6
Figure 6. Figure 6: Comparison of processing time between CellProfiler [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. DyMorph-B2I: Dynamic and Morphology-Guided Binary-to-Instance Segmentation for Renal Pathology

    cs.CV 2025-08 reject novelty 4.0 of 10

    DyMorph-B2I combines watershed, skeletonization, and morphology with per-image tuning to separate renal structures from binary masks, but its evaluation is circular because tuning uses the ground-truth object counts.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages · cited by 1 Pith paper

  1. [4]

    J. Chen, Y . Wang, R. Deng, Q. Liu, C. Cui, T. Yao, Y . Liu, J. Zhong, A. B. Fogo, H. Yang, et al., ”Spatial pathomics toolkit for quantitative analysis of podocyte nuclei with histology and spatial transcriptomics data in renal pathology,” inMedical Imaging 2024: Digital and Com- putational Pathology, vol. 12933, pp. 252–260, 2024, SPIE

    work page 2024

  2. [1]

    Gupta, T

    R. Gupta, T. Kurc, A. Sharma, J. S. Almeida, and J. Saltz, ”The emergence of pathomics,” Current Pathobiology Reports, vol. 7, pp. 73–84, 2019, Springer

    work page 2019

  3. [2]

    D. R. Stirling, M. J. Swain-Bowden, A. M. Lucas, A. E. Carpenter, B. A. Cimini, and A. Goodman, ”CellProfiler 4: improvements in speed, utility and usability,” BMC Bioinformatics, vol. 22, pp. 1–11, 2021, Springer

    work page 2021

  4. [3]

    D. L. H ¨olscher, N. Bouteldja, M. Joodaki, M. L. Russo, Y .-C. Lan, A. V . Sadr, M. Cheng, V . Tesar, S. V . Stillfried, B. M. Klinkhammer, et al., ”Next-Generation Morphometry for pathomics-data mining in histopathology,”Nature Communications, vol. 14, no. 1, p. 470, 2023, Nature Publishing Group UK London

    work page 2023

  5. [5]

    Gilley, K

    P. Gilley, K. Zhang, N. Abdoli, Y . Sadri, L. Adhikari, K.-M. Fung, and Y . Qiu, ”Utilizing a pathomics biomarker to predict the effectiveness of bevacizumab in ovarian cancer treatment,”Bioengineering, vol. 11, no. 7, p. 678, 2024

    work page 2024

  6. [6]

    M. P. Humphries, P. Maxwell, and M. Salto-Tellez, ”QuPath: The global impact of an open source digital pathology system,” Compu- tational and Structural Biotechnology Journal, vol. 19, pp. 852–859, 2021, Elsevier

    work page 2021

  7. [7]

    Bankhead, M

    P. Bankhead, M. B. Loughrey, J. A. Fern ´andez, Y . Dombrowski, D. G. McArt, P. D. Dunne, S. McQuaid, R. T. Gray, L. J. Murray, H. G. Coleman, et al., ”QuPath: Open source software for digital pathol- ogy image analysis,” Scientific Reports, vol. 7, no. 1, pp. 1–7, 2017, Nature Publishing Group

    work page 2017

  8. [8]

    Kothari, J

    S. Kothari, J. H. Phan, T. H. Stokes, and M. D. Wang, ”Pathology imaging informatics for quantitative analysis of whole-slide images,” Journal of the American Medical Informatics Association , vol. 20, no. 6, pp. 1099–1108, 2013, BMJ Publishing Group

    work page 2013

Show all 13 references
  1. [9]

    Chakroun, N

    I. Chakroun, N. Michiels, and R. Wuyts, ”GPU-accelerated CellPro- filer,” in 2018 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), pp. 321–326, 2018, IEEE

    2018

  2. [10]

    Rslan, H

    E. Rslan, H. A. Hameed, and E. Ezzat, ”Spatial R-tree index based on grid division for query processing,” Int. J. Database Manag. Syst. (IJDMS), vol. 9, no. 6, pp. 25–36, 2017

    2017

  3. [11]

    Zhang, J.-P

    Z.-B. Zhang, J.-P. Zhang, J. Yang, and Y . Yang, ”A new approach to creating spatial index with R-tree,” in 2007 International Conference on Machine Learning and Cybernetics, vol. 5, pp. 2645–2648, 2007, IEEE

    2007

  4. [12]

    A. B. Fogo, ”Learning from deep learning and pathomics,” Kidney International, vol. 104, no. 6, pp. 1050–1053, 2023, Elsevier

    2023

  5. [13]

    Kidney precision medicine project data

    “Kidney precision medicine project data.” https://www.kpmp. org. Accessed Month Day, Year. Funded by the National Institute of Diabetes and Digestive and Kidney Diseases (Grant numbers: U01DK133081, U01DK133091, U01DK133092, U01DK133093, U01DK133095, U01DK133097, U01DK114866, ...

Pith tools

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