
==== Front
BMC Med Imaging
BMC Med Imaging
BMC Medical Imaging
1471-2342
BioMed Central London

39300334
1423
10.1186/s12880-024-01423-0
Research
Three-dimensional numerical schemes for the segmentation of the psoas muscle in X-ray computed tomography images
Paolucci Giulio 1
Cama Isabella 1
Campi Cristina 12
Piana Michele piana@dima.unige.it

12
1 https://ror.org/0107c5v14 grid.5606.5 0000 0001 2151 3065 MIDA, Dipartimento di Matematica, Università di Genova, via Dodecaneso 35, Genova, 16145 Italy
2 https://ror.org/04d7es448 grid.410345.7 0000 0004 1756 7871 IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi 10, Genova, 16132 Italy
19 9 2024
19 9 2024
2024
24 2516 1 2024
6 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
The analysis of the psoas muscle in morphological and functional imaging has proved to be an accurate approach to assess sarcopenia, i.e. a systemic loss of skeletal muscle mass and function that may be correlated to multifactorial etiological aspects. The inclusion of sarcopenia assessment into a radiological workflow would need the implementation of computational pipelines for image processing that guarantee segmentation reliability and a significant degree of automation. The present study utilizes three-dimensional numerical schemes for psoas segmentation in low-dose X-ray computed tomography images. Specifically, here we focused on the level set methodology and compared the performances of two standard approaches, a classical evolution model and a three-dimension geodesic model, with the performances of an original first-order modification of this latter one. The results of this analysis show that these gradient-based schemes guarantee reliability with respect to manual segmentation and that the first-order scheme requires a computational burden that is significantly smaller than the one needed by the second-order approach.

Keywords

Image segmentation
X-ray Computed Tomography (CT)
Three-dimensional level set methods
Sarcopenia
PROGRAMMA OPERATIVO REGIONE LIGURIA FONDO SOCIALE EUROPEO 2014-2020RLOF18ASSRIC/58/1 Hub Life Science – Digital Health (LSH-DH) PNC-E3-2022-23683267 - Progetto DHEAL-COMD33C22001980001 D33C22001980001 issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

Sarcopenia [1] is a pathological condition of the skeletal muscle characterized by loss of mass, and depletion of strength and physical performance. In the last decade, several studies have shown that sarcopenia is associated to specific chronic diseases. Just as examples, and not exhaustively, low skeletal muscle mass is associated with post-transplantation death [2, 3]; it is a risk factor for poor performances and toxicity of chemotherapy intervention in cancer patients [4–8]; it has a high probability of co-occurrence with several neurodegenerative conditions such as Parkinson’s disease and multiple sclerosis [9, 10].

Despite its relevance as a frailty biomarker, sarcopenia is difficult to quantitatively determine, most studies in this field currently relying on estimates of variation of the mass of the psoas muscle determined by post-processing specific slices of low resolution X-ray computed tomography (CT) scans [11]. Specifically, in a typical pipeline utilizing clinical data, abdominal CT images of the patient are downloaded from the hospital Picture Archiving and Communication System (PACS); a cross-sectional image in correspondence with the third lumbar vertebra (L3) is selected; the psoas is identified in this image; and the corresponding psoas area and density are computed.

Different approaches to automated two-dimensional psoas segmentation have been proposed. For example, anatomy-based manually extracted prior shapes have been utilized to initialize energy-based iterative schemes [12–14]. More recently, artificial intelligence (AI) relied on convolutional neural networks (CNNs) and Generative Adversarial Networks (GANs) to realize CT images segmentation specifically tailored to psoas identification at L3 [15, 16]. Of course, deep learning algorithms generally designed for muscles and internal organs segmentation [17, 18] can be applied to the specific task of psoas identification.

This image-based approach to sarcopenia quantification relies on the controversial assumption that a single CT slice is representative of the whole psoas condition [19–21]. Very recently [20, 22], two studies conceived in our (extended) research group accounted for this ambiguity and showed that a complete three-dimensional segmentation of the psoas muscle in CT images of oncological and neurological patients may be exploited to obtain metabolic information about the diseased tissue. Both studies utilized hybrid data made of low-dose CT and positron emission tomography (PET) data of patients injected with 18F-Fluorodeoxyglucose (FDG-PET/CT), and, in both cases, the computational paradigm was based on the following three steps [23–25]: A segmentation algorithm was applied to co-registered CT data in order to automatically identify the psoas muscle.

A binary mask was constructed, where each voxel within the image region corresponding to the muscle was given value 1 and all other voxels were given value 0.

The binary mask was voxel-wise multiplied with the FDG-PET volumes in order to extract the metabolic information associated to the psoas tissue.

Using this scheme, paper [20] introduced a novel image-driven biomarker, the Attenuation Metabolic Index (AMI), that can be applied as a prognostic tool in metastatic castration-resistant prostate cancer patients, while paper [22] proved that there is a common underlying mechanism for skeletal muscle and spinal cord hypermetabolism in amyotrophic lateral sclerosis (ALS).

Of course, the reliability of these metabolic outcomes strongly depends on the accuracy with which the segmentation step 1 above is realized. Specifically, in the first study [20] the segmentation approach was based on an extended version of the Hough transform [26], which can be used only if a realistic parametric model for the muscle is at disposal; the second study [22] utilized a rather heuristic approach based on histogram equalization and an α-shape algorithm to identify the region corresponding to the inner muscle. However, despite the increasing potential relevance of this kind of analysis for diagnostic/prognostic applications, no systematic study has been performed so far on the effectiveness of three dimensional numerical schemes for image segmentation of the psoas muscle in low-resolution CT data.

The objective of the present paper is to address the computational problem of segmenting the psoas muscle in clinical low-dose CT data by means of three dimensional numerical schemes originating from the level set methodology [27–31]. The use of three-dimensional methods is particularly appropriate for assessing sarcopenia. Indeed, volumetric segmentation naturally exploits the information coming from adjacent regions, improving the accuracy with which the tissue density is computed [32, 33]. Specifically, we focused on gradient-based methods that search for an optimal solution via a search orientation based on the minimization of an objective function [34]. These methods can incorporate in their objective function weighted curve length and variation of intensity in the segmented regions [35], and can be combined with deep learning techniques [36]. In the present study we considered a classical evolution model, a standard second-order geodesic model, and a first-order modification of the geodesic model that is formulated for the first time in the present paper. In all three cases the model discretization has been performed by means of a standard finite-difference scheme. As a result, we obtained three segmentation algorithms, all with the nice property that they require the users to select just one point for each psoas in the data volume to initialize the segmentation process. The application to CT images contained in hybrid PET/CT volumes showed a notable segmentation accuracy, although with rather significant differences in the required computational burden. However, due to the low spatial resolution of these clinical CT data, in some cases the segmented shapes presented some artifacts at the end of the prescribed iterations. Therefore, we have formulated and implemented an algorithm for removing these false positives, which can be applied at the end of the evolution process when necessary.

We finally point out that three-dimensional segmentation schemes allow the leveraging of information from adjacent slices [37]. However, in the case of medical imaging applications, the use of three-dimensional deep-learning-based algorithms may be hampered by the computational burden and, sometimes, by the lack of sufficiently large data sets, necessary in the training phase of the networks [32]. Therefore, the introduction of three-dimensional deterministic numerical schemes may be particularly supportive for this kind of applications.

The plan of the paper is as follows. “Evolution models and their discretization” section describes the level set models employed in the study, together with the numerical scheme utilized for their discretization. “Pre-processing and post-processing” section discusses the more heuristic approaches exploited to pre-processing the CT images and post-processing the outcomes of the numerical evolution process. “Numerical results” section compares the results provided by the different numerical procedures when applied to clinical CT data. Our conclusions are offered in “Comments and conclusions” section.

Evolution models and their discretization

The general framework considered in the present study assumes that segmentation is realized by letting a three-dimensional wave-front1 γt:={x∈R3:v(t,x)=0}

evolve from an initial point inside the shape to segment until it reaches its boundary. Following the level set method [27], the field v=v(t,x) is assumed to be the solution of the Cauchy problem2 vt+g(x)|Dv|=0,(x,t)∈R3×(0,T)v(0,x)=v0(x)x∈R3,

where D is the three-dimensional gradient operator, g(x) is the speed in the normal direction and v0 is a proper representation of the front γ0 (i.e. denoting by Ω0 the region enclosed by the front, v0 has to be equal to 0 on γ0, positive outside Ω0 and negative inside). It is well known that, under appropriate assumptions on the velocity field g(x) and the initial condition v0, the viscosity solution of (2) exists and is unique (see, e.g., [38]). In the following we refer to the eikonal Eq. (2) as the classical model. This approach can be generalized to a second order Hamilton-Jacobi equation of the form3 vt+H(x,Dv,D2v)=0,

which includes more sophisticated definitions of the velocity field g, possibly dependent on the curvature of the surface, as for example in the Mean Curvature Motion. Here we consider the version proposed by Caselles et al. in [39] denominated geodesic model4 H(x,Dv,D2v)=-g(x)|Dv|ε∇·Dv|Dv|-μ-ηDg·Dv,

where ε,μ,η are positive parameters to appropriately tune. To give a better insight on the choices made in the numerical section below, we recall that these parameters are used to properly weigh the different components of the model. For instance, η is related to the transport driven term, which attracts the front towards the boundaries of the image, and ε is related to the curvature part that regularizes the solution and tends to shorten the front. Finally, the constant velocity μ is added to balance the shrinking effect of the curvature term and allows outward evolution, as in our case, where the front expands from an initial sphere towards the surface of the muscle. We decided to give a normalized value to the constant μ=1 and adapted the other parameters accordingly. More precisely, the higher the value of η, the slower the expansion toward the boundary, whereas high values for the curvature term enforce the use of very restrictive Courant-Friedrichs-Lewy (CFL) conditions to keep the finite difference scheme stable.

We also propose a first-order modification of Eq. (4), i.e.5 H(x,Dv)=μg(x)|Dv|-ηDg·Dv,

which is characterized by a lower computational burden. In order to complete the description of the segmentation model we have to define the velocity g(x) such that it forces the front towards the relevant boundary and then it stops its evolution. A typical definition in gradient-based methods is the following6 g(x)=11+|∇(G∗I(x))|p,p≥1,

where I represents the intensity values of the image, and G a Gaussian filter of deviation σ, used to regularize noisy images. This function has the desired properties for edge detection; in fact, it takes values in [0, 1], and it is decreasing with respect to the gradient of the image I. If a spike in the gradient is present, then g is close to 0; otherwise it is equal to 1 in smooth regions. In practice, the filtering process is done once before the evolution, and the values of the smoothed image G⋆I are used instead of I in the computation of the gradient. High values for the parameter p may enlarge the region where the function is close to 0, reducing the overall accuracy of the scheme. This is why here we used (6) with p=2 as in [39].

We now provide some details about the numerical schemes utilized for the discretization of the first order Eqs. (2) and (5), which can also be adapted to Eq. (4) with minor changes. For sake of simplicity the description of these schemes is provided in the two-dimensional setting, although the codes we used for data analysis have been implemented for three-dimensional segmentation.

We realized the numerical discretization of all three models by means of the standard difference form7 ui,jn+1=SΔ(un)i,j:=ui,jn-Δthxj,yi,Dx-ui,jn,Dx+ui,jn,Dy-ui,jn,Dy+ui,jn,

where the numerical Hamiltonian is the local Lax-Friedrichs Hamiltonian8 h(x,y,p-,p+,q-,q+):=Hx,y,p++p-2,q++q-2-αx(p-,p+)2(p+-p-)-αy(q-,q+)2(q+-q-),

with9 αx(p-,p+):=maxx,y,q,p∈I(p-,p+)Hp(x,y,p,q),αy(q-,q+):=maxx,y,p,q∈I(q-,q+)Hq(x,y,p,q),

with I(a,b):=[min(a,b),max(a,b)]. We point out that h(x,y,p-,p+,q-,q+) is a Lipschitz continuous function, with10 Dx±ui,jn:=±ui,j±1n-ui,jnΔx,

and11 Dy±ui,jn:=±ui±1,jn-ui,jnΔy.

Further, in this context consistency is imposed via condition12 h(x,a,a,b,b)=H(a,a,b),

whereas monotonicity requires the numerical Hamiltonian h to be non-decreasing with respect to its second and fourth argument and non-increasing with respect to the third and the fifth ones. The resulting scheme (7)-(8) is monotone under the Courant-Friedrichs-Lewy (CFL) type condition ΔtΔx·αx+ΔtΔy·αy≤1. Moreover, as suggested in [40], this scheme can also be used to solve Eq. (4) by first splitting the Hamiltonian into two componentsH(x,y,Dv,D2v)=H0(x,y,Dv)+H1(k(x,y),Dv),

where H0 includes the eikonal and advection term, whereas H1 accounts for the dependence on the curvature k(x, y), which is computed ask(x,y)=∇·Dv(x,y)|Dv(x,y)|=vxxvy2-2vxvyvxy+vyyvx2(vx2+vy2)3/2.

Then, it is enough to discretize the first-order part using the Lax-Friedrichs scheme and H1 with simple centered finite difference approximations. In this case, the stability of the scheme requires a parabolic-type CFL condition of the form Δt∼Δx2.

Pre-processing and post-processing

Before solving the evolution equations, we pre-processed the CT data according to the following process: Hounsfield normalization: given that normal bone tissue in CT has a Hounsfield unit (HU) of around 160, all HU values above 120 have been set to 0, in order to reduce the relative weight of the bones [41].

For image enhancement purposes, a histogram equalization algorithm has been applied, which performed a re-mapping of the image based on the probability distribution of the input gray levels [42].

A Gaussian filter with standard deviation σ has been applied to improve the signal-to-noise ratio of each image [43].

This process was executed slice-wise in a two dimensional manner. Finally, we created a mask for all HU values inside the interval [5, 120] to save the position of the non-muscular tissue, which is needed to improve the stability of the scheme for practical application.

The initialization of the scheme requires the user to select two points, one for each muscle. These points become the centers of the two spheres of radius r=5Δx used as initial condition for the evolution, which is executed separately for each psoas. For all simulations the spatial discretization (the length of the voxel side) is set to Δx=0.1, whereas the maximum number of iterations Nmax is heuristically chosen case by case, depending on the number of slices in the CT volume and the overall quality of the images.

As previously remarked, clinical low-dose CT images as the ones acquired by PET/CT scanners, are usually noisy and present a lower level of resolution with respect to high-resolution CT acquisitions. Further, the complexity of the human body in the psoas region often implies the presence of sections where different organs or muscles, which have close HU values, have numerous contact points. Therefore, the definition of a clear boundary of the psoas becomes an intricate issue, which impacts the reliability of the numerical scheme in automatically segmenting the image. Specifically, in some cases it may happen that the wave-front leaves the correct boundary and starts detecting spurious objects outside the muscle of interest, thus making the usual stopping criterion ineffective, with unreliable final results.

In order to solve this issue, we devised a post-processing routine able to clean the result of (most of) the external objects and eventually give an acceptable final result. The procedure is composed of two steps and has a very low computational cost, which makes it efficient even if more than one iteration of the second step is needed. The routine works separately for each slice, but makes use of the information given by the closest slice in order to exploit the intrinsic continuity of the surrounding problem. We have summarized the procedure in Algorithm 1, while more details on the implementation are given in the following items.First step. The main aim of this step is to reduce the number of connected components (CCs) in each slice, detecting the one belonging to the psoas and discarding the remaining ones. To do so, the algorithm starts from the slice containing the center chosen by the user and selects the correct CC by trivially excluding the ones that do not include the chosen point. Then, the remaining slices are scanned toward the top and the bottom, separately, always starting from the center. The main loop consists in first computing the centroid of the CC detected in the previous slice, and then, in the current slice, in selecting just the CC containing such point. The final result is a segmented muscle in which only the spurious objects directly connected to the boundary may be still present.

Second step. In this second step the starting slice and the previous one, in the direction of the scanning, are assumed not to need the cleaning procedure. With this assumption, a vector tol is initialized with length equal to P, the set of segmented voxels in the starting slice, where the i-th entry is the sum of minimum distance of the voxel i∈P from the set of voxels Pprev in the previous slice with a fixed constant (0.75). Then the algorithm loops over the slice, setting Pprev=P and first computing the difference set Pdiff=P-Pprev, where P and Pprev are the sets of voxels in the current and previous slices, respectively; then, for each point j∈Pdiff, the algorithm computes the set distance dist(j,Pprev) with respect to Pprev. This distance is obtained in correspondence of a voxel imin∈Pprev. If dist(j,Pprev)>tol(imin), voxel j is excluded from P and hence removed from the segmentation. After all voxels in Pdiff have been scanned and removed if needed, a new tolerance vector tol is computed using the current cleaned set of voxels P and Pprev.

The second step, in its initialization phase, assumes that the starting slice and the previous one (according to the scanning direction) do not need any cleaning. If this is not the case, a second iteration of this step is needed, using a different starting point.

We finally point out that, on the one hand, the pre-processing procedures described at the beginning of this section, which are rather standard, can be realized by means of algorithms different than the ones employed in this study. For example, image enhancement can be obtained by means of methods based on integral transforms of AI-based techniques (see [44] for a review of these methods); and noise reduction can be obtained by means of adaptive or anisotropic Gaussian filtering [45, 46]. On the other hand, the post-processing procedure has been tailored on the three-dimensional schemes introduced in this study, although both steps address issues that are typical of image segmentation performed by means of level set methods. At a more general level, we validated the outcomes of the overall segmentation process with respect to modifications of the input parameters in the pre-processing steps, and found out a notable robustness of the outcomes. However, for a few subjects such outcomes (before post-processing) are characterized by a lower quality (see Fig. 3 below for some examples), which is due to the topographical complexity of the district close to the psoas and to the poor resolution of the low-dose CT volumes. For such cases, the second step of the post-processing phase typically requires more than one iteration. In any case, we think that a more systematic investigation of the impact of different post-processing procedures on the segmentation quality is worthwhile consideration for future studies.

Algorithm 1 Cleaning procedure

Numerical results

The pipeline made of pre-processing, numerical segmentation, and post-processing has been applied to the CT images (see Fig. 1) acquired for nine subjects belonging to a retrospective data set of patients recruited at the IRCCS Ospedale Policlinico San Martino, Genova, Italy. All subjects were submitted to FDG-PET/CT and the study was performed according to the Declaration of Helsinki, Good Clinical Practice, and local ethics regulation. All enrolled patients signed a written informed consent at the time of FDG-PET/CT, econmpassing the use of anonymized data for retrospective research purposes.

For numerical segmentation we used the three numerical schemes introduced in “Evolution models and their discretization” section, i.e., the finite difference scheme (7)-(11) applied to (2) (‘classical’ from now on), and to the two geodesic schemes (4) and (5) (‘GMFD 1ord’ and ‘GMFD 2ord’ from now on, respectively). For both the two finite difference goedesic model schemes we set μ=1 and ν=0.25, while ϵ=0.05 for the second order model. Further, ’GMFD 1ord’ exploits the linear CFL condition Δt=λΔx, with13 λ=(3‖g‖+‖Dxg‖+‖Dyg‖+‖Dzg‖)-1,

and ‘GMFD 2ord’ utilizes the parabolic CFL condition Δt=0.25Δx2.Fig. 1 Coronal view of the psoas muscles (in yellow) and the corresponding axial views at different levels of the abdomen

Table 1 contains the number of iterations and CPU time employed to realize the segmentation by the three numerical schemes when applied to the volumes of all nine patients. Figure 2 contains the results of the segmentation provided by the three methods in the case of six subjects and Fig. 3 describes the performances of the post-processing algorithm in the case of the remaining three subjects. Specifically, the three cases in this latter figure are the ones for which the outcomes of the numerical schemes present a particularly significant amount of artifacts. However, the post-processing step can be applied to all outcomes, including the ones described in Fig. 2. The segmentations in Figs. 2 and 3 show that the run of the numerical schemes without any post-processing may imply either the persistence of some artifacts that typically show up in the region between the two psoas muscles, or a slight shrinking of the psoas volume. The presence of false positive or shrinking may impact the assessment of sarcopenia realized by means of the PET/CT-based process described in the Introduction. In fact, an inaccurate evaluation of the psoas dimension implies the construction of inaccurate binary maps that are used to extract metabolic information from the PET data, and may result in an over- or under-estimated determination of the FDG uptake of the psoas tissue. Table 1 Number of iterations, CPU time, values of p and σ, and number of slices for the three segmentation algorithms for the nine CT volumes considered in this study

Method	Patient	Nmax	CPU	p	σ	Slices	
classical	CTTO1	400	178.48	2	1	49	
GMFD 1ord	CTTO1	400	197.59	2	1	49	
GMFD 2ord	CTTO1	1900	1955.84	2	1	49	
classical	CTTO02	350	133.19	2	0	45	
GMFD 1ord	CTTO02	350	147.56	2	0	45	
GMFD 2ord	CTTO02	2000	1870.16	2	0	45	
classical	CTTO03	450	232.91	2	1	57	
GMFD 1ord	CTTO03	450	250.65	2	1	57	
GMFD 2ord	CTTO03	2700	3298.88	2	1	57	
classical	PZGE01	700	229.96	2	2	33	
GMFD 1ord	PZGE01	700	234.94	2	2	33	
GMFD 2ord	PZGE01	3500	2654.22	2	2	33	
classical	PZGE02	700	232.67	2	2	35	
GMFD 1ord	PZGE02	700	246.20	2	2	35	
GMFD 2ord	PZGE02	3500	2787.62	2	2	35	
classical	PZGE03	700	207.53	2	2	31	
GMFD 1ord	PZGE03	700	224.22	2	2	31	
GMFD 2ord	PZGE03	3500	2499.52	2	2	31	
classical	PZTO01	500	258.00	2	0	57	
GMFD 1ord	PZTO01	500	275.68	2	0	57	
GMFD 2ord	PZTO01	2500	2991.82	2	0	57	
classical	PZTO03	500	229.58	2	0	52	
GMFD 1ord	PZTO03	500	250.62	2	0	52	
GMFD 2ord	PZTO03	2500	2758.94	2	0	52	
classical	PZTO04	450	207.79	2	1	51	
GMFD 1ord	PZTO04	450	230.97	2	1	51	
GMFD 2ord	PZTO04	2300	2540.63	2	1	51	

We compared the segmentation results we obtained with the three methods described in the previous sections with the ones provided by manually drawn profiles within the OsiriX software package [47] and by the 3D Slicer extension TotalSegmentator (TS) [48], a fully automated total body segmentation deep learning routine based on nnU-Net [49], which has been applied with reliable results for the segmentation of abdominal tissues even in the three-dimensional setting [50–52]. In our analysis, TS has been used in both fast (TS-FM) and full-resolution (TS-FRM) modes. In order to evaluate the agreement between the manually drawn psoas profiles and the segmentations obtained with the three numerical schemes and the TotalSegmentator tool we used six metrics: the Dice similarity coefficient [53], the Jaccard index [54], the Hausdorff distance [55], the Average Symmetric Surface Distance (ASSD) [56], the Surface Dice (Surface Dice) and the Robust Hausdorff distance (Robust Hausdorff). Surface Dice and Robust Hausdorff have been computed via the packages provided by https://github.com/google-deepmind/surface-distance/tree/master.

The Dice similarity coefficient and the Jaccard index vary between 0 and 1 where 0 represents an empty intersection between the two segmentations while 1 a perfect overlap. The Hausdorff distance and ASSD are indeed distances, so smaller values correspond to closer objects. Specifically, the Hausdorff distance is the greatest distance from a point in one segmentation mask to the closest point in the other mask, while ASSD measures the average of all distances from points on the boundary of the first mask to the boundary of the second one, and viceversa. Similarly to the classical volumetric Dice, Surface Dice measures the overlap of two surfaces, where a surface element is counted as overlapping when the closest distance to the other surface is less than or equal to the specified tolerance, fixed at 3 mm. The range of the Surface Dice is [0,1]. Robust Hausdorff is computed like the classical Hausdorff distance but it uses a fixed percentile of the distances instead of the maximum distance. We choose the 75th percentile. In Table 2 we show the average values of these metrics computed across the nine subjects, together with the corresponding standard deviations.Fig. 2 Segmentations provided by the classical evolution model (2) (‘classical’, first column), the first-order geodesic model (5) (‘GMFD 1ord’, second column), and the second order geodesic model (4) (‘GMFD 2ord’, third column). First row: subject CTTO02 (45 slices) with p=2 and σ=0. Second row: subject PZGE01 (33 slices) with p=2 and σ=2. Third row: subject PZGE02 (35 slices) with p=2 and σ=2. Fourth row: subject PZGE03 (31 slices) with p=2 and σ=2. Fifth row: subject PZTO03 (52 slices) with p=2 and σ=0. Sixth row: subject PZTO04 (51 slices) with p=2 and σ=1

Fig. 3 Post-processing pipeline for the patients CTTO01, CTTO03, and PZTO01 when applied to the outcomes of ‘GMFD 1ord’. First column: the starting data. Second colum: the outcome of the first step (starting slice: 20 for CTTO01, 44 for CTTO03, 25 for PZTO01). Third column: the outcome of the second step

Table 2 Reliability assessment for the three three-dimensional numerical schemes when compared to manual segmentation. The metrics utilized for this analysis are the Dice coefficient (optimal value: 1), the Surface Dice coefficient (optimal value: 1), the Jaccard coefficient (optimal value: 1), the Hausdorff distance, the Robust Hausdorff distance, and ASSD distance. The last two columns report the metrics values in the case of two standard software tools utilized in medical imaging, i.e., the TotalSegmentator 3D extension of Slicer in both fast and full-resolution mode (TS-FM and TS-FRM, respectively)

Metric	Classical	GMFD 1ord	GMFD 2ord	TS-FM	TS-FRM	
Dice	0.86 ± 0.03	0.87 ± 0.02	0.87 ± 0.03	0.83 ± 0.02	0.85 ± 0.05	
Surface Dice	0.88 ± 0.05	0.89 ± 0.02	0.88 ± 0.04	0.83 ± 0.03	0.84 ± 0.07	
Jaccard	0.76 ± 0.05	0.77 ± 0.04	0.77 ± 0.04	0.72 ± 0.03	0.75 ± 0.07	
Hausdorff	12.15 ± 3.85	12.95 ± 2.96	12.64 ± 3.76	29.18 ± 25.95	29.55 ± 15.58	
Robust Hausdorff	2.1 ± 0.7	2 ± 0.4	2 ± 0.5	2.8 ± 0.4	2.7 ± 0.9	
ASSD	0.40 ± 0.10	0.40 ± 0.04	0.42 ± 0.08	0.64 ± 0.35	0.66 ± 0.38	

We also tested the above described techniques on data downloaded from the Cancer Imaging Archive at https://www.cancerimagingarchive.net/. Specifically, we used six subjects collected in the Cancer Moonshot Biobank - Prostate Cancer Collection (CMB-PCA) [57, 58]. Figure 4 shows the segmentations of a dicom image of the psoas for one of the six subject in the dataset, while the average values of the metrics computed across the six subjects are shown in Table 3.

We finally point out that, on the one hand, the computational burden required by the ‘classical’ and ‘GMFD 1ord’ schemes is one order of magnitude smaller than the one required by ‘GMFD 2ord’. On the other hand, deep learning segmentation in fast mode (TS-FM) is typically performed within a computational time that is almost one order of magnitude smaller than the ones employed by the two fastest level set schemes (although with a rather significantly lower segmentation reliability).Fig. 4 Different segmentations for one of the six subjects in the CMB-PCA dataset. First row, from left to right: ground truth, TS-FM and TS-FRM segmentations. Second row, from left to right: ‘classical’, ‘GMFD 1ord’ and ‘GMFD 2ord’ segmentations

Table 3 Reliability assessment for the three three-dimensional numerical schemes when compared to manual segmentation on the six subjects in the CMB-PCA dataset. The metrics utilized for this analysis are the Dice coefficient (optimal value: 1), the Surface Dice coefficient (optimal value: 1), the Jaccard coefficient (optimal value: 1), the Hausdorff distance, the Robust Hausdorff distance, and ASSD distance. The last two columns report the metrics values in the case of two standard software tools utilized in medical imaging, i.e., the TotalSegmentator 3D extension of Slicer in both fast and full-resolution mode (TS-FM and TS-FRM, respectively)

Metric	Classical	GMFD 1ord	GMFD 2ord	TS-FM	TS-FRM	
Dice	0.84 ± 0.01	0.84 ± 0.01	0.84 ± 0.02	0.86 ± 0.02	0.86 ± 0.02	
Surface Dice	0.84 ± 0.05	0.84 ± 0.04	0.84 ± 0.07	0.87 ± 0.06	0.88 ± 0.05	
Jaccard	0.72 ± 0.2	0.73 ± 0.02	0.72 ± 0.06	0.75 ± 0.03	0.76 ± 0.03	
Hausdorff	20.56 ± 9.6	20.56 ± 9.6	21.77 ± 3.76	46.03 ± 15.37	48.21 ± 15.16	
Robust Hausdorff	2.4 ± 0.6	2.3 ± 0.6	2.4 ± 0.9	2.2 ± 0.7	2.2 ± 0.6	
ASSD	0.54 ± 0.18	0.53 ± 0.17	0.56 ± 0.25	0.88 ± 0.41	0.93 ± 0.43	

Comments and conclusions

This paper has proposed a 3D approach to the segmentation of the psoas muscle based on level set models. These numerical schemes provide a higher degree of automation with respect to the approach utilized in [22], which is notably heuristic, and a higher degree of generality with respect to the approach utilized in [20], whose performances rely on the availability of prior geometrical models for the psoas. For all numerical tests the geodesic models provide stable results, whereas the classical scheme may produce over-resolved shape, as it is particularly evident in the cases of patient CTTO02 (Fig. 2, first column, first row) and patient PZGE02 (Fig. 2, first column, third row). This effect is not evident for the two geodesic models, probably because the transport term in (4)-(5) pushes the front toward the external boundary of the muscle. Further, this same transport term usually helps the front stick to the right boundary, reducing the overflow outside the muscle, which happens when some complex structure is present in some portion of the data, as shown in Fig. 2, fifth row. As far as the reliability of the segmentation is concerned, the average Hausdorff and ASSD distances computed between the masks produced by our schemes and the manually identified ground-truth are, respectively, more than two times and more than 1.5 times lower than the distances obtained using the nnUNet-based TS masks (see Tables 2 and 3). This fact highlights that, even if the Jaccard and the Dice metrics indicate comparable overlap (of both volume and surface), our segmentation masks present less outlier. These results are confirmed by the use of the Surface Dice modification, which is 4-5% higher with respect to TS, and of the Robust Haussdorff distance, which is 22% lower, on average. In the case of the images downloaded from the CMB-PCA dataset, TS and the three level set numerical schemes present a similar degree of reliability if Dice, Surface Dice and Jaccard indicators are applied, while the deterministic numerical methods have a significantly better behaviour than deep learning according to the Haussdorff and ASSD metrics (see Table 3).

The first order version of the geodesic model has performances with comparable reliability to the ones characterizing the second order model. Indeed, the metrics and distances used in Table 2 have similar mean values and uncertainties. These same uncertainties are clearly smaller than the ones associated to the classical model. Further, the computational cost required by ’GMFD 1ord’ to achieve its metrics and distances values is one order of magnitude smaller than the one typical of the second order model (see Table 1). The higher computational burden of ’GMFD 2ord’ is mainly due to the parabolic CFL condition, which slows down the evolution and forces the algorithm to utilize a greater number of iterations to reach convergence.

It is clear that the simple post-processing procedure devised to clean the results provided by the numerical schemes is able to remove all spurious objects in most cases, still keeping the volume of the segmented muscle mostly untouched with respect to the results of segmentation. From a heuristic viewpoint, we could test that, if the outcome given by the scheme is already accurate on its own, one iteration of the procedure is enough to remove most artifacts. On the other hand, if the spurious objects are many and, in particular, some of them intersect the slice containing the center of the initial condition, a second run with a (possibly) different starting point is necessary. In any case, if the latter point is chosen correctly, the final result can be satisfactory and only in few cases some small objects still persist.

From a computational viewpoint, we observe that higher order schemes, as defined for example in [59], could be considered and adapted to the modified model (4) to increase the accuracy of the segmentation at the expense of longer computational times, but this goes beyond the scope of the present work and will be considered in future research.

The clinical impact of the present study can be of potential interest. In fact, in general, a systematic application of three-dimensional approaches for the segmentation of psoas can contribute to assess the limitations of the currently widespread clinical protocol, which is based on the analysis of a single cross-sectional image in correspondence to L3. More specifically, the use of deterministic numerical scheme like the ones based on level set methods may have some advantage with respect to AI-based models, since machine/deep learning may suffer data scarcity in the training phase.

In conclusion, on the one hand, this approach to segmentation of medical imaging data can be of course extended to the case of other muscular districts, although the psoas dimension and shape make this muscle particularly appropriate for automated segmentation in the case of sarcopenia assessment. On the other hand, the application of this approach in the case of other imaging modalities like MRI and PET can be problematic, due to the lower inter-organ contrast typical of MRI data and low spatial resolution typical of PET data. Finally, future work could include exploring adaptive parametric selection in the numerical schemes, integration with deep learning models, and applications to prospective clinical studies.

Abbreviations

CT X-ray computed tomography

PACS Picture Archiving and Communication System

AI Artificial intelligence

CNN Convolutional neural network

GAN Generative Adversarial Networks

FDG 18F-Fluorodeoxyglucose

PET Positron emission tomography

AMI Attenuation Metabolic Index

ALS Amyotrophic lateral sclerosis

CFL Courant-Friedrichs-Lewy

HU Hounsfield Unit

TS Total Segmentator

ASSD Average Symmetric Surface Distance

CMB-PCA Cancer Moonshot Biobank - Prostate Cancer Collection

Acknowledgements

G.P. and M.P. would like to dedicate this paper to the memory of Maurizio Falcone.

Authors' contributions

GP formulated the three numerical schemes and implemented the corresponding codes with the support of CC. CC applied the segmentation schemes to both the data set involved in the retrospective study at the San Martino Hospital and the data set downloaded from the Cancer Imaging Archive. IC computed the reliability coefficients in Tables 2 and 3. MP conceived the study, designed all the tests used for validation, and wrote the manuscript’s text with the support of all other authors.

Funding

CC and MP acknowledge the financial support of the “Hub Life Science - Digital Health (LSH-DH) PNC-E3-2022-23683267 - Progetto DHEAL-COM - CUP: D33C22001980001”, granted by the Italian Ministero della Salute within the framework of the Piano Nazionale Complementare to the “PNRR Ecosistema Innovativo della Salute - Codice univoco investimento: PNC-E.3”. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

The data utilized for this study can be provided by the corresponding author via email after motivated request. The MATLAB software implementing the three numerical schemes and the cleaning procedure is available at https://github.com/theMIDAgroup/psoas.

Declarations

Ethics approval and consent to participate

The CT images acquired for the nine subjects included in this paper belong to a retrospective data set of patients recruited at the IRCCS Ospedale Policlinico San Martino, Genova, Italy. All subjects were submitted to FDG-PET/CT. The study was approved by the Ethics Committees of IRCCS Ospedale Policlinico San Martino and performed according to the Declaration of Helsinki, Good Clinical Practice, and local ethics regulation. All enrolled patients signed a written informed consent at the time of FDG-PET/CT, econmpassing the use of anonymized data for retrospective research purposes.

Consent for publication

All enrolled patients signed a written informed consent for publication at the time of FDG-PET/CT, econmpassing the use of anonymized data for retrospective research purposes.

Competing interests

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Isabella Cama, Cristina Campi and Michele Piana contributed equally to this work.
==== Refs
References

1. Cruz-Jentoft AJ Bahat G Bauer J Boirie Y Bruyère O Cederholm T Sarcopenia: revised Europsensus on definition and diagnosis Age Ageing. 2019 48 1 16 31 10.1093/ageing/afy169 30312372
Cruz-Jentoft AJ, Bahat G, Bauer J, Boirie Y, Bruyère O, Cederholm T, et al. Sarcopenia: revised Europsensus on definition and diagnosis. Age Ageing. 2019;48(1):16–31.30312372
2. Kaido T Ogawa K Fujimoto Y Ogura Y Hata K Ito T Impact of sarcopenia on survival in patients undergoing living donor liver transplantation Am J Transplant. 2013 13 6 1549 1556 10.1111/ajt.12221 23601159
Kaido T, Ogawa K, Fujimoto Y, Ogura Y, Hata K, Ito T, et al. Impact of sarcopenia on survival in patients undergoing living donor liver transplantation. Am J Transplant. 2013;13(6):1549–56.23601159
3. Masuda T Shirabe K Ikegami T Harimoto N Yoshizumi T Soejima Y Sarcopenia is a prognostic factor in living donor liver transplantation Liver Transplant. 2014 20 4 401 407 10.1002/lt.23811
Masuda T, Shirabe K, Ikegami T, Harimoto N, Yoshizumi T, Soejima Y, et al. Sarcopenia is a prognostic factor in living donor liver transplantation. Liver Transplant. 2014;20(4):401–7.
4. Vergara-Fernandez O Trejo-Avila M Salgado-Nesme N Sarcopenia in patients with colorectal cancer: a comprehensive review World J Clin Cases. 2020 8 7 1188 10.12998/wjcc.v8.i7.1188 32337193
Vergara-Fernandez O, Trejo-Avila M, Salgado-Nesme N. Sarcopenia in patients with colorectal cancer: a comprehensive review. World J Clin Cases. 2020;8(7):1188.32337193
5. Chindapasirt J Sarcopenia in cancer patients Asian Pac J Cancer Prev. 2015 16 18 8075 7 10.7314/APJCP.2015.16.18.8075 26745041
Chindapasirt J. Sarcopenia in cancer patients. Asian Pac J Cancer Prev. 2015;16(18):8075–7.26745041
6. Pamoukdjian F Bouillet T Levy V Soussan M Zelek L Paillaud E Prevalence and predictive value of pre-therapeutic sarcopenia in cancer patients: a systematic review Clin Nutr. 2018 37 4 1101 1113 10.1016/j.clnu.2017.07.010 28734552
Pamoukdjian F, Bouillet T, Levy V, Soussan M, Zelek L, Paillaud E. Prevalence and predictive value of pre-therapeutic sarcopenia in cancer patients: a systematic review. Clin Nutr. 2018;37(4):1101–13.28734552
7. Collins J Noble S Chester J Coles B Byrne A The assessment and impact of sarcopenia in lung cancer: a systematic literature review BMJ Open. 2014 4 1 e003697 10.1136/bmjopen-2013-003697 24384894
Collins J, Noble S, Chester J, Coles B, Byrne A. The assessment and impact of sarcopenia in lung cancer: a systematic literature review. BMJ Open. 2014;4(1):e003697.24384894
8. Villasenor A Ballard-Barbash R Baumgartner K Baumgartner R Bernstein L McTiernan A Prevalence and prognostic effect of sarcopenia in breast cancer survivors: the HEAL Study J Cancer Survivorship. 2012 6 398 406 10.1007/s11764-012-0234-x
Villasenor A, Ballard-Barbash R, Baumgartner K, Baumgartner R, Bernstein L, McTiernan A, et al. Prevalence and prognostic effect of sarcopenia in breast cancer survivors: the HEAL Study. J Cancer Survivorship. 2012;6:398–406.
9. Drey M Hasmann SE Krenovsky JP Hobert MA Straub S Elshehabi M Associations between early markers of Parkinson’s disease and sarcopenia Front Aging Neurosci. 2017 9 53 10.3389/fnagi.2017.00053 28326036
Drey M, Hasmann SE, Krenovsky JP, Hobert MA, Straub S, Elshehabi M, et al. Associations between early markers of Parkinson’s disease and sarcopenia. Front Aging Neurosci. 2017;9:53.28326036
10. Yuksel H Balaban M Tan OO Mungan S Sarcopenia in patients with multiple sclerosis Mult Scler Relat Disord. 2022 58 103471 10.1016/j.msard.2021.103471 34998245
Yuksel H, Balaban M, Tan OO, Mungan S. Sarcopenia in patients with multiple sclerosis. Mult Scler Relat Disord. 2022;58:103471.34998245
11. Derstine BA Holcombe SA Ross BE Wang NC Su GL Wang SC Skeletal muscle cutoff values for sarcopenia diagnosis using T10 to L5 measurements in a healthy US population Sci Rep. 2018 8 1 11369 10.1038/s41598-018-29825-5 30054580
Derstine BA, Holcombe SA, Ross BE, Wang NC, Su GL, Wang SC. Skeletal muscle cutoff values for sarcopenia diagnosis using T10 to L5 measurements in a healthy US population. Sci Rep. 2018;8(1):11369.30054580
12. Inoue T Kitamura Y Li Y Ito W Ishikawa H Menze B Langs G Montillo A Kelm M Muller H Zhang S Psoas Major Muscle Segmentation Using Higher-Order Shape Prior Medical Computer Vision: Algorithms for Big Data 2016 Cham Springer International Publishing 116 124
Inoue T, Kitamura Y, Li Y, Ito W, Ishikawa H, et al. Psoas Major Muscle Segmentation Using Higher-Order Shape Prior. In: Menze B, Langs G, Montillo A, Kelm M, Muller H, Zhang S, et al., editors. Medical Computer Vision: Algorithms for Big Data. Cham: Springer International Publishing; 2016. p. 116–24.
13. Kamiya N Zhou X Chen H Muramatsu C Hara T Yokoyama R Automated segmentation of psoas major muscle in X-ray CT images by use of a shape model: preliminary study Radiol Phys Technol. 2012 5 5 14 10.1007/s12194-011-0127-0 21755416
Kamiya N, Zhou X, Chen H, Muramatsu C, Hara T, Yokoyama R, et al. Automated segmentation of psoas major muscle in X-ray CT images by use of a shape model: preliminary study. Radiol Phys Technol. 2012;5:5–14.21755416
14. Chen B Huang S Liang Z Chen W Pan B A fractional order derivative based active contour model for inhomogeneous image segmentation Appl Math Model. 2019 65 120 136 10.1016/j.apm.2018.08.009
Chen B, Huang S, Liang Z, Chen W, Pan B. A fractional order derivative based active contour model for inhomogeneous image segmentation. Appl Math Model. 2019;65:120–36.
15. Hashimoto F Kakimoto A Ota N Ito S Nishizawa S Automated segmentation of 2D low-dose CT images of the psoas-major muscle using deep convolutional neural networks Radiol Phys Technol. 2019 12 210 215 10.1007/s12194-019-00512-y 30937726
Hashimoto F, Kakimoto A, Ota N, Ito S, Nishizawa S. Automated segmentation of 2D low-dose CT images of the psoas-major muscle using deep convolutional neural networks. Radiol Phys Technol. 2019;12:210–5.30937726
16. Duong F Gadermayr M Merhof D Kuhl C Bruners P Loosen SH Automated major psoas muscle volumetry in computed tomography using machine learning algorithms Int J CARS. 2022 17 355 361 10.1007/s11548-021-02539-2
Duong F, Gadermayr M, Merhof D, Kuhl C, Bruners P, Loosen SH, et al. Automated major psoas muscle volumetry in computed tomography using machine learning algorithms. Int J CARS. 2022;17:355–61.
17. Kamiya N Li J Kume M Fujita H Shen D Zheng G Fully automatic segmentation of paraspinal muscles from 3D torso CT images via multi-scale iterative random forest classifications Int J CARS. 2018 13 1697 1706 10.1007/s11548-018-1852-1
Kamiya N, Li J, Kume M, Fujita H, Shen D, Zheng G. Fully automatic segmentation of paraspinal muscles from 3D torso CT images via multi-scale iterative random forest classifications. Int J CARS. 2018;13:1697–706.
18. Villarini B, Asaturyan H, Kurugol S, Afacan O, Bell JD, Thomas EL. 3D Deep Learning for Anatomical Structure Segmentation in Multiple Imaging Modalities. In: 2021 IEEE 34th International Symposium on Computer-Based Medical Systems (CBMS). New York City: IEEE; 2021. pp. 166–171.
19. Manabe T Ogawa C Takuma K Nakahara M Oura K Tadokoro T Usefulness of the Measurement of Psoas Muscle Volume for Sarcopenia Diagnosis in Patients with Liver Disease Diagnostics. 2023 13 7 1245 10.3390/diagnostics13071245 37046463
Manabe T, Ogawa C, Takuma K, Nakahara M, Oura K, Tadokoro T, et al. Usefulness of the Measurement of Psoas Muscle Volume for Sarcopenia Diagnosis in Patients with Liver Disease. Diagnostics. 2023;13(7):1245.37046463
20. Bauckneht M Lai R D’Amico F Miceli A Donegani MI Campi C Opportunistic skeletal muscle metrics as prognostic tools in metastatic castration-resistant prostate cancer patients candidates to receive Radium-223 Ann Nucl Med. 2022 36 4 373 383 10.1007/s12149-022-01716-w 35044592
Bauckneht M, Lai R, D’Amico F, Miceli A, Donegani MI, Campi C, et al. Opportunistic skeletal muscle metrics as prognostic tools in metastatic castration-resistant prostate cancer patients candidates to receive Radium-223. Ann Nucl Med. 2022;36(4):373–83.35044592
21. Zopfs D Theurich S Grosse Hokamp N Knuever J Gerecht L Borggrefe J Single-slice CT measurements allow for accurate assessment of sarcopenia and body composition Eur Radiol. 2020 30 1701 1708 10.1007/s00330-019-06526-9 31776743
Zopfs D, Theurich S, Grosse Hokamp N, Knuever J, Gerecht L, Borggrefe J, et al. Single-slice CT measurements allow for accurate assessment of sarcopenia and body composition. Eur Radiol. 2020;30:1701–8.31776743
22. Bauckneht M Lai R Miceli A Schenone D Cossu V Donegani MI Spinal cord hypermetabolism extends to skeletal muscle in amyotrophic lateral sclerosis: a computational approach to [18F]-fluorodeoxyglucose PET/CT images EJNMMI Res. 2020 10 1 1 10 10.1186/s13550-020-0607-5 31900594
Bauckneht M, Lai R, Miceli A, Schenone D, Cossu V, Donegani MI, et al. Spinal cord hypermetabolism extends to skeletal muscle in amyotrophic lateral sclerosis: a computational approach to [18F]-fluorodeoxyglucose PET/CT images. EJNMMI Res. 2020;10(1):1–10.31900594
23. Sambuceti G Brignone M Marini C Massollo M Fiz F Morbelli S Estimating the whole bone-marrow asset in humans by a computational approach to integrated PET/CT imaging Eur J Nucl Med Mol Imaging. 2012 39 1326 1338 10.1007/s00259-012-2141-9 22639281
Sambuceti G, Brignone M, Marini C, Massollo M, Fiz F, Morbelli S, et al. Estimating the whole bone-marrow asset in humans by a computational approach to integrated PET/CT imaging. Eur J Nucl Med Mol Imaging. 2012;39:1326–38.22639281
24. Fiz F Marini C Campi C Massone AM Podestà M Bottoni G Allogeneic cell transplant expands bone marrow distribution by colonizing previously abandoned areas: an FDG PET/CT analysis Blood J Am Soc Hematol. 2015 125 26 4095 4102
Fiz F, Marini C, Campi C, Massone AM, Podestà M, Bottoni G, et al. Allogeneic cell transplant expands bone marrow distribution by colonizing previously abandoned areas: an FDG PET/CT analysis. Blood J Am Soc Hematol. 2015;125(26):4095–102.
25. Marini C Morbelli S Cistaro A Campi C Caponnetto C Bauckneht M Interplay between spinal cord and cerebral cortex metabolism in amyotrophic lateral sclerosis Brain. 2018 141 8 2272 2279 10.1093/brain/awy152 30730551
Marini C, Morbelli S, Cistaro A, Campi C, Caponnetto C, Bauckneht M, et al. Interplay between spinal cord and cerebral cortex metabolism in amyotrophic lateral sclerosis. Brain. 2018;141(8):2272–9.30730551
26. Beltrametti MC Massone AM Piana M Hough transform of special classes of curves SIAM J Imaging Sci. 2013 6 1 391 412 10.1137/120863794
Beltrametti MC, Massone AM, Piana M. Hough transform of special classes of curves. SIAM J Imaging Sci. 2013;6(1):391–412.
27. Osher S Fedkiw RP Level set methods: an overview and some recent results J Comput Phys. 2001 169 2 463 502 10.1006/jcph.2000.6636
Osher S, Fedkiw RP. Level set methods: an overview and some recent results. J Comput Phys. 2001;169(2):463–502.
28. Li H Li P Gao L Zhang L Wu T A level set method for topological shape optimization of 3D structures with extrusion constraints Comput Methods Appl Mech Eng. 2015 283 615 635 10.1016/j.cma.2014.10.006
Li H, Li P, Gao L, Zhang L, Wu T. A level set method for topological shape optimization of 3D structures with extrusion constraints. Comput Methods Appl Mech Eng. 2015;283:615–35.
29. Rumpf M Preusser T A level set method for anisotropic geometric diffusion in 3D image processing SIAM J Appl Math. 2002 62 5 1772 1793 10.1137/S0036139901384662
Rumpf M, Preusser T. A level set method for anisotropic geometric diffusion in 3D image processing. SIAM J Appl Math. 2002;62(5):1772–93.
30. Pan S, Dawant BM. Automatic 3D segmentation of the liver from abdominal CT images: a level-set approach. In: Medical Imaging 2001: Image Processing, vol. 4322. SPIE; 2001. pp. 128–138.
31. Luo X, Chen J, Song T, Wang G. Semi-supervised medical image segmentation through dual-task consistency. In: Proceedings of the AAAI conference on artificial intelligence, vol. 35. Washington, DC: Association for the Advancement of Artificial Intelligence (AAAI); 2021. pp. 8801–8809.
32. Crespi L, Loiacono D, Sartori P. Are 3D better than 2D Convolutional Neural Networks for Medical Imaging Semantic Segmentation? In: 2022 International Joint Conference on Neural Networks (IJCNN). IEEE; 2022. pp. 1–8.
33. Mai DVC Drami I Pring ET Gould LE Lung P Popuri K A systematic review of automated segmentation of 3D computed-tomography scans for volumetric body composition analysis J Cachex Sarcopenia Muscle. 2023 14 5 1973 1986 10.1002/jcsm.13310
Mai DVC, Drami I, Pring ET, Gould LE, Lung P, Popuri K, et al. A systematic review of automated segmentation of 3D computed-tomography scans for volumetric body composition analysis. J Cachex Sarcopenia Muscle. 2023;14(5):1973–86.
34. Daoud MS Shehab M Al-Mimi HM Abualigah L Zitar RA Shambour MKY Gradient-based optimizer (GBO): a review, theory, variants, and applications Arch Comput Methods Eng. 2023 30 4 2431 2449 10.1007/s11831-022-09872-y 36597494
Daoud MS, Shehab M, Al-Mimi HM, Abualigah L, Zitar RA, Shambour MKY. Gradient-based optimizer (GBO): a review, theory, variants, and applications. Arch Comput Methods Eng. 2023;30(4):2431–49.36597494
35. Hell B Kassubeck M Bauszat P Eisemann M Magnor M An approach toward fast gradient-based image segmentation IEEE Trans Image Process. 2015 24 9 2633 2645 10.1109/TIP.2015.2419078 25838522
Hell B, Kassubeck M, Bauszat P, Eisemann M, Magnor M. An approach toward fast gradient-based image segmentation. IEEE Trans Image Process. 2015;24(9):2633–45.25838522
36. Clough JR Byrne N Oksuz I Zimmer VA Schnabel JA King AP A topological loss function for deep-learning based image segmentation using persistent homology IEEE Trans Pattern Anal Mach Intel. 2020 44 12 8766 8778 10.1109/TPAMI.2020.3013679
Clough JR, Byrne N, Oksuz I, Zimmer VA, Schnabel JA, King AP. A topological loss function for deep-learning based image segmentation using persistent homology. IEEE Trans Pattern Anal Mach Intel. 2020;44(12):8766–78.
37. Kawamoto M, Kamiya N, Zhou X, Kato H, Hara T, Fujita H. Simultaneous Learning of Erector Spinae Muscles for Automatic Segmentation of Site-Specific Skeletal Muscles in Body CT Images. IEEE Access. 2024;12:15468–76.
38. Crandall MG Lions PL Viscosity solutions of Hamilton-Jacobi equations Trans Am Math Soc. 1983 277 1 1 42 10.1090/S0002-9947-1983-0690039-8
Crandall MG, Lions PL. Viscosity solutions of Hamilton-Jacobi equations. Trans Am Math Soc. 1983;277(1):1–42.
39. Caselles V Kimmel R Sapiro G Geodesic active contours Int J Comput Vis. 1997 22 61 79 10.1023/A:1007979827043
Caselles V, Kimmel R, Sapiro G. Geodesic active contours. Int J Comput Vis. 1997;22:61–79.
40. Osher S Sethian JA Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations J Comput Phys. 1988 79 1 12 49 10.1016/0021-9991(88)90002-2
Osher S, Sethian JA. Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations. J Comput Phys. 1988;79(1):12–49.
41. Podgorsak EB. Radiation oncology physics: A handbook for teachers and students. International Atomic Energy Agency (IAEA). 2005.
42. Kaur M, Kaur J, Kaur J. Survey of contrast enhancement techniques based on histogram equalization. Int J Adv Comput Sci Appl. 2011;2(7):137–41.
43. Nixon M, Aguado A. Feature extraction and image processing for computer vision. Academic Press; 2019.
44. Qi Y, Yang Z, Sun W, Lou M, Lian J, Zhao W, et al. A comprehensive overview of image enhancement techniques. Arch Comput Methods Eng. 2022;29:583–607.
45. Deng G, Cahill L, An adaptive Gaussian filter for noise reduction and edge detection. In: 1993 IEEE conference record nuclear science symposium and medical imaging conference. IEEE; 1993. pp. 1615–9.
46. Geusebroek JM Smeulders AW Van De Weijer J Fast anisotropic gauss filtering IEEE Trans Image Process. 2003 12 8 938 943 10.1109/TIP.2003.812429 18237967
Geusebroek JM, Smeulders AW, Van De Weijer J. Fast anisotropic gauss filtering. IEEE Trans Image Process. 2003;12(8):938–43.18237967
47. Rosset A Spadola L Ratib O OsiriX: An Open-Source Software for Navigating in Multidimensional DICOM Images J Digit Imaging. 2004 17 3 205 216 10.1007/s10278-004-1014-6 15534753
Rosset A, Spadola L, Ratib O. OsiriX: An Open-Source Software for Navigating in Multidimensional DICOM Images. J Digit Imaging. 2004;17(3):205–16.15534753
48. Wasserthal J Breit HC Meyer MT Pradella M Hinck D Sauter AW TotalSegmentator: Robust Segmentation of 104 Anatomic Structures in CT Images Radiol Artif Intell. 2023 5 5 e230024 10.1148/ryai.230024 37795137
Wasserthal J, Breit HC, Meyer MT, Pradella M, Hinck D, Sauter AW, et al. TotalSegmentator: Robust Segmentation of 104 Anatomic Structures in CT Images. Radiol Artif Intell. 2023;5(5):e230024.37795137
49. Isensee F Jaeger PF Kohl SAA Petersen J Maier Hein KH nnU-Net: a self-configuring method for deep learning-based biomedical image segmentation Nat Methods. 2021 18 2 203 211 10.1038/s41592-020-01008-z 33288961
Isensee F, Jaeger PF, Kohl SAA, Petersen J, Maier Hein KH. nnU-Net: a self-configuring method for deep learning-based biomedical image segmentation. Nat Methods. 2021;18(2):203–11. 10.1038/s41592-020-01008-z.33288961
50. Rozynek M Tabor Z Klek S Wojciechowski W Body composition radiomic features as a predictor of survival in patients with non-small cellular lung carcinoma: A multicenter retrospective study Nutrition. 2024 120 112336 10.1016/j.nut.2023.112336 38237479
Rozynek M, Tabor Z, Klek S, Wojciechowski W. Body composition radiomic features as a predictor of survival in patients with non-small cellular lung carcinoma: A multicenter retrospective study. Nutrition. 2024;120:112336.38237479
51. Rozynek M, Gut D, Kucybala I, Strzalkowska-Kominiak E, Tabor Z, Urbanik A, et al. Fully automated 3D body composition analysis and its association with overall survival in head and neck squamous cell carcinoma patients. Front Oncol. 2023;13:1–8.
52. Gut D Tabor Z Szymkowski M Rozynek M Kucybala I Wojciechowski W Benchmarking of deep architectures for segmentation of medical images IEEE Trans Med Imaging. 2022 41 11 3231 3241 10.1109/TMI.2022.3180435 35666795
Gut D, Tabor Z, Szymkowski M, Rozynek M, Kucybala I, Wojciechowski W. Benchmarking of deep architectures for segmentation of medical images. IEEE Trans Med Imaging. 2022;41(11):3231–41.35666795
53. Dice LR Measures of the Amount of Ecologic Association Between Species Ecology. 1945 26 3 297 302 10.2307/1932409
Dice LR. Measures of the Amount of Ecologic Association Between Species. Ecology. 1945;26(3):297–302.
54. Jaccard P. The distribution of the flora in the alpine zone.1. New Phytol. 1912;11(2):37–50.
55. Hausdorff F Grundzuge der Mengenlehre 1914 Leipzig Aufl
Hausdorff F. Grundzuge der Mengenlehre. Leipzig: Aufl; 1914.
56. Heimann T van Ginneken B Styner MA Arzhaeva Y Aurich V Bauer C Comparison and Evaluation of Methods for Liver Segmentation From CT Datasets IEEE Trans Med Imaging. 2009 28 8 1251 1265 10.1109/TMI.2009.2013851 19211338
Heimann T, van Ginneken B, Styner MA, Arzhaeva Y, Aurich V, Bauer C, et al. Comparison and evaluation of methods for liver segmentation from CT datasets. IEEE Trans Med Imaging. 2009;28(8):1251–65.19211338
57. Cancer Moonshot Biobank - Prostate Cancer Collection (CMB-PCA) (Version 5). Cancer Imaging Arch. 2022. 10.7937/25T7-6Y12.
58. Clark K Vendt B Smith K Freymann J Kirby J Koppel P The Cancer Imaging Archive (TCIA): Maintaining and Operating a Public Information Repository J Digit Imaging. 2013 26 6 1045 1057 10.1007/s10278-013-9622-7 23884657
Clark K, Vendt B, Smith K, Freymann J, Kirby J, Koppel P, et al. The Cancer Imaging Archive (TCIA): Maintaining and Operating a Public Information Repository. J Digit Imaging. 2013;26(6):1045–57. 10.1007/s10278-013-9622-7.23884657
59. Falcone M Paolucci G Tozza S A high-order scheme for image segmentation via a modified level-set method SIAM J Imaging Sci. 2020 13 1 497 534 10.1137/18M1231432
Falcone M, Paolucci G, Tozza S. A high-order scheme for image segmentation via a modified level-set method. SIAM J Imaging Sci. 2020;13(1):497–534.
