Nicola Rossberg
a,
Simon Corrie
b,
Lisbeth Grøndahl
c and
Imanda Jayawardena
*d
aTaighde Éireann – Research Ireland Center for Research Training in Artificial Intelligence, School of Computer Science & Information Technology, University College Cork, Co. Cork T12 K8AF, Ireland
bChemical and Biological Engineering Department, ARC Centre of Excellence in Carbon Science and Innovation, Monash University, Clayton, VIC 3800, Australia
cSchool of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, QLD 4072, Australia
dTyndall National Institute, Lee Maltings Complex, Dyke Parade, Cork, Co. Cork T12R5CP, Ireland. E-mail: imanda.jayawardena@tyndall.ie
First published on 10th July 2025
Quantitative assessment of pore size and morphology is crucial in biomaterials design and evaluation, particularly hydrogels and scaffolds used in tissue engineering and drug delivery. In recent years, a growing number of studies have proposed or adopted automated image analysis tools to evaluate pore characteristics; however, the absence of standardised protocols, validation criteria, and consistent reporting practices has limited reproducibility and cross-study comparability. This perspective, for the first time, examines recent trends in automated pore size analysis in biomaterials research, highlighting commonly used algorithms, their implementation in image-based workflows, and their ability to resolve pore geometries in disordered materials. We discuss the influence of imaging dimension, resolution, algorithm assumptions, and image pre-processing on outcomes and highlight common challenges such as over-segmentation, user bias, and the misidentification of irregularly shaped pores. By drawing on selected examples from the literature, we illustrate both the strengths and limitations of current approaches and emphasise the need for transparent, standardised methodologies in the field.
In hydrogel-based materials, mesh size (the space between crosslinking sites in polymer chains) and mesh radius (size of largest spherical solute that could move through a mesh portal) governs mechanical properties and permeability of the material, whereas pore size (the voids within e.g. hydrogel network or scaffold) plays a crucial role in cellular infiltration, nutrient diffusion, and tissue regeneration.1,5,7,12–17 Hydrogels with small mesh sizes provide higher mechanical stability but may restrict cellular infiltration, whereas larger pores enhance cell migration and diffusion but can compromise mechanical integrity.1,7,18 Similarly, scaffold-based biomaterials rely on optimised pore architectures to balance mechanical support with biological functionality.1,5,8,19
Traditional methods for characterising pore structures, such as manual analysis from images are time-intensive, and are prone to over or under estimation of parameters, subjectivity and bias.1,3,5,16,20 These methods often fail to provide comprehensive characterisation, particularly for hydrogels, which have hydrated and dynamic porous networks. While semi-automated approaches have been developed, limitations such as bias during user input, low sample sizes during verification steps and lack of accessibility to both software and sample measurement techniques has necessitated automation.3,5,20,21 As a result, there has been an increase in adopting automated pore analysis techniques that utilise advanced image processing and artificial intelligence (AI) algorithms.1,10,20,22–24 Automated segmentation and feature extraction enable high-throughput, reproducible, and quantitative analysis of pore architectures across various length scales.25 Furthermore, automated models can enhance image contrast, remove noise, and classify pore structures with minimal human intervention.7
Numerous studies in the literature develop and/or use automated methods for pore size analysis. However, in the absence of standardised methodologies, reporting guidelines, or benchmarking criteria, comparison across different materials becomes challenging. It is therefore important to evaluate studies from the current literature to illustrate such limitations and strengths. This perspective is, to our knowledge, the first to provide a focused and critical evaluation of automated pore size analysis methods as applied to biomaterials such as hydrogels and porous scaffolds. We explore various automated approaches used to quantify pore characteristics from images recorded using direct two-dimensional (2D) and three-dimensional (3D) imaging methods, discuss the challenges associated with analysing disordered porous materials, and provide recommendations for future developments in automated pore analysis.
In polymeric biomaterials, such as hydrogels, mesh size refers to the distance between crosslinking sites in polymer chains, mesh radius28 describes the voids within a polymeric network, while pore size2 is the area of void space within solid material where the structures separating the void spaces (e.g. pore walls) consists of polymer bundles. The latter description relates to the pores that can be visualised using images captured by various techniques.
![]() | ||
Fig. 1 2D and 3D pore architectures in biomaterials. All images have been reproduced with permission. (A) Left: How a pore is defined in a 2D image.2 Right: Graphical representation of pore size calculation and slices used in a 3D reconstruction.45 (B) Left: 3D visualisation of pores detected as spherical bubbles: orange depicts segmented collagen fibers and blue represents pores. Right: 2D visualisation (circles) of detected 3D bubbles in the left image, in an exemplary 2D image slice. Black depicts segmented collagen fibers, blue represents determined pores of a single analysis process, orange represents detected pores of a second residual analysis.16 (C) Top left corner and centre: scanning electron microscopy (SEM) images of salt leached and gas foamed scaffolds.3 Top right centre: Scanning electron micrograph of the microfluidic foaming poly(vinyl alcohol) (PVA) scaffold.45 Top right corner: field-emission microscopy (FESM) images of electrospun poly(L-lactic acid)-co-poly(ε-caprolactone)–gelatin (PLACL/Gel) nanofibers.37 Bottom left corner: SEM micrographs of a PCL scaffolds produced by bioextrusion.15 Bottom left centre: Micro-CT 2D section of a collagen-based composite scaffold.46 Bottom right centre: Confocal fluorescence microscopy image of a collagen network.47 Bottom right corner: SEM image of an agarose hydrogel.2 (D) 3D bubble method applied to an in silico fibrin gel. Left: Red spheres represent the largest 3D bubbles that can be optimally fit in the pore zones of the gel and produce their maximum filling. Right: Zoom of a single sphere touching four different fibers.17 (E) Left: Synchroton micro-CT image of a bioactive glass scaffold. Right: Pores within the scaffold identified by image analysis algorithms.48 |
Alternatively, biomaterials may have a mixture of different pore architecture and/or irregular arrangement of pores, and as such, can be categorised as a disordered structure. It should be noted that not only the pore structure but also the pore walls may be variable across a disordered structure. Disordered structures such as disordered scaffolds,8,16,18,19,37 hydrogels2,7,13,38–40 and fibrous networks8,17,24,25,39,41–44 are examples of this and are comprised of inherently irregular porous architecture, making automated characterisation challenging. Automated methods, such as micro-computed tomography (micro-CT), 3D image reconstruction and AI-based segmentation models, have been used to analyse these architectures.10,26
This perspective primarily focuses on disordered porous materials, such as hydrogels and scaffolds, where variability in pore shape, size, and connectivity complicates standard pore characterisation techniques.38 However, understanding how automated methods have been applied to ordered porous structures, such as mesoporous silica and MOFs, whether successfully or with limitations, facilitates their adaptation for more complex materials.1,16,29,35
A fundamental challenge of 2D pore analysis is its inability to capture the true three-dimensional pore architecture. For instance, scanning electron microscopy (SEM) images (literature examples provided in Fig. 1C, 2F and 3C) captured using a high resolution technique have inherent limitations; including restricted planes of view due to fractured surfaces examined by SEM, minor pores being unclear or easily overlooked, challenges when imaging multilayered biomaterials, and lack of standardised protocols for improved robustness and limiting user-bias during pore quantification.19,43,46,47 Circular or elliptical pore cross-sections in 2D may misrepresent actual 3D pore connectivity, requiring stereological or computational corrections.17,19,20,50 Moreover, orientation dependence and sampling bias introduce further inaccuracies, particularly in fibrous or interconnected scaffold structures.17,46 Pore shape is crucial, because unless specifically engineered to achieve a highly ordered shape, pores are irregularly shaped (Fig. 1), thus physical or virtual sectioning anywhere other than where the pore diameter is maximum (Fig. 1A) leads to an underestimation of pore size.2,28,46 This is also influenced by the threshold chosen for the sample and imaging artefacts/tilts/curvatures, requiring filtering and segmentation to improve the accuracy of pore boundary determination.50 While statistical sampling methods, such as unbiased stereology are used to approximate 3D characteristics, they may not be suitable for highly interconnected porous networks.17
![]() | ||
Fig. 2 Examples of segmentation methods. All images have been reproduced with permission. (A) 3D collagen network binarisation method. Left: Representative image cube of a 5(6)-carboxytetramethylrhodamine N-succinimidylester (TAMRA)-labelled collagen scaffold. Right: A representative 2D slice of original image data, Otsu's thresholding, adaptive local thresholding of a denoised and the final segmentation result.16 (B) An unthresholded, undespeckled section of a micro-CT image of a collagen scaffold, demonstrating the noise which can be present in this type of data. Top right: Shows collagen in white, with the noise manually removed. Black represents void. Bottom right: Shows the isolated noise, seen as light grey specks.57 (C): Overview of how the DiameterJ algorithm analyses fiber diameter and other scaffold properties.56 (D) Top left: CLSM images of nanomicrofiber scaffolds. Top right: Filtered with median 3D. Bottom left: Binarised. Bottom right: 3D reconstructed with Avizo Fire software.43 (E) Top: Binarised image of a PVA scaffold. Bottom: Separation of pores using the watershed algorithm.45 (F) Effect of pixel intensity threshold, K, on pore identification and resulting effect on pore size of scaffolds.3 |
![]() | ||
Fig. 3 Examples of automated pore size determination methods. All images have been reproduced with permission. (A) Quantification of 3D pore networks of sol–gel derived bioactive glass foams from micro-CT data: (a) 2D slice of raw data, (b) 2D slice showing application of the dilation algorithm, (c) 2D representation of pores derived from the watershed algorithm, (d) 3D image of identified pores, (e) 3D image of the interconnects obtained from the top down algorithm, and (f) demonstration of the bounding box method of measuring the interconnect length.19 (B) Illustration of pore size determination of alginate hydrogel microcapsules using TEM images. Image A depicts a hydrogel as it is typically observed using transmission electron microscopy. Image B shows the results of the image segmentation after binarisation. Image C shows the result of a Euclidean distance transformation. Image D gives an overlay of the pore region image skeleton (red lines) with the original image. Image skeletons are one-pixel wide center axes. They are defined via the set of inner pore pixels. The set is defined via local distance maxima with respect to alginate segments.38 (C) Fiber network diameter determination of scaffolds (opposite of pore size determination, using similar algorithms). Top left: Starting SEM image. Top centre: Image histogram equalsation followed by 3 by 3 median filtering. Top right: Local thresholding through Otsu method. Middle left: Thinning, smoothing and removal of isolated pixel areas through a cascade of different morphological operators. Middle center: Skeletonization. Middle right–bottom left: Binary filters for Delaunay network refinement. Bottom center: Modified Delaunay network associated to the real fiber network. Bottom right: Final network and fiber diameters detected.25 (D) Simplified bubble analysis of the pore space of a random biopolymer network in two steps. Left: First, the Euclidean distance map (EDM) of the fluid space of the network structure is computed (shades of gray (colours online) indicate the distance of each fluid pixel to the nearest fiber pixel). (black) Fibers. Second, the local maxima of the EDM (white crosses) determine the centres of all 2D bubbles. To avoid bubbles of similar size in close proximity, the EDM was smoothed with a 5 × 5 Gaussian kernel with a sigma of one before the local maxima were determined. (inset) A local maximum of the EDM (red circle) is a pixel whose eight neighbours all have smaller values. Right: Resulting 2D bubbles (black circles) fit into the pore zones of the fiber structure (black). (red crosses) Bubble centers.102 (E) Comparison of detected pores with cutoff values of 85 vs. 120 in chitosan-gelatin cryogels using PoreVision software. Red outlines are pores outside the analysis boundary, blue outlines are pores removed for being too small (most likely dust, cracks, or folds), and green outlines are identified pores.20 (F) Illustration of varying results provided by micro-CT 2D pore size analysis of collagen-based scaffolds. Pores (in gray) of 3 differing shapes (left, centre, right) were evaluated by means of 4 micro-CT 2D parameters (MT – mean thickness, MD – major diameter, BICD – biggest inner circle diameter, AECD – area-equivalent circle diameter) and their values are presented in panels below the images (in mm). The results tend to differ with increasing shape irregularity.46 (G) Colour map and pixel intensity based measurement of macro pore size distribution in a polymeric scaffold.21 |
3D imaging techniques, such as micro-CT, confocal laser scanning microscopy (CLSM), transmission electron microscopy (TEM) tomography reconstruction and focused ion beam SEM (FIB-SEM), have been used for obtaining non-destructive, high-resolution information regarding pore morphology28,43,47 (Fig. 1B, D, E, 2A, B, D and 3A, B). The convoluted data and naturally complex biological networks make manual data extraction from confocal images challenging, driving automation.17 Furthermore, a resolution and volume interdependence for TEM tomography and FIB-SEM have been reported.51 Micro-CT (Fig. 1E, 2B and 3A), in particular, enables quantitative assessment of porosity and connectivity, linking these parameters to material function.20,47 Further, the technique allows virtual sectioning, visualisation in different planes and colour coded visualisation for improved comprehension.46 However, the accuracy of 3D segmentation relies on grayscale differentiation and thresholding, which are often subjective and affected by environmental factors such as lighting and computer monitor resolution and personal factors such as fatigue.20,27 3D methods often acquire repeated ‘slices’ of a 3D material, process them in a 2D format and aggregate to reconstruct the 3D image. Furthermore, the high cost and low accessibility of micro-CT, especially for nanometre level visualisation, computational intensity and data scarcity for machine learning-based 3D segmentation further complicate widespread adoption.17,20,38,43,52 As noted by Mickel et al. ‘a generic definition (for pore size) and a robust method to extract pore sizes from experimental three-dimensional microscopy data sets have been lacking’.53
Regardless of imaging method, 2D or 3D, quality of the source image is a key determinant of the resultant pore size metrics, and in manual image analysis methods in particular, there is a trade-off between accuracy and time and effort.47,52 Studies directly comparing 2D vs. 3D pore analysis have shown systematic discrepancies (% discrepancy) in extracted values, with 2D methods often underestimating porosity and connectivity.38,54 Given these challenges, hybrid approaches, combining stereological corrections with machine learning-assisted 3D reconstruction, are emerging as promising solutions.52
As we focus on direct 2D and 3D imaging techniques, the accuracy of pore quantification relies on automated image processing techniques that extract meaningful information from images. This section discusses the segmentation process (to isolate pores from the surrounding material) in Tables 1, 2 and Box 1, and pore analysis techniques (to quantify pore characteristics) in Box 2 and Table 3. Given the complexity of porous materials, different approaches, ranging from thresholding-based methods to machine learning-driven segmentation,10 have been developed to improve accuracy and reproducibility in pore characterisation.1 It is important to note that the suitability of algorithms will vary as a function of the considered biomaterial and imaging modality. Depending on available equipment, signal-to-noise ratio and contrast, separation of background and specimen and image resolution will vary, changing the task difficulty, accuracy and consequentially the choice of segmentation and analysis approach.47,54,55
Pre-processing modality | Purpose | Software for implementation |
---|---|---|
Edge detection | Allows the separation of touching objects without specifying a subjective threshold.52 | ImageJ/MatLab/OpenCV |
Contrast maximisation | Maximises the contrast between pores and walls or pores and the background to ease segmentation and analysis.21 | ImageJ/MatLab/OpenCV |
Despeckling | Despeckling sweep removes all but the largest objects in a space. Despeckling sieve removes objects below a certain area or volume. Adaptable for 2D or 3D space.57 | ImageJ |
Morphological operations | Operations that can analyse and modify image shapes and structures. Can also be implemented for noise removal and feature extraction.18,19,21,24,45,52,54 | Image J/MatLab/Quanfima python package24 |
Gaussian blur | A technique to smooth an image by averaging pixel values within a Gaussian window, effectively reducing noise.58,59 | ImageJ/MatLab/OpenCV/Scikit – image package for python |
Anisotropic diffusion filtering | Reduces noise while preserving edges by allowing diffusion to occur primarily in homogenous regions.60 | ImageJ/MatLab/MedPy Python Package |
Mean/median filtering | Reduces noise by replacing each pixel's value with the mean or median value of the surrounding pixels.12,19,25 | ImageJ/MatLab/SciPy Python Package/OpenCV |
Software | Built-in segmentation |
---|---|
https://www.imagemet.com/products_/spip/features/particle_analysis_in_images/82–84 | Intensity-based thresholding. |
https://porespy.org/index.html85–88 | Watershed segmentation, deep-learning for pixel-wise classification. |
https://examples.itk.org/80,89 | Otsu and binary thresholding, watershed segmentation. |
https://www.thermofisher.com/in/en/home/electron-microscopy/products/software-em-3d-vis/porometric-software.html90 | Gradient-based and marker-controlled watershed methods for segmentation refinement. |
https://www.manula.com/manuals/mipar/user-manual/latest/en/topic/getting-started42,91 | Multiple threshold, watershed and deep learning-based methods. |
https://www.wolfram.com/?source=nav38,92,93 | Adaptive and histogram threshold, various watershed analyses. |
ImageJ based software | Global and local thresholding, and watershed segmentation. Compatible with segmentation plugins including BoneJ and MorphoLibJ. |
Software | Built-in analysis |
---|---|
https://www.imagemet.com/products_/spip/features/particle_analysis_in_images/82–84 | (1) Pore size distribution |
https://www.manula.com/manuals/mipar/user-manual/latest/en/topic/getting-started42,91 | (2) Interconnectivity, and transport pathways |
https://gades-solutions.com/project/sigmascan-pro/96,104 | (3) Quantification using geometric fitting methods |
https://www.materialise.com/en/healthcare/mimics/ai-enabled-segmentation65 | (4) Skeletonisation |
https://www.bruker.com/en/products-and-solutions/preclinical-imaging/micro-ct/3d-suite-software.html46,105,106 | (5) Voxel-based analysis |
https://quanfima.readthedocs.io/en/latest/24 | |
https://mediacy.com/image-pro/107–110 | |
https://porespy.org/index.html79–82 | Enhances pore network characterisation with skeletonisation-based metrics. |
https://www.math2market.com/index.html44 | Covering radius transform (CRT) and micro-CT-based pore size distribution (PSD) reconstruction. |
https://www.thermofisher.com/in/en/home/electron-microscopy/products/software-em-3d-vis/amira-software/cell-biology.html?cid=msd_vds_ls_none_amr_123456_gl_pso_gaw_escchz&gad_source=119,111 | Skeletonisation. |
https://imagej.net/plugins/diameterj24,41,56 | Measure porosity of fibre networks. |
https://github.com/sedmorteza/ND112 | Measure porosity of scaffolds. |
https://imagej.net/plugins/bonej99 | Measure porosity of bone structures. |
https://www.mdpi.com/2310-2861/11/2/132#:~:text=Two%2Ddimensional%20image%2Dbased%20analyses,obtained%20from%20commonly%20available%20microscopy20 | Pore size analysis including measurement, distribution and range. Morphological analysis. |
In Box 1, we outline a range of (semi-)automated segmentation approaches. Thresholding is the simplest and most widely used segmentation technique, where pixels are classified as either pore or material, based on intensity values (Fig. 2). Thresholding-based segmentation is computationally efficient but may struggle with detecting intricate pore networks, particularly in fibrous or highly porous materials.43 Additionally, it may yield variable results due to user-defined parameters and a lack of benchmarking.20,27 Machine learning (ML) techniques, both shallow learning (e.g., support vector machines, decision trees) and deep learning (e.g., convolutional neural networks (CNN), U-Net models), can be used for pixel-wise classification in pore segmentation.5,20,39,61–63 However, machine learning is data-hungry, thus requiring larger sample sizes and manual ‘ground-truth’ labelling which is labour intensive and may introduce user-bias into the segmentation process.
Box 1: Segmentation approaches![]() • Global thresholding: a single threshold is applied across the entire image, making it suitable for materials with uniform contrast but less effective for heterogeneous structures.1,3,8,21,38,43,47,52,56,64–67 These techniques are also more sensitive to noise.68 Some examples of global thresholding methods are given below, more methods and an in-depth discussion on thresholding could be found through Rajagopalan et al.27 ○ Otsu's global thresholding: selects an optimal segmentation threshold by maximising inter-class variance. This method is effective for materials with bimodal intensity distributions, where pores and solid structures have distinct grayscale values.19,67,69 ○ Entropy-based thresholding: the foreground and background are treated as separate sources and the optimal threshold is chosen at the maximum of the sum of the two class entropies. This is useful for images with complex histograms.24,27,45 ○ Histogram thresholding: analysing the concavity points of the image's histogram's convex hull. Any valley may be considered as a potential threshold, with the deepest concavities being favoured.27,70 • Local (adaptive) thresholding: the threshold varies across different regions, allowing segmentation in images, in which contrast changes due to uneven illumination or varying material densities. This is useful for heterogenous structures.16,20,25,54,71 Several local thresholding techniques exist as outlined by Rajagopalan et al.27 • Contrast enhancement with Fourier transform: Fourier-based contrast enhancement has been integrated to improve edge detection in biomaterial samples, particularly in scaffold imaging.5 • Manual thresholding in ImageJ: adjusting the threshold manually to match original image features, as seen in some studies.1,72–76 This technique can introduce various problems associated with room lighting, fatigue of the operator and limited grey-scale shade perception.27 • MATLAB-optimised thresholding: ○ Region-growing and edge-detection algorithms: allows for adaptive thresholding, enhancing pore segmentation in fibrous and highly interconnected structures.52,77 ○ Gradient filters: improve thresholding accuracy by refining intensity distributions, minimising over-segmentation errors78,79 • Watershed-based segmentation: a region-based segmentation technique that treats an image as a topographic surface, where intensity variations correspond to elevations. It is useful for segmenting complex pore structures in porous materials by identifying watershed lines that separate adjacent regions. Following initial segmentation, the watershed algorithm refines pore boundaries by treating intensity gradients as a height map. By simulating the flooding of an image, it effectively separates adjacent pores that may have been grouped together in thresholding-based segmentation. While this method enhances segmentation accuracy, it may require post-processing to address over-segmentation.18,19,43,45,48,66,80 MATLAB-based implementations of the watershed algorithm further improve pore segmentation accuracy through gradient-based seed point refinement.43 ○ Gradient-based watershed: uses the gradient magnitude image to detect high-intensity ridges that define the segmented boundaries.24 ○ Marker-controlled watershed: introduces predefined seed points to prevent over-segmentation, which is beneficial for fibrous or highly interconnected porous materials.20,43,66 ○ 3D watershed for pore connectivity: applied in voxel-based 3D reconstructions to quantify pore interconnectivity in porous scaffolds.18,19,43,48,66,80,81
• Shallow ML models: shallow learners struggle to classify raw image input and depend on efficient preprocessing and predefined, representative characteristics such as intensity and edges.27 • Deep learning models: deep learning models have demonstrated enhanced performance in segmenting irregular and disordered pore architectures, particularly in hydrogels and scaffolds.10,62,63 CNNs and U-Net architectures learn hierarchical representations from large datasets, allowing robust segmentation of complex pore structures.5,20 These model are capable to automatic feature selection and extraction, enhancing contrast in low-resolution porous material and removing the need for manual tuning.62 Furthermore, advanced machine learning techniques improve segmentation across multiple imaging modalities, making them adaptable to various pore architectures.61 One problem with deep learning models is their data hungry nature, requiring large, annotated datasets for successful training and prevention of overfitting. ○ CNN-based approaches: the convolutional nature of CNN's enables multi-scale feature extraction, which improves segmentation accuracy, enhances pore boundary detection and reduces errors in low-contrast regions.63 The hierarchical nature of CNNs refines pore morphology, effectively distinguishing pores from material phases. Pretrained CNN models, such as VGG16 and ResNet, further accelerate adaptation by reducing the need for large, manually annotated training sets, while transfer learning enhances model generalisability across different biomaterial types, ensuring consistent segmentation accuracy.61,62 ○ U-Net architecture approaches: provide pixel-wise segmentation, making them particularly effective for highly interconnected porous biomaterial. Their encoder–decoder structure allows for detailed feature extraction, while skip connections preserve fine-grained pore structures, significantly improving segmentation accuracy. Furthermore, data augmentation techniques enhance model robustness, compensating for limited training datasets and increasing adaptability across various biomaterial imaging conditions5,61 |
Software given in Table 2 include modalities for different segmentation techniques. It is important to note that many studies utilise in-house, tailor-made code to conduct segmentation, that best befits their requirements.
Both threshold and machine learning-based approaches offer unique advantages to image segmentation with neither demonstrating consistently better performance albeit machine learning offers faster computation time. The choice of segmentation technique directly influences measures of porosity and pore structures, and the careful validation of image-based pore data remains crucial.45,64,94
Box 2: Pore analysis approaches![]() • Pore-to-material ratio: the fraction of pore area relative to material area, expressed as a decimal or percentage. This is calculated using tools such as ImageJ or FiJi, which is useful for assessing porosity, however, lacks structural context.24,64,67,77,95,96 • Pore count: the total number of detected pores; limited when pores are irregular or overlapping. Suited for homogenous structures.5,21 • Pore area distribution: captures variability in pore sizes; however, segmentation errors can skew results.6,21,45,59,67,77
Pixel-based approaches • Exponential decay fit to pore pixel spacing distribution: quantifies the spatial arrangement of pores by fitting an exponential decay model to the pore-pixel spacing distribution.27 This approach helps distinguish between homogeneous (evenly distributed pores) and heterogeneous (randomly distributed or clustered pores) porous networks. • Fourier transform analysis: used to detect periodicity in pore spacing, Fourier analysis can reveal structural anisotropy within porous scaffolds and hydrogels.36,71,77,97 • Skeletonisation-based analysis: by reducing pore structures to their one-dimensional (1D) medial axis, this approach allows for measuring pore branching, interconnectivity, tortuosity and transport pathways in porous biomaterials scaffold.19,38,56,98 Geometric transform-based methods • Maximum covering radius transform: the CRT method determines the largest inscribed circle within each pore, quantifying local pore size variations across a sample.46,55,66,99,100 • Morphological opening and closing transformations: these operations refine pore boundaries by removing small artifacts and enhancing true pore structures, aiding in more accurate segmentation.1,16,18,19,47,54,65,66 For instance, connected pores can be separated using morphological dilation followed by erosion—dilation shrinks the pores by expanding scaffold boundaries, and erosion then restores pore size while maintaining separation.1 • Voxel-based pore size distribution reconstruction: 3D voxel reconstruction has been used to quantify PSDs, for improved interconnectivity analysis.4,38,56,66,67 • Fourier-based feature extraction for pore shape: fourier-based shape analysis has been integrated into voxel-based reconstructions to improve pore morphology quantification.5 Pore fitting methods • Sphere/circle fitting: pore fitting methods are based on fitting a circle or a sphere to a pore, based on criteria such major diameter (MD, major diameter of analysed pore), mean thickness (MT, based on circle-fitting algorithm similar to sphere-fitting method), biggest inner circle diameter (BICD, diameter of the biggest circle that fits the pore) and area-equivalent circle diameter (AECD, diameter of the circle with an equivalent area to that of the pore).12,45,46,48,65–67,97,99,101 Ideally well suited for circular pores.1 • Bubble analysis: a geometric method that mimics fluid dynamics in fibrous networks and is based on the largest possible circle inside a pore that touches three surrounding fibers,16,17,46,53 however, according to Fischer et al. may not fully account for the residual fluid volume.16 In an article-response dynamic, Molteni et al. proposed an algorithm that randomly seeds a pore and expands its boundary until it tangentially contacts three fibers,17 for which Münster & Fabry developed a simplified approach using Euclidean distance mapping to detect local maxima, identifying the largest pore regions.102 However, Molteni et al. argued that this method lacks filtering and requires additional post-processing to remove overlapping pores and those that do not satisfy the three-fiber contact criterion.103 Applicable to fibrous networks such as bundled F-actin, fibrin, cytoskeletal filament networks, given individual fibers can be resolved.53 • Destroy and rebuild method for PSD analysis: This micro-CT-based technique reconstructs 3D pore networks from 2D slice data, offering an alternative to direct segmentation-based methods. ImageJ's particle analyser (PA) function extracts mass centre coordinates (Xc, Yc) and pore section areas for each binarised slice. An algorithm groups sections by spatial proximity, assuming pores maintain a spherical shape. Edge-affected pores are removed, and equivalent volume and radius are calculated for each pore to determine the PSD.5,45,56 In addition, MATLAB-based PSD reconstruction methods have been applied to improve accuracy in volume estimations of irregularly shaped pores.77 |
Imaging quality—including resolution, magnification, and field of view—strongly influences the accuracy of segmentation and measurement methods. Poor image quality or insufficient resolution often leads to inconsistent or inaccurate pore size detection, particularly in fibrous or highly irregular networks.41,47
Circle-fitting algorithms can be useful for regular, isotropic materials with well-defined, circular pores, but often over- or underestimate pore sizes in more complex biomaterials. Many image analysis tools still depend on manual thresholding, subjective filtering, or operator bias, particularly when the software allows for interactive selection or exclusion of pores during post-processing.20 While software such as DiameterJ and PoreSpy offer user-friendly platforms, their results are sensitive to image parameters, such as pixel size and contrast depth, and sometimes require subjective validation.41
Importantly, there is no universal standard or pipeline for PSD analysis that suits all materials. Algorithms often need to be tailored to specific scaffold morphologies. Nevertheless, some methods demonstrate cross-material applicability, successfully analysing electrospun meshes, hydrogels, and decellularised tissues with the same underlying algorithm.25
A growing concern in the literature is the lack of detailed reporting on how PSD is computed. Some studies present PSD values without any or very limited mention of segmentation methods, threshold settings, or image pre-processing steps, and many studies cite ‘in-house algorithms’ without sufficient detail or provided code.18,20,24,41 This hampers reproducibility and impedes method comparison while inconsistency in terminology increases the entry barrier to the field. Moreover, claims of “fully automated” analysis often mask semi-automated steps that introduce bias, such as filtering out overlapping or irregular pores without objective criteria.41 Closely related to this is the inconsistency across reported metrics, with many studies failing to report objective evaluations of the PSD methodology.
Most current approaches rely on unsupervised learning methods without a reference standard or ‘ground-truth’, which makes performance evaluation difficult. Many studies rely on qualitative validation (visual matching) rather than reporting consistent quantitative metrics such as measurement error or precision. When validation is provided, it is often inconsistent or non-comparable across studies.20 This may be partially attributed to the lack of available ‘ground-truth’ values, which makes the implementation of supervised learning challenging. However, while the manual measurement of pore size is difficult and includes human error, studies should aim to validate proposed automated methods using ‘ground-truth’ measurements to ascertain the reliability of the selected methods.
Closely related is a need for reporting objective evaluation metrics of automated methods. While the lack of ‘ground-truth’ complicates the implementation of supervised learning, steps should be taken to evaluate algorithms in the absence of a ‘true’ value. While some studies have implemented algorithm validation using standard reference images,1 others compare results between methods to establish mean performance and improve reliability of findings. While several supervised and unsupervised evaluation methods thus exist, many studies continue to accept pore size estimation at face value without considering the associated bias of the measurement tool. Future studies should aim to conduct critical evaluation of their estimation techniques and in the absence of an objective ‘ground-truth’ should at the very least compare the performances of different algorithms.
There is also limited ability among automated algorithms to handle irregularly shaped or anisotropic pores. Many algorithms interpret elongated or eccentric pores as two or more circular pores, which misrepresents the actual pore size distribution.35 While some algorithms have introduced shape classification features, their accuracy for small or faint shapes remains limited.
Looking ahead, a promising direction is the integration of machine learning and deep learning approaches trained on annotated datasets to improve pore recognition and classification accuracy. This includes the potential for algorithms that are not only shape-aware but also context-aware—able to distinguish between artefacts and true pores based on their position, orientation, or relationship with surrounding structures.
We recommend the following:
• Standardised reporting: all studies should include detailed information on segmentation algorithms, image processing steps, and evaluation metrics used in PSD analysis.
• Material-specific validation: algorithms should be validated for specific scaffold types and pore morphologies, not assumed to be generalisable.
• Open-source development: wider adoption of open-source platforms would promote reproducibility, comparability, and collaborative development.
• Integration of AI and ML: supervised machine learning and AI-enhanced approaches could overcome current limitations in pore shape recognition, classification, and error quantification.
• Consistency in evaluation: future studies should report standard error metrics and comparison benchmarks to assess accuracy, ideally using both synthetic and real-world datasets.
Ultimately, while automated PSD analysis is advancing rapidly, further work is needed to balance ease of use, reproducibility, and the ability to capture complex pore features across diverse biomaterials, including dynamic hydrogels and 4D scaffolds.
This journal is © The Royal Society of Chemistry 2025 |