$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Diese Studie wurde vom Institutional Review Board des Shenzhen University General Hospital genehmigt (IRB-Zulassungsnr. KYLL-2026-077-1). Die Anforderung einer schriftlichen informierten Einwilligung wurde aufgrund des retrospektiven Designs der Studie aufgehoben.
1. Studienkohorte
Patienten mit akutem ischämischen Schlaganfall (AIS), die sich einer multiparametrischen MRT unterzogen hatten, wurden rückwirkend aus der institutionellen klinischen Datenbank identifiziert. Der gesamte Studienablauf, einschließlich Datensatzaufteilung, Entwicklung von Bildmodellen, klinischer Modellkonstruktion und multimodaler Fusion, ist in Abbildung 1 dargestellt.

Abbildung 1: Arbeitsablauf des multimodalen Ergebnisvorhersage-Frameworks. (A) Datenerfassung. Patienten mit akutem ischämischen Schlaganfall (AIS), die die vordefinierten Einschlusskriterien erfüllten, wurden rückwirkend identifiziert. Der Datensatz wurde in eine Trainings-Validierungskohorte (n = 250), eine interne, unabhängige Testkohorte (n = 50) und eine externe Testkohorte (n = 37) unterteilt. Innerhalb der Trainings-Validierungskohorte wurde eine geschichtete fünffache Kreuzvalidierung durchgeführt. (B) Entwicklung von Bildgebungsmodellen. Multiparametrische MRT-Sequenzen, einschließlich diffusionsgewichteter Bildgebung (DWI), scheinbarer Diffusionskoeffizienten (ADC) und T2-fluid-attenuated inversions-Recovery (T2-FLAIR), wurden mit einer hybriden 3-dimensionalen konvolutionellen neuronalen Netzwerk-Vision Transformer (CNN-ViT)-Architektur zur Ergebnisvorhersage verarbeitet. (C) Entwicklung klinischer Modelle. Strukturierte klinische Variablen wurden verwendet, um maschinelle Lernmodelle zur Vorhersage des 90-Tage-Funktionsergebnisses zu trainieren. (D) Multimodale Fusionsstrategie. Die von den bildgebenden und klinischen Modellen generierten Vorhersagen wurden mithilfe eines Stacked Logistic Regression Meta-Learners integriert, um Endergebnisvorhersagen zu erstellen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
MRT-Untersuchungen wurden mit klinischen MRT-Systemen durchgeführt. Für die Bildaufnahme wurde eine dedizierte Phased-Array-Kopfspule verwendet. Alle Sequenzen wurden in der Axialebene mit konsistenter Schnittpositionierung über verschiedene Modalitäten hinweg aufgenommen. Die diffusionsgewichtete Bildgebung (DWI) wurde mittels einer Einzelaufnahme mit Spin-Echo-Echo-Echo-Bildabbildung mit Wiederholungszeit (TR) von 3.000–5.000 ms und Echozeit (TE) von 80–90 ms aufgenommen. Die Diffusionssensibilisierung wurde mit b-Werten von 0 und 800–1.000 s/mm2 in mindestens drei orthogonalen Richtungen angewendet. Das Gesichtsfeld reichte von 220–240 mm mit einer Matrixgröße von 128 × 128. Die Scheibendicke betrug 5–6 mm mit einem Zwischenspalt von 1–1,5 mm. Es wurden zwei bis vier Signaldurchschnitte ermittelt.
Apparent Diffusion Coefficient (ADC)-Karten wurden automatisch aus DWI-Daten auf der Scanner-Arbeitsstation unter Verwendung von monoexponentieller Anpassung basierend auf den erfassten b-Werten erzeugt. Die ADC-Werte wurden voxelweise berechnet und zur quantitativen Analyse exportiert. T2-Aufnahmen der fluid-attenuated inversion recovery (T2-FLAIR) wurden mit einer Inversionswiederherstellungssequenz mit TR von 8.000–10.000 ms, TE von 80–140 ms und einer Inversionszeit von 2.200–2.600 ms aufgenommen. Das Sichtfeld reichte von 220–240 mm mit einer Matrixgröße von 192 × 192 bis 256 × 256. Die Scheibendicke betrug 4–5 mm mit einem Zwischenspalt von 1–1,5 mm.
2. Datenvorverarbeitung
Klinische Daten wurden aus strukturierten Tabellenkalkulationen importiert, die Fallkennungen, mRS-Outcome-Scores, Datensatz-Split-Labels sowie demografische und klinische Variablen enthielten. Fallkennungen wurden standardisiert, um die Konsistenz mit den Dateinamen der Bildgebung zu gewährleisten. Klinische Variablen umfassten demografische Merkmale, die Schwere des Schlaganfalls mit der National Institutes of Health Stroke Scale (NIHSS), vaskuläre Risikofaktoren und bei der Aufnahme erfasste Begleiterkrankungen.
Das primäre Ergebnis war der funktionelle Status 90 Tage nach Schlaganfallbeginn, gemessen mit der modifizierten Rankin-Skala (mRS). Ein günstiges Ergebnis wurde als mRS ≤ 2 definiert und ein ungünstiges als mRS > 2. Der Datensatz wurde zufällig in Trainings-Validierungs- und interne unabhängige Testkohorten unterteilt, wobei stratifizierte Stichproben auf Basis der mRS-Verteilung verwendet wurden, um das Ergebnisgleichgewicht zu erhalten. Eine externe Testkohorte wurde separat verarbeitet und während der Modellentwicklung nicht verwendet. Innerhalb der Trainings-Validierungskohorte wurde geschichtete k-fache Kreuzvalidierung angewendet, um konsistente Ergebnisverteilungen über Falten hinweg aufrechtzuerhalten.
Alle DWI-, ADC- und T2-FLAIR-Bände wurden vor dem Modelltraining vorbearbeitet, um räumliche und numerische Konsistenz zwischen den Modalitäten sicherzustellen. Die gleiche Preprocessing-Pipeline wurde auf die Trainings-Validierung, den internen Test und die externen Testkohorten ohne Änderungen angewendet. Die Bilder wurden in die kanonische RAS-Orientierung umorientiert und in Gleitkomma-Arrays umgewandelt. Die räumliche Auflösung in der Ebene wurde mittels linearer Interpolation auf 256 × 256 neu besetzt. Die Durch-Ebenen-Dimension wurde auf 20 Schnitte standardisiert, indem eine zenterbasierte Strategie verwendet wurde: Volumina, die mehr als 20 Schnitte enthielten, wurden mittelbeschnitten, während Bände mit weniger als 20 Scheiben symmetrisch nullgepolstert wurden. Die endgültige Bandgröße betrug 256 × 256 × 20.
Die Intensitätsnormalisierung wurde unabhängig für jedes Volumen mittels z-Score-Normalisierung durchgeführt:

wobei x die Voxelintensität bezeichnet, μ die mittlere Intensität des Volumens und σ die Standardabweichung ist. Wenn σ = 0, wurde keine Normalisierung angewendet, um numerische Instabilität zu verhindern. Verarbeitete Volumes wurden im NIfTI-Format mit einer standardisierten affinen Matrix für die nachgelagerte Deep-Learning-Analyse gespeichert.
3. Entwicklung klinischer Modelle
Strukturierte klinische Variablen wurden verwendet, um Machine-Learning-Modelle zur Outcome-Vorhersage zu entwickeln. Zu den Kandidatenprädiktoren gehörten demografische Merkmale, vaskuläre Risikofaktoren, Schlaganfallätiologie und klinische Schweregradmaße zur Ausgangsbasis. Kontinuierliche Variablen wurden mittels Medianwerten imputiert und standardisiert. Kategoriale Variablen wurden mit der häufigsten Kategorie imputiert und mit One-Hot-Codierung codiert.
Mehrere maschinelle Lernalgorithmen wurden bewertet, darunter logistische Regression, Random Forest, Gradient Boosting, Support Vector Machine (SVM), extreme Gradient Boosting und Light Gradient Boosting Machine-Modelle. Die Modellentwicklung folgte einem geschichteten, fünffachen Cross-Validation-Rahmen innerhalb der Trainings-Validierungskohorte zur Schätzung der Generalisierungsleistung. Endgültige Vorhersagen für die internen und externen Testkohorten wurden durch Durchschnittsvorhersagen aus Modellen generiert, die auf jeder Kreuzvalidierungsfalte trainiert wurden. Während der Kreuzvalidierung oder Modellauswahl wurden keine externen Testproben verwendet.
4. Deep-Learning-Modellarchitektur
Eine dreidimensionale CNN-ViT-Hybridarchitektur wurde zur Ergebnisvorhersage aus multimodalen MRT-Volumina implementiert. Das Netzwerk wurde entwickelt, um lokale räumliche Merkmalsextraktion mit globaler kontextueller Modellierung innerhalb eines einheitlichen Rahmens zu kombinieren. Die Eingangsvolumina bestanden aus mehrkanaligen dreidimensionalen Bildern, die End-to-End-Verfahren verarbeitet wurden.
Die Merkmalsextraktion erfolgte zunächst mit einem hierarchischen dreidimensionalen Faltungsrückgrat, das aus 4 Stufen bestand. Jede Stufe enthielt 2 Faltungsschichten mit einer Kerngröße von 3 × 3 × 3 und Auffüllung von 1 Voxel, gefolgt von Batch-Normalisierung und rectifizierter linearer Einheitsaktivierung. Die räumliche Auflösung wurde schrittweise durch dreidimensionale Max-Pooling-Schichten reduziert, die nach den ersten drei Stufen angewendet wurden, während die Feature-Kanaltiefe auf jeder Ebene erhöht wurde, um höhere semantische Darstellungen zu erfassen. Die Dropout-Regularisierung (Dropout-Rate = 0,1) wurde nach der letzten Faltungsstufe angewendet, um Überanpassung zu reduzieren. Das Faltungsrückgrat transformierte das Eingangsvolumen von Größe C × D × H × W in eine kompakte, hochrangige Merkmalsrepräsentation mit reduzierten räumlichen Dimensionen.
Die resultierende Merkmalskarte wurde in eine Folge von Tokens umgeformt, indem räumliche Dimensionen abgeflacht wurden, sodass die Anzahl der Tokens entsprach:
N = D' x H' x W'
Globale kontextuelle Beziehungen zwischen Token wurden mit Transformer-Encoder-Schichten modelliert, die aus Multi-Head Self-Attention- und Feedforward-Netzwerken bestehen. Selbstaufmerksamkeit wurde berechnet wie folgt:

wobei Q, K und V die Abfrage-, Schlüssel- und Wertmatrizen bezeichnen und d die Einbettungsdimension darstellt. Schichtnormalisierung und Dropout wurden innerhalb jeder Encoder-Schicht angewendet, um die Trainingsstabilität zu verbessern. Das Transformatormodul bestand aus 3 Encoder-Schichten mit 8 Aufmerksamkeitsköpfen und einer Einbettungsdimension von 256.
Nach der Transformer-Codierung wurde die Darstellung des Klassifikationstokens extrahiert und normalisiert. Eine vollständig verbundene lineare Schicht erzeugte eine einzelne Logit-Ausgabe für die binäre Klassifikation.
Die vorgeschlagene hybride CNN-ViT-Architektur wurde bewusst als leichtes und parametereffizientes Modell konzipiert, um Repräsentationskapazität und Überanpassungsrisiko auszubalancieren. Das Modell umfasste 3,79 Millionen trainierbare Parameter (etwa 14,4 MB in fp32-Präzision), darunter 1,38 Millionen im Faltungsrückgrat und 2,37 Millionen im Transformator-Encoder.
5. Bildgebungsmodelltraining
Das Training des Bildgebenden Modells erfolgte mit einem geschichteten fünffachen Kreuzvalidierungsrahmen, um die Ergebnisverteilung über Falten zu erhalten und gleichzeitig die Leistungsabschätzung zu verbessern. Der Datensatz wurde in eine Trainings-Validierungskohorte und eine interne unabhängige Testkohorte unterteilt, wobei stratifizierte Stichproben basierend auf der Ergebnisverteilung verwendet wurden. Innerhalb der Trainings-Validierungskohorte wurde eine geschichtete fünffache Kreuzvalidierung angewendet. Für jede Falte wurde das Bildgebungsmodell mit der Trainingsteilmenge trainiert und mit der entsprechenden Validierungsteilmenge ausgewertet, während die abgehaltene Testkohorte ausschließlich für die endgültige Leistungsbewertung reserviert war.
Die Modelloptimierung erfolgte mit einem Deep-Learning-Framework auf einer GPU-ausgestatteten Workstation mit dem AdamW-Optimierer, Lernrate von 3 × 10⁻5 und einem Gewichtsabfall von 3 × 10⁻4. Die Ausbildung erfolgte mit einer Batchgröße von 8 Personen für bis zu 200 Epochen. Die binäre Kreuzentropie mit Logits wurde als Verlustfunktion verwendet.
Um das Klassenungleichgewicht zu beheben, wurde für jede Faltung ein positiver Klassengewichtungsfaktor berechnet, basierend auf dem Verhältnis negativer zu positiver Stichproben, und in die Verlustfunktion einbezogen. Gradientennorm-Clipping mit einer maximalen Norm von 0,5 wurde angewendet, um die numerische Stabilität während der Optimierung zu verbessern. Automatisches Mixed-Precision-Training wurde ermöglicht, um die Recheneffizienz zu verbessern.
Frühes Stoppen wurde eingeführt, wenn sich die Validierungsleistung über 30 aufeinanderfolgende Epochen nicht um mindestens 1 × 10⁻4 verbesserte. Der leistungsstärkste Modell-Checkpoint aus jedem Falt wurde beibehalten. Nach Abschluss aller Folds wurden Vorhersagen für die internen und externen Testkohorten mit jedem fold-spezifischen Modell erstellt, und die endgültigen Wahrscheinlichkeiten wurden durch Durchschnitt der Vorhersagen über die fünf Modelle hinweg ermittelt, um Ensemble-Ergebnisse zu erzeugen.
6. Multimodales Fusionsmodell
Eine gestapelte Fusionsstrategie wurde implementiert, um bildgebende Vorhersagen mit strukturierten klinischen Informationen zu integrieren. Das Deep-Learning-Bildgebungsmodell und das klinische Vorhersagemodell dienten als Basislerner, und ihre vorhergesagten Wahrscheinlichkeiten wurden als Eingabefunktionen für einen logistischen Regressions-Meta-Lerner verwendet. Zusätzliche Interaktionsmerkmale, einschließlich Produkt- und absoluter Differenz der vorhergesagten Wahrscheinlichkeiten, wurden integriert, um ergänzende Informationen zwischen Bildgebung und klinischer Vorhersagen zu erfassen.
Um Informationslecks zu verhindern, wurde der Meta-Lernende mit außerhalb der Falt vorhergesagten Wahrscheinlichkeiten aus der Trainings-Validierungskohorte trainiert. Kreuzvalidierte Out-of-Fold-Wahrscheinlichkeiten aus den Bildgebungs- und klinischen Modellen wurden anhand der Patientenidentifikator zusammengeführt, um den Meta-Lerner-Trainingsdatensatz zu erstellen. Für die internen und externen Testkohorten wurden entsprechende Testset-Wahrscheinlichkeiten aus den bildgebenden und klinischen Modellen als Eingaben für den trainierten Meta-Lernenden verwendet, um fusionierte Wahrscheinlichkeiten zu erzeugen. Der ausgebildete Meta-Lerner wurde auf beide Testkohorten angewendet, ohne neu anzupassen.
7. Ablationsstudie
Ablationsexperimente wurden durchgeführt, um die Beiträge einzelner MRT-Sequenzen und der Vision Transformer-Komponente unter Verwendung identischer Trainings- und Bewertungssettings wie das primäre Modell zu bewerten. Sequenzablationen umfassten Einzelsequenzmodelle, Leave-One-Sequence-out-Modelle und das vollständige multiparametrische Modell. Um den Beitrag des Vision Transformer-Moduls zu bewerten, wurde die vorgeschlagene hybride Architektur zusätzlich mit einer ausschließlich CNN-basierten Basislinie verglichen, bei der der Transformer-Encoder entfernt und das gleiche Faltungsrückgrat erhalten blieb.
8. Statistische Analyse
Die Modellleistung wurde separat in der internen unabhängigen Testkohorte und der externen Testkohorte unter Verwendung der Fläche unter der Empfänger-Betriebscharakteristiekkurve (AUC) als primäre Diskriminierungsmetrik bewertet. Empfänger-Betriebscharakteristik (ROC)-Kurven wurden mit den von jedem Modell generierten vorhergesagten Wahrscheinlichkeiten konstruiert. Sensitivität, Spezifität und Gesamtgenauigkeit wurden zusätzlich berechnet, um die Klassifikationsleistung zu charakterisieren.
Binäre Ergebnisse wurden anhand von Schwellenwerten generiert, die gemäß Youdens Index bestimmt wurden. Sensitivität, Spezifität und Genauigkeit wurden anschließend am optimalen Schwellenwert berechnet. Alle statistischen Analysen, die Entwicklung von maschinellen Lernmodellen und das Training von Deep-Learning-Modellen wurden mit standardisierten wissenschaftlichen Computing- und Machine-Learning-Softwarepaketen umgesetzt.