Изобретение относится к медицине, а именно к хирургии, и касается диагностики заболеваний вен нижних конечностей с применением математической экспертной системы для обработки данных комбинированной термометрии.
Целью изобретения является разработка способа диагностики заболеваний вен нижних конечностей, который позволит классифицировать заболевания вен нижних конечностей по данным комбинированной термометрии и провести численную верификацию диагноза.
Наиболее близок к заявляемому нами и служит основным прототипом способ диагностики рака молочной железы по данным микроволновой радиометрии с применением математической экспертной системы. Авторы - Вайсблат А.В. и Бурдина Л.М. (1999), проанализировав большое количество термограмм больных раком молочной железы, выделили 6 характерных признаков этого заболевания, которые были формализованы:
1. Повышенное максимальное значение температуры соска по сравнению со средней температурой молочной железы.
2. Повышенная разница температур сосков правой и левой молочной железы.
3. Повышенное максимальное значение разности температур между одноименными точками правой и левой молочной железы.
4. Повышенное среднеквадратичное значение разностей температур между одноименными точками правой и левой молочной железы.
5. Повышенное среднеквадратичное значение разброса температур в одной из молочных желез вычисляется отдельно по точкам для каждой из молочных желез, рассматривается большее из двух значений.
Принцип компьютерной диагностики рака молочной железы состоит в следующем: термограммы всех верифицированных случаев рака вводятся в память компьютера. Далее в процессе диагностики программа сравнивает, насколько термограмма обследуемой пациентки, по выделенным признакам, близка к термограммам пациенток с установленным диагнозом рака молочной железы.
Однако анатомические и функциональные особенности молочной железы не позволяют применить эту экспертную систему для способа диагностики патологии вен нижних конечностей.
Поставленная нами цель достигается путем проведения последовательного измерения кожной (инфракрасной) температуры и глубокой микроволновой термометрии (РТМ) в 12 симметричных точках, расположенных по задней поверхности обеих голеней пациента. Для этого мы используем диагностический комплекс РТМ-01-РЭС, позволяющий регистрировать как инфракрасную температуру с поверхности тела, так и внутреннюю температуру тканей на глубине от 4 до 6 см по их естественному тепловому излучению в микроволновом диапазоне. В состав комплекса входят: антенна-аппликатор для регистрации микроволнового излучения, датчик регистрации инфракрасной температуры, блок обработки информации, персональный компьютер (ПЭВМ). Антенна-аппликатор устанавливается на соответствующую точку на поверхности голени плотно, без зазора между кожей и плоскостью аппликатора, без дополнительного нажима. Температура измеряется последовательно попеременно в 12-ти точках по задней поверхности правой и левой голени. Первые три точки расположены: на вершине латеральной головки икроножной мышцы (1), в подколенной ямке (2), на вершине медиальной головки икроножной мышцы (3). Второй ряд точек расположен: по центру латеральной головки икроножной мышцы (4), между головками икроножной мышцы (5) и по центру медиальной головки икроножной мышцы (6). Третий ряд точек расположен в нижней части икроножной мышцы - латерально (7), в центре (8) и медиально (9). Последние точки измерения температуры находятся на наружной стороне ахилового сухожилия в области латеральной лодыжки (10), по центру ахилового сухожилия (11) и по его внутренней поверхности в области медиальной лодыжки (12). Измерение температуры в указанных точках проводится последовательно радиодатчиком (РТМ) и датчиком температуры кожи (ИК) в положении пациента «лежа на животе» и «стоя». Результаты измерений поступают в виде постоянного напряжения на блок обработки и далее на ПЭВМ, где осуществляется фиксация данных измерений температуры, привязанных к позиции датчиков. Данные по температуре обрабатываются и могут быть отображены на мониторе или принтере в виде термограммы или в виде цветового поля температур. Все обследуемые больные имеют предварительный диагноз, поставленный на основании ультразвукового доплеровского обследования на ультразвуковом доплерографе БИОСС-1 датчиком 8 МГц с помощью проб Вальсальвы и манжеточной пробы (Market A., Meissner M., Manzo a., et al., 1994), а также ультразвукового сканирования состояния подкожных и глубоких вен в положении стоя на ультразвуковом сканере "Лоджик а-100" фирмы «Дженерал электрик».
Проведенные нами исследования здоровых конечностей у лиц различного возраста выявили определенный диапазон и закономерность распределения температур на голени здорового человека (Табл. 1, 2).
Из представленных данных видно, что для голени здорового человека характерен определенный термографический профиль, который нарушается при появлении у него острого венозного тромбоза вен (Табл. 3, 4), посттромботической болезни (Табл. 5, 6) или варикозной болезни (Табл. 7, 8). При статистической обработке данных термометрии учитывались изменения вен каждой ноги каждого пациента. Т.е., если пациент имеет правую здоровую ногу, а левую - с варикозной болезнью, то правая нога обсчитывалась в группе здоровых конечностей, а левая - в группе больных варикозной болезнью.
Для математической верификации диагноза термометрии мы представили наши исходные данные как упорядоченный набор из 48 признаков-координат, каждый из которых представляет собой температуру, измеренную в определенной точке правой и левой голени пациента в положении лежа и стоя микроволновым и инфракрасным датчиком.
В наших исследованиях мы использовали метод классификации пациентов с болезнями вен нижних конечностей на основе непараметрического оценивания плотности распределения вероятности с помощью многомерного ядра Пуассона и два алгоритма дискриминантного анализа, основанные на отношении плотностей распределения двух классов, оценки которых построены методами, описанными ниже. Наши исходные данные мы представили как упорядоченный набор из р признаков-координат (р=48), т.е. в виде матрицы:
где уi выступает в качестве номера болезни 1-го наблюдения.
Обозначим X=(x1, …, xР) - новое наблюдение, которое необходимо отнести к одному из имеющихся классов и причислить номер класса уi.
За Hj тогда обозначим соответствующую гипотезу о принадлежности нового наблюдения к Xj-му классу.
Fj(…)=F(…|Hj) - функция распределения вектора X, принадлежащего j-му классу.
fj(…)=fj(…|Hj) - плотности распределения вероятностей вектора X, принадлежащего j-му классу, j=1, 4.
Основная задача дискриминантного анализа - задача отнесения нового наблюдения к одному из двух классов (частный случай нашей задачи). Она тесно связана с классической задачей проверки простой гипотезы Н1 против альтернативной гипотезы Н2
H1 : X ∈ F1
Если пространство всех возможных исходов Х обозначить как Ω, то чтобы проверить какую-либо гипотезу H1 исходя из выборки наблюдений, мы должны разделить выборочное пространство Ω на две области. Если наблюдаемая выборочная точка Х попадает в одну из этих областей, скажем ω, то гипотеза H1 отвергается, если же Х попадает в дополнительную область Ω/ω, то гипотеза принимается. Область ω называется критической областью критерия, a Ω/ω - областью принятия гипотезы. Если известно распределение вероятностей наблюдений, соответствующее проверяемой гипотезе H1, то можно определить ω так, чтобы при выполнении гипотезы H1 вероятность отвергнуть эту гипотезу (т.е. вероятность попадания Х в ω) была равна заранее заданной величине α, т.е.
P{X∈ω|H1}=α
Число α называют уровнем значимости критерия. Ошибки, возникающие при проверке статистической гипотезы, могут быть двух типов:
1) можно ошибочно отвергнуть гипотезу, когда она верна;
2) можно ее ошибочно принять, когда она неверна.
Эти ошибки называются соответственно ошибками первого и второго рода. Вероятность ошибки первого рода равна уровню значимости. Вероятность ошибки второго рода, обозначаемая β, зависит от рассматриваемой альтернативной гипотезы Н2. Таким образом,
Р{Х∈(Ω/ω)|H2}=β или P{X∈ω|H2}=1-β.
Вероятность (1-β) называется мощностью критерия для проверки гипотезы против альтернативной гипотезы Н2.
Функция правдоподобия выборки x1, x2, …, определяются как
Среди всех критериев, различающих две альтернативные гипотезы H1 и Н2 с уровнем значимости а, наиболее мощным является критерий отношения правдоподобия, основанный на статистике
где L - функция правдоподобия. При γ(X)≤Cα, принимается гипотеза Н1, при γ(X)>Cα - соответственно Н2.
Введя априорные вероятности гипотез πj, которые вычисляются как отношение числа наблюдений Nj из j-го класса к общему числу наблюдений N, рассмотрим частный случай описанного критерия.
Пусть
Если γ(Х)>1, принимается гипотеза Н2, если γ(Х)≤1, соответственно принимается гипотеза H1.
Формула 2 представляет собой байесовский классификатор, минимизирующий вероятность принятия неверного решения по поводу принадлежности нового наблюдения к той или иной болезни.
Известно, что большинство используемых на практике алгоритмов классификации строится исходя из формулы 1. При этом возможно оценивание неизвестных параметров θ предполагаемых теоретических распределений и замена их оценками θ в плотности распределения. Далее вычисляется оценка отношения :
Это параметрический метод построения алгоритмов классификации.
Если нельзя сделать предположение о структуре многомерной плотности распределения вероятности, то для данной точки сразу, минуя оценку параметров, строится оценка отношения
Это непараметрические методы. В нашем изобретении мы использовали непараметрический метод оценивания плотности распределения, основанный на приближении δ-образной последовательностью многомерной δ-функции Дирака.
Тогда оценка плотности для j-класса представима в виде:
Где - произвольные параметры сглаживания, s=1, … р.
Если различна мера разброса р признаков-координат, то предлагается в качестве параметров сглаживания bj, j=1, … 4 рассматривать минимальное, среднее или максимальное расстояние , … от s-й координаты нового наблюдения
Х=(x1, …xs, …xp) до s-й координаты наблюдений , , … из j-го класса:
Здесь Nj - число наблюдений в j-м классе,
- наблюдения из базы i=1, …Nj для j-й болезни.
На практике подбор подходящего параметра сглаживания bj производится путем численного эксперимента.
Проверка проводилась двумя алгоритмами классификации: на основе байесовского классификатора и с использованием апостериорных вероятностей.
Каждого человека из базы рассматривали как нового больного и классифицировали с помощью экспертной системы относительно остальных наблюдений. Далее сравнивали предварительный диагноз (по данным доплеровского исследования и комбинированной термометрии) с поставленным посредством разработанного способа классификации.
При проведении численных экспериментов были использованы два способа задания параметров сглаживания : один, описанный формулой (5), другой представлял собой случай, когда , ∀j, ∀s.
Пример 1.
Первый алгоритм основан на байесовском классификаторе и на методе объединения больных в общий класс с последующим поэтапным вычленением каждой болезни.
а. На первом шаге, вычисляя для здоровых оценку плотности вероятности fi(X) и для объединенного класса всех больных находим оценку отношения
Обозначим априорные вероятности для класса здоровых π1, для класса больных -
π234.
Если ,
то принимается гипотеза Н234, говорящая о том, что новый пациент принадлежит классу больных. В ином случае выполнена гипотеза Н1 - человек из класса здоровых.
b. В случае, если новое наблюдение на первом шаге не было классифицировано как здоровый, вычисляем для класса ВБ оценку плотности вероятности , априорную вероятность π2, для объединенного класса двух других болезней (ОВТ и ПТБ) общую оценку плотности распределения и априорную вероятность π34. Снова находим оценку: .
Если , то принимается гипотеза
Н34, говорящая о том, что новый пациент принадлежит классу ОВТ и ПТБ-больных. В ином случае выполнена гипотеза Н2 - человек из класса ВБ.
с. Далее, если вновь пришедший больной не классифицирован как здоровый или ВБ-больной, то вычисляем для класса ПТБ-больных оценку плотности распределения , априорную вероятность
π3 и сравниваем их с оценкой плотности распределения класса ОВ-больных , вычислив априорную вероятность ОВТ-класса π4:
При γ(Х)≥π3/π4 принимается гипотеза Н4 (болен ОВТ), в ином случае - гипотеза Н3 (болен ПТБ).
Пример 2.
Второй алгоритм базируется на вычислении апостериорных вероятностей.
Одним из основных определений, которое будет использоваться в дальнейшем, является понятие условной (апостериорной) вероятности отнесения наблюдения к данному классу при заданном значении вектора признаков X. Используя формулу Байеса, запишем условную вероятность в виде:
В ситуации с полной информацией решающее правило, минимизирующее условную вероятность ошибки, может быть определено в виде: объект относится к классу H1, если выполняется неравенство
P[H1|x]>P[H2|x],
и к классу Н2 в противном случае.
Обобщая на случай четырех болезней можно сказать, что тот класс, которому соответствует наибольшее значение апостериорной вероятности P[Hj|X], и будет искомым классом для наблюдения X.
Для каждого класса строится оценка
где P[Hj|x] - условная (апостериорная) вероятность появления j-го класса при заданном значении вектора признаков X, и полученные апостериорные вероятности затем сравниваются между собой.
Результаты численных экспериментов.
1. Результаты численных экспериментов для классификации "Здоров - Болен" сравнивались с базой данных диагнозов заболеваний вен нижних конечностей, полученных на предварительном обследовании пациентов: 58 человек больные, 19 человек здоровые.
1.1. Классификатор, основанный на апостериорных вероятностях со средними в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 25 человек.
Классифицированные как Больные - 52 человека.
Ошибка классификации - 13 человек.
1.2. Классификатор, основанный на апостериорных вероятностях с единицей в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 11 человек.
Классифицированные как Больные - 66 человек.
Ошибка классификации - 8 человек.
1.3. Байесовский классификатор со средними в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 31 человек.
Классифицированные как Больные - 46 человек.
Ошибка классификации - 16 человек.
1.4. Байесовский классификатор с единицей в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 16 человек.
Классифицированные как Больные - 61 человек.
Ошибка классификации - 3 человека.
Наименьшую ошибку классификации на два класса (3.89%) имеет методика анализа, построенная на байесовском классификаторе с параметром сглаживания b, равным единице.
Вычислим чувствительность, точность и специфичность разработанного метода с наименьшей ошибкой классификации.
По результатам данной методики, заболевания вен нижних конечностей были подтверждены у 58 человек (истинно положительный диагноз). Обозначим это число Nip.
Ни один человек с верифицированным диагнозом болезней вен нижних конечностей не был классифицирован как здоровый (ложноотрицательный диагноз). Обозначается это число NLQ.
У 16 человек (NIQ) с неварикозными заболеваниями их характер был подтвержден (истинноотрицательный диагноз).
У 3 человек (NLp) с неварикозными заболеваниями были найдены признаки варикоза (ложноположительный диагноз).
Обозначим общее число людей через N.
Чувствительность метода (выявляемость заболеваний вен), вычисляемая по формуле для разработанного метода составила 100%.
Точность метода, вычислим как
Для нашей методики она составила 96.1%.
Специфичность метода, вычисляемая как
в нашем случае равна 84.2%.
Приведем для сравнения чувствительность, точность и специфичность экспертной системы по дифференциальной диагностике рака молочной железы у женщин: чувствительность обследования составляет 89.7%; точность - 86.4%; специфичность метода - 81.8%.
2. Результаты решения задачи, когда необходимо поставить точный диагноз (более широкая задача), класс "Больные" разбиваются на три болезни: ВБ, ПТБ и ОВТ
Результаты также сравнивались с клиническими данными: 19 здоровых, 20 человек с диагнозом ВБ, 22 человека - ПТБ, 16 человек - ОВТ.
2.1. Классификатор, основанный на апостериорных вероятностях со средними в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 28 человек.
Классифицированные как ВБ - 18 человек.
Классифицированные как ПТБ - 20 человек.
Классифицированные как ОВТ - 11 человек.
Ошибка
классификации - 16 человек.
2.2. Классификатор, основанный на апостериорных вероятностях с единицей в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 14 человек.
Классифицированные как ВБ - 21 человек.
Классифицированные как ПТБ - 25 человек.
Классифицированные как ОВТ - 17 человек.
Ошибка
классификации - 5 человек.
2-3. Байесовский классификатор со средними в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 31 человек.
Классифицированные как ВБ - 16 человек.
Классифицированные как ПТБ - 0 человек.
Классифицированные как ОВТ - 30 человек.
Ошибка
классификации - 36 человек.
2.4. Байесовский классификатор с единицей в качестве параметра сглаживания b, общее число людей: 77.
Классифицированные как Здоровые - 16 человек.
Классифицированные как ВБ - 21 человек.
Классифицированные как ПТБ - 10 человек.
Классифицированные как ОВТ - 30 человек.
Ошибка классификации - 25 человек.
Наименьшую ошибку классификации на четыре класса (6.5%) имеет методика анализа, основанная на апостериорных вероятностях с единицей в качестве параметра сглаживания b.
Нами получены следующие результаты.
1. В случае, когда все признаки-координаты измеряются в одной температурной шкале (в нашем случае эта шкала температур в С°), наилучший результат классификации достигается, если в качестве параметров сглаживания выступают единицы (, ∀j, ∀s).
2. Предложенная методика, основанная на байесовском классификаторе, с 3.89% ошибкой классифицирует новое наблюдение к одному из двух классов: "Здоров" или "Болен". Чувствительность данной методики составляет 100%, точность и специфичность превосходит запатентованную и разрешенную к применению в медицинской практике экспертную систему по дифференциальной диагностике рака молочной железы у женщин на 8.7% и 2.4% соответственно.
3. Применение разработанных методик для более точной классификации с указанием конкретного диагноза пациента дает немного худшие результаты (ошибка классификации 6.4%). Экспертная система в этом случае базируется на вычислении апостериорных вероятностей, т.к байесовский классификатор выдает ошибку классификации 33.7%.
Таким образом, предлагаемый нами способ диагностики на основе байесовского классификатора и апостериорных вероятностей с единицей в качестве параметра сглаживания по данным комбинированной термометрии имеет высокую точность (максимальная ошибка классификации 6,4%) и может быть предложен в качестве базы для разработки интеллектуальной экспертной компьютерной системы диагностики заболеваний вен нижних конечностей.
Использованная литература
1. Кендалл М., Стьюарт А. "Статистические выводы и связи." М.: Наука. 1973 г.
2. Айвазян С.А., Бухштабер В.М., Епюков И.С. Мешалкин Л.Д. "Прикладная статистика: Классификация и снижение размерности: Справ, изд." Под ред. Айвазяна С.А. - М.: Финансы и статистика, 1989 г.
3. Колмогоров А.Н., Фомин СВ. "Элементы теории функций и функционального анализа." М.: Наука, 1972.
4. Деврой Л., Дьерфи Л. "Непараметрическое оценивание плотностей. L1-подход". М.: Мир, 1988 г.
5. Дуда Р., Хард П. "Распознавание образов и анализ сцен." Изд. "Мир", Москва, 1976 г.
6. Фукуиага К. "Введение в статистическую теорию распознавания образов". М.: Наука, 1979 г.
7. Бурдина Л.М., Хайленко В.А., Кижаев Е.В., Легков А.А., Пинхосевич Е.Г., Мустафип Ч.К., Вайсблат А.В., Веснин С.Г., Тихомирова Н.Н. "Применение радиотермометра диагностического компьютеризированного интегральной глубинной температуры ткани для диагностики рака молочной железы". - Пособие для врачей. - М.: 1999 г.
8. Вайсблат А.В., Веснин С.Г., Конкин М.А. «Использование микроволновой радиотермометрии в диагностике рака молочной железы». - www.resltd.ru
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у здоровых лиц в положении лежа
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у здоровых лиц в положении стоя
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у больных с ОВТ в положении лежа
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у больных с ОВТ в положении стоя
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у больных посттромботической болезнью в положении лежа
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у больных посттромботической болезнью в положении стоя
Среднее значение инфракрасного (ИК) и микроволнового (РТМ) датчиков температуры в различных точках голени у больных варикозной болезнью в положении лежа
Среднее значение инфракрасного (ИК) и микроволнового (MB) датчиков температуры в различных точках голени у больных варикозной болезнью в положении стоя
Изобретение относится к медицине и может быть использовано для верификации диагнозов болезней вен нижних конечностей по данным комбинированной термометрии на основе байесовского классификатора и апостериорных вероятностей. Для этого измеряют поверхностные и глубокие температуры тканей голени пациентов. Затем для классификации пациентов применяют математическую экспертную систему на основе байесовского классификатора и апостериорных вероятностей с единицей в качестве параметра сглаживания. На основе полученных данных диагностируют заболевания вен нижних конечностей. Заявленный способ позволяет достигнуть высокой точности верификации диагнозов болезней вен нижних конечностей и может быть использован в качестве базы для разработки интеллектуальной экспертной компьютерной системы диагностики заболеваний вен нижних конечностей. 8 табл.
Способ диагностики заболеваний вен нижних конечностей по изменению поверхностной и глубокой температуры тканей голени пациентов, отличающийся применением для классификации пациентов математической экспертной системы на основе байесовского классификатора и апостериорных вероятностей с единицей в качестве параметра сглаживания.
СПОСОБ ДИАГНОСТИКИ ТРОМБОЗА ГЛУБОКИХ ВЕН ГОЛЕНИ И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ | 2003 |
|
RU2262883C2 |
ПРИМЕНЕНИЕ ЛАВАНДОВОГО МАСЛА ДЛЯ ПРОФИЛАКТИКИ И ЛЕЧЕНИЯ НЕВРАСТЕНИИ, СОМАТОФОРМНОГО РАССТРОЙСТВА И ДРУГИХ БОЛЕЗНЕЙ, СВЯЗАННЫХ СО СТРЕССОМ | 2005 |
|
RU2406521C2 |
КОЗЛОВ В.Б | |||
Применение тепловизионного исследования для выявления недостаточности соединительных вен при варикозной болезни нижних конечностей | |||
// Клин | |||
хирургия, 1989, №7, с.73, 74 | |||
МАРКЕЛОВ С.И | |||
Диагностика заболеваний вен нижних конечностей с применением термографии жидкими |
Авторы
Даты
2009-06-20—Публикация
2006-12-19—Подача