Метод оценки имплантационной динамики компонента "Литар" на основе сжатого частотного спектра рентгеновских снимков
Метод обработки изображений, представленных рентгеновским снимком, отражающим текущее состояние воздействия имплантированного компонента "ЛитАр" на поврежденную ткань человека. Оценка взаимодействия имплантируемого компонента с окружающей его тканью.
Рубрика | Медицина |
Вид | статья |
Язык | русский |
Дата добавления | 30.01.2020 |
Размер файла | 37,4 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Метод оценки имплантационной динамики компонента "Литар" на основе сжатого частотного спектра рентгеновских снимков
Л.С. Бекасов
В статье рассматривается метод обработки изображений, представленных рентгеновским снимком, отражающим текущее состояние воздействия имплантированного компонента «ЛитАр» на поврежденную ткань человека. Предлагаемый метод нелинейного спектрального преобразования позволяет получить сжатый образ изображения с высокой разрешающей способностью и тем самым количественно оценить взаимодействие имплантируемого компонента с окружающей его тканью. Метод не требует относительно больших затрат времени и характеризуется инвариантностью по отношению к сдвигу данных в пределах анализируемого окна. рентгеновский снимок имплантируемый ткань
Ключевые слова: рентгеновский снимок; имплантируемый материал; импульсная система базисных функций; сжатый спектр.
Введение
В последнее десятилетие вопрос применения имплантационных материалов, способных восстанавливать пациентам натуральную ткань в области дефекта, становится наиболее распространенным. Особый интерес представляют материалы, которые имеют высокую скорость биодеградации и универсальный характер действия, то есть могут восстанавливать биологические ткани разного типа. Материалом, чьи свойства наиболее приближены к желаемым, является имплантационный биодеградируемый полимер-солевой композит «ЛитАр», производство которого осуществляется в г. Самаре организацией, носящей то же имя [1-5].
Для отслеживания динамики процесса регенерации ткани используют рентгенологические методы. Рентгенограммы обрабатываются методом апостериорного компьютерного анализа, позволяющего более точно выявить качество восстановления ткани. В частности, области дефекта на разных этапах регенерации представляются в виде нормализованных гистограмм яркости.
До операции на гистограммах примерно в равной степени присутствовали уровни яркости, соответствующие двум основным составляющим кости: «мягкотканевой» - коллагену и «солевой» - минеральному компоненту. Через 20 дней после операции на гистограмме, соответствующей области дефекта, определялось большее присутствие и большая интенсивность уровней яркости, близких к плотностной характеристике костной ткани. Дальнейшие наблюдения подтвердили формирование костной ткани, так как плотность в зоне дефекта увеличивалась [6]. Слабый контраст является наиболее распространенным свойством рентгеновских изображений и существенно затрудняет их визуальное восприятие. Это обусловлено тем, что диапазон изменения яркости изображения относительно усредненного значения очень мал. При этом яркость меняется не от черного до белого, а в небольшом диапазоне от темно-серого до светло-серого. Более того, рентгеновский снимок по своей природе несет информацию в аналоговой форме.
В данной статье предлагается новый метод оценки динамики восстановления тканей на основе рентгеновских снимков, поскольку такие снимки являются простейшими безболевыми средствами контроля процесса лечения внутренних органов человека. Затронуты также некоторые вопросы методики получения сжатого образа для рентгеновского изображения, важными достоинствами которой являются высокая разрешающая способность (в пределах одного пиксела), инвариантность по отношению к сдвигу и изменению масштаба анализируемого сигнала, а также повышенное быстродействие.
Исходные данные. В качестве исходных данных рассматривается фрагмент (рис. 1) рентгеновского снимка, характеризующего плотность ткани, в которую имплантировано средство «ЛитАр». Очевидно, что в таком виде информация о действии имплантата в силу ее аналоговой природы воспринимается медицинским персоналом только в качественной форме. Устранение этого негативного момента требует преобразования данных в дискретную форму, что достигается при помощи процедуры Read (чтение) программы Matcad. Применив эту процедуру к рассматриваемому фрагменту, получаем массив пикселов в виде числовой таблицы, содержащей данные с интенсивностями свечения каждого пиксела. Каждый из пикселов имеет соответствующее значение яркости, от белого до черного (255-0 пикс.). Очевидно, что в таком виде информация о текущем состоянии поврежденной ткани трудно воспринимается медицинским персоналом в силу большого объема числовых данных.
Предположим, что исследуемый фрагмент рентгеновского снимка после преобразования программой Matcad представлен табл. 1 двумерным массивом.
Таблица 1
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
||
0 |
1,27 |
1,23 |
1,25 |
1,26 |
1,12 |
1,09 |
1,07 |
1,06 |
|
1 |
1,27 |
1,24 |
1,2 |
1,18 |
1,17 |
1,16 |
1,15 |
1,10 |
|
2 |
1,25 |
1,24 |
1,19 |
1,18 |
1,16 |
1,14 |
1,10 |
1,10 |
|
3 |
1,20 |
1,22 |
1,18 |
1,15 |
1,14 |
1,08 |
1,06 |
1,05 |
|
4 |
1,19 |
1,21 |
1,25 |
1,23 |
1,12 |
1,09 |
1,06 |
1,04 |
|
5 |
1,08 |
1,2 |
1,23 |
1,17 |
1,05 |
1,07 |
1,05 |
1,02 |
|
6 |
1,05 |
1,18 |
1,06 |
1,08 |
1,00 |
1,05 |
1,03 |
0,95 |
|
7 |
0,97 |
1,12 |
1,09 |
1,12 |
0,97 |
1,01 |
1,00 |
0,93 |
Далее этот массив формально рассматривается как дискретное множество , где l и k (l = k) - его координаты.
Методы и средства. Одним из свойств преобразования Фурье в экспоненциальном базисе является свойство инвариантности амплитудно-частотного спектра к сдвигам сигнала , благодаря которому значительно упрощается проблема сопоставления различных спектров [7, 8]. Кроме того, переход в частотную область позволяет эффективно обрабатывать сигнал с целью защиты от помех и сжатия информации. Несмотря на достоинства такого подхода, необходимо отметить, что при получении спектра с помощью экспоненциального базиса теряется в явном виде структурная информация о сигнале. Действительно, для определения значения любой k-той гармоники фактически вычисляется (интегрируется) взаимная энергия между сигналом и k-той базисной функцией, в результате величины значения соответствующих взаимных энергий могут быть одинаковыми для сигналов различной формы. Возможны случаи, когда одинаковые значения спектральных составляющих соответствуют разным по форме сигналам. Другими словами, при использовании экспоненциального базиса в явном виде отсутствует информация о локальных особенностях , например о резких всплесках. Следует также добавить, что выполнение рассмотренных преобразований в экспоненциальном базисе требует значительных временных затрат. При большом k анализ полученного спектра осложняется, особенно если данные имеют сложную структуру, а при отбрасывании части гармоник теряется информация о сигнале.
Очевидно, что в случае со сложноструктурированными данными необходимо сжатие, причем оно должно учитывать самые небольшие изменения сигнала. В связи с этим интерес представляет возможность использования ортогональных систем кусочно-постоянных базисных функций, позволяющих исключить такие временные операции, как вычисление тригонометрических функций и умножение вещественных величин на соответствующие значения анализируемого сигнала. Необходимо подчеркнуть, что эта проблема приобретает особую актуальность при обработке больших массивов цифровых данных. Среди различных систем кусочно-постоянных базисных функций в настоящее время наиболее изученной является система функций Уолша и подобных ей, но при применении данной системы функций при преобразовании Фурье не выполняется теорема о сдвиге сигнала, т. е. полученный спектр находится в зависимости от сдвига сигнала. Таким образом, в зависимости от базисных функций при использовании преобразования Фурье можно получать спектральные характеристики сигнала, которые зависят либо от сдвига, либо от частоты. Поэтому в первом случае отсутствует информация о частотном характере, а во втором - информация о структуре сигнала.
Избежать данных недостатков можно, применив вейвлет-преобразование [9]. Семейство вейвлет-функций
(1)
генерируется одной главной функцией путем растяжения и сдвига, где и определяют масштаб и сдвиг соответственно. Главной проблемой вейвлет-преобразования, несмотря на его преимущества, является необходимость подбора конкретного вейвлета в выражении (1) и низкое быстродействие метода, что отрицательно скажется при обработке больших массивов данных.
Для решения перечисленных проблем предложена базисная комплексная система импульсных функций.
Из всех известных линейных спектральных форм представления сигнала в силу ряда причин наибольшее распространение получило его разложение по ортогональной системе экспоненциальных базисных функций [7, 8, 9]. Однако процедуры спектрального преобразования сигналов, представленных в двухкоординатном (плоскостном) формате, требуют значительных временных затрат, но не устраняют их негативные свойства.
Для решения перечисленных проблем предлагается комплексная система импульсных функций, с помощью которой можно осуществить нелинейное спектральное преобразование, отвечающее как требованию инвариантности к сдвигу, так и требованию быстродействия. Для этого рассмотрим систему комплексных функций [8]:
. (2)
Предложенная ортогональная система функций определяется на дискретном множестве
,
где , - число вертикальных и горизонтальных отсчетов, составляющих период дискретного сигнала , который определяется из условий:
.
Функции и формируются следующим образом:
, для ; при и
,
где представляет собой единичный импульс, задержанный на d и b сдвигов:
В соответствии с формулой Эйлера система комплексных функций (1) может быть выражена через экспоненциальные функции:
(3)
где k изменяется от 0 до с шагом , l изменяется от 0 до с шагом .
Для получения спектра, инвариантного к сдвигам сигнала , необходимо производить циклические сдвиги функций (2) относительно этого сигнала. С учетом сдвигов система функций (3) может быть записана в виде:
где , .
Функции и формируются в результате сдвигов на и подынтервалов функций и вдоль горизонтальной и (или) вертикальной осей.
В отличие от одномерных функций [11] при вычислении спектра данные функции сдвигаются не только в горизонтальном направлении, но и в вертикальном и диагональном направлениях.
В результате этих сдвигов вычисляются частные скалярные произведения между каждой базисной функцией и и всеми возможными значениями сигнала. На заключительном этапе двумерный амплитудно-частотный спектр анализируемого сигнала может быть получен согласно следующему выражению:
, (4)
Где
В данных выражениях представляет собой значение анализируемого сигнала в точках с координатами , где и , а л и м - целочисленные переменные, изменяющиеся в диапазонах от 0 до и от 0 до соответственно. Соотношение (4) указывает также на одно важное свойство, а именно: результат не зависит от масштаба переменных l и k.
Результаты. Спектр исходного фрагмента, полученный с помощью предложенного метода, представлен табл. 2.
Таблица 2
№ гармоники |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
Значение гармоники |
71,93 |
1,6753 |
46,413 |
2,4666 |
46,4931 |
0.2460 |
0,49201 |
0,57801 |
0,28160 |
Пусть с течением времени произошли некоторые изменения ткани в окрестности имплантата, например отсчет с координатами 6/3 (шестой столбец, 3-я строка в табл. 1) поменял свое значение с 1,06 на 1,05. Тогда значения спектральных компонент будут такими, как показано в табл. 3.
Таблица 3
№ гармоники |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
Значение гармоники |
71,92 |
1,6749 |
46,408 |
2,4650 |
46,4833 |
0,24158 |
0,494 |
0,58603 |
0,29155 |
Выводы
Предложенная методика позволяет преобразовывать большие объемы данных в компактную форму, позволяющую количественно и качественно оценивать динамику имплантата, т. е. отслеживать малейшие изменения в биологических тканях путем анализа рентгеновских снимков.
При формировании двумерного спектра исходного фрагмента согласно преобразованию (3) фактически удалось исключить времяемкую операцию умножения, так как импульсные функции принимают только три значения: -1, 0 и +1. Поэтому на основе разработанной системы импульсных функций можно значительно быстрее сформировать сжатый спектр анализируемого сигнала, чем, например, при использовании экспоненциальной системы комплексных функций. Кроме того, замечательным свойством данного спектра является свойство его инвариантности по отношению к смещениям анализируемого сигнала.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Леонтьев В.К., Литвинов С.Д., Судакова Т.В. Имплантационные материалы для замещения дефектов костной и хрящевой ткани // Российский вестник дентальной имплантологии. - 2003. -№2. - С. 10-19.
2. Краснов А.Ф., Литвинов С.Д. Медицинская практика применения материала «ЛитАр»: история и реальность // Ортопедия, травматология и протезирование (Харьков). - 2003. - №3. - С. 136-142.
3. Марков И.И., Литвинов С.Д., Марков А.И. Имплантационный материал «ЛитАр» индуцирует ангиогенез // Морфологические ведомости. - №1-2. - 2003. - С. 74-76.
4. Куликов А.Н., Литвинов С.Д. Имплантационный материал «ЛитАр» и сочетанная травма. Современные подходы к диагностике, лечению и реабилитации пострадавших с сочетанными повреждениями: сб. материалов Межд. конф., г. Москва, 25-26 мая 2006 г. - М.: Объединенная редакция МВД России, 2006. - С. 43-48.
5. Тарасов Ю.В., Литвинов С.Д., Серегин А.С. Применение материала «ЛитАр» в детской оториноларингологии: мат-лы XVII съезда оториноларингологов России (тезисы). - Н. Новгород, 7-9 июня 2006 г. - СПб, 2006. - С. 490.
6. Литвинов С.Д. Наноразмерный композитный материал «ЛитАр» - универсальный имплантат: монография. - Самара: Книга, 2008. - 250 с.
7. Гоноровский И.С. Радиотехнические цепи и сигналы. - М.: Советское радио, 1977. - 608 с.
8. Баскаков С.И. Радиотехнические цепи и сигналы. - М.: Высшая школа, 1983. - 536 с.
9. Переберин А.В. О систематизации вейвлет-преобразований // Вычислительные методы и программирование. - 2001. - Т. 2. - С. 15-40.
10. Бахрушин А.П. Спектральный анализ одномерных и двумерных сигналов на основе различных базисных функций. - Владивосток: Дальнаука, 2004. - 170 с.
Размещено на Allbest.ru
...Подобные документы
Ряд способов получения МР-томограмм. Получение изображения тонких слоев тела человека в любом сечении. Оценка магнитно-резонансной томографии спинного мозга. Отграничение кистозного опухолевого компонента от сопутствующих сирингомиелических изменений.
презентация [282,8 K], добавлен 29.03.2015Появление идеи создания особого биопластического материала, способного восстанавливать дефекты покровных тканей. Получение биокожи - биополимера на основе гиалуроновой кислоты - основного межклеточного компонента кожи. Лечение с помощью биоматрикса.
презентация [5,3 M], добавлен 27.05.2016Понятие и принципы организации гормональной контрацепции как наиболее эффективного метода предохранения от нежелательной беременности. Их разновидности и оценка практической эффективности: по дозировке эстрогенного компонента, прогестагену и фазности.
презентация [3,8 M], добавлен 12.04.2015Метод рациональной психотерапии. Метод отвлечения. Метод ролевого аутотренинга. Метод молниеносной мышечной релаксации. Метод пассивного очищающего дыхания. Метод динамического мышечного напряжения. Метод добрых дел.
монография [26,7 K], добавлен 21.08.2007Значение пищеварения - главного компонента функциональной системы питания - для жизнедеятельности организма. Строение пищеварительной системы, ее функции. Морфологические и физиологические особенности пищеварительных органов ребёнка и взрослого человека.
курсовая работа [1022,8 K], добавлен 19.04.2016Четыре компонента повреждения сердца, составляющие тетраду Фалло: подаортальный дефект межжелудочковой перегородки, обструкция выхода из правого желудочка, гипертрофия его миокарда. Диагностика данной патологии на основе рентгенографии грудной клетки.
презентация [1,3 M], добавлен 21.04.2016Оценка физического развития методом сигмальных отклонений. Метод оценки физического развития по шкалам регрессии. Комплексная оценка физического развития. Центильный метод оценки физического развития. Основные параметры морфологического статуса.
презентация [894,9 K], добавлен 23.04.2015Оценка риска для здоровья человека. Характеристика вредных эффектов, способных развиться в результате воздействия факторов окружающей среды на группу людей. Передача информации о риске. Анализ продолжительности воздействия факторов риска на человека.
презентация [211,5 K], добавлен 01.10.2014Патопсихологические и нейропсихологические особенности при органических патология сосудистого генеза. Виды патологий памяти. Нарушение динамики мнестической деятельности, опосредованной памяти, мотивационного компонента памяти, интеллекта и внимания.
дипломная работа [425,4 K], добавлен 06.05.2011Понятие ДВС-синдрома, его формы, признаки и стадии. Основы проблемы свертывания крови и агрегации тромбоцитов. Рассмотрение процесса активации плазменного компонента гемостаза, распространение тромбоза, проявление развернутого геморрагического синдрома.
презентация [186,7 K], добавлен 26.05.2014Высокая патогенетическая значимость иммунного компонента в развитии атеросклероза и его клинических проявлений. Роль модифицированных ЛПНП в инициации аутоиммунного ответа. Прогрессирование атеросклероза и нестабильности атеросклеротической бляшки.
реферат [27,6 K], добавлен 20.03.2009Четыре категории общей неотложной стоматологической патологии. Анатомия зуба, два основных компонента нормального пародонта. Боль одонтогенной природы при прорезывании зубов, наличии околозубного абсцесса. Пароксизмальная боль неврологической этиологии.
реферат [26,8 K], добавлен 08.06.2009Отек как увеличение внесосудистого (интерстициального) компонента внеклеточной жидкости. Патофизиологические причины отеков, их признаки, последствия. Внезапный отек одной из нижних конечностей. Лимфатические и статические отеки. Тромбоз подмышечной вены.
презентация [540,9 K], добавлен 21.01.2014Причины возникновения пародонтоза, его характерые признаки. Патогенетическая классификация заболеваний пародонта, причины возникновения воспалительного процесса. Факторы, влияющие на развитие зубного налета. Влияние аллергического компонента на пародонт.
презентация [4,2 M], добавлен 29.05.2014Аутосомно наследуемый дефицит или аномальная структура синтезируемого в эндотелии сосудов макромолекулярного компонента. Этиология и патогенез болезни Виллебранда, его типы и диагностические критерии. Клиническая картина заболевания, методика лечения.
презентация [174,8 K], добавлен 16.02.2015Методы планирования семьи, направленные на снижение числа незапланированных беременностей. Анализ состояния коагуляционного компонента гемостаза у женщин, использующих эстроген-гестагенный препарат для контрацепции с эндовагинальным путем введения.
статья [27,9 K], добавлен 02.08.2013Современное состояние гирудотерапии и ее краткая история как метода лечения. Основные физиологические механизмы воздействия медицинской пиявки на организм человека. Отвлекающее и противовоспалительное действие гирудотерапии, ее лечебный эффект.
реферат [86,2 K], добавлен 17.05.2009Стеклоиономерный цемент как это материал, который состоит из основного компонента стекла и кислотного полимера и отвердевает посредством кислотно-основной реакции между этими компонентами. Его свойства и применение в терапевтической стоматологии.
презентация [1,6 M], добавлен 10.01.2016Основные трудности диагностики и терапии мигрени. Феномен мигренозной боли. Преобладающие аффективные компоненты при феномене мигренозных болей. Оптимизация различных методик, направленных на адаптацию и модификацию существующих терапевтических моделей.
презентация [257,6 K], добавлен 20.11.2013Использование одноразового инструментария и шприцов. Металлические конструкции для остеосинтеза. Самостоятельное изготовление шовного материала. Возникновение гнойно-септических осложнений. Основные направления профилактики имплантационной инфекции.
презентация [458,0 K], добавлен 29.10.2014