Abstract
التوافقي الحركة التصوير بالموجات فوق الصوتية المركزة ل(HMIFU) هو الاسلوب الذي يمكن أن تؤدي ومراقبة عالية الكثافة (HIFU) الاجتثاث ركزت الموجات فوق الصوتية. يتم إنشاء حركة متذبذبة في محور 93 عنصر و 4.5 ميغاهيرتز تردد مركز HIFU محول عن طريق تطبيق 25 هرتز، التضمين السعة إشارة باستخدام مولد وظيفة. يتم وضع 64 عنصر و 2.5 ميغاهيرتز محول التصوير مع 68kPa ذروة الضغط confocally في وسط محول HIFU للحصول على الترددات الراديوية (RF) بيانات القناة. في هذا البروتوكول، وصفت رصد في الوقت الحقيقي الاجتثاث الحرارية باستخدام HIFU مع القوة الصوتية من 7 W على كبد الكلاب في المختبر. يتم تطبيق العلاج HIFU على النسيج خلال 2 دقيقة وتصوير المنطقة ذاب في الوقت الحقيقي باستخدام التصوير مباعدة أو موجة الطائرة تصل إلى 1000 لقطة / ثانية. يتم ضرب مصفوفة بيانات القناة RF بواسطة مصفوفة متفرق لإعادة الإعمار الصورة. مجال أعيد النظر هو 90 درجة لمدة متباينة واهاء و 20 مم للتصوير موجة الطائرة ويتم أخذ عينات البيانات في 80 ميغاهرتز. يتم تنفيذ إعادة الإعمار على وحدة المعالجة الرسومية (GPU) من أجل الصورة في الوقت الحقيقي بمعدل 4.5 إطار العرض. ويستخدم 1-D تطبيع عبر الارتباط للبيانات RF أعيد بناؤها لتقدير التشريد المحوري في المنطقة المحورية. حجم النزوح الذروة إلى الذروة في العمق البؤري يتناقص خلال الاجتثاث الحرارية التي تدل على تشنج الأنسجة نتيجة لتشكيل الآفة. وكانت نسبة النزوح إشارة إلى الضوضاء (SNR د) في منطقة محورية للموجة مستوى أعلى 1.4 مرة من لتباين موجة تبين أن التصوير موجة الطائرة تظهر لإنتاج النزوح أفضل خرائط الجودة لHMIFU من تباين التصوير الموجة.
Protocol
وقد وافق هذا البروتوكول من قبل لجنة رعاية واستخدام الحيوان المؤسسي من جامعة كولومبيا. تم تنفيذ جميع الحصول على البيانات ومعالجتها باستخدام بيئة ماتلاب.
1. انشاء التجريبية
- ديغا في الجسم الحي الكلاب عينة الكبد السابقين خلال 90 دقيقة. وضع عينة الكبد في خزان مملوء نزع الغاز الفوسفات الحل (الشكل 1). إصلاح العينة الكبد على امتصاص الصوتية مع الإبر في الأطراف من الكبد.
- إدراج 64 عنصرا، 0.32 ملم الملعب، 2.5 ميغاهيرتز تردد مركز مجموعة مراحل (التصوير) من خلال ثقب دائري يقع في وسط ل93-عنصر صفيف نصف كروية HIFU محول (العلاجي) في مركز تردد 4.5 ميغاهيرتز، 70 مم العمق البؤري و 1.7 مم × 0.4 مم حجم البؤري 11. محاذاة كل من محولات شارك في محوريا وإصلاح محول التصوير في محول العلاجي مع مسامير التكيف.
- تغطية محول HIFU مع مركباتolume تسيطر عليها غشاء من مادة البولي يوريثين مليئة المياه المتدفقة نزع الغاز لتبريده. تركيب محول التجمع على 3-D مناور الكمبيوتر التي تسيطر عليها.
- قم بتوصيل محول HIFU إلى مولد وظيفة إرسال 25 هرتز السعة التضمين الموجي شرط مع 500 فولت كحد أقصى السعة. قم بتوصيل محول التصوير لنظام الموجات فوق الصوتية للبرمجة بشكل كامل باستخدام برنامج ماتلاب.
ملاحظة: هناك برامج المرتبطة بنظام الموجات فوق الصوتية واستخدام بيئة ماتلاب أن تكون مثبتة على جهاز كمبيوتر متصل إلى النظام. يجب أن توضع A 50 ديسيبل مكبر للصوت RF وشبكة المطابقة بين محول HIFU ومولد وظيفة لتضخيم التوالي القوة وتطابق مقاومة. - إنشاء شبكة القطبية، وذلك باستخدام ماتلاب، بدءا من 50 مم من سطح مجموعة و 40 مم العميقة في الاتجاه شعاعي مع خطوة المكانية من 9.625 ميكرون و90 درجة في الاتجاه السمتي مع 128 خطوط والتي الأصل هو FOسلطات الجمارك من موجة متباينة. تحديد مصدر متباينة موجة 10.24 مم (نصف حجم الفتحة) وراء السطح من مجموعة، وتركزت في الاتجاه الجانبي.
- إنشاء شبكة الديكارتية، وذلك باستخدام ماتلاب، بدءا من 50 مم من سطح مجموعة و 40 مم العميقة في الاتجاه المحوري مع خطوة المكانية من 9.625 ميكرون و 20 ملم واسعة في الاتجاه الجانبي مع 64 خطوط لموجة الطائرة. تحديد مصدر موجة الطائرة على سطح المصفوفة. لكل شبكة، حساب الوقت من المصدر إلى كل نقطة من الشبكة والعودة إلى كل عنصر من عناصر المصفوفة.
- أدخل "ReconMat_DW" لتباين التصوير موجة أو "ReconMat_PW" للتصوير موجة الطائرة في إطار الأوامر مطلب واضغط على "دخول" لإنشاء مصفوفة إعادة الإعمار المرتبطة خوارزمية تأخير ومبلغ القياسية لكل شبكة. تطبيق خوارزمية تأخير ومحصلتها إلى كل متجه أساس معيار واسترداد elemen غير الأصفارالخبر من الناتج مصفوفة 11. تخصيص عناصر غير الصفرية التي تم الحصول عليها من المصفوفة الناتجة إلى مصفوفة متفرق في الموقع المقابل. حفظ مصفوفة إعادة الإعمار على القرص الصلب لجهاز الكمبيوتر.
ملاحظة: المتباينة وأساليب موجة الطائرة استخدام اثنين من مصفوفات إعادة الإعمار متميزة.- يلقي مصفوفة إعادة الإعمار إلى مصفوفة GPU. أدخل "SetUpP4_2Flash_4B_streaming_DW" لتباين التصوير موجة أو "SetUpP4_2Flash_4B_streaming_PW" للتصوير موجة الطائرة في إطار الأوامر ماتلاب ثم اضغط على "دخول" لإنشاء ملف الإعداد للحصول على البيانات قناة الموجات فوق الصوتية باستخدام البرنامج النصي المرتبطة مع مجموعة مراحل، وقدمت من قبل الشركة المصنعة من نظام الموجات فوق الصوتية. اسم ملف الإعداد "P4-2Flash_DivergingWave.mat" لتباين التصوير موجة و "P4-2Flash_PlaneWave.mat" للتصوير موجة الطائرة.
ملاحظة: حزمة البرمجيات التجارية يجب أن تكون مثبتة على جهاز الكمبيوتر رس يلقي إعادة الإعمار المصفوفة المتناثرة لمصفوفة GPU.
- يلقي مصفوفة إعادة الإعمار إلى مصفوفة GPU. أدخل "SetUpP4_2Flash_4B_streaming_DW" لتباين التصوير موجة أو "SetUpP4_2Flash_4B_streaming_PW" للتصوير موجة الطائرة في إطار الأوامر ماتلاب ثم اضغط على "دخول" لإنشاء ملف الإعداد للحصول على البيانات قناة الموجات فوق الصوتية باستخدام البرنامج النصي المرتبطة مع مجموعة مراحل، وقدمت من قبل الشركة المصنعة من نظام الموجات فوق الصوتية. اسم ملف الإعداد "P4-2Flash_DivergingWave.mat" لتباين التصوير موجة و "P4-2Flash_PlaneWave.mat" للتصوير موجة الطائرة.
- مزامنة نظام الموجات فوق الصوتية مع مولد وظيفة باستخدام الزناد الخارجي بحيث الحصول على البيانات ارتفاع معدل الإطار الموجات فوق الصوتية للكبد يبدأ في نفس الوقت كما HIFU.
- فتح ماتلاب. تشغيل البرنامج النصي الإعداد "SetUpP4_2Flash_4B.m" المقدمة من قبل الشركة المصنعة نظام الموجات فوق الصوتية لاستخدام التصوير B-الوضع. اسم ملف الإعداد التي تم إنشاؤها: "P4-2Flash_4B_Bmode.mat". استخدام "VSX" القيادة وعندما "اسم ملف .mat على العملية:" إن المطالبة، أدخل اسم ملف الإعداد "P4-2Flash_4B_Bmode.mat". نقل كل محولات الطاقة واستخدام شاشة عرض B-الوضع الذي ظهر على شاشة الكمبيوتر لوضعها في المنطقة المستهدفة من الكبد إلى يجتذ. استهداف المنطقة حوالي 1 سم تحت سطح الكبد لتجنب ارتفاع توهين الموجات فوق الصوتية بسبب الامتصاص. حفظ التقليدية B-وضع صورة للكبد على الكمبيوتر.
ملاحظة:نحن هنا يقوم ablations HIFU في 11 مواقع مختلفة في عينتين الكبد عن طريق تحريك محولات مع مناور 3-D لكل الاجتثاث.
2. الموجات فوق الصوتية الحصول على البيانات
- فتح ماتلاب. استخدام "VSX" القيادة وعندما "اسم ملف .mat على العملية:" إن المطالبة، أدخل اسم ملف الإعداد "P4-2Flash_DivergingWave.mat" لتباين التصوير موجة أو "P4-2Flash_PlaneWave.mat" لموجة الطائرة التصوير. بدء HIFU وتطبيقه خلال 2min إلى المنطقة المستهدفة.
- الحصول على بيانات القناة RF في 1000 لقطة في الثانية خلال 2 دقيقة باستخدام موجات متباينة. بدلا من ذلك، الحصول على بيانات القناة RF في 1000 لقطة في الثانية خلال 2 دقيقة باستخدام موجات الطائرة.
- نقل البيانات إلى الكمبيوتر المضيف كل 200 لقطة عن طريق كابل صريح PCI. بدلا من ذلك، في الوقت الحقيقي الجري، الحصول على بيانات القناة RF عند 167 لقطة في الثانية خلال 2 دقيقة باستخدام موجات الطائرة وtransfإيه البيانات إلى الكمبيوتر المضيف كل 2 الإطارات.
ملاحظة: وسائل التصوير مع مجموعة من 200 لقطة يوفر القرار الزماني عالية ضمن مجموعة ولكن خلق كل الفجوات بين كل مجموعة ومناسبة للمعالجة خارج الخط. طريقة التصوير في 167 إطارا في الثانية لديها من القرار الزمني أقل ولكن لا يخلق أي ثغرات في الوقت الاجتثاث بأكمله و المناسب في الوقت الحقيقي الجري. - يلقي مصفوفة بيانات القناة RF إلى واحدة مصفوفة الدقة GPU مع مطلب. ضرب مصفوفة بيانات القناة RF بالمصفوفة إعادة الإعمار للحصول على البيانات RF أعيد بناؤها 11.
3. التشرد التصوير
- إنشاء 6 أجل عشر بتروورث مرشح تمرير منخفض في 4 ميغاهيرتز تردد قطع باستخدام DSP نظام أدوات للمطلب. تطبيق منخفض هذا مرشح تمرير للبيانات RF أعيد بناؤها لتصفية عنصر ميغاهيرتز HIFU 4.5.
- تقدير انزياح بين إطارات متتالية باستخدام 1-D تطبيع عبر الارتباطمع طول 3.1mm نافذة و 90٪ التداخل.
- إنشاء 6 أجل عشر بتروورث مرشح تمرير منخفض عند 100 هرتز تردد قطع باستخدام DSP نظام أدوات للمطلب. تطبيق منخفض هذا مرشح تمرير للبيانات النزوح الزمنية باستخدام ماتلاب لاسترداد المكون تردد 50 هرتز، متذبذبة.
- تحديد المنطقة ذات الاهتمام (ROI) كمنطقة محورية في -6 ديسيبل (1.7 X 0.4 مم في الماء) وتقع 70 ملم بعيدا عن سطح محول. استخراج البيانات التشريد في هذا ROI. تقدير تشريد نسبة الإشارة إلى الضوضاء (SNR د) في منطقة الوصل بعد 2 دقيقة من الاجتثاث بأنها النسبة بين يعني التشرد والانحراف المعياري للنزوح في العائد على الاستثمار.
- استخراج 50 هرتز إشارة النزوح الزمني في التركيز من البيانات تشريد المصفوفة. تحويل الإشارات النزوح الزمنية في التركيز إلى صوت مسموع باستخدام ماتلاب.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
في الوقت الحقيقي يتدفقون من HMI النزوح خلال HIFU الاجتثاث يمكن الحصول عليها باستخدام مباعدة وطائرة التصوير موجة الشكل 2 هو القبض على الشاشة فيديو يظهر الوقت الحقيقي العرض من قوة الإشعاع الصوتية النزوح الناجم عن استخدام التصوير موجة الطائرة في في المختبر الكلاب كبد خلال HIFU الاجتثاث . هي تدفق من النزوح في الوقت الحقيقي على شاشة الكمبيوتر في معدل الإطار عرض 4.5 هرتز. وتظهر التشريد إيجابية في نزوح الحمراء والسلبية في الزرقاء. تم تسليم الآفات بنجاح باستخدام HIFU الاجتثاث. ويبين الشكل 3 الآفة التي تم الحصول عليها في الكبد بعد استئصال الموافق الشكل 2.
بانخفاض قدره HMI ذروة النزوح سلبي السعة خلال HIFU الاجتثاث ولا يمكن تصوير مع كل مباعدة وطائرة التصوير موجة الشكل 4 يظهر نزوح سلبي HMI الذروة في مرحلة مختلفة من الاجتثاث مع المتباينة وموجة الطائرةالتصوير. وقد أظهرت نزوح السلبية الذروة دون ومع تراكب على B-وضع لرؤية أكثر وضوحا في نمط التشرد ورؤية المنطقة المستهدفة في الكبد على التوالي كلا. تأسست في 50 هرتز HMI الصوت الإزاحة المقابلة لالاجتثاث رصدها مع موجة الطائرة (الشكل 4C) إلى الفيديو. انخفاض HMI سعة الإزاحة بسبب الاجتثاث ويمكن سماع الذي يوفر أداة رصد إضافية. الشكل 4 يشير ايضا الى ان حجم المنطقة ولع زيادات HIFU خلال الاجتثاث. 5A الشكل و5B يظهر نزوح HMI في منطقة الوصل خلال الاجتثاث عن موجة متباينة والطائرة على التوالي. انخفاض HMI حجم النزوح واضحة للعيان سواء بالنسبة المتباينة والتصوير موجة الطائرة الشكل 6 يوضح انخفاض النزوح الذروة إلى الذروة لجميع المواقع المستهدفة في الكبد سواء بالنسبة للتباين (الشكل 6A) وطائرة (الشكل 6B) التصوير الموجة. انخفاض النزوح الذروة إلى الذروة لموجة الطائرة لا يختلف كثيرا عن واحد تم الحصول عليها متباينة الموجة.
تم العثور التصوير موجة الطائرة لديها SNR العالي د في التركيز من تباين التصوير موجة الشكل 7 يبين SNR د في العائد على الاستثمار لجميع المواقف الآفة في الكبد للتباين (الشكل 7A) وطائرة (الشكل 7B) التصوير موجة . وكان متوسط SNR د لمستوى أعلى 1.7 مرة من لتباين التصوير الموجة.
الشكل 1. التجريبية مجموعة المتابعة. (A) تمثيل نظام HMIFU. (B) صورة انشاء التجريبية.OM / ملفات / ftp_upload / 53050 / 53050fig1large.jpg "الهدف =" _ فارغة "> الرجاء انقر هنا لعرض نسخة أكبر من هذا الرقم.
الشكل 2. في الوقت الحقيقي نزوح HMI. التقاط شاشة الكمبيوتر تظهر الوقت الحقيقي يتدفقون من التشريد HMI مع التصوير موجة الطائرة خلال HIFU الاجتثاث من الكبد الكلاب في 4.5Hz معدل الإطار العرض. معارض اليسرى لوحة الجانب والتشريد HMI التي تمت تصفيتها ويظهر الحق على لوحة الجانب والتشريد HMI تصفيتها مضافين على ما قبل الاجتثاث B-وضع الكبد. الرجاء انقر هنا لعرض نسخة أكبر من هذا الرقم.
الشكل 3. الآفات التي تسببها مرحباFU. صورة المتوسطة المقطع العرضي للآفة بعد العلاج HIFU. الرجاء النقر هنا لمشاهدة نسخة أكبر من هذا الرقم.
الرقم 4. متباينة والتصوير موجة الطائرة من النزوح. الذروة سلبي HMI التصوير النزوح خلال HIFU الاجتثاث من الكبد الكلاب باستخدام متباينة موجة من دون غطاء B-وضع (A)، مع B-وضع تراكب (B)، وذلك باستخدام التصوير موجة الطائرة مع عدم وجود B-وضع تراكب (C) ومع B-وضع تراكب (D). 50 هرتز HMI الصوت الإزاحة المقابلة لالاجتثاث رصدها مع موجة الطائرة (الشكل 4C) تأسست على الفيديو. Pتأجير انقر هنا لعرض نسخة أكبر من هذا الرقم.
الرقم 5. النزوح البؤري HMI. HMI التهجير في المنطقة خلال التنسيق HIFU الاجتثاث باستخدام متباينة (A) وطائرة (B) التصوير الموجة. الرجاء انقر هنا لعرض نسخة أكبر من هذا الرقم.
الرقم 6. انخفاض النزوح الذروة إلى الذروة. انخفاض النزوح الذروة إلى الذروة في المنطقة خلال التنسيق HIFU الاجتثاث باستخدام متباينة (A) وطائرة (B) التصوير الموجة. الرجاءانقر هنا لمشاهدة نسخة أكبر من هذا الرقم.
الرقم 7. المهجرين نسبة الإشارة إلى الضوضاء. نسبة النزوح إشارة إلى الضوضاء في تركيز متباينة (A) وطائرة (B) التصوير موجة لموقف الاجتثاث مختلفة. الرجاء انقر هنا لعرض نسخة أكبر من هذا الرقم .
Subscription Required. Please recommend JoVE to your librarian.
Discussion
في الوقت الحقيقي رصد الآفات HIFU مهم لضمان السليم والفعال تسليم الآفة. كما الأشكال الآفة، والأنسجة صمود والسعة تلقاء تحت الإثارة النقصان. تطبيق HIFU في المنطقة من نتائج الأنسجة في قوة الإشعاع الصوتية أن يدفع النزوح الأنسجة. التغير النسبي في النزوح هو بديل من التغير النسبي في تصلب الأنسجة. هذه التقنية توفر ميزة مراقبة HIFU الآفة دون توقف العلاج على النقيض من الأساليب الأخرى التي تتخذ من الموجات فوق الصوتية. وكان هذا القرار الزمني للرصد في الوقت الحقيقي في هذه الدراسة (4.5 هرتز) أعلى من تلك التي حصلت في MR الموجهة HIFU الاجتثاث (1 هرتز).
معالجة سريعة للبيانات RF الموجات فوق الصوتية هو خطوة حاسمة في الوقت الحقيقي يتدفقون من النزوح. إعادة بناء صورة هو أبطأ خطوة من المعالجة. في هذا البروتوكول، وكان محسن سرعة إعادة بناء صورة عن طريق الحصول على الإطار بأكمله باليودنانوغرام عملية واحدة. وتتألف هذه العملية في ضرب بيانات القناة RF بواسطة مصفوفة. تم تخصيص العناصر غير الصفر من المصفوفة فقط لتحسين حساب الوقت و
تم تنفيذ الضرب على GPU. تم استخدام سريع 1-D تطبيع طريقة عبر ارتباط لتقدير التشريد. A التداخل نافذة 80٪ يسمح للمفاضلة جيدة بين وقت الحساب والقرار المحوري للصور التشرد.
طريقة الإرسال تكوين الشعاع يمكن أن تؤثر أيضا على جودة الصورة النزوح. تم العثور على SNR د أن تكون أقل من ذلك بكثير لمتباينة من لطائرة موجة التصوير باستخدام عينة اثنين اختبار (ت). كان حجم النزوح أيضا أقل لمتباينة من لتصوير موجة الطائرة. ويمكن تفسير ذلك من خلال حقيقة أن الاتجاه المحوري للموجة متباينة لا تتماشى مع شعاع HIFU في ROI بأكمله نظرا لطبيعة المتباينة للموجة فيعلى النقيض من موجة الطائرة. انخفاض النزوح أقل من الذروة إلى الذروة تم العثور عليها ل الآفة رقم 3 لتباين التصوير موجة يمكن أن يكون بسبب وجود سفينة في مركز الآفة لوحظ بعد أمراض الإجمالي. في أقل SNR تم العثور عليها ل الآفة رقم 4 للتصوير موجة الطائرة د يمكن أن يكون نتيجة لقربها من التركيز على سطح الكبد. كما أن لديها تجدر الإشارة إلى أن توهين بسبب نشر هندسية في الاتجاه الأفقي يحدث للموجات متباينة وليس لموجة الطائرة التي يمكن أن تؤثر على نوعية تقدير الحركة. ومع ذلك، عند استخدام نفس محول بالموجات فوق الصوتية، والتصوير موجة متباينة يقدم أكبر مجال الرؤية من التصوير موجة الطائرة وهي تهم باستمرار صورة أكبر جزء من منطقة إلى يجتذ.
في هذا البروتوكول، استخدمت مجموعة مراحل لصورة النزوح بحيث تم تصويرها فقط قطاعا عريضا من المنطقة ذاب. ويمكن استخدام A 2-D مجموعة محول لصورة رانه كامل حجم المنطقة ذاب. وقد تحقق الاجتثاث في مواقع مختلفة من الكبد عن طريق نقل محول المتعلقة الكبد. يمكن أن يؤديها شعاع القيادة مع لجنة التحقيق HIFU لاستهداف مواقع مختلفة من المنطقة المراد علاجها للسماح لمزيد من الاستهداف المناسب. إلى جانب التحسينات التقنية المذكورة آنفا، وتشمل التوجهات المستقبلية الترجمة السريرية لهذا الأسلوب.
Subscription Required. Please recommend JoVE to your librarian.
Materials
| Name | Company | Catalog Number | Comments |
| P4-2 Phased array | ATL | ||
| H-178 HIFU transducer | Sonic Concepts | ||
| 3-D positioner | Velmex Inc. | ||
| AT33522A function generator | Agilent Technologies | ||
| V-1 ultrasound system | Verasonics | ||
| 3100L RF amplifier | ENI | ||
| Matching network | Sonic Concepts | ||
| Degasing system | Sonic Concepts | ||
| Programming software | Matlab | ||
| Jacket software package | Accelereyes |
References
- Al-Bataineh, O., Jenne, J., Huber, P. Clinical and future applications of high intensity focused ultrasound in cancer. Cancer Treat Rev. 38, 346-353 (2012).
- Dewhirst, M. W., Viglianti, B. L., Lora-Michiels, M., Hanson, M., Hoopes, P. J. Basic principles of thermal dosimetry and thermal thresholds for tissue damage from hyperthermia. Int J Hyperthermia. 19, 267-294 (2003).
- Napoli, A., et al. MR-guided high-intensity focused ultrasound: current status of an emerging technology. Cardiovasc Intervent Radiol. 36, 1190-1203 (2013).
- Gudur, M. S., Kumon, R. E., Zhou, Y., Deng, C. X. High-frequency rapid B-mode ultrasound imaging for real-time monitoring of lesion formation and gas body activity during high-intensity focused ultrasound ablation. IEEE Trans Ultrason Ferroelectr Freq Control. 59, 1687-1699 (2012).
- Jensen, C. R., Cleveland, R. O., Coussios, C. C. Real-time temperature estimation and monitoring of HIFU ablation through a combined modeling and passive acoustic mapping approach. Phys Med Biol. 58, 5833-5850 (2013).
- Mariani, A., et al. Real time shear waves elastography monitoring of thermal ablation: in vivo evaluation in pig livers. J Surg Res. 188, 37-43 (2014).
- Bing, K. F., Rouze, N. C., Palmeri, M. L., Rotemberg, V. M., Nightingale, K. R. Combined ultrasonic thermal ablation with interleaved ARFI image monitoring using a single diagnostic curvilinear array: a feasibility study. Ultrason Imaging. 33, 217-232 (2011).
- Athanasiou, A., et al. Breast lesions: quantitative elastography with supersonic shear imaging--preliminary results., Radiology. 256, 297-303 (2010).
- Maleke, C., Konofagou, E. E. Harmonic motion imaging for focused ultrasound (HMIFU): a fully integrated technique for sonication and monitoring of thermal ablation in tissues. Phys Med Biol. 53, 1773-1793 (2008).
- Maleke, C., Konofagou, E. E. In vivo feasibility of real-time monitoring of focused ultrasound surgery (FUS) using harmonic motion imaging (HMI). IEEE Trans Biomed Eng. 57, 7-11 (2010).
- Hou, G. Y., et al. Sparse matrix beamforming and image reconstruction for 2-D HIFU monitoring using harmonic motion imaging for focused ultrasound (HMIFU) with in vitro validation. IEEE Trans Med Imaging. 33, 2107-2117 (2014).
