OSCR

Leveraging experimental vasculature data for high-resolution brain tumor simulations.

Code ↔ Paper

11 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 11 matches · 2 of them tie a paragraph to a whole file without given lines: weak matches, whose lines are not tinted
  1. [1] § Materials and methods › Model parameters › Cell motility magnitude and persistence ↔ src/lib/action/cellorientation.cpp, lines 31 to 93 · score 0.83 · persistent random walk, Wiener process, random movement, random vector, magnitude, motility · Report a mismatch
  2. [2] § Results and discussion › Vascular data processing for large-scale simulations › Microscopy stack grouping ↔ data/local_density_analysis/extract_local_vessel_properties.py, lines 63 to 126 · score 0.79 · removing single voxels, binary mask, vessel voxels, network length, convolution, erosion · Report a mismatch
  3. [3] § Materials and methods › Model parameters › Cell motility magnitude and persistence ↔ src/lib/action/cellorientation.h, lines 22 to 73 · score 0.77 · persistent random walk, Wiener process, random movement, magnitude, motility, cell · Report a mismatch
  4. [4] § Results and discussion › Influence of vessel network topology on tumor growth › Tumor growth behavior ↔ data/local_density_analysis/extract_local_vessel_properties.py, lines 63 to 126 · score 0.68 · vessel properties, tumor center, local density, network length, radius, simulations · Report a mismatch
  5. [5] § Materials and methods › Model description ↔ src/lib/functions/energy/energyfactory.h, lines 31 to 120 · score 0.64 · cellular Potts model, individual cell, algorithms, stencil · Report a mismatch
  6. [6] § Results and discussion › Influence of vessel network topology on tumor growth › Tumor growth behavior ↔ generate_figures.py, lines 269 to 313 · score 0.64 · network length fraction, local network length, local density, tumor cell, logistic, curve · Report a mismatch
  7. [7] § Materials and methods › Model parameters › Cell division rate ↔ src/lib/action/celldivision.h, lines 34 to 88 · score 0.61 · randomly oriented, cell division, probability, surface, ID, properties · Report a mismatch
  8. [8] § Materials and methods › Model parameters › Cell division rate ↔ src/lib/action/celldivision.cpp, lines 26 to 82 · score 0.59 · cell divides, cell division, oriented, probability, voxels, properties · Report a mismatch
  9. [9] § Results and discussion › Vascular data processing for large-scale simulations › Processing into 3D volumes ↔ mousebrainsimulator/src/mousebrainsimulator/main.py, the whole file · a weak match · score 0.57 · binary mask, blood vessel, filter, locations, substacks, resolution · Report a mismatch
  10. [10] § Materials and methods › Loading of blood vessel data into CiS ↔ mousebrainsimulator/src/mousebrainsimulator/main.py, the whole file · a weak match · score 0.56 · binary mask, blood vessel, workers, spaces, split, rank · Report a mismatch
  11. [11] § Results and discussion › Influence of vessel network topology on tumor growth › Tumor growth behavior ↔ generate_figures.py, lines 269 to 313 · score 0.54 · network length fraction, local network length, tumor cell, logistic, curve, density · Report a mismatch

Paper

Loaded from Europe PMC by your browser; OSCR does not store it: doi.org · Europe PMC

The paper is loaded when this pane is shown.

The authors' code

Python · 126 lines · 4.7 KB · no license · 2 matches

  1. import glob
  2. import os
  3. import shutil
  4. import sys
  5. import numpy as np
  6. import json
  7. import jsbeautifier
  8. import pyvista as pv
  9. from copy import deepcopy,copy
  10. from skimage.morphology import skeletonize_3d,binary_dilation,binary_erosion,binary_opening
  11. import matplotlib.pyplot as plt
  12. from scipy.signal import convolve
  13. def loadjson(filename):
  14. with open(filename, "r") as f:
  15. config = json.load(f)
  16. return config
  17. def processed_vtk_to_numpy(pv_obj):
  18. dimensions = pv_obj.dimensions
  19. return np.array(pv_obj["cells"]).reshape(dimensions)
  20. def grid_3d_to_vti(data, filename="test_file.vti"):
  21. data = data.astype(np.uint32)
  22. grid = pv.ImageData()
  23. grid.dimensions = data.shape
  24. grid.origin = (0, 0, 0)
  25. grid.spacing = (1, 1, 1) # (0.65, 0.65, 2)
  26. grid["cells"] = data.ravel(order="C")
  27. grid.save(filename, binary=True)
  28. def mask_3d_array_within_radius(array, position, radius):
  29. # Get the shape of the input array
  30. x_size, y_size, z_size = array.shape
  31. # Create a grid of coordinates for each index in the array
  32. x, y, z = np.indices((x_size, y_size, z_size))
  33. # Calculate the distance from the position (x0, y0, z0)
  34. x0, y0, z0 = position
  35. distance = (x - x0)**2
  36. distance += (y - y0)**2 #np.sqrt((x - x0)**2 + (y - y0)**2 + (z - z0)**2)
  37. distance += (z - z0)**2
  38. distance = np.sqrt(distance)
  39. # Create a mask where the distance is less than or equal to the radius
  40. mask = distance <= radius
  41. # Apply the mask to the array, values outside the radius will be 0
  42. result = np.where(mask, array, 0)
  43. return result, np.sum(mask)
  44. #Define dirs
  45. path_to_vti_dir = "../binary_masks"
  46. simu_dir = "../simulations"
  47. path_to_simus_run_file = "../chunks_to_simulate.txt"
  48. N_per_dim = 2
  49. N_simu_folders = int(N_per_dim**3)
  50. box_padding=250#µm
  51. with open(f"local_vessel_properties.csv","w") as f, open(path_to_simus_run_file, "r") as g:
  52. f.write("#chunk name,folder nbr,vessel voxel count,skeleton voxel count,volume\n")
  53. for line in g:
  54. chunk_name = line[:-5]#Ignore "\n" at EoL
  55. path_to_vti = f"{path_to_vti_dir}/{chunk_name}.vti"
  56. for folder_nbr in range(N_simu_folders):
  57. print(f"Analyzing simu {folder_nbr} of {chunk_name}")
  58. path_to_simu_folder = f"{simu_dir}/{chunk_name}/{folder_nbr:05d}"
  59. config = loadjson(f"{path_to_simu_folder}/config.json")
  60. tumor_entry = config["Filling"]["cells"][1]
  61. vessel_voxel_count = -1
  62. network_length = -1
  63. volume_of_surrounding = -1
  64. tumor_center = copy(tumor_entry["center"])
  65. tumor_center = (tumor_center[0],tumor_center[1],tumor_center[2]) ####<---- CAREFUL
  66. box_start = [tumor_center[i] - box_padding for i in range(3)]
  67. box_end = [tumor_center[i] + box_padding for i in range(3)]
  68. print(box_start,box_end)
  69. whole_vessel_structure = processed_vtk_to_numpy(pv.read(path_to_vti)) == 1
  70. print(whole_vessel_structure.shape)
  71. print("Getting surrounding")
  72. tumor_surrounding = whole_vessel_structure[box_start[0]:box_end[0],box_start[1]:box_end[1],box_start[2]:box_end[2]]
  73. tumor_surrounding_tmp, volume_of_surrounding = mask_3d_array_within_radius(tumor_surrounding, [tumor_center[i]-box_start[i] for i in range(3)], box_padding)
  74. vessel_voxel_count = np.sum(tumor_surrounding_tmp)
  75. print("Skeletonizing surrounding")
  76. # Remove single voxels
  77. print("Opening")
  78. tumor_surrounding = binary_opening(tumor_surrounding)
  79. for erosion_count in range(2):
  80. print(f"Erosion nbr {erosion_count}")
  81. tumor_surrounding = binary_erosion(tumor_surrounding)
  82. tumor_surrounding = binary_opening(tumor_surrounding)
  83. print("Building skeleton...")
  84. tumor_surrounding_skeleton = skeletonize_3d(tumor_surrounding).astype(bool)
  85. print("Skeleton built")
  86. # Remove some remoining artifacts
  87. print("Convoluting...")
  88. kernel = np.ones((3,3,3),dtype= int)
  89. convolution = convolve(tumor_surrounding_skeleton.astype(int), kernel)[1:-1,1:-1,1:-1] * tumor_surrounding_skeleton
  90. tumor_surrounding_skeleton = (convolution == 2) | (convolution == 3)
  91. print("Done")
  92. tumor_surrounding_skeleton, _ = mask_3d_array_within_radius(tumor_surrounding_skeleton, [tumor_center[i]-box_start[i] for i in range(3)], box_padding)
  93. network_length = np.sum(tumor_surrounding_skeleton)
  94. f.write(f"{chunk_name},{folder_nbr:05d},{vessel_voxel_count:1.5f},{network_length:1.5f},{volume_of_surrounding:1.5f}\n")

extract_local_vessel_properties.py at commit c9d7e6d, no license · at the source

Overview

Authors: Eric Behle1, Julian Herold2, Alexander Schug1
ORCID iDs: Alexander Schug
  1. Jülich Supercomputing Centre, Jülich Research Centre, Jülich, North Rhine-Westphalia, Germany
  2. Scientific Computing Center, Karlsruhe Institute of Technology, Eggenstein-Leopoldshafen, Germany
Journal: Biophysical journal, volume 125, issue 7, pages 1701-1712
Dates: received 25 June 2025; accepted 3 March 2026; published online 25 March 2026; in print 7 April 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1016/j.bpj.2026.03.005 · PMID 41887211 · PMCID PMC13351720 · OpenAlex W7140274101
Open access: hybrid, a free copy (OpenAlex)
Status: code verified (how the code is found and verified)
Categories: mouse (organism), other condition (population), computational (subfield)
MeSH: Blood Vessels*, Brain Neoplasms*, Computer Simulation*, Models, Biological*, Neovascularization, Pathologic*, Animals, Brain, Cell Proliferation, Mice (* major topic)
Topic: Angiogenesis and VEGF in Cancer (Molecular Biology, Biochemistry, Genetics and Molecular Biology), according to OpenAlex
Funding: Gauss Centre for Supercomputing; Helmholtz Association; National Institute of Corrections
Citations: not cited yet (Europe PMC); 42 references in the paper

Abstract

Cancer remains a leading cause of mortality. Multidisciplinary studies probe its complex pathology to increase treatment options. Computational modeling of tumor growth on high-performance computing resources offers microscopic insight into its progress and a valuable avenue for advancing our understanding. However, the effective initialization and parameterization of the underlying models require high-resolution data from real tissue structures. Here, we leveraged high-performance computing resources and a massive data set of a mouse brain’s entire vascular network. We processed these image stacks into detailed three-dimensional representations, identified brain regions of interest, and conducted a series of large-scale simulations to investigate how tumor growth is influenced by local vascular network characteristics. By simulating tumor growth with subcellular resolution, we can probe to which extent vessel density and vessel network length influence tumor growth. We determined that vessel density is more likely to affect the growth rate than the network length. Finally, our results allowed us to extrapolate tumor cell growth predictions for the entire mouse brain, highlighting the critical role of vascular topology in tumor progression. Such increasingly realistic simulations of cancer cells and their microenvironment enable researchers to increasingly bridge the gap between basic biology and clinical practice, ultimately supporting the development of more effective personalized cancer therapies.

Reproduced under the paper's license (CC BY), from the paper cited above.

Repositories

Its files are read in the Code ↔ Paper reader above, with 11 matches between paragraphs and lines of code.

nastja/nastja

License: MPL-2.0
State: the link answers, verified on 30 September 2026
Evidence: files inventoried
Commit: d621cd35c5ce10d2cca11d72d91d29e2fc0a9e27, 5 March 2026
Languages: C/C++ (231), C++ (133), Shell (28), Python (17), CUDA (5), C (5)
Size: 679 files, 419 scripts
Software Heritage: archived
Found in: “Data and code availability”
Holds: README, license file, tests, continuous integration, documentation
Not found: CITATION.cff, environment file
Tools: NumPy (9 files), Matplotlib (2 files), pandas (2 files), SciPy (1 file)
Availability: 1 check, the latest on 30 September 2026: the link answers
  • 30 September 2026: the link answers
421 files

gitlab.jsc.fz-juelich.de/behle2/mousebrainvasculaturepaper

License: none: the authors keep all their rights
State: the link answers, verified on 30 September 2026
Evidence: files inventoried
Commit: c9d7e6dfc3f1fd4f31748f59011bf540ee83a9d9, 26 May 2025
Languages: Python (14), Shell (5)
Size: 76 files, 19 scripts
Software Heritage: not archived
Found in: “Data and code availability”
Holds: README, environment (mousebrainsimulator/pyproject.toml, mousebrainsimulator/setup.cfg, mousebrainsimulator/setup.py)
Not found: license file, CITATION.cff, tests, continuous integration, documentation
Tools: NumPy (9 files), Matplotlib (5 files), scikit-image (4 files), SciPy (4 files), pandas (2 files), scikit-learn (1 file)
Availability: 1 check, the latest on 30 September 2026: the link answers
  • 30 September 2026: the link answers
20 files

The paper's code and data availability statement is in the Data section.

Tracing map

Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.

What the map holds:

  • 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 438 scripts, each with its path and the digest of its content;
  • 11 matches between paragraphs of the paper and lines of the code (method lexical-v1);
  • neither the text of the paper nor the code itself.

Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.

Data

No dataset and no data link were found in the paper.

Data and code availability

The code of CiS can be found at https://gitlab.com/nastja/nastja.

The data and analysis scripts required for reproducing our analysis can be found at https://gitlab.jsc.fz-juelich.de/behle2/mousebrainvasculaturepaper.

Due to the size of the entire simulation data (roughly 1.1 TB), we only include the last output file of each simulation within this repository.

The full data are available upon reasonable request.

Reproduced under the paper's license (CC BY), from the paper cited above.

Versions

The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.

Version 1, 30 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 3 authors, 9 MeSH terms, 3 funders, 29 references.

Cite

This paper

Behle, E., Herold, J., & Schug, A. (2026). Leveraging experimental vasculature data for high-resolution brain tumor simulations. Biophysical journal, 125(7), 1701-1712. https://doi.org/10.1016/j.bpj.2026.03.005

BibTeX

@article{behle2026leveraging,
author = {Behle, Eric and Herold, Julian and Schug, Alexander},
title = {{Leveraging experimental vasculature data for high-resolution brain tumor simulations}},
journal = {Biophysical journal},
year = {2026},
month = mar,
volume = {125},
number = {7},
pages = {1701--1712},
publisher = {The Biophysical Society},
issn = {0006-3495},
doi = {10.1016/j.bpj.2026.03.005},
url = {https://doi.org/10.1016/j.bpj.2026.03.005},
pmid = {41887211},
pmcid = {PMC13351720}
}

RIS

TY - JOUR
AU - Behle, Eric
AU - Herold, Julian
AU - Schug, Alexander
TI - Leveraging experimental vasculature data for high-resolution brain tumor simulations
T2 - Biophysical journal
J2 - Biophys J
PY - 2026
DA - 2026/03/25
VL - 125
IS - 7
SP - 1701
EP - 1712
SN - 0006-3495
PB - The Biophysical Society
DO - 10.1016/j.bpj.2026.03.005
UR - https://doi.org/10.1016/j.bpj.2026.03.005
LA - en
ER -

CSL-JSON

{
"id": "10.1016/j.bpj.2026.03.005",
"type": "article-journal",
"title": "Leveraging experimental vasculature data for high-resolution brain tumor simulations",
"container-title": "Biophysical journal",
"author": [
{
"family": "Behle",
"given": "Eric"
},
{
"family": "Herold",
"given": "Julian"
},
{
"family": "Schug",
"given": "Alexander"
}
],
"container-title-short": "Biophys J",
"volume": "125",
"issue": "7",
"page": "1701-1712",
"DOI": "10.1016/j.bpj.2026.03.005",
"PMID": "41887211",
"PMCID": "PMC13351720",
"ISSN": "0006-3495",
"publisher": "The Biophysical Society",
"URL": "https://doi.org/10.1016/j.bpj.2026.03.005",
"language": "en",
"issued": {
"date-parts": [
[
2026,
3,
25
]
]
}
}

The tracing map gets a citation of its own once an author has validated it and it has a DOI.

Similar papers

The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.

[1] doi:10.1038/s41467-022-32848-2 [code]
An analysis modality for vascular structures combining tissue-clearing technology and topological data analysis.
Journal: Nature communications
In common: scikit-image, scikit-learn, pandas, 3 other tools, mouse, 1 reference
[2] doi:10.1126/sciadv.aeb4205 [code]
Learning stochastic dynamics and cell-fate landscapes from single-cell snapshots via optimal transport.
Journal: Science advances
In common: scikit-image, scikit-learn, pandas, 3 other tools, 1 reference
[3] doi:10.1038/s41586-024-07780-8 [code]
Inhibitory specificity from a connectomic census of mouse visual cortex.
Journal: Nature
In common: scikit-image, scikit-learn, pandas, 3 other tools, mouse, 1 reference
[4] doi:10.7554/elife.82376 [code]
Sex, strain, and lateral differences in brain cytoarchitecture across a large mouse population.
Journal: eLife
In common: scikit-image, scikit-learn, pandas, 3 other tools, mouse, 1 reference
[5] doi:10.1016/j.crmeth.2023.100436 [code]
Rapid and fully automated blood vasculature analysis in 3D light-sheet image volumes of different organs.
Journal: Cell reports methods
In common: scikit-image, scikit-learn, SciPy, 2 other tools, mouse, 1 reference
[6] doi:10.1002/mrm.30629 [code]
Fast MR signal simulations of microvascular and diffusion contributions using histogram-based approximation and recurrent neural networks.
Journal: Magnetic resonance in medicine
In common: scikit-learn, pandas, SciPy, 2 other tools, computational, 1 reference
[7] doi:10.1038/s44303-025-00081-3 [code]
Temperature dependent chemical shifts of pyruvate and lactate enable in vivo hyperpolarized &lt;sup&gt;13&lt;/sup&gt;C MRSI thermometry.
Journal: Npj imaging
In common: scikit-image, pandas, SciPy, 2 other tools, 1 reference
[8] doi:10.1126/sciadv.abm6340 [code]
Neural network learning defines glioblastoma features to be of neural crest perivascular or radial glia lineages.
Journal: Science advances
In common: scikit-learn, pandas, SciPy, 2 other tools, other condition, 1 reference
[9] doi:10.7554/elife.83792 [code]
Dynamic fibronectin assembly and remodeling by leader neural crest cells prevents jamming in collective cell migration.
Journal: eLife
In common: scikit-learn, pandas, SciPy, 2 other tools, 1 reference
[10] doi:10.1101/2025.07.16.665209 [code]
Dual-feature selectivity enables bidirectional coding in visual cortical neurons
Journal: bioRxiv (preprint)
In common: scikit-image, pandas, SciPy, 2 other tools, 1 reference

Contribute

The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.

Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.

Add the badge to the code's README

OSCR: the badge of this paper

If you are an author or a maintainer of this paper's code, you can add this badge. It says what OSCR knows of the paper (its code, its validated tracing map) and links the code to this page. Paste one of these snippets under the README's title. Keep the two marks around it, so that the badge can be replaced in place and is never added twice. OSCR changes nothing in your repository.

Markdown
<!-- oscr-badge:start -->
[![OSCR](https://openscicode.org/badge/doi/10.1016/j.bpj.2026.03.005.svg)](https://openscicode.org/badge/doi/10.1016/j.bpj.2026.03.005)
<!-- oscr-badge:end -->
reStructuredText
.. oscr-badge:start

.. image:: https://openscicode.org/badge/doi/10.1016/j.bpj.2026.03.005.svg
   :target: https://openscicode.org/badge/doi/10.1016/j.bpj.2026.03.005
   :alt: OSCR

.. oscr-badge:end
HTML
<!-- oscr-badge:start -->
<a href="https://openscicode.org/badge/doi/10.1016/j.bpj.2026.03.005"><img src="https://openscicode.org/badge/doi/10.1016/j.bpj.2026.03.005.svg" alt="OSCR"></a>
<!-- oscr-badge:end -->

If you choose to, you can have the badge kept up to date for you, with oscr badge --write from the command line or with the GitHub Action in your repository's workflows. Every other form and option (a card, AsciiDoc, LaTeX, shields.io) is on the badge page.

Request its removal

To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).

Discussion, reproductions, activity

Discussion: code errors and mismatches between the code and the paper about this paper and its code, from readers and its authors, beside the paper and its tracing map.

Report a mismatch between the code and the paper · Report a code error (signed in)

Reproductions: reports from readers who ran the authors' code and did not get the paper's result, what they ran, with which environment, commit and data.

Report a reproduction failure (signed in)

Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.