Method Article
כאן, אנו מציגים בפירוט את פרוטוקול המדידה של המכשיר האולטראסוני של שידור צירי דו-כיווני (BDAT) ובודקים אותו במחקר שחזור, בהתחשב ב-14 משתתפים בריאים ו-3 מפעילים. המהימנות, שנמדדה עם מקדמי המתאם התוך-מחלקתיים (ICC), הייתה טובה עד מצוינת עבור ארבעה פרמטרים מעניינים.
שברי שבירות עדיין מהווים נטל בריאותי עולמי בהקשר של הזדקנות האוכלוסייה. בפרט, המספר העולמי של שברים בירך צפוי להכפיל את עצמו בין 2020 ל-2050. לכן, חיוני לאתר חולים בסיכון לשבר שברירי בקנה מידה של אוכלוסייה. תקן הזהב הנוכחי הוא ספיגת קרני רנטגן כפולה (DXA), המספקת את צפיפות המינרלים של העצם השטחית (aBMD). מכשירים קוליים, בדרך כלל ניידים וזולים יותר ממכשירי רנטגן, מייצגים חלופות DXA מעניינות ככלי סינון. עם זאת, תלות במפעיל מוכרת בדרך כלל כחיסרון העיקרי שלהם. במחקר זה מוצג בפירוט פרוטוקול המדידה של מכשיר אולטרסאונד שידור צירי דו-כיווני (BDAT). הגשושית האולטראסונית הייעודית ממוקמת ברדיוס הדיסטלי של שליש של האמה הלא דומיננטית באמצעות ג'ל צימוד קונבנציונאלי. הממשק המונחה מספק בזמן אמת (כ-2 הרץ) ארבעה פרמטרים מעניינים: מהירויות האות המגיע הראשון (vFAS) ומצב A0 (vA0), עובי קליפת המוח (Ct.Th) ונקבוביות (Ct.Po), כמו גם ארבעה פרמטרים של איכות. המפעיל מזיז את הגשושית באיטיות באתר המדידה, תוך התבוננות מדוקדקת במשוב שמספק הממשק עד למציאת מיקום יציב והתחלת סדרה של 10 רכישות. כאשר מתקבלות לפחות ארבע סדרות עקביות, המדידה מסתיימת ומופק דוח אוטומטי. המדידה אורכת בדרך כלל כ-5 דקות. כדי לקבוע את החוסן של פרוטוקול זה, נערך מחקר שחזור בקרב 3 מפעילים (מומחה אחד ושני טירונים) ו-14 משתתפים בריאים (6 נשים, 8 גברים, בגילאי 21-53). מקדמי המתאם התוך-מחלקתיים (ICC) נמצאו טובים עבור vA0 (0.76), Ct.Po (0.80) או מצוינים עבור Ct.Th (0.87) ו-vFAS (0.91). סטיות התקן נמצאו פחות מ-10% מכלל הטווחים בפרקטיקה הקלינית.
אוסטאופורוזיס ושברים שבריריים נלווים עדיין מהווים בעיה מרכזית בבריאות הציבור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 רכישות. כאשר מתקבלות לפחות ארבע סדרות עקביות, המדידה מסתיימת ומופק דוח אוטומטי.
המחקר אושר על ידי הוועדה האתית של אוניברסיטת ולפראיסו, צ'ילה, תחת פרוטוקול מספר CEC213-20. ניתנה הסכמה מדעת בכתב למשתתפים. נערך ראיון טלפוני כדי לקבוע הכללה/אי הכללה של המשתתפים. המחקר נרשם תחת האסמכתא הבאה: NCT05424536.
1. הגדרת המכשיר
2. התקנת משתתף
3. מחפש עמדה יציבה
הערה: ה-HMI מציג ארבעה פרמטרים מעניינים: שתי מהירויות, vFAS ו-vA0, ושני ערכי בעיה הפוכים, עובי קליפת המוח (Ct.Th) ונקבוביות קליפת המוח (Ct.Po). ה-HMI מציג גם ארבעה פרמטרי איכות המסומנים אלפא, lowk, max ו-diff. פרמטרים אלה מתוארים בפירוט בהקדמה.
4. רכישת נתונים
מחקר שחזור נערך בהתבוננות ב-3 מפעילים (מומחה אחד, שני טירונים) ו-14 משתתפים בריאים (6 נשים, 8 גברים, בגילאי 21-53). מפעילים מתחילים הוכשרו במשך כ-3 שעות להבין ולתרגל את פרוטוקול הרכישה. לאחר מכן, המשתתפים נמדדו במהלך שבועיים באוגוסט 2023. כל מדידה בוצעה באופן עצמאי. כל המפעילים היו עיוורים, כלומר, מפעיל אחד לא ידע את התוצאות שהשיגו השניים האחרים.
יכולת חזרה תוך-אופרטור
איור 5 ממחיש את יכולת החזרה התוך-אופרטור עבור מקרה עקבי: 4 סדרות של 10 רכישות שנרכשו על משתתף אחד ומפעיל אחד. העמודה הראשונה מתאימה לפונקציות הבעיה ההפוכה, בעוד שהעמודה השנייה מציגה את פיזור המצב המונחה הניסיוני בהשוואה למודל המתאים ביותר. כל קו איור מתאים לסדרה מוצלחת. מספר הרכישות שנשמרו מופיע בכותרת: 7 מתוך 10 עבור העונה הראשונה ו-10 מתוך 10 עבור הסדרות הבאות. ממוצעים וסטיות תקן וארבעת הפרמטרים (vFAS, vA0, Ct.Th ו-Ct.Po) מסופקים עבור כל סדרה. בנוסף, מוצגים גם פרמטרי איכות: זווית דו-כיוונית (אלפא), מקסימום של פונקציית הבעיה ההפוכה (מקסימום), ההפרש המוחלט עם המקסימום השני (diff), והערך הממוצע של החלק הנמוך ביותר של ה-GWSI (lowk).
סטיות התקן התוך-סדרתיות נמוכות, כ-0.02 מ"מ עבור עובי קליפת המוח, פחות מ-0.5% עבור נקבוביות קליפת המוח, ופחות מ-20 m∙s-1 עבור שתי המהירויות, מה שמצביע על כך שמיקומי בדיקת ההגעה היו יציבים. לאחר מכן, ניתן לראות שהערכים הממוצעים המתקבלים עבור כל סדרה קרובים מאוד, במיוחד עבור ערכי העובי הנעים בין 3.4 מ"מ ל-3.5 מ"מ וערכי ה-vFAS הנעים בין 4040 m∙s-1 ל-4070 m∙s-1. שימו לב שהפרש של 40 m∙s-1 מתאים להבדל של 1% לערך ממוצע של 4000 m∙s-1. וריאציות גדולות יותר נצפות עבור נקבוביות קליפת המוח, הנעות בין 8% ל-12%, ומהירות vA0, הנעה בין 1700 m∙s-1 ל-1740 m∙s-1. במקרה עקבי זה, כמעט כל הרכישות עקביות, כלומר קרובות זו לזו. אין כמעט עמימות בתוצאה הסופית של ארבעת הפרמטרים המעניינים.
מקרה שני מתואר באיור 6 עבור מקרה פחות רגיל. במקרה זה, שלושה פרמטרים, vFAS, vA0 ו-Ct.Po יציבים מאוד, עם ערכים בהתאמה כ-4120 m∙s-1, 1700 m∙s-1 ו-10%. המקרה של עובי קליפת המוח קשה יותר מכיוון ששני פתרונות מעורפלים, 2.6 מ"מ ו-3.5 מ"מ, נצפים בהסכמה עם ערך הפרש קטן (הפרש של פחות מ-0.5%) בין שני המקסימום הראשונים של פונקציית הבעיה ההפוכה. במקרה הרגיל הקודם, הבדל זה נע בין 1% ל-3%. העמימות מוסרת על ידי ניתוח מומחים, במקרה זה, תוך התבוננות בהסכמה בין מצבים מונחים ניסיוניים ותיאורטיים (עמודה ימנית). במקרה של העובי הנמוך ביותר, ההסכמה טובה יותר בחלק הנמוך ביותר של הספקטרום (שתי סדרות ראשונות). עבור שתי הסדרות האחרונות, יש מצב תיאורטי עם מעט מאוד נקודות ניסוי, בסביבות 0.5 מגה-הרץ, מה שמצביע על הסכמה גרועה יותר בהשוואה לסדרה הקודמת. יתר על כן, פרמטר ההבדל (0.1%) קטן מהערכים של שתי הסדרות הראשונות (0.4% ו-0.2%). במקרה זה, הבחירה בתמיסה השמורה (2.6 מ"מ) עדיין אינה אוטומטית, ועדיין יש צורך במומחה. עם זאת, שלושת המפעילים התמודדו עם בעיות דומות ובחרו בפתרונות דומים, קרוב ל-2.6 מ"מ.
אמינות מחוץ למפעיל
כל התוצאות עבור 4 הפרמטרים המעניינים שהתקבלו על ידי 3 האופרטורים עם 14 המשתתפים מוצגות באיור 7. מקדמי המתאם התוך-מחלקתיים (ICC) חושבו בעקבות הנוסחה וקוד Matlab שפורסם בעבר 37,38,39. ICC משמש בדרך כלל להערכת האמינות של סולמות מדידה, במיוחד עבור יישומים ביו-רפואיים. ערכי ICC הנעים בין 0.75 ל-0.9 קשורים בדרך כלל לאמינות טובה, בעוד שערכי ICC מעל 0.9 נחשבים לאמינות מצוינת. הערך הנמוך ביותר עבור ICC (0.76) התקבל עבור הפרמטר vA0. סטיית התקן הייתה שווה ל-17 m∙s-1, שהם כ-7% מטווח המדידה בסדר גודל של 250 m∙s-1. ערכים דומים נצפו עבור Ct.Po עם ICC שווה ל-0.80, וסטיית תקן של 1.1%, כ-10% מהטווח. אמינות מצוינת (ICC כ-0.9) התקבלה עבור שני הפרמטרים האחרים, Ct.Th ו-Ct.Po, עם סטיית תקן נמוכה מ-10% מהטווח.
איור 1: חלקים שונים של ההתקן האולטרה-סוני עם שידור צירי דו-כיווני (BDAT). אב הטיפוס כולל בידוד אלקטרוני (1), מתג דוושה (2), שני בדיקות (3.1 ו-3.2), מודול אלקטרוני (4), מחשב (5) וסרגל (6). אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
איור 2: רכישה טיפוסית על האמה in vivo. שתי פונקציות הנורמה (אחת לכל כיוון התפשטות) מוצגות בשורה העליונה של התמונה. הם מסומנים גם תמונת ספקטרום גל מודרך (GWSI). מהמקסימום של תמונות אלה, ניתן לחלץ את המצבים המונחים הניסיוניים (נקודות כחולות ואדומות) בהשוואה למודל המתאים ביותר (תת-תמונה משמאל למטה). המודל המתאים ביותר מסווג על ידי שני ערכים, עובי קליפת המוח (Ct.Th) ונקבוביות (Ct.Po), המתאימים למיקום המקסימלי של פונקציית הבעיה ההפוכה (תת-תמונה ימנית תחתונה). הערכים שלהם מוצגים בכותרת של כל לוח. התאמת vA0 מוצגת כקו מקווקו (משמאל). תמונות ספקטרום הגל המודרך ותמונת הבעיה ההפוכה מנורמלות (כלומר, ערך הפיקסלים נע בין 0 ל-1) על ידי בנייה19,20. הערכים של שתי המהירויות הנמדדות, vFAS ו-vA0, מצוינים בכותרת של תמונת המשנה הימנית העליונה. אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
איור 3: ממשק מכונה אנושית (HMI). ה-HMI מציג במציאות את שני ה-GWSIs (אחד לכל כיוון התפשטות), את תמונת הבעיה ההפוכה, את הפרמטרים המעניינים ואת פרמטרי האיכות. אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
איור 4: דוגמה לדוח אוטומטי. הדוח ציין את נתוני המשתתף והמפעיל וכן את הערכים הסופיים של פרמטרי העניין ופרמטרי האיכות. אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
איור 5: דוגמה לדוח השני עבור מקרה עקבי, זהה לדוח האוטומטי המוצג באיור 4. האיור מציג 4 סדרות של משתתף אחד ומפעיל אחד: תמונות בעיות הפוכות (עמודה שמאלית) ומספרי גל ניסיוניים בהשוואה למודל המתאים ביותר (עמודה ימנית). הערכים המוצגים תואמים את הממוצע וסטיית התקן של הרכישות שנשמרו בסדרת 10 הרכישות. מספר הרכישות שנשמרו מצוין בכותרת העמודה הימנית, למשל, (7/10) עבור הסדרה הראשונה. אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
איור 6: דוגמה לדוח השני למקרה מעורפל של עובי קליפת המוח. האיור מציג 4 סדרות של משתתף אחד ומפעיל אחד: תמונות בעיות הפוכות (עמודה שמאלית) ומספרי גל ניסיוניים בהשוואה למודל המתאים ביותר (עמודה ימנית). הערכים המוצגים תואמים את הממוצע וסטיית התקן על פני הרכישות שנשמרו בסדרת 10 הרכישות. מספר הרכישות שנשמרו מצוין בכותרת העמודה הימנית, למשל, (10/10) עבור הסדרה הראשונה. אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
איור 7: אמינות בין-מפעילים. התוצאות מוצגות עבור 4 הפרמטרים המעניינים שהתקבלו על ידי 3 האופרטורים עם 14 משתתפים. הערכים המתקבלים על ידי 3 האופרטורים (ציר y) מושווים לערך הממוצע של 3 האופרטורים (ציר x). מקדמי המתאם התוך-מחלקה (ICC) מצוינים בכותרות. סטיית התקן σ ומקדם השונות CV מצוינים גם הם. אנא לחץ כאן לצפייה בגרסה גדולה יותר של איור זה.
הנקודה הקריטית להצלחת המדידה היא מיקום הבדיקה הנכון ביחס לעצם. המיקום ביחס למשטח העצם נפתר בהצלחה על ידי שילוב עם הממשק המנחה את התיקון הדו-כיווני16,17. הוכח כי ללא תיקון זה, ניתן להגיע לשגיאה של כמה אחוזים עבור פרמטרי vFAS16. תיקון זה נמצא חיוני להבחנה מדויקת בין חולים שכן כל טווח המדידה (כ-3800-4200 m∙s-1) הוא כ-10% מהערך הממוצע17, כ-4000 m∙s-1. יכולת השחזור שנמצאה במחקר זה (סטיית תקן של 24 m∙s-1) הייתה דומה למחקר הקודם17, והצביעה על סטיית תקן של 20 m∙s-1 ומקדם שונות (CV) של כ-0.5%. מקדם שונות דומה (0.5%) נמצא בתדר נמוך יותר, כלומר 0.3 מגה-הרץ40. שים לב שהאמינות המצוינת של פרמטר vFAS (ICC = 0.91) נמצאה גבוהה מזו שהושגה לאחרונה עם התקן AT אחר באוכלוסייה מבוגרת יותר (ICC = 0.77)41. יש לציין כי האוכלוסייה במקרה זה הייתה לא רק מבוגרת יותר אלא גם כל הנשים וכמה נבדקים עם סימנים קליניים אחרים של נקבוביות עצם מוגברת41.
האתגר השני הוא גם למקם את הגשושית בצורה נכונה ביחס לציר העצם ולשני הממשקים על מנת להשיג את המצבים המודרכים, בפרט אלה של מהירות פאזה גבוהה או מספר גל נמוך. מצבים אלה קרובים לתהודה וקשורים לנתיבי השתקפות מרובים. אם היישור אינו נכון, הגלים יתפזרו מחוץ לטווח המדידה המרחבי של הגשושית. נהפוך הוא, אם הגשושית ממוקמת נכון, מצבי מהירות פאזה גבוהה אלה יופיעו כענפים רציפים ב-ROI3. עבור הגרסה הראשונה של ה-HMI הנוכחי, יישור זה נפתר על ידי התבוננות ב-GWSI20 בזמן אמת. עם זאת, גישה זו נמצאה לא מספקת כדי להגיע למדידה חזקה: שיעור הכישלון נמצא כ-20% במחקר הפיילוט הקליני20. הכללת פרמטרים איכותיים בזמן אמת מאפשרת להפחית את הכישלון לכ-10% במחקר קליני שני23. לאחר מכן, הכללת תמונת הבעיה ההפוכה ב-HMI "בזמן אמת" מאפשרת שיפור נוסף, עם שיעור כשל נוכחי למדידת רדיוס של כ-5%25. שימו לב שהכשל הראשוני במכשיר הראשון שמדד vFAS בלבד נמצא כ-15%17.
מהימנות על עובי קליפת המוח נמצאה דומה למחקר קודם22. בגרסה קודמת של אותו אב טיפוס, נמצא גם ICC קרוב ל-0.9, עם סטיית תקן של כ-0.1 מ"מ ו-CV של כ-3%. עם זאת, נצפה שיפור ברור בנקבוביות קליפת המוח: ICC עלה מ-0.622 ל-0.8, וסטיית התקן ירדה מ-1.5% ל-1%. הפרמטר הקשה ביותר הוא vA0 בשל הקרבה של מהירותו למהירות הרקמה הרכה, כ-1500 m∙s-1. אם הגשושית מיושרת כהלכה, מצב A0 נראה ייחודי ורציף. להיפך, הוא נראה לא רציף ו/או מרובה עקב צימוד עם הרקמות הרכות שמסביב. השפעה זו חזקה במיוחד עבור אנשים עם עובי רקמה רכה קטן (פחות מ-4 מ"מ) או גדול (יותר מ-10 מ"מ) (מסומן ST.Th).
כפי שנחשף קודם, נקודת המפתח היא למצוא את מיקום הבדיקה הנכון. עם זאת, הפוזיציה צריכה להיות לא רק נכונה אלא גם יציבה על מנת לבצע ממוצע וסטיית תקן עבור סדרה של 10 רכישות. ברוב המקרים, מציאת אורווה אורכת פחות מדקה אחת, והמדידה המלאה של המטופל תימשך כ-5 דקות. גם אם הפרוטוקול הסטנדרטי המתואר במחקר זה מתאים היטב לרוב המטופלים, יש אנשים שקשה יותר למדוד אותם; קשה אך אפשרי למצוא עמדות נכונות, אך כמעט בלתי אפשרי למצוא עמדה יציבה. במקרה זה, המפעיל יכול לבחור להקליט סדרות ארוכות יותר, עד 200 רכישות. הרכישות הטובות ביותר נקבעות מאוחר יותר במצב לא מקוון באמצעות פרמטרי איכות. סינון זה אמור להיות מיושם בזמן אמת בעתיד. בפועל, המפעיל רושם יותר סדרות ו/או רכישות מאשר 4 סדרות של 10 רכישות של פרוטוקול המצב האידיאלי. עם זאת, זמן המדידה הרגיל נשאר כ-5 דקות עבור אתר אחד, האמה או הרגל. אם לאחר מספר דקות ה- HMI אינו מזהה מיקום נכון, המדידה נעצרת ונחשבת לכישלון. עם המכשיר והפרוטוקול הנוכחיים, שיעור הכשל נמצא נמוך מ-5%25.
המכשיר הנוכחי מתמודד עם מגבלות שונות:
(ט)גודל ומשקל: מכשיר ה-BDAT הנוכחי נייד: הוא מתאים למזוודות סטנדרטיות ומשקלו כ-25 ק"ג. עם זאת, משקל זה גדול ביחס למכשירים האולטראסונים החדשים ביותר. ניתן לשקול עיצוב אלקטרוני חדש, כאשר הגשושית ועיבוד האותות יישארו זהים. עם זאת, ניתן להזיז את המכשיר הנוכחי, במיוחד עבור חולים במיטה או בבית עם ניידות מוגבלת.
(ב)מהירות רכישה: קצב הפריימים הנוכחי הוא כ-2-4 הרץ, כלומר חקר המדידה איטי בהשוואה לרכישה בזמן אמת בפועל, כלומר גבוה מ-25 הרץ. ניתן לשפר זאת בעתיד, תוך התחשבות במחשבים מהירים יותר, ניתוח נתונים מהיר יותר והעברה בין האלקטרוניקה למחשב. עלייה במהירות הרכישה תשפר את הידידותיות למשתמש של המדידה, במיוחד את החיפוש אחר מיקום נכון של הגשושית.
(ג)עובי רקמות רכות: הגישה הנוכחית מוגבלת על ידי שכבת רקמה רכה גדולה, בדרך כלל עדיפה על 10 מ"מ. במקרה זה, האות הראשון שמגיע מקושר לנתיב הרקמות הרכות ולא לעצם קליפת המוח. לפיכך, לא ניתן להשתמש ב-vFAS ובזווית הדו-כיוונית הקשורה אליו. כמו כן, קשה מאוד למדוד vA0 עבור שכבות גדולות של רקמות רכות. ללא שתי מהירויות אלו, לא ניתן לבצע את הבעיה ההפוכה. בעתיד ניתן יהיה ליישם שיטות אחרות לתיקון דו-כיווני באמצעות טכניקות הדמיה למשל. חולים עם שכבות גדולות של רקמות רכות קשורים בדרך כלל להשמנת יתר ומדד מסת גוף (BMI) גדול מ-30 ק"ג-מ"ר.
(ד)סדירות מוליך גל: גישת הבעיה ההפוכה מניחה מוליך גל רגיל עם מסלולי התפשטות מרובים. עבור חולים אוסטאופורוטיים, הממשק הפנימי בקליפת המוח יכול להיות לא סדיר ולכן מרמז על תמונות ספקטרום גלים מונחות גרועות, במיוחד בחלק התחתון. חולים אלה קשורים בדרך כלל לעמימות גבוהה בתמיסה. אם הרקמה הרכה או המיקום הלקוי אינם יכולים להיחשב כמקור לתמונת ספקטרום גרועה, ואם ערך הפרמטר lowk נמוך, מוליך הגל כביכול לא סדיר, ונחשב פתרון העובי הנמוך ביותר. ניתןלהשתמש גם בגישות המבוססות על למידת מכונה, שאינה זקוקה למודלים פיזיים.
כפי שנדון בהקדמה, תקן הזהב הנוכחי לאיתור חולים בסיכון לשבר שברירי הוא DXA, שמתמודד עם כמה מגבלות: גודלו הגדול, זמינותו הירודה באזורים מסוימים, עלותו הגבוהה יחסית ויעילותו המתונה יחסית. ניתן להקל על המגבלות הראשונות באמצעות מכשירים קוליים, הידועים בניידות ובעלויות האטרקטיביות שלהם. עם זאת, היכולת לזהות ביעילות חולים בסיכון צריכה להיות לפחות שווה ערך ל-DXA. במציאות, לפעמים הוא צפוי להיות גבוה יותר מ-DXA כדי להצדיק את ההתאמה של רוב ההפניות (החלטות רפואיות, טיפול, עלויות, חדרים וכו') המקושרות לתקן הזהב. זו הסיבה שחלק מהמכשירים האולטראסונים מציעים תחליפי aBMD 10,11,42. עם זאת, חיסרון אחד של פרמטרים קליניים, כגון aBMD וגם vFAS, הוא שילוב של תכונות שונות של עצם קליפת המוח. לכן נקודות מבט משלימות מוצעות על ידי מכשירים קוליים אחרים, כולל BDAT, המציעים פרמטרים שניתן לפרש ביתר קלות על ידי הרופא והמטופל, כגון נקבוביות קליפת המוח, עובי או התפלגות גודל נקבוביות9. פרמטרים אלה משקפים תכונות גיאומטריות וחומריות: ניתן להעריך את העצם בקליפת המוח במונחים של וריאציות בלתי תלויות בכמות או באיכות. נקודת מבט זו יכולה להיות מועילה מאוד בחקירת גורמים אפשריים שונים לשבריריות העצם. לדוגמה, שברים שבריריים של הירך בתוך או מחוץ לקפסולה, כלומר אתרי שבר בצוואר הירך או בטרוכנטר, אמורים להיות בעלי מקורות רפואיים שונים43. כמו כן, ניתן לעקוב אחר תרופות שונות המכוונות להשפעות שונות על עצם קליפת המוח, כמו גם מבחינת כמות או איכות3.
שימו לב שהדיוק המתקבל עם BDAT לעובי קליפת המוח (0.1 מ"מ) טוב יותר משיטות אולטרסאונד אחרות, בדרך כלל מעל 0.25 מ"מ44. הבדל זה נובע בחלקו מהעובדה שבעיית ה-BDAT ההפוכה לוקחת בחשבון שינויים גיאומטריים וחומריים משולבים. כמה גישות אחרות, כמו בהד דופק, שוקלות תכונות ייחודיות של חומר עצם עבור כל החולים10,44. ערך דיוק זה, כ-0.1 מ"מ (CV כ-3%), הוא אכן חיוני להבחנה עדינה בין מטופלים מכיוון שטווח העובי הוא פחות מ-2 מ"מ. הדיוק בנקבוביות קליפת המוח (1%, CV כ-14%) עדיין לא טוב כמו בעובי. עם זאת, שיפורים משמעותיים כבר נצפו ביחס למחקר השחזור הקודם22. ניתן לצפות כי שיפורים דומים יושגו בעתיד הקרוב הודות לשיפורים עתידיים ב-HMI, במיוחד במונחים של קצב פריימים קרוב יותר לזמן אמת.
ניתן להשתמש ב-BDAT בקנה מידה גדול לסינון אוכלוסייה באזורים שבהם DXA אינו זמין באופן נרחב. יתר על כן, התוצאות הקליניות האחרונות הראו את הפוטנציאל שה-BDAT יכול להיות יעיל אפילו יותר מ-DXA. עם זאת, יש לאשר תוצאות אלה על ידי הכללת חולים נוספים. האתגר הבא צריך להיות מחקרים רב-מרכזיים ו/או מחקרי אורך11,12. עם זאת, מכשיר ה-BDAT הוא עדיין אב טיפוס שניתן לפגוע בו לשיתוף פעולה מדעי מכיוון שהוא כבר נעשה בגרמניה23 ובבריטניה24. יש צורך במאמצים לקראת תיעוש הדור הבא של מכשירי BDAT, שבוודאי יהיה מהיר ונייד יותר.
המחברים אינם מצהירים על ניגודי עניינים.
מחקר זה מומן על ידי הפרויקטים הצ'יליאניים ANID / Fondecyt / Regular 1201311 ו-1241091. המחברים רוצים להודות לארגון BleuSolid על תמיכתו במהלך השיפורים האחרונים ב-HMI, ולג'רופוליס על חדר הצילום.
Name | Company | Catalog Number | Comments |
Computer | Notebook HP | mod Zbook (16 Go RAM, Inrel Core i7) | to receive the sampled signals and applying the signal processing steps. Results are displayed in quasi real time (up to 4 per second) through a dedicated Human Machine Interface (HMI, BleuSolid, Pomponne, France) allowing the measurement guidance; |
Electric insulation transformer | REOMED, Solingen, Germany | IEC / EN 60601-1 | to protect the device, the patient and other devices from any electric hazard |
Electronic module | Althaïs, Tours, France | in-house | send excitation signals (half a period of negative voltage of 170 V) and discretize received signals (1024 time points per receiver at 20 MHz sampling frequency and 12 bit dynamic) before sending them to the computer. Delays and amplification can be adapted using linear laws in order to optimize data sampling within the accessible ranges. |
Human Machine Interface | BleuSolid, Pomponne, France | N/A | HMI |
Pedal switch | Scythe, Germany | USB Foot Switch 2 | to start an acquisition series |
Ruler | Westcott, USA | 10417 | to locate the measurement site |
Ultrasonic probe radius | Vermon, Tours, France | in-house | 1 MHz central frequency, 24 receivers with 0.8 mm pitch and two blocks of 5 transmitters with 1 mm pitch. |
Ultrasonic probe tibia | Vermon, Tours, France | in-house | 0.5 MHz central frequency, 24 receivers with 1.2 mm pitch and two block of 5 transmitters with 1.5 mm pitch |
Ultrasonic probes | designed according to the bidirectional geometry: a single receiver array surrounded by two transmitter arrays. The three arrays are aligned, mechanically and electrically isolated in order to minimize coupling signals. The probes are adapted to two different sites, one third distal radius and mid tibia. |
Request permission to reuse the text or figures of this JoVE article
Request PermissionThis article has been published
Video Coming Soon
Copyright © 2025 MyJoVE Corporation. All rights reserved