Аннотация 4
1. Введение 5
1.1. Актуальность 5
1.2. Обзор литературы и ключевых работ в предметной области 6
1.2.1. Travel-time томография 6
1.2.2. Compressive Sensing 9
1.3. Предварительный пример 13
1.4. Краткий обзор результатов работы 14
1.5. Структура работы 15
2. Постановка задачи 16
3. Основной результат 17
3.1. Исследование разреженности изображения 17
3.2. Ожидания по масштабированию задачи 19
3.3. Проектирование универсальной матрицы
измерений 22
3.4. Алгоритм реконструкции при использовании
Compressive Sensing 24
3.5. Архитектура программного решения 24
4. Моделирование 26
4.1. Результаты эксперимента на маломасштабной модели . . 27
5. Заключение 30
Список литературы
Ультразвуковая томография по качеству и разрешающей способности достигла сопоставимого с МРТ уровня и активно применяется для исследования мягких тканей [7]. Преимуществами УЗИ являются относительно низкая стоимость оборудования и обслуживания, безопасность для организма, неинвазивность техники исследования.
Для повышения качества и разрешающей способности получаемого в ходе исследования изображения требуется увеличение количества используемых датчиков и частоты дискретизации сигнала. Всё это ведёт к значительному увеличению объема передаваемых данных, что усложняет как производственный процесс, так и проведение диагностики. Кроме того особенности проведения исследования допускают наличие помех и искажений [28].
Опухоли преимущественно имеют более высокую скорость прохождения ультразвука, чем окружающие ткани. Это делает возможным реконструкцию плотностей тканей в исследуемой зоне с помощью уравнений с участием предполагаемых путей распространения сигнала и временем его прибытия на датчики кругового массива (travel time tomography)[25]. Современные системы основаны именно на таком подходе и получили применение в диагностике рака молочных желез [26][30].
Для метода томографии travel-time типична квадратичная от числа сенсоров зависимость получаемых “сырых” данных для анализа: последовательно с каждого датчика пускается ультразвуковой импульс, который принимают остальные к — 1 сенсоров. Это приводит к серьезному повышению требований к вычислительной части устройства томографа, а также к увеличению времени обработки: современные методы реконструкции работают итеративно с временной сложностью итерации O(N log N) [11].
На настоящий момент использование ультразвуковой томографии осложнено высокими требованиями к вычислительному комплексу и временем на обработку [24]. Уменьшение объема данных позволит эффективно использовать FPGA вычислители для задач определения времени прибытия сигнала на датчик и последующей реконструкции изображения. Compressive Sensing также позволит сократить время на сбор данных путем уменьшения количества производимых проекций, что частично сократит зашумление данных из-за колебаний пациента в процессе обследования.
В ходе работы были достигнуты следующие результаты:
• Исследована возможность применения техники Compressive Sen¬sing в задаче ультразвуковой томографии.
• Разработан и реализован на языке MATLAB рандомизированный алгоритм для сбора и обработки данных ультразвуковой томографии.
• Проведены эксперименты на программной симуляции.
В этой работе была исследована возможность применения техники Опознания со сжатием, или “Compressive Sensing”, в области ультразвуковой томографии “времени прибытия”: изучены свойства разреженности результатов полученных изображений, а также влияние увеличения разрешения снимка на числовую характеристику разреженности s. Пренебрежимое увеличение этой характеристики позволяет эффективно проводить обработку данных, даже при увеличении их исходного количества на много порядков.
Также была разработана и апробирована на искусственных моделях модификация алгоритма реконструкции томографических снимков с применением техники опознания со сжатием. Эксперименты показали, что возможна реконструкция изображения без существенной потери качества при использовании неполного объема данных. Кроме того, некоторые варианты построения матриц измерений продемонстрировали положительное влияние на шумоустойчивость алгоритма.
[1] Граничин ОН, Поляк БТ. Рандомизированные алгоритмы оценивания и оптимизации при почти произвольных помехах. — Наука М., 2003.
[2] Граничин Олег Николаевич, Павленко Дмитрий Валентинович. Рандомизация получения данных и 1-оптимизация (опознание со сжатием) // Автоматика и телемеханика. — 2010. — no. 11. — P. 3¬28.
[3] A simple proof of the restricted isometry property for random matrices / Richard Baraniuk, Mark Davenport, Ronald DeVore, Michael Wakin // Constructive Approximation. — 2008. — Vol. 28, no. 3. — P. 253-263.
[4] An improved automatic time-of-flight picker for medical ultrasound tomography / Cuiping Li, Lianjie Huang, Nebojsa Duric et al. // Ultrasonics. — 2009. — Vol. 49, no. 1. — P. 61-72.
[5] Berinde Radu, Indyk Piotr, Ruzic M. Practical near-optimal sparse recovery in the l1 norm // Communication, Control, and Computing, 2008 46th Annual Allerton Conference on / IEEE. — 2008. — P. 198¬205.
[6] Breast imaging using ultrasound tomography: From clinical requirements to system design / Olivier Roy, Signe Schmidt, Cuiping Li et al. // Ultrasonics Symposium (IUS), 2013 IEEE International / IEEE. -- 2013. -- P. 1174-1177.
[7] Breast imaging with 3D ultrasound computer tomography: results of a first in-vivo study in comparison to MRI images / Torsten Hopp, Lukas Sroba, Michael Zapf et al. // Breast Imaging. — Springer,
2014. — P. 72-79.
[8] Cande Emmanuel J, Wakin Michael B. An introduction to compressive sampling // Signal Processing Magazine, IEEE. — 2008. — Vol. 25, no. 2. — P. 21-30.
[9] Candes Emmanuel, Romberg Justin. Sparsity and incoherence in compressive sampling // Inverse problems. — 2007. — Vol. 23, no. 3. — P. 969.
[10] Candes Emmanuel J, Romberg Justin, Tao Terence. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information // Information Theory, IEEE Transactions on. — 2006. — Vol. 52, no. 2. — P. 489-509.
[11] Chen Chen, Huang Junzhou. Compressive sensing MRI with wavelet tree sparsity // Advances in neural information processing systems. — 2012. —P. 1115-1123.
[12] Chen Scott Shaobing, Donoho David L, Saunders Michael A. Atomic decomposition by basis pursuit // SIAM review. — 2001.— Vol. 43, no. 1. —P. 129-159.
[13] Compressed sensing MRI / Michael Lustig, David L Donoho, Juan M Santos, John M Pauly // Signal Processing Magazine, IEEE. — 2008. — Vol. 25, no. 2. — P. 72-82.
[14] Compressed sensing of ultrasound images: sampling of spatial and frequency domains / Celine Quinsac, Adrian Basarab, Jean¬Marc Girault, Denis Kouame // Signal Processing Systems (SIPS), 2010 IEEE Workshop on / IEEE. — 2010. — P. 231-236.
[15] Detection of breast cancer with ultrasound tomography: First results with the Computed Ultrasound Risk Evaluation (CURE) prototype / Nebojsa Duric, Peter Littrup, Lou Poulo et al. // Medical physics. — 2007. - Vol. 34, no. 2. - P. 773-785.
[16] Dines Kris A, Lytle R Jeffrey. Computerized geophysical tomography // Proceedings of the IEEE. — 1979. — Vol. 67, no. 7. — P. 1065-1073.
[17] Donoho David L. Compressed sensing // Information Theory, IEEE Transactions on. — 2006. — Vol. 52, no. 4. — P. 1289-1306.
[18] Donoho David L. For most large underdetermined systems of linear equations the minimal L1-norm solution is also the sparsest solution // Communications on pure and applied mathematics. -- 2006. -- Vol. 59, no. 6. -- P. 797-829.
[19] Granichin Oleg. Linear regression and filtering under nonstandard assumptions (Arbitrary noise) // Automatic Control, IEEE Transactions on. — 2004. — Vol. 49, no. 10. — P. 1830-1837.
[20] Leonid Kunyansky. Fast reconstruction algorithms for the thermoacoustic tomography in certain domains with cylindrical or spherical symmetries //Inverse Problems and Imaging.— 2012.— Vol. 6, no. 1. —P. 111-131.
[21] Lustig Michael, Donoho David, Pauly John M. Sparse MRI: The application of compressed sensing for rapid MR imaging // Magnetic resonance in medicine. — 2007. — Vol. 58, no. 6. — P. 1182-1195.
[22] Mallat Stephane G, Zhang Zhifeng. Matching pursuits with time¬frequency dictionaries // Signal Processing, IEEE Transactions on. -¬1993. —Vol. 41, no. 12. —P. 3397-3415.
[23] Natarajan Balas Kausik. Sparse approximate solutions to linear systems // SIAM journal on computing. — 1995.— Vol. 24, no. 2.— P. 227-234.
[24] Ozmen N. Ultrasound Imaging Methods for Breast Cancer Detection. — TU Delft, Delft University of Technology, 2014.
[25] Quan Youli, Huang Lianjie. Sound-speed tomography using first¬arrival transmission ultrasound for a ring array // Medical Imaging / International Society for Optics and Photonics. -- 2007. -- P. 651306¬651306.
[26] Robust ultrasound travel-time tomography using the bent ray model / Ali Hormati, Ivana Jovanovic, Olivier Roy, Martin Vetterli // Proc. SPIE. — Vol. 7629. — 2010. — P. 76290I.
[27] Schuster Arthur. An introduction to the theory of optics. — E. Arnold, 1904.
[28] Shannon Claude E. Communication in the presence of noise // Proceedings of the IRE. — 1949. — Vol. 37, no. 1. — P. 10-21.
[29] Treeby Bradley E, Cox Benjamin T. k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields // Journal of biomedical optics. — 2010. — Vol. 15, no. 2. — P. 021314-021314.
[30] Ultrasound transmission computed tomography of the breast. / JS Schreiman, JJ Gisvold, James F Greenleaf, RC Bahn // Radiology. — 1984. — Vol. 150, no. 2. — P. 523-530.