stringtranslate.com

МРТ в реальном времени

МРТ сердца человека в реальном времени (2-камерная проекция) с разрешением 22 мс [1]
МРТ речевого тракта в реальном времени во время пения с разрешением 40 мс.

Магнитно-резонансная томография в реальном времени ( RT-MRI ) относится к непрерывному мониторингу («съемке») движущихся объектов в реальном времени. Поскольку МРТ основана на трудоемком сканировании k-пространства , МРТ в реальном времени была возможна только с низким качеством изображения или низким временным разрешением. С помощью алгоритма итеративной реконструкции эти ограничения недавно были устранены: новый метод МРТ в реальном времени обеспечивает временное разрешение от 20 до 30 миллисекунд для изображений с разрешением в плоскости от 1,5 до 2,0 мм. [2] МРТ в реальном времени обещает добавить важную информацию о заболеваниях суставов и сердца . Во многих случаях МРТ-обследования могут стать для пациентов проще и комфортнее.

История

1977/1978 - Раймонд Дамадиан построил первый сканер МРТ и провел первое МРТ-сканирование здорового человеческого тела (1977) с целью диагностики рака. [3] Кроме того, Питер Мэнсфилд разрабатывает эхо-планарную технику, позволяющую получать изображения за считанные секунды и становящуюся основой для быстрой МРТ. [4]

1983 — Введение k-пространства Д.Б. Твигом [5].

1987 г. – разработана первая МРТ сердца в реальном времени [6].

1997 - Д.К. Содиксон представил параллельную визуализацию с помощью массива радиочастотных катушек [7].

1999 г. - К.П. Прюсманн представил реконструкцию изображений SENSE [8].

2002 — Марк Грисволд представил реконструкцию изображений GRAPPA [9].

Физическая основа

Обзор

Уравнение для расчета скорости нарастания и прилагаемая диаграмма

В целом, МРТ в реальном времени опирается на последовательности градиентного эха , эффективную выборку k-пространства и методы быстрой реконструкции для ускорения процесса получения изображения. [10] Последовательности градиентного эха обеспечивают более короткое время эхо-сигнала, поскольку для каждой последовательности требуется только один РЧ- импульс. [11] Современные градиентные катушки с быстрым переключением также требуют увеличения скорости нарастания напряжения , что позволяет быстрее изменять последовательности градиентного эхо и уменьшать время повторения . [12]

Прямолинейная выборка k-пространства

выборка в k-пространстве

Другие траектории выборки k-пространства

Эффективная выборка k-пространства также сокращает время сбора данных. Прямолинейное сканирование стало стандартным методом отбора проб в k-пространстве для МРТ. [13] Однако этот процесс занимает относительно много времени, поскольку он равномерно оцифровывает все k-пространство. Из-за этой задержки для захвата движения в реальном времени используются другие методы выборки. Планарная визуализация одиночного эха — это чрезвычайно быстрый метод выборки, при котором все данные для МР-изображения собираются из одного радиочастотного импульса. [14] Однако важно отметить, что метод EPI по-прежнему является декартовым методом выборки, как и прямолинейное сканирование, одинаково осуществляющим выборку всего k-пространства. Спиральная выборка, как и EPI, требует только одного радиочастотного импульса для выборки всего k-пространства. Радиальная и спиральная выборка также используются в качестве методов эффективной выборки k-пространства, при этом спираль также требует только одного радиочастотного импульса для выборки k-пространства. Как радиальная, так и спиральная выборка более эффективны, чем декартовы методы, поскольку они передискретизируют низкие частоты, что позволяет обеспечить общий захват движения и лучшую реконструкцию изображения в реальном времени. [10] Таким образом, радиальная или спиральная выборка k-пространства в настоящее время является предпочтительным методом МРТ-реконструкции в реальном времени.

Параллельная визуализация

Коллекция катушек параллельного формирования изображений

Параллельная визуализация предполагает добавление нескольких катушек, окружающих цель, при этом каждая катушка получает часть общего изображения. Поскольку современные графические процессоры обладают возможностями параллельной обработки, они могут реконструировать каждую часть изображения одновременно. Таким образом, чем больше катушек используется, тем быстрее происходит получение МР-изображений. [15]

Последовательности градиентного эха

ФЛЭШ-МРТ

В то время как ранние приложения были основаны на эхопланарной визуализации, которая нашла важное применение в функциональной МРТ в реальном времени (rt-fMRI), [16] недавний прогресс основан на итеративной реконструкции и FLASH MRI . [17] [18] Метод визуализации в реальном времени, предложенный Юкером и его коллегами [2], сочетает в себе радиальную FLASH MRI, [19] которая обеспечивает быстрый и непрерывный сбор данных, устойчивость к движению и устойчивость к недостаточной выборке, с итеративным методом реконструкции изображения. основан на постановке реконструкции изображения как нелинейной обратной задачи . [20] [21] Путем интеграции данных с нескольких приемных катушек (т.е. параллельного МРТ) и использования избыточности во временных рядах изображений с использованием регуляризации и фильтрации этот подход увеличивает возможную степень недодискретизации данных на один порядок. величину, так что высококачественные изображения могут быть получены всего лишь из 5–10% данных, необходимых для нормальной реконструкции изображения.

Из-за очень короткого времени эхо-сигнала (например, от 1 до 2 миллисекунд ) метод не страдает от нерезонансных эффектов, так что изображения не демонстрируют артефактов чувствительности и не зависят от подавления жира. В то время как испорченные последовательности FLASH обеспечивают плотность вращения или контраст T1, версии с перефокусированными или полностью сбалансированными градиентами обеспечивают доступ к контрасту T2/T1. Выбор времени градиентного эхо (например, синфазные или противофазные условия) дополнительно изменяет представление сигналов воды и жира на изображениях и позволяет создавать отдельные фильмы воды/жира.

Сбалансированная устойчивая свободная прецессия

Другая последовательность GRE, обычно используемая в RT-MRI, - это сбалансированная стационарная свободная прецессия (bSSFP), как упоминалось выше, со сбалансированными градиентами. [10] Свободная прецессия в устойчивом состоянии предполагает время повторения (TR), которое короче T2. Это предотвращает полное затухание магнитного сигнала до подачи следующего радиочастотного импульса, что затем с течением времени устанавливает устойчивый сигнал. [22] Короткий TR также делает bSSFP идеальным для RT-MRI.

Уравнение для пикового сигнала MR в bSSFP имеет вид:

Где – начальная намагниченность, и .

Таким образом, сигнал МР пропорционален T2/T1. Материалы с одинаковыми T1 и T2, такие как жидкости и жир, демонстрируют высокий контраст T2/T1 и могут иметь интенсивность сигнала до

Сигнал bSSFP также превышает сигнал FLASH в несколько раз.

. [22]

Due to this strong fluid/tissue contrast, RT-MRI with bSSFP lends itself to cardiac imaging and visualizing blood flow.[22]

Image reconstruction

SENSitivity Encoding (SENSE)

Certain image reconstruction algorithms used alongside parallel imaging address the potential issues that can arise from undersampling the k-space. SENSitivity Encoding (SENSE) is a method that reconstructs the partial k-space data from each coil and combines the partial images into the final scan in the spatial domain.[23] Coil sensitivities must first be acquired either before the actual imaging or during the imaging process. During the rest of imaging, the k-space is undersampled to skip every other line, resulting in a ½ FOV.

As a two-point example, pixels on the original aliased images can be “unfolded” through the following equations to give the final scan:

for two points, and , in the final image. and denote the image signal for the aliased image. and are the sensitivity values for coil 1 at points and , respectively, and and are the sensitivity values for coil 2 at points and , respectively.[23]

GeneRalized Autocalibrating Partial Parallel Acquisition (GRAPPA)

Another reconstruction algorithm used is GeneRalized Autocalibrating Partial Parallel Acquisition (GRAPPA). GRAPPA fills in the undersampled k-space data in the k-space domain before reconstructing the final image.[24] Lines through the center of the k-space are fully sampled, typically alongside the actual image, to give the autocalibration signal (ACS) region. Weighing factors are calculated using the ACS, and these factors reflect the coil-specific distortions that each coil applies on the full field-of-view frequency domain. Then, the filled-in k-space data undergoes the inverse Fourier transform to construct the partial, non-aliased images. These images are then simply combined directly in the spatial domain.[24]

If the k-space data is non-Cartesian, reconstruction is computationally more difficult, since the fast Fourier transform (FFT) requires Cartesian values. Typically, k-space data must be resampled into Cartesian coordinates before applying the FFT. GRAPPA can address these issues by obtaining large quantities of calibration data; however, the fastest reconstructions will generally require Cartesian data.[10]

Signal-to-noise ratio

Lastly, within parallel image reconstruction there is another factor to consider, which is the signal to noise ratio (SNR). The SNR for parallel imaging can be calculated using the following equation:

[25]

Где – коэффициент ускорения, а – пространственно-зависимый геометрический фактор (пропорциональный количеству используемых катушек или взаимодействию между катушками). Следовательно, чем больше катушек используется, тем быстрее процесс визуализации и тем больше взаимодействий между катушками; следовательно, тем ниже SNR. [25]

Приложения

(a,c) CINE и (b,d) RT CMR при 1,5 Тл у пациента с фибрилляцией предсердий в проекции SA среднего желудочка во время конечной диастолы (a,b) и конца систолы (c,d). CINE-маги имеют пониженное диагностическое качество (оценка 1), тогда как RT-изображения имеют хорошее диагностическое качество.

МРТ сердца

Хотя применения МРТ в реальном времени охватывают широкий спектр — от немедицинских исследований турбулентного потока [26] до неинвазивного мониторинга интервенционных (хирургических) процедур, наиболее важным применением новых возможностей является визуализация сердечно-сосудистой системы. [1] Предыдущая МРТ сердца (ММР) использовала кинотехнику для фиксации периодических движений сердца. Однако это неосуществимо для пациентов с аритмией , у которых сердечный цикл непредсказуем. [27] С помощью нового метода можно получать фильмы бьющегося сердца в реальном времени со скоростью до 50 кадров в секунду при свободном дыхании и без необходимости синхронизации с электрокардиограммой . [28] Исследование, проведенное Laubrock et al. [27] продемонстрировали, что RT-MRI дает изображения более высокого качества с более высоким SNR, чем киноCMR с последовательностью bSSFP и радиальной выборкой в ​​k-пространстве. RT-MRI также устраняет необходимость задержки дыхания во время визуализации, что также повышает комфорт пациента. [27]

МРТ опорно-двигательного аппарата

Помимо МРТ сердца, другие приложения, работающие в режиме реального времени, занимаются функциональными исследованиями кинетики суставов (например, височно-нижнечелюстного сустава , [29] колена и запястья [30] ) или изучают скоординированную динамику артикуляторов, таких как губы, язык, мягкое небо и голосовые складки во время разговора ( артикуляционная фонетика ) [31] или глотания . [32] Визуализация скелетно-мышечной системы особенно полезна при наблюдении в реальном времени. Исследователи из Медицинской школы имени Гроссмана Нью-Йоркского университета [33] разработали перчатку RT-MRI для визуализации движений руки. В перчатке используются катушки с высоким импедансом для предотвращения генерации вихревых токов из-за быстро меняющихся магнитных полей и bSSFP для быстрого получения изображений. Катушки с высоким импедансом устраняют необходимость в специальной конформации катушек и активном градиентном экранировании. [34]

Инвазивные процедуры под контролем МРТ

Применение в интервенционной МРТ , которая относится к мониторингу минимально инвазивных хирургических процедур, возможно путем интерактивного изменения таких параметров, как положение и ориентация изображения. Это приложение особенно полезно, когда во время операции требуется трехмерное изображение ткани. [10] Для этого требуется наличие дисплея в кабинете, который врач сможет использовать во время процедуры, а также использование хирургических инструментов, безопасных для МРТ. К ним относятся керамика, пластик или титан, который является парамагнитным металлом. Благодаря использованию bSSFP и параллельной визуализации с несколькими катушками была достигнута частота кадров 5–10 кадров в секунду, что позволяет визуализировать сердечные процедуры. [10]

Будущие направления

Параллельная визуализация

Катушки для параллельной визуализации доступны для визуализации туловища и сердца, но они еще не стандартизированы для других частей тела. Динамические катушки для визуализации речи и скелетно-мышечной системы являются ключевыми областями текущих исследований. [10]

Стандартная реконструкция (слева) и реконструкция с DL (справа) для Т2-взвешенного турбоспин-эхосканирования предстательной железы

Машинное обучение

Реконструкция изображений в RT-MRI выигрывает от машинного обучения (ML) или глубокого обучения (DL). Нелинейное ядро ​​или функция отображения может быть разработано на основе ACS для заполнения данных k-пространства и создания окончательного изображения. [35] Этот процесс в целом значительно ускоряет процесс МРТ. Сегментация изображений или идентификация повреждений могут быть достигнуты с помощью машинного обучения. При глубоком обучении с помощью сверточной нейронной сети функция отображения может быть задана сетью. ML и DL улучшают разрешение изображения, а также скорость обработки изображений. [35]

Высокопроизводительные сканеры с низким полем зрения

Высокопроизводительные МРТ-сканеры с низким полем также находятся в области разработки. [10] Эти сканеры работают при относительно низкой напряженности магнитного поля, например 0,35 Тл или 0,55 Тл. Многие последовательности сбора данных RT-MRI, такие как bSSFP, испытывают значительные внерезонансные эффекты. Внерезонансные эффекты линейно возрастают с увеличением напряженности поля B0, поэтому минимизация B0 также сводит к минимуму эти эффекты, которые могут привести к артефактам и искажению изображения. [36] Это позволяет использовать более длинные TR, что затем открывает двери для более широкого спектра методов выборки в k-пространстве и конструкций последовательностей. [10] Наконец, МРТ-сканеры меньшей мощности снизят опасность, связанную с нагревом металлических имплантатов, и снизят стоимость МРТ. [37]

Рекомендации

  1. ^ ab S Zhang, M Uecker, D Voit, KD Merboldt, J Frahm (2010a) Сердечно-сосудистый магнитный резонанс в реальном времени с высоким временным разрешением: радиальная FLASH с нелинейной обратной реконструкцией. J Cardiovasc Magn Reson 12, 39, [1] doi : 10.1186/1532-429X-12-39
  2. ^ ab M Uecker, S Zhang, D Voit, A Karaus, KD Merboldt, J Frahm (2010a) МРТ в реальном времени с разрешением 20 мс. ЯМР Biomed 23: 986-994, [2] doi : 10.1002/nbm.1585.
  3. ^ «Кто создал Америку? | Новаторы | Раймонд Дамадиан». www.pbs.org . Проверено 9 декабря 2022 г.
  4. ^ "Призывник Национального Зала славы изобретателей и изобретатель МРТ Питер Мэнсфилд" . www.invent.org . Проверено 9 декабря 2022 г.
  5. ^ Твиг, Д.Б. (1983). «Формулировка k-траектории процесса ЯМР-визуализации с применением в анализе и синтезе методов визуализации». Медицинская физика . 10 (5): 610–621. Бибкод : 1983MedPh..10..610T. дои : 10.1118/1.595331. ISSN  0094-2405. ПМИД  6646065.
  6. ^ «Краткая история МРТ | Два взгляда» . два-views.com . Проверено 9 декабря 2022 г.
  7. ^ Содиксон, ДК; Мэннинг, WJ (октябрь 1997 г.). «Одновременное получение пространственных гармоник (SMASH): быстрое получение изображений с помощью массивов радиочастотных катушек». Магнитный резонанс в медицине . 38 (4): 591–603. дои : 10.1002/mrm.1910380414. ISSN  0740-3194. PMID  9324327. S2CID  17505246.
  8. ^ Прюсманн, КП; Вайгер, М.; Шайдеггер, МБ; Бозигер, П. (ноябрь 1999 г.). «SENSE: кодирование чувствительности для быстрой МРТ». Магнитный резонанс в медицине . 42 (5): 952–962. doi : 10.1002/(SICI)1522-2594(199911)42:5<952::AID-MRM16>3.0.CO;2-S . ISSN  0740-3194. PMID  10542355. S2CID  16046989.
  9. ^ Грисволд, Марк А.; Якоб, Питер М.; Хайдеманн, Робин М.; Ниттка, Матиас; Желлюс, Владимир; Ван, Цзяньминь; Кифер, Бертольд; Хаазе, Аксель (июнь 2002 г.). «Обобщенная автокалибровка, частично параллельный сбор данных (GRAPPA)». Магнитный резонанс в медицине . 47 (6): 1202–1210. дои : 10.1002/mrm.10171 . ISSN  0740-3194. PMID  12111967. S2CID  14724155.
  10. ^ abcdefghi Наяк, Кришна С.; Лим, Юнван; Кэмпбелл-Уошберн, Эдриен Э.; Стиден, Дженнифер (январь 2022 г.). «Магнитно-резонансная томография в реальном времени». Журнал магнитно-резонансной томографии . 55 (1): 81–99. дои : 10.1002/jmri.27411. ISSN  1053-1807. ПМЦ 8435094 . ПМИД  33295674. 
  11. ^ «Градиентное эхо (GRE)» . Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  12. ^ «Характеристики градиента». Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  13. ^ «Радиальная выборка». Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  14. ^ Пущи-Амин, Мехди; Мировиц, Скотт А.; Браун, Джеффри Дж.; МакКинстри, Роберт С.; Ли, Тао (1 мая 2001 г.). «Принципы и применение эхопланарной визуализации: обзор для врача-рентгенолога». Радиографика . 21 (3): 767–779. doi : 10.1148/radiographics.21.3.g01ma23767. ISSN  0271-5333. ПМИД  11353123.
  15. ^ «Параллельная визуализация (PI)» . Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  16. ^ Коэн М.С. (2001). «Функциональная магнитно-резонансная томография в реальном времени». Методы . 25 (2): 201–220. дои : 10.1006/meth.2001.1235. ПМИД  11812206.
  17. ^ Дж. Фрам , А. Хаазе, В. Хенике, К. Д. Мербольдт, Д. Маттеи (1985) Hochfrequenz-Impuls und Gradienten-Impuls-Verfahren zur Aufnahme von schnellen NMR-Tomogrammen unter Benutzung von Gradientenechos. Заявка на патент Германии P 35 04 734.8, 12 февраля 1985 г.
  18. ^ Дж. Фрам , А. Хаазе, Д. Маттеи (1986) Быстрая ЯМР-визуализация динамических процессов с использованием техники FLASH. Magn Reson Med 3:321-327 [3] doi :10.1002/mrm.1910030217
  19. ^ С. Чжан, КТ Блок КТ, Дж. Фрам (2010b) Магнитно-резонансная томография в реальном времени: достижения в использовании радиальной FLASH. J Magn Reson Imag 31: 101-109, [4] doi :10.1002/jmri.21987
  20. ^ М. Юкер, Т. Хохаге, К. Т. Блок, Дж. Фрам (2008) Реконструкция изображения посредством регуляризованной нелинейной инверсии - Совместная оценка чувствительности катушки и содержания изображения. Magn Reson Med 60: 674-682, [5] doi : 10.1002/mrm.21691.
  21. ^ М. Юкер, С. Чжан, Дж. Фрам (2010b) Нелинейная обратная реконструкция для МРТ человеческого сердца в реальном времени с использованием радиальной FLASH с недостаточной выборкой. Magn Reson Med 63: 1456-1462, [6] doi :10.1002/mrm.22453
  22. ^ abc Биери, Оливер; Шеффлер, Клаус (июль 2013 г.). «Основы сбалансированной стационарной МРТ со свободной прецессией: основы сбалансированной МРТ SSFP». Журнал магнитно-резонансной томографии . 38 (1): 2–11. дои : 10.1002/jmri.24163 . PMID  23633246. S2CID  23442745.
  23. ^ ab "ЧУВСТВО/АКТИВ?". Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  24. ^ ab "GRAPPA/ARC?". Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  25. ^ ab «Параллельная визуализация (PI): шум». Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  26. ^ В. Ведин; Кроули; Р. Вайскофф; Дж. Холмванг и М. С. Коэн (1990). «МРТ потока структурированной жидкости в реальном времени». Общество магнитного резонанса в медицине : 164.
  27. ^ abc Лауброк, Керстин; фон Леш, Тассило; Стейнмец, Майкл; Лотц, Иоахим; Фрам, Йенс; Юкер, Мартин; Унтерберг-Бухвальд, Кристина (1 января 2022 г.). «Визуализация аритмии: магнитно-резонансная томография сердца в режиме реального времени при фибрилляции предсердий». Европейский журнал радиологии Open . 9 : 100404. doi : 10.1016/j.ejro.2022.100404. ISSN  2352-0477. ПМЦ 8899235 . ПМИД  35265735. 
  28. ^ И. Уяник, П. Линднер, Д. Шах, Н. Цекос и Павлидис (2013) Применение метода набора уровней для разрешения физиологических движений при МРТ сердца со свободным дыханием и без ворот. FIMH, 2013, «Лаборатория вычислительной физиологии» (PDF) . Архивировано из оригинала (PDF) 22 июля 2018 г. Проверено 1 октября 2013 г.
  29. ^ С. Чжан, Н. Герсдорф, Дж. Фрам (2011) Магнитно-резонансная томография динамики височно-нижнечелюстного сустава в реальном времени. Открытый журнал медицинской визуализации, 2011, 5, 1-7, «Магнитно-резонансная томография динамики височно-нижнечелюстного сустава в реальном времени» (PDF) . Архивировано из оригинала (PDF) 26 сентября 2012 г. Проверено 16 сентября 2011 г.
  30. ^ Бутин Р.Д., Буонокоре М.Х., Иммерман И., Эшвелл З., Сонико Г.Дж., Сабо Р.М. и Чаудхари А.Дж. (2013) Магнитно-резонансная томография (МРТ) в реальном времени во время активного движения запястья — первоначальные наблюдения. PLoS ONE 8(12): e84004. doi:10.1371/journal.pone.0084004
  31. ^ Нибергалл А., Чжан С., Кунай Э., Кейдана Г., Джоб М. и др. МРТ речи в реальном времени с разрешением 33 мс: радиальная FLASH с недостаточной выборкой и нелинейной обратной реконструкцией. Magn Reson Med 2010, номер документа : 10.1002/mrm.24276.
  32. ^ Чжан С., Олтхофф А. и Фрам Дж. Магнитно-резонансная томография нормального глотания в реальном времени. J Magn Reson Imaging 2011; 35: 1372-1379. дои : 10.1002/jmri.23591.
  33. ^ Чжан, Бэй; Содиксон, Дэниел К.; Клоос, Мартейн А. (август 2018 г.). «Высокоимпедансная детекторная перчатка для магнитно-резонансной томографии руки». Природная биомедицинская инженерия . 2 (8): 570–577. дои : 10.1038/s41551-018-0233-y. ISSN  2157-846X. ПМК 6405230 . ПМИД  30854251. 
  34. ^ «Активно экранированные градиенты» . Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  35. ^ аб Гассенмайер, Себастьян; Кюстнер, Томас; Никель, Доминик; Херрманн, Юдит; Хоффманн, Рюдигер; Альмансур, Хайдара; Афат, Саиф; Николау, Константин; Осман, Ахмед Э. (24 ноября 2021 г.). «Приложения глубокого обучения в магнитно-резонансной томографии: стало ли будущее настоящим?». Диагностика . 11 (12): 2181. doi : 10.3390/diagnostics11122181 . ISSN  2075-4418. ПМК 8700442 . ПМИД  34943418. 
  36. ^ «Внерезонансные импульсы». Вопросы и ответы по МРТ . Проверено 9 декабря 2022 г.
  37. ^ Кэмпбелл-Уошберн, Адриенн Э.; Рамасавми, Раджив; Рестиво, Мэтью С.; Бхаттачарья, Ипшита; Басар, Бурку; Герцка, Дэниел А.; Хансен, Майкл С.; Роджерс, Тоби; Бандеттини, В. Патрисия; МакГирт, Делани Р.; Манчини, Кристина; Гродский, Дэвид; Шнайдер, Райнер; Маджид, Вакас; Бхат, Химаншу (ноябрь 2019 г.). «Возможности интервенционной и диагностической визуализации с использованием высокоэффективной МРТ с низкой напряженностью поля». Радиология . 293 (2): 384–393. дои : 10.1148/radiol.2019190452. ISSN  1527-1315. ПМК 6823617 . ПМИД  31573398. 

Внешние ссылки