CC BY
38
Fedotov Aleksandr Aleksandrovich - e-mail: fedoaleks@yandex.ru, 303-11, Kirova avenue, Samara, 443114, Russia; phone: +78469277668; the department of radioengineering and biomedical engineering; engineer; post raduate student.
Lebedeva Ekaterina Petrovna - Samara State Medical University; e-mail: milaorion@mail.ru; 89, Chapaeva street, Samara, 443099; Russia; phone +78462674549: intern.
УДК 621.317.332: 616-073.7
А.В. Киреев
ОБ УСЛОВИЯХ СЪЁМА ПЕРВИЧНОЙ ИЗМЕРИТЕЛЬНОЙ ИНФОРМАЦИИ, ИСПОЛЬЗУЕМОЙ ПРИ РЕКОНСТРУКЦИИ ПРОСТРАНСТВЕННОГО РАСПРЕДЕЛЕНИЯ ЭЛЕКТРОПРОВОДНОСТИ БИОЛОГИЧЕСКИХ ТКАНЕЙ
Сформулированы требования, предъявляемые к точности измерения параметров импеданса на поверхности тела с позиции визуализации пространственного распределения электропроводности внутренних тканей и органов. Выделены основные факторы, ограничивающие точность измерения импеданса и установлена их взаимосвязь с результирующей погрешностью измерений. Получена оценка предельно достижимой точности измерения импеданса корреляционным методом и определены оптимальные параметры измерительной системы, обеспечивающие её достижение. Предложен способ учёта нестационарно-сти импеданса с помощью сглаживающих окон полиномиального скользящего среднего, позволяющий дополнительно повысить точность измерения импеданса.
Томография; импеданс; погрешности; нестационарность; сглаживание.
A.V. Kireev
ABOUT CONDITIONS OF OBTAINING PRIMARY MEASURING INFORMATION USED AT RECONSTRUCTION OF SPATIAL DISTRIBUTION OF AN ELECTRIC IMPEDANCE OF BIOLOGICAL FABRICS
The requirements shown to accuracy of measurement of parameters of an impedance on a surface of a body from a position of visualization of spatial distribution of an elektroprovodnost of internal fabrics and bodies are formulated. The major factors limiting accuracy of measurement of an impedance are allocated and their interrelation with a resultant error of measurements is established. The mark of extremely achievable accuracy of measurement of an impedance is received by a correlation method and the optimum parameters of measuring system providing its achievement are defined. The way of the accounting of a nestatsionarnost of an impedance by means of smoothing windows of the polinomialny sliding average is offered, allowing in addition to increase accuracy of measurement of an impedance.
Tomography; impedance; errors; nonstationarity; smoothing.
Идея получения пространственного распределения электропроводности биологических тканей по результатам измерения потенциалов на поверхности тела, возбуждаемых внешним источником электрического тока, весьма привлекательна, так как она является основой для создания новых, абсолютно безвредных для организма, диагностических приборов. Сравнительно низкая себестоимость в сочетании с их удобством и оперативностью потенциально делают область применения таких приборов очень широкой. Однако ряд объективных препятствий, связанных с достижением приемлемой информативности измерений, а также отсутствие алгоритмов, дающих удовлетворительные результаты реконструкции распределения
электропроводности при использовании данных, полученных при измерениях на реальном человеческом теле, сильно затрудняют осуществление этой идеи на уровне, представляющем интерес с позиции медицинской диагностики [1].
Информативность картины пространственного распределения электропроводности, так же как и информативность любого томографического изображения определяется двумя показателями - пространственным разрешением и чувствительностью к изменению картируемого параметра. В настоящее время достигнуто пространственное разрешение распределения электропроводности не более 5мм, и это только на незначительной глубине от поверхности. На сравнительно большой (100 мм и более) глубине от поверхности это значение обычно составляет не менее 20 мм. Достигнутая чувствительность к изменению значения электропроводности составляет не более 3 %. Причём, учитывая плохую обусловленность прямой задачи импедансной томографии, с увеличением пространственного разрешения картины распределения электропроводности следует ожидать снижения чувствительности и как следствие - сужения динамического диапазона изображения. По этой причине электроимпедансная томография предъявляет особенно жесткие требования к точности первичной измерительной информации.
Данная статья посвящена анализу факторов, ограничивающих точность измерения электрического импеданса биологических тканей и определению оптимальных условий съёма первичной измерительной информации.
При измерении импеданса обычно рассматривается линейная модель взаимосвязи возбуждающего тока с откликом напряжения
t
U (t ) = J к(т-1 )I (T)dr+ e(t), (1)
0
где I(t) - возбуждающий ток, U(t) - отклик напряжения, e(t) - шум, h(t) - импульсная реакция токопроводящего объекта.
В частотной области эта зависимость записывается как
U * (a>) = G(a>) • I * (И) + e* (и) , (2)
где G(ffl) - частотная зависимость импеданса, U* (ю) - Фурье (изображение сигнала отклика напряжения); I* (ю) - Фурье (изображение сигнала возбуждающего тока напряжения); e (t) - Фурье (изображение шума).
В рамках такой модели основными факторами, ограничивающими точность измерения импеданса, являются: нестационарность и нелинейность биообъекта; погрешности измерения отклика напряжения и задания тока; частотная нестабильность элементов измерительной системы; внутренний импеданс измерительной системы; наводки, шумы и помехи.
При измерении импеданса на фиксированной частоте наиболее целесообразно применение возбуждающего тока синусоидальной формы
I(t) = 10 cos(fflt) + eI (t), (3)
где eI (t) - абсолютная ошибка задания возбуждающего тока.
Сигнал отклика напряжения на действие такого тока определяется как
U(t) = U0 cos (и + ф0) + ^Ui cos ([i + 1]®i + <pi) +eu (), (4)
i=1
где eu(t) - абсолютная ошибка измерения отклика напряжения, - начальные фа-
зы соответствующих гармонических составляющих отклика напряжения.
Второе слагаемое выражения (4) описывает высшие гармоники отклика напряжения, возникающие в нём вследствие нелинейных искажений возбуждающего тока реальным проводящим объектом. Погрешность гармонической линеаризации проводящего объекта при этом определяется как
1 cos (imt + (pi-1)
Д/n =---------- ---------------------------------------------------. (5)
U0 cos (cat + p0) + I—. 4 cos (imt + р1Л)
i=2
Вследствие нестационарности электропроводящего объекта амплитуда отклика напряжения может изменяться. Наиболее просто такие изменения описываются моделью линейного тренда, имеющей следующий вид
Ui = Ui0 (1 ±в*), (6)
где Ui0 - начальная амплитуда, 6i - максимальная относительная скорость изменения импеданса на частоте i-й гармоники, т - длительность измерения.
Во время проведения измерений частота возбуждающего тока и частота тактового генератора АЦП, используемого для оцифровки сигнала отклика напряжения также могут изменяться. При этом диапазон частот измерительной системы в первом приближении определяется как
a = а0 ± Да , (7)
где ю0 - средняя частота, Дю - частотная нестабильность.
Влияние наводок, шумов и помех учитывается путём аддитивного добавления шума e(t) к сигналу отклика напряжения (4). Их можно учесть также при до-
бавлении аддитивного шума к току, но параметры шума напряжения в целом более стабильны, и их значительно проще нормировать.
Учитывая принятые допущения, корреляционная оценка действительной части импеданса определяется как
2 U
Re[G(m)] = ^т Ju(t) • I(t)dt = -0-LtJ h
sin(Дaт + p) + sin(2m + p)
+ —0
—L вт+-1- Z—.
I I T^
in ¿n*- • i
sin (im* + p)
И
Дат 2ат
e—(m) + e(m) + U 0
(8)
1 2 eI
I 0и I 0И
(И
Анализ выражения (8) позволяет выделить основные составляющие погрешности измерения импеданса. Относительная погрешность измерений
= вг + С1 -^™(ДЮТГ| + ^п(2ют) , 81 +8и +8ип +_^ (9)
Д сот ) 2ют сот 10ют
где 5/ - относительная погрешность задания тока, 5и - относительная погрешность измерения напряжения, 5цп - относительная погрешность линеаризации импеданса, се - среднеквадратическая амплитуда шумов напряжения, приходящаяся на единицу импеданса. Второе слагаемое выражения (9) может быть сведено к нулю выбором продолжительности измерения т = л1 ю, погрешность линеаризации таким способом устранить нельзя, так как начальные фазы высших гармоник требуют для компенсации разных т.
При рабочей частоте ю = 50 кГц и типовой кратковременной нестабильности частоте кварцевого генератора 10-10 второго слагаемое (9) имеет порядок 10-7 для продолжительности измерений более 100 секунд. Следовательно, погрешность, связанная с нестабильностью частоты очень незначительна.
При 0 = 0,05 с-1, вычисленной по типовой реограмме, для обеспечения приемлемой чувствительности время измерения не должно превышать 0,1 с. В противном случае погрешность, связанная с нестационарностью, не позволит регистрировать пульсовые колебания импеданса тела. При использовании измеренных данных для реконструкции пространственного распределения электропроводности это время должно быть, по крайней мере, на 2 порядка меньше, с тем, чтобы обеспечить приемлемый динамический диапазон томографического изображения.
Таким образом, при частоте 50 КГц продолжительность измерений должна составлять около 1 мс, т.е. включать в себя 50 периодов зондирующего тока. В данном случае для оцифровки сигналов будет достаточно использовать 8-битные ЦАП-АЦП, обеспечивающие относительную погрешность измерения импеданса (четвёртое слагаемое в 1) на уровне 2,5 10-5.
Для обеспечения требуемого уровня точности погрешность линеаризации должна быть на уровне 0,4 %. Согласно данным [2], при использовании для съёма потенциала отдельных неполяризующихся электродов это требование выполняется при амплитудах отклика напряжения не более 30 мВ.
При базовом импедансе тела 100 Ом амплитуда зондирующего тока должна составлять не более 300 мкА, а среднеквадратическая амплитуда шумов не должна превышать 150 мкВ. В этом случае результирующая относительная погрешность измерения импеданса будет на уровне 0,01 %.
Таким образом, в результате проведённого анализа основных составляющих погрешности измерения импеданса установлено, что основным фактором, ограничивающим точность измерения импеданса является нестационарность биообъекта. Сокращение продолжительности измерений позволяет снизить обусловленную ей погрешность измерения импеданса до 0,01 %, при условии обеспечения очень эффективного шумоподавления. Дальнейшего повышения точности измерений можно достичь повышая частоту зондирующего тока. Так на частоте 1МГ ц, при длительности измерений 0,2 мс, разрядности ЦАП-АЦП 12 бит и амплитуде зондирующего тока 20-30 мкА, относительная погрешность измерения импеданса снизится до уровня 0,002 %, при максимально допустимом уровне посторонних шумов напряжения на входе измерительной системы не более 500 мкВ.
Дальнейшее увеличение частоты зондирующего тока приводит к дополнительным сложностям, связанным с усилением излучения и как следствие - к снижению точности задания возбуждающего тока, а так же с усилением влияния конфигурации измерительных цепей на результаты измерений. По этой причине измерения биоимпеданса на частотах более 1 МГц обычно не проводят [3].
Очевидно, что дальнейшего повышения точности измерения импеданса возможно только при явном учёте нестационарности. Проще всего это достигается посредством оконного сглаживания спектра. Корреляционная оценка импеданса в этом случае определяется как
О(ю) = ^ ¡¡к (( -£) ■и(()[собЮ) + ] Бт(ю()й(йТ, (10)
0 0 0
где К(() - оконная функция, Т - ширина окна, Q - нормирующий коеффициент.
В качестве оконной функции наиболее целесообразно использовать полиномиальные скользящие средние, которые строятся из условия наилучшей аппроксимации локального участка тренда полиномом заданной степени.
Весовые коэффициенты этого окна определяются первой строкой матрицы
К = (1-1 ■ .т), (11)
к
где 1(и, V) = ^ ги+”-2, 3(ы, V) = (у - к +1)“-1, и - номер строки, V - номер столбца
г=-к
(число строк/ столбцов квадратной матрицы здесь равно ширине сглаживающего окна), к - степень сглаживающего полинома.
На рис. 1 представлены результаты сглаживания ступенчатого тренда (наи-хутший вариант, при котором в случае обычного усреднения, на некоторых участках погрешность может достигать 100 % от величины скачка) полиномиальными скользящими средними разной степени.
Рис. 1. Результаты сглаживания ступенчатого тренда полиномиальными скользящими средними разной степени
Из рис. 1 следует, что даже в случае использования модели линейного тренда (к = 1) ошибка, связанная с нестационарностью уменьшается в 4 раза. С увеличением степени сглаживающего полинома погрешность, связанная с нестационарно-стью уменьшается всё сильнее. Так при к = 15 она составляет менее 10 %.
При использовании сглаживающих окон, однако, необходимо учитывать некоторые особенности. При вычислении ошибки, связанной с нестационарностью проводящего объекта здесь не применима формула (9), постулирующая линейный тренд, полностью устраняемый при к=1. С увеличением порядка полинома шумоподавление при корреляционной оценке по (10) несколько уменьшается, по сравнению с оценками без применения сглаживающего окна. Для учёта этой особенности, вместо длительности измерения т в формуле (9) следует использовать эффективную ширину сглаживающего окна, которая уменьшается с увеличением порядка полинома. Кроме того, для полиномов большой степени и большой ширине сглаживающего окна матрица I, используемая для вычисления весовых коеффици-ентов, становится плохо обусловленной, что приводит к большим ошибкам в оценках весовых коэффициентов и падению эффективности сглаживающего окна. Кроме того, необходимо учитывать, что полиномное сглаживающее окно не является казуальной системой, так как использует в своей работе будущие отсчёты сигнала. Это приводит к неизбежной задержке измерения параметров импеданса относительно реального времени, увеличивающейся с расширением окна.
Указанные особенности ограничивают порядок как сглаживающего полинома, так и ширину сглаживающего окна, и как следствие затрудняют адекватный учёт нестационарностей сложной формы. Тем не менее по ориентировочным оценкам применение сглаживающих окон позволяет снизить погрешность измерения импеданса, связанную с его нестационарностью ниже 0,001 %. В этих условиях на первый план выступают погрешности, связанные с нелинейностью импеданса и снижение их негативного влияния становится основной задачей.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Бамбер Дж., Тристам М., Лич М., Джонс К., Уэбб С, Хилл К., Бентли Р. Физика визуализации изображений в медицине: В 2-х томах. Т. 2: Пер. с англ. / Под ред. С. Уэбба. - М.: Мир, 1991.
2. Киреев А.В. О нелинейной природе импеданса биологических тканей // Известия ЮФУ. Технические науки. - 2010. - № 8 (109). - C. 115.
3. Мартиросов Э.Г., Николаев Д.В., Руднев С.Г. Технологии и методы определения состава тела человека. - М.: Наука, 2006. - C. 68.
Статью рекомендовал к опубликованию д.м.н., профессор А.И. Сафронов.
Киреев Андрей Владимирович - ГОУ ВПО «Пензенская государственная технологическая академия»; e-mail: kireewska@mail.ru; 440039, г. Пенза, пр. Байдукова, ул. Гагарина, 1а/11; тел.: 88412496155; кафедра информационных технологий и менеджмента в медицинских и биотехнических системах; к.т.н.; доцент.
Kireev Andrey Vladimirovich - State Educational Institution of the Higher Vocational Training "Penza State Technological Academy”; e-mail: kireewska@mail.ru; 1a/11, Jour.
Baydukova/Gagarina, Penza, 440039, Russia; phone: +78412496155; the department of information technology and management in medical and biotechnical systems; cand. of eng. sc.; associate professor.
УДК 51-76: 612.88
А.В. Киреев, Р.Ф. Акжигитов
МОДЕЛИРОВАНИЕ ОСНОВНЫХ ЗАКОНОМЕРНОСТЕЙ ФОРМИРОВАНИЯ СТАБИЛОГРАФИЧЕСКИХ СИГНАЛОВ В РАМКАХ КОНЦЕПЦИИ ПЕРЕВЁРНУТОГО МАЯТНИКА
Описана модель формирования стабилографического сигнала основанная на представлении системы равновесия в виде композиции линейной динамической и нелинейной статической подсистем. В процессе поддержания равновесия первая из них выполняет функцию грубого управления и описывает плавные перемещения проекции центра масс большой величины. Вторая подсистема выполняет функцию тонкого управления и описывает малые хаотические отклонения проекции центра масс от точки покоя. В рамках разработанной модели выполнен анализ устойчивости системы равновесия по Ляпунову, определены границы устойчивости для свободных параметров модели, намечены пути их интерпретации и предложен обобщенный показатель запаса устойчивости системы равновесия.
Стабилограмма; тренд; линеаризация; устойчивость; авторегрессия.
A.V. Kireev, R.F. Akzhigitov
SIMULATION OF THE MAIN REGULARITIES OF FORMATION STABILOGRAM WITHIN THE CONCEPT OF THE OVERTURNED
PENDULUM
The model of formation of a stabilogram based on representation of system of stability in the form of composition of the linear dynamic and non-linear static subsystems is described. In the course of stability maintenance the first of them executes function of rough control and describes the smooth relocation of a projection of center of masses of big value. The second subsystem executes function of thin control and describes small chaotic deviations of a projection of center of masses from a rest point. Within the developed model it is analysed stability of system of equilibrium according to Lyapunov, stability boundaries for the free model parameters are defined, ways of their interpretation are scheduled and the generalized index of an inventory of stability of system of equilibrium is offered.
Stabilogram; trend; linearization; stability; autoregression.
Стабилометрия является одним из наиболее молодых методов функциональной диагностики, активно использующимся не более 20 последних лет. Благодаря достаточно высокой чувствительности параметров этого метода к состоянию большого числа функциональных систем организма, при абсолютной безвредности для организма и простоте использования, в последнее время стабилометрия привлекает к себе всё больше внимания со стороны различных исследователей.
