Орбита форм и образов превращается в метрическое пространство путем наведения метрики на группу диффеоморфизмов. Изучение метрик на группах диффеоморфизмов и изучение метрик между многообразиями и поверхностями стало областью значительных исследований. [1] [2] [3] [4 ] [ 5] [6] [7] [8] [9] В вычислительной анатомии метрика диффеоморфометрии измеряет, насколько близко и далеко друг от друга находятся две формы или изображения. Неформально метрика строится путем определения потока диффеоморфизмов , которые соединяют элементы группы друг с другом, так что для тогда . Метрика между двумя системами координат или диффеоморфизмами тогда является кратчайшей длиной или геодезическим потоком , соединяющим их. Метрика на пространстве, связанном с геодезическими, задается как . Метрики на орбитах наследуются от метрики, индуцированной на группе диффеоморфизмов.
Часто знакомая евклидова метрика не применима напрямую, поскольку шаблоны форм и изображений не образуют векторное пространство . В модели римановой орбиты вычислительной анатомии диффеоморфизмы, действующие на формы, не действуют линейно. Существует много способов определения метрик, и для множеств, связанных с формами, метрика Хаусдорфа является другим. Метод, используемый для индукции римановой метрики , заключается в индукции метрики на орбите форм путем определения ее в терминах метрической длины между диффеоморфными преобразованиями систем координат потоков. Измерение длин геодезического потока между системами координат на орбите форм называется диффеоморфометрией .
Группа диффеоморфизмов, генерируемая как лагранжевы и эйлеровы потоки
с эйлеровыми векторными полями в для . Обратное для потока задается как
и матрица Якоби для потоков в задается как
Для обеспечения гладких потоков диффеоморфизмов с обратными векторными полями они должны быть по крайней мере однократно непрерывно дифференцируемыми в пространстве [10] [11] , которое моделируется как элементы гильбертова пространства с использованием теорем вложения Соболева, так что каждый элемент имеет производные, интегрируемые в квадрате 3, что подразумевает гладкие вложения в однократно непрерывно дифференцируемые функции. [10] [11] Группа диффеоморфизмов — это потоки с векторными полями, абсолютно интегрируемыми в норме Соболева:
( Группа Диффеоморфизмов )
Модель римановой орбиты
Формы в вычислительной анатомии (CA) изучаются с помощью диффеоморфного отображения для установления соответствий между анатомическими системами координат. В этой настройке трехмерные медицинские изображения моделируются как диффеоморфные преобразования некоторого образца, называемого шаблоном , в результате чего наблюдаемые изображения становятся элементами модели случайных орбит CA . Для изображений они определяются как , а для диаграмм, представляющих подмногообразия, обозначаются как .
Риманова метрика
Орбита форм и образов в Computational Anatomy генерируется действием группы , . Они превращаются в римановы орбиты путем введения метрики, связанной с каждой точкой и соответствующим касательным пространством. Для этого метрика определяется на группе, которая индуцирует метрику на орбите. Возьмем в качестве метрики для Computational anatomy в каждом элементе касательного пространства в группе диффеоморфизмов
с векторными полями, смоделированными в гильбертовом пространстве с нормой в гильбертовом пространстве . Мы моделируем как воспроизводящее ядро гильбертово пространство (RKHS), определяемое дифференциальным оператором 1-1 , где — дуальное пространство. В общем случае — обобщенная функция или распределение, линейная форма, связанная со скалярным произведением, и норма для обобщенных функций интерпретируется путем интегрирования по частям в соответствии с для ,
Когда , векторная плотность,
Дифференциальный оператор выбирается так, чтобы ядро Грина, связанное с обратным оператором, было достаточно гладким, чтобы векторные поля поддерживали 1-непрерывную производную . Аргументы теоремы вложения Соболева были приведены для демонстрации того, что 1-непрерывная производная требуется для гладких потоков. Оператор Грина , сгенерированный из функции Грина (скалярный случай), связанный с дифференциальным оператором, сглаживает.
Для правильного выбора тогда есть RKHS с оператором . Ядра Грина, связанные с дифференциальным оператором, сглаживают, поскольку для управления достаточным количеством производных в квадратично-интегральном смысле ядро непрерывно дифференцируемо по обеим переменным, что подразумевает
Диффеоморфометрия пространства форм и образов
Правоинвариантная метрика на диффеоморфизмах
Метрика на группе диффеоморфизмов определяется расстоянием, определенным на парах элементов в группе диффеоморфизмов согласно
( метрические диффеоморфизмы )
Это расстояние обеспечивает правоинвариантную метрику диффеоморфометрии, [12] [13] [14] инвариантную к перепараметризации пространства, поскольку для всех ,
Метрика форм и образов
Расстояние на изображениях, [15] ,
( метрические-формы-формы )
Расстояние по формам и видам, [16] ,
( метрические-формы-формы )
Метрика геодезических потоков ориентиров, поверхностей и объемов внутри орбиты
Для вычисления метрики геодезические представляют собой динамическую систему, поток координат и управление векторным полем связаны посредством Гамильтоново представление [17] [18] [19] [20] [21] перепараметризует распределение импульса в терминах гамильтонова импульса , множителя Лагранжа, ограничивающего лагранжеву скорость . Соответственно:
Принцип максимума Понтрягина [17] дает гамильтониан
Оптимизирующее векторное поле с динамикой . Вдоль геодезической гамильтониан постоянен: [22] . Метрическое расстояние между системами координат, связанными геодезической, определяется индуцированным расстоянием между единичным и групповым элементом:
Геодезические ориентиры или точки
Для ориентиров , , гамильтонов импульс
с динамикой Гамильтона, принимающей форму
с
Метрика между ориентирами
Динамика, связанная с этими геодезическими, показана на прилагаемом рисунке.
Геодезические исследования поверхности
Для поверхностей гамильтонов импульс определяется по всей поверхности, имеет гамильтониан
и динамика
Метрика между координатами поверхности
Объемная геодезическая
Для объемов гамильтониан
с динамикой
Метрика между объемами
Программное обеспечение для диффеоморфного отображения
Программные пакеты , содержащие различные алгоритмы диффеоморфного отображения, включают в себя следующее:
^ Миллер, MI; Юнес, Л. (2001-01-01). «Групповые действия, гомеоморфизмы и сопоставление: общая структура». International Journal of Computer Vision . 41 (1–2): 61–84. doi :10.1023/A:1011161132514. ISSN 0920-5691. S2CID 15423783.
^ Younes, L. (1998-04-01). «Вычислимые упругие расстояния между формами». Журнал SIAM по прикладной математике . 58 (2): 565–586. CiteSeerX 10.1.1.45.503 . doi :10.1137/S0036139995287685.
^ Mio, Вашингтон; Srivastava, Anuj; Joshi, Shantanu (2006-09-25). «О форме плоских упругих кривых». International Journal of Computer Vision . 73 (3): 307–324. CiteSeerX 10.1.1.138.2219 . doi :10.1007/s11263-006-9968-0. S2CID 15202271.
^ Michor, Peter W.; Mumford, David; Shah, Jayant; Younes, Laurent (2008). «Метрика в пространстве форм с явными геодезическими». Rend. Lincei Mat. Appl. () . 9 (2008): 25–57. arXiv : 0706.4299 . Bibcode :2007arXiv0706.4299M.
^ Michor, Peter W.; Mumford, David (2007). «Обзор римановых метрик на пространствах кривых с использованием гамильтоновского подхода». Applied and Computational Harmonic Analysis . 23 (1): 74–113. arXiv : math/0605009 . doi :10.1016/j.acha.2006.07.004. S2CID 732281.
^ Куртек, Себастьян; Классен, Эрик; Гор, Джон К.; Дин, Чжаохуа; Шривастава, Анудж (2012-09-01). «Эластичные геодезические пути в пространстве форм параметризованных поверхностей». Труды IEEE по анализу шаблонов и машинному интеллекту . 34 (9): 1717–1730. doi :10.1109/TPAMI.2011.233. PMID 22144521. S2CID 7178535.
^ Шривастава, Анудж; Классен, Эрик; Джоши, Шантану Х.; Джермин, Ян Х. (2011). «Анализ формы упругих кривых в евклидовых пространствах». Труды IEEE по анализу шаблонов и машинному интеллекту . 33 (7): 1415–1428. doi :10.1109/TPAMI.2010.184. ISSN 1939-3539. PMID 20921581. S2CID 12578618.
^ Jermyn, Ian H.; Kurtek, Sebastian; Klassen, Eric; Srivastava, Anuj (2012), Fitzgibbon, Andrew; Lazebnik, Svetlana; Perona, Pietro; Sato, Yoichi (ред.), "Elastic Shape Matching of Parameterized Surfaces Using Square Root Normal Fields", Computer Vision – ECCV 2012 , т. 7576, Берлин, Гейдельберг: Springer Berlin Heidelberg, стр. 804–817, doi : 10.1007/978-3-642-33715-4_58 , ISBN978-3-642-33714-7
^ Джермин, Ян Х.; Куртек, Себастьян; Лага, Хамид; Шривастава, Анудж (15.09.2017). «Анализ упругой формы трехмерных объектов». Synthesis Lectures on Computer Vision . 7 (3): 1–185. doi :10.2200/s00785ed1v01y201707cov012. ISSN 2153-1056. S2CID 52096321.
^ ab P. Dupuis, U. Grenander, MI Miller, Существование решений для потоков диффеоморфизмов, Quarterly of Applied Math, 1997.
^ аб А. Труве. Действия группы бесконечных измерений и разведка форм. CR Acad Sci Paris Ser I Math, 321(8):1031–1034, 1995.
^ Миллер, MI; Юнес, Л. (2001-01-01). «Групповые действия, гомеоморфизмы и сопоставление: общая структура». Международный журнал компьютерного зрения . 41 : 61–84. CiteSeerX 10.1.1.37.4816 . doi :10.1023/A:1011161132514. S2CID 15423783.
^ Миллер, М. И.; Юнес, Л.; Труве, А. (2014). «Диффеоморфометрия и геодезические системы позиционирования для анатомии человека». Технологии . 2 (1): 36. doi :10.1142/S2339547814500010. PMC 4041578. PMID 24904924 .
^ Миллер, Майкл И.; Труве, Ален; Юнес, Лоран (2015-01-01). «Гамильтоновы системы и оптимальное управление в вычислительной анатомии: 100 лет со дня Д'Арси Томпсона». Ежегодный обзор биомедицинской инженерии . 17 (1): 447–509. doi :10.1146/annurev-bioeng-071114-040601. PMID 26643025.
^ Миллер, MI; Юнес, Л. (2001-01-01). «Групповые действия, гомеоморфизмы и сопоставление: общая структура». Международный журнал компьютерного зрения . 41 : 61–84. CiteSeerX 10.1.1.37.4816 . doi :10.1023/A:1011161132514. S2CID 15423783.
^ Миллер, Майкл И.; Юнес, Лоран; Труве, Ален (март 2014 г.). «Диффеоморфометрия и геодезические системы позиционирования для анатомии человека». Технологии . 2 (1): 36. doi :10.1142/S2339547814500010. ISSN 2339-5478. PMC 4041578. PMID 24904924 .
^ ab Миллер, Майкл И.; Труве, Ален; Юнес, Лоран (2015-01-01). «Гамильтоновы системы и оптимальное управление в вычислительной анатомии: 100 лет со дня Дарси Томпсона». Ежегодный обзор биомедицинской инженерии . 17 (1): 447–509. doi :10.1146/annurev-bioeng-071114-040601. PMID 26643025.
^ Glaunès J, Trouvé A, Younes L. 2006. Моделирование изменения плоской формы с помощью гамильтоновых потоков кривых. В Statistics and Analysis of Shapes, ред. H Krim, A Yezzi Jr, стр. 335–61. Model. Simul. Sci. Eng. Technol. Boston: Birkhauser
^ Аргиллер С., Трела Э., Труве А., Юнес Л. 2014. Анализ деформации формы с точки зрения оптимального управления. arXiv : 1401.0661 [math.OC]
^ Миллер, MI; Юнес, L; Труве, A (2014). «Диффеоморфометрия и геодезические системы позиционирования для анатомии человека». Technology (Singap World Sci) . 2 (1): 36. doi :10.1142/S2339547814500010. PMC 4041578. PMID 24904924 .
^ Michor, Peter W.; Mumford, David (2007-07-01). "Обзор римановых метрик на пространствах кривых с использованием гамильтоновского подхода". Applied and Computational Harmonic Analysis . Специальный выпуск по математической визуализации. 23 (1): 74–113. arXiv : math/0605009 . doi :10.1016/j.acha.2006.07.004. S2CID 732281.
^ Миллер, Майкл И.; Труве, Ален; Юнес, Лоран (2015-01-01). «Гамильтоновы системы и оптимальное управление в вычислительной анатомии: 100 лет со дня Д'Арси Томпсона». Ежегодный обзор биомедицинской инженерии . 17 (1): 447–509. doi :10.1146/annurev-bioeng-071114-040601. PMID 26643025.
^ Программное обеспечение - Стэнли Дуррлеман (Отчет).
^ Avants, Brian B.; Tustison, Nicholas J.; Song, Gang; Cook, Philip A.; Klein, Arno; Gee, James C. (2011-02-01). «Воспроизводимая оценка показателей подобия ANT при регистрации изображений мозга». NeuroImage . 54 (3): 2033–2044. doi :10.1016/j.neuroimage.2010.09.025. ISSN 1053-8119. PMC 3065962 . PMID 20851191.
^ Эшбёрнер, Джон (15.10.2007). «Быстрый алгоритм регистрации диффеоморфных изображений». NeuroImage . 38 (1): 95–113. doi :10.1016/j.neuroimage.2007.07.007. PMID 17761438. S2CID 545830.
^ "Программное обеспечение - Том Веркаутерен" . сайты.google.com . Проверено 11 декабря 2015 г.
^ Бег, М. Фейсал; Миллер, Майкл И.; Труве, Ален; Юнес, Лоран (2005-02-01). «Вычисление метрических отображений больших деформаций с помощью геодезических потоков диффеоморфизмов». Международный журнал компьютерного зрения . 61 (2): 139–157. doi :10.1023/B:VISI.0000043755.93987.aa. ISSN 0920-5691. S2CID 17772076.