$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
מחקר זה אושר על ידי ועדת הביקורת המוסדית של בית החולים הכללי של אוניברסיטת שנזן (אישור IRB מס' KYLL-2026-077-1). הדרישה להסכמה מדעת בכתב בוטלה בשל עיצוב הרטרוספקטיבי של המחקר.
1. קבוצת לימוד
מטופלים עם שבץ איסכמי אקוטי (AIS) שעברו MRI רב-פרמטרי זוהו בדיעבד מתוך מאגר הנתונים הקליני המוסדי. תהליך העבודה הכולל של המחקר, הכולל חלוקת מערכי נתונים, פיתוח מודלים דימותיים, בניית מודלים קליניים ומיזוג רב-מודלי, מוצג באיור 1.

איור 1: זרימת העבודה של מסגרת חיזוי התוצאות הרב-מודלית. (א) רכישת נתונים. חולים עם שבץ איסכמי אקוטי (AIS) שעומדים בקריטריוני הכללה מוגדרים מראש זוהו בדיעבד. מערך הנתונים חולק לקוהורט אימון-אימות (n = 250), קבוצת בדיקות פנימית עצמאית (n = 50), וקוהורט בדיקה חיצוני (n = 37). בוצע אימות צולב מדורג בחמישה שלבים בתוך קבוצת ההכשרה-אימות. (ב) פיתוח מודלים של דימות. רצפי MRI רב-פרמטריים, כולל דימות משוקלל דיפוזיה (DWI), מקדם דיפוזיה נראה (ADC) ושחזור היפוך נוזלי T2 (T2-FLAIR), עובדו באמצעות ארכיטקטורת רשת עצבית קונבולוציונית תלת-ממדית היברידית עם משמר ראייה (CNN-ViT) לחיזוי התוצאות. (ג) פיתוח מודל קליני. משתנים קליניים מובנים שימשו לאימון מודלים של למידת מכונה לחיזוי תוצאה פונקציונלית של 90 יום. (ד) אסטרטגיית מיזוג רב-מודלית. התחזיות שנוצרו על ידי מודלי ההדמיה והקליניים שולבו באמצעות מטא-לומד רגרסיה לוגיסטית מוערם להפקת תחזיות תוצאות סופיות. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.
בדיקות MRI בוצעו באמצעות מערכות MRI קליניות. סליל ראש ייעודי עם מערך פאזה שימש לרכישת תמונה. כל הרצפים נרכשו במישור הצירי באמצעות מיקום חתך עקבי בין מודאליות. הדמיה משוקללת דיפוזיה (DWI) נרכשה באמצעות רצף דימות אקו-מישורי ספין-אקו חד-פעמי, עם זמן חזרה (TR) של 3,000–5,000 מילישניות וזמן הד (TE) של 80–90 מילישניות. רגישות דיפוזיה יושמה באמצעות ערכי b של 0 ו-800–1,000 s/mm2 לפחות בשלושה כיוונים אורתוגונליים. שדה הראייה נע בין 220–240 מ"מ עם גודל מטריצה של 128 × 128. עובי הפרוסה היה 5–6 מ"מ עם רווח בין הפרוסות של 1–1.5 מ"מ. נרכשו שניים עד ארבעה ממוצעים של אותות.
מפות מקדם דיפוזיה נראה (ADC) נוצרו אוטומטית מנתוני DWI בתחנת העבודה של הסורק באמצעות התאמת מונו-אקספוננציאלית המבוססת על ערכי b שנרכשו. ערכי ADC חושבו לפי ווקסל ויוצאו לניתוח כמותי. תמונות T2-FLAIR (T2-FLAIR) של שחזור היפוך מוחלש בנוזל T2 צולמו באמצעות רצף התאוששות היפוך עם TR של 8,000–10,000 מילישניות, TE של 80–140 מילישניות, וזמן היפוך של 2,200–2,600 מילישניות. שדה הראייה נע בין 220–240 מ"מ עם גודל מטריצה של 192 × 192 עד 256 × 256. עובי הפרוסה היה 4–5 מ"מ עם רווח בין הפרוסות של 1–1.5 מ"מ.
2. קדם-עיבוד נתונים
הנתונים הקליניים יובאו מגיליונות אלקטרוניים מובנים המכילים מזהי מקרים, ציוני תוצאות mRS, תוויות מפוצלות במערכי נתונים, ומשתנים דמוגרפיים וקליניים. מזהי המקרים הותקנו כדי להבטיח עקביות עם שמות הקבצים של הדמיה. המשתנים הקליניים כללו מאפיינים דמוגרפיים, חומרת שבץ שנמדדה באמצעות סולם השבץ של המכון הלאומי לבריאות (NIHSS), גורמי סיכון וסקולריים ומצבים נלווים שנרשמו בעת האשפוז.
התוצאה העיקרית הייתה מצב תפקודי ב-90 ימים לאחר תחילת השבץ, שנמדד באמצעות סולם רנקין (mRS) המותאם למדידה. תוצאה חיובית הוגדרה כ-mRS ≤ 2 ותוצאה שלילית כ-mRS > 2. מערך הנתונים חולק באקראי לקוהורטות אימות-אימון וקבוצות בדיקות פנימיות בלתי תלויות, תוך שימוש בדגימה שכבתית המבוססת על התפלגות mRS לשמירה על איזון התוצאות. קבוצת בדיקה חיצונית עובדה בנפרד ולא שימשה במהלך פיתוח המודל. בתוך קבוצת האימון-אימות, יושמה אימות צולב מדורג של k-fold כדי לשמור על התפלגויות תוצאות עקביות בין הקיפולים.
כל נפחי DWI, ADC ו-T2-FLAIR עובדו מראש לפני אימון המודלים כדי להבטיח עקביות מרחבית ומספרית בין מודאליות. אותו צינור עיבוד מוקדם יושם על קבוצות האימון-וולידציה, בדיקות פנימיות ובדיקות חיצוניות ללא שינוי. התמונות הוסבו לכיוון RAS קנוני והוסבו למערכי נקודה צפה. הרזולוציה המרחבית במישור דוגמה מחדש ל-256 × 256 באמצעות אינטרפולציה ליניארית. ממד המישור הסטנדרטי ל-20 פרוסות באמצעות אסטרטגיה מבוססת מרכז: נפחים המכילים יותר מ-20 פרוסות נחתכו במרכז, בעוד נפחים המכילים פחות מ-20 פרוסות היו מרופדים סימטרית באפס. גודל הכרך הסופי היה 256 × 256 × 20.
נרמול עוצמה בוצע באופן עצמאי לכל נפח באמצעות נרמול ציון z:

כאשר x מציין עוצמת ווקסלים, μ היא עוצמת הנפח הממוצעת, ו-σ היא סטיית התקן (סטיית תקן). אם σ = 0, לא יושמה נרמול כדי למנוע חוסר יציבות נומרית. נפחי עיבוד נשמרו בפורמט NIfTI עם מטריצה אפינית סטנדרטית לניתוח למידה עמוקה במורד הזרם.
3. פיתוח מודל קליני
משתנים קליניים מובנים שימשו לפיתוח מודלים של למידת מכונה לחיזוי תוצאות. המנבאים המועמדים כללו מאפיינים דמוגרפיים, גורמי סיכון וסקולריים, אטיולוגיה לשבץ ומדדי חומרה קלינית בסיסיים. משתנים רציפים הוערכו באמצעות ערכי חציון וסטנדרטים. משתנים קטגוריים הונחו באמצעות הקטגוריה הנפוצה ביותר וקודדו באמצעות קידוד חד-חם.
הוערכו מספר אלגוריתמים של למידת מכונה, כולל רגרסיה לוגיסטית, יער אקראי, הגברת גרדיאנט, מכונות וקטוריות תמיכה (SVM), הגברת גרדיאנט קיצוני ומודלים של מכונות להגברת גרדיאנט אור. פיתוח מודלים התבצע לפי מסגרת אימות צולב מדורגת בת חמישה שלבים בתוך קבוצת ההכשרה-אימות כדי להעריך ביצועי הכללה. תחזיות סופיות לקוהורטות הבדיקה הפנימיות והחיצוניות נוצרו על ידי ממוצע תחזיות ממודלים שאומנו בכל קיפול אימות צולב. לא נעשה שימוש בדגימות בדיקה חיצוניות במהלך אימות צולב או בחירת מודל.
4. ארכיטקטורת מודלים של למידה עמוקה
יושמה ארכיטקטורה היברידית תלת-ממדית של CNN-ViT לחיזוי תוצאות מנפחי MRI רב-מודליים. הרשת תוכננה לשלב חילוץ תכונות מרחבי מקומי עם מידול הקשר גלובלי במסגרת אחידה. נפחי הקלט כללו תמונות תלת-ממדיות רב-ערוציות שעובדו מקצה לקצה.
חילוץ התכונות בוצע בתחילה באמצעות עמוד שדרה קונבולוציוני תלת-ממדי היררכי המורכב מארבעה שלבים. כל שלב כלל 2 שכבות קונבולוציוניות בגודל גרעין של 3 × 3 × 3 וריפוי של ווקסל אחד, ואחריו נרמול אצווה והפעלה ליניארית של יחידה. הרזולוציה המרחבית הופחתה בהדרגה באמצעות שכבות מקסימום איגוד תלת-ממדיות שהוחלו לאחר שלושת השלבים הראשונים, בעוד שעומק ערוץ התכונה גדל בכל רמה כדי ללכוד ייצוגים סמנטיים ברמה גבוהה יותר. רגולריזציה של דרופאאוט (קצב נפילות = 0.1) יושמה לאחר שלב הקונבולוציוני הסופי כדי להפחית התאמת יתר. עמוד השדרה הקונבולוציוני הפך את נפח הקלט בגודל C × D × H × W לייצוג תכונות קומפקטי ברמה גבוהה עם ממדים מרחביים מוקטנים.
מפת התכונות שהתקבלה עוצבה מחדש לרצף של טוקנים על ידי השטחת ממדים מרחביים כך שמספר הטוקנים תואם ל:
N = D' x H' x W'
מערכות יחסים קונטקסטואליות גלובליות בין טוקנים דומנו באמצעות שכבות מקודד טרנספורמר המורכבות מרשתות ריבוי ראשים של תשומת לב עצמית ורשתות הזנה קדימה. תשומת הלב העצמית חושבה כך:

כאשר Q, K ו-V מציינים מטריצות שאילתה, מפתח וערכים, בהתאמה, ו-d מייצג ממד הטמעה. נורמליזציה של שכבות ונפילות יושמו בכל שכבת מקודד לשיפור יציבות האימון. מודול הטרנספורמר כלל 3 שכבות מקודד עם 8 ראשי תשומת לב ומימד הטמעה של 256.
לאחר קידוד טרנספורמר, הייצוג המתאים לטוקן הסיווג הופק וננורמל. שכבה ליניארית מחוברת לחלוטין יצרה פלט לוגיט יחיד לסיווג בינארי.
הארכיטקטורה ההיברידית המוצעת של CNN-ViT תוכננה במכוון כמודל קל משקל ויעיל מבחינת פרמטרים, כדי לאזן בין קיבולת ייצוג לסיכוני התאמת יתר. המודל כלל 3.79 מיליון פרמטרים ניתנים לאימון (כ-14.4 מגה-בייט בדיוק fp32), כולל 1.38 מיליון בעמוד השדרה הקונבולוציוני ו-2.37 מיליון במקודד הטרנספורמר.
5. אימון מודל הדמיה
הכשרה של מודל ההדמיה בוצעה באמצעות מסגרת אימות צולבת חמש-שכבתית לשמירה על התפלגות התוצאה בין הקיפולים תוך שיפור עמידות הערכת הביצועים. מערך הנתונים חולק לקוהורט אימון-אימות וקוהורט בדיקות פנימי עצמאי תוך שימוש בדגימה שכבתית המבוססת על התפלגות התוצאות. בתוך קבוצת האימון-אימות, בוצעה אימות צולב מדורג בחמישה שכבות. לכל קיפול, מודל ההדמיה אומן באמצעות תת-קבוצת האימון והוערך באמצעות תת-הקבוצה המתאימה, בעוד שקבוצת הבדיקה שהוחזקה הוקדשה אך ורק להערכת ביצועים סופית.
אופטימיזציה של מודלים בוצעה באמצעות מסגרת למידה עמוקה על תחנת עבודה מצוידת ב-GPU עם אופטימייזר AdamW, קצב למידה של 3 × 10⁻5, ודעיכת משקל של 3 × 10⁻4. ההכשרה התבצעה באמצעות גודל אצווה של 8 יחידות עד 200 אפוקים. אנטרופיה חוצה בינארית עם לוגיטים שימשה כפונקציית הפסד.
כדי להתמודד עם חוסר איזון במחלקה, חושב גורם משקל חיובי-מחלקה לכל קיפול בהתבסס על היחס בין מדגמים שליליים לחיוביים ושולב בפונקציית ההפסד. חיתוך נורמות גרדיאנט עם נורמה מקסימלית של 0.5 יושם לשיפור היציבות המספרית במהלך האופטימיזציה. הכשרה אוטומטית עם דיוק מעורב אפשרה לשפר את היעילות החישובית.
עצירה מוקדמת יושמה כאשר ביצועי האימות לא השתפרו לפחות ב-1 × 10⁻4 במשך 30 תקופות רצופות. נקודת הביקורת המודל הטובה ביותר מכל קיפול נשמרה. לאחר סיום כל הקיפולים, נוצרו תחזיות עבור קוהורטות הבדיקה הפנימיות והחיצוניות באמצעות כל מודל ספציפי לקיפול, וההסתברויות הסופיות התקבלו על ידי ממוצע תחזיות בין חמשת המודלים ליצירת פלטי אנסמבל.
6. מודל היתוך רב-מודלי
יושמה אסטרטגיית מיזוג מוערמת לשילוב תחזיות שמקורם בהדמיה עם מידע קליני מובנה. מודל הדמיית הלמידה העמוקה ומודל החיזוי הקליני שימשו כבסיס ללמידה, וההסתברויות החזויות שימשו כמאפייני קלט למטה-לומד רגרסיה לוגיסטית. תכונות אינטראקציה נוספות, כולל התוצר וההפרש המוחלט בין ההסתברויות החזויות, שולבו כדי ללכוד מידע משלים בין הדמיה לתחזיות קליניות.
כדי למנוע דליפת מידע, המטה-לומד עבר הכשרה באמצעות הסתברויות חזויות מחוץ לקפל שנוצרו מקבוצת האימון-אימות. הסתברויות מחוץ לאזור מאומתות ממודלי ההדמיה והקליניים אוחדו על ידי מזהה המטופל לבניית מערך הנתונים לאימון המטה-לומדים. עבור קבוצות הבדיקה הפנימיות והחיצוניות, השתמשו בהסתברויות מקבוצה מתאימה מהמודלים של ההדמיה והקליניים כקלטים למטה-לומד המיומן ליצירת הסתברויות מאוחדות. הלומד המטא-מיומן הוחל על שתי קבוצות המבחן ללא התאמה.
7. מחקר אבלציה
נערכו ניסויי אבלציה להערכת תרומתם של רצפי MRI בודדים ורכיב Vision Transformer תוך שימוש באותן הגדרות אימון והערכה כמודל הראשי. אבלציות רצף כללו מודלים של רצף יחיד, מודלים של השארת רצף אחד החוצה, ומודל רב-פרמטרי מלא. כדי להעריך את תרומת מודול Vision Transformor, הארכיטקטורה ההיברידית המוצעת הושוותה בנוסף לבסיס CNN בלבד, שבו מקודד הטרנספורמר הוסר תוך שמירה על אותו עמוד שדרה קונבולוציונלי.
8. ניתוח סטטיסטי
ביצועי המודל הוערכו בנפרד בקוהורט הבדיקות הפנימי העצמאי ובקוהורט הבדיקה החיצוני, תוך שימוש בשטח מתחת לעקומת מאפיין ההפעלה של המקלט (AUC) כמדד ההבחנה הראשוני. עקומות מאפיין הפעלה של מקלט (ROC) נבנו באמצעות הסתברויות חזויות שנוצרו על ידי כל מודל. רגישות, ספציפיות ודיוק כולל חושבו גם כדי לאפיין את ביצועי הסיווג.
תוצאות בינאריות נוצרו באמצעות ספים שנקבעו לפי מדד יודן. רגישות, ספציפיות ודיוק חושבו לאחר מכן בסף האופטימלי. כל הניתוחים הסטטיסטיים, פיתוח מודלים של למידת מכונה והדרכת מודלים ללמידה עמוקה יושמו באמצעות חבילות תוכנה סטנדרטיות למחשוב מדעי ולמידת מכונה.