Цифровая технология для диагностики ретинопатии недоношенных

ЛИТЕРАТУРА

1. Асташева И. Б., Сидоренко Е. И., Аксенова И. И. Лазеркоа-гуляция в лечении различных форм ретинопатии недоношенных // Вестник офтальмологии. - 2005. - № 2. - С. 31-34.

2. Коротких С. А., Степанова Е. А., Кулакова М. В. Злокачественная ретинопатия недоношенных: особенности перинатального анамнеза, клиники, эффективность лазерного лечения // Российская педиатрическая офтальмология. - 2007. -№ 4. - С. 13-14.

3. Терещенко А. В., Белый Ю. А., Терещенкова М. С., Володин П. Л. Особенности диагностики и лечения задней агрессивной ретинопатии недоношенных // Офтальмохирургия. - 2010. - № 1. -С. 29-33.

4. Терещенко А. В., Белый Ю. А., Трифаненкова И. Г., Терещенкова М. С. Анализ состояния сосудов сетчатки в прогнозировании течения ретинопатии недоношенных // Офтальмохирургия. -2006. - № 3. - С. 37-40.

5. Altinsoy H., Mutlu F., Gungor R. Combination of laser photocoagulation and intravitreal bevacizumab in aggressive posterior retinopathy of prematurity // Ophthalmic surg lasers imaging. - 2010. -Vol. 9. - Р. 1-5.

6. Axer-Siegel R., Maharshak I., Snir M. Diode laser treatment of retinopathy of prematurity: anatomical and refractive outcomes // Retina. - 2008. - Vol. 28. №. 6. - Р. 839-846.

7. Blumenkranz M., Yellachich D., Andersen D. et al. Semiautomated patterned scanning laser for retinal photocoagulation // Retina. -2006. - Vol. 26. - Р. 370-376.

8. Buksh M., Dai S., Kuschel C. AP-ROP in an infant with minimal oxygen exposure // J. paediatr child health. - 2008. - Vol. 44. №. 4. -Р. 228-230.

9. Hiraoka M., Nishina S., Nakagawa A. Case of aggressive posterior retinopathy of prematurity with atypical neovascular growth // Jpn. j. ophthalmol. - 2008. - Vol. 52. №. 5. - Р. 417-419.

10. Honda S., Hirabayashi H., Tsukahara Y. Acute contraction of the proliferative membrane after an intravitreal injection of bevacizumab for advanced retinopathy of prematurity // Graefes arch. clin. exp. ophthalmol. - 2008. - Vol. 246. № 7. - Р. 1061-1063.

11. Kim M., Kim S., Yu Y. The risk for retinal detachment associated with hemorrhages pre- and postlaser treatment in retinopathy of prematurity // Retina. - 2008. - Vol. 28. № 10. - Р. 1451-1457.

12. Law J., Recchia F., Morrison D. Intravitreal bevacizumab as adjunctive treatment for retinopathy of prematurity // J AAPOS. - 2010. -Vol. 14. № 1. - Р. 6-10.

13. Sanghi G., Dogra M., Das P. Aggressive posterior retinopathy of prematurity in Asian Indian babies: spectrum of disease and outcome after laser treatment // Retina. - 2009. - Vol. 29. № 9. - Р. 1335-1339.

14. Sanghvi C., Mclauchlan R., Delgado C. Initial experience with the Pascal photocoagulator: a pilot study of 75 procedures // Br. j. ophthalmol. - 2008. - Vol. 92. - P. 1061-1064.

15. The International classification of retinopathy of prematurity revisited // Arch. ophthalmol. - 2005. - Vol. 123. № 7. -Р. 991-999.

Поступила 25.09.2010

А. В. ТЕРЕЩЕНКО, Ю. А. БЕЛЫЙ, М. С. ТЕРЕЩЕНКОВА, И. Г. ТРИФАНЕНКОВА, Ю. А. ЮДИНА

ЦИФРОВАЯ ТЕХНОЛОГИЯ ДЛЯ ДИАГНОСТИКИ РЕТИНОПАТИИ НЕДОНОШЕННЫХ

Калужский филиал ФГУ «МНТК «Микрохирургия глаза» им. академика С. Н. Федорова Росмедтехнологии»,

Россия, 248007, г. Калуга, ул. Вишневского, 1а. E-mail: nauka@mntk.kaluga.ru

Разработано программное обеспечение для количественной оценки состояния сетчатки и ретинальных сосудов при ретинопатии недоношенных (РН). Программа обеспечивает выполнение следующих функций: проецирование двухмерных цифровых изображений глазного дна на поверхность виртуальной трехмерной сферы; выполнение заданных пользователем перемещений; вычисление диаметра и извитости ретинальных сосудов, угла между сосудистыми аркадами, площади васку-ляризированной и аваскулярной сетчатки. Получаемые данные необходимы для достоверного прогноза течения активных стадий РН и оценки эффективности проводимого лечения.

Ключевые слова: ретинопатия недоношенных, трехмерная модель глазного яблока, количественный анализ сетчатки и ретинальных сосудов.

А. V. TERESHCHENKO, Yu. А. BELYY, М. S. TERESHCHENKOVA,

I. G. TRIFANENKOVA, Yu. А. YUDINA

DIGITAL TECHNOLOGY FOR RETINOPATHY OF PREMATURITY DIAGNOSTICS

The S. Fyodorov Eye microsurgery Federal state institution,

Russia, 248007, Kaluga, Vishnevskogo street, 1а. E-mail: nauka@mntk.kaluga.ru

The soft for quantification of retinopathy of prematurity (ROP) via retina and retinal vessels parameters was designed. The soft provides the next functions: two-dimensional digital fundus photos projection on the virtual sphere; any photos movements; vessels diameter and tortuosity calculations, an angle between the vascular arcades and vascularized and avascular retina square calculations. The data that the soft provides is necessary for reliable prognosis of ROP active stages course and to assess an efficacy of the treatment.

Key words: retinopathy of prematurity, three-dimension eye model, retina and retinal vessels quantification.

Кубанский научный медицинский вестник № 1 (124) 2011 УДК 617-753-089

Кубанский научный медицинский вестник № 1 (124) 2011

Развитие высокоинформативных методов диагностики ретинопатии недоношенных (РН), без которых невозможно проведение своевременного и адекватного лечения, связано с разработкой и внедрением в клиническую практику современных систем для получения высококачественных изображений глазного дна.

В 1998 году была создана широкопольная ретинальная цифровая педиатрическая видеосистема «RetCam-120» («Massie Research Laboratories Inc.», США), позволяющая проводить цифровую ретиноско-пию и регистрировать в динамике все стадии и клинические проявления РН [1].

В 1999 году в Германии появилось первое сообщение об использовании «RetCam-120» для скрининга недоношенных младенцев и, кроме того, представлены возможности ретинальной системы в получении панорамного изображения сетчатки [8].

В настоящее время установлено, что совокупность таких факторов, как площадь аваскулярной зоны сетчатки, расширение и извитость ретинальных сосудов, а также величина угла между магистральными сосудами височной аркады, является прогностическими факторами в развитии РН.

Объективно определить степень этих нарушений для оценки течения заболевания и выбора оптимальной тактики ведения и лечения пациента возможно при проведении их достоверной количественной оценки, основанной на комплексном анализе цифровых изображений глазного дна.

До сих пор все попытки количественно определить параметры сетчатки при РН ограничивались подсчетами калибра и извитости сосудов на двухмерной фотографии заднего полюса глаза [2, 9, 10].

Цель работы - разработка программного обеспечения количественной оценки состояния сетчатки и ретинальных сосудов при ретинопатии недоношенных на базе трехмерной модели глазного яблока.

Материалы и методы

В основу программы положен метод проецирования цифровых фотографий глазного дна, полученных с помощью ретинальной камеры «RetCam120» с линзой 130°, на сферическую поверхность модели глаза.

Все измерения (площади заданных пользователем замкнутых контуров, длины прямолинейного и извитости криволинейного отрезка, угла между прямолинейными отрезками), необходимые для определения диагностических признаков (площади аваскулярной зоны сетчатки, диаметра и извитости ретинальных сосудов, величины угла между магистральными сосудами височной аркады), базируются на использовании криволинейных отрезков на сфере, соединяющих некоторые лежащие на ней точки A и B по кратчайшему пути. Такой отрезок всегда является дугой большого круга [4], причём его центр и радиус всегда такие же, как у сферы, а концы дуги совпадают с точками A и B. Кроме того, эта дуга полностью принадлежит плоскости AOB, проходящей через центр и концы дуги.

Исходя из указанных принципов, мы разработали математическую модель поверхности сетчатки глаза, а также методики расчета длин отрезков, углов между отрезками, извитости кривой и площадей участков произвольной формы на сфере, которые затем были реализованы в компьютерной программе.

математическая модель

поверхности сетчатки глаза

Процесс построения графической модели сетчатки сопряжен с определёнными трудностями. Известно, что не существует способа получения развертки сферы на плоскость. Например, в картографии [5, 7] для создания карты мира (земного шара) применяют несколько способов построения проекций, при этом в одном из них сохраняются углы между произвольными отрезками (равноугольные проекции), в других сохраняются площади (равновеликие проекции).

После проведения предварительных исследований нами было принято решение моделировать трёхмерный рисунок поверхности сетчатки, «растянутый» на идеальной сфере. Такой подход позволяет использовать математический аппарат сферической тригонометрии [3].

Топология внутренней поверхности глазного яблока может быть чрезвычайно сложной, и смоделировать её, основываясь на одних лишь снимках глазного дна, невозможно. Кроме того, современная медицина не располагает неинвазивными способами оцифровки поверхности сетчатки, которые давали бы возможность получения рисунка сосудов. Тем не менее сферическая модель позволяет достичь необходимой точности расчётов, а в абсолютном рельефе сетчатки нет необходимости.

Нами был разработан и реализован метод трехмерного моделирования поверхности сетчатки, основанный на наложении двухмерных цифровых фотографий глазного дна и учитывающий свойства оптических сред глаза, который имеет следующие допущения:

• поверхность сетчатки образует идеальную сферу - предполагается, что пациент (грудной ребёнок) находится под наркозом, вследствие чего глазодвигательные мышцы расслаблены, и глазное яблоко стремится к сферической форме;

• фотографии, полученные с помощью ретинальной камеры <^еГСат-120», по ширине охватывают телесный угол 130° (величина угла заявлена производителем прибора) на поверхности сетчатки при таком расположении камеры, когда её ось совпадает с оптической осью глаза;

• деформации роговицы под действием усилия со стороны камеры незначительны;

• ответственность за корректное расположение и масштабирование фотографий глазного дна полностью лежит на пользователе программы;

• параметры оптической системы глаза заданы постоянными значениями и равны среднестатистическим.

При создании модели поверхности сетчатки производится проецирование фотографии на сферическую поверхность модели глаза (обратный процесс по отношению к происходящему во время съёмки фотографии). Прямоугольник (с центром в точке А), содержащий фотографию, располагается на расстоянии АВ от точки В (лежащей на сфере), обозначающей центр хрусталика. При этом прямоугольник всегда расположен перпендикулярно отрезку АВ. Луч света, исходящий из точки А, проходит в точку С. Аналогично луч света из некоторой точки D фотографии проходит по прямой через точку В в точку Е. Следует отметить, что радиус сферы R и расстояние АВ выбирают таким образом, чтобы обеспечить п. 2 перечисленных выше допущений.

Полученная таким образом проекция фотографии на сферу представляет собой рисунок сетчатки на сфере без учета оптических сил прозрачных сред гла-

за, что в некоторых случаях может обеспечить более высокую точность вычислений.

Улучшенная расчётная модель учитывает коэффициенты преломления роговицы, жидкости в передней камере, хрусталика и стекловидного тела, а также расположение границ этих сред, исходя из среднестатистических величин и некоторых параметров, задаваемых вручную (размеров передней камеры, переднезаднего отрезка, диаметра роговицы и т. п.).

метод расчёта площадей на сфере

Некоторый участок на сфере описывается с помощью ограничивающего его контура. Контур определяется массивом радиус-векторов, последовательность которых даёт непрерывную цепочку дуг, при этом замыкающая дуга проходит от последнего вектора массива к первому.

Лежащим внутри контура считается тот участок сферы, который располагается по правую сторону от дуг при движении вдоль них от первого вектора массива к последнему.

Заданный полигон подвергается дроблению на примитивы - эйлеровы треугольники [6], при этом используется алгоритм, основанный на условии видимости вершин. Так, для каждого полученного треугольника должно выполняться следующее условие: стороны треугольника (дуги) не должны пересекаться с дугами заданного контура. Ряд дополнительных условий позволяет провести корректную разбивку контура на примитивы в автоматическом режиме.

Площадь заданной области (S) определяется суммой площадей сферических треугольников, вычисляемых по формуле:

S = (а + ß + у — п) • R2

где а, ß, у - углы при вершинах треугольника,

R - радиус сферы.

метод расчёта длин отрезков

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

Отрезок задаётся двумя радиус-векторами - A и B. Соединяющая их дуга лежит в плоскости AOB, лежит на сфере и имеет её радиус R. Тогда длина дуги L выражается формулой:

-в- —> —>

L = а *R = A B*R = R* arccos

метод расчёта углов между отрезками

Для определения величины угла между некоторыми отрезками AB и BC указываются три радиус-вектора A, B и C, лежащих на поверхности сферы. Поскольку дуги AB и BC располагаются аналогично тому, как описано в методике расчета длин отрезков, то угол а между ними определяется как угол между нормалями к плоскостям, в которых они лежат:

ШІ

, = arccos (пАов • Пвос).

метод расчёта извитости кривой

Произвольная кривая определяется массивом радиус-векторов, описывающих некоторый интересующий пользователя сосуд. Коэффициент извитости (к) определяется отношением суммарной длины всех дуг,

последовательно соединяющих векторы массива, к длине дуги, проходящей от первого до последнего вектора:

т У R•arccos

к = —=^

L

R •arccos

Результаты и обсуждение

Разработанный нами математический аппарат был реализован в компьютерной программе (рис. 1).

Программа содержит меню «Вставка изображений», «Параметры глаза», «Контуры», «Вычисления», «Морфометрия» и обеспечивает выполнение следующих функций: проецирование двухмерных цифровых изображений глазного дна на поверхность виртуальной трехмерной сферы с нормализацией размеров плоского изображения; выполнение заданных пользователем перемещений спроецированных на сферу изображений в произвольных направлениях с наложением; вычисление диаметра и извитости ретинальных сосудов, угла между сосудистыми аркадами, площади аваскулярной и васкуляризирован-ной сетчатки.

Проецирование двухмерных цифровых изображений глазного дна на поверхность виртуального трехмерного глазного яблока.

Для получения трехмерной картины глазного дна пациента с ретинопатией недоношенных выбирают не менее 5 цифровых изображений сетчатки: центральная область сетчатки, верхнее, нижнее, височное и носовое поле сетчатки.

Каждое изображение помещают в отдельную вкладку «Изображение п» в меню «Вставка изображений», после чего программа автоматически проецирует его на сферу.

Совмещение изображений сетчатки на сфере (рис. 2) осуществляют путем последовательного наложения соседних изображений друг на друга таким образом, чтобы две однозначно идентифицируемые точки, например, точки дихотомического деления сосудов, на одном изображении совпали с этими же точками на соседнем изображении. Фотографии свободно перемещаются по поверхности сферы в произвольных направлениях с возможностью вращения и масштабирования.

Совместив таким образом все фотографии сетчатки, получают ее трехмерное изображение (рис. 3), которое можно вращать, обеспечивая просмотр в любой проекции.

Расчет площади сетчатки и ее зон

Для диагностики и прогнозирования течения РН особый интерес представляет площадь аваскулярной сетчатки. Разработанная нами программа автоматически рассчитывает данную величину. Для этого вначале в меню «Параметры глаза» вводят размеры передней камеры, переднезаднего отрезка, диаметра роговицы. Затем в меню «Контуры» обозначают границы аваскулярной зоны сетчатки замкнутой ломаной линией. После чего в меню «Вычисления» по команде «Вычислить» в соответствующем окне в абсолютном и относительном выражении отображаются значения площади аваскулярной и васкуляризированной зон сетчатки, а также общая площадь сетчатки глаза пациента с РН (рис. 4).

а

Кубанский научный медицинский вестник № 1 (124) 2011

Кубанский научный медицинский вестник № 1 (124) 2011

Рис. 1. Интерфейс программы для количественной оценки состояния сетчатки и ретинальных сосудов при РН

г  d  П Т П Ш

Рис. 2. Совмещение в программе четырех изображений сетчатки

Расчет диаметра и извитости сосудов, угла между магистральными сосудами височной аркады.

Расчет данных величин осуществляют на фотографии центральной области сетчатки, спроецированной на сферу. Для этого в меню «Морфо-метрия» при помощи команд «Отобразить линию», «Отобразить кривую», «Отобразить угол» интересующие диаметр сосуда, его извитость, угол между магистральными сосудами височной аркады обозначают отрезком, кривой и двумя лучами соответственно. Затем по команде «Вычислить» программа выводит на экран в окне вычислений искомые значения.

Заключение

Разработанная нами компьютерная программа количественной оценки состояния сетчатки и ретинальных сосудов при РН позволяет проводить индивидуальный расчет совокупности объективных и информативных количественных показателей: площади аваскулярной и васкуляризированной сетчатки, диаметра сосудов, их извитости, угла между магистральными сосудами височной аркады.

Рис. 3. Трехмерное изображение сетчатки пациента с РН с обозначенной ломаной линией границей аваскулярной зоны сетчатки

Рис. 4. Результаты вычислений площадей аваскулярной и васкуляризированной зон сетчатки и общей площади сетчатки глаза пациента с РН

Получаемые с помощью программы данные необходимы для достоверного прогноза течения активных стадий РН и оценки эффективности проводимого лечения.

литература

1. Harrison S. A., McCannell C. A., Brill E. R., Schwarts S. Retinopathy of prematurity using the 120 digital fundus camera // Invest. ophthalmol. vis. sci. - 1998. - Vol. 39. № 4. - Р. 591.

2. Heneghan C., Flynn J., O'Keefe M., Cahill M. Characterization of changes in blood vessel width and tortuosity in retinopathy of prematurity using image analysis // Med. image anal. - 2002. - Vol. 6. № 4. - P. 407-429.

3. Lorenz B., Bock M., Muller H., Massie N. Telemedecine based screening of infants at risk of retinopathy of prematurity // Stad. heals technol. inform. - 1999. - Vol. 64. - Р. 155-163.

4. Wallace D. K., Zhao Z., Freedman S. F. A pilot study using «ROPtool» to quantify plus disease in retinopathy of prematurity // J. AAPOS. - 2007. - Vol. 11. № 4. - P. 381-387.

5. Wallace D. K. Computer-assisted quantification of vascular tortuosity in retinopathy of prematurity // Trans am. ophthalmol. soc. -2007. - Vol. 105. - P. 594-615.

Поступила 25.09.2010