способ получения томографического изображения тела и электроимпедансный томограф
Классы МПК: | A61B5/05 измерение с помощью электрического тока или магнитных полей для диагностических целей |
Автор(ы): | Корженевский А.В., Культиасов Ю.С., Черепенин В.А. |
Патентообладатель(и): | Корженевский Александр Владимирович, Черепенин Владимир Алексеевич |
Приоритеты: |
подача заявки:
1996-12-11 публикация патента:
10.03.1999 |
Изобретение представляет способ получения томографического изображения тела и электроимпедансный томограф. Серию контактных электродов размещают на поверхности тела и поочередно подключают источник тока к произвольным парам электродов. На каждой паре остальных электродов измеряют разность потенциалов, определяют опорные значения разностей путем аппроксимаций измеренного распределения. Реконструкцию изображения пространственного распределения проводимости производят путем обратного проецирования вдоль эквипотенциалей. Томограф содержит систему контактных электродов, устройства возбуждения тока, измерения разностей потенциалов с дифференциальным усилителем, блок микропроцессора, аналоговые коммутаторы, цепь компенсации синфазной составляющей напряжений на каждой паре электродов и цепь компенсаций контактной разности потенциалов. Способ позволяет визуализировать распределение абсолютной проводимости в сечении человеческого тела. 2 c. и 5 з.п. ф-лы, 10 ил.
Рисунок 1, Рисунок 2, Рисунок 3, Рисунок 4, Рисунок 5, Рисунок 6, Рисунок 7, Рисунок 8, Рисунок 9, Рисунок 10
Формула изобретения
1. Способ получения томографического изображения тела, включающий размещение серии контактных электродов на поверхности тела, поочередное подключение источника тока к парам электродов, измерение разности потенциалов Uim(j) на каждой паре остальных электродов, определение опорных значений разностей потенциалов Uir(j) путем аппроксимации измеренного распределения разностей потенциалов Uim(j) в соответствии с выражениемUir(j) = Ci1fi1(j) + Ci2fi2(j) + Ci3,
отличающийся тем, что осуществляют подключение источника тока к произвольным парам электродов, а реконструкцию изображения пространственного распределения проводимости тела путем обратного проецирования проводимости тела путем обратного проецирования вдоль эквипотенциалей, производят в соответствии с выражениями:
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127075/2127075-4t.gif)
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127011/955.gif)
где Wet, Wzt - весовые факторы, определенные в соответствии с процедурой обратного проецирования в направлении от "левого" и от "правого" пересечения эквипотенциали с границей тела, соответственно:
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127043/931.gif)
Ulmt,rt - напряжения, измеренные на левом и правом "концах" эквипотенциальной линии, проходящей через данную точку реконструируемого поперечного сечения;
Ulrt,rt - опорные разности потенциалов, соответствующие телу с однородной проводимостью;
i - номер возбуждающей пары электродов;
j - номер измеряющей пары электродов;
f1 (j) - заданное распределение напряжения между соседними электродами вдоль границы эталонного объекта;
f2 (j) - сигналы, обусловленные паразитными связями;
Ci
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127010/945.gif)
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127010/945.gif)
2. Способ по п.1, отличающийся тем, что сигналы fi2(j), обусловленные паразитными связями, определяют путем проведения измерений. 3. Способ по п.1, отличающийся тем, что для диагностики органов с изменяющейся во времени проводимостью производят автоматически серию измерений разностей потенциалов последовательно во времени, выполняют спектральное Фурье - преобразование временных зависимостей полученных результатов и производят реконструкцию изображений пространственного распределения проводимости органов для каждой частотной компоненты. 4. Способ по п.1, отличающийся тем, что производят реконструкцию изображения пространственного распределения абсолютной проводимости тела путем нормировки полученных значений проводимости, исходя из того, что наименьшее значение проводимости соответствует проводимости костных тканей, а наибольшее значение проводимости соответствует проводимости крови. 5. Электроимпедансный томограф, содержащий систему контактных электродов, устройство возбуждения тока, устройство измерения разностей потенциалов с дифференциальным усилителем на входе, микропроцессорную схему управления, входные и выходные аналоговые коммутаторы, входы которых подключены к контактным электродам, выходы входных коммутаторов - к входу дифференциального усилителя, а выходы выходных коммутаторов - к входу дифференциального усилителя, а выходы выходных коммутаторов - к устройству возбуждения напряжения, отличающийся тем, что томограф содержит цепь компенсации синфазной составляющей напряжений на каждой паре электродов, цепь компенсаций контактной разности потенциалов, при этом цепь компенсации синфазной составляющей напряжений выполнена в виде цепи обратной связи, содержащей операционный усилитель, выход синфазного сигнала дифференциального усилителя подключен к инвертирующему входу операционного усилителя, выход которого подключен к устройству возбуждения напряжения. 6. Томограф по п.5, отличающийся тем, что усилитель устройства возбуждения напряжения выполнен компаратором, выход которого подключен к блоку микропроцессора. 7. Томограф по п.5, отличающийся тем, что цепь компенсации контактной разности потенциалов выполнена в виде цепи обратной связи, включающей последовательно соединенные аналоговый ключ и интегратор, при этом выход интегратора подключен к входу установки нуля дифференциального усилителя, а вход аналогового ключа - к выходу усилителя устройства измерения разностей потенциалов.
Описание изобретения к патенту
Изобретение относится к области медицины, а более конкретно к способам диагностики с помощью приборов для получения томографического изображения тела пациента. Известны способы получения томографического изображения тела пациента, основанные на измерении пространственного распределения физического поля или излучения, пронизывающего объект, и последующей реконструкции изображения пространственного распределения измеренного параметра математическими методами свертки и обратной проекции (Физика визуализации изображений в медицине./Под ред. С. Уэбба. М.: Мир, 1991 г. с. 105 - 216.)Известны томографы, основанные на использовании рентгеновского излучения или ядерного магнитного резонанса (ЯМР) (Физика визуализации изображений в медицине./Под ред. С. Уэбба. М.: Мир, 1991 г. с. 105 - 216.)
Известные способы обеспечивают высокую разрешающую способность. Однако использованные для диагностики сложные установки рентгеновского излучения или ядерного магнитного резонанса дорогостоящи и сложны в эксплуатации, процедура обследования достаточно длительна, кроме того, излучение, пронизывающее тело, не является безвредным для пациента и обслуживающего персонала. Известен способ получения томографического изображения тела для диагностики в медицине, основанный на использовании электрического тока в качестве средства, зондирующего исследуемый объект - электроимпедансная томография (Патент Великобритании N 2119520 A, кл. A 61 B 5/05, 1983 г.) В известном способе на поверхности тела пациента располагают серию контактных электродов, последовательно подключают источник тока к парам электродов, производят измерения разностей потенциалов (напряжения) между каждой парой остальных электродов, возникающих из-за протекания тока через объект, определяют опорные значения напряжений в предположении об однородности электрической проводимости объекта или измеряют их на том же объекте в другой момент времени, если проводимость объекта меняется, и производят реконструкцию изображения пространственного распределения проводимости тела или измерений его проводимости, произошедших между двумя измерениями путем обратного проецирования относительных разностей измеренных и опорных напряжений вдоль эквипотенциалей электрического поля. Установлено, что проводимость биоткани зависит от ее физиологических свойств. Изображение распределения проводимости тела дает возможность увидеть кости, мягкие ткани, сосуды. Способ осуществления на электроимпедансном томографе, содержащем систему контактных электродов, устройство возбуждения тока, устройство измерения разностей потенциалов с дифференциальным усилителем на входе, блок микропроцессора, входные и выходные аналоговые коммутаторы, входы которых подключены к контактным электродам, а выходы - к устройству возбуждения тока и ко входу дифференциального усилителя (Патент Великобритании N 2119520 A, кл A 61 B 5/05, 1983 г.). Однако использование способа в клинической практике до настоящего времени тормозилось, в частности, нерешенностью проблемы получения абсолютных или "статических" изображений удовлетворительного качества при проведении измерений на человеческом теле. Существующие томографы позволяют получать только динамические томограммы - изображения измерений проводимости, произошедшие между двумя измерениями, малоинформативные для медицинских приложений. Отсутствие визуализации статических объектов объясняется невозможностью полного решения обратной задачи реконструкции проводимости из-за трудности получения опорных значений разностей потенциалов при неизвестных точно геометрии граничной поверхности исследуемого объекта и положения на ней измерительных электродов. Получение возможности визуализации распределения абсолютной проводимости в сечении человеческого тела с высокой скоростью сбора данных стало возможно с помощью компактного томографа с управлением всеми функциями его измерительной системы с персонального компьютера, осуществляющего обработку, визуализацию и хранение данных. (В. А. Черепенин, А. В. Корженевский и др. Электроимпедансный томограф: новые возможности. IX Международная конференция по электрическому биоимпедансу. Гейдельберг, Германия, 1995, с. 430 - 433. - Cherepenin V.A., Korjenevsky A.V et al. The electrical impedance tomograph: new capabilities. IX International Conference on Electrical Bio-Impedance, Proceedings. - Heidelberg, 1995, p. 430 - 433). В способе получения томографического изображения тела, включающем размещение серии контактных электродов на поверхности тела, поочередное дипольное подключение источника тока к парам соседних электродов, измерение разности потенциалов (напряжений) между каждой парой остальных электродов, определение опорных значений разностей потенциалов и реконструкцию изображения пространственного распределения проводимости тела путем обратного проецирования вдоль эквипотенциалей взвешенных относительных разностей опорных и измеренных напряжений, опорные значения разностей потенциалов uir(j) определяют путем аппроксимации измеренного распределения разностей потенциалов uin(j) в соответствии с выражением:
uir(j) = ci1fi1(j) + ci2fi2(j) + ci3 (1),
где
i - номер возбуждающей пары электродов;
j - номер измеряющей пары электродов,
f1(j) - заданное распределение напряжения между соседними электродами вдоль границы эталонного объекта,
fi2(j) - сигналы, обусловленные паразитными связями,
ci
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127010/945.gif)
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127010/945.gif)
В описанном решении, используя аппроксимацию измеренных данных uim(j) гладкими зависимостями из набора простых линейно независимых функций, удается построить опорный набор данных, не содержащий информацию о внутренней структуре объекта. Этот набор совместно с исходным набором, включающим вариации, отражающие внутреннюю структуру объекта, используется для реконструкции абсолютной проводимости объекта. Измеренные разности потенциалов могут содержать значительные систематические погрешности, вызванные, в основном, паразитным проникновением сигналов из канала в канал в интегральных коммутаторах и входных цепях томографа. В известных решениях при дифференциальной томографии они не оказывают существенного влияния на реконструируемое изображение, поскольку паразитные составляющие имеются как в текущих, так и в опорных данных, и взаимно компенсируются. При реконструкции распределения абсолютной проводимости эти помехи вызывают появление артефактов и существенно снижают качество изображения. Для устранения их влияния в состав базовых функций, используемых для аппроксимации входных данных, может быть включен набор паразитных сигналов, которые могут быть, например, измерены с помощью эталонного объекта, для которого известны (вычислены) результаты идеально точных измерений. Наилучшие результаты получаются при использовании комбинации трех функций, приведенных в выражении (1). В качестве распределения f1(j) - задается распределение напряжения между соседними электродами вдоль границы цилиндрического объекта с однородной проводимостью при подключении источника тока к паре смежных электродов. Разработанный алгоритм реконструкции распределения проводимости позволяет получать с помощью томографа привычные в медицине и весьма информативные "статические" изображения, характеризующие физиологические состояния органов и тканей. Однако низкая чувствительность при невысокой разрешающей способности ограничивает область применения этого метода. В настоящем изобретении предлагается способ получения томографического изображения тела, позволяющий получить качественное изображение проводимости объема тела с высокой чувствительностью и удовлетворительной разрешающей способностью, динамично отражающее состояние внутренних органов и тканей с высокой достоверностью. Способ позволяет более чем на порядок увеличить отношение сигнал-шум и тем самым повысить чувствительность и разрешающую способность реальных устройств, что позволяет найти предложенным устройствам широкое применение в медицинской практике и при клинических исследованиях. Новые свойства предложенного способа позволяют выявлять и измерять структуры и процессы, определение которых с помощью рентгеновских лучей или ядерного магнитного резонанса затруднено или невозможно. Способ позволяет обеспечить простое и безопасное исследование пациентов не только в клинических условиях, но и обычных кабинетах и лабораториях. Задача повышения качества томографического изображения решается за счет разработки нового эффективного метода измерения разностей потенциалов и реконструкции изображения пространственного распределения проводимости тела. Поскольку при проведении измерений на теле человека величина инжектируемого тока ограничена условием безопасности воздействия, величина измеряемых прибором сигналов оказывается весьма небольшой и серьезное влияние на качество реконструируемых изображений оказывает отношение амплитуды измеряемых сигналов к величине собственных шумов прибора и внешних электрических помех на частоте измерения. В настоящем изобретении за счет инжекции тока через два необязательно смежные электрода, достигается увеличение отношения сигнал-шум, что позволяет повысить чувствительность и во многих случаях разрешающую способность реальных устройств. Для этого более общего случая, по сравнению с дипольной инжекцией, как в прототипе, необходимо скорректировать метод реконструкции изображения пространственного распределения проводимости. Исходными данными являются разности потенциалов между соседними электродами, закрепленными на коже пациента вдоль замкнутого контура, охватывающего тело, при подключении каких-либо двух электродов к источнику тока. При использовании N электродов имеем N профилей, соответствующих каждому варианту подключения источника тока, содержащих по (N - 4) отсчетов разностей потенциалов между свободными парами электродов. В отличие от дипольного метода инжекции процедура обратного проецирования для произвольного метода инжекции тока описывается следующим образом:
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127075/2127075-2t.gif)
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127011/955.gif)
где
Wlt, Wn - весовые факторы, определенные в соответствии с процедурой обратного проецирования в направлении от "левого" и от "правого" пересечения эквипотенциали с границей тела, соответственно,
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127043/931.gif)
ulmt,rt - напряжения, измеренные на левом и правом "концах" эквипотенциальной линии, проходящей через данную точку реконструируемого поперечного сечения,
ulrt,rt - опорные разности потенциалов, соответствующие телу с однородной проводимостью. Кроме того, уменьшение ошибок измерения в изобретении достигается введением цепи контроля качества контактов. Изобретение иллюстрируется чертежом, где на фиг. 1 изображена блок-схема измерительной системы томографа, на фиг. 2 показана зависимость проводимости в области сердца от времени, полученная при измерении с дипольной инжекцией, как в известных решениях, а на фиг. 3 - та же зависимость, полученная при измерении с полярной инжекцией, согласно изобретению, на фиг. 4 - спектр этой зависимости, на фиг. 5 показана импедансная томограмма грудной области пациента Y, на фиг. 6 приведено пространственное распределение спектральной амплитуды проводимости в сечении грудной области на частоте, соответствующей частоте сердечных сокращения пациента Y. На фиг. 7 показана импедансная томограмма грудной области пациента X, на фиг. 8 - обзорная рентгенограмма грудной области пациента X, на фиг. 9 приведена импедансная томограмма нижней части левой ноги пациента Z, фиг. 10 приведена импедансная томограмма верхней части левой ноги пациента Z. На блок-схеме томографа (фиг. 1) обозначены: контактные электроды 1 - 16, формирователь 17 импульсных напряжений, преобразователь напряжение-ток 18, компаратор 19, фазоимпульсный модулятор 20, выходные аналоговые коммутаторы 21 и входные аналоговые коммутаторы 22, каждый из которых имеет по N входов (N - количество электродов) для подключения к электродам и по одному выходу, дифференциальный усилитель 23, операционный усилитель 24, интегратор 25, аналоговый ключ 26, усилитель 27, синхронный детектор 28, ключ 29, интегратор 30, ключ 31, аналогово-цифровой преобразователь 32, микропроцессор 33, формирователь 34. Работа системы сбора данных электроимпедансного томографа сводится к измерениям потенциалов на поверхности тела человека с помощью контактных электродов 1 - 16 при протекании через него слабого электрического тока. Возбуждающая пара электродов подключается к источнику тока, а на остальных парах последовательно измеряются напряжения, обусловленные действием тока в цепи возбуждающих электродов. Описываемый вариант электроимпедансного томографа построен по одноканальной схеме и рассчитан на работу с 16-ю электродами. Время выполнения каждого измерения составляет 350 мс , а частота выполнения полных циклов измерений - 11 кадров в секунду. Измерения выполняются при импульсном воздействии сигнала специальной формы со средней частотой 8 кГц. Использование такого сигнала позволило упростить аппаратуру по сравнению с известными системами с синусоидальным возбуждением при сохранении высокой точности и быстродействия. Устройство возбуждения напряжения содержит прецизионный формирователь 17 импульсных напряжений, усилитель 18 с токовым выходом и фазоимпульсный модулятор 21. Формирователь 17 вырабатывает импульсное напряжение, представляющее собой 2 периода меандра, сдвинутые друг относительно друга на 180o, с регулируемой амплитудой, которая устанавливается 8-разрядным кодом с блока микропроцессора 34. Использование импульсного тока в томографе позволило отказаться от прецизионного формирователя синусоидального сигнала. Преобразователь напряжение-ток 18 с токовым выходом собран на двух операционных усилителях. Он обладает стабильным коэффициентом усиления, имеет большое выходное сопротивление и обеспечивает независимость величины выходного тока от сопротивления цепи возбуждающей пары электродов. В процессе измерений может ухудшаться качество электрического контакта электродов с кожей человека. Это приводит к ошибкам измерений, возникающим из-за нарушения работы источника тока при достижении максимально допустимого выходного напряжения (порядка напряжения источника питания) а также увеличения наводок и погрешностей во входных цепях устройства измерения разностей потенциалов. Поэтому важно иметь возможность контролировать в процессе измерений качество контакта электродов с кожей и определять, какой именно контакт имеет повышенное электрическое сопротивление и требует вмешательства. Для осуществления этих функций в схему прибора введено устройство контроля качества контактов, содержащее компаратор 19, сравнивающий выходное напряжение устройства возбуждения напряжения с заданными постоянными положительным +U и отрицательным -U напряжениями, равными максимально допускаемому положительному и отрицательному выходному напряжению источника тока. Сигнал с выхода компаратора 19, соответствующий выходу напряжения за допустимые пределы, поступает на микропроцессорный блок и вместе с текущими адресами активных контактов передается на персональный компьютер PC. Это позволяет в процессе измерений вывести на экран сообщение о наличии плохих контактов и во многих случаях точно установить номер электрода, требующего внимания. Фазоимпульсный модулятор 20 собран на аналоговых ключах и один раз за период одного измерения переворачивает на 180o фазу выходного сигнала возбудителя тока. Выбранная пара электродов подключается к возбудителю двумя выходными аналоговыми коммутаторами 21 с общими адресными шинами управления. Напряжения с приемных пар электродов 1 - 16 контактной системы через входные аналоговые коммутаторы 22 подаются на устройство измерения разностей потенциалов, состоящее из дифференциального усилителя 23, усилителя 27, аналоговых ключей 29 и 31, синхронного детектора 28, интегратора 30 и аналого-цифрового преобразователя 32. Полезная информация заключена в дифференциальной составляющей напряжений на паре приемных электродов. Синфазная составляющая этих напряжений должна подавляться дифференциальным усилителем 23. Поскольку полного подавления синфазной помехи добиться не удается, в томограф введена цепь компенсации синфазной составляющей с использованием операционного усилителя 24. На инвертирующий вход этого усилителя поступает синфазная составляющая пары измеряемых напряжений, а его выход связан с возбуждающей парой через модулятор 20. На измеряемый дифференциальный сигнал накладывается контактная разность потенциалов, которая для электродов из нержавеющей стали, используемых в наших лабораторных образцах, составляет
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127020/177.gif)
![способ получения томографического изображения тела и электроимпедансный томограф, патент № 2127075](/images/patents/344/2127075/2127075-3t.gif)
Класс A61B5/05 измерение с помощью электрического тока или магнитных полей для диагностических целей