$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Deze studie werd goedgekeurd door de Institutionele Beoordelingscommissie van het Shenzhen University General Hospital (IRB-goedkeuring nr. KYLL-2026-077-1). De eis van schriftelijke geïnformeerde toestemming werd opgeheven vanwege het retrospectieve ontwerp van de studie.
1. Studiecohort
Patiënten met acute ischemische beroerte (AIS) die een multiparametrische MRI ondergingen, werden retrospectief geïdentificeerd uit de institutionele klinische database. De algemene studieworkflow, inclusief dataset-partitionering, ontwikkeling van beeldvormingsmodellen, klinische modelconstructie en multimodale fusie, wordt geïllustreerd in Figuur 1.

Figuur 1: Werkstroom van het multimodale uitkomstvoorspellingskader. (A) Dataverwerving. Patiënten met een acute ischemische beroerte (AIS) die voldeden aan vooraf gedefinieerde inclusiecriteria, werden retrospectief geïdentificeerd. De dataset werd verdeeld in een trainings-validatiecohort (n = 250), een interne onafhankelijke testcohort (n = 50) en een externe testcohort (n = 37). Gestratificeerde vijfvoudige kruisvalidatie werd uitgevoerd binnen de training-validatie cohort. (B) Ontwikkeling van beeldvormingsmodellen. Multiparametrische MRI-sequenties, waaronder diffusiegewogen beeldvorming (DWI), schijnbare diffusiecoëfficiënt (ADC) en T2 vloeistofverzwakte inversie recovery (T2-FLAIR), werden verwerkt met behulp van een hybride 3-dimensionale convolutioneel neuraal netwerk-Vision Transformer (CNN-ViT) architectuur voor uitkomstvoorspelling. (C) Ontwikkeling van klinisch model. Gestructureerde klinische variabelen werden gebruikt om machine learning-modellen te trainen voor de voorspelling van functionele uitkomsten na 90 dagen. (D) Multimodale fusiestrategie. Voorspellingen die door de beeldvormings- en klinische modellen werden gegenereerd, werden geïntegreerd met behulp van een gestapelde logistische regressie-meta-learner om uiteindelijke uitkomstvoorspellingen te maken. Klik hier om een grotere versie van deze figuur te bekijken.
MRI-onderzoeken werden uitgevoerd met behulp van klinische MRI-systemen. Voor beeldverwerving werd een speciale fase-array kopspoel gebruikt. Alle sequenties werden in het axiale vlak verkregen met consistente slice-positionering over modaliteiten. Diffusiegewogen beeldvorming (DWI) werd verkregen met behulp van een enkelvoudige spin-echo echo-planaire beeldvorming met herhalingstijd (TR) van 3.000–5.000 ms en echotijd (TE) van 80–90 ms. Diffusiesensitisatie werd toegepast met b-waarden van 0 en 800–1.000 s/mm2 in ten minste 3 orthogonale richtingen. Het gezichtsveld varieerde van 220–240 mm met een matrixgrootte van 128 × 128. De sneeddikte was 5–6 mm met een tussenspleet van 1–1,5 mm. Er werden twee tot vier signaalgemiddelden verkregen.
Schijnbare diffusiecoëfficiënten (ADC) kaarten werden automatisch gegenereerd uit DWI-gegevens op het scannerwerkstation met behulp van mono-exponentiële fitting op basis van de verkregen b-waarden. ADC-waarden werden voxel-gewijs berekend en geëxporteerd voor kwantitatieve analyse. T2 vloeistofverzwakte inversieherstel (T2-FLAIR) beelden werden gemaakt met behulp van een inversieherstelsequentie met TR van 8.000–10.000 ms, TE van 80–140 ms en een inversietijd van 2.200–2.600 ms. Het gezichtsveld varieerde van 220–240 mm met een matrixgrootte van 192 × 192 tot 256 × 256. De sneeddikte was 4–5 mm met een tussenspleet van 1–1,5 mm.
2. Gegevensvoorverwerking
Klinische gegevens werden geïmporteerd uit gestructureerde spreadsheets met case-identificaties, mRS-uitkomstscores, dataset-splitlabels en demografische en klinische variabelen. Zaakidentificaties werden gestandaardiseerd om consistentie met beeldbestandsnamen te waarborgen. Klinische variabelen omvatten demografische kenmerken, de ernst van een beroerte gemeten met de National Institutes of Health Stroke Scale (NIHSS), vasculaire risicofactoren en comorbide aandoeningen die bij opname werden geregistreerd.
De primaire uitkomst was de functionele status 90 dagen na het begin van de beroerte, gemeten met de aangepaste Rankin Scale (mRS). Een gunstig resultaat werd gedefinieerd als mRS ≤ 2 en een ongunstig resultaat als mRS > 2. De dataset werd willekeurig verdeeld in trainingsvalidatie- en interne onafhankelijke testcohorten met behulp van gestratificeerde steekproeven op basis van mRS-verdeling om het uitkomstevenwicht te behouden. Een externe testcohort werd apart verwerkt en werd niet gebruikt tijdens modelontwikkeling. Binnen de training-validatiecohort werd gestratificeerde k-voudige kruisvalidatie toegepast om consistente uitkomstverdelingen over folds te behouden.
Alle DWI-, ADC- en T2-FLAIR-volumes werden vooraf verwerkt vóór modeltraining om ruimtelijke en numerieke consistentie tussen modaliteiten te waarborgen. Dezelfde preprocessing-pijplijn werd toegepast op de trainings-validatie, interne test en externe testcohorten zonder aanpassing. Beelden werden heroriënterend naar canonieke RAS-oriëntatie en omgezet naar floating-point arrays. De ruimtelijke resolutie in het vlak werd opnieuw bemonsterd naar 256 × 256 met behulp van lineaire interpolatie. De door-vlak dimensie werd gestandaardiseerd tot 20 slices met een center-gebaseerde strategie: volumes met meer dan 20 slices werden center-cropped, terwijl volumes met minder dan 20 slices symmetrisch zero-padding waren. De uiteindelijke volumegrootte was 256 × 256 × 20.
Intensiteitsnormalisatie werd onafhankelijk uitgevoerd voor elk volume met behulp van z-score normalisatie:

waarbij x de voxelintensiteit aanduidt, μ de gemiddelde intensiteit van het volume is, en σ de standaarddeviatie. Als σ = 0, werd normalisatie niet toegepast om numerieke instabiliteit te voorkomen. Verwerkte volumes werden opgeslagen in NIfTI-formaat met een gestandaardiseerde affiene matrix voor downstream deep learning-analyse.
3. Ontwikkeling van klinisch model
Gestructureerde klinische variabelen werden gebruikt om machine learning-modellen te ontwikkelen voor uitkomstvoorspelling. Kandidaatvoorspellers waren onder andere demografische kenmerken, vasculaire risicofactoren, beroerteetiologie en klinische ernstmetingen van de basislijn. Continue variabelen werden geïimputeerd met behulp van mediane waarden en gestandaardiseerd. Categorische variabelen werden geïimputeerd met de meest voorkomende categorie en gecodeerd met one-hot encoding.
Meerdere machine learning-algoritmen werden geëvalueerd, waaronder logistische regressie, random forest, gradient boosting, support vector machine (SVM), extreme gradient boosting en light gradient boosting machinemodellen. De ontwikkeling van het model volgde een gestratificeerd vijfvoudig crossvalidatiekader binnen de training-validatie cohort om de generalisatieprestaties te schatten. De definitieve voorspellingen voor de interne en externe testcohorten werden gegenereerd door het gemiddelde van voorspellingen van modellen die op elke crossvalidatiefolder waren getraind. Er werden geen externe testmonsters gebruikt tijdens kruisvalidatie of modelselectie.
4. Deep learning modelarchitectuur
Er werd een driedimensionale CNN-ViT hybride architectuur geïmplementeerd voor uitkomstvoorspelling van multimodale MRI-volumes. Het netwerk is ontworpen om lokale ruimtelijke feature-extractie te combineren met globale contextuele modellering binnen een uniform kader. Invoervolumes bestonden uit meerkanaals driedimensionale beelden die end-to-end werden verwerkt.
Feature-extractie werd aanvankelijk uitgevoerd met behulp van een hiërarchische driedimensionale convolutionele ruggengraat bestaande uit 4 fasen. Elke fase bevatte 2 convolutionele lagen met een kerngrootte van 3 × 3 × 3 en opvulling van 1 voxel, gevolgd door batchnormalisatie en rectificeerde lineaire eenheidsactivatie. De ruimtelijke resolutie werd geleidelijk verminderd met behulp van 3-dimensionale max-pooling-lagen die na de eerste 3 fasen werden toegepast, terwijl de diepte van het feature-kanaal op elk niveau toenam om semantische representaties op hoger niveau vast te leggen. Dropout-regularisatie (dropoutpercentage = 0,1) werd toegepast na de laatste convolutionele fase om overfitting te verminderen. De convolutionele ruggengraat transformeerde het invoervolume van grootte C × D × H × W in een compacte hoog-niveau kentekenrepresentatie met verminderde ruimtelijke afmetingen.
De resulterende feature map werd hervormd tot een reeks tokens door ruimtelijke dimensies af te vlakken zodat het aantal tokens overeenkwam met:
N = D' x H' x W'
Globale contextuele relaties tussen tokens werden gemodelleerd met behulp van Transformer-encoderlagen die bestaan uit multi-head self-attention en feedforward-netwerken. Zelfaandacht werd berekend als:

waarbij Q, K en V respectievelijk query-, key- en waardematrices aanduiden, en d de embedding-dimensie vertegenwoordigt. Laagnormalisatie en dropout werden binnen elke encoderlaag toegepast om de trainingsstabiliteit te verbeteren. De transformatormodule bestond uit 3 encoderlagen met 8 aandachtshoofden en een embedding-dimensie van 256.
Na Transformer-codering werd de representatie die overeenkwam met het classificatietoken geëxtraheerd en genormaliseerd. Een volledig verbonden lineaire laag produceerde één enkele logit-output voor binaire classificatie.
De voorgestelde hybride CNN-ViT-architectuur is bewust ontworpen als een lichtgewicht en parameter-efficiënt model om representatiecapaciteit en overfittingrisico in balans te brengen. Het model bestond uit 3,79 miljoen trainbare parameters (ongeveer 14,4 MB in fp32-precisie), waaronder 1,38 miljoen in de convolutionele ruggengraat en 2,37 miljoen in de transformator-encoder.
5. Beeldvormingsmodeltraining
De training van het beeldvormingsmodel werd uitgevoerd met behulp van een gestratificeerd vijfvoudig cross-validatiekader om de uitkomstverdeling over vouwen te behouden en tegelijkertijd de robuustheid van prestatieschatting te verbeteren. De dataset werd verdeeld in een trainingsvalidatiecohort en een interne onafhankelijke testcohort met behulp van gestratificeerde steekproeven op basis van uitkomstverdeling. Binnen de training-validatie cohort werd gestratificeerde vijfvoudige kruisvalidatie toegepast. Voor elke vouw werd het beeldvormingsmodel getraind met de trainingssubset en geëvalueerd met de bijbehorende validatiesubset, terwijl de vastgehouden testcohort exclusief werd gereserveerd voor de definitieve prestatie-evaluatie.
Modeloptimalisatie werd uitgevoerd met een deep learning-framework op een met GPU uitgeruste werkstation met de AdamW-optimizer, leersnelheid van 3 × 10⁻5 en gewichtsafname van 3 × 10⁻4. De training werd gegeven met een batchgrootte van 8 personen tot wel 200 tijdperken. Binaire kruis-entropie met logits werd gebruikt als verliesfunctie.
Om klasse-onevenwicht aan te pakken, werd voor elke vouw een positieve klassewegingsfactor berekend op basis van de verhouding van negatieve tot positieve monsters en opgenomen in de verliesfunctie. Gradient norm clipping met een maximale norm van 0,5 werd toegepast om de numerieke stabiliteit tijdens optimalisatie te verbeteren. Automatische mixed-precision training werd mogelijk gemaakt om de rekenkundige efficiëntie te verbeteren.
Vroegtijdig stoppen werd ingevoerd wanneer validatieprestaties niet verbeterden met minstens 1 × 10⁻4 over 30 opeenvolgende epochs. Het best presterende modelcontrolepunt van elke vouw bleef behouden. Na voltooiing van alle folds werden voorspellingen voor de interne en externe testcohorten gegenereerd met elk fold-specifiek model, en werden de uiteindelijke kansen verkregen door voorspellingen over de 5 modellen te middelen, zodat ensemble-output werd gegenereerd.
6. Multimodaal fusiemodel
Er werd een gestapelde fusiestrategie geïmplementeerd om beeldvormingsvoorspellingen te integreren met gestructureerde klinische informatie. Het deep learning imagingmodel en het klinische voorspellingsmodel dienden als basislers, en hun voorspelde kansen werden gebruikt als invoerfuncties voor een logistische regressie meta-learner. Extra interactiekenmerken, waaronder het product en het absolute verschil van voorspelde waarschijnlijkheden, werden toegevoegd om aanvullende informatie tussen beeldvorming en klinische voorspellingen vast te leggen.
Om informatielek te voorkomen werd de meta-leerling getraind met behulp van buiten de vouw voorspelde kansen die werden gegenereerd uit de training-validatiecohort. Kruisgevalideerde out-of-fold kansen uit de beeldvormings- en klinische modellen werden samengevoegd op patiëntidentificatie om de meta-leerling trainingsdataset te construeren. Voor de interne en externe testcohorten werden overeenkomstige testset-kansen uit de beeldvormings- en klinische modellen gebruikt als input voor de getrainde meta-leerling om gefuseerde kansen te genereren. De getrainde meta-leerling werd toegepast op beide testcohorten zonder aanpassing.
7. Ablatiestudie
Ablatie-experimenten werden uitgevoerd om de bijdragen van individuele MRI-sequenties en de Vision Transformer-component te beoordelen, waarbij gebruikmaken van identieke trainings- en evaluatie-instellingen als het primaire model. Sequentie-ablaties omvatten enkel-sequentiemodellen, leave-one-sequence-out modellen en het volledige multiparametrische model. Om de bijdrage van de Vision Transformer-module te evalueren, werd de voorgestelde hybride architectuur bovendien vergeleken met een CNN-only baseline waarbij de Transformer-encoder werd verwijderd terwijl dezelfde convolutionele ruggengraat behouden bleef.
8. Statistische analyse
De modelprestaties werden afzonderlijk geëvalueerd in de interne onafhankelijke testcohort en de externe testcohort, waarbij het gebied onder de receiver operating characteristic curve (AUC) de primaire discriminatiemaatstaf was. Receiver operating characteristic (ROC) curves werden geconstrueerd met voorspelde kansen die door elk model zijn gegenereerd. Gevoeligheid, specificiteit en algehele nauwkeurigheid werden bovendien berekend om classificatieprestaties te karakteriseren.
Binaire uitkomsten werden gegenereerd met behulp van drempels die volgens Youden's index werden bepaald. Gevoeligheid, specificiteit en nauwkeurigheid werden vervolgens berekend op de optimale drempel. Alle statistische analyses, ontwikkeling van machine learning-modellen en training van deep learning-modellen werden uitgevoerd met behulp van standaard wetenschappelijke computing- en machine learning-softwarepakketten.