OSCR

Quantification of Ki-67 labeling index in pediatric brain tumor immunohistochemistry images.

Code ↔ Paper

2 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 2 matches
  1. [1] § RESULTS › Density maps ↔ summary_ratios.py, lines 106–146 · score 0.61 · QuPath, LI map, positive cell density, negative cell density, ratio, slide
  2. [2] § METHODS › Postprocessing in Python ↔ summary_ratios.py, lines 106–146 · score 0.60 · QuPath, LI map, cell density maps, ratio, negative cell, Ki

Paper

Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC

The paper is loaded when this pane is shown.

The authors' code

Python · 389 lines · 16 KB · GPL-3.0 · 2 matches

  1. # %% IMPORTS
  2. import json
  3. import os
  4. import numpy as np
  5. import math
  6. import pandas as pd
  7. import argparse
  8. import matplotlib
  9. import matplotlib.pyplot as plt
  10. import seaborn as sns
  11. from PIL import Image
  12. from mpl_toolkits.axes_grid1 import make_axes_locatable
  13. # %% ARGUMENT PARSING
  14. parser = argparse.ArgumentParser(description='Create summary ratios and density maps.')
  15. parser.add_argument('--maps_dir', type=str, default=False, help='Directory of the cell density maps generated by QuPath.')
  16. parser.add_argument('--data_dir', type=str, default=False, help='Directory of the project data')
  17. parser.add_argument('--area_path', type=str, default=False, help='Path to Area.txt which is saved under the Results folder.')
  18. parser.add_argument('--csv_path', type=str, required=True, help='Path to the WSI dataframe CSV file')
  19. parser.add_argument('--WSIs_dir', type=str, default=False, help='Directory of the folder where the WSIs exist.')
  20. parser.add_argument('--norm_maps_dir', type=str, default=False, help='Directory to save the normalized cell density maps.')
  21. parser.add_argument('--result_dir', type=str, default=False, help='Directory to save the summary tables and graphs.')
  22. args = parser.parse_args()
  23. # %% DIRECTORIES & FILES NEEDED FOR THE ANALYSIS
  24. maps_dir = args.maps_dir
  25. data_dir = args.data_dir
  26. csv_path = args.csv_path
  27. # check if the CSV file exists
  28. if not os.path.isfile(csv_path):
  29. raise FileNotFoundError(f"The specified CSV file does not exist: {csv_path}")
  30. area_path = args.area_path
  31. # check if the Area.txt file exists
  32. if not os.path.isfile(area_path):
  33. raise FileNotFoundError(f"The specified Area.txt file does not exist: {area_path}")
  34. WSI_df = pd.read_csv(csv_path)
  35. WSIs_dir = args.WSIs_dir
  36. norm_maps_dir = args.norm_maps_dir
  37. # create the directory to save the normalised maps if it does not exist
  38. if not os.path.exists(norm_maps_dir):
  39. os.makedirs(norm_maps_dir)
  40. result_dir = args.result_dir
  41. # create the directory to save the results if it does not exist
  42. if not os.path.exists(result_dir):
  43. os.makedirs(result_dir)
  44. # get the area and total detections
  45. annotation_df = pd.read_csv(area_path,names=['name', 'detections', 'area'],delimiter=";",index_col=False)
  46. # %% RETRIEVE INFORMATION FROM JSON FILES
  47. file_list = [f for f in os.listdir(data_dir) if os.path.isdir(os.path.join(data_dir,f))]
  48. WSIs = os.listdir(WSIs_dir)
  49. # data frame to store the results
  50. results_df = pd.DataFrame(columns = ['case_id', 'slide_id', 'label', 'Positive_Cell_Count', 'Positive_Cell_Density', 'Negative_Cell_Count', 'Negative_Cell_Density', 'Ki-67 LI', 'Cell_Density'])
  51. proc_l = len(file_list)
  52. proc = -1
  53. failed_list = []
  54. failed_ann = []
  55. for folder in file_list:
  56. proc += 1
  57. proc_perc = round(proc/proc_l*100, 2)
  58. print('Progress: ' + str(proc_perc) + '%')
  59. print('Processing image ' + str(proc) + ' of ' + str(proc_l))
  60. with open(data_dir + '/' + folder + '/server.json', 'r') as s_file:
  61. server_data = json.load(s_file)
  62. slide_id = server_data['builder']['metadata']['name'][:-4]
  63. case_id = slide_id.split('___')[0]
  64. # Convert slide_id to int to match CSV data type
  65. try:
  66. slide_id_lookup = int(slide_id)
  67. except ValueError:
  68. # If conversion fails, keep as string (for cases with non-numeric IDs)
  69. slide_id_lookup = slide_id
  70. label = WSI_df['label'].loc[WSI_df['slide_id']==slide_id_lookup].values[0]
  71. print('Reading data from image ' + slide_id + ':')
  72. with open(data_dir + '/' + folder + '/summary.json', 'r') as file:
  73. data = json.load(file)
  74. detection_data = data['hierarchy']['detectionClassificationCounts']
  75. if 'Positive' not in detection_data.keys() and 'Negative' not in detection_data.keys():
  76. print('The slide has not been properly segmented yet, skipping.')
  77. failed_list.append(slide_id)
  78. continue
  79. # calculate the Ki-67 LI
  80. pos_cell_count = detection_data.get('Positive', 0)
  81. neg_cell_count = detection_data.get('Negative', 0)
  82. # Check if there are any cells detected
  83. if pos_cell_count + neg_cell_count == 0:
  84. print('No cells detected in this slide, skipping.')
  85. failed_list.append(slide_id)
  86. continue
  87. # round ratio to 2 decimal places
  88. ki67_li = round(pos_cell_count / (pos_cell_count + neg_cell_count), 2)
  89. for WSI in WSIs:
  90. WSI_name = WSI.split('.')[0]
  91. if slide_id in WSI_name:
  92. slide_file_name = WSI
  93. print('Found the corresponding WSI file: ' + slide_file_name)
  94. break
  95. try:
  96. annAreamm = annotation_df['area'].loc[annotation_df['name'] == slide_file_name].values[0]*0.000001
  97. numdet = annotation_df['detections'].loc[annotation_df['name'] == slide_file_name].values[0]
  98. cells_smm = numdet/annAreamm
  99. pos_cell_dens = pos_cell_count/annAreamm
  100. neg_cell_dens = neg_cell_count/annAreamm
  101. except:
  102. print('Could not find any annotation data.')
  103. failed_ann.append(slide_file_name)
  104. else:
  105. print("Area of the segmented tissue: " + str(annAreamm))
  106. print("Total number of detections: " + str(numdet))
  107. print('Number of detections per mm^2: ' + str(cells_smm))
  108. results_df.loc[len(results_df)] = [case_id, slide_id, label, pos_cell_count, pos_cell_dens, neg_cell_count, neg_cell_dens, ki67_li*100, cells_smm]
  109. print('Positive cell count: ' + str(pos_cell_count))
  110. print('Negative cell count: ' + str(neg_cell_count))
  111. print('Positive to negative ratio: ' + str(ki67_li))
  112. # %% GENERATE THE NORMALISED AND RATIO CELL DENSITY MAPS
  113. # get the cell density maps generated by QuPath
  114. HM_file_neg = maps_dir + '/' + slide_id + '_NegDMap.tif'
  115. HM_file_pos = maps_dir + '/' + slide_id + '_PosDMap.tif'
  116. failed_cell_density_list = []
  117. # check if both cell density map files exist before proceeding
  118. if os.path.isfile(HM_file_neg) and os.path.isfile(HM_file_pos):
  119. neg_pp_im = norm_maps_dir + '/' + slide_id + '_Neg_Cell_Dens_Map.png'
  120. pos_pp_im = norm_maps_dir + '/' + slide_id + '_Pos_Cell_Dens_Map.png'
  121. ratio_pp_im = norm_maps_dir + '/' + slide_id + '_Ki-67_LI_Map.png'
  122. # generate the normalised negative cell density map
  123. if not os.path.isfile(neg_pp_im):
  124. print("Generating normalized negative cell density map.")
  125. # get the raw grayscale negative cell density map
  126. im_neg = Image.open(HM_file_neg)
  127. width, height = im_neg.size
  128. # set color bar orientation depending on image dimensions:
  129. if width < height:
  130. orient = 'vertical'
  131. else:
  132. orient = 'horizontal'
  133. # normalize
  134. imarray_neg = np.array(im_neg)
  135. cm_lim_neg = np.max(imarray_neg)
  136. new_imarray_neg = imarray_neg/cm_lim_neg
  137. # keep values inside normalized range for a correct comparison between heat maps
  138. new_imarray_neg[new_imarray_neg>1] = 1
  139. # colorise the negative cell density map
  140. cmapQ = plt.get_cmap('jet')
  141. A_neg = cmapQ(new_imarray_neg)
  142. # convert dark blue background to white
  143. A_neg[new_imarray_neg == 0] = [1, 1, 1, 1]
  144. # plot and save normalised negative cell density map
  145. fig, ax = plt.subplots()
  146. im = ax.imshow(A_neg)
  147. ax.set_xticks([])
  148. ax.set_yticks([])
  149. ax.set_title('Negative Cell Density Map')
  150. divider = make_axes_locatable(ax)
  151. if orient == 'vertical':
  152. cax = divider.append_axes("right", size="5%", pad=0.20)
  153. else:
  154. cax = divider.append_axes("bottom", size="5%", pad=0.20)
  155. norm = matplotlib.colors.Normalize(vmin=0, vmax=neg_cell_dens)
  156. cb = plt.colorbar(plt.cm.ScalarMappable(norm=norm, cmap='jet'), cax=cax, orientation=orient, shrink=1.0)
  157. plt.title('Negative Cell Density Map')
  158. plt.savefig(neg_pp_im)
  159. cb.remove()
  160. plt.close("all")
  161. # mask of the negative cell density map
  162. A_neg_blue_mask = A_neg
  163. A_neg_blue_mask[np.any(A_neg_blue_mask[:, :, :3] != [1, 1, 1], axis=2)] = [0, 0, 1, 1]
  164. else:
  165. print("Normalized negative cell density map already exists, skipping.")
  166. # generate the normalised positive cell density map
  167. if not os.path.isfile(pos_pp_im):
  168. print("Generating normalized positive cell density map.")
  169. # get the raw grayscale positive cell density map
  170. im_pos = Image.open(HM_file_pos)
  171. width, height = im_pos.size
  172. # set color bar orientation depending on image dimensions:
  173. if width < height:
  174. orient = 'vertical'
  175. else:
  176. orient = 'horizontal'
  177. # normalise
  178. imarray_pos = np.array(im_pos)
  179. cm_lim_pos = np.max(imarray_pos)
  180. new_imarray_pos = imarray_pos/cm_lim_pos
  181. # keep values inside normalized range for a correct comparison between heat maps
  182. new_imarray_pos[new_imarray_pos>1] = 1
  183. # colorise
  184. cmapQ = plt.get_cmap('jet')
  185. A_pos = cmapQ(new_imarray_pos)
  186. # convert dark blue background to white
  187. A_pos[new_imarray_pos == 0] = [1, 1, 1, 1]
  188. new_A_pos = A_pos
  189. white_mask = np.all(new_A_pos[:, :, :3] == [1, 1, 1], axis=2)
  190. new_A_pos[white_mask] = A_neg_blue_mask[white_mask]
  191. # plot and save normalised positive cell density map
  192. fig, ax = plt.subplots()
  193. im = ax.imshow(new_A_pos)
  194. ax.set_xticks([])
  195. ax.set_yticks([])
  196. ax.set_title('Positive Cell Density Map')
  197. divider = make_axes_locatable(ax)
  198. if orient == 'vertical':
  199. cax = divider.append_axes("right", size="5%", pad=0.20)
  200. else:
  201. cax = divider.append_axes("bottom", size="5%", pad=0.20)
  202. norm = matplotlib.colors.Normalize(vmin=0, vmax=pos_cell_dens)
  203. cb = plt.colorbar(plt.cm.ScalarMappable(norm=norm, cmap='jet'), cax=cax, orientation=orient, shrink=1.0)
  204. plt.savefig(pos_pp_im)
  205. cb.remove()
  206. plt.close("all")
  207. else:
  208. print("Normalized positive cell density map already exists, skipping.")
  209. # generate the Ki-67 LI map
  210. if not os.path.isfile(ratio_pp_im):
  211. print("Generating ratio cell density map.")
  212. neg_im = Image.open(HM_file_neg)
  213. width, height = neg_im.size
  214. # set color bar orientation depending on image dimensions:
  215. if width < height:
  216. orient = 'vertical'
  217. else:
  218. orient = 'horizontal'
  219. neg_imarray = np.array(neg_im)
  220. pos_im = Image.open(HM_file_pos)
  221. pos_imarray = np.array(pos_im)
  222. for k in range(np.shape(neg_imarray)[0]):
  223. for l in range(np.shape(neg_imarray)[1]):
  224. denominator = neg_imarray[k, l] + pos_imarray[k, l]
  225. # if denominator is 0, assign a default value
  226. if denominator == 0:
  227. neg_imarray[k, l] = 0.0001
  228. else:
  229. neg_imarray[k, l] = pos_imarray[k, l] / denominator
  230. if neg_imarray[k, l] < 0.0001:
  231. neg_imarray[k, l] = 0.0001
  232. elif neg_imarray[k, l] > 1:
  233. neg_imarray[k, l] = 1
  234. elif math.isnan(neg_imarray[k, l]):
  235. neg_imarray[k, l] = 0.0001
  236. # only take into account foreground pixels
  237. relevant_values = neg_imarray[neg_imarray > 0]
  238. n_relval = len(relevant_values)
  239. adjusted_n_relval = round(n_relval*0.95)
  240. sorted_relval = sorted(relevant_values)
  241. min_cm_lim = 0.03
  242. cm_lim = max(max(sorted_relval[:adjusted_n_relval]), min_cm_lim)
  243. # normalize with obtained value
  244. neg_imarray = neg_imarray/cm_lim
  245. # keep values inside normalized range for a correct comparison between heat maps
  246. neg_imarray[neg_imarray>1] = 1
  247. # colorise
  248. cmapQ = plt.get_cmap('jet')
  249. A_ki67_li = cmapQ(neg_imarray)
  250. # convert dark blue background to white
  251. dark_blue = [0, 0, 0.5, 1] # RGBA values for dark blue
  252. white = [1, 1, 1, 1] # RGBA values for white
  253. A_ki67_li[
  254. (A_ki67_li[:, :, 0] == dark_blue[0]) &
  255. (A_ki67_li[:, :, 1] == dark_blue[1]) &
  256. (A_ki67_li[:, :, 2] == dark_blue[2]) &
  257. (A_ki67_li[:, :, 3] == dark_blue[3])
  258. ] = white
  259. new_A_ki67_li = A_ki67_li
  260. white_mask = np.all(new_A_ki67_li[:, :, :3] == [1, 1, 1], axis=2)
  261. new_A_ki67_li[white_mask] = A_neg_blue_mask[white_mask]
  262. # plot and save the ki-67 li map
  263. fig, ax = plt.subplots()
  264. im = ax.imshow(new_A_ki67_li)
  265. ax.set_xticks([])
  266. ax.set_yticks([])
  267. ax.set_title('Ki-67 LI Map')
  268. divider = make_axes_locatable(ax)
  269. if orient == 'vertical':
  270. cax = divider.append_axes("right", size="5%", pad=0.20)
  271. else:
  272. cax = divider.append_axes("bottom", size="5%", pad=0.20)
  273. norm = matplotlib.colors.Normalize(vmin=0, vmax=ki67_li)
  274. cb = plt.colorbar(plt.cm.ScalarMappable(norm=norm, cmap='jet'), cax=cax, orientation=orient, shrink=1.0)
  275. plt.savefig(ratio_pp_im)
  276. cb.remove()
  277. plt.close("all")
  278. else:
  279. print("Normalized Ki-67 LI map already exists, skipping.")
  280. else:
  281. print(f"Required cell density map files not found for slide ID: {slide_id}. Skipping.")
  282. # %% CREATE SUMMARY STATISTICS
  283. # save the results to a csv file
  284. results_df.to_csv(os.path.join(result_dir, 'summary_results.csv'), index=False)
  285. statistics_df = pd.DataFrame(columns = ['variable', 'label', 'mean', 'sd', 'median', 'minimum', 'maximum'])
  286. summary_points = ['Positive_Cell_Count', 'Positive_Cell_Density', 'Negative_Cell_Count', 'Negative_Cell_Density', 'Ki-67 LI', 'Cell_Density']
  287. y_units = {'Positive_Cell_Count': 'Number of positive cells',
  288. 'Positive_Cell_Density': 'Number of positive cells per mm^2',
  289. 'Negative_Cell_Count': 'Number of negative cells',
  290. 'Negative_Cell_Density': 'Number of negative cells per mm^2',
  291. 'Ki-67 LI': 'Ki-67 LI',
  292. 'Cell_Density': 'Number of cells per mm^2'}
  293. res_df = pd.read_csv(os.path.join(result_dir, 'summary_results.csv'))
  294. label_list = np.unique(res_df['label'].tolist())
  295. # summary table and boxplots
  296. for sp in summary_points:
  297. box_data = []
  298. for label in label_list:
  299. sp_array = res_df[sp].loc[res_df['label'] == label]
  300. sp_array = sp_array[~np.isnan(sp_array)]
  301. box_data.append(sp_array)
  302. sp_mean = round(np.mean(sp_array), 2)
  303. sp_sd = round(np.std(sp_array), 2)
  304. sp_median = round(np.median(sp_array), 2)
  305. sp_min = round(np.min(sp_array), 2)
  306. sp_max = round(np.max(sp_array), 2)
  307. # summary statistics
  308. statistics_df.loc[len(statistics_df)] = [sp, label, sp_mean, sp_sd, sp_median, sp_min, sp_max]
  309. # boxplots
  310. palette = sns.color_palette("deep", 10)
  311. custom_colors = [palette[9]]
  312. palette = sns.color_palette(custom_colors * len(label_list))
  313. sns.set_style("whitegrid", {'axes.grid': False})
  314. plt.figure(figsize=(10, 6))
  315. plt.gca().set_facecolor('white')
  316. flierprops = dict(marker='D', markerfacecolor='darkgrey', markersize=5, linestyle='none')
  317. boxplot = sns.boxplot(
  318. x='label',
  319. y=sp,
  320. data=res_df,
  321. flierprops=flierprops,
  322. palette=palette,
  323. linewidth=2
  324. )
  325. plt.title(f'{sp} across all diagnoses')
  326. plt.xlabel('Diagnosis')
  327. plt.ylabel(y_units[sp])
  328. plt.tight_layout()
  329. plt.savefig(os.path.join(result_dir,sp + '_BoxPlot.png'))
  330. plt.close()
  331. # summary statistics
  332. statistics_df.to_csv(os.path.join(result_dir, 'summary_statistics.csv'), index=False)
  333. # %%

summary_ratios.py at commit d1efa7c, under GPL-3.0 · at the source

Overview

Authors: Christoforos Spyretos1,2, Juan Manuel Pardo Ladino1, Hakon Andersen Blomstrand3,4, Per Nyman2,5,6, Oscar Snödahl2,6,7, Alia Shamikh8,9, Nils Elander4,10, Neda Haj-Hosseini1,2
  1. Department of Biomedical Engineering, Linköping University, Linköping, Sweden
  2. Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden
  3. Clinical Department of Clinical Pathology, Region Östergötland, Linköping, Sweden
  4. Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
  5. Crown Princess Victoria Children's Hospital, Region Östergötland, Linköping, Sweden
  6. Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
  7. Clinical Department of Radiology in Linköping, Region Östergötland, Linköping, Sweden
  8. Department of Clinical Pathology and Cancer Diagnostics, Karolinska University Hospital, Solna, Sweden
  9. Department of Oncology-Pathology, Karolinska Institute, Solna, Sweden
  10. Clinical Department of Oncology in Linköping, Region Östergötland, Linköping, Sweden
Journal: Journal of neuropathology and experimental neurology, volume 85, issue 5, pages 475-486
Dates: published online 10 March 2026; in print May 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1093/jnen/nlaf163 · PMID 41806389 · PMCID PMC13127889 · OpenAlex W7134976525
Open access: hybrid, a free copy (OpenAlex)
Status: code verified
Categories: histology / microscopy (modality), human (organism), other condition (population)
Methods: Connectivity, Statistics
Keywords: pediatric brain tumor, histopathology, immunohistochemistry, Ki-67 labeling index, deep learning, QuPath, cell segmentation, cell classification, cell density map
MeSH: Brain Neoplasms*, Ki-67 Antigen*, Biomarkers, Tumor, Child, Child, Preschool, Female, Glioma, Humans, Immunohistochemistry, Infant, Male (* major topic)
Topic: Glioma Diagnosis and Treatment (Genetics, Medicine), according to OpenAlex
Funding: Linköping University's Cancer Strength Area (2024), and the Medical Research Council of Southeast Sweden (FORSS-1011571); Children's Brain Tumor Network; Children’s Brain Tumor Tissue Consortium; Medical Research Council of Southeast Sweden (FORSS-1011571); Joanna Cocozza's Foundation (2025-2026); Swedish Childhood Cancer Foundation (MT2021-0011, MT2022-0013); Linköping University’s Cancer Strength Area; Children's Brain Tumor Tissue Consortium/the Children's Brain Tumor Network
Citations: cited by 1 paper (Europe PMC); 40 references in the paper

Abstract

Quantification of the Kiel 67 (Ki-67) labeling index (LI) is critical for assessing proliferation and prognosis in tumors but manual scoring remains a common practice. We present an automated framework for Ki-67 scoring in whole slide images (WSIs) developed for research settings using an Apache Groovy code script for QuPath and complemented by a Python postprocessing script that provides cell density maps and summary tables. Tissue segmentation is performed by pixel classifiers and cell segmentation is conducted using StarDist, a deep learning model, followed by adaptive thresholding to classify Ki-67 positive and negative nuclei. The pipeline was applied to a cohort of 632 pediatric brain tumor cases with 734 Ki-67 WSIs from the Children’s Brain Tumor Network. Medulloblastomas showed the highest Ki-67 LI (median: 19.84), followed by atypical teratoid rhabdoid tumors (median: 19.36), brainstem glioma-diffuse intrinsic pontine gliomas (median: 11.50), high-grade gliomas (grades 3, 4) (median: 9.50), and ependymomas (median: 5.88). Lower indices were found in meningiomas (median: 1.84) and the lowest were seen in low-grade gliomas (grades 1, 2) (median: 0.85), dysembryoplastic neuroepithelial tumors (median: 0.63), and gangliogliomas (median: 0.50). The results demonstrate a significant correlation (P < .05) in Ki-67 LI across most of the tumor families/types aligning with neuro-oncology and neuropathology consensus.

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

Repository

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

Christoforos-Spyretos/QuPath-Automatic-Cell-Detection-for-Ki-67-WSIs

License: GPL-3.0
State: the link answers, verified on 28 September 2026
Evidence: files inventoried
Commit: d1efa7cff94c115107a0e4d9790d606b0eba0134, 21 May 2026
Languages: Python (1)
Size: 13 files, 1 script
Software Heritage: not archived
Found in: “DATA AVAILABILITY”
Holds: README, license file, environment (requirements.txt)
Not found: CITATION.cff, tests, continuous integration, documentation
Tools: Matplotlib (1 file), NumPy (1 file), pandas (1 file), Pillow (1 file), seaborn (1 file)
Availability: 1 check, the latest on 28 September 2026: the link answers
  • 28 September 2026: the link answers
3 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:

  • 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 1 script, each with its path and the digest of its content;
  • 2 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

Datasets cited

Data availability

The open-access dataset used in this study was obtained from the Children’s Brain Tumor Network (https://cbtn.org). The dataset for neuroendocrine and testicular seminoma are open-access and are available at https://github.com/cialab/neuroendocrine_ and https://zenodo.org/records/11218961, respectively. The code and detailed instructions for execution linked to this study are available in the GitHub repository: https://github.com/Christoforos-Spyretos/QuPath-Automatic-Cell-Detection-for-Ki-67-WSIs

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, 28 September 2026: the first record

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

Cite

This paper

Spyretos, C., Pardo Ladino, J. M., Andersen Blomstrand, H., Nyman, P., Snödahl, O., Shamikh, A., Elander, N., & Haj-Hosseini, N. (2026). Quantification of Ki-67 labeling index in pediatric brain tumor immunohistochemistry images. Journal of neuropathology and experimental neurology, 85(5), 475-486. https://doi.org/10.1093/jnen/nlaf163

BibTeX

@article{spyretos2026quantification,
author = {Spyretos, Christoforos and Pardo Ladino, Juan Manuel and Andersen Blomstrand, Hakon and Nyman, Per and Snödahl, Oscar and Shamikh, Alia and Elander, Nils and Haj-Hosseini, Neda},
title = {{Quantification of Ki-67 labeling index in pediatric brain tumor immunohistochemistry images}},
journal = {Journal of neuropathology and experimental neurology},
year = {2026},
month = may,
volume = {85},
number = {5},
pages = {475--486},
publisher = {Oxford University Press},
issn = {0022-3069},
doi = {10.1093/jnen/nlaf163},
url = {https://doi.org/10.1093/jnen/nlaf163},
pmid = {41806389},
pmcid = {PMC13127889}
}

RIS

TY - JOUR
AU - Spyretos, Christoforos
AU - Pardo Ladino, Juan Manuel
AU - Andersen Blomstrand, Hakon
AU - Nyman, Per
AU - Snödahl, Oscar
AU - Shamikh, Alia
AU - Elander, Nils
AU - Haj-Hosseini, Neda
TI - Quantification of Ki-67 labeling index in pediatric brain tumor immunohistochemistry images
T2 - Journal of neuropathology and experimental neurology
J2 - J Neuropathol Exp Neurol
PY - 2026
DA - 2026/05/01
VL - 85
IS - 5
SP - 475
EP - 486
SN - 0022-3069
PB - Oxford University Press
DO - 10.1093/jnen/nlaf163
UR - https://doi.org/10.1093/jnen/nlaf163
LA - en
ER -

CSL-JSON

{
"id": "10.1093/jnen/nlaf163",
"type": "article-journal",
"title": "Quantification of Ki-67 labeling index in pediatric brain tumor immunohistochemistry images",
"container-title": "Journal of neuropathology and experimental neurology",
"author": [
{
"family": "Spyretos",
"given": "Christoforos"
},
{
"family": "Pardo Ladino",
"given": "Juan Manuel"
},
{
"family": "Andersen Blomstrand",
"given": "Hakon"
},
{
"family": "Nyman",
"given": "Per"
},
{
"family": "Snödahl",
"given": "Oscar"
},
{
"family": "Shamikh",
"given": "Alia"
},
{
"family": "Elander",
"given": "Nils"
},
{
"family": "Haj-Hosseini",
"given": "Neda"
}
],
"container-title-short": "J Neuropathol Exp Neurol",
"volume": "85",
"issue": "5",
"page": "475-486",
"DOI": "10.1093/jnen/nlaf163",
"PMID": "41806389",
"PMCID": "PMC13127889",
"ISSN": "0022-3069",
"publisher": "Oxford University Press",
"URL": "https://doi.org/10.1093/jnen/nlaf163",
"language": "en",
"issued": {
"date-parts": [
[
2026,
5,
1
]
]
}
}

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.1002/alz.71649 [code]
Postmortem brain MRI reveals differential associations of subcortical and limbic volumes with cortical thinning and neurodegenerative pathologies.
Journal: Alzheimer's & dementia : the journal of the Alzheimer's Association
In common: Pillow, seaborn, pandas, 2 other tools, histology / microscopy, other condition
[2] doi:10.1371/journal.pcbi.1014571 [code]
SynAPSeg: A novel dataset and image analysis framework for deep learning-based synapse detection and quantification.
Journal: PLoS computational biology
In common: Pillow, seaborn, pandas, 2 other tools, 1 reference
[3] doi:10.1038/s41593-026-02388-9 [code]
Hippocampal CA3 connectomics reveals a gradient of mossy fiber inputs and selective feedforward inhibition onto pyramidal cells.
Journal: Nature neuroscience
In common: Pillow, seaborn, pandas, 2 other tools, histology / microscopy
[4] doi:10.1016/j.xpro.2026.104659 [code]
Protocol for simultaneous in vivo two-photon imaging and locomotion quantification during olfactory stimulation and pharmacology in walking Drosophila.
Journal: STAR protocols
In common: Pillow, seaborn, pandas, 2 other tools, histology / microscopy
[5] doi:10.1038/s41467-026-72709-w [code]
An epifluorescence microscope design for naturalistic behavior and cellular activity in freely moving Caenorhabditis elegans.
Journal: Nature communications
In common: Pillow, seaborn, pandas, 2 other tools, histology / microscopy
[6] doi:10.1186/s13244-026-02365-7 [code]
Super-resolution MRI and 2.5D deep learning for intratumoral-peritumoral radiomics in preoperative prediction of rectal cancer perineural invasion.
Journal: Insights into imaging
In common: Pillow, seaborn, pandas, 2 other tools, other condition
[7] doi:10.1038/s41467-026-76837-1 [code]
Drug screen and machine learning predict neuroprotective agents in a preclinical human model of childhood dementia.
Journal: Nature communications
In common: Pillow, seaborn, pandas, 2 other tools, other condition
[8] doi:10.1523/eneuro.0041-26.2026 [code]
Ocular Speech Tracking Persists in Blindness, but Its Dynamics and Oculo-Cerebral Connectivity Depend on Visual Status.
Journal: eNeuro
In common: Pillow, seaborn, pandas, 2 other tools, other condition
[9] doi:10.1162/imag.a.1264 [code]
Data-driven subtyping and staging of ALS: A multicenter, longitudinal, deformation-based morphometry study.
Journal: Imaging neuroscience (Cambridge, Mass.)
In common: Pillow, seaborn, pandas, 2 other tools, other condition
[10] doi:10.1371/journal.pmed.1004809 [code]
Brain morphology in Anorexia Nervosa and its subtypes: A multi-cohort study of individual participant data.
Journal: PLoS medicine
In common: Pillow, seaborn, pandas, 2 other tools, other condition

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.

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: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.

Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.

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