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

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

Гайдель А.В.1, Переушкин С .С2 1 Самарский государственный аэрокосмический униеерситетимени академика С. П. Королёва (национальный исследовательский университет), 2 Самарский государственный медицинский университет

Аннотация

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

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

Введение

Остеопороз - это системное заболевание, поражающее все кости скелета, сопровождающееся снижением плотности и прочности костной ткани, а также нарушением её структуры, что приводит к высокому риску переломов даже при небольшой травме. Остеопения - это лёгкая форма остеопоро-за. Для предупреждения переломов важно как можно раньше поставить правильный диагноз, чтобы назначить лекарственные средства, замедляющие развитие заболевания [1, 2].

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

С другой стороны, имеется метод субъективной диагностики остеопороза по рентгеновским снимкам определённых областей интереса, на которых хорошо заметны изменения в так называемой трабе-кулярной структуре костной ткани [3]. Например, на рис. 1а и 1б видно, как меняется структура шейки бедра, поражённой остеопорозом: на втором изображении структура выглядит более прозрачной, волокнистой, редкой.

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

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

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

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

В зарубежных изданиях описаны некоторые попытки использовать отдельные группы признаков, например, маски Лоу [4], признаки Харалика [5], спектральные характеристики [6], вейвлет-ана-

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

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

Разрешение таких изображений до сканирования, по сути, зависит от разрешения рентгеночувстви-тельной плёнки. В данном исследовании используется плёнка KodakMedical X-Ray Film General Purpose Green (MXG). Её разрешающая способность составляет примерно 16 линий на мм или 400 точек на дюйм.

С нашей точки зрения, изображение на плёнке можно рассматривать как функцию ло (t1, t2), возвращающую яркость изображения в точке с координатами (t1, t2). Изображение ограничено, то есть 0 < t1 < l1 и 0 < t2 < l2. Оно монохромное, то есть имеет только один канал яркости - оттенки серого.

Полученные физические фотографии подвергаются сканированию, причём разрешение при этом не теряется. Можно считать, что в процессе оцифровки непрерывное изображение л0 (t1 , t2 ) последовательно подвергается двум преобразованиям: дискретизации и квантованию, в результате чего получается двумерное дискретное изображение л(m, n), состоящее из M х N отсчётов (0 < m < M и 0 < n < N). Каждый отсчёт характеризует уровень яркости и принимает целочисленные значения от 0 до L -1 = 255 включительно.

В порядке предварительной обработки изображения ориентируются, на них локализуется область интереса. От дополнительных преобразований, таких как контрастирование и фильтрация шумов, было решено отказаться, чтобы не потерять информацию, потенциально содержащуюся на изображении. Кроме того, на изображении x(m,n) специалистом выделяется полигональная область Dx, включающая только трабекулярную структуру.

2. Текстурные признаки

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

В данном исследовании были использованы следующие признаки на основе моментов функции яркости: средняя яркость х, средняя энергия (1), сред-неквадратическое отклонение о, коэффициент асимметрии, коэффициент эксцесса.

= D Z m n).

(i)

(m,n)e Dx

Корреляционные характеристики Часто в качестве текстурных признаков используют отсчёты выборочной нормированной корреляционной функции [10]:

X х (Я х (1 + т J + п)

Я(т,п) = (1'(т,п) ,-,-,

зШх (т, п)

где Ох(тп) = {(1,7')е Д I (1 + m,] + п)е Пх} .

Во многих источниках рекомендуется использовать именно такую корреляционную функцию без вычитания математического ожидания, так как колебания непосредственно функции яркости более информативны, чем её сдвиг относительно среднего значения, который может быть вызван, например, неравномерной освещённостью изображения [10, 11]. В ходе исследований использовались 8 отсчётов такой функции: для двух расстояний и четырёх направлений.

ПризнакиГабора Двумерная функция Габора определяется как

G (t1, t2; t0, t0, w0, ю0 ) = e

Mf+Mf

20]2 2g2

-i(W0ti+w0tz)

. (2)

Произвольный двумерный сигнал х0 (10) может

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

а (¿°, < ю02 ) =

+¥ +¥ (3)

= 11 х (1,*о)с(1,1°,ю0,ю0)а^.

Фурье-образ функции Габора (2) выражается в виде [10]:

^ 02 ) =

2 (w-w0 )2 +s2 (w-m°2 f j -;((ш1-ш0 +(w2-w0 )t0) '

G (w,, w2; t0, t0, ш0, ш2, о,

= 0яо1о0е

Используя Фурье-образ х0 (11,10) сигнала х0 (11,10), можно переписать разложение (3) в частотной области [10]:

а(^ t°,ш°,) =

= Л Л (t,, t2) G(t,, t2; t1, tl, w0, w0)dt,dt2

Это преобразование можно интерпретировать как поточечное умножение Фурье-образа сигнала на

s

гауссовское окно с последующим обратным преобразованием Фурье. Именно эта идея используется на практике.

Рассмотрим дополненный нулями дискретный сигнал

X (щ, щ ) =

Сх(п + п0,п2 + п0), (п1 + п0,п2 + п0)е Ц; (4) (О, (п + п0,п2 + п0)е Бх.

Здесь п° и п0 - левая и нижняя границы множества Ц соответственно. Функция X(п1, п2) определена

для п1 е{0,1,_, ММ -1} и п2 е {0,1,___, N-1},

где

М -1 и N -1 - правая и нижняя границы множества Ц соответственно.

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

- (т-т (к))2 (т-то (к))2 ^

2о2 (к) + 2о2 (к)

% (^ т ) =е -

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

Ук (т ) = х(т1, Щ ) ( %к (ml, т2 ) + %к (ml, т2 )) ,

где окно Шк (т1,т2) имеет центр (-т° (к),-т2° (к)),

чтобы полученный спектр всё ещё соответствовал вещественным изображениям. Признаками будем считать средние энергии этих изображений [10], аналогичные (1):

1 Ь 2 J Ь 2 J

ёк = X . X -Ук(^т)(ml,т2).

Разобьём частотную область на прямоугольные участки, как показано на схеме (рис. 2).

13 14 7 8 15 16

9 10 2 и 12

3 4 5 6

1 0 1

6 5 2 4 3

12 11 8 7 10 9

16 15 14 13

Центры гауссовских окон (т°(к), т^к)) поместим в центры соответствующих участков, а параметры о1(к) и о2(к) зададим как половины соответствующих линейных размеров этих участков. Такое разбиение позволяет лучше отличать высокочастотные текстуры [10].Таким образом, имеем 17 признаков Габора.

Признаки Харалика Рассмотрим р0(д 7) - вероятность того, что два отсчёта функции яркости, находящиеся на расстоянии ё в направлении 8 друг от друга, имеют цвета 1 и 7. Положим

(1, j ) =

= {(т, п)е Ц (с11, с12) |{х(т, п), х(т + с11, п + с12)} =

= {1,7}} ,

тогда эти вероятности могут быть оценены как

8 = 0;

Р«я( и 7 ) =

Ц (0, ё )| К (1,7 )1

Ц ( л , ё )|

Щ*,- л (1,7 )|

8 = —;

8 = +—;

8 = -—. 4

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

Ц (+Л,-Л )|

Здесь 1 и 7 могут принимать значения от 0 до Ь-1 включительно. Для фиксированных ё и 8 (далее эти индексы опущены) могут быть рассчитаны 14 признаков Харалика [12]: второй угловой момент, контраст Харалика, корреляция Харалика, дисперсия Харалика, обратный разностный момент, суммарное среднее, суммарная дисперсия, суммарная энтропия, энтропия, разностная дисперсия, разностная энтропия, первая информационная мера корреляции, вторая информационная мера корреляции, максимальный коэффициент корреляции.

В данном исследовании использовалось два расстояния и четыре направления, так что всего рассчитывалось 112 признаков Харалика. Дополнительные сведения об этих признаках и способах их быстрого вычисления можно найти в [13].

Признаки Тамуры

Имеются 6 признаков, признанных существенными для зрительного восприятия в результате психологического эксперимента [14]. Все они вычисляются по эвристическим процедурам и никак не масштабируются.

- Зернистость - это признак, связанный с расстоянием между заметными пространственными колебания-

2

4

М-1

N-1

+

+

т

2

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

- Контраст Тамуры - это мера того, насколько сильно и резко может меняться цвет на изображении.

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

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

- Регулярность - это общая изменчивость первых четырёх признаков между различными частями изображения.

- Грубость - субъективная оценка грубости переходов на изображении.

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

3. Экспериментальные исследования За время выполнения исследования в клиники Самарского государственного медицинского университета с подозрением на остеопороз поступило 50 пациентов. Каждому с помощью денситометрии был поставлен один из трёх диагнозов: остеопороз, остео-пения или без патологий. Кроме того, были получены рентгеновские снимки шейки бедра каждого пациента, они были оцифрованы и использовались как входные данные, наряду с известными диагнозами.

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

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

Пусть имеется обучающая выборка и е ККхУ, где и. - значение 1 -го признака для у -го вектора из обучающей выборки. Относительно каждого вектора-столбца признаков и7 е Ж.К известен его класс Ф (иу). Вычислим значения оценок внутриклассовых дисперсий

*)=V I и - к))2,

У1 у ф( и, )=1

где О)( к ) - дисперсия к -го признака внутри 1 -го

класса ( 1 принимает значения 0,1,2 для нормы, ос-теопении и остеопороза соответственно),

{уе {1,2,..., V} Ф(и7) = 1} - количество векторов класса 1 в обучающей выборке, М(1)(к) = — I и7 - среднее значение признака

^ у Ф(ц,)=1 к внутри 1 -го класса.

V =

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

1 2

к ) = 1 Ц ^)( к) .

3 1=0

(5)

Кроме того, определим значение оценки дисперсии смеси распределений:

(6)

IV %2

ОД = VIи - м<») ,

' 7=1

1 V

где М(к) = — 1ику - среднее значение к -го

V у=1

признака по всем элементам обучающей выборки.

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

О (к)

1 (к ) =

к)

(7)

после чего оставить только несколько первых.

В табл. 1 приведены десять лучших по критерию (7) признаков. Видно, что в данной задаче наиболее эффективными оказались корреляционные характеристики.

Таблица 1. Результаты дисперсионного анализа

№ Признак 1(к) е

1 Ближняя побочнодиагональная корреляция 3,71 0,32

2 Ближняя горизонтальная корреляция 3,60 0,52

3 Ближняя вертикальная корреляция 2,82 0,48

4 Дальняя побочнодиагональная корреляция 2,79 0,32

5 Ближняя главнодиагональная корреляция 2,50 0,56

6 Дальняя вертикальная корреляция 2,32 0,52

7 Дальняя горизонтальная корреляция 2,31 0,52

8 Ближний горизонтальный контраст Харалика 2,02 0,56

9 Дальний побочнодиагональный контраст Харалика 1,98 0,68

10 Дальняя побочнодиагональная разностная дисперсия 1,96 0,76

Рассмотрим контрольную выборку И е МКх"

объёма V. Настроенный с помощью обучающей выборки и классификатор реализует отображение Ф(к), которое каждому вектору признаков Оу ставит в соответствие класс, к которому его относит классификатор. Кроме того, относительно каждого вектора Оу известен его реальный класс Ф(11у).

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

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

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

1,005 1,000

0,995

0,990

0,985

0,980

* *• D

* * □ без патологий в остеопения * остеопороз

0,986 0,990 0,994 0,998 1,002 Ближняя побочнодиагональная корелляция

Рис. 3. Векторы из контрольной выборки

в пространстве из двух лучших признаков

Учитывая, что априорные вероятности появления объектов из различных классов неизвестны, имеем несмещённую оценку вероятности ошибочной классификации [17]:

1

£ = — V

{ jе {l,2,..., V} | Ф(Üj) (Üy)}

(8)

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

Для отбора лучшей группы признаков воспользуемся аналогом процедуры, описанной в [16]. Однако в отличие от [16] будем добавлять в текущую группу очередной признак в порядке убывания критерия (7) и оценивать ошибку классификации (8) векторов из контрольной выборки, произведённой с помощью данной группы.

Таблица 2. Результаты отбора группы признаков

Группа признаков e

1 0,32

1, 2 0,32

1, 2, 3 0,44

1, 2, 3, 4 0,36

1, 2, 3, 4, 5 0,40

1, 2, 3, 4, 5, 6 0,36

1, 2, 3, 4, 5, 6, 7 0,28

1, 2, 3, 4, 5, 6, 7, 8 0,20

1, 2, 3, 4, 5, 6, 7, 8, 9 0,24

1, 2, 3, 4, 5, 6, 7, 8, 9, 10 0,28

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

корреляционных и один признак Харалика. Оценка вероятности ошибки составляет 0,2.

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

Заключение

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

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

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

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

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

Литература

1. Руководство по остеопорозу / под ред. Л. И. Беневоленской. - М.: БИНОМ. Лаборатория знаний, 2003.

- 524 с.

2. Остеопороз / под ред. О. М. Лесняк, Л. И. Беневоленской. - М.: ГЭОТАР-Медиа, 2009. - 272 с.

3. Рейнберг, С. А. Рентгенодиагностика заболеваний костей и суставов / С.А. Рейнберг. - М.: Медицина, 1964.

- Т. 1. - 530 с.

4. Smyth, P.P . Application of computer texture analysis to the Singh Index / P.P. Smyth, J.E. Adams, R.W. White-house, C.J. Taylor // The British Journal of Radiology. -1997. - March, 70. - P. 242-247.

5. Lespessailles, E. Clinical interest of bone texture analysis in osteoporosis: a case control multicenter study / E. Lespessailles [et al.] // Osteoporosis International. -2008. - Vol. 19.7. - P. 1019-1028.

6. Brunet-Imbault, B. A new anisotropy index on trabecular bone radiographic imagesusing the fast Fourier transform /

B. Brunet-Imbault [et al.] // BMC Medical Imaging. -2005. - Vol. 5:4. - 11 p.

7. Mengko, T.R. Texture analysis of radiographs in the assessment of osteoporosis / T.R. Mengko, J.T. Pramudito // IAPR Workshop on Machine Vision Applications. - 2002.

- Dec. 11-13. - P. 184-187.

8. Хеннан, Э. Многомерные временные ряды / Э. Хеннан; пер. с англ. - М.: Мир, 1974. - 576 с.

9. Волков, И.К. Случайные процессы / И.К. Волков,

C.М. Зуев, Г.М. Цветкова. - М.: Изд-во МГТУ им. Н.Э. Баумана, 1999. - 448 с.

10. Petrou, M. Image Processing: Dealing with Texture / M. Petrou, P. Garcia Sevilla. - Chichester, UK: John Wiley & Sons, Ltd., 2006. - 618 p.

11. Марпл-мл., С.Л. Цифровой спектральный анализ и его приложения / С.Л. Марпл-мл.; пер. с англ. - М.: Мир, 1990. - 584 с.

12. Haralick, R.M. Textural features for image classification / R.M. Haralick, K. Shanmugam, Dinstein // IEEE Transactions on Systems, Man, and Cybernetics. - 1973. - November, Vol. SMC-3. - P. 610-621.

13. Gipp, M. Haralick's Texture Features Computed by GPUs for Biological Applications / M. Gipp, G. Marcus, N. Harder, A. Suratanee, K. Rohr, R. König, R. Männer // IAENG International Journal of Computer Science. - 2009. - February, Vol. IJCS_36_1_09. - 10 p.

14. Tamura, H. Textural Features Corresponding to Visual Perception / H. Tamura, Sh. Mori, T. Yamawaki // IEEE Transactions on Systems, Man, and Cybernetics. - 1978. -June, Vol. SMC-8. - P. 460-473.

15. Методы компьютерной обработки изображений / под ред. В.А. Сойфера. - М.: ФИЗМАТЛИТ, 2003. - 784 с.

16. Глумов, Н.И. Метод отбора информативных признаков на цифровых изображениях / Н.И. Глумов, Е.В. Мясников // Компьютерная оптика. - 2007. - Т. 31, № 3. -С. 73-76.

17. Фукунага, К. Введение в статистическую теорию распознавания образов / К. Фукунага; пер. с англ. - М.: Наука, 1979. - 368 с.

References

1. Osteoporosis Handbook / edited by L.I. Benevo-lenskaya. Moscow: "BINOM. Knowledge Laboratory" Publisher, 2003. - 524 p. - (In Russian).

2. Osteoporosis / edited by O.M. Lesnyak and L.I. Benevo-lenskaya. - Moscow: "GEOTAR-Media" Publisher, 2009.

- 272 p. - (In Russian).

3. Reinberg, S.A. X-ray diagnostics of bones and joints diseases / S.A. Reinberg. - Moscow: "Medicine" Publisher, 1964. - V. 1. - 530 p. - (In Russian).

4. Smyth, P.P . Application of computer texture analysis to the Singh Index / P.P. Smyth, J.E. Adams, R.W. White-house, C.J. Taylor // The British Journal of Radiology. -1997. - March, 70. - P. 242-247.

5. Lespessailles, E. Clinical interest of bone texture analysis in osteoporosis: a case control multicenter study / E. Lespessailles, [et all.] // Osteoporosis International. -2008. - Vol. 19.7. - P. 1019-1028.

6. Brunet-Imbault, B. A new anisotropy index on trabecular bone radiographic imagesusing the fast Fourier transform / B. Brunet-Imbault, [et all.] // BMC Medical Imaging. -2005. - Vol. 5:4. - 11 p.

7. Mengko, T.R.Texture analysis of radiographs in the assessment of osteoporosis / T.R. Mengko, J.T. Pramudito // IAPR Workshop on Machine Vision Applications. - 2002. - Dec. 11-13. - P. 184-187.

8. Hannan, E.J. Multipletimeseries. - N.Y.: John Wiley & Sons, Inc., 1970.

9. Volkov, I.K. Stochastic processes / I.K. Volkov, S.M. Zu-yev, G.M. Tsvetkova. - Moscow: "Bauman Moscow State Technical University" Publisher, 1999. - 448 p. - (In Russian).

10. Petrou, M. Image Processing: Dealing with Texture / M. Petrou, P. Garcia Sevilla. - Chichester, UK: John Wiley & Sons, Ltd., 2006. - 618 p.

11. Marple, S.L., Jr. Digital spectral analysis with applications / S.L. Marple, Jr. - Prentice-Hall, Inc. Englewood Cliffs, New Jersey, 1987.

12. Haralick, R.M. Textural features for image classification / R.M. Haralick, K. Shanmugam, Dinstein // IEEE Transactions on Systems, Man, and Cybernetics. - 1973. - November, Vol. SMC-3. - P. 610-621.

13. Gipp, M. Haralick's Texture Features Computed by GPUs for Biological Applications / M. Gipp, G. Marcus, N. Harder, A. Suratanee, K. Rohr, R. König, R. Männer // IAENG International Journal of Computer Science. - 2009. - February, Vol. IJCS_36_1_09. - 10 p.

14. Tamura, H. Textural Features Corresponding to Visual Perception / H. Tamura, Sh. Mori, T. Yamawaki // IEEE Transactions on Systems, Man, and Cybernetics. - 1978. -June, Vol. SMC-8. - P. 460-473.

15. Methods of the computer image processing / edited by V.A. Soifer. - Moscow: "Fizmatlit" Publisher, 2003. -784 p. - (In Russian).

16. Glumov, N.I. Method of the informative features selectio-non the digital images / N.I. Glumov, E.V. Myasnikov // Computer Optics. - 2007. - V. 31, N 3. - P. 73-76. - ISSN 0134-2452. - (In Russian).

17. Fukunaga, K. Introduction to statistical pattern recognition. - N.Y.: Academic Press, 1972.

RESEARCH OF THE TEXTURAL FEATURES FOR THE BONY TISSUE DISEASES DIAGNOSTICS

USING THE ROENTGENOGRAMS

A. V. Gaidel1, S.S. Pervushkin2

1S. P. Korolyov Samara State Aerospace University (National Reseach University),

2Samara State Medical University

Abstract

In the paper method of the musculoskeletal system osteoporotic diseases computer diagnostics with a help of the femoral neck roentgenograms texture an alysisis proposed and researched. Problem of the texture features providing acceptable recognition qualityselectionissolved. Method work is re-

searched on the real roentgenograms set with the known diagnoses. Researches results demonstrating a possibility to use proposed method int he clinical practice are adduced. During the work computer program designed both for proposed method research and for clinical diagnostics is developed.

Key words: image processing, texture analysis, pattern recognition, diagnostics, Gabor filter, Haralick'sfeatures, Tamura's features, features selection, analysis of variance.

Сведения об авторах Гайдель Андрей Викторович, 1989 года рождения, в 2012 году с отличием окончил Самарский государственный аэрокосмический университет имени академика С. П. Королёва по направлению «Прикладная математика и информатика». В настоящее время аспирант кафедры технической кибернетики СГАУ. Область научных интересов: компьютерная обработка изображений, теория распознавания образов, математическое моделирование, теория алгоритмов. E-mail: andrey.gaidel@gmail. com .

Andrey Viktorovich Gaidel(b. 1989) graduated with honors from Samara State Aerospace University in 2012, majoring in Applied Mathematics and Informatics. Currently he is a postgraduate student at Technical Cybernetics Department, SSAU. His research interests are currently focused on computer image processing, pattern recognition, mathematical modeling and theory of computation.

Первушкин Сергей Сергеевич, 1985 года рождения, в 2009 году окончил лечебный факультет Самарского государственного медицинского университета, в 2011 окончил ординатуру по специальности «Лучевая диагностика», аспирант кафедры лучевой диагностики и лучевой терапии с курсом медицинской информатики Самарского государственного медицинского университета. Область научных интересов: лучевая диагностика, ос-теопороз, компьютерная обработка медицинских изображений.

E-mail: sergey.pervushkin @gmail. com .

Sergey Sergeevich Pervushkin (b. 1985) graduated fromthe Department of General Medicine, Samara State Medical University in 2009 and from residency of Radiology in 2011, postgraduate student at Department of Radiology and Radiation Therapy with a Course of Medical tate Medical University. Research interests: radiology, osteoporosis, computer processing of

nocmynwa e peдaкцuю 1 дeкaбpя 2012 г.

mm

Informatics, Samara S medical images.