Научная статья на тему 'Компьютерная обработка папиллярных изображений при диагностике сахарного диабета'

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

CC BY
95
47
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Компьютерная обработка папиллярных изображений при диагностике сахарного диабета»

циальностеи в вузе; п^^ - количество курсов (потоков) на 1-й специальности; Ц^уи - количество студентов на ]-м потоке 1-й специальности, которые должны использовать 1-е ЛО в к-м семестре в течение часов; пло,норм - нормативное количество студентов, одновременно использующих единицу ЛО (пЛо,норм >1); Тлаб,ы - допустимый (потенциальный) объем лабораторных занятий с использованием 1-го ЛО в к-м семестре по вузу:

Алаб,к| К*~нед,к ^^факт,| ^ ^Ir , (2)

где кнед,к - количество недель в k-м семестре; Кфакгг,| - фактическое количество единиц l-го ЛО; кдн/нед ,к — количество учебных дней в неделю; tir — количество часов работы лаборатории с l-м оборудованием в r-й день.

Если неравенство (1) не выполняется, требуемое количество единиц лабораторного оборудования l-го типа для k-го семестра рассматриваемого (планируемого) учебного года можно оценить по формуле

K и £

план,к1

T

ло,к1

T

(3)

Если требуемое количество единиц К план,1> >К факг,ь необходимо приобрести или изготовить Лкло,и единиц ЛО 1-го типа.

Если выполняется условие К план,к1<К факт,и, ЛО 1-го типа достаточно для выполнения учебного процесса; при этом, если Кло план,к1<Кло факт,к1, следует перейти к процедуре выработки решения по использованию излишних единиц оборудования.

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

- применение возможно, но требуется переход на новое ЛО;

- оборудование соответствует современным требованиям.

При оценке предложений по приобретению нового ЛО могут использоваться следующие формулировки:

- не приобретать;

- приобретение возможно, допустимо;

- приобретать как современное и перспективное.

Оценки экспертов можно оставить в лингвистической форме или использовать дополнительные оценки в баллах. Модель данных должна обеспечивать хранение индивидуальных и согласованных (принятых) экспертных оценок (рис. 3). Процедура формирования согласованных оценок зависит от принятых для реализации методов проведения экспертиз [3].

Разработанные функциональные и инфологи-ческие модели процесса управления МТО учебного процесса являются основой для создания соответствующего информационного, программного и методического обеспечения ИАС вуза, что, в свою очередь, позволит выполнять мониторинг состояния МТО, повысить объективность и оперативность принятия управленческих решений по формированию образовательных ресурсов надлежащего качества.

Литература

1. Соломенцев Ю.М., Позднеев Б.М., Солдатов А.В. Эффективное управление ресурсами вуза. URL: http://tqm.stan-kin.ru/articles/2003/18/index.htm (дата обращения: 15.11.2010).

2. Васин Д.А., Пикулин В.В. Информационная система планирования ресурсов образовательных программ // Программные продукты и системы. 2011. № 1. С. 114-118.

3. Мартемьянов Ю.Ф., Лазарева Т.Я. Экспертные методы принятия решений: учеб. пособие. Тамбов: Изд-во Тамб. гос. технич. ун-та, 2010. 80 с.

k

УДК 004:616-07

КОМПЬЮТЕРНАЯ ОБРАБОТКА ПАПИЛЛЯРНЫХ ИЗОБРАЖЕНИЙ ПРИ ДИАГНОСТИКЕ САХАРНОГО ДИАБЕТА

А.А. Азази; Г.А. Дмитриев, д.т.н.; Н.А. Семенов, д.т.н.

(Тверской государственный технический университет, dmitrievtstu@maHru)

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

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

Широко распространенные заболевания (такие как сахарный диабет, гипертоническая болезнь,

коронарная болезнь, бронхиальная астма, шизофрения и др.) являются одновременным проявлени-

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

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

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

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

Процесс дерматоглифического анализа включает три основные стадии обработки дактилокар-ты:

1) сканирование и обработка изображения;

2) анализ и распознавание типа изображения с индивидуальными особенностями (выявление дерматоглифических признаков);

3) постановка диагноза на основании выявленных дерматоглифических признаков.

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

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

— дефекты, обусловленные особенностями папиллярных гребней гребешковой кожи человека (белые линии, дисплазия папиллярных гребней); низкое качество дерматоглифических изображений, обусловленное этими дефектами, при-

водит к значительным ошибкам при определении классификационных признаков.

Для улучшения изображений используются методы обработки в пространственной области (пространственные методы) и методы обработки в частотной области (частотные методы). Пространственные методы основаны на прямом манипулировании элементами изображения. Процессы пространственной обработки описываются уравнением g(x, у)=Т[1"(х, у)], где 1"(х, у) - входное изображение; g(x, у) - обработанное изображение; Т - оператор над 1"(х,у), определенный в некоторой окрестности точки (х, у).

Зашумленное изображение 1"(х, у) можно представить как сумму двух компонент: 1"(х, y)=g(x, у)+И(х, у), где И(х, у) - аддитивный шум. В этом случае g(x, у) выступает некоторым идеальным изображением. Главную роль при улучшении изображений и сегментации играет информация, извлекаемая из гистограмм яркости. Типичная гистограмма пальцевого узора приведена на рисунке 1.

Рис. 1. Исходное изображение и гистограмма

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

Частотные методы базируются на двухмерном дискретном преобразовании Фурье: =

j2n (ux/M+ vy/N)

при u=0, 1, 2, ..., M-1 и

x=0 y=0

у=0, 1, 2, ..., N-1. Преобразование Фурье обеспечивает значительную гибкость при разработке и реализации алгоритмов фильтрации при решении задач улучшения изображений. Визуальный анализ спектра существенно помогает при использовании инструментов анализа изображений в частотной области. На рисунке 2 приведен спектр Фурье для исходного изображения, приведенного на рисунке 1.

Кроме того, улучшения изображений можно добиться с помощью его морфологической обра-

M1 N1

Рис. 2. Реализация случайной функции и спектр Фурье

ботки, используя операции замыкания и размыкания.

Следующим этапом обработки является бинаризация изображения, сущность которой заключается в разбиении полутоновых изображений на две области, одна из которых содержит все пиксели со значением ниже некоторого порога, а другая - выше этого порога. Для этого используется гистограмма яркости пикселей вида ^(г^п^ где г - з-й уровень яркости из интервала [1, S]; пз - число пикселей изображения, уровень яркости которых равен гз. Поскольку общее число пикселей изображения равно п, вероятность появления уровня интенсивности г в данном изображении определяется выражением рз=^(гз)/п.

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

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

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

Матрица Сг задается выражением 1 к

с* К-™.^-^ , где 2к - вектор-

к -1 к=1

столбец, составленный из яркостей одноименных

пикселей g

родственных

изображений. При

K=MN: m = — 2 z

z XT ^^

K k=1

Преобразование главных компонент определяется формулой y=A(z-mz).

Строками матрицы A служат нормированные собственные векторы матрицы Ck. Это преобразование необходимо для приближения исходного

л A T

вектора z = АтУ- mz, когда используется лишь q

собственных векторов. Ошибка приближенного

восстановления вектора составляет величину g q

а = 2^ -, где Xj - собственные числа мат-j=i j=i

рицы Cz.

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

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

Эффективным методом решения данной задачи является согласованная фильтрация, основанная на корреляционном сопоставлении. Корреляционное сопоставление заключается в нахождении позиций на изображении f(x, y), которые лучше всего соответствуют заданному эталону w(x, y). Для этого используется корреляционная метрика - нормированный коэффициент корреля-

ттт™. „ 2 ((f(x,y) - f")(w(x,y) - w°)

ции: f = . I = ,

42((f(x,y) -f°)W2(w(x,y) - w°)

где f°, w° - средние значения интенсивности для исходного изображения и маски соответственно.

Если дано n эталонных изображений {wj}, i=1, ..., n, каждое из которых соответствует i-му классу, то обнаружение фрагмента изображения f(x, y) основано на решении оптимизационной задачи вида j = argmax(Rf,w).

Это есть задача поиска максимального значения нормированного коэффициента корреляции на множестве {wi}.

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

f(x,y)°w(x,y)»F(u,v)W(u,v). где f(x, у) о w(x, y) - пространственная корреляция. Таким образом, пространственную корреляцию можно получить с помощью обратного пре-

Оцифрованное изображение

k J J

Д 1 1>

Изображения

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

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

К) =-оЛСГ^-

справедливо для любого распределения случайной величины

Байесовская процедура классификации состоит в отнесении вектора наблюдений £ к если Р№ | С)>Р(^21 С), и к если Р(^11 С)<Р(^21 С), где qi - априорная вероятность того, что индивидуум принадлежит к популяции i=1, 2 (предполагается, что сумма априорных вероятностей q1+q2 равна 1); Р(С№) - условная вероятность получения некоторого вектора наблюдений если известно, что объект принадлежит к популяции

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Величины и Р(¥2|£) являются апосте-

риорными вероятностями.

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

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

- определение качества изображения и его корректировка;

- бинаризация изображения;

- ориентация изображения и его нормализация;

- морфологический анализ и выделение отдельных элементов изображения (ладонных линий);

- истончение (скелетизация) изображения для определения количественных показателей, таких как гребневый счет;

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

Модуль предварительной обработки

I

Генератор эталонных изображений

d

БД изо- БД

бражений эталонов

<Z. Z5 d

Модуль распознавания изображений

Генератор дескрипторов

Классификатор

Рис. 3. Структура информационной системы обработки дерматоглифических изображений и ранней диагностики сахарного диабета

Последовательность обработки изображений определяется сценарием дерматоглифического анализа и управляется программой-супервизором, составляющей ядро информационной системы. Модули обработки изображения, формирования дескрипторов и диагностики подключаются последовательно по мере выполнения алгоритма анализа. Кроме того, система содержит вспомогательные модули, обеспечивающие управление выполнением заданий, описание сценариев обработки, визуализацию данных, модуль, обрабатывающий некоторые возникающие ошибки. Этап обработки изображений заканчивается генерацией классификационных признаков, на их основе делается вывод о предрасположенности пациента к заболеванию сахарным диабетом. Значения классификационных признаков передаются в модуль-классификатор для дальнейшей обработки и формирования заключения. БД системы содержит сведения о пациентах и цифровые изображения отпечатков ладоней правой и левой руки, а также цифровые изображения эталонов (масок) для поиска необходимых фрагментов. Система реализована в среде MatLab с использованием пакетов расширения Image Processing Toolbox и Database Toolbox. Для классификации в данном случае использовался пакет статистического анализа Statistica. В дальнейшем для этой цели будет использоваться MatLab. Система предназначена для специалистов, врачей, проводящих дерматогли-фические исследования, которые являются необходимой частью клинического осмотра в генетической клинике. В настоящее время система проходит апробацию на реальных данных.

Литература

1. Рогачева Е.А. Особенности дерматоглифики у больных сахарным диабетом 1 -го типа // Успехи современного естествознания. 2009. № 6. С. 18-19.

2. Якушева М.Ю. [и др.]. Выявление генетической предрасположенности к возникновению сахарного диабета 2-го типа // Эфферентная терапия. 2007. Т. 13. N° 4. С. 9-13.

3. Гонсалес Р., Вудс Р., Эддинс С. Цифровая обработка изображений в среде MatLab. М.: Техносфера, 2006. 616 с.

4. Визельтер Ю.В. [и др.]. Обработка и анализ цифровых изображений. М.: ДМК Пресс, 2007. 464 с.

i Надоели баннеры? Вы всегда можете отключить рекламу.