Background This study evaluates the performance of MedGemma-4B, a specialized medical vision-language model (VLM), on six fundamental radiological image classification tasks. We hypothesize that, despite their diagnostic capabilities, such models may underperform on clinically essential perceptual tasks that may be underrepresented in training data. Methods MedGemma-4B was assessed using 600 multicenter radiological images, divided equally into six classification tasks: modality, body part, orientation, contrast, image mode, and organ. Using standardized visual question answering (VQA) prompts, model outputs were categorized by two expert radiologists as correct, incorrect, or imprecise. Performance was evaluated using accuracy, sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and Matthews correlation coefficient (MCC), with 95% confidence intervals and χ2 tests. Results Performance varied substantially across tasks. The model demonstrated robust accuracy in modality (97%) and body part identification (80%), but critically underperformed in orientation (41%), contrast phase recognition (32%), and image mode classification (48%). Organ identification achieved moderate accuracy (74%). MCC values showed strong correlation for modality (96.2%) and body part (78.7%), but poor correlation for orientation (33.1%). Statistical analysis confirmed significant performance disparities (p
Artificial intelligence (AI) has rapidly transitioned from a futuristic concept to an integral part of everyday life, with AI-driven algorithms now embedded in a broad range of technologies.1 Modern AI systems rely on machine learning trained on vast datasets, enabling specialized tasks with high efficiency.2 Notable examples include AlphaFold for protein prediction and OpenAI’s GPT models for natural language processing.3–5 In medicine, AI shows promise in enhancing diagnostics and workflow automation, with models like BioBERT and LYNA demonstrating expert-level performance in specific tasks.6–8
Subsequently, numerous models have been developed with medical imaging as their central focus. For instance, CheXNet was designed to detect pneumonia, and BioViL took it a step further by allowing multiple diagnoses on chest X-rays.9,10
Among these, vision language models (VLMs) show considerable promise in analyzing biomedical images.11 They are AI systems designed to process both images and text together. They combine a visual encoder (like a CNN or Vision Transformer) that processes images with a language model (such as BERT or GPT) that interprets text. Trained on large datasets of image-text pairs, VLMs can perform tasks such as image captioning, visual question answering, and multimodal search. Despite their promise, they can struggle with tasks outside their training data and require careful validation before clinical use.12
One of these recent innovations is Google’s Med-PaLM family of models, which includes MedGemma-4B, adapted for use in medical contexts.13–15 MedGemma builds upon the foundation of general-purpose AI models, refining them through fine-tuning on medical imaging datasets. Recent studies have already shown that, fine-tuned for chest X-ray analysis, it ranked above all its competitors in terms of diagnostic accuracy.16
While such models are marketed as being able to solve complex diagnostic tasks, performance variability may occur in tasks underrepresented in training data, regardless of perceived task simplicity, such as classification or multimodal comprehension.12,17 Therefore, we hypothesized that even medically specialized AI models, such as MedGemma-4B, may exhibit significant limitations when applied to tasks that, despite being clinically routine, may be underrepresented in their training data, and therefore cannot be assumed to be simple for the model. To explore this hypothesis, we submitted a series of publicly accessible radiology images to the MedGemma-4B model. Our goal was not to evaluate its diagnostic capability, but rather to assess its performance on image interpretation tasks that clinicians consider fundamental, but whose simplicity for a VLM cannot be assumed, as model performance may be strongly influenced by training data coverage rather than clinical task complexity.18,19 This study aims to provide insights into the generalization capabilities and limitations of emerging medical AI tools. It raises important questions about their future readiness for integration into everyday clinical practice.
The evaluated AI model, MedGemma-4B-it, is an advanced architecture combining a vision encoder with a large language model (LLM) decoder, featuring 4 billion parameters. The model (MedGemma-4B-it version 1.0.0, release date 20/05/2025) was accessed and utilized via the Hugging Face Transformers library (https://huggingface.co/google/medgemma-4b-it , accessed on 16/06/2025), enabling seamless integration with the underlying architecture.20 The model was deployed locally via an integrated development environment (IDE) setup using Visual Studio Code (VS Code), v1.101.0 (Universal), coupled with Python, v3.10.1. This local deployment ensured full control over data privacy and model execution. No additional fine-tuning was performed on the model; it was employed in its original pre-trained state without further adaptation. The model inference and experimentation were conducted on a personal, non-commercial workstation, specifically a 2021 MacBook Pro equipped with an M1 Pro chip and 32 GB of RAM, as processing time optimization was not an issue in this academic setting.
The dataset comprised a total of 600 medical images, systematically divided into six distinct classification tasks, each containing 100 images ( Figure 1). All images were sourced from the open-access platform Radiopaedia.org, which is not part of MedGemma’s disclosed training datasets, made available under the Creative Commons Non-Commercial Attribution CC-NC-BY-SA 3.0 license.21 Proper citation of each case was ensured in accordance with the licensing requirements. Detailed information on the images used, including the DOI of each source case and the corresponding image URL, is provided in Extended data 1. All images were fed into the model in their native resolution and content, without resizing, cropping, or other pre-processing modifications. Images were carefully selected to exclude any identifying information, such as embedded patient data, facial features, or institutional labels, in accordance with ethical and privacy standards.
A total of 600 medical images were selected and divided into six distinct groups of 100 images each, corresponding to specific classification tasks. Each group was independently submitted to the vision-language model (MedGemma-4B-it) for analysis and performance assessment.
The six datasets were designed to evaluate the VLM’s capabilities across fundamental radiological competencies: (1) Modality Dataset: testing recognition of basic radiological imaging modalities techniques (e.g., X-ray, CT, MRI, etc., 5 subgroups in total); (2) Body part Dataset: testing anatomical knowledge of the broader body part imaged in the study (e.g., brain, thorax, abdomen, etc., 5 subgroups in total, further divided into 16 sub-classes such as hand, wrist, elbow etc.); (3) Orientation Dataset: testing specific spatial orientation of the image (e.g., axial, coronal, sagittal, etc. 8 subgroups in total); (4) Contrast Dataset: testing the absence or presence of contrast agents and specific contrast phases (e.g., non-contrast, arterial, venous, etc., 8 subgroups in total); (5) Image mode Dataset: testing specific image setting such as CT window (e.g., mediastinal window, lung window, etc.), MRI sequence (e.g., T1, T2, etc.) or US modes (e.g. B-mode, Doppler) (10 subgroups in total); (6) Organ Dataset: testing the recognition of the main/central organ in the image (e.g. heart, thyroid, etc., 7 subgroups in total). The datasets were conceived to cover a broad spectrum of imaging scenarios. For each dataset, 5 to 10 subgroups of 10 to 20 images were created to ensure diversity and facilitate the identification of areas of high or low performance on specific subtasks. Dataset groups and subgroups are summarized in Table 1.
The body part group is further divided into subclasses for the brain/neck, upper limbs and lower limbs subgroups. Complete subgroup definitions and individual images are provided in Extended data 1.
For each dataset, a maximum of three images per patient case were selected, avoiding similar images, focusing on normal or typical pathological cases while avoiding extreme cases (e.g., unilateral lung, distortions, blurriness, etc.). Paediatric cases were excluded from the datasets. Within each subcategory, images were varied according to orientation, contrast, and modality to ensure a robust evaluation framework.
The tasks assigned to the VLM were carefully structured using a Visual Question Answering (VQA) framework. For every medical image input presented to the model, a standardized series of questions was formulated to assess the model’s capability to identify the key features related to the specific dataset task, in accordance with recent best practices in prompt engineering.22,23 This involved explicitly defining the role that the model was expected to take on, namely, that of an imaging specialist, providing context for required expertise. The prompts clearly outlined the nature of the task, needing precise identification of one specific aspect per question. The expected response format was constrained to a single word to maintain uniformity in the model’s answers. Examples of acceptable answers were incorporated within the prompt design, guiding the model towards the desired output. The complete set of prompts, customized for each dataset used in this study, is available in Extended data 2. For all tasks, the model’s top textual output was directly taken as the prediction, without applying any probabilistic thresholding.
For each question, model outputs were categorized as correct, incorrect, or imprecise, by two expert radiologists (4 and 8 years of experience). A response was considered correct if it was fully accurate or provided either a synonym or a more specific yet valid term (e.g., labelling a chest X-ray as ‘antero-posterior’ instead of ‘frontal’). An output was classified as incorrect if it used a categorically wrong term (e.g., referring to an MRI scan as CT). A response was deemed imprecise if it was ambiguous or partially accurate (e.g., ‘limb’ instead of ‘lower limb’; or labelling a ‘brain’ image as ‘neck’ due to marginal cervical structures). Cases where the classification of an answer was ambiguous were determined through agreement.
Statistical analyses were performed using R (v4.3.2) through RStudio (v2024.09.0 + 375). Accuracy was calculated as the proportion of correct responses, with 95% confidence intervals computed for each of the six datasets. Chance-levels (i.e., expected accuracy under random guessing) were reported, based on the number of possible categories. Other performance metrics included sensitivity (SE), specificity (SP), positive predictive value (PPV), negative predictive value (NPV), and Matthews correlation coefficient (MCC). 95% confidence intervals were computed for SE, SP, PPV, and NPV using the Wilson method, and for MCC using bootstrap resampling (10,000 iterations, percentile method). These metrics were selected as they are standard and complementary in AI diagnostic studies: accuracy offers an intuitive overall performance measure, SE reflects the model’s ability to correctly identify each class, SP its ability to reject non-members, PPV and NPV the reliability of positive and negative predictions respectively, and MCC provides a single balanced summary that is robust to class imbalance.
For each multi-class task (“modality”, “body part”, “orientation”, etc.), each class (“axial”, “sagittal”, etc.) was treated in a one-vs-rest framework. Predicted labels matching the ground-truth class in positive instances were counted as true positives (TP). False negatives (FN) corresponded to cases where the ground-truth class was not correctly identified, including instances assigned to another class. For each class, all samples not belonging to the class of interest formed the negative set: predictions incorrectly assigned to the class were counted as false positives (FP), while correctly rejected negatives were counted as true negatives (TN). Regarding “imprecise” responses, these were systematically treated as incorrect predictions for the purpose of metric computation, reflecting a conservative evaluation strategy. To obtain task-level metrics, TP, FN, FP, and TN counts were summed across all classes in the task to create an aggregated confusion matrix for the entire task. A one-vs-rest decomposition was used to derive class-specific performance metrics in the multi-class setting. This approach is standard for multi-class problems without a natural binary structure, and provides a consistent and interpretable framework for computing sensitivity, specificity, predictive values, and MCC across tasks with varying numbers of categories. Although minor class imbalance was present (maximum class ratio 2:1), it was limited by design, and the use of MCC reduces the impact of imbalance on performance estimation.
Comparative analyses were performed both between tasks (inter-group) and within subgroups of each task (intra-group) using a Chi-squared (χ2) test, followed by post-hoc analysis with Holm correction to identify differences. A significance threshold of p < 0.05 was adopted. Reproducibility was assessed by resubmitting a subset of test cases (specifically, the 1st, 11th, 21st, etc., image from each dataset) and compared with a McNemar test.
This study was conducted and reported in accordance with the CLAIM (Checklist for Artificial Intelligence in Medical Imaging) 2024 guidelines, available in Extended data 3.24
MedGemma’s classification performance was evaluated across six radiological image datasets, including modality, body part, orientation, contrast, image mode, and organ. The VLM model demonstrated an overall classification accuracy of 62% (95% CI: 58.12%–65.88%), with important performance variation across different diagnostic tasks. Overall model performances are summarized in Figure 2 and Tables 2 and 3.
Correct answers are highlighted in green, imprecise answers in yellow and incorrect answers in orange. Values are expressed as percentages of the total for each group.
For each category, the number of images correctly classified (Correct), misclassified (Incorrect), and imprecisely classified (Imprecise) are reported, along with the total number of images evaluated (Total, N = 100 per category). Accuracy is presented as a percentage, with the corresponding 95% confidence interval (CI95). Full results with subgroup performances of each dataset are available in Extended data 4.
For each category, sensitivity (SE), specificity (SP), positive predictive value (PPV), negative predictive value (NPV), and Matthews correlation coefficient (MCC) are reported. All values are expressed as percentages rounded to one decimal place. Metrics were computed using a one-vs-rest approach for each class; task-level values were derived from aggregated confusion matrices obtained by summing TP, FP, TN, and FN across classes. 95% CIs for SE, SP, PPV, and NPV were computed using the Wilson method; 95% CIs for MCC were derived by bootstrap resampling (10,000 iterations, percentile method).
The model was the most accurate in correctly identifying the image modality, achieving a correct classification rate of 97%, with only 3 errors and no imprecise outputs (95% CI: 90.99%–99.25%, chance level: 20%). Among the three errors, one involved misclassifying an axial abdominal T2-weighted MRI slice as a CT scan, while the remaining two errors resulted from mistaking angiography images for conventional X-rays. All remaining categories (X-ray, CT, and ultrasound) were classified with perfect accuracy (100%).
Performance was also high in body part identification, with an accuracy of 80% (95% CI: 70.93%–86.59%, chance level: 20%). While most predictions were correct, 11% of outputs were deemed imprecise, typically involving ambiguous or overlapping anatomical regions (e.g, “wrist” was predicted instead of “hand”). Additionally, 9% were incorrect, with the majority arising from musculoskeletal misclassifications (e.g., knee instead of hand).
By contrast, performance dropped in the remaining tasks. In identifying image orientation, the model achieved only 41% accuracy (95% CI: 32.06%–51.10%, chance level: 12.5%), with 59% of images misclassified and no imprecise cases. Similarly, contrast usage yielded an accuracy of just 32% (95% CI: 23.62%–41.84%, chance level: 12.5%). This category also had the highest rate of imprecise answers (17%). Classification by image mode yielded an accuracy of 48% (95% CI: 38.41%–57.75%, chance level: 10%). Nearly as many predictions were incorrect (42%) as correct, with 9% of outputs classified as imprecise.
Finally, in the organ classification task, MedGemma achieved a moderate accuracy of 74% (95% CI: 64.81%–81.75%, chance level: 14.3%). The model demonstrated proficiency in distinguishing among common organs but occasionally struggled with images that featured multiple structures or had subtle anatomical boundaries. Detailed subgroup performance metrics for each dataset are provided in Figure 3 and Extended data 4. Furthermore, examples of typical misclassifications across the different datasets are illustrated in Figure 4.
Stacked bar charts illustrate the distribution of prediction outcomes for each subgroup within the six datasets: (A) modality classification (e.g., X-ray, CT, MRI), (B) body part classification (e.g., brain/neck, chest, abdomen), (C) orientation classification (e.g., axial, coronal, sagittal views), (D) contrast classification (e.g., arterial CT, venous MRI), (E) image mode classification (e.g., T1/T2 MRI, B-mode ultrasound), and (F) organ classification (e.g., heart, liver, kidney). Each horizontal bar represents the proportion of correct (lightest shade) and imprecise (darker shade) predictions; the unfilled portion of each bar corresponds to incorrect classifications (e.g., a fully empty bar indicates 100% incorrect predictions). Values are expressed as percentages of the total for each subgroup.
The errors include: modality (MRI predicted as CT, A); anatomical region (Neck as Brain, B); imaging plane (Sagittal as Axial, C); contrast (Arterial as Non-contrast, D); image mode (FLAIR as T1, E) and organ (Spleen as Liver, F).
Sensitivity and specificity varied markedly across tasks. Modality and body part classification achieved near-perfect sensitivity (97.0% and 80.0%) and specificity (99.2% and 98.7%), whereas orientation and contrast tasks showed low sensitivity (41.0% and 32.0%) despite high specificity (91.9% and 91.0%). MCC values confirmed these disparities, with strong correlation for modality (96.2%) and body part (78.7%), but moderate-to-poor correlation for orientation (33.1%) and contrast (23.5%).
A chi-squared test revealed a statistically significant difference in classification accuracy across the six tasks (p < 0.001). Pairwise post hoc comparisons with Holm-adjusted p-values confirmed that the modality group exhibited significantly higher accuracy than all other groups (p < 0.001 for all comparisons). The body part group also outperformed orientation, contrast and image-mode (p < 0.001), but did not differ significantly from the organ group (p = 0.79). The organ group, in turn, showed significantly higher accuracy than orientation, contrast and image-mode groups (p < 0.01), and was statistically comparable to the body part group. No significant differences were observed between orientation, contrast and image-mode groups (all p > 0.12).
The reproducibility of MedGemma’s predictions across repeated runs was high, with a mean reproducibility of 91.6% across all tasks. Perfect consistency (100%) was observed for both the modality and organ classifications. Body part, orientation, and contrast tasks each showed 90% reproducibility, reflecting minor variability. The lowest reproducibility was noted in the image mode task (80%). All accuracy results reported in this study reflect the model’s first-run predictions. McNemar’s test showed no statistically significant difference between initial and repeated predictions (χ2 = 0.80, p = 0.37), indicating that observed discrepancies likely occurred by chance and do not reflect systematic variability.
Our evaluation of MedGemma-4B revealed significant variability in its performance across fundamental radiological tasks. The model demonstrated near-perfect accuracy in identifying imaging modalities (97%) and strong performance in anatomical localization (80% for body parts), suggesting robust capabilities in high-level broad categorizations. However, its proficiency dropped markedly for tasks requiring finer-grained distinctions: orientation (41%), contrast (32%), and image modes (48%) tasks were notably weaker, although they remained above chance. Organ identification achieved moderate accuracy (74%) but is more debatable for the definition of ground truth compared to the other tasks, since multiple other organs can sometimes be partially visible in the images. This disparity in performances underscores a critical gap in the model’s generalization abilities, excelling in broad categorizations possibly better represented in its training data, but faltering in tasks whose fine-grained labels may be less consistently or comprehensively captured in large-scale medical datasets.
These findings may appear counterintuitive to clinicians with limited familiarity with AI development, who may intuitively assume that tasks routinely mastered by radiology trainees should be straightforwardly solved by an advanced VLM. However, this assumption conflates human-defined task simplicity with model performance, overlooking a critical factor: training data dependency. The model’s performance is likely influenced by the scope and quality of its training dataset. While this study focuses on MedGemma-4B, similar evaluations on other VLMs are needed to confirm generalizability, as performance strongly depends on training data coverage and model design.
Thus, the observed patterns raise important questions about the origins of these limitations. For example, the model’s high accuracy for modalities such as X-ray (a frequently encountered imaging type in training datasets) contrasts with its poor performance in underrepresented tasks (e.g., contrast phase identification, or image mode, which may lack consistent annotation in large-scale, curated radiology datasets). This suggests that gaps in training data diversity may partially explain the results; however, inherent challenges in visual ambiguity (e.g., distinguishing between arterial and venous contrast phases on static images) likely compound the issue.
As a newly released model in the AI ecosystem, MedGemma has yet to be extensively evaluated in independent studies.25 Prior to MedGemma, several medical VLMs were developed, but none combined its three key advantages: open accessibility, training dataset scale/diversity, and multimodal versatility. For example, Med3DVLM tackles 3D medical vision-language alignment, achieving 61% R@1 retrieval on volumetric data, but was trained on a single M3D dataset, lacking the multi-source, multi-domain training foundation of models like Med-Gemma.26 Meanwhile, MedVLM-R1 showed efficient reasoning capabilities for radiology QA through reinforcement learning (78.22% accuracy with just 2B parameters), but its training data remains partially restricted, and its application is limited to closed-set questions.27 While many of these models report strong performance across various standardized, closed-benchmark systems, only a limited number have undergone clinical validation to date.28
Multiple benchmarks have been proposed for medical VLM in the literature.12,29–32 However, while such evaluations are essential for standardization and reproducibility, they also present notable limitations. Models tend to over-optimize for benchmark tasks, potentially at the expense of real-world clinical applicability. Smaller-scale, clinically grounded evaluations, although more limited, are necessary to bridge this gap; yet, they remain markedly underrepresented in the current academic landscape.33
The current work focused on fundamental radiological classification tasks, areas that have remained relatively underexplored in the existing literature. For instance, a 2019 study developed a classifier achieving 94.8% accuracy in identifying imaging modality and anatomical region.34 Similarly, a 2020 study developed a modality classification system with ResNet-50, achieving 87.9% accuracy across 31 imaging classes, demonstrating the potential of such systems to support clinical case retrieval from large-scale archives.35 Furthermore, a 2021 study developed a model capable of classifying both imaging modality and anatomical location with over 99.5% accuracy on converted JPEG images, demonstrating classification feasibility even when DICOM metadata is unavailable.36 While these results are promising, they remain insufficient given the vast number of models yet to be evaluated and the immense diversity of potential clinical tasks in radiology.
The findings of this study highlight both the potential of MedGemma-4B and other possible future medical VLMs in real-world radiological workflows. The model’s high accuracy in modality and body part classification suggests its utility for automated pre-sorting of medical images, streamlining the organization of large-scale PACS (Picture Archiving and Communication Systems) databases. Future iterations or other models could achieve robust performance across finer-grained tasks. For instance, reliable contrast phase recognition (e.g., arterial vs. venous) could automate protocol verification, reducing human error in dynamic imaging studies. Similarly, precise MRI sequence identification might enable real-time quality control, flagging mislabeled sequences or suboptimal acquisitions before interpretation. Beyond classification, advanced models could standardize ambiguous annotations (e.g., “oblique” projections or “portal phase” timing) by learning consensus patterns from multi-institutional datasets, bridging variability in radiological lexicons. Such capabilities would unlock applications like intelligent PACS structuring, where AI pre-annotates studies with metadata (e.g., “contrast-enhanced CT, arterial phase, axial liver”), accelerating retrieval and triage. Moreover, the seamless integration of multimodal reasoning, which combines image features with clinical context from reports, could support the creation of automated preliminary summaries, thereby streamlining radiologist workflows. While current models remain imperfect, their trajectory points toward a future where AI handles routine perceptual tasks, allowing clinicians to focus on complex diagnostics and patient care.
This study has several limitations that should be considered when interpreting the results. First, the dataset size for specific subcategories was relatively small, which may limit the statistical power to detect subtle performance differences. While our image source provided diverse, multicentric data, the images may not fully represent the variability seen in single-institution datasets, where acquisition protocols and artifacts can differ. Additionally, model responses were manually categorized (correct, incorrect, or imprecise) by expert radiologists, introducing potential subjectivity in borderline cases. While the introduction of the imprecise category was intended to reduce forced binary judgments, it also creates an additional boundary where subjectivity may arise, and future studies could benefit from formal inter-rater reliability assessment. Furthermore, the one-vs-rest aggregation strategy used to compute task-level metrics inherently expands the true negative pool, which can inflate specificity and NPV. The MCC was included to counterbalance this, as it remains robust to such class imbalance. Furthermore, all evaluations were conducted on single static images, consistent with MedGemma’s design for image-level interpretation. However, modalities such as ultrasound or fluoroscopy are typically assessed dynamically, so this may not fully represent real-world conditions. Finally, the evaluation employed a fixed prompt structure without exploring variations in phrasing or context, which may have influenced MedGemma’s performance. Future studies could benefit from larger, prospectively collected datasets and systematic prompt optimization to better assess real-world applicability.
MedGemma-4B demonstrated strong performance on certain foundational radiological tasks, such as modality and anatomical region identification, while exhibiting variable or lower performance on other tasks considered clinically straightforward, such as image contrast or orientation classification. This pattern indicates that tasks perceived as simple by clinicians are not necessarily easy for AI models, whose performance may be strongly influenced by the coverage and quality of their training data rather than human-perceived task simplicity.
This study utilised publicly available anonymised medical images from Radiopaedia.org, which are distributed under the Creative Commons Non-Commercial Attribution CC-NC-BY-SA 3.0 license. No patient identifying information was included. As no human participants were directly involved and no identifiable data were used, formal ethical approval was not required.
The data and code supporting the findings of this study are openly available in the Zenodo repository: https://doi.org/10.5281/zenodo.21456064.37
The Zenodo repository contains the following files supporting this study:
• Extended data 1: input dataset containing image URLs, source DOIs, and case information.
• Extended data 2: complete prompt designs for all six evaluation tasks.
• Extended data 3: CLAIM 2024 checklist.
• Extended data 4: detailed subgroup performance metrics, including the values underlying the reported analyses and figures.
• medgemma-online.py: body part classification script.
• medgemma-reproducibility.py: reproducibility script.
The deposited files are available under the Creative Commons Attribution 4.0 International license (CC-BY 4.0).
AI-based language tools (ChatGPT by OpenAI) were used to improve the clarity and style of the manuscript.