$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
אוסטאופורוזיס ושברים שבריריים נלווים עדיין מהווים בעיה מרכזית בבריאות הציבור1. בפרט, המספר העולמי של שברים בירך צפוי להכפיל את עצמו עד 20502. שבריריות העצם נובעת מתהליך איטי ושקט של דה-מינרליזציה ואובדן עצם ללא סימני התראה משמעותיים לפני אירוע השבר השברירי. תקן הזהב הנוכחי לאיתור חולים בסיכון לשבר שברירי הוא ספיגת קרני רנטגן כפולה (DXA), המספקת תמונת רנטגן דו-ממדית ברזולוציה נמוכה עם פיקסל אפורמכויל 3. מתמונה זו ניתן לחלץ את צפיפות המינרלים של העצם השטחית (aBMD ב-g.cm-2) באזורי עניין שונים הקשורים לאתרי השבר השבריריים העיקריים: עמוד השדרה, שורש כף היד והירך. ערך ה-aBMD יורד ככל שקצב השבר השברירי עולה3. יתרה מכך, נורמליזציה של ציון T, ביחס לאוכלוסייה בריאה רגילה, מאפשרת השוואה של חולים שנמדדו עם מכשירים המוצעים על ידי יצרנים שונים. ציון ה-T של DXA הוצע על ידי ארגון הבריאות העולמי כדי להגדיר את אבחון האוסטאופורוזיס בשלושה שלבים: נורמלי (ציון T < -1), אוסטאופני (-1 < ציון T <-2.5) ואוסטאופורוטי (ציון T <-2.5)4.
DXA מציג מספר מגבלות: גודלו, עלותו הגבוהה יחסית, הצורך בחדר ייעודי ויכולתו להבחין בין שבור ללא שבירה, כמו גם זמינותו במדינות רבות, כמו באמריקה הלטינית, שניהם מתונים5. לפיכך, יש צורך בחלופות DXA ככלי סינון להערכת סיכון לשברים שבריריים6. עם זאת, כמה חלופות DXA, כגון טומוגרפיה ממוחשבת כמותית ונגזרותיה7, הדמיית תהודה מגנטית (MRI)8, הן גם מגושמות ואינן זמינות באופן נרחב. אולטרסאונד כמותי (QUS) מציג את הפוטנציאל למכשירי סינון ניידים, חזקים וקלים לשימוש. פותחו מכשירים שונים להערכת עצם קליפת המוח, הקשורים לתדרים שונים הנעים בין כמה קילו-הרץ לכמה מגה-הרץ ומיקום מתמר שונה בשידור, דיפוזיה רטרו9, הד דופק10 ושידור צירי שבו המתמרים מיושרים עם ציר עצם ארוכה כגון רדיוס ושוקה. מכשירים מסוימים מספקים חלופות aBMD11, בעוד שאחרים מספקים פרמטרים קוליים "קלאסיים" כגון מהירויות12 או מקדמי הנחתה9 ואפילו פרמטרים גיאומטריים וחומריים, למשל, עובי קליפת המוח, נקבוביות או התפלגות גודל נקבוביות9. עם זאת, עד היום, QUS עדיין לא הצליח להיות בשימוש נרחב בפרקטיקה הקלינית להערכת עצם, בין היתר בשל היעדר הומוגניזציה בין מכשירים ותלות במפעיל13.
בין טכנולוגיות QUS המוצעות כחלופות DXA, לשידור צירי (AT) יש יתרון שניתן לבצע את המדידה באמה, אתר (i) נגיש בקלות ו-(ii) קרוב לאחד האתרים העיקריים של שברים שבריריים, כלומר שורש כף היד. פרמטר ה-AT המוצע הראשון תלוי במהירות ההתפשטות האולטראסונית בשכבת קליפת המוח, המסומנת במהירות הקול (SOS) או במהירות אות ההגעה הראשון (vFAS), בהתאם לעיבוד האותות והמכשירים, חלקם מסחריים12,14 ואבות טיפוס מעבדתיים אחרים15,16. פרמטר זה הצליח להבחין בין קבוצות חולים עם או בלי שברים שבריריים עם ביצועים דומים ל-BMD במספר מחקרים קליניים מאז סוף שנות ה-9014,15. הוא גם יושם בהצלחה למחקרי אורך רב-מרכזיים, המדגים את היישום הקליני והחוסן שלו12. דיוק ה-vFAS שופר על ידי שילוב שני כיווני ההתפשטות המנוגדים על מנת להפחית את ההטיה עקב הזווית בין הגשושית למשטח העצם16,17. נקודת מבט זו סומנה דו-כיוונית AT (BDAT).
גם אם vFAS הראה עניין קליני, החיסרון העיקרי שלו, בדומה ל-BMD, הוא שהוא משלב תכונות מפתח שונות של עצם קליפת המוח כגון תכונות גיאומטריות וחומריות, מה שהופך את הפרשנות הקלינית שלו לא פשוטה. לכן הוצעה נקודת המבט של הגל המודרך, בהתחשב בפוטנציאל שלה בשל הרגישות העדינה של גלים מונחים לתכונות מוליך הגל. גישה זו צריכה לשלב עיבוד אותות, מידול מוליכי גל ובעיות הפוכות והיא משמשת בעיקר בבדיקות לא הרסניות בהתחשב למשל במוליכי גל מתכתיים, כגון לוחות או צינורות18. לפיכך, מכשיר BDAT מהדור השני פותח צעד אחר צעד מאז 2010, מרפאים מחקי עצם19 לאימות ex vivo 20 ומדידות in vivo 21. המכשיר נבדק בהצלחה במחקרים קליניים בצרפת22, גרמניה23, בריטניה24 וצ'ילה25, והוא הראה תוצאות משופרות מבחינת אחוזי הצלחה ואפליית מטופלים.
מחקר זה נועד לחקור את יכולת השחזור של המכשיר האולטראסוני הנוכחי BDAT. ראשית, המכשיר ופרוטוקול המדידה יפורטו. תוצאות שהתקבלו עם 14 משתתפים ו-3 מפעילים יוצגו וידונו במונחים של סינון אוכלוסייה לאיתור חולים בסיכון לשברים שבריריים.
עקרון מדידה: עיבוד אותות, פרמטרים מעניינים ופרמטרים של איכות
מכשיר השידור הצירי הדו-כיווני (BDAT) מורכב מחלקים שונים, העיקריים הם הגשושית האולטראסוני, המודול האלקטרוני והמחשב. הרשימה המלאה מפורטת בטבלת החומרים ומומחשת באיור 1. להלן מתוארים הפרמטרים המעניינים, פרמטרי איכות המדידה ופרוטוקול המדידה.
vFAS
לאחר שהאותות שנדגמו מתקבלים על ידי המחשב, הם מעובדים לפי שלבים שונים. השלב הראשון מורכב מעיבוד אותות בתחום הזמן, זיהוי FAS באמצעות הפרוטוקול שתואר קודם לכן16,17. לאחר קבלת זמן ההגעה עבור כל מקלט, ניתן לקבוע את מהירות ה-FAS, המסומנת מאוחר יותר vFAS, שהיא הממוצע ההרמוני של המהירויות המתקבלות בשני כיווני ההתפשטות. בשילוב המידע משני כיווני ההתפשטות, ניתן לקבל את זווית הערך בין כיווני הגשושית ומשטח העצם ולהפיק ערך vFAS בלתי מוטה16. זווית דו-כיוונית זו מסומנת מאוחר יותר אלפא ומשמשת כפרמטר לאיכות המדידה. עיבוד זמני זה מאפשר גם להעריך את עובי הרקמה הרכה בין משטח העצם לבדיקה, המסומן ST.Th26.
תמונת ספקטרום גלים מודרכת
השלב השני מורכב מעיבוד אותות בתחום פורייה, בהתחשב בתדרים הזמניים והמרחביים, המסומנים f ו-k. הגישה היא שיטה מבוססת SVD, המאפשרת להפוך את האותות המרחביים-זמניים לפונקציית הנורמה, המצוינת גם בתמונת ספקטרום גל מודרך (GWSI), כפי שמוצג באיור 2 עבור זרוע in vivo 19. השיטה משלבת שני טרנספורמציות פורייה (זמן ומרחב) ופירוק ערך יחיד (SVD), המאפשר לדמיין את קצב הנוכחות באותות המתקבלים (בסולם 0-1) של המצבים המונחים על ידי שכבת העצם בקליפת המוח. ניתן לפרש את ה-GWSI כשיפור של התמרת פורייה המרחבית-זמנית, כאשר כל פיקסל משויך למישור עצמאי של תדר f ומספר גל k. שימו לב שהגישה שופרה על מנת לקחת בחשבון את ההשפעה של הנחתת חומר27 ושינוי עובי ליניארי28.
תשומת לב מיוחדת תינתן לחלק העליון של הספקטרום, הקשור למצב A0, וגם לחלק הנמוך ביותר, הקשור לערכי מהירות הפאזה הגבוהים ביותר, כלומר, גדולים מ-4 מ"מ·μs-1. חלק זה מתאים לאזור העניין 3 (ROI 3)29. הערך הממוצע של ROI 3, שסומן מאוחר יותר lowk, משמש גם כפרמטר איכות. ערך גדול מתאים למוליך גל רגיל, המאפשר השתקפויות גל ברורות בממשקי העצם. אם הערך יורד, זה יכול להיות בגלל מוליך גל לא סדיר או בדיקה שאינה במקומה.
דגם מוליך גל
פיזור הגל המודרך, או הווריאציה של מהירות הפאזה של כל מצב מונחה ביחס לתדר, תלוי הן בתכונות החומר והן בתכונות הגיאומטריות של מוליך הגל. לפיכך, ניתן פוטנציאלית לאחזר מאפיינים אלה באמצעות עיבוד אותות ייעודי, מידול מוליכי גל וסכמות בעיות הפוכות. במקרה של BDAT, מודל מוליך הגל מתאים לצלחת חופשית איזוטרופית רוחבית דו-ממדית, בהתאם לחומר מוליך הגל ופרמטר גיאומטרי אחד, העובי30. חומר העצם בקליפת המוח הומוגני בהתחשב בפרמטרים קבועים למטריצת העצם ונקבוביות משתנה31. לפיכך, הבעיה ההפוכה תלויה בשני פרמטרים, עובי קליפת המוח המסומן (Ct.Th) ונקבוביות קליפת המוח (Ct.Po). ההשפעות של ספיגת חומרים, עקמומיות מוליך הגל והרקמות הרכות שמסביב אינן נלקחות בחשבון במודל, גם אם הן משפיעות על המדידה. עם זאת, משקלם על תוצאת הבעיה ההפוכה לא נמצא כמכריע, כלומר המצבים בשני אזורי העניין העיקריים (A0 והחלק התחתון) אינם משתנים באופן משמעותי על ידי העקמומיות והרקמות הרכות32.
בעיה הפוכה
בתחילה, הבעיה ההפוכה חולקה לשני שלבים: ראשית, לחלץ את פיזור הגל המונחה הניסיוני, ושנית, להשוות למודל מוליך הגל. נקודת מבט זו הוגבלה על ידי תיוג רעש ומצב30,32. לפיכך, הוצעה גישה ייעודית להתגברות על מגבלות אלה כהרחבה של נקודת המבט של פונקציית הנורמה. במקום לשקול כל גל מישורי בנפרד, נלקחים בחשבון רק הגלים המונחים האפשריים המסופקים על ידי מודל מוליך הגל20. זה מוביל לתמונת הבעיה ההפוכה, המתבטאת בתחום הפרמטרים של המודל, כלומר, מישור Ct.th - Ct.Po (איור 2 למטה מימין). המודל המתאים ביותר מקבל את המיקום המקסימלי, בעוד ששיאים משניים בסופו של דבר (המסומנים על ידי תמונות הבעיה ההפוכה עם נקודה אפורה) תואמים לפתרונות דו-משמעיים, המצוינים בהשוואה f-k עם מצבי ניסוי עם קווים אפורים בהירים. כמו קודם, ערך הפיקסלים מנורמל על ידי בנייה ומשקף, במקרה זה, את נוכחותו של מודל מוליך גל מסוים באותות שהתקבלו. הערך המקסימלי (המסומן מקסימום) וההפרש עם המקסימום השני (מסומן הבדל) משמשים גם כפרמטרי איכות.
הבעיה ההפוכה הוצעה במקור לחישוב לא מקוון, כלומר לאחר רכישת האותות, תוך שימוש בערכים המדויקים של מספרי גלי המודל. גישה זו אומתה הן עבור אתרי רדיוס והן עבור אתרי שוק בהתחשב במחקרי ex vivo20,33 ו-in vivo21,34,35. על מנת לכלול חישובים אלה בממשק אדם-מכונה (HMI), הוצעה גרסה משוערת, התואמת ליישום בזמן אמת, תוך שימוש בנקודת מבט מטריצה דלילה36.
vA0
מה-GWSI, ניתן גם לחלץ את המהירות של המצב המונחה האיטי ביותר, הקשור למצב האנטי-סימטרי הראשון A0 של הצלחת החופשית או דגם הכבש33,35. ניתן להעריך את החלק העליון של ספקטרום הגלים המונחים באופן ליניארי, כאשר השיפוע מספק את ערך המהירות vA0 (איור 2 משמאל למטה).
סיכום פרמטרים:
לבסוף, נמדדים ארבעה פרמטרים מעניינים: (i) vFAS: מהירות האות המגיע הראשון (m·s-1); (ii) vA0: מהירות המצב המונחה האיטי ביותר (m·s-1); (iii) Ct.Th: עובי קליפת המוח (מ"מ); ו-(iv) Ct.Po: נקבוביות קליפת המוח (%).
ארבעה פרמטרים איכותיים נחשבים: (i) אלפא: זווית דו-כיוונית (°); (ii) lowk: הערך הממוצע של החלק הנמוך ביותר של ה-GWSI (ערך מנורמל בין 0 ל-1); (iii) מקסימום: מקסימום של פונקציית הבעיה ההפוכה (ערך מנורמל בין 0 ל-1); ו-(iv) diff: ההפרש בין המקסימום הראשון לשני של פונקציית הבעיה ההפוכה (ערך מנורמל בין 0 ל-100).
כל הפרמטרים הללו, כמו גם שתי תמונות ספקטרום הגלים המונחים (כיוון התפשטות אחד תואם) ותמונת הבעיה ההפוכה, מוצגים "בזמן אמת" על ידי ה-HMI, עם קצב פריימים של כ-2 הרץ. דוגמה טיפוסית מומחשת באיור 3. בסעיף הבא מתוארת בפירוט שיטת השימוש בפרמטרים אלה. הרעיון המרכזי הוא שהמפעיל מזיז את הגשושית לאט באתר המדידה, תוך התבוננות קפדנית במשוב המסופק על ידי חלקי הממשק השונים עד למציאת מיקום יציב והתחלת סדרה של 10 רכישות. כאשר מתקבלות לפחות ארבע סדרות עקביות, המדידה מסתיימת ומופק דוח אוטומטי.