Next Article in Journal
Langmuir Trough Study of the Interactions of Tear Mimetic Eyedrop Formulation with Human Meibum Films
Next Article in Special Issue
2D Winograd CNN Chip for COVID-19 and Pneumonia Detection
Previous Article in Journal
Modified Artificial Bee Colony Based Feature Optimized Federated Learning for Heart Disease Diagnosis in Healthcare
Previous Article in Special Issue
Artificial Intelligence, Sensors and Vital Health Signs: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Survey on Recent Trends in Medical Image Classification Using Semi-Supervised Learning

by
Zahra Solatidehkordi
* and
Imran Zualkernan
Department of Computer Science and Engineering, American University of Sharjah, Sharjah 26666, United Arab Emirates
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(23), 12094; https://doi.org/10.3390/app122312094
Submission received: 6 November 2022 / Revised: 18 November 2022 / Accepted: 23 November 2022 / Published: 25 November 2022

Abstract

:
Training machine learning and deep learning models for medical image classification is a challenging task due to a lack of large, high-quality labeled datasets. As the labeling of medical images requires considerable time and effort from medical experts, models need to be specifically designed to train on low amounts of labeled data. Therefore, an application of semi-supervised learning (SSL) methods provides one potential solution. SSL methods use a combination of a small number of labeled datasets with a much larger number of unlabeled datasets to achieve successful predictions by leveraging the information gained through unsupervised learning to improve the supervised model. This paper provides a comprehensive survey of the latest SSL methods proposed for medical image classification tasks.

1. Introduction

Recent years have seen an increase in the digitalization of health records as well as a gradual increase in computational power, enabling the use of machine learning and deep learning algorithms in various areas of medicine such as diagnosis, outcome prediction, and clinical decision-making [1,2]. Medical images are a major component of patients’ electronic health records, making medical image analysis an active field of research for machine learning and deep learning [2]. Currently, the extraction of information from medical images is performed manually by radiologists and other professionals in the healthcare field, making the process dependent on the individual’s experience and knowledge, vulnerable to human error, and time-consuming. In order to address these issues, extensive research has been performed on the usage of machine learning and deep learning in medical image analysis to achieve more efficient and precise results [3]. The tasks studied in the context of medical image analysis include segmentation, detection, de-noising, reconstruction, and classification [4]. This survey paper focuses on recent studies in image classification only.
Image classification involves the assignment of one or more labels to an image based on the features extracted from it. The process consists of various steps, including preprocessing, feature extraction, feature selection, and classification. The extracted features can be low-level such as color, shape, intensity, texture, boundary, and position information, as well as high-level such as bag-of-words, scale-invariant feature transform, and fisher vector [5,6]. In traditional machine learning methods, features are generally extracted separately and then fed into the classification model. On the other hand, deep learning algorithms allow the feature extraction and classification steps to be combined in one network [5]. In addition, deep learning models apply end-to-end learning, in which the model receives a dataset of labeled images and extracts descriptive, hierarchical, and highly representative features with respect to each label which the model then uses for the classification task [5,7]. Deep learning methods can combine complex and low-level features and are less vulnerable to human errors [6]. Research has shown that deep learning models frequently outperform traditional machine learning algorithms in image classification tasks [5,8,9,10]. However, deep learning methods have their own limitations, such as the requirement for more time and higher computing power. Furthermore, traditional machine learning methods may achieve better results on smaller datasets [11,12,13].
Medical image classification plays an essential role in the analysis and diagnosis of diseases as well as clinical care and treatment. Medical images can be of different modalities based on the biomedical devices used to generate them. These modalities include Projectional Radiography (X-rays), Computed Tomography (CT), Magnetic Resonance Imaging (MRI), and Ultrasound Imaging [14]. X-ray is the most common medical imaging technique in which electromagnetic waves are used to create a two-dimensional representation of the body’s internal structure based on the varying rates of wave absorption in tissues with different densities. CT is a computerized procedure used to generate cross-sectional X-ray images, which can be combined to create a three-dimensional image, providing a more detailed representation of the region compared to conventional X-rays [15]. MRI generates three-dimensional anatomical images based on the behavior of the atomic nuclei in response to powerful magnetic fields and radiofrequency pulses [16]. Ultrasound imaging uses high-frequency sound waves to generate two-dimensional images of organs and tissues based on the reflection of the sound waves from different body structures [17]. In addition to these, dermoscopic images are also often used for medical image classification. These are two-dimensional images visualizing the subsurface skin structures. The benchmark datasets most used for research on medical image classification are the National Institutes of Health (NIH) Chest X-ray14 dataset and the International Skin Imaging Collaboration 2018 (ISIC2018) Skin Lesion Analysis dataset. The ChestX-ray14 dataset consists of 112,120 chest x-rays of 30,805 patients. 53.9% of the images are normal, and 46.1% are abnormal. The abnormal images are labeled with 14 different pathologies. Each image can have one or more labels, with 40.1% of abnormal images being multi-labeled, making the dataset suitable for multi-label classification. The ISIC2018 dataset consists of 10,015 skin lesion dermoscopy images labeled with seven types of skin lesions. Each image belongs to one class, making the dataset suitable for multi-class classification.
One major issue facing medical image analysis is the lack of large datasets of high-quality labeled samples because manual labeling of images is a labor-intensive, expensive task [18]. However, obtaining unlabeled medical data from clinical sites is relatively easier [19]. Therefore, semi-supervised learning (SSL) is appropriate for medical image analysis because SSL uses a combination of both labeled and unlabeled data for learning. The goal of SSL is to use unsupervised learning with unlabeled data to improve the performance of the supervised model developed from labeled data [20]. This paper provides a survey of some of the latest SSL methods developed for the purpose of medical image classification.
This paper is organized as follows: Section 2 discusses commonly used SSL approaches, Section 3 explains the survey methodology, Section 4 presents the review of recent SSL approaches, Section 5 offers a summary of the results, and Section 6 discusses gaps and future work.

2. Background

This section provides background information about commonly used SSL approaches. Most current SSL methods rely on three assumptions. First, the smoothness assumption which states that if two points in a high-density region are close, their output labels must also be close. Second, the cluster assumption states that the points belonging to one class tend to form a cluster, meaning if certain points are in the same cluster, they likely belong to the same class. Lastly, the manifold assumption states that samples that are close in high-dimensional space are also close in a low-dimensional manifold [21,22]. The state-of-the-art (SOTA) semi-supervised methods for medical image classification can be divided into three broad categories of consistency-based, graph-based and adversarial [19,23]. Each of these is described next.

2.1. Consistency-Based Methods

In consistency-based methods, a model is encouraged to give consistent results on perturbed versions of the same inputs, such as input images with Gaussian noise added [21]. In other words, if an input x i is classified as class c, then the perturbed input x i must also be classified as class c. As previous research has shown, an ensemble of neural networks can considerably improve the ability of the network to generalize [24,25,26]. One of the most dominant consistency-based SSL models is the Mean Teacher (MT) [27], which is based on the Temporal Ensembling [24] method. This method involves self-ensembling in neural networks, wherein the results of the same network over different epochs are combined during training to infer new training targets from the unlabeled data. Each training epoch has the same network but different regularization and input augmentation conditions to enforce consistency. Due to the use of dropout regularization, the prediction generated from each epoch is from an individual sub-network, making the final result an ensemble prediction of many different sub-networks. An exponential moving average (EMA) prediction is maintained for each training instance and updated after each training step with the new predictions. As a result, the EMA prediction is the ensemble result of the model generated after the latest training step and all the models generated in the previous steps. The predictions obtained through the temporal self-ensembling method are closer to the real labels of the unlabeled part of the dataset, and as such, these inferred labels can be used as training targets for unlabeled data. The Mean Teacher (MT) [27] builds on the temporal self-ensembling method. The model takes on the role of a student and a teacher; as a student, the model normally learns, while as a teacher, the model creates targets for learning which the student then uses to train. The teacher is an ensemble of the student model, similar to the previous method. This model uses the average weights of the models as opposed to the average of the predictions. Thus, the teacher maintains the EMA weights of the student model. At training step t, the EMA( θ t ) of successive θ weights is defined as:
θ t = α θ t 1 + 1 α θ t
where α is a smoothing coefficient hyperparameter. Mean Teacher maintains a total cost function which is a weighted sum of the consistency cost and classification cost. The consistency cost is calculated by comparing the outputs of the teacher and student, while the classification cost involves comparing the student’s prediction and the target label. The Mean Teacher method allows the model to work with large datasets and improves classification accuracy.

2.2. Adversarial Methods

We will discuss two commonly used methods called Virtual Adversarial Training (VAT) and Generative Adversarial Network (GAN) methods in this category.
Adversarial training methods involve the presence of “adversarial attacks,” which cause various errors in prediction [28]. These attacks are small perturbations to the input data applied in the direction where the model’s label assignment is most sensitive [29]. In VAT, adversarial attacks are generated and incorporated into the training of the model in order to increase its robustness. The aim of this method is to achieve high smoothness at each data point in order to prevent small perturbations from causing classification errors. The loss function used is the local distributional smoothness (LDS) loss, a measure of the local distributional smoothness which indicates the robustness of the model against local perturbation around each input data point. The LDS is calculated as follows:
L D S x * , θ : = D p ( y | x * , θ ^ ) , p y | x * + r v a d v , θ
wherein x * is a labeled or unlabeled data point and D p , p is a function that evaluates the divergence between the distributions p and p . p ( y | x * , θ ^ ) is an approximation of the true distribution of the output label calculated based on probabilistically generated “virtual” labels. Thus, the goal is to minimize the divergence between the distribution of x * and its transformed version x * + r v a d v , making the model more robust to adversarial perturbations. The full objective function of VAT is as follows:
D l , θ + α R v a d v D l , D u l , θ
wherein D l , θ is the negative log-likelihood of the labeled dataset, and α is the regularization coefficient which controls the relative balance between D l , θ and R v a d v . R v a d v is the regularization term which is defined as the average of the loss L D S x , θ over all input data points.
In GAN-based methods, the adversary is a generator network whose role is to create fake images that fit the distribution of the real data as closely as possible. This generator competes against a discriminator network whose role is to distinguish the fake data from the real data. As the two networks compete, they both gradually improve. The optimal state of the model is reached when the generator has learned to create perfectly realistic samples which the discriminator cannot distinguish from the real data [30]. This ability to synthesize samples with high levels of realism can be used to solve the issue of scarcity of labeled data. Additionally, the classifier training can also be embedded directly into the GAN framework for improved generalization and robustness [4].
Generative Adversarial Nets (GANs) [31] are deep neural network architectures involving two networks; a generator model (G) and a discriminator model (D). The name “adversarial” comes from the fact that G and D are pitted against each other; G generates samples while D learns to distinguish between the samples generated by G and samples from the real data. The competition between the two networks leads to both networks improving over time, with G generating more realistic images increasingly until they are undistinguishable from the real samples.
Deep Convolutional Generative Adversarial Networks (DCGAN) [32] is a type of GAN that uses Convolutional Neural Networks (CNN) in its generator in order to create a higher quality generator model and achieve stable training across different datasets. Some of the architecture guidelines of DCGAN include the replacement of pooling layers with stride convolutions in D and fractional-strided convolutions in G, the usage of Batch Normalization in both G and D, the elimination of fully connected hidden layers for deeper architectures, the usage of the ReLU activation in G for all layers except the output which uses the Tanh function, and lastly, usage of LeakyReLU activation in D for all layers.
Auxiliary Classifier Generative Adversarial Networks (ACGAN) [33] is a model based on GAN that uses label conditioning in image synthesis, resulting in higher quality samples with increased global coherence. In this method, every sample created by G has a class label in addition to noise. In turn, D has to estimate both the likelihood of the image being correct and the likelihood of the label is correct. The addition of label conditioning in the model results in higher-quality generated images and more stabilized training.

2.3. Graph-Based Methods

In graph-based methods, graphs are constructed based on the data, where nodes represent samples (labeled or unlabeled), and edges represent neighborhood relations between the samples [22]. Generally, graph-based methods in SSL work by propagating labels from labeled samples to unlabeled samples based on the information gained from the graph [21]. Graph-based techniques are popular in recent works as they tend to be convex, scalable, and effective [22]. Additionally, their strong mathematical properties allow for fast computation [34]. A graph is represented by the notation G = V ,   E with V denoting the nodes of the graph and E represents edges between two nodes. Each node is a sample from the data and each edge represents the relationship between two samples. Graphs are called directed if the edges have directions and undirected if they do not. Graphs can also be weighted in a weighted graph G = V ,   E ,   W , w i j represents the weight of the edge between nodes i and j, and W is a matrix of size N × N which contains the weights of all edges. If w i j = w j i , the W matrix is symmetric. If w i j = 0 , this signifies that the nodes i and j are not connected [22,35]. The degree of a vertex v i is defined as d i = j = 1 n w i j . The degree matrix D is a diagonal matrix in which the values on the diagonal are d 1 , d n [35].
The unnormalized Graph Laplacian matrix is defined as L = D W , while the normalized graph Laplacian matrix is defined as:
L sym : = D 1 / 2 L D 1 / 2 = I D 1 / 2 W D 1 / 2
where I is an identity matrix.
The Graph p-Laplacian is a nonlinear generalization of the graph Laplacian, which aims to improve the representation of the local structure information [36]. The unnormalized and normalized matrix forms of the graph p-Laplacian operator Δ p are as follows:
Δ p ( u ) i = j V M i j φ p f i f j
Δ p ( u ) i = 1 d i j V M i j φ p f i f j
where φ p x = | x | p 1 sign x , d i is the degree function of the graph, M i j is the weight of an edge, f is the eigenvector of the graph p-Laplacian operator, and p 1 is a parameter set by the user.

2.4. Other Methods

In addition to the three main categories, we also review papers implementing methods that are less common in the recent literature. For example, one such method is pseudo-labeling, where a supervised learning model is trained on the labeled data and then used to classify a part of the unlabeled data. The pseudo-labels generated with high confidence are then added to the labeled data, and the model has trained again [37]. Traditional pseudo-labeling methods have several disadvantages, such as their tendency to be biased toward the majority class and their inability to be seamlessly adapted to multi-label and multi-class problems. Confidence-based pseudo-labeling leads to the selection of samples from the majority class as they tend to have higher confidence, resulting in a majority-biased classification model [37]. As such, pseudo-labeling methods are usually used in conjunction with consistency-based, adversarial, or graph-based methods in order to mitigate these disadvantages. In addition to pseudo-labeling, papers addressing Federated Learning, Learning by Association, Contrastive Learning, and MixMatch will also be reviewed.
Federated Learning (FL) is a decentralized system that allows training to be conducted on distributed data sources such as different hospitals separately. The parameters of the locally trained models can be used to update the main algorithm without exchanging private information, making FL an appropriate privacy-reserving method of data collaboration between different institutions. In semi-supervised FL, each of the local datasets can be labeled or unlabeled [38].
The Learning by Association (LA) [39] model assumes that good embeddings that belong to the same class will have high similarity. It builds association cycles through a “walk” from the embedding vector A obtained from a batch of labeled data and embedding vector B obtained from a batch of unlabeled data. The ”walker” goes from A to B and then back to A based on mutual similarities. If walker ends up in the same class that it started from, the cycle is considered correct. Thus, the goal is to maximize the possibility of a correct association cycle. The “walker” loss is defined as:
L walker   = H T , P a b a
with H representing the cross-entropy between the association cycle probabilities P a b a and the uniform target distribution of correct association cycles T . Additionally, Haeusser et al. [39] argue that instead of making associations only among “easy” samples, it is beneficial to “visit” all unlabeled samples in order to make the best use of the dataset and encourage embeddings that generalize better. Thus, the visit loss L v i s i t is proposed, which is the cross-entropy between the visit probabilities P v i s i t and the uniform target distribution V:
L v i s i t = H V , P v i s i t
Contrastive Learning (CL) [40] is a form of self-supervised learning which involves generating augmentations of datapoints to learn the higher-level features of a dataset from unlabeled samples. The augmentations of the same datapoint are considered positive samples and encouraged to be close to each other, while the augmentations of different datapoints are considered negative samples and encouraged to be far apart. The goal of this process is to create a model with higher generalization capability across different data distributions. Joint Contrastive Learning (JCL) [41] is a form of CL in which an infinite number of key-query pairs are generated from unlabeled data. As opposed to each key-query pair being compared independently as in the original CL method, multiple key-query pairs containing positive keys paired with the same shared query q are encouraged to be close. This leads to more consistent representations in each instance-specific class.
MixMatch is a recent SSL method proposed by Berthelot et al. [42] that aims to unify three different loss terms, such as entropy minimization loss, which encourages the generation of labels with high confidence. The consistency regularization loss encourages output consistency in the presence of perturbations, and the generic regularization loss encourages better generalization. Consistency regularization is achieved by augmenting both labeled and unlabeled images with standard augmentation techniques (random horizontal flips and crops for images) and enforcing consistency in their results. Entropy minimization is achieved by sharpening the target distribution for unlabeled data. Lastly, traditional regularization is applied using a weight decay that penalizes the L2 norm of the model parameters. Additionally, the Mixed Sample Data Augmentation (MSDA) method MixUp [43] is used for both regularization and semi-supervised learning. MSDA methods involve mixing different samples and their labels based on a particular policy to create an augmented dataset. The MixUp approach involves training the model on convex combinations of pairs of images and their labels to encourage linear behavior in-between training samples. Thus, the algorithm works as follows. A batch of labeled data X and a batch of unlabeled data U is fed to the model, and augmentations are applied to both. The label predictions of all augmentations of U are computed and averaged, and sharpening is applied to the average prediction. The augmented labeled dataset X ^ and the augmented unlabeled dataset with guessed labels U ^ are combined and shuffled to obtain set W. Then, MixUp is applied to the labeled dataset X ^ and entries from W to obtain the processed labeled dataset X ; in the same way, MixUp is applied to the unlabeled dataset U ^ and the rest of W to obtain the processed unlabeled dataset U . After the processed labeled and unlabeled (with guessed labels) batches X and U are obtained, they are used for semi-supervised training with the loss function L X + λ U L U which combines the supervised loss L X and the unsupervised loss L U with the hyperparameter λ U . L X is the cross-entropy loss obtained from training on X and L U is the squared L2 loss obtained from training on U′. FixMatch [44] is an SSL method similar to MixMatch that combines different existing SSL approaches, including consistency regularization and pseudo-labeling. FixMatch encourages label consistency between the pseudo-label of a weakly-augmented version of an input image and the pseudo-label of a strongly augmented version of the same image. Comatch [45] is another`` SSL method that combines entropy minimization, consistency regularization, contrastive learning, and graph-based SSL. In this method, the training data is represented by a class probability and a low-dimensional embedding which evolve in a co-training framework. A smoothness constraint is applied to the class probabilities to improve pseudo-labels which are then used to train the classification head and projection head with a cross-entropy loss and graph-based contrastive loss respectively.

3. Survey Methodology

We used a variety of sources to search for journal articles and conference papers. These sources include the following databases and digital libraries: IEEE Xplore, ACM Digital Library, Springer, Elsevier, ScienceDirect, ResearchGate, Google Scholar, Scopus, MDPI online database, and PubMed. The distribution of papers read from different sources can be seen in Figure 1. The main keywords used for our searches were medical, image, and semi-supervised. Other searches included the addition of the terms “analysis” and “classification.” Additionally, for each category, the keywords related to that category were added to the search; consistency, adversarial, graph, GANs. Results were filtered for articles published between the years 2019 and 2022, and sorting was conducted based on relevance and number of citations when possible.
The exclusion and inclusion of research articles were decided first by a preliminary abstract analysis, then by a full review of the article. The research papers which focused only on medical image segmentation and did not contain a dedicated section on classification were excluded. The inclusion criteria for our research were as follows:
  • The study must have a primary focus on semi-supervised learning.
  • The study must include a clear report on the performance of the classification model.
  • The study must present an in-depth description of the model architecture.
The exclusion criteria of our research were as follows:
  • The study is not peer-reviewed or indexed in a reputable database.
  • The study does not propose any significant addition or change to previously existing deep learning models.
  • The study presents vague descriptions of the experimentation and classification results.
For the popular benchmark datasets X-ray14 and ISIC2018, the models which did not achieve results comparable to existing state-of-the-art models were excluded.
As the latest survey on this topic covered papers published in 2018 and prior [46], our paper presents research published between 2019 and 2022. Additionally, while [46] offered a general overview of unsupervised and semi-supervised methods in medical image analysis, this paper is specifically about semi-supervised learning for classification tasks and provides a more in-depth description of each of the models presented.

4. Survey Results

We have categorized the methods that were reviewed into five types of SSL approaches: consistency-based, adversarial, graph-based, other approaches, and hybrid methods. Figure 2 shows a diagram of the methods belonging to each category as well as the various datasets used for experimentation. The application of each method is discussed next.

4.1. Consistency-Based Methods

Relation-driven Self-Ensembling Model (SRC-MT) [19]: This method employs the Mean Teacher framework with a consistency-enforcing approach. In addition, this method explores the intrinsic relation between the images, which is generally overlooked in consistency-based methods, including Mean Teacher. Studying the relation between the images in an unsupervised fashion helps extract useful information from the unlabeled data [47,48]. The SRC-MT method proposes a novel Sample Relation Consistency (SRC) paradigm which enforces the consistency of the relation between two images after perturbation; if two images have a similarity before perturbation, this relation should be preserved after they are perturbed. In other words, the relationship between input samples x 1 and x 2 must be the same as the relationship between perturbed samples x 1 and x 2 . Thus, this method achieves consistency in both labeling and relations after perturbation. The total objective functions of the framework are as follows:
L = L s + λ L u ,   with   L u = L c + β L s r c
wherein L s is the supervised objective and L u the unsupervised objective. L u is made up of the standard consistency loss L c , and the relation consistency loss L s r c . λ is the trade-off weight between unsupervised and supervised loss, and β is a hyperparameter to balance L c and L s r c .
SRC-MT was tested on two public benchmark medical image datasets; thorax disease classification with ChestX-ray14 for multilabel classification and ISIC2018: Skin Lesion Analysis for single label classification. Perturbations such as random rotation, horizontal and vertical flips, and translation were applied to the input images for consistency enforcing. For the skin legion dataset, the images were all resized to 224 × 224. The dataset was divided into 70% training, 10% validation, and 20% testing. DenseNet121 [23], pre-trained on ImageNet [49], was used as the network backbone. For the ChestX-ray14 dataset, the images were resized to 384 x 384, and the same data split as the previous dataset was used. DenseNet169 [23], pre-trained on ImageNet, was used as the network backbone. For comparison, the upper bound and baseline performances were calculated; the upper bound was the performance of the fully supervised model trained on 100% labeled data, and the baseline was the performance of the fully supervised model trained on 20% of the data. Additionally, the model was compared to other semi-supervised learning methods, including the self-training-based method [50], SS-DCGAN [51], the Temporal Ensembling method, and the Mean Teacher. The training set used for the proposed method consisted of 20% labeled and 80% unlabeled data. The metrics used to gauge the performance of the model for single-label classification was AUC, sensitivity, specificity, accuracy, and F1 scores. For the single-label classification task on the skin legion dataset, SRC-MT achieved the highest AUC, sensitivity, accuracy, and F1 scores compared to the baseline and other semi-supervised models. SRC-MT also achieved the second-highest specificity. Additionally, the AUC score of SRC-MT (93.58) was 3% higher than the baseline and 2% less than the upper bound and the F1 score (60.68) was 8% more than the baseline and 10% less than the upper bound. For the multi-labeled classification on the Chest X-ray14 dataset, the metric used was AUC, and for comparison the upper bound and baseline as well as the GraphXNET [34] method were used. The AUC scores were calculated for each of the 14 labels and averaged. The SRC-MT achieved a score of 79.23, approximately 0.4% higher than the Mean Teacher, 1.2% higher than GraphXNET, 2% higher than the baseline, and 2% lower than the upper bound.
NoTeacher (NoT) [52]: In the Mean Teacher, the consistency target, which is the teacher, relies on the EMA of the student. In other words, the teacher’s weight is an ensemble of student weights. As such, the weights of the student and teacher are tightly coupled. This leads to the model being vulnerable to confirmation bias [53]; the teacher continues to reinforce what it already believes. NoTeacher is a consistency-based SSL framework that addressed this issue by employing two independent networks, thus removing the need for the time-averaging component. The NoTeacher framework works as follows: two random augmentations are applied to an input x , which can be a 2D or 3D sample. These augmentations result in samples x 1 and x 2 , which becomes the input of F 1 and F 2 , two networks with similar architectures. The outputs are named f 1 L and f 2 L for labeled inputs and f 1 U and f 2 U for unlabeled inputs. Next, a loss function is calculated to enforce prediction consistency between F 1 and F 2 . This loss function consists of a consistency loss and a supervised cross-entropy loss; if x 1 and x 2 are augmented versions of the same input x , then f 1 resulting from x 1 and f 2 resulting from x 2 must be similar. Additionally, if x is a labeled input, then the outputs of both networks should match the ground truth target y . Lastly, the total loss is backpropagated to update network parameters. The NoTeacher method is like the Mean Teacher as they both employ two networks with similar architectures. However, the NoTeacher method introduces two differences; first, it removes the EMA, making the networks entirely independent, and second, NoTeacher bases its loss function on a graphical model. The graphical model involves a graph consisting of f 1 , f 2 and y as the nodes, each connected to a consensus function f c 0 , 1 , which enforces mutual agreement of the outputs on labeled and unlabeled data. The three datasets used for testing this model were Chest X-ray14, Radiological Society of North America (RSNA) Brain CT, and the Knee MRNet dataset containing knee MRI exams. The backbone networks used for each dataset were DenseNet121 [23], DenseNet169 [23], and MRNet [54], respectively, all pre-trained on ImageNet. All inputs were normalized based on ImageNet statistics before training. For all the datasets, the experiment was repeated six or seven times with different percentages of the dataset being labeled. The Chest X-ray14 dataset was used for multi-label classification. It was divided into 70% training, 10% validation, and 20% testing. The AUROC scores of the method were compared to the supervised baseline as well as other semi-supervised methods, including Mean Teacher, VAT, Pseudo-labeling, and the previously discussed SRC-MT. The NoTeacher method consistently achieved higher scores than the previous methods. It achieved a score of 82.10 at 50% labeled data, 1% less than the supervised model at 100% labeled data. The SRC-MT method attained lower scores than MT at all labeling budgets below 20%, but a higher score at 20%. The RSNA Brain CT dataset was used for multi labeled classification. It includes 19,520 CT brain scans from the RSNA 2019 challenge. Each slice in these scans is labeled with five binary labels corresponding to the absence or presence of five different forms of hemorrhage. 14.3% of the images include the presence of at least one form of hemorrhage, with 30.1% including multiple forms. 85.7% of the data are samples with no abnormalities. The dataset was divided into 60% training, 20% validation and 20% testing. NoTeacher achieved consistently higher scores than both the supervised baseline and MT. It also achieved overall higher scores than the other semi-supervised methods. On the other hand, MT provided lower scores than the supervised baseline at higher percentages of labeled data. Lastly, the Knee MRNet dataset consists of 1370 data points labeled with the presence or absence of abnormality, leading to a single labeled classification task. 80.6% of the data points contain abnormalities. 1130 training scans and 120 testing scans were used, with 20% of the training scans being randomly sampled for the validation set. Similar to the previous experiments, NoTeacher achieved higher scores than both the supervised baseline and the other semi-supervised methods. Additionally, at the labeled data percentage of 50%, NoTeacher managed to achieve a score higher than the supervised baseline at 100% labeled data.
Self-supervised Mean Teacher for Semi-supervised (S2MTS2) [55]: S2MTS2 is a consistency-based SSL method for chest x-ray classification, which uses the Mean Teacher framework. This method consists of two learning stages. In the first stage, the student-teacher model is pre-trained in a self-supervised manner on the labeled and unlabeled data using JCL [41]. JCL involves learning an infinite number of key-query pairs generated from unlabeled data with the goal of enforcing dependencies between different key-query pairs containing the same query, which leads to more consistent representations in each instance-specific class [41]. Thus, each query q i paired simultaneously with multiple of its positive keys k i , m + must return a low loss value. The loss of each q i , k i , m + pair is defined as follows:
L i , m = l o g e x p 1 τ q i k i , m + e x p 1 τ q i k i , m + + j = 1 K e x p 1 τ q i k i , j
where τ is the temperature hyperparameter, k i , m + is the m t h positive key of q i and k i , j is the j t h negative key of q i . Thus, the total JCL loss can be calculated as follows:
p D X , θ 2 , θ 2 = 1 D X 1 M i = 1 D X m = 1 M L i , m
where D X is the set of labeled and unlabeled images, and M is the number of positive keys. In the second stage, the pre-trained student-teacher model is fine-tuned using the original Mean Teacher approach by maintaining an EMA following Equation (1). The datasets used for testing were the Chest X-ray14, ISIC2018 Skin Lesion Analysis, and CheXpert. The CheXpert dataset contains 220,000 images labeled with 14 diseases, wherein each image can have more than one label. The backbone model used for all datasets was DenseNet121 [23], with the two-layer multi-layer perceptron (MLP) projection head being replaced with a three-layer MLP. The Chest X-rat14 dataset was resized to 512 × 512, and the ISIC2018 and CheXpert datasets were resized to 224 × 244. The data augmentations used were random resize, random crop, random horizontal flipping, and random rotation. The Chest X-ray14 dataset was evaluated on different percentages of labeled data (2%, 5%, 10%, 15%, 20%), resulting in AUC scores of 74.69, 78.96, 79.90, 80.31, 81.06, and 82.51. For the CheXpert dataset, the model achieved AUC scores of 66.15, 67.85, 70.83, 71.37, and 71.58 for 100, 200, 300, 400, and 500 labeled samples, respectively. Lastly, the results obtained on the ISIC2018 dataset with 20% labeled samples were an AUC score of 94.71% and an F1 score of 62.67.
A summary of the methods discussed in this section is provided in Table 1. Three consistency-based methods were reviewed: SRC-MT, NoT, and S2MTS2. For each model, the table contains the datasets used for testing as well as the best performance score obtained on each dataset.

4.2. Adversarial Methods

Uncertainty-Guided Virtual Adversarial Training With Batch Nuclear-Norm Optimization [56]: This method incorporates batch-nuclear norm (BNN) optimization [57] to prevent overfitting to the labeled data and to improve the diversity and discriminability of the model. As proposed by Cui et al. [57], the nuclear-norm P θ * of the m × n prediction matrix P θ is calculated as follows:
P ( θ ) * = i = 1 m i n ( m , n ) σ i ( P ( θ ) )
where σ i P θ represents the ith largest singular value of P θ . As mentioned previously, the goal of incorporating BNN optimization is to improve generalization and avoid overfitting the labeled data. This is achieved by minimizing the BNN of the labeled data and maximizing the BNN of the unlabeled data batch. Thus, the labeled BNN loss L l B N N and unlabeled BNN loss L u B N N are defined as follows:
L l B N N = 1 B l P l ( θ ) *
L u B N N = 1 B u P u ( θ ) *
B l and B u represent the size of the labeled and unlabeled datasets, while P l ( θ ) * and P u ( θ ) * represent the nuclear norm of the labeled and unlabeled prediction matrices, respectively. In addition to BNN, the model incorporates uncertainty guidance during VAT loss computation to filter out unlabeled samples close to the decision boundary. The uncertainty U i is calculated for each unlabeled sample x u i in a batch, and predictions with high uncertainty are filtered out to ensure that only reliable targets are used for learning. U i is calculated as follows:
U i = j = 1 c p u i , j l o g p u i , j . i 1 B u
with p u i , j representing the predicted probability of x u i for the jth class and C denoting the number of classes. Thus the model is trained using the following losses: the cross-entropy loss of the supervised model ( L c l s ), the VAT loss calculated from labeled data ( L l v a t ), the uncertainty-guided VAT loss computed from the unlabeled data ( L ˜ u v a t ), and the BNN losses L l B N N and L u B N N . The overall loss for the labeled data is the combination of all losses calculated over the labeled data:
L l = L c l s + λ v a t L l v a t + λ l B N N L l B N N
Similarly, the unlabeled loss can be calculated as follows:
L u = λ v a t L ˜ u v a t + λ u B N N L u B N N
wherein λ v a t , λ l B N N and λ u B N N are weighting coefficients. The overall objective function is the addition of the supervised and unsupervised losses L l + L u . The model was tested on the ISIC2018 and Chest X-ray14 datasets as well as an in-house collected hip x-ray dataset composed of 26,075 x-ray images labeled into four categories. The ISIC2018 and hip x-ray images were resized into 224 × 224 images. DenseNet121 [23], pre-trained on ImageNet, was used as the backbone network for both datasets. The X-ray14 dataset was resized into 256 × 256 images with DenseNet169 [23] as the backbone network. The model achieved an AUC score of 96.04 and an F1 score of 69.67 on the ISIC2018 dataset with 20% labeled data. On the Chest X-ray14 dataset, the method achieved AUC scores of 69.75, 74.50, 77.52, 79.49, and 80.69 on labeled percentages of 2%, 5%, 10%, 15%, and 20%. Lastly, the hip x-ray dataset generated an AUC score of 92.43 and an F1 score of 72.39 with 20% labeled data. Other state-of-the-art methods such as MT, SRC-MT, and VAT were also tested on this dataset achieving AUC scores of 92.00, 91.40, 91.66, and F1 scores of 70.52, 70.90, and 70.57, respectively.
Semi-supervised Adversarial Classification (SSAC) [58]: this GAN-based method involves a reconstruction network R, and a supervised classification network C. Learnable transition (T) layers are used to transfer the image representation ability learned by network R to C. R is an adversarial autoencoder-based unsupervised network consisting of a generator G and a discriminator D. G contains an encoder and a decoder. The encoder, which has an architecture similar to ResNet [59], takes 64 × 64 patches, and the decoder generates reconstructed patches of the same size. D is a deep convolutional neural network containing four convolutional layers. C contains two parts, one with the same architecture as the encoder of R, and one a fully connected layer with two neurons, separated by a global average pooling (GAP) layer. There is no parameter sharing between R and C; the learnable T layers, with each layer consisting of a 1 × 1 convolutional layer, transfers the feature maps obtained by R to the corresponding blocks in C. In the experimentation, C was pre-trained on ImageNet, and R was pre-trained on both the labeled and unlabeled data. The loss function was defined as follows:
L S S A C X m = λ 1 l m s e G X m , X m l b c D G X m , 0 + l b c D X m , 1 + l b c C X m , Y m
where X m is the mth input sample and λ 1 is a weighting parameter. The three terms in the function represent the mean square reconstruction loss of G, the cross-entropy adversarial loss of D, and the supervised classification loss, respectively. Two datasets were used for experimentation. The first dataset was a database of lung nodules on CT-scans from the lung image database consortium (LIDC), and the image database resource initiative (IDRI) called LIDC-IDRI. The second dataset was the Tianchi Lung Nodule Detection dataset. The LIDC-IDRI dataset contains 1018 chest CT scans from which the nodules of diameter 3mm to 30 mm were identified and annotated. From the Tianchi Lung Nodule Detection dataset, 1000 patients’ data was used to extract 1227 unlabeled nodules. In total, the combination of the two datasets resulted in 644 malignant, 1301 benign, and 1839 unlabeled samples, making this a single-label classification task. For this task, the MK-SSAC model was constructed, which consists of nine knowledge-based collaborative SSAC models (KBC-SSAC) trained on the patches extracted from the lung nodules whose predictions were subsequently combined by two output neurons. The model was tested five independent times with 10-fold cross-validation. The metrics calculated to gauge the results were accuracy, sensitivity, specificity, and AUC. The performance of the model was compared to other SSL methods (MK-CatGAN, MK-AAE, MK-Ladder Network), hand-crafted feature-based methods (3D GLCM + SVM, MVF + SVM), and DCNN-based methods (Fuse-TSD, MV-KBC). The average scores achieved by the model were 92.53 accuracies, 84.94 sensitivity, 96.28 specificities, and 95.81 AUC. The model obtained the highest accuracy and AUC scores as well as the second-highest sensitivity and specificity among all the methods compared. Additionally, the MK-SSAC model was tested with different percentages of the labeled data (100%, 80%, 60%, 40%, and 20%). From 20% to 100%, the scores increased by about 2% for accuracy and 3% for AUC. In comparison, the supervised MK-C model was also tested with the same percentages of labeled data and showed an increase in 7% in accuracy and 5% in AUC. These results show that the MK-SSAC model has relative robustness to variation in the number of labeled data.
Retinal Image Synthesis and Semi-Supervised Learning for Glaucoma Assessment (SS-DCGAN) [51]: SS-DCGAN is a model created to synthesize retinal images and determine the presence or absence of glaucoma. This model is based on the DCGAN architecture with a modification in the last output layer of D: one neuron for synthesis and three neurons for training the glaucoma classifier. D thus becomes a classifier that has to label each sample as either Normal, Glaucoma, or Synthetic. The loss function of the method was as follows:
L = L supervised   + L unsupervised
where L supervised for K classes is the cross-entropy loss function:
L supervised   = E x , y p data   x , y l o g p model   ( y | x , y < K + 1 )
and L unsupervised is GAN’s two-player minmax game:
L unsupervised   = E x p data   x l o g D x + E z p z z l o g 1 D G z
where D x is the probability of x being from the real data and G z is the probability of z being from the generator. The dataset used for experimentation was a combination of fourteen public databases with a total of 86,926 images. Samples with glaucoma or normal labels were divided into a 70% training set (669 glaucomas and 981 normal) and a 30% testing set (287 glaucomas and 420 normal). The labeled training set, as well as all unlabeled images (84,569 samples), were used to train the model. The metrics calculated were AUC, F1-score, specificity, and sensitivity. The final results were compared to three previous models: ResNet50 [59] as well as the CNNs proposed by Alghamdi et al. [60] and Chen et al. [61]. SS-DCGAN achieved the highest AUC (0.90) and F1-score (0.84), as well as the second-highest specificity and sensitivity scores among the models.
Bi-Modality Medical Image Synthesis Using Semi-Supervised Sequential Generative Adversarial Networks [62]: multi-modality imaging involves the simultaneous usage of two or more imaging modalities for one examination [63]. For example, single photon emission computed tomography (SPECT), magnetic resonance imaging (MRI), and positron emission tomography (PET) can be used to combine optical, magnetic, and radioactive reporters to detect abnormalities in the epileptic brain. Common examples of bi-modality images are PET-SPECT and PET-CT [63]. The model proposed by Yang et al. was intended to synthesize high-quality bi-modality medical images using GANs by creating two sequential generative networks, each corresponding to one modality. The first modality in the sequence is the one with the least complexity, measured automatically by a complexity measurer algorithm. The goal of this is to decrease the difficulty of generating the second, more complex modality by conditioning it on the first. The proposed generator network is trained via semi-supervised learning in order to generate realistic images with a large diversity. The supervised learning model learns the joint distribution of the different modalities, while the unsupervised learning model learns the marginal distribution of the different modalities via adversarial learning. The architecture of the generator is as follows: first, a real image of modality 1 is encoded into a low-dimensional latent vector, which is then mapped to a fake image of the same modality by a decoder. For the second image, an image-to-image translator is used to create a fake image of modality 2 with the help of information gained from the previously generated fake image of modality 1. In the supervised training, the original paired images are provided, and so for each pair of fake images generated by the generator, the corresponding true pair in the original data can be found. Thus, the loss function of the supervised training is the pixel-wise reconstruction loss:
L 1 = E I 1 , I 2 p I 1 , I 2 I 1 I ^ 1 + I 2 I ^ 2
where I ^ 1 and I ^ 2 represent the fake images while I 1 and I 2 represent the real images and x x ^ represents the average Manhattan-distance calculated pixel-wise between intensities of images x and x ^ . In order to counter a severe overfitting issue present in the supervised learning model due to the small number of labeled images available, an unsupervised learning model is also deployed. In the unsupervised model, the generator is trained using unpaired images and noise vectors rather than encodings. The model aims to minimize the Wasserstein (W) distances between the fake and real images. As such, the loss function of the unsupervised generator is:
L u n s u p = W 1 + W 2
where W 1 is the W distance between the real and fake images of modality 1 and W 2 is the W distance between the real and fake images of modality 2. The semi-supervised training is carried out as follows: in one iteration, the generator is trained in a semi-supervised manner with a set of paired training images. In the next iteration, the decoder and image translator are trained in an unsupervised manner using a set of unpaired images. The supervised and unsupervised training are alternated for 40,000 iterations. The usage of supervised learning allows for the creation of correctly paired images, while unsupervised training allows for higher realism and diversity. The images generated by this model were used as real training data in a single-label prostate cancer classification task in which each pair of images is classified as clinically significant (CS) or non-clinically significant (nonCS). The classifier was trained using 483 synthetic multimodal images of type ADC-T2w, and the testing set contained 50 real CS images and 50 real nonCS images. The evaluation metrics used were the Inception Score (IS) and Fréchet Inception Distance (FID) for quality assessment of the synthetic images, Mutual Information Distance (MID) to assess the correctness of the pairings between ADC and T2w synthetic images, and finally, the accuracy of the classification task. When compared with four state-of-the-arts GAN-based image synthesis methods (Costa et al. [64], CoGAN [65], CycleGAN [66], and pix2pix [67]), the proposed model achieved better results in all metrics, with a classification accuracy of 93%.
A summary of the methods discussed in this section is provided in Table 2. Three adversarial methods have been reviewed: Uncertainty-guided VAT with BNN, SSAC, and Bi-Modality Image Synthesis. For each model, the table contains the specific adversarial approach (VAT-based or GAN-based), the datasets used for testing, as well as the best performance score obtained on each dataset.

4.3. Graph-Based Methods

GraphXNET: GraphXNET is a semi-supervised graph-based framework that performs a classification task with an extremely small number of labeled samples and a large number of unlabeled samples. The function used in this model is Δ 1 u = W D 1 u , based on the normalized graph p-Laplacian with p = 1. The algorithm is as follows: for each class k , the model assumes that there exists a set of labeled nodes I k 1 l . For each class k , a variable u k is chosen that has values in all nodes of the graph. Assume the total number of classes is L . For all unlabeled nodes i > l , the L chosen variables are coupled with the following constraint: k = 1 L u i k = 0 , i > l . Using a small positive number ϵ , the following constraint is applied:
u i k ϵ                   if   i I k u i k ϵ           if   i I k   and   k k
The goal of the model is then to minimize the normalized ratios k Δ 1 ( u k ) | u k | . This model was evaluated using the ChestX-ray14 dataset. The results of the classification task were compared to two previously existing state-of-the-art deep learning methods for X-ray classification proposed by Wang et al. [68] and Yao et al. [69]. For both methods, 70% of the data was labeled, whereas, for the GraphXNET method, only 20% was labeled. The metric used for comparison was the average of the AUC score over all 14 classes. GraphXNET achieved the highest average score (0.78) despite using a much lower number of labeled samples. Additionally, when tested with three different, random partitions of data, GraphXNET proved to be more stable to variations in the data partition.
Graph-Embedded Random Forest [70]: This method makes the standard random forest algorithm more reliable when dealing with a low number of labeled samples. With the traditional method, insufficient training data leads to a limited depth, an inaccurate prediction model for leaf nodes, and a sub-optimal splitting strategy [71]. In the graph-based model, Gu et al. [70] improved the splitting strategy by replacing the information gain algorithm with a graph-embedded entropy. The goal is to utilize the local structure of unlabeled data to achieve higher reliability when using a small number of labeled data while maintaining the advantages of the random forest, such as low computational burden and robustness towards overfitting. The loss function is a sum of the supervised loss and a graph Laplacian regularization term. First, a graph G = V ,   E ,   W is constructed using labeled and unlabeled data. The nodes represent the training samples and W is the symmetric weight matrix, calculated as follows:
W i j = e x p x i x j 2 2 σ 2   if   x i , x j   are   neighbors 0 if   x i , x j   are   not   neighbors
Based on the label information for unlabeled samples gained from the graph embedding, the new information gained is as follows:
G m w , τ , X l , Y l , X u = G m S S l G m S l + S u G m S u / S
where S is the node, S l is the left child node, S u is the right child node, τ is a threshold, X l and Y l represent labeled instances and the class labels, respectively, and X u represents unlabeled instances. The medical image datasets used for the classification task were the Digital Retinal Images for Vessel Extraction (DRIVE) 2D retinal vessel dataset and the Big Neuron 3D neuronal dataset. For both datasets, 40,000 samples were randomly selected from the training and testing sets, respectively, each containing 20,000 positive and 20,000 negative samples. The results were compared with the standard random forest trained on labeled data only, as well as the RobustNode method [71] trained on both labeled and unlabeled data. The results showed considerable improvement over both the standard random forest and the RobustNode method. For the DRIVE dataset, training with a labeled dataset size of 400 achieved an accuracy of 79.42%. For the Big Neuron dataset, training with 1500 labeled data achieved an accuracy of 74.16%.
Colorectal Cancer Tissue Classification Using Semi-Supervised Hypergraph Convolutional Network [72]: Bakht et al. [72] propose a hypergraph-based approach for colorectal cancer (CRC) classification. Hypergraphs allow more complex relationships between nodes compared to standard graphs by allowing one edge to join multiple nodes. The images used for classification were CRC Whole Slide Images (WSIs), high-resolution images captured from a microscope slide representing tissue structures that can be used to identify malignancy. First, patches of size 224 × 224 are extracted from the images, and a VGG-19 [73] model in feed-forward mode is used to extract the feature matrix X from the set of n patches. The hypergraph G V , E , W is then constructed from X by connecting each vertex with its k nearest neighbors. The hypergraph is represented by a vertex-edge probabilistic incidence matrix H of size n × n as follows:
h n , e = e x p d p m a x d avg ,       if   n e 0 , if   n   e
where d is the Euclidean distance between the current node and the neighbor node, d avg is the average Euclidean distance between the k neighbors, and p m a x is the maximum probability. The degrees of each vertex v V and edge e E were calculated as follows:
d ( v ) = v V h ( v , e ) , d ( e ) = e E h ( v , e )
The combination of all node and edge degrees results in the diagonal matrices D v and D e . In the classification step, X and H are fed to a hypergraph neural network (HGNN) consisting of three hidden convolutional layers and a SoftMax classification layer. Spectral graph convolution is used for representation learning as follows:
X L + 1 = σ D v 1 / 2 H W D e 1 H T D v 1 / 2 X L θ L
where σ is the activation function, θ is a parameter to be learned during training, and W is one’s diagonal matrix. Layer L outputs X L + 1 which is received as an input by layer L + 1. The model was tested on a dataset of CRC WSIs containing 4995 patches of size 224 × 224 representing 7 different tissue types. 70% of the dataset was used for training and 30% for testing. The model achieved an F1 score of 0.94 for the seven-class classification task.
A summary of the methods discussed in this section is provided in Table 3. Three graph-based methods have been reviewed: GraphXNET, Graph-Embedded RF, and CRC with HGNN. For each model, the table contains the datasets used for testing as well as the best performance score obtained on each dataset.

4.4. Other Methods

The following section provides reviews of methods that do not neatly fit into the previous three categories.
Anti-curriculum Pseudo-labeling for Semi-supervised Medical Image Classification (ACPL) [37]: Liu et al. [37] propose a pseudo-labeling based image classification method for the Chest X-ray14 and ISIC2018 Skin Lesion Analysis datasets. The ACPL approach aims to address the disadvantages of traditional pseudo-labeling and develop a model capable of achieving state-of-the-art results on par with consistency-based methods. ACPL relies on the existence of a distribution shift between labeled and unlabeled data. The unlabeled samples chosen for pseudo-labeling are the ones located as far as possible from the labeled data’s distribution, making them more likely to belong to the minority class and, thus, more informative, resulting in a more balanced training process. The informativeness of a sample is determined using the cross-distribution sample informativeness (CDSI) measure, which computes the closeness of the unlabeled sample to the set of most informative labeled samples called the anchor set ( D A ). The CDSI was calculated as follows:
h f θ x , D A = 1 ,           p γ ζ = h i g h | x , D A > τ 0 ,           otherwise
in which ζ denotes the information content random variable (low, medium, high), γ represents the parameters of the Gaussian Mixture Model (GMM) and τ = m a x p γ ζ =   low   | x , D A , p γ ζ =   medium   | x , D A . Once the set of most informative unlabeled samples is selected, the informative mixup (IM) method is used for pseudo-labeling. This method mixes the labels obtained from the K-nearest neighbor (KNN) classification with the labels obtained from the model p θ . , where p θ x = σ f θ x wherein f θ x is the input image feature and σ . is the final activation function to calculate an output in [ 0 , 1 ] Y . The IM method carries out the pseudo-labeling process by calculating a linear combination of the model prediction p θ x and the KNN prediction is weighted by the density score. Following the pseudo-labeling stage, the Anchor Set Purification (ASP) algorithm is used to select the most informative pseudo-labeled samples to be added to the anchor set. The backbone model used for experimentation was DenseNet-121 [23] for both datasets. The images in the X-ray14 dataset were resized to 512 × 512. The training was performed with batch size 16 and a learning rate of 0.05 for 20 epochs in the first step and then 50 epochs, where the anchor set was updated with new pseudo-labeled samples every 10 epochs. The model achieved AUC scores of 74.69, 79.96, 79.90, 80.31, and 81.06 for labeled dataset sizes of 2%, 5%, 10%, 15%, and 20%, respectively. For the ISIC2018 dataset, images were resized to 224 × 224. The training was conducted with batch size 32 and a learning rate of 0.001 for 40 epochs in the warm-up stage, followed by 100 epochs, with the anchor set being updated with ASP every 20 epochs. The model achieved AUC and F1-scores of 94.36 and 62.23 at a labeled set size of 20%.
Federated Semi-supervised Medical Image Classification via Inter-client Relation Matching [74]: Liu et al. [74] propose an SSL method based on FL. Existing SSL methods cannot be used reliably for FL as they depend on the labeled dataset being accessible during training. However, in the case of FL, the local data at a particular location may be entirely unlabeled. As such, Liu et al. [74] build an interaction between the learning of labeled and unlabeled data by leveraging inherent disease relationships which are independent of a specific local dataset. This client-invariant disease relation information can be extracted from the labeled data and used to supervise the learning of unlabeled data to ensure it captures similar disease relationships. An uncertainty-based algorithm is used to filter out unreliable pseudo-labels. The backbone FL framework of this method follows the standard FL paradigm in which the central server sends the global parameters θ to every client k. The model is trained at each k with a local objective function L k and the updated local parameters are sent back to the central server and aggregated to update the global model. In the proposed model, local parameters are aggregated using weights proportional to the size of their database following the federated averaging algorithm (FedAvg) method [75]. The local learning objectives use the cross-entropy loss at labeled clients and the consistency regularization mechanism at unlabeled clients with transformations such as rotation, translation, and flip. The inter-client relation matching (IRM) method is introduced for disease relation estimation to assist learning at unlabeled clients, in which the disease relation matrix M l is derived from the data’s pre-SoftMax features at each labeled client and averaged to obtain the matrix M representing the general release relation. The model was tested on two datasets; ISIC2018 and the RSNA 2019 Brain CT dataset. From the RSNA dataset, 25,000 slices were sampled and divided into 70%, 10%, and 20% splits for training, validation, and testing. The ISIC2018 images were resized to 224 × 224 and split similarly. The FL setting was simulated by partitioning the training set into 10 local clients, of which 20% (two clients) were labeled. DenseNet121 was used as the backbone model. On the RSNA dataset, the model achieved an AUC score of 87.56 and F1-score of 59.86. On the ISIC2018 dataset, it achieved an AUC score of 92.46 and F1-score of 55.81.
Semi-HIC [76]: Su et al. [76] propose an SSL method called Semi-HIC for the classification of histopathological images, which are representations of biopsy samples captured by a microscope. This model used a patch-wise CNN to learn embeddings of labeled and unlabeled images, wherein three convolutional blocks are used to extract low-level detail information, and four cascaded inception blocks are used to create discriminative representations for patches. Semi-HIC uses a novel loss function obtained from combining an association cycle consistency (ACC) loss and a maximal conditional association (MCA) loss. These losses are intended to address the two major challenges in histopathological image classification, which are inter-class similarity and intra-class variation. They are based on the walker loss and visit loss proposed by Haeusser et al. [39] in their LA method. The original visit loss is vulnerable to the inter-class similarity in histopathological images, as a cancerous labeled embedding E i l and unlabeled embedding E k u with the underlying non-cancerous class may have high similarity. Su et al.’s MCA loss aims to address this by separately calculating the visit probability of non-cancerous embeddings p k v n and the visit probability of cancerous embeddings p k v c . If the unlabeled embedding E k u has an underlying cancerous class, it is expected to have a high p k v c and low p k v n and vice versa. Thus, the MCA loss L M C A is defined as the conditional entropy between p k v n and p k v c :
L M C A = 1 M k = 1 M p k v c l o g p k v c + p k v n l o g p k v n
where M is the size of E u . The goal is to maximize the difference between p k v c and p k v n for the model to learn discriminative representations for unlabeled patches. The ACC loss, on the other hand, applies modifications to the walker loss in order to make the model robust to intra-class variation. This is achieved by penalizing association cycles that start and end at labeled embeddings belonging to the same image. The ACC loss was calculated as follows:
L A C C = 1 N i = 1 N j = 1 N φ h i , h j P i j T l o g P i j l u l
representing the cross-entropy between the association cycle probabilities P l u l and uniform target distribution P T . The penalty factor φ h i , h j is 1 if the association cycle does not start and end at the same image ( h i h j ) and 5 if it does. Finally, the total loss function of the model is:
L total   = λ L B C E + γ L A C C + β L M C A
wherein L B C E is the binary cross-entropy loss of the CNN model. λ , γ , and β are the weights of each loss which are set to 1, 1, and 0.75, respectively. This model was evaluated on the Bioimaging2015 and the Grand Challenge on BreAst Cancer Histology (BACH) datasets. The Bioimaging2015 dataset contains 249 training and 36 testing images of cancerous and non-cancerous histopathological images. The training set is split into 222 training and 27 validation images. The BACH dataset contains 400 cancerous and non-cancerous histopathological images split with the ratio 70:10:20 into the training, validation, and testing sets. After patch extraction, the size of the dataset is 8420, 1163, and 2369 patches for training, validation, and testing, respectively. 5-fold cross-validation was carried out on the BACH dataset, while the average of 5 runs was computed for the Bioimaging2015 dataset. For the BACH dataset, the model obtained error rates of 30.38, 24.87, and 15.46 for 20, 40, and 80 labeled images per class, respectively. For the Bioimaging2015 dataset, the model achieved error rates of 30.13, 28.39, and 25.23 for 20, 40, and 80 labeled images per class, respectively. The model proved to be more successful than state-of-the-art supervised models for histopathology image classification.
A summary of the methods discussed in this section is provided in Table 4. Three methods have been reviewed: ACPL, FL with IRM, and Semi-HIC. For each model, the table contains the type of approach to SSL, the datasets used for testing, as well as the best performance score obtained on each dataset.

4.5. Hybrid Methods

This section covers a few models which combine two or more of the previously reviewed approaches.
Local and Global Consistency Regularized Method [77]: This method is both consistency-based and graph-based. This method uses the Mean Teacher framework and also enforces the local and global consistency of the data [78]. Local consistency is the tendency of instances from the same class to be in the same area in the feature space, and global consistency is the tendency of instances from the same global structure to have the same label. This method integrates label propagation (LP) to enforce local and global consistency. LP is a semi-supervised learning algorithm that involves the propagation of labels from labeled samples to unlabeled samples based on their closeness, which is determined by the affinity matrix. For an unlabeled instance x , the label is calculated by taking the weighted average of the labeled instances close to x ; then, the new label of x can be propagated to other neighboring unlabeled data. Next, a graph is constructed based on the labels created by the LP algorithm as well as the ground truth labels:
A i j = 1 ,           if   y i = y j 0 ,           otherwise
in which y i and y j represent the labels of the data. Contrastive Siamese loss [79] is used in order to enforce local and global consistencies by pulling the instances from the same class closer and pushing the ones from different classes farther apart:
L s = z i z j 2 , if   A i j = 1 m a x 0 , m z i z j 2 ,         if   A i j = 0
where m is a hyperparameter and z is the feature vector extracted from the intermediate layers of the student network. The final loss function is as follows:
L total = Loss m t + w ( τ ) λ g 1 x i , x j X l L s 1 + λ g 2 x i X l , x j X u L s 2
where the total loss is a weighted addition of the Mean Teacher loss ( Loss m t ) and the graph-based loss ( L s 1 and L s 2 ). λ g 1 is the weight of the loss calculated on the labeled instances and λ g 2 is the weight of the loss calculated on both labeled and unlabeled instances. The loss on the unlabeled instances is not calculated separately as the labels predicted by the LP algorithm are noisy, and their inclusion may harm the performance of the method. Testing was conducted using two datasets, the Multi-organ Nucleus Segmentation (MoNuSeg) dataset containing 22,462 samples and the Ki-67 nucleus dataset containing 17,516 samples. In both datasets, each sample is labeled with one of four types, leading to a single labeled classification task. For both datasets, 80% of the data was used for training and 20% for testing. From the training data, 20% was used for validation. The model architecture is a 13-layer convolutional neural network similar to the one proposed in the original Mean Teacher method [27]. The metric used to gauge the method’s performance is F1. Testing was conducted at different percentages of labeled data (5, 10, 25, 50, 100 for MoNuSeg and 1, 5, 10, 100 for Ki-67). For the MoNuSeg dataset, the proposed method achieved higher scores than both the supervised baseline and the Mean Teacher for all percentages, notably in the 5% and 10% percentages, where it obtained scores approximately 2% higher than the Mean Teacher. The highest score obtained was 76.89 for 50% of labeled data, and the lowest was 75.02 for 5% of labeled data. In the case of the Ki-67 dataset, the proposed method also achieved a higher score than the supervised baseline and Mean Teacher, the highest increase from Mean Teacher being approximately 2.5% at the 5% labeled percentage. The highest score was 79.79 for 10% of labeled data, and the lowest was 74.9 for 1% of labeled data.
Deep virtual adversarial self-training with consistency regularization [80]: This model is both adversarial and consistency-based. It consists of three parts: self-training, consistency regularization, and virtual adversarial training. As previously mentioned, self-training involves using the model itself to generate labels for unlabeled samples, which will then become part of the labeled training set in the next iteration. The model outputs a probability distribution over all classes, and the generated label is only kept when the highest probability is above a predetermined threshold. Consistency enforcement is applied to both labeled and unlabeled samples. For the labeled samples, weak augmentation is applied, and the generated label is encouraged to be consistent with the ground truth label. For unlabeled samples, weak augmentation is applied, then a pseudo-label is generated by the self-training model; the same input image is then strongly augmented, and the new prediction must be consistent with the previously obtained pseudo-label. Lastly, virtual adversarial training is applied to improve the robustness of the model and strengthen its generalization capability. Thus, the loss function of the model becomes the weighted sum of three losses: the supervised cross-entropy loss s for labeled data, the regularization loss r for unlabeled data and the virtual adversarial training loss vat applied to both labeled and unlabeled data.
L = s + α · r + β · vat
where α and β are the weighting coefficients. Experimentation was performed using two datasets; a breast ultrasound (BRUS) dataset and an optical coherence tomography (OCT) dataset. The BRUS dataset contains 39,904 images, 22,026 of which are labeled as malignant (14,557) or benign (7469). The other 17,878 images are unlabeled. The labeled data were randomly split into 34% training (7662), 34% validation (7579), and 32% testing (6785). The OCT dataset contains 109,309 images, each labeled with four different types (three diseases and one normal). There were 20,000 images taken for the validation and testing datasets, respectively, and 1% of the remaining data (691 images) was taken for the labeled dataset, whereas the rest was considered unlabeled. The network architecture used for this experimentation is DenseNet-121. The weak augmentations include horizontal and/or vertical flips with a 50% probability, as well as randomly applied horizontal and/or vertical translations by less than 20% of the image width or height. Strong augmentation includes transformations such as color channel reduction, contrast maximization, solarization, rotation, and posterization. The results are compared with Mean Teacher, VAT, MixMatch [42], and FixMatch [44], as well as the fully supervised DenseNet-121 as the baseline. The metrics calculated were the accuracy, F1, AUC, and Kappa scores. The proposed method achieved higher scores than all the compared methods on both datasets, with an F1-score of 0.88 on the BRUS dataset and a Macro F1-score of 0.91 on the OCT dataset.
Pseudo-labeling generative adversarial networks (PLGAN) [81]: This model proposed by Mao et al. [81] is based on pseudo-labeling and GANs. It additionally incorporates CL and MixMatch. The training process of this method is comprised of 4 steps; (1) pretraining, (2) image generation, (3) finetuning, and (4) pseudo-labeling. In step 1, the feature layer of ResNet50 [59] is pre-trained using CL to extract key image features. In step 2, GANs are used to generate images that simulate the real distribution of the labeled images with random Gaussian noise. In step 3, the cross-entropy loss is used to classify the images created by the generator in order to finetune the discriminator for classification. The model contains two classifiers, a global and a local classifier, both of which use two convolutional blocks to extract feature information. In step 4, the MixMatch approach is used for pseudo-labeling. The trained generator is used to create additional unlabeled samples in order to expand the unlabeled dataset, and these generated images are combined with the original dataset. In the same fashion, the pseudo-labels created for both the original and generated samples are combined to create the complete set of pseudo-labels. The overall loss function of the model contains four loss functions for each of the 4 steps. The loss function for the first step is the infoNCE [82] loss for CL and the reconstruction loss for intermediate vectors to avoid CL pattern collapse. For the second step, the loss function is the least squares loss. The loss of the semi-supervised finetuning step is the addition of the supervised cross-entropy loss and the unsupervised loss. Lastly, the fourth step uses the loss function of MixMatch. The model was tested on a dataset of optical coherence tomography (OCT) images for the classification of retinal degeneration, with each sample belonging to one of four categories (three disease labels and one normal label). Additionally, it was tested on a COVID-19 dataset, a brain tumor MRI dataset, and a chest x-ray dataset. For the OCT and chest x-ray datasets, 100 labeled and 1000 unlabeled images were used. In the MRI dataset, 100 labeled images and 480 unlabeled images were selected. Lastly, 200 labeled and 200 unlabeled images were used for the COVID-19 dataset. The results achieved by the model were classification accuracies of 87.06, 89.50%, 80.50%, and 96.80% for the OCT, COVID-19, MRI, and x-ray datasets, respectively. In addition to PLGAN, the authors developed PLGAN+, which uses Comatch [45] for pseudo-labeling. This model achieved accuracy scores of 97.10%, 89.30%, 86.80%, and 97.50% for the OCT, COVID-19, MRI, and x-ray datasets, respectively. These results were compared to the results gained by several other state-of-the-art SSL models, including Comatch, MixMatch, FixMatch, and VAT. Overall, PLGAN+ achieved more favorable results than all other models.
Semi-supervised medical image classification based on CamMix [83]: Guo et al. [83] proposed an SSL model for medical image classification similar to MixMatch, which combines several SSL approaches. The consistency-based approach is used on unlabeled data to generate pseudo-labels that are robust to various augmentations. The authors argue that the MixUp method, which mixes samples by linearly interpolating the input samples and labels, results in unnatural mixed samples. Thus, they proposed CamMix, a novel MSDA method that mixes pairs of input samples and labels based on the class activation mask generated from the predictions of both labeled and unlabeled samples. As in MixMatch, entropy minimization is achieved by sharpening the target distribution for unlabeled data. The class activation map is obtained as follows for each batch b at each epoch:
GradMaxCam b = max ReLU k w k b A k
where w k b is the weight of feature map k for batch b and A. w k b is calculated as follows:
w k b = 1 Z i j Y b A i j k
where A i j k is the pixel value at location (i, j) of k, Z represents the number of pixels in k and Y b denotes the maximum prediction score of batch b generated by the classification model. Random threshold λ 0 , 1 is applied to the gray-level G r a d M a x C a m b to obtain the binary mask CamMask. The pixels higher than 1 − λ are set to 1, and all others are set to 0. The CamMix algorithm takes a batch of labeled and unlabeled data and their corresponding predictions, including weak and strong augmentations of the data, and generates a mixed batch of original samples and shuffled samples, whose corresponding labels are mixed based on the number of pixels in CamMask. Thus, considering the original samples i n p u t 1 and the shuffled samples i n p u t 2 = i n p u t r a n d _ i n d e x as well as the corresponding label targets t a r g e t 1 and t a r g e t 2 , the CamMask is obtained by calculating GradMaxCam i n p u t 1   , 1 λ , and the parameter l a m is calculated based on the number of pixels in the CamMask as follows:
l a m = s u m CamMask   = = 1 / CamMask . size 0 CamMask . size 1
The mixed batch mixed _ input is obtained from i n p u t 1 and i n p u t 2 by calculating the following:
m i x e d i n p u t = i n p u t 1 CamMask + i n p u t 2 1 CamMask
Lastly, the overall loss of the model is defined as:
l = c r i t e r i o n ( logits ,   t a r g e t 1 )   *   lam + c r i t e r i o n ( logits ,   t a r g e t 2 )   *   ( 1   lam )
wherein logits = model(mixed_input). The model was tested with a ResNet-18 on two datasets: The Interstitial Lung Disease (ILD) dataset and the ISIC2018 dataset. Instead of the F1-score, the metric used for evaluation was F a v g , the average F1-score over the different classes. The ILD dataset contains 109 CT scans with an average of 25 slices per case for a total of 2795 slices. All regions of interest were cropped into patches of size 32 × 32, resulting in approximately 7000 labeled patches and 19,000 unlabeled patches. Six typical lung patterns were used for the classification task. By using CamMix, the model achieved an F a v g score of 95.34, higher than previous state-of-the-art MSDA methods for the ILD dataset, such as MixUp, CutMix [84], and Fmix [85]. For the ISIC2018 dataset, the model was trained with 20% labeled and 80% unlabeled data, obtaining an AUC score of 94.04 and F a v g of 78.15. These scores were higher than the scores obtained by the previously mentioned MSDA methods. The SRC-MT model was also used for comparison and achieved a F a v g score of 78.15, higher than the one obtained by CamMix. On the other hand, it obtained a lower AUC score.
A summary of methods discussed in this section is provided in Table 5. Three hybrid methods have been reviewed: Local and Global Consistency Regularized, PLGAN, SSL with CamMix, and Wan et al.’s method. We have reported the SSL methods which have been combined to create each hybrid method, as well as the datasets used for testing each model and the best performance score obtained on each dataset.

5. Summary of Results

Table 6 shows a comparison of the models which were tested on the Chest X-ray14 dataset for multi-label classification and/or the ISIC2018 for multi-class classification. For both datasets, all models used 20% of the data as the labeled dataset and 80% as the unlabeled dataset. Among them, the highest score on the Chest X-ray14 dataset was achieved by the S2MTS2 method (AUC of 82.51), followed by ACPL (AUC of 81.77). On the ISIC2018 dataset, the highest score was achieved by Uncertainty-guided VAT with BNN with an F1 score of 69.67 and AUC of 96.04, considerably higher than all other models. It is worth noting that ACPL obtained relatively high scores on both datasets, showing that pseudo-labeling can indeed achieve state-of-the-art results by improving the measure used for the selection of informative pseudo-labels.
Figure 3 showcases the results obtained by different methods for each dataset.
Among the GAN-based methods, the bi-modality image synthesis model is notable as it allows a more comprehensive look at a certain medical problem by combining different types of scans of the same area, leading to a more successful diagnosis. Despite the usefulness of bi-modality and multi-modality medical image diagnosis, labeled data in this category is even more scarce, making the semi-supervised image synthesis method to create more labeled inputs a promising method.
Hybrid models which use a combination of different SSL approaches are popular in recent research and show much promise. Most of these models incorporate consistency regularization, including the MixMatch method, which has gained considerable attention. Hybrid models inspired by or incorporating MixMatch, as well as DSNA approaches such as MixUp, are common. Among the hybrid methods, the consistency-based adversarial method proposed by Wang et al. achieved impressive results with a labeled dataset of only 691 images. However, this study assumes that the labeled and unlabeled data have the same distribution. Therefore, further experimentation with differently distributed data is needed.

6. Gaps and Future Work

As one would expect, the performance of the models varied with different percentages of labeled data. Some methods, such as SRC-MT and NoTeacher, revealed up to a 10% difference in score from 1% to 10% labeled data. GraphXNET showed an increase in more than 20%, with an AUC score of 0.53 and 0.78 for 2% and 20% labeled datasets, respectively. Other methods had smaller differences: the SSAC model’s AUC score increased by 3% from 20% to 100% labeled data, and the bi-modality synthesis model increased by 4% from 1% to 10% labeled data. Further research can be conducted on accomplishing higher robustness to changes in the percentage of labeled data and achieving successful training on lower numbers of labeled data.
As mentioned earlier, the areas of bi and multi-modality image synthesis is an areas that can be studied further. Additionally, given the promise shown by the ACPL method, it may be beneficial to explore other strategies for measuring confidence and uncertainty in order to achieve improved results with the traditional pseudo-labeling method. Another area of learning that may be worth exploring is active learning. In this method, instead of using an existing dataset, a few valuable unlabeled samples are carefully selected, and a medical expert is asked to label them. As the samples to be labeled are chosen specifically based on the information they offer, this method results in a smaller amount of labeled data needed and, therefore, less time and effort in labeling from experts.
Lastly, more thorough experimentation is needed to gauge the performance of these models in a real-world clinical setting and with larger patient cohorts.

Author Contributions

Conceptualization: Z.S. and I.Z.; investigation: Z.S.; writing—original draft preparation: Z.S.; writing—review and editing: I.Z. and Z.S.; supervision: I.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sidey-Gibbons, J.A.M.; Sidey-Gibbons, C.J. Machine learning in medicine: A practical introduction. BMC Med. Res. Methodol. 2019, 19, 64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ker, J.; Wang, L.; Rao, J.; Lim, T. Deep Learning Applications in Medical Image Analysis. IEEE Access 2018, 6, 9375–9389. [Google Scholar] [CrossRef]
  3. AlAmir, M.; AlGhamdi, M. The Role of Generative Adversarial Network in Medical Image Analysis: An in-depth survey. ACM Comput. Surv. 2022. [Google Scholar] [CrossRef]
  4. Kazeminia, S.; Baur, C.; Kuijper, A.; van Ginneken, B.; Navab, N.; Albarqouni, S.; Mukhopadhyay, A. GANs for medical image analysis. Artif. Intell. Med. 2020, 109, 101938. [Google Scholar] [CrossRef]
  5. Wang, W.; Liang, D.; Chen, Q.; Iwamoto, Y.; Han, X.-H.; Zhang, Q.; Hu, H.; Lin, L.; Chen, Y.-W. Medical Image Classification Using Deep Learning. In Deep Learning in Healthcare: Paradigms and Applications; Springer International Publishing: Cham, Switzerland, 2020; pp. 33–51. [Google Scholar]
  6. Swati, Z.N.K.; Zhao, Q.; Kabir, M.; Ali, F.; Ali, Z.; Ahmed, S.; Lu, J. Brain tumor classification for MR images using transfer learning and fine-tuning. Comput. Med. Imaging Graph. 2019, 75, 34–46. [Google Scholar] [CrossRef]
  7. O’Mahony, N.; Campbell, S.; Carvalho, A.; Harapanahalli, S.; Hernandez, G.V.; Krpalkova, L.; Riordan, D.; Walsh, J. Deep Learning vs. Traditional Computer Vision; Springer: Cham, Switzerland, 2020; pp. 128–144. [Google Scholar]
  8. Loussaief, S.; Abdelkrim, A. Deep learning vs. bag of features in machine learning for image classification. In Proceedings of the 2018 International Conference on Advanced Systems and Electric Technologies (IC_ASET), Hammamet, Tunisia, 22–25 March 2018; pp. 6–10. [Google Scholar]
  9. Boumaraf, S.; Liu, X.; Wan, Y.; Zheng, Z.; Ferkous, C.; Ma, X.; Li, Z.; Bardou, D. Conventional Machine Learning versus Deep Learning for Magnification Dependent Histopathological Breast Cancer Image Classification: A Comparative Study with Visual Explanation. Diagnostics 2021, 11, 528. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, Z.; Tang, C.; Sima, X.; Zhang, L. Research on Application of Deep Learning Algorithm in Image Classification. In Proceedings of the 2021 IEEE Asia-Pacific Conference on Image Processing, Electronics and Computers (IPEC), Dalian, China, 14–16 April 2021; pp. 1122–1125. [Google Scholar]
  11. Liu, P.; Choo, K.-K.R.; Wang, L.; Huang, F. SVM or deep learning? A comparative study on remote sensing image classification. Soft Comput. 2017, 21, 7053–7065. [Google Scholar] [CrossRef]
  12. Wang, P.; Fan, E.; Wang, P. Comparative analysis of image classification algorithms based on traditional machine learning and deep learning. Pattern Recognit. Lett. 2021, 141, 61–67. [Google Scholar] [CrossRef]
  13. Devi, M.R.S.S.; Kumar, V.R.V.; Sivakumar, P. A Review of image Classification and Object Detection on Machine learning and Deep Learning Techniques. In Proceedings of the 2021 5th International Conference on Electronics, Communication and Aerospace Technology (ICECA), Coimbatore, India, 2–4 December 2021; pp. 1–8. [Google Scholar]
  14. Miranda, E.; Aryuni, M.; Irwansyah, E. A survey of medical image classification techniques. In Proceedings of the 2016 International Conference on Information Management and Technology (ICIMTech), Jakarta, Indonesia, 16–18 November 2016; pp. 56–61. [Google Scholar]
  15. Withers, P.J.; Bouman, C.; Carmignato, S.; Cnudde, V.; Grimaldi, D.; Hagen, C.K.; Maire, E.; Manley, M.; Du Plessis, A.; Stock, S.R. X-ray computed tomography. Nat. Rev. Methods Prim. 2021, 1, 18. [Google Scholar] [CrossRef]
  16. Grover, V.P.B.; Tognarelli, J.M.; Crossey, M.M.E.; Cox, I.J.; Taylor-Robinson, S.D.; McPhail, M.J.W. Magnetic Resonance Imaging: Principles and Techniques: Lessons for Clinicians. J. Clin. Exp. Hepatol. 2015, 5, 246–255. [Google Scholar] [CrossRef]
  17. Köse, G.; Darguzyte, M.; Kiessling, F. Molecular Ultrasound Imaging. Nanomaterials 2020, 10, 1935. [Google Scholar] [CrossRef] [PubMed]
  18. Litjens, G.; Kooi, T.; Bejnordi, B.E.; Setio, A.A.A.; Ciompi, F.; Ghafoorian, M.; Van Der Laak, J.A.; Van Ginneken, B.; Sánchez, C.I. A survey on deep learning in medical image analysis. Med. Image Anal. 2017, 42, 60–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Liu, Q.; Yu, L.; Luo, L.; Dou, Q.; Heng, P.A. Semi-Supervised Medical Image Classification with Relation-Driven Self-Ensembling Model. IEEE Trans. Med. Imaging 2020, 39, 3429–3440. [Google Scholar] [CrossRef]
  20. van Engelen, J.E.; Hoos, H.H. A survey on semi-supervised learning. Mach. Learn. 2020, 109, 373–440. [Google Scholar] [CrossRef] [Green Version]
  21. Chapelle, O.; Scholkopf, B.; Zien, E.A. Semi-Supervised Learning (Chapelle, O. et al., Eds.; 2006) [Book reviews]. IEEE Trans. Neural Netw. 2009, 20, 542. [Google Scholar] [CrossRef]
  22. Chong, Y.; Ding, Y.; Yan, Q.; Pan, S. Graph-based semi-supervised learning: A review. Neurocomputing 2020, 408, 216–230. [Google Scholar] [CrossRef]
  23. Huang, G.; Liu, Z.; Weinberger, K. Densely Connected Convolutional Networks. arXiv 2016, arXiv:1608.06993. [Google Scholar]
  24. Laine, S.; Aila, T. Temporal ensembling for semi-supervised learning. arXiv 2016, arXiv:1610.02242. [Google Scholar]
  25. Maqsood, I.; Khan, M.R.; Abraham, A. An ensemble of neural networks for weather forecasting. Neural Comput. Appl. 2004, 13, 112–122. [Google Scholar] [CrossRef]
  26. Alam, K.M.R.; Siddique, N.; Adeli, H. A dynamic ensemble learning algorithm for neural networks. Neural Comput. Appl. 2020, 32, 8675–8690. [Google Scholar] [CrossRef] [Green Version]
  27. Tarvainen, A.; Valpola, H. Mean teachers are better role models: Weight-averaged consistency targets improve semi-supervised deep learning results. arXiv 2017, arXiv:1703.01780. [Google Scholar]
  28. Allen-Zhu, Z.; Li, Y. Feature Purification: How Adversarial Training Performs Robust Deep Learning. In Proceedings of the 2021 IEEE 62nd Annual Symposium on Foundations of Computer Science (FOCS), Denver, CO, USA, 7–10 February 2022; pp. 977–988. [Google Scholar]
  29. Miyato, T.; Maeda, S.; Koyama, M.; Ishii, S. Virtual Adversarial Training: A Regularization Method for Supervised and Semi-Supervised Learning. IEEE Trans. Pattern Anal. Mach. Intell. 2019, 41, 1979–1993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Pan, Z.; Yu, W.; Yi, X.; Khan, A.; Yuan, F.; Zheng, Y. Recent Progress on Generative Adversarial Networks (GANs): A Survey. IEEE Access 2019, 7, 36322–36333. [Google Scholar] [CrossRef]
  31. Goodfellow, I.; Pouget-Abadie, J.; Mirza, M.; Xu, B.; Warde-Farley, D.; Ozair, S.; Courville, A.; Bengio, Y. Generative adversarial networks. Commun. ACM 2020, 63, 139–144. [Google Scholar] [CrossRef]
  32. Radford, A.; Metz, L.; Chintala, S. Unsupervised representation learning with deep convolutional generative adversarial networks. arXiv 2015, arXiv:1511.06434. [Google Scholar]
  33. Odena, A.; Olah, C.; Shlens, J. Conditional image synthesis with auxiliary classifier gans. In Proceedings of the International Conference on Machine Learning, Sydney, Australia, 6 August 2017; pp. 2642–2651. [Google Scholar]
  34. Aviles-Rivero, A.I.; Papadakis, N.; Li, R.; Sellars, P.; Fan, Q.; Tan, R.T.; Schönlieb, C.-B. GraphX $$^\mathbf {\small NET}-$$ Chest X-ray Classification under Extreme Minimal Supervision. In Proceedings of the International Conference on Medical Image Computing and Computer-Assisted Intervention, Shenzhen, China, 13–17 October 2019; pp. 504–512. [Google Scholar]
  35. Von Luxburg, U. A tutorial on spectral clustering. Stat. Comput. 2007, 17, 395–416. [Google Scholar] [CrossRef]
  36. Fu, S.; Liu, W.; Zhang, K.; Zhou, Y.; Tao, D. Semi-supervised classification by graph p-Laplacian convolutional networks. Inf. Sci. 2021, 560, 92–106. [Google Scholar] [CrossRef]
  37. Liu, F.; Tian, Y.; Chen, Y.; Liu, Y.; Belagiannis, V.; Carneiro, G. ACPL: Anti-Curriculum Pseudo-Labelling for Semi-Supervised Medical Image Classification. In Proceedings of the Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, New Orleans, LA, USA, 21–24 June 2022; pp. 20697–20706. [Google Scholar]
  38. Kaissis, G.A.; Makowski, M.R.; Rückert, D.; Braren, R.F. Secure, privacy-preserving and federated machine learning in medical imaging. Nat. Mach. Intell. 2020, 2, 305–311. [Google Scholar] [CrossRef]
  39. Haeusser, P.; Mordvintsev, A.; Cremers, D. Learning by association—A versatile semi-supervised training method for neural networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 89–98. [Google Scholar]
  40. Jaiswal, A.; Babu, A.R.; Zadeh, M.Z.; Banerjee, D.; Makedon, F. A survey on contrastive self-supervised learning. Technologies 2020, 9, 2. [Google Scholar] [CrossRef]
  41. Cai, Q.; Wang, Y.; Pan, Y.; Yao, T.; Mei, T. Joint contrastive learning with infinite possibilities. Adv. Neural Inf. Process. Syst. 2020, 33, 12638–12648. [Google Scholar]
  42. Berthelot, D.; Carlini, N.; Goodfellow, I.; Papernot, N.; Oliver, A.; Raffel, C.A. Mixmatch: A holistic approach to semi-supervised learning. arXiv 2019, arXiv:1905.02249. [Google Scholar]
  43. Zhang, H.; Cisse, M.; Dauphin, Y.N.; Lopez-Paz, D. mixup: Beyond empirical risk minimization. arXiv 2017, arXiv:1710.09412. [Google Scholar]
  44. Sohn, K.; Berthelot, D.; Carlini, N.; Zhang, Z.; Zhang, H.; Raffel, C.A.; Cubuk, E.D.; Kurakin, A.; Li, C.-L. Fixmatch: Simplifying semi-supervised learning with consistency and confidence. Adv. Neural Inf. Process. Syst. 2020, 33, 596–608. [Google Scholar]
  45. Li, J.; Xiong, C.; Hoi, S.C. Comatch: Semi-supervised learning with contrastive graph regularization. In Proceedings of the IEEE/CVF International Conference on Computer Vision, Montreal, QC, Canada, 11–17 October 2021; pp. 9475–9484. [Google Scholar]
  46. Cheplygina, V.; de Bruijne, M.; Pluim, J.P. Not-so-supervised: A survey of semi-supervised, multi-instance, and transfer learning in medical image analysis. Med. Image Anal. 2019, 54, 280–296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Liu, Y.; Cao, J.; Li, B.; Yuan, C.; Hu, W.; Li, Y.; Duan, Y. Knowledge Distillation via Instance Relationship Graph. In Proceedings of the 2019 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR), Long Beach, CA, USA, 15–20 June 2019; pp. 7089–7097. [Google Scholar]
  48. Battaglia, P.W.; Hamrick, J.B.; Bapst, V.; Sanchez-Gonzalez, A.; Zambaldi, V.; Malinowski, M.; Tacchetti, A.; Raposo, D.; Santoro, A.; Faulkner, R. Relational inductive biases, deep learning, and graph networks. arXiv 2018, arXiv:1806.01261. [Google Scholar]
  49. Deng, J.; Dong, W.; Socher, R.; Li, L.-J.; Li, K.; Fei-Fei, L. Imagenet: A large-scale hierarchical image database. In Proceedings of the 2009 IEEE Conference on Computer Vision and Pattern Recognition, Miami, FL, USA, 20–25 June 2009; pp. 248–255. [Google Scholar]
  50. Bai, W.; Oktay, O.; Sinclair, M.; Suzuki, H.; Rajchl, M.; Tarroni, G.; Glocker, B.; King, A.; Matthews, P.M.; Rueckert, D. Semi-supervised learning for network-based cardiac MR image segmentation. In Proceedings of the International Conference on Medical Image Computing and Computer-Assisted Intervention, Quebec City, QC, Canada, 11–13 September 2017; pp. 253–260. [Google Scholar]
  51. Diaz-Pinto, A.; Colomer, A.; Naranjo, V.; Morales, S.; Xu, Y.; Frangi, A.F. Retinal Image Synthesis and Semi-Supervised Learning for Glaucoma Assessment. IEEE Trans. Med. Imaging 2019, 38, 2211–2218. [Google Scholar] [CrossRef] [PubMed]
  52. Unnikrishnan, B.; Nguyen, C.; Balaram, S.; Li, C.; Foo, C.S.; Krishnaswamy, P. Semi-supervised classification of radiology images with NoTeacher: A teacher that is not mean. Med. Image Anal. 2021, 73, 102148. [Google Scholar] [CrossRef]
  53. Ke, Z.; Wang, D.; Yan, Q.; Ren, J.; Lau, R.W. Dual student: Breaking the limits of the teacher in semi-supervised learning. In Proceedings of the IEEE/CVF International Conference on Computer Vision, Seoul, Republic of Korea, 27 October–2 November 2019; pp. 6728–6736. [Google Scholar]
  54. Bien, N.; Rajpurkar, P.; Ball, R.L.; Irvin, J.; Park, A.; Jones, E.; Bereket, M.; Patel, B.N.; Yeom, K.W.; Shpanskaya, K. Deep-learning-assisted diagnosis for knee magnetic resonance imaging: Development and retrospective validation of MRNet. PLoS Med. 2018, 15, e1002699. [Google Scholar] [CrossRef] [Green Version]
  55. Liu, F.; Tian, Y.; Cordeiro, F.R.; Belagiannis, V.; Reid, I.; Carneiro, G. Self-Supervised Mean Teacher for Semi-Supervised Chest X-ray Classification; Springer: Cham, Switzerland, 2021; pp. 426–436. [Google Scholar]
  56. Liu, P.; Zheng, G. Handling Imbalanced Data: Uncertainty-Guided Virtual Adversarial Training with Batch Nuclear-Norm Optimization for Semi-Supervised Medical Image Classification. IEEE J. Biomed. Health Inform. 2022, 26, 2983–2994. [Google Scholar] [CrossRef]
  57. Cui, S.; Wang, S.; Zhuo, J.; Li, L.; Huang, Q.; Tian, Q. Towards discriminability and diversity: Batch nuclear-norm maximization under label insufficient situations. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, Seattle, DC, USA, 13–19 June 2020; pp. 3941–3950. [Google Scholar]
  58. Xie, Y.; Zhang, J.; Xia, Y. Semi-supervised adversarial model for benign–malignant lung nodule classification on chest CT. Med. Image Anal. 2019, 57, 237–248. [Google Scholar] [CrossRef]
  59. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep residual learning for image recognition. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. [Google Scholar]
  60. Alghamdi, H.S.; Tang, H.L.; Waheeb, S.A.; Peto, T. Automatic optic disc abnormality detection in fundus images: A deep learning approach. In Proceedings of the Ophthalmic Medical Image Analysis MICCAI 2016 Workshop (MICCAI 2016), Athens, Greece, 17–21 August 2016; pp. 10–17. [Google Scholar]
  61. Chen, X.; Xu, Y.; Wong, D.W.K.; Wong, T.Y.; Liu, J. Glaucoma detection based on deep convolutional neural network. In Proceedings of the 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Milan, Italy, 25–29 August 2015; pp. 715–718. [Google Scholar]
  62. Yang, X.; Lin, Y.; Wang, Z.; Li, X.; Cheng, K.T. Bi-Modality Medical Image Synthesis Using Semi-Supervised Sequential Generative Adversarial Networks. IEEE J. Biomed. Health Inform. 2020, 24, 855–865. [Google Scholar] [CrossRef] [PubMed]
  63. Moseley, M.; Donnan, G. Multimodality imaging: Introduction. Stroke 2004, 35, 2632–2634. [Google Scholar] [CrossRef] [Green Version]
  64. Costa, P.; Galdran, A.; Meyer, M.I.; Niemeijer, M.; Abràmoff, M.; Mendonça, A.M.; Campilho, A. End-to-end adversarial retinal image synthesis. IEEE Trans. Med. Imaging 2017, 37, 781–791. [Google Scholar] [CrossRef] [PubMed]
  65. Liu, M.-Y.; Tuzel, O. Coupled generative adversarial networks. arXiv 2016, arXiv:1606.07536. [Google Scholar]
  66. Zhu, J.-Y.; Park, T.; Isola, P.; Efros, A.A. Unpaired image-to-image translation using cycle-consistent adversarial networks. In Proceedings of the IEEE International Conference on Computer Vision, Venice, Italy, 22–29 October 2017; pp. 2223–2232. [Google Scholar]
  67. Isola, P.; Zhu, J.-Y.; Zhou, T.; Efros, A.A. Image-to-image translation with conditional adversarial networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 1125–1134. [Google Scholar]
  68. Wang, X.; Peng, Y.; Lu, L.; Lu, Z.; Bagheri, M.; Summers, R.M. Chestx-ray8: Hospital-scale chest X-ray database and benchmarks on weakly-supervised classification and localization of common thorax diseases. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 2097–2106. [Google Scholar]
  69. Yao, L.; Prosky, J.; Poblenz, E.; Covington, B.; Lyman, K. Weakly supervised medical diagnosis and localization from multiple resolutions. arXiv 2018, arXiv:1803.07703. [Google Scholar]
  70. Gu, L.; Zhang, X.; You, S.; Zhao, S.; Liu, Z.; Harada, T. Semi-supervised learning in medical images through graph-embedded random forest. Front. Neuroinformatics 2020, 14, 49. [Google Scholar] [CrossRef] [PubMed]
  71. Liu, X.; Song, M.; Tao, D.; Liu, Z.; Zhang, L.; Chen, C.; Bu, J. Random Forest Construction with Robust Semisupervised Node Splitting. IEEE Trans. Image Process. 2015, 24, 471–483. [Google Scholar] [CrossRef] [PubMed]
  72. Bakht, A.B.; Javed, S.; AlMarzouqi, H.; Khandoker, A.; Werghi, N. Colorectal Cancer Tissue Classification Using Semi-Supervised Hypergraph Convolutional Network. In Proceedings of the 2021 IEEE 18th International Symposium on Biomedical Imaging (ISBI), Nice, France, 13–16 April 2021; pp. 1306–1309. [Google Scholar]
  73. Francesco, P.; Enrico, M.; Ficarra, E. Colorectal Cancer Classification using Deep Convolutional Networks. An Experimental Study. In Proceedings of the 11th International Joint Conference on Biomedical Engineering Systems and Technologies (BIOIMAGING 2018), Madeira, Portugal, 19–21 January 2018; pp. 58–66. [Google Scholar]
  74. Liu, Q.; Yang, H.; Dou, Q.; Heng, P.-A. Federated Semi-Supervised Medical Image Classification via Inter-Client Relation Matching; Springer: Cham, Switzerland, 2021; pp. 325–335. [Google Scholar]
  75. McMahan, B.; Moore, E.; Ramage, D.; Hampson, S.; Arcas, B.A.y. Communication-efficient learning of deep networks from decentralized data. In Proceedings of the Artificial Intelligence and Statistics, Fort Lauderdale, FL, USA, 20–22 April 2017; pp. 1273–1282. [Google Scholar]
  76. Su, L.; Liu, Y.; Wang, M.; Li, A. Semi-HIC: A novel semi-supervised deep learning method for histopathological image classification. Comput. Biol. Med. 2021, 137, 104788. [Google Scholar] [CrossRef]
  77. Su, H.; Shi, X.; Cai, J.; Yang, L. Local and global consistency regularized mean teacher for semi-supervised nuclei classification. In Proceedings of the International Conference on Medical Image Computing and Computer-Assisted Intervention, Shenzhen, China, 13–17 October 2019; pp. 559–567. [Google Scholar]
  78. Zhou, D.; Bousquet, O.; Lal, T.; Weston, J.; Schölkopf, B. Learning with local and global consistency. In Proceedings of the 16th International Conference on Neural Information Processing Systems, Whistler, BC, Canada, 9–11 December 2003. [Google Scholar]
  79. Bromley, J.; Guyon, I.; LeCun, Y.; Säckinger, E.; Shah, R. Signature verification using a “siamese” time delay neural network. In Proceedings of the 6th International Conference on Neural Information Processing Systems, Denver, CO, USA, 29 November–2 December 1993. [Google Scholar]
  80. Wang, X.; Chen, H.; Xiang, H.; Lin, H.; Lin, X.; Heng, P.-A. Deep virtual adversarial self-training with consistency regularization for semi-supervised medical image classification. Med. Image Anal. 2021, 70, 102010. [Google Scholar] [CrossRef]
  81. Mao, J.; Yin, X.; Zhang, G.; Chen, B.; Chang, Y.; Chen, W.; Yu, J.; Wang, Y. Pseudo-labeling generative adversarial networks for medical image classification. Comput. Biol. Med. 2022, 147, 105729. [Google Scholar] [CrossRef]
  82. Oord, A.v.d.; Li, Y.; Vinyals, O. Representation learning with contrastive predictive coding. arXiv 2018, arXiv:1807.03748. [Google Scholar]
  83. Guo, L.; Wang, C.; Zhang, D.; Xu, K.; Huang, Z.; Luo, L.; Peng, Y. Semi-supervised medical image classification based on CamMix. In Proceedings of the 2021 International Joint Conference on Neural Networks (IJCNN), Shenzhen, China, 18–22 July 2021; pp. 1–7. [Google Scholar]
  84. Yun, S.; Han, D.; Oh, S.J.; Chun, S.; Choe, J.; Yoo, Y. Cutmix: Regularization strategy to train strong classifiers with localizable features. In Proceedings of the IEEE/CVF International Conference on Computer Vision, Seoul, Republic of Korea, 27–28 October 2019; pp. 6023–6032. [Google Scholar]
  85. Harris, E.; Marcu, A.; Painter, M.; Niranjan, M.; Prügel-Bennett, A.; Hare, J. Fmix: Enhancing mixed sample data augmentation. arXiv 2020, arXiv:2002.12047. [Google Scholar]
Figure 1. Percentage of studies read from different sources.
Figure 1. Percentage of studies read from different sources.
Applsci 12 12094 g001
Figure 2. Diagram of the types of SSL approaches and the datasets used to test each method.
Figure 2. Diagram of the types of SSL approaches and the datasets used to test each method.
Applsci 12 12094 g002
Figure 3. The performance of the methods on all datasets. The best result for each dataset is bolded.
Figure 3. The performance of the methods on all datasets. The best result for each dataset is bolded.
Applsci 12 12094 g003
Table 1. Comparison of consistency-based methods.
Table 1. Comparison of consistency-based methods.
MethodTypeDataset UsedBest Results
SRC-MT [19]Consistency-basedChest X-ray14 and ISIC2018X-ray14: AUC 79.23
ISIC2018: AUC 93.58, F1 60.68
NoT [52]Consistency-basedChest X-ray14, RSNA Brain CT, Knee MRNetX-ray14: AUC 79.49
RSNA: AUC 93.31
MRNet: AUC 91.49
S2MTS2 [55]Consistency-basedChest X-ray14, ISIC2018, CheXpertX-ray14: AUC 82.51
ISIC2018: AUC 94.71, F1 62.67
CheXpert: AUC 71.58
Table 2. Comparison of adversarial methods.
Table 2. Comparison of adversarial methods.
MethodTypeDataset UsedBest Results
Uncertainty-guided VAT with BNN [56]Adversarial (VAT-based)Chest X-ray14, ISIC2018 Skin Lesion Analysis, Hip X-rayX-ray14: 80.69 AUC
ISIC2018: AUC 96.04, F1 69.67
Hip X-ray: AUC 92.43, F1 72.39
SSAC [58]Adversarial (GAN-based)LIDC-IDRI and Tianchi Lung Nodule dataset combinedAUC 92.82, Accuracy 90.58
Bi-Modality Image Synthesis [62]Adversarial (GAN-based)Locally collected dataset and PROSTATEx combinedAccuracy 93
Table 3. Comparison of graph-based methods.
Table 3. Comparison of graph-based methods.
MethodTypeDataset UsedBest Results
GraphXNET [34]Graph-basedChest X-ray14AUC 0.78
Graph-Embedded RF [70]Graph-basedDRIVE 2D retinal vessel and Big Neuron 3D neuronalDRIVE: accuracy 79.42
Big Neuron: accuracy 74.16
CRC with HGNN [72]Graph-basedCRC WSI datasetF1 0.94
Table 4. Comparison of other methods.
Table 4. Comparison of other methods.
MethodTypeDataset UsedBest Results
ACPL [37]Pseudo-labelingChest X-ray14 and ISIC2018 Skin Lesion AnalysisChest X-ray: AUC 81.77
ISIC2018: AUC 94.36, F1 62.23
FL with IRM [74]Federated learningRSNA Brain CT and ISIC2018 Skin Lesion AnalysisRSNA: AUC 87.56, F1 59.86
ISIC2018: AUC 92.46, F1 55.81
Semi-HIC [76]Learning by associationBioimaging2015 and BACHBioimaging2015: error rate 25.23
BACH: error rate 15.46
Table 5. Comparison of hybrid methods.
Table 5. Comparison of hybrid methods.
MethodTypeDataset UsedBest Results
Local and Global Consistency Regularized [77]Consistency-based + Graph-basedMoNuSeg and Ki-67MoNuSeg: F1 76.89
Ki-67: F1 79.79
Wang et al. [80]Consistency-based + Adversarial (VAT)Breast Ultrasound (BRUS) and OCTBRUS: F1 0.88
OCT: Macro F1 0.91
PLGAN [81]Pseudo-labeling + Adversarial (GANs)OCT, COVID-19, Brain Tumor MRI, Chest X-rayAccuracies:
OCT: 97.10
COVID-19: 89.30
MRI: 86.80
X-ray: 97.50
SSL with CamMix [83]Consistency-based + CamMix (MSDA)ILD dataset, ISIC2018ILD: AUC 95.34
ISIC2018: AUC 94.04
Table 6. Comparison of methods tested on the Chest X-ray14 and ISIC2018 datasets, at 20% labeled data for all models.
Table 6. Comparison of methods tested on the Chest X-ray14 and ISIC2018 datasets, at 20% labeled data for all models.
ModelTypeChest X-ray14ISIC2018
AUCAUCF1
SRC-MT [19]Consistency-based79.2393.5860.68
NoT [52]Consistency-based79.49--
S2MTS2 [55]Consistency-based82.5194.7162.67
Uncertainty-guided VAT with BNN [56]Adversarial (VAT-based)80.6996.0469.67
GraphXNET [34]Graph-based78.00--
ACPL [37]Pseudo-labeling81.7794.3662.23
FL with IRM [74]Federated learning-92.4655.81
SSL with CamMix [83]Consistency-based + CamMix (MSDA)-94.04-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Solatidehkordi, Z.; Zualkernan, I. Survey on Recent Trends in Medical Image Classification Using Semi-Supervised Learning. Appl. Sci. 2022, 12, 12094. https://doi.org/10.3390/app122312094

AMA Style

Solatidehkordi Z, Zualkernan I. Survey on Recent Trends in Medical Image Classification Using Semi-Supervised Learning. Applied Sciences. 2022; 12(23):12094. https://doi.org/10.3390/app122312094

Chicago/Turabian Style

Solatidehkordi, Zahra, and Imran Zualkernan. 2022. "Survey on Recent Trends in Medical Image Classification Using Semi-Supervised Learning" Applied Sciences 12, no. 23: 12094. https://doi.org/10.3390/app122312094

APA Style

Solatidehkordi, Z., & Zualkernan, I. (2022). Survey on Recent Trends in Medical Image Classification Using Semi-Supervised Learning. Applied Sciences, 12(23), 12094. https://doi.org/10.3390/app122312094

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