מוצגת שיטה לרישום משותף של תמונה רב-מודאלית שאינה דורשת סמנים פידוקיאליים לסריקות ייצור. הגישה מבוססת הפנטום יוצרת טרנספורמציה דיפרנציאלית בין מערכות הקואורדינטות של שני אופני הדמיה.
שארית של היתוך ואימות הטרנספורמציה הדיפרנציאלית
בעת חישוב הטרנספורמציה הדיפרנציאלית, התוכנה מציגה שארית של היתוך במילימטרים, המייצגת את ממוצע השורש שגיאהריבועית 19 של הטרנספורמציה. אם שארית זו חורגת מסדר הגודל של גודל הווקסל, מומלץ לבדוק את מערכי הנתונים עבור בעיות כלליות. עם זאת, מכיוון שלכל התמונות יש עיוותים קלים, השארית אינה יכולה להיות קטנה באופן שרירותי; זה רק משקף את ההתאמה של הסמנים המשמשים. לדוגמה, רישום משותף עם שלושה סמנים עשוי לגרום לשאריות קטנות יותר באותם מערכי נתונים מאשר המרה עם ארבעה סמנים מבוזרים היטב. זה קורה כי הסמנים עצמם עשויים להיות overfitted כאשר פחות fiducials מועסקים. הדיוק בכל מערך הנתונים משתפר עם מספר גדול יותר של סמנים.
הדיוק הכמותי של השיטה תלוי בזוג המכשירים הספציפיים המשמשים. ניתן לאמת את הטרנספורמציה הדיפרנציאלית המחושבת בין מערכות הקואורדינטות של שני מכשירים באמצעות השלבים הבאים: היצמדות לשלב 4 של הפרוטוקול, אך שימוש בפנטום עם סמנים פידוקיאליים כ"מדגם" שוב. הצבת הפנטום בכל מיקום, וידוא שהוא שונה מזה המשמש להערכת הטרנספורמציה הדיפרנציאלית. ניתן גם להשתמש בפנטום אחר המתאים לשיטות המתאימות, אם קיים כזה. לאחר מכן, יישום הטרנספורמציה הדיפרנציאלית שנקבעה קודם לכן (שלב 4.2.5) כדי ליישר את שני האופנים. לאחר מכן, הצבת סמנים על התמונות משני האופנים לפי שלב 3.2 של הפרוטוקול. כדי לחשב את שאריות ההיתוך עבור סמנים אלה, לחץ על תפריט היתוך > כיסוי רישום לשכבה תחתונה > המציג ציון שיורי.
השגיאה השיורית מתארת את המיקום השגוי הממוצע של האות וצריכה להיות בסדר גודל הווקסל. ספי קבלה קונקרטיים תלויים ביישום ועשויים להיות תלויים במספר גורמים, כגון הקשיחות והדיוק של מערכות ההדמיה, אך יכולים להיות מושפעים גם מתוצרי שחזור תמונה.
פתרון בעיות של עקביות עצמית
לעתים קרובות, קשיים עם עקביות עצמית נובעים מיקום אמין. שגיאה נפוצה היא הצבת המוביל במצב הפוך לרוחב. באופן אידיאלי, זה צריך להיות מוכנס מכנית לתוך מכשיר הדמיה רק בכיוון אחד. אם הדבר אינו אפשרי, יש להוסיף סימונים מובנים עבור המשתמש. נושא שכיח נוסף הוא האפשרות של תנועה בציר האורך, מה שהופך את המיקום הצירי לא אמין. מומלץ להשתמש בספייסר שניתן לחבר בקצה אחד כדי לאבטח את מיטת העכבר במקומה. ספייסרים מותאמים אישית יכולים, למשל, להיווצר במהירות ובקלות על ידי הדפסה תלת ממדית שלהם. עם זאת, התקנים מסוימים אינם יכולים לספק עקביות עצמית עם שדות ראייה משתנים. במקרים כאלה, מומלץ ליצור קשר עם הספק, אשר אמור לאשר את אי התאימות ואולי לטפל בה בעדכון עתידי. אחרת, השיטה נשארת אמינה אם נשמר שדה ראייה זהה עבור כל הסריקות, כולל כיול והדמיית ייצור.
עבור סריקות ייצור מסוימות עם מיקום חריג, טרנספורמציה למיקום המכויל אפשרית, אם ניתן להבחין במבנה נשא מספיק. עבור הדמיה in vivo , בעל החיים המרומם חייב להישאר במנשא אחד, ובניית מנשא יחיד שמתאים היטב בשני המכשירים אינה תמיד ברת השגה. לעתים קרובות, מיטת עכבר עבור מודל מבוסס נותב משמש, ולאחר מכן המיקום מאולתר לתוך מכשיר CT. לדוגמה, באיור 5A, מיטת עכבר MPI הונחה על מיטת עכבר CT עקב אילוצים מכניים. מרחב תמרון אקסיאלי ואפשרות גלגול הופכים את המיקום הזה ללא אמין. במקרים כאלה, מומלץ לתכנן מתאם המחליף את מיטת העכבר התחתונה ומאפשר התאמה משתלבת. הוא עשוי, למשל, להשתמש בטרוניונים המחוברים לחלק התחתון ובחורים נוספים בתחתית מיטת העכבר העליונה.
עם זאת, תיקון רטרוספקטיבי לתמונות קיימות אפשרי, שכן מיטת העכבר ניתנת לזיהוי בתמונת CT. הפרוטוקול מחייב סריקות כיול, ולאחריהן חישוב טרנספורמציה דיפרנציאלית של שכבת העל לשכבה התחתונה. ההליך דומה אך חייב גם למפות כל סריקת CT ייצור בודדת לסריקת הכיול, תוך שימוש במבני מיטת העכבר כפידוקיאלים.

איור 5: מיקום פתרון בעיות. (A) מיטת עכבר MPI מונחת על גבי מיטת עכבר CT. לפיכך, המיקום ב- CT לא ניתן לשחזר באופן אמין. עקביות עצמית יכולה להיות מושגת על ידי איחוי כל תמונת CT לתמונת CT הייחוס המשמשת להערכת השינוי הדיפרנציאלי. (ב-ד) מפושט לדו-ממדי. (B) כל תמונת CT הפקה נטענת כשכבת על ונרשמת בתמונת CT הייחוס (שכבה תחתונה) באמצעות מבנים של מיטת העכבר הנראים ב- CT. תמונת CT הייצור המתוקנת עקבית כעת עם CT הייחוס וניתן להשתמש בה עם הטרנספורמציה הדיפרנציאלית T. (C) שכבת MPI רשומה לתמונת CT הייחוס באמצעות הסמנים הפידוקיאליים של פנטום. (D) התמונות הרב-מודאליות מורכבות. לשם כך, כל תמונת CT ממופה למיקום הייחוס עם הטרנספורמציה הדיפרנציאלית האינדיבידואלית שלה. לאחר מכן, שכבת ה- MPI נרשמת גם למיקום הייחוס באמצעות טרנספורמציה דיפרנציאלית, שתקפה לכל התמונות של המכשיר. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.
כדי למפות את סריקות CT הייצור לסריקת הכיול, עיין בסעיף 3 של הפרוטוקול, המשלב את השינויים הבאים. לשם הבהרה, התיאור ממשיך להשתמש בדוגמה של שכבת CT ושכבת על של MPI: בשלב 3.1, טען את סריקת כיול ה- CT (תמונה A) כשכבה התחתונה ואת סריקת ה- CT לתיקון כשכבת-על. השתמש במבנים של מיטת העכבר MPI כסמנים לשלב 3.2 או כהפניות חזותיות לשלב 3.3. עקוף את שלב 3.4, אך שמירת שכבת העל מייצגת את אמצעי האחסון CT המתוקן (Menu File > Overlay > Saving Overlay as). בתיבת הדו-שיח הבאה, הזן שם חדש ולחץ על שמור. סגור את שכבת העל על-ידי ניווט אל File > Menu File > Closing Overlay. טען את סריקת ה- CT הבאה הדורשת תיקון כשכבת העל וחדש את ההליך משלב 3.2 של הפרוטוקול. הרעיון העומד בבסיס שלב זה מומחש באיור 5B.
מיטת העכבר מיושרת כעת באופן כמעט זהה לסריקת הכיול בכל נפחי ה- CT שנשמרו לאחרונה. כחלק מההליך הסטנדרטי, סריקת הכיול נרשמת לתמונות MPI באמצעות טרנספורמציה דיפרנציאלית T (איור 5C). כדי למזג לאחר מכן את תמונת ה- CT עם MPI, השתמש תמיד בנפח ה- CT המתוקן (איור 5D).
פתרון בעיות של תמונות הפוכות ושינוי קנה מידה
שיטת הרישום המוצגת כאן מניחה איכות תמונה מדויקת למדי ומתאימה רק את הסיבוב והתרגום. הוא אינו מתקן תמונות הפוכות או שינוי גודל שגוי. עם זאת, ניתן לטפל בשתי בעיות אלה באופן ידני לפני חישוב השינוי הדיפרנציאלי.
חוסר עקביות בין תבניות נתונים של יצרנים שונים עלול לגרום לערכות נתונים מסוימות, במיוחד אלה בתבנית DICOM, להיות מוצגות כהיפוך שיקוף בתוכנה. מכיוון שמיטות פנטום ועכבר הן לעתים קרובות סימטריות, ייתכן שבעיה זו לא תהיה ברורה באופן מיידי. זיהוי תמונות הפוכות קל יותר כאשר הסריקה מכילה אותיות ניתנות לזיהוי במודאליות המתאימה, כגון האותיות המוגבהות בכיוון הנכון שנראה בפנטום באיור 3H. בדוגמה המודגמת באיור 6, נתוני CT נטענים כשכבה התחתונה, ונתוני MPI נטענים כשכבת-על. זוהי סריקה in vivo של עכבר הממוקם במיטת עכבר MPI עם סמנים fiducial מחוברים. מיטת העכבר MPI ממוקמת מעל מיטת עכבר μCT (איור 6A). על-ידי היצמדות לפרוטוקול וסימון הפידוקיאלים הן בשכבה התחתונה והן בשכבת-העל בכיוון עקבי של סיבוב, מתקבלת תוצאה בלתי הולמת בעליל (איור 6B). עם זאת, בבדיקה מדוקדקת יותר ניתן לזהות את הבעיה. הפידוקיאלים יוצרים משולש א-סימטרי. בהתבוננות בצלעות המשולש במבט צירי (איור 6C, D) מהקצר ביותר לאמצעי לארוך ביותר, ניתן לראות סיבוב בכיוון השעון בנתוני ה-CT, בעוד סיבוב נגד כיוון השעון ניכר בנתוני ה-MPI. זה מדגים שאחת התמונות הפוכה לרוחב. במקרה זה, אנו מניחים שנתוני ה-CT מדויקים. כדי לתקן את שכבת העל של MPI, התמונה מתהפכת: לשם כך, העבר את השכבה שנבחרה לכיסוי ולחץ על Menu Edit > Flip > Flip X. השינוי הדיפרנציאלי המחושב על ידי התוכנה מקיף את כל הסיבובים הדרושים, ולכן "Flip X" מספיק גם אם התמונה נראית הפוכה בכיוון אחר.

איור 6: פתרון בעיות טרנספורמציה. נתוני CT נטענים כשכבה תחתונה בגודל ווקסל של 0.240 מ"מ, ונתוני MPI כשכבת על בגודל ווקסל של 0.249 מ"מ. מיטת העכבר מכילה סמנים פידוקיאליים. (A) תצוגה תלת-ממדית של תמונת הכיסוי שלא תוקנה. הפידוקיאלים בשכבת ה-CT מסומנים בחצים. הפידוקיאלים בשכבת העל של MPI נראים ככדורים בטבלת הצבעים של NIH. (B) תוצאה לא תואמת של טרנספורמציה שבוצעה ללא תיקונים מתאימים. שארית של היתוך = 6.94 מ"מ. (C) מדידת המרחקים בין הפידוקיאלים ב-CT. סיבוב בכיוון השעון מהמרחק הקצר ביותר לארוך ביותר. (ד) מדידת המרחקים בין הפידוקיאלים ב-MPI. סיבוב נגד כיוון השעון מהמרחק הקצר ביותר לארוך ביותר. השוואה למדידות CT מביאה למקדם קנה מידה של 0.928774. (E) שכבת על מתוקנת לאחר היפוך ושינוי קנה מידה. (F) שינוי צורה עם תוצאות תואמות בתצוגה תלת-ממדית. (G) טרנספורמציה עם תוצאות תואמות בתצוגה צירית. שארית של היתוך = 0.528 מ"מ. אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.
ניתן לתקן ידנית גם ערכות נתונים עם גדלי ווקסל שגויים. מכיוון שממדי הפנטום צריכים להיות ידועים, ניתן לאמת זאת בתמונה. השיטה הפשוטה ביותר היא שימוש בקצה באורך ידוע. לחץ על [Ctrl + לחצן העכבר הימני] בקצה אחד של קצה, ותוך כדי לחיצה ממושכת על הלחצן, הזז את מצביע העכבר לקצה השני של הקצה ושחרר את הלחצן. בתיבת הדו-שיח הבאה, התוכנה מציגה את אורך המרחק הנמדד בתמונה. בדוגמה המודגמת באיור 6 ניכר כי הגדלים אינם תואמים כאשר משווים את המרחקים בין הפידוקיאלים בשני האופנים (איור 6C,D). שוב, ההנחה היא שנתוני ה-CT מדויקים. כדי לשנות את קנה המידה, מחושב גורם קנה מידה (SF). מכיוון שהיחס בין האורכים (CT/MPI) אינו זהה במדויק עבור כל צד של המשולש, המנה הממוצעת מחושבת: SF = ((l1CT/l1MPI) + (l2CT/l2MPI) + (l2CT/l2MPI)) / 3.
לאחר מכן, התאם את גודל הווקסל של שכבת העל על ידי הכפלת כל ממד ב- SF. כדי להשיג זאת, החליפו את השכבה שנבחרה לשכבת-על ופתחו את Menu Edit > Change Voxel Sizes. חשב כל ממד, הזן את הערך ולאחר מכן לחץ על אישור. התוצאה של שני התיקונים מוצגת באיור 6E. לאחר מכן, שכבת העל נרשמת לשכבה התחתונה על פי הפרוטוקול. היישור המתקבל מוצג באיור 6F,G. בעוד שהדבר מספק פתרון מהיר לתיקון סריקה קיימת, אנו ממליצים לכייל את מכשיר ההדמיה לשימוש בייצור.
מגבלות
שיטה זו מוגבלת לרישום משותף מרחבי של נתונים נפחיים קיימים המורכבים מווקסלים בצורת קובייה. הוא אינו כולל תהליך שחזור שמחשב את הנפח מנתונים גולמיים המופקים על ידי מכשיר ההדמיה (למשל, תחזיות ב- CT). טכניקות שונות לשיפור התמונה קשורות לשלב זה, כגון שיטות איטרטיביות20,21 ויישום בינה מלאכותית21. למרות שהשיטה המתוארת ישימה באופן עקרוני לכל השיטות המייצרות תמונות תלת-ממדיות עם ווקסלים בצורת קובייה, לא ניתן להשתמש בה לאיחוי נתונים תלת-ממדיים עם נתונים דו-ממדיים, כמו נפח MRI בשילוב עם תרמוגרפיה אינפרא-אדומהדו-ממדית 22 או הדמיה פלואורסצנטית, שעשויה להיות רלוונטית ביישומי ניתוח מונחה תמונה. הרישום של נתונים תלת-ממדיים אינו מתקן עיוותים, כגון אלה המתרחשים בתמונות MRI בקצה הסליל. אמנם לא חובה, אך תוצאות אופטימליות מושגות כאשר מתקנים עיוותים במהלך תהליך השיקום. ההמרה האוטומטית גם אינה מטפלת בתמונות הפוכות או בשינוי קנה מידה שגוי. עם זאת, ניתן לפתור שתי בעיות אלה באופן ידני כמתואר בסעיף פתרון בעיות.
משמעות השיטה
השיטה המוצעת מבטלת את הצורך בסמנים פידוקיאליים בסריקות ייצור, ומציעה מספר יתרונות. זה מועיל לשיטות שעבורן נדרשת תחזוקת סמן או החלפה תכופה. לדוגמה, רוב סמני ה-MRI מבוססים על לחות אך נוטים להתייבש עם הזמן, וסמני PET רדיואקטיביים דועכים. על ידי הסרת הצורך בפידוקיאלים בסריקות ייצור, ניתן לצמצם את שדה הראייה, מה שמוביל לזמני רכישה קצרים יותר. זה מועיל בהגדרות תפוקה גבוהה כדי להפחית עלויות ולמזער את מינון הרנטגן בסריקת CT. מינון מופחת רצוי מכיוון שקרינה יכולה להשפיע על המסלולים הביולוגיים של חיות ניסוי במחקרי הדמיה אורכיים23.
יתר על כן, השיטה אינה מוגבלת לשיטות ספציפיות. הפשרה עבור רב-תכליתיות זו היא שפחות שלבים הם אוטומטיים. שיטה שפורסמה בעבר להתכת נתוני μCT ו-FMT משתמשת בסמנים מובנים במיטת עכבר בכל סריקה ויכולה לבצע זיהוי סמן אוטומטי ותיקון עיוות במהלך שחזור24. שיטות אחרות מבטלות את הצורך בסמנים על ידי שימוש בדמיון בין תמונות. בעוד שגישה זו מניבה תוצאות טובות ויכולה גם לתקן עיוותים25, היא ישימה רק אם שתי השיטות מספקות תמונות דומות מספיק. זה בדרך כלל לא המקרה בשילוב של מודאליות מפורטת אנטומית ומודל מבוסס נותב. עם זאת, שילובים אלה נחוצים להערכת הפרמקוקינטיקה של סוכנים ממוקדים26, שיש להם יישומים בתחומים כגון ננותרפיה אנטי סרטנית 27,28.
מכיוון שבקרת האיכות פחות קפדנית ביישומים פרה-קליניים בהשוואה ליישומים קליניים, חוסר התאמה של מכשירי הדמיה משולבים הוא בעיה מוכרת29. נתונים המושפעים מחוסר התאמה זה יכולים להשתפר בדיעבד על ידי סריקת פנטום וקביעת השינוי הדיפרנציאלי, מה שעשוי להפחית עלויות ולמזער את הפגיעה בבעלי חיים. בנוסף לשיטה המודגמת המשתמשת בסמנים פידוקיאליים לחישוב טרנספורמציה דיפרנציאלית, אשר מוחלת לאחר מכן על סריקות ייצור, אפשרויות נוספות למיזוג תמונה מתוארות ומשמשות. סקירה כללית, הכוללת הפניות לתוכנות זמינות שונות, ניתן למצוא ב- Birkfellner et al.30.
לסיכום, השיטה המוצגת מציעה פתרון יעיל לרישום משותף של תמונות רב-מודאליות. הפרוטוקול ניתן להתאמה בקלות לשיטות הדמיה שונות, וטכניקות פתרון הבעיות המסופקות משפרות את עמידות השיטה מפני בעיות אופייניות.