Next Article in Journal
The Effect of Interventions on Preventing Musculoskeletal Injuries Related to Nurses Work: Systematic Review
Previous Article in Journal
Multi-Objective Genetic Algorithm for Cluster Analysis of Single-Cell Transcriptomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deep Learning-Based Segmentation and Volume Calculation of Pediatric Lymphoma on Contrast-Enhanced Computed Tomographies

by
Michał Klimont
1,2,3,*,
Agnieszka Oronowicz-Jaśkowiak
2,4,
Mateusz Flieger
2,
Jacek Rzeszutek
2,
Robert Juszkat
1 and
Katarzyna Jończyk-Potoczna
3
1
Department of Radiology, Poznań University of Medical Sciences, Dluga 1/2, 61-848 Poznan, Poland
2
Fast-Radiology, Poland
3
Department of Pediatric Radiology, Institute of Pediatrics, Poznań University of Medical Sciences, Szpitalna 27/33, 60-572 Poznan, Poland
4
1st Department of Radiology, National Institute of Oncology, W.K. Roentgena 5, 02-781 Warsaw, Poland
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2023, 13(2), 184; https://doi.org/10.3390/jpm13020184
Submission received: 15 November 2022 / Revised: 16 January 2023 / Accepted: 17 January 2023 / Published: 20 January 2023

Abstract

:
Lymphomas are the ninth most common malignant neoplasms as of 2020 and the most common blood malignancies in the developed world. There are multiple approaches to lymphoma staging and monitoring, but all of the currently available ones, generally based either on 2-dimensional measurements performed on CT scans or metabolic assessment on FDG PET/CT, have some disadvantages, including high inter- and intraobserver variability and lack of clear cut-off points. The aim of this paper was to present a novel approach to fully automated segmentation of thoracic lymphoma in pediatric patients. Manual segmentations of 30 CT scans from 30 different were prepared by the authors. nnU-Net, an open-source deep learning-based segmentation method, was used for the automatic segmentation. The highest Dice score achieved by the model was 0.81 (SD = 0.17) on the test set, which proves the potential feasibility of the method, albeit it must be underlined that studies on larger datasets and featuring external validation are required. The trained model, along with training and test data, is shared publicly to facilitate further research on the topic.

1. Introduction

Lymphomas are the most common blood malignancies in the developed world [1]. The two main categories of lymphomas are non-Hodgkin lymphomas (NHL) and Hodgkin lymphomas (HL) [1]. Worldwide, lymphomas are the ninth most common malignant neoplasms as of 2020 and were diagnosed in 627,439 persons, and caused 283,169 deaths [2].
Standardized staging and response criteria are key to successfully managing patients; what is more, they are essential to compare results and endpoints between studies when comparing treatment efficacy in a population.
There are multiple staging systems for both HL and NHL based on various criteria, including the anatomic disease extent and involvement of extra-nodal sites, as well as clinical and biochemical parameters [3]. The classification that is most commonly used in clinical practice is the Lugano staging classification, introduced in 2011 [4]. Factors taken into account include the number of lymph node regions involved, the presence of the disease on one or both sides of the diaphragm, the involvement of extranodal organs, the presence of systemic symptoms, and the presence of a bulky manifestation [4]. Evaluation of response, according to the Lugano classification, can be based on computed tomography (CT) alone or 18F-fluorodeoxyglucose (FDG) positron emission tomography/computed tomography (PET/CT) [4]. In PET/CT FDG, uptake is commonly assessed against regions of increased physiological activity using the Deauville score [5] and compared to baseline; the disease is then classified as complete, partial, or no metabolic response. In the Lugano criteria for FDG PET/CT, the volume of the metabolically active tumor tissue is not taken into account. CT assessment, on the other hand, is based on two-dimensional tumor measurement of up to six target lesions [4]. Tumor dimensions of each lymph node are multiplied, and these values are then added for all lesions and compared with the baseline [4]. Therefore, it can be argued that the value that is received is proportional to the surface area in the axial plane, which, in turn, reflects the volume of the lesion. The sum of the product for selected target lesions should be an approximate reflection of total tumor volume; however, it may be argued that it is less intuitive than volume for a human reader.
While relatively straightforward and easily applicable, the Lugano response criteria were primarily based on expert opinion, and their development has not been supported by large-scale data analysis [6]. What is more, the criteria for computed tomography restaging are elaborate, and the measurements are difficult to replicate, which results in a relatively high inter- and intra-observer variability [7]. Criteria based on FDG-avidity, on the other hand, are seen as controversial in lymphomas with variable FDG avidity, such as marginal zone lymphomas [4]. Additionally, in some cases, FDG-avid areas constituted only a small portion of the tumor mass (approximately 25%), and that monitoring based on FDG-avid areas potentially limits the prediction of the treatment sensitivity of the whole tumor mass to this small region [8]. In spite of those pitfalls, many lymphoma clinical trials continue to use the Lugano criteria as the best option that is available, albeit many introduce some modifications [9,10].
Because of these challenges in using Lugano and other classifications for staging, restaging, and follow-up of lymphomas, current guidelines recommend the use of FDG PET/CT [11], with the standardized glucose uptake value (SUV), used as a threshold to determine the metabolic tumor volume (MTV), defined as a total tumor volume demonstrating high glucose metabolism [12].
The use of a volume metric has numerous advantages over two-dimensional measurements, as proposed in the Lugano criteria. It eliminates the need for a subjective selection of a lesion and determining the best way of measuring it and, in the case of multiple tumors, the tedious process of identifying and comparing previously selected lesions. It can be hypothesized that using automatically determined tumor volume as a metric could result in a higher accuracy and a lower reading time compared to other methods and measurements performed by human readers. Indeed, it has been shown that PET/CT pretherapy metabolic tumor volume may be an independent prognostic factor in patients with some types of lymphoma, for example, large B-cell lymphoma [13]. It has also been shown that PET/CT can be used in both HL and NHL for monitoring the disease, although there is currently no consensus on cut-off points for both volume [3] and SUV threshold [14], and the positive influence of altering treatment on the basis of PET/CT results alone is not clear [3].
However, it is not always practical to perform PET/CT as a lymphoma follow-up. As mentioned above, some less common lymphoma types have varying FDG avidity, although pediatric lymphomas are more often FDG-avid than adult lymphomas. More importantly, PET/CT examinations may result in increased radiation exposure when compared to a stand-alone CT examination, as the effective dose is a combination of the dose from PET and CT [15], although it is worth noting that new hybrid PET/low dose CT scans do not have to translate to a higher radiation dose than CT. This is especially important when multiple follow-up examinations are required. What is more, PET/CT examination is significantly more expensive and less available, and the use of this examination for a regular follow-up for the sake of a somewhat better accuracy cannot be justified in many cases.
As shown above, both methods based on two-dimensional measurements taken by a radiologist or automatically computed MTV can be useful for staging and follow-up, but they have numerous limitations. From a practical standpoint, it seems that it would be ideal to have a method that would be cheaper and expose patients to less radiation than FDG PET/CT while at the same time being less prone to inter- and intra-observer variability than two-dimensional measurements taken subjectively by a radiologist on a regular CT scan. An alternative use of whole-body MRI is proposed; however, its application seems to be focused more on the initial staging, as diffusion-weighted imaging rather than facilitates finding suspected lymph nodes, than it helps in the treatment evaluation [16]. Thus, it seems appropriate to develop a volume based approach that would be based on CT scans.
It could be argued that calculating volume on plain CT scans does not take into account the metabolic state of the mass and may therefore be less reliable. However, it has been shown that the evolution of MTV and total tumor volume during treatment is similar for HL and that both FDG PET/CT avid and non-avid areas shrink at a similar rate [8].
With all the challenges described above, researchers are searching for alternative methods of volume measurements for lymphoma. To this date, researchers mostly focused their approach on a search for quick and precise estimates of the volume performed by the reporting radiologist that was not computer-aided [17]. However, with the advancement in machine learning, segmentation tasks in medical imaging can now be approached with fully automated methods. There are already solutions developed for anatomical structures segmentation, such as colon [18], lung [19], or even more robust multi-organ segmentation, such as the work by Wasserthal et al., which focused on segmenting 104 anatomical structures [20]. To our knowledge, automated volume calculation based on machine learning has not been used to calculate the volume of chest lymphoma.
The aim of this paper is to present a solution based on machine learning for automatic segmentation and volume calculation of chest lymphoma. Such a tool could theoretically assist both radiologists and clinicians in staging and follow-up of lymphoma in a more objective way than currently available methods. However, it should be emphasized that volume is one of many prognostic and staging factors. Importantly, the aim of this paper is also to release the dataset of segmented pediatric lymphomas for other researchers.

2. Materials and Methods

In total, 30 CT scans from 30 different pediatric patients diagnosed with any type of lymphoma were collected. Five patients were randomly selected as a test set and the rest of the cases were used for training and validation. All patients were hospitalized at Karol Jonscher University Hospital, Poznan, Poland between the years 2013 and 2020. Only one CT scan per patient was included to maximize potential variance in the data and increase model generalizability. Scans from the following two CT scanners were used: Siemens SOMATOM Definition AS+ and Siemens SOMATOM Force (Siemens Healthineers, Erlangen, Germany). In all cases, Visipaque (GE Healthcare, Chicago, IL, USA) was used as the contrast agent. The dosage of the contrast agent was calculated using the following formula: for children under 35 kg, a dose of 1 mL per kg of body weight was used; for children over 35 kg, a dose of 40–50 mL was used. The injection rate was between 1.5 and 2.5 mL/s. All scans were performed in portal venous phase.
There were 17 boys and 13 girls in the study group, and the mean age was 12.8 (SD 4.1, minimum age 2, maximum age 17). Patients’ age distribution is demonstrated in Figure 1. Median scan size was 296 × 512 × 512 voxels. Full dataset is publicly available at [21], model achieving the best Dice score is available for download at [22], and code for inference is available at [23].
The study, as well as making the scanning data publicly available, was approved by the local Bioethics Committee and the head of the Radiology Department. Furthermore, the patients and/-or their legal guardians consented to the data being used for retrospective research purposes. All data were fully anonymized.
Subsequently, two radiology residents manually segmented the thoracic lymphoma manifestations in 3D Slicer [24]. While performing the segmentations, multiple challenges were encountered, and it was attempted to create a set of rules about how to determine whether or not to include a lesion in the segmentation (Table 1).
To perform automated segmentations, we decided to explore the feasibility of applying deep learning. First developed by Ronneberger et al. in 2015 [25], the U-Net architecture has become a popular approach to segmentation tasks. Its impressive, robust performance has encouraged many researchers to develop modifications, further improving the results [26]. This has led to the development of a variety of new configurations (e.g., extending from two-dimensional to three-dimensional input), however, the parameters were generally preselected for a specific task and not generalizing well to every other problem. This makes parameter selection a time- and resource-consuming task. In 2018, Fabian Isensee et al. [27] developed nnU-Net (‘no-new-UNet’), which is a self-adapting framework that attempts to automatically prepare a well-performing configuration based on specific dataset properties. Not only does it automatically select certain parameters for the researcher, but it also streamlines the research process by providing tools for preparing and performing cross-validation, obtaining predictions for validation set, determining whether to apply post-processing, and deciding upon which, model (2-dimensional, 3-dimensional, 3-dimensional “cascade”, or a combination of them-i.e., ensemble), performs best. It is a popular approach to segmentation in biomedical tasks [26], and considering its simplicity and robustness, serves as a good reference baseline.
Great generalizability of nnU-Net across diverse datasets was demonstrated by Isensee et al. [27]. The framework was developed on 53 various segmentation tasks. Input data modalities included magnetic resonance imaging, computed tomography, electron microscopy, and fluorescence microscopy. The input format was either two- or three-dimensional. With such a broad scope of challenges, Isensee et al. were able to achieve new state-of-the-art results in 33 out of 53 tasks and presented results comparable to the top of the leaderboard for the remaining 20 tasks.
In this research, preprocessing and parameter selection followed the standard nnU-Net approach, which selects all hyperparameters automatically. The method of selection differs depending on the type of parameters, which can be divided into the following three groups: fixed, rule-based, and empirical. Fixed parameters, i.e., the parameters that remain the same for all applications, include learning rate, loss function, architecture template, optimizer, data augmentation, training procedure, and inference procedure. The rule-based parameters, i.e., parameters selected based on dataset properties, include the following: intensity normalization, image target spacing, network topology, patch size, batch size, trigger of 3D U-Net Cascade, configuration of low-resolution 3D U-Net [27]. Some of the rule-based parameters are co-dependent, as they have to meet GPU memory constraints, e.g., image size affects patch size, which later limits batch size. Empirical rules, i.e., parameters that are tested during cross-validation, include configuration of postprocessing and ensemble selection. Key parameters selected for our task are listed in Table 2.
The nnU-Net framework explores 3 different U-Net configurations to find the best result. These 3 configurations include a 2D U-Net, 3D U-Net, and 3D U-Net Cascade. In case of the 2D U-Net, each slice of the 3D image is fed into the network separately. The 3D U-Net operates on full 3D images by cropping 3D patches that cover voxel space. The 3D U-Net Cascade operates in two stages. During the first stage, segmentations are made on a downsampled version of the image. In the second stage, these segmentations are improved using both the segmentation obtained in the first stage and the full-resolution image as input. Cascade is omitted if the 3D U-Net patch size can accommodate a large part of the input image.
Performance of each of the configurations is explored with the 5-fold cross-validation, during which the dataset is consistently split 5 times into training and validation sets and then tested on the test (unseen during the training) data. Based on the cross-validation results, nnU-Net compares the results of each of the configurations and ensembles of them and decides which performed best. The inference is performed by averaging over predictions from models trained during different folds of cross-validation. The sliding window method is used with a window size equal to the training patch size and Gaussian importance weighting to reduce stitching artifacts. The model selection as well as the inference are performed automatically and do not require additional user interaction [27].
Obtained training time was oscillating around 2 days per fold per model. Computations were run using 10 CPUs, a single GPU card Tesla V100 (Nvidia, Santa Clara, CA, USA). Version 1.7.0 of nnU-Net package was used for performing experiments. The Dice coefficient was used as the evaluation metric. The Dice coefficient, also known as Sørensen-Dice coefficient or F1 score, is one of the most commonly used evaluation and validation metrices in medical imaging machine learning segmentation tasks [28]. The formula for the Dice coefficient is as follows:
2 T P 2 T P + F P + F N
where TP, FP, and FN represent the number of true positives, false positives, and false negatives, respectively. The Dice coefficient ranges from 0 to 1, with 1 meaning that two samples are identical, and 0 meaning that there are no mutual true positive data points.
In addition to being the standard metric for segmentation in medical imaging, the Dice coefficient was ideal for our application, as it is more focused on the total volume of the segmented lesion rather than a perfect reproduction of lesion borders [29]. As mentioned above, lymphoma borders could not always be identified in our dataset. What is more, the focus on volume reflects our suggested application of the model, which is comparison of segmentation volume in time.
Lymphoid volume for each test set segmentation was calculated using 3D Slicer software [24].

3. Results

The results of evaluating nnU-Net performance on each available model configuration can be found in Table 3. The highest Dice coefficient was obtained for the 3-D U-Net (without Cascade) model (0.7262 for the validation set). However, other model configurations achieved results that were comparable.
Table 4 presents the Dice coefficient for patients in the test set. The mean Dice coefficient for all the patients in the test set was 0.81 (SD = 0.17).
It should be noted that the model performed reasonably well on 4 out of 5 test cases, achieving a Dice score ranging from 0.73 to 0.95. One of the CT scans in the test set had a noticeably lower Dice score, i.e., 0.55. This was the result of the model interpreting part of the brain as a tumor-reference volume for this case was 146.19 cm3, and the predicted volume was 331.21 cm3. The volume of the mistakenly segmented brain tissue was larger than the reference volume of the tumor tissue, which translated to a low Dice score in this case. This was due to an unusually broad scanning range, and this is an example of the data coming from outside of the original distribution. Simple postprocessing by limiting the scanning range for that specific case would increase the Dice score greatly (to 0.92), but we refrained from that as the test set, in our opinion, should not be edited in any way to improve the results.
Exemplary segmentations are presented in Figure 2. It can be noticed that compared with manual segmentation, automatic segmentations have smoother borders. This reflects the challenges that are inherently associated with manual segmentation and clearly shows the potential advantage of machine learning-based approaches.

4. Discussion

Deep learning algorithms might be a solution to the growing number of imaging studies that greatly outpaces the number of radiologists being certified each year. With the aid of machine learning models, researchers develop ways of increasing the pace of reading imaging studies without compromising the quality, all the while reducing interobserver variability. For example, Lim et al. [30] recently published a study that applied deep learning algorithms for lumbar spine MRI assessment, which reduced reading time and improved the consistency of the stenosis rating between radiologists.
However, deep learning segmentation methods have their limitations that are, for the most part, independent of the research problem. Typical challenges include obtaining the ground truth segmentations, as this process is most often performed manually. This process is very time-consuming and may be prone to interobserver variability, as the borders of the lesions frequently cannot be discerned. Usually, a compromise needs to be made between the quality of the segmentation and the size of the dataset. In the following paragraphs, specific challenges encountered in our study and their consequences on our results are reported.
While evaluating the performance of nnU-Net on our dataset, some limitations were noticed that were mostly attributable to the limited size and diversity of the dataset. Firstly, there were cases of the algorithm mistaking parts of a solid organ, such as the brain or iliopsoas muscle, for tumor tissue. This issue occurred in patients with a broader than usual scanning scope. We believe this could be prevented by either including segmented whole-body CT scans from healthy individuals in the dataset or alternatively, by cropping the scope of the scans for existing cases to some standardized range. Although we are aware that not implementing these changes limits the use of our algorithm, it is beyond the scope of our project.
Another caveat of using human-made segmentations as the ground truth is that they can be quite far from the objective ground truth, in spite of best efforts. Perhaps the best illustration of the problem is presented by Tingelhoff et al. [31] as follows: in their experiment, 21 participants (10 ENT surgeons, 10 medical students, and one engineer) were asked to segment maxillary and ethmoid sinuses on the same CT scan. The total volume of the segmentation for the same patient varied between 30.9 cm3 and 47.1 cm3. This proves that the practical quality of even a relatively well-defined segmentation (as one would expect, would be the case for sinuses, being clearly limited by bones) is not as objective as it may first appear. What is more, manual segmentation is extremely time-consuming, and small improvements in quality may come with a significant increase in time per segmentation.
In addition to the ingrained inter-observer variability associated with any segmentation attempts, our team has encountered some further practical difficulties during the data preparation stage of the project, as described in the Materials and Methods section. This introduced further uncertainty to our human-made ground truth segmentation. The biggest challenges included the distinction between thymus and lymphoid tissue, as well as the presence of pericardial and pleural effusion and liquefactive necrosis within the tumor. However, it must be noted that even experienced radiologists might find the differentiation difficult while reviewing computed tomography. These issues translate into the model handling poorly some edge cases during cross-validation. Depending on the fold (and, therefore, on the training set subset), the model tended to be more likely to include or exclude effusion, liquefactive necrosis, and normal thymus in the automatic segmentation. Nonetheless, as the output of our model consists of voxel-wise segmentations, it can be edited by a reviewing radiologist to exclude any regions that were erroneously included. Providing editable, three-dimensional output segmentations might provide at least a partial solution to this problem.
It is important to underline that other researchers encounter similar difficulties. However, the positive impact of publishing imperfect training datasets on which other researchers can advance their research is worth emphasizing and promoting.
Our dataset included 30 CT scans from 30 different patients. To the best of our knowledge, all the previously published studies attempting automatic lymph node segmentation were also performed on small datasets, with the notable exception of the paper by Roth et al. [32], who published a large dataset of 176 abdominal CT scans with segmented lymph nodes. However, on manual inspection, it could be noted that many of the lymph nodes in the region of interest were not segmented, which is understandable, given the size of the dataset and the fact that, on average, there are 230–250 abdominal and pelvic lymph nodes [33]. To the best of our knowledge, there is no publicly available large dataset of segmented CT scans of patients with lymphoma. It can be hypothesized that the reason for this is a very large time investment on behalf of scientists attempting to create such a dataset, as one case may take several hours to segment. While it is clear that the scientific community would greatly benefit from large, publicly available, high-quality datasets with segmentations of various tissues, this project attempted to create a valuable dataset that would enable to train a reasonably well-performing model under practical constraints. It should be underlined that increasing the size of the dataset could decrease the variance in performance on the test set (including the model interpreting part of the brain tissue as a tumor, as seen on Patient 4 of the test set), and therefore, contribute to its overall better performance. What is more, the lack of appropriate publicly available datasets made it impossible to perform external validation on our model. We hope international scientific cooperation can help to alleviate some of these problems in the future.
In spite of being generally performed on small datasets, multiple attempts of automatic lymph node segmentation on CT scans have been reported. Early attempts employed numerical methods analyzing the shape, gray value, and borders of the nodes that generally required placement of a marker inside a lymph node that was to be segmented [34], although some models were also able to detect lymph nodes [35] without placing the marker. More recently, Iuga et al. [36] proposed a CNN for the detection and segmentation of thoracic lymph nodes in patients with possible lymph nodes metastasis. However, their approach differed from ours, as multiple small, non-pathological lymph nodes were generally assessed, and their primary benchmark was the number of detected lymph nodes and not their volume. The total detection rate was 69.9% for the validation dataset.
The development of new deep learning architectures has allowed a significant advancement in many applications of machine learning in medical imaging. For this project, the nnU-Net framework was used. This decision was based not only on the application of a well-designed deep learning model but also on a toolkit for pre-processing of the data, increasing the efficacy of hardware use, training process, cross-validation, supporting reproducibility, and more. Good deep learning practices that are an integral part of the nnU-Net guarantee that many common errors will be avoided and resulting models will be of high quality.

5. Conclusions

In this paper, a novel approach to a potentially important clinical task—volume calculation of thoracic lymphoma in pediatric patients—is presented. As discussed in the introduction, lymphoma volume correlates well with disease prognosis; therefore, accurate evaluation directly affects clinical management. We suggest that such tasks should be fully automated using recent advancements in imaging processing with machine learning methods. This study proves the feasibility of this approach and can serve as a building block for the further development of fine-tuned methods. It should be emphasized that there are some limitations of this study, including the size of the dataset and lack of external validation, both connected with the lack of public datasets. The nnU-Net, which was used for this project, is a freely available software package. The dataset used for the training and evaluation is published as a supplementary material with open access, and we encourage other researchers to use and modify it. To our best knowledge, it is the first, open, labeled dataset for pediatric lymphoma segmentation. As pointed out by Varoquaux et al. [37], dataset availability greatly influences areas of focus in a given field of research. Therefore, we expect other researchers to use this dataset freely to develop new machine learning algorithms, which will eventually bring benefit to the patients affected by this condition.
What is more, further prospective studies verifying the practical value of calculating the total tumor volume are needed, and the relationship between total volume change and clinical progression/remission should be documented.

Author Contributions

Conceptualization, M.K., A.O.-J., M.F., J.R. and K.J.-P.; Methodology, M.K., A.O.-J., M.F. and J.R.; Data collection, M.K. and K.J.-P.; Data analysis, M.K., A.O.-J., M.F. and J.R.; Writing—original draft preparation, M.K., A.O.-J., M.F. and J.R.; Review and Editing, M.K., A.O.-J., M.F. and J.R.; Supervision, K.J.-P. and R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Polish National Science Centre grant 2020/36/T/NZ5/00515.

Institutional Review Board Statement

The study as well as making the scanning data publicly available was approved by the local Bioethics Committee and the head of the Radiology Department. Furthermore, the patients and/-or their legal guardians consented for the data being used for retrospective research purposes. All data were fully anonymized.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are openly available in Zenodo at [21,22].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Morton, L.M.; Wang, S.S.; Devesa, S.S.; Hartge, P.; Weisenburger, D.D.; Linet, M.S. Lymphoma Incidence Patterns by WHO Subtype in the United States, 1992–2001. Blood 2006, 107, 265–276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  3. Kostakoglu, L.; Cheson, B.D. Current Role of FDG PET/CT in Lymphoma. Eur. J. Nucl. Med. Mol. Imaging 2014, 41, 1004–1027. [Google Scholar] [CrossRef] [PubMed]
  4. Cheson, B.D.; Fisher, R.I.; Barrington, S.F.; Cavalli, F.; Schwartz, L.H.; Zucca, E.; Lister, T.A.; Alliance, Australasian Leukaemia and Lymphoma Group; Eastern Cooperative Oncology Group; European Mantle Cell Lymphoma Consortium; et al. Recommendations for Initial Evaluation, Staging, and Response Assessment of Hodgkin and Non-Hodgkin Lymphoma: The Lugano Classification. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2014, 32, 3059–3068. [Google Scholar] [CrossRef]
  5. Barrington, S.F.; Mikhaeel, N.G.; Kostakoglu, L.; Meignan, M.; Hutchings, M.; Müeller, S.P.; Schwartz, L.H.; Zucca, E.; Fisher, R.I.; Trotman, J.; et al. Role of Imaging in the Staging and Response Assessment of Lymphoma: Consensus of the International Conference on Malignant Lymphomas Imaging Working Group. J. Clin. Oncol. 2014, 32, 3048–3058. [Google Scholar] [CrossRef] [PubMed]
  6. Younes, A.; Hilden, P.; Coiffier, B.; Hagenbeek, A.; Salles, G.; Wilson, W.; Seymour, J.F.; Kelly, K.; Gribben, J.; Pfreunschuh, M.; et al. International Working Group Consensus Response Evaluation Criteria in Lymphoma (RECIL 2017). Ann. Oncol. 2017, 28, 1436–1447. [Google Scholar] [CrossRef]
  7. Möller, J.; Steyn, T.; Combrinck, N.; Joubert, G.; Sherriff, A.; Rensburg, J.J.V. Inter-Observer Variability Influences the Lugano Classification When Restaging Lymphoma. SA J. Radiol. 2018, 22, 1357. [Google Scholar] [CrossRef]
  8. Girinsky, T.; Ghalibafian, M.; Bonniaud, G.; Bayla, A.; Magne, N.; Ferreira, I.; Lumbroso, J. Is FDG-PET Scan in Patients with Early Stage Hodgkin Lymphoma of Any Value in the Implementation of the Involved-Node Radiotherapy Concept and Dose Painting? Radiother. Oncol. 2007, 85, 178–186. [Google Scholar] [CrossRef]
  9. Batlevi, C.L.; De Frank, S.; Stewart, C.; Hamlin, P.A.; Matasar, M.J.; Gerecitano, J.F.; Moskowitz, A.J.; Hamilton, A.M.; Zelenetz, A.D.; Drullinsky, P.; et al. Phase I/II Clinical Trial of Ibrutinib and Buparlisib in Relapsed/Refractory Diffuse Large B-Cell Lymphoma, Mantle Cell Lymphoma, and Follicular Lymphoma. J. Clin. Oncol. 2018, 36, 7520. [Google Scholar] [CrossRef]
  10. Munakata, W.; Terauchi, T.; Maruyama, D.; Nagai, H. Revised Staging System for Malignant Lymphoma Based on the Lugano Classification. Jpn. J. Clin. Oncol. 2019, 49, 895–900. [Google Scholar] [CrossRef]
  11. Juweid, M.E. FDG-PET/CT in Lymphoma. In Positron Emission Tomography, 1st ed.; Juweid, M.E., Hoekstra, O.S., Eds.; Humana Press: Totowa, NJ, USA, 2011; pp. 1–19. ISBN 978-1-61779-062-1. [Google Scholar]
  12. Im, H.-J.; Bradshaw, T.; Solaiyappan, M.; Cho, S.Y. Current Methods to Define Metabolic Tumor Volume in Positron Emission Tomography: Which One Is Better? Nucl. Med. Mol. Imaging 2018, 52, 5–15. [Google Scholar] [CrossRef] [PubMed]
  13. Sasanelli, M.; Meignan, M.; Haioun, C.; Berriolo-Riedinger, A.; Casasnovas, R.-O.; Biggi, A.; Gallamini, A.; Siegel, B.A.; Cashen, A.F.; Véra, P.; et al. Pretherapy Metabolic Tumour Volume Is an Independent Predictor of Outcome in Patients with Diffuse Large B-Cell Lymphoma. Eur. J. Nucl. Med. Mol. Imaging 2014, 41, 2017–2022. [Google Scholar] [CrossRef] [PubMed]
  14. Schöder, H.; Moskowitz, C.H. Metabolic Tumor Volume in Lymphoma: Hype or Hope? J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2016, 34, 3591–3594. [Google Scholar] [CrossRef] [PubMed]
  15. Huang, B.; Law, M.W.-M.; Khong, P.-L. Whole-Body PET/CT Scanning: Estimation of Radiation Dose and Cancer Risk. Radiology 2009, 251, 166–174. [Google Scholar] [CrossRef]
  16. Littooij, A.S.; Kwee, T.C.; Barber, I.; Granata, C.; Vermoolen, M.A.; Enríquez, G.; Zsíros, J.; Soh, S.Y.; de Keizer, B.; Beek, F.J.A.; et al. Whole-Body MRI for Initial Staging of Paediatric Lymphoma: Prospective Comparison to an FDG-PET/CT-Based Reference Standard. Eur. Radiol. 2014, 24, 1153–1165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Baues, C.; Görgen, H.; Semrau, R.; Nast-Kolb, B.; Assenmacher, K.; Celik, E.; Morgenthaler, J.; Rosenbrock, J.; Trommer, M.; Houbois, C.; et al. Volumetric Assessment of Mediastinal Lymphoma Masses in Hodgkin Lymphoma. Leuk. Lymphoma 2019, 60, 3244–3250. [Google Scholar] [CrossRef]
  18. Gayathri Devi, K.; Radhakrishnan, R. Automatic Segmentation of Colon in 3D CT Images and Removal of Opacified Fluid Using Cascade Feed Forward Neural Network. Comput. Math. Methods Med. 2015, 2015, 670739. [Google Scholar] [CrossRef]
  19. Hofmanninger, J.; Prayer, F.; Pan, J.; Röhrich, S.; Prosch, H.; Langs, G. Automatic Lung Segmentation in Routine Imaging Is Primarily a Data Diversity Problem, Not a Methodology Problem. Eur. Radiol. Exp. 2020, 4, 50. [Google Scholar] [CrossRef]
  20. Wasserthal, J.; Meyer, M.; Breit, H.-C.; Cyriac, J.; Yang, S.; Segeroth, M. TotalSegmentator: Robust Segmentation of 104 Anatomical Structures in CT Images 2022. arXiv 2022, arXiv:2208.05868. [Google Scholar]
  21. Klimont, M.; Oronowicz-Jaśkowiak, A.; Flieger, M.; Rzeszutek, J.; Jończyk-Potoczna, K. Pediatric Chest Lymphoma with Segmentation. Zenodo 2022. [Google Scholar] [CrossRef]
  22. Klimont, M.; Oronowicz-Jaśkowiak, A.; Flieger, M.; Rzeszutek, J.; Jończyk-Potoczna, K. Pretrained Model for Pediatric Chest Lymphoma Segmentation. Zenodo 2022. [Google Scholar] [CrossRef]
  23. Deep Learning-Based Segmentation and Volume Calculation of Pediatric Lymphoma as an Alternative for Staging and Monitoring. Available online: https://github.com/fast-radiology/lymphoma (accessed on 13 November 2022).
  24. Fedorov, A.; Beichel, R.; Kalpathy-Cramer, J.; Finet, J.; Fillion-Robin, J.-C.; Pujol, S.; Bauer, C.; Jennings, D.; Fennessy, F.; Sonka, M.; et al. 3D Slicer as an Image Computing Platform for the Quantitative Imaging Network. Magn. Reson. Imaging 2012, 30, 1323–1341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ronneberger, O.; Fischer, P.; Brox, T. U-Net: Convolutional Networks for Biomedical Image Segmentation. In Proceedings of the International Conference on Medical Image Computing and Computer-Assisted Intervention, Munich, Germany, 5–10 October 2015; pp. 234–241. [Google Scholar]
  26. Siddique, N.; Paheding, S.; Elkin, C.P.; Devabhaktuni, V. U-Net and Its Variants for Medical Image Segmentation: A Review of Theory and Applications. IEEE Access 2021, 9, 82031–82057. [Google Scholar] [CrossRef]
  27. Isensee, F.; Jaeger, P.F.; Kohl, S.A.A.; Petersen, J.; Maier-Hein, K.H. NnU-Net: A Self-Configuring Method for Deep Learning-Based Biomedical Image Segmentation. Nat. Methods 2021, 18, 203–211. [Google Scholar] [CrossRef] [PubMed]
  28. Fastai. Available online: https://docs.fast.ai/index.html (accessed on 1 November 2022).
  29. Taha, A.A.; Hanbury, A. Metrics for Evaluating 3D Medical Image Segmentation: Analysis, Selection, and Tool. BMC Med. Imaging 2015, 15, 29. [Google Scholar] [CrossRef] [Green Version]
  30. Lim, D.S.W.; Makmur, A.; Zhu, L.; Zhang, W.; Cheng, A.J.L.; Sia, D.S.Y.; Eide, S.E.; Ong, H.Y.; Jagmohan, P.; Tan, W.C.; et al. Improved Productivity Using Deep Learning–Assisted Reporting for Lumbar Spine MRI. Radiology 2022, 305, 220076. [Google Scholar] [CrossRef]
  31. Tingelhoff, K.; Eichhorn, K.W.G.; Wagner, I.; Kunkel, M.E.; Moral, A.I.; Rilk, M.E.; Wahl, F.M.; Bootz, F. Analysis of Manual Segmentation in Paranasal CT Images. Eur. Arch. Otorhinolaryngol. 2008, 265, 1061–1070. [Google Scholar] [CrossRef]
  32. Roth, H.; Lu, L.; Seff, A.; Cherry, K.M.; Hoffman, J.; Wang, S.; Liu, J.; Turkbey, E.; Summers, R.M. A New 2.5 D Representation for Lymph Node Detection in CT. Med. Image Comput. Comput. Assist. Interv. 2014, 17, 520–527. [Google Scholar]
  33. Einstein, D.M.; Singer, A.A.; Chilcote, W.A.; Desai, R.K. Abdominal Lymphadenopathy: Spectrum of CT Findings. Radiographics 1991, 11, 457–472. [Google Scholar] [CrossRef] [Green Version]
  34. Dornheim, J.; Seim, H.; Preim, B.; Hertel, I.; Strauss, G. Segmentation of Neck Lymph Nodes in CT Datasets with Stable 3D Mass-Spring Models: Segmentation of Neck Lymph Nodes. Acad. Radiol. 2007, 14, 1389–1399. [Google Scholar] [CrossRef]
  35. Feulner, J.; Kevin Zhou, S.; Hammon, M.; Hornegger, J.; Comaniciu, D. Lymph Node Detection and Segmentation in Chest CT Data Using Discriminative Learning and a Spatial Prior. Med. Image Anal. 2013, 17, 254–270. [Google Scholar] [CrossRef] [PubMed]
  36. Iuga, A.-I.; Carolus, H.; Höink, A.J.; Brosch, T.; Klinder, T.; Maintz, D.; Persigehl, T.; Baeßler, B.; Püsken, M. Automated Detection and Segmentation of Thoracic Lymph Nodes from CT Using 3D Foveal Fully Convolutional Neural Networks. BMC Med. Imaging 2021, 21, 69. [Google Scholar] [CrossRef] [PubMed]
  37. Varoquaux, G.; Cheplygina, V. Machine Learning for Medical Imaging: Methodological Failures and Recommendations for the Future. Npj Digit. Med. 2022, 5, 48. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Age distribution of the patients included in the training and testing dataset.
Figure 1. Age distribution of the patients included in the training and testing dataset.
Jpm 13 00184 g001
Figure 2. Three-dimensional renderings of the segmentations performed by a human reader (red) and from the best performing nnU-Net model (green). Smooth edge segmentation can be noted in the segmentations performed by the nnU-Net when compared to the human reader.
Figure 2. Three-dimensional renderings of the segmentations performed by a human reader (red) and from the best performing nnU-Net model (green). Smooth edge segmentation can be noted in the segmentations performed by the nnU-Net when compared to the human reader.
Jpm 13 00184 g002
Table 1. Challenges encountered by the radiologists during segmentation and solutions used in the database.
Table 1. Challenges encountered by the radiologists during segmentation and solutions used in the database.
Region/IssueSolution
Cervical lymph nodes not always distinguishable from surrounding tissuesInclude cervical lymph nodes whenever possible
Unsharp border between lymphoma and thymic tissueExclude thymus from segmentation only when a clear border between lymphoma and thymus is visible; include thymus in segmentation when no clear border is visible
Unsharp borders between lymphoma/liquefactive necrosis and fluid in pericardium and pleural cavitiesTry to exclude any pericardial and pleural effusion and include liquefactive necrosis in the segmentation (difficult in some cases)
Abdominal lymph nodesDo not include in the segmentation
Table 2. Key parameters of the nnU-Net framework. SGD-stochastic gradient descent.
Table 2. Key parameters of the nnU-Net framework. SGD-stochastic gradient descent.
ParameterValue
Batch size2D: 12
3D: 2
Float precision 16-bitYes
Max number of epochs * 1000
Number of batches per epoch * 250
Number of input channels 1
Initial learning rate * 0.01
Momentum * 0.99
Optimizer * SGD
Patch size 2D: 512 × 512
3D: 96 × 160 × 160
Weight decay * 0.00003
* fixed parameters.
Table 3. Model performance measured by the average Dice coefficient during cross-validation. The model with the highest dice score is highlighted.
Table 3. Model performance measured by the average Dice coefficient during cross-validation. The model with the highest dice score is highlighted.
ModelAverage Dice Coefficient
2D U-Net0.7065
3D U-Net0.7262
3D U-Net Cascade0.7024
2D U-Net + 3D U-Net0.7221
2D U-Net + 3D U-Net Cascade0.7203
3D U-Net + 3D U-Net Cascade0.7148
Table 4. The test set results for the 3D U-Net model. The dice coefficient, reference volume based on the manual segmentation, volume predicted by the model, and volume difference are presented for each patient.
Table 4. The test set results for the 3D U-Net model. The dice coefficient, reference volume based on the manual segmentation, volume predicted by the model, and volume difference are presented for each patient.
PatientDiceManual Segmentation [cm3]Automatic Segmentation [cm3]Volume Difference [cm3]
Patient 10.88288.79257.6831.11
Patient 20.73631.34865.01−233.67
Patient 30.92776.99686.1490.85
Patient 40.55146.19331.21−185.02
Patient 50.95354.63352.092.54
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Klimont, M.; Oronowicz-Jaśkowiak, A.; Flieger, M.; Rzeszutek, J.; Juszkat, R.; Jończyk-Potoczna, K. Deep Learning-Based Segmentation and Volume Calculation of Pediatric Lymphoma on Contrast-Enhanced Computed Tomographies. J. Pers. Med. 2023, 13, 184. https://doi.org/10.3390/jpm13020184

AMA Style

Klimont M, Oronowicz-Jaśkowiak A, Flieger M, Rzeszutek J, Juszkat R, Jończyk-Potoczna K. Deep Learning-Based Segmentation and Volume Calculation of Pediatric Lymphoma on Contrast-Enhanced Computed Tomographies. Journal of Personalized Medicine. 2023; 13(2):184. https://doi.org/10.3390/jpm13020184

Chicago/Turabian Style

Klimont, Michał, Agnieszka Oronowicz-Jaśkowiak, Mateusz Flieger, Jacek Rzeszutek, Robert Juszkat, and Katarzyna Jończyk-Potoczna. 2023. "Deep Learning-Based Segmentation and Volume Calculation of Pediatric Lymphoma on Contrast-Enhanced Computed Tomographies" Journal of Personalized Medicine 13, no. 2: 184. https://doi.org/10.3390/jpm13020184

APA Style

Klimont, M., Oronowicz-Jaśkowiak, A., Flieger, M., Rzeszutek, J., Juszkat, R., & Jończyk-Potoczna, K. (2023). Deep Learning-Based Segmentation and Volume Calculation of Pediatric Lymphoma on Contrast-Enhanced Computed Tomographies. Journal of Personalized Medicine, 13(2), 184. https://doi.org/10.3390/jpm13020184

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop