CC BY
63
Радіоелектроніка біомедичних технологій
РАДІОЕЛЕКТРОНІКА БІОМЕДИЧНИХ ТЕХНОЛОГІЙ
УДК 621.372.061
СРАВНЕНИЕ КЛАССИЧЕСКОГО МЕТОДА РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧИ ИМПЕДАНСНОЙ ТОМОГРАФИИ С МЕТОДОМ «ЗОН» ПРОВОДИМОСТИ
Сушко И.А., аспирантка; Рыбин А.И., д.т.н. профессор
Национальный технический университет Украины «Киевский политехнический институт», г. Киев, Украина
Введение
Методы импедансной томографии находят всё большее применение в медицине, технике, физике, геодезии [1 - 5] и т.д. Методы привлекательны тем, что в отличие от других («линейных» — рентгеновских, эмиссионных, ядерного магнитного резонанса и т.д.) методов томографии технические средства для реализации этих методов характеризуются относительной дешевизной и простотой реализации. Однако математические методы решения задачи реконструкции образа по проекциям отличаются (в сравнении с «линейными») значительно большей сложностью и громоздкостью [6 - 18]. Это обусловлено тем, что траектории сигнала носителя информации (линии равного напряжения или равного тока) в импедансной томографии являются нелинейными даже для однородного объекта исследования [13] и их форма сильно усложняется в случае наличия неоднородностей. Поэтому отдельной и важнейшей задачей импедансной томографии является реконструкция (обратная задача или задача синтеза) образа по проекциям. При этом считается известным начальное приближение, т.е. распределение поверхностных проводимостей областей неоднородностей внутри томографического сечения (фантома). Также, в частности, начальное приближение фантома может быть и равномерным. Для решения обратной задачи используют два метода: метод обратной проекции и итерационный метод вычислений.
Метод обратной проекции отличается простотой и высокой вероятностью получения неверного результата (в зависимости от картины неоднородностей, появившихся в начальном фантоме). При использовании этого метода проекции (измеренные передаточные сопротивления) проецируются на фантом в секторы между линиями равного напряжения (вычисленными для начального приближения), начинающимися на измерительных электродах. Метод обратной проекции является «калькой» широко распространённых методов «линейной» томографии (рентгеновской, эмиссионной, ядерного магнитного резонанса, ультразвуковой [13, 16 - 22 ] и т.д.). При наличии одной дополнительной (к фантому первого приближения)
166
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
неоднородности метод используется, например, при диагностике легких [3,
4].
При реконструкции образа итерационными методами необходимо решить две задачи:
1. Задачу анализа или вычисление по заданному фантому передаточных сопротивлений по его обводу (по начальному приближению распределения поверхностных сопротивлений областей неоднородности).
2. Задачу синтеза или реконструкции картины (образа) внутреннего распределения неоднородностей и их поверхностных сопротивлений, при которой результаты анализа и результаты измерений с заданной точностью совпадают.
Задача реконструкции является весьма трудоёмкой и громоздкой, включающей, к тому же, на каждой итерации решение задачи анализа. Для решения как прямой, так и обратной задачи необходимо представить фантом в виде его математической модели. Для этого томографическое сечение с непрерывным распределением поверхностной проводимости следует в соответствии с методом конечных элементов [23] представить в виде совокупности конечных элементов (для определённости — в виде квадратов), представленных их эквивалентными схемами (в виде квадрата [14] сопротивлений с двумя диагоналями). Как прямую, так и обратную задачи импедансной томографии целесообразно решать методами теории цепей [15 - 18].
Основные теоретические положения
Решение задачи реконструкции фантома по проекциям в итерационной процедуре, например, методом Ньютона-Рафсона заключается в составлении и решении системы уравнений вида
Jdo2 ... д ^1 пер .зыч /5%1 ГД ad
dz
1 пер .зыч
,/да1 dz
1 пер.зыч
dz■ dz
2 пер.зы1-:
Jdai dz-
2 пер.зы*-:
,/do2
dzn
пер,
Jdab
IV пер.зыч
Jda^ dz
IV пер ,зыч
Jdiт2
dz
N пер.зыч
Jda
ivJ
= пер выч /dak]x[Aak] = -[ДZJ, (1)
где dzj пер ,3M~/dad — производная от у-го передаточного сопротивления,
вычисленного по фантому по поверхностной проводимости d-ro конечного элемента; Аак — искомое корректирующее значение приращения поверхностной проводимости -ГО конечного элемента; ^тпер.изм. — Z І7Ї пер ,зыч, = Дzm — разница между измеренным и вычисленным по фантому значениями m-го передаточного сопротивления; N — количество конечных эле-
х
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
167
Радіоелектроніка біомедичних технологій
ментов в фантоме.
Таким образом, для решения задачи реконструкции необходимо многократно решать задачу анализа, т.е. задачу вычисления обратной матрицы (на каждой итерации для текучих значений поверхностных проводимостей неоднородностей фантома), по значениям которой находить передаточные
сопротивления пер.зыч.. По вычисленным передаточным сопротивлениям далее следует находить разности для оценки соответствия картины распределения неоднородности внутри фантома истинному распределению в исследуемом сечении. При неудовлетворительно большой невязке, оцениваемой по [Д7т] (1), процедуру вычислений по (1) следует по-
вторить.
Если задача анализа может быть решена достаточно просто, учитывая разреженность матрицы проводимостей системы уравнений узловых напряжений фантома методом модификаций [24 —26], то решение системы (1) представляет значительную проблему.
Так, количество конечных элементов N должно быть достаточно большим, чтобы полученное в результате реконструкции изображение визуально воспринималось. Пусть N = 1000. Тогда порядок матрицы [5Zm пер зыч /также равен 1000, т.е. матрица содержит 10б элементов,
и эта матрица, в общем случае, не содержит нулевых элементов. Её структура также априори неизвестна. При наличии обратной матрицы текущего фантома задача вычисления всех 10б производных dZj пер,зыу/d<jd не представляет трудности [27 - 32]. Однако, при этом матрица [&Zт пер ЗЬ[Т,/5су,]
порядка 1000 всегда будет плохо обусловленной. Это связано с тем, что производные от передаточных сопротивлений по поверхностной проводимости для двух соседних конечных элементов, расположенных в центре фантома, будут весьма малыми (в сравнении со значением передаточной проводимости) и достаточно близкими по величине. Последнее является достаточным условием равенства нулю определителя этой матрицы. Так, например, для фантома рис.1 при нормированных значениях поверхностной проводимости каждого конечного элемента, равных 1, были вычислены передаточные сопротивления от тока источника между узлами, расположенными по горизонтали (рис. 1... .рис.4) и напряжения на измерительных электродах, расположенных в основании равных дуг длиной 2тт/16 (точки по обводу контура, пронумерованные от 0 (крайняя точка слева) по часовой стрелке до 14, причём правый крайний узел не пронумерован), а также соответствующие производные.
168
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
Результаты машинного эксперимента
Рис.3 Рис.4
Для равномерного фантома рис.1 значения передаточных сопротивле-
нии
=0.42526, = 0.67047, =0.71303,
Z;L"rJ=0.52703,
=0.52703, ZT^=0.57956, Z'f^-
=0.71303,
=0.67047,
(nep]_
(пер]
'll
=0.79461,
=0.62754,
(nep]_
zt:
пер]
=0.62754, = 0.79461, =0.57956,
Zf^=0.42526.
(2)
Производные dZ'f^ jdak от передаточных сопротивлений по поверхностным проводимостям и1га2 конечных элементов, расположенных в центре фантома (показаных на рис.1) соответственно равны
(пер]
/дах = -2.9624 ■ ю“е, dZ'f^/da, = -1.6187- ю“%
(пер]
/да-у = -4.S33S ■ іо_% jba-y = -7.166 ■ іо“% jddy = -8.6158 ■ 10“%
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
169
Радіоелектроніка біомедичних технологій
3Z,
(пер]
/даг = -8.0258 ■ іо % dZ'^0ep}/da1 = -7.1906 ■ ю
dZ'f/^/da1 = -5.8206 ■ ю“% dZ'^^/da, = -3.9089 ■ io“%
(пер ]
dZ'f^/da, = -2,0821 ■ 10
13 / л-’-' ,
(пеР] /д^ _ .„лг Д7'СпеР]
dZ^'/дъ = -4.5114 ■ 10 е
dZ.
]/да2 = -2.2325 ■ 10“% dZf*P]/да2 = -1.2933 ■ 10“%
/гтаи 1
(пер ]
/да2 = /да2 =
-4.2759- 10'
-6.8085 ■ 10'
/да2 = -8.4817 ■ 10'
дҐ,
пер]
/да2 = -7.7147- 10“% дҐ///р]/да2 = -6.7114 ■
10'
dZf^/da2 = -5.
1998 ■ 10'
flZj“ep]/3fj2 = -3.305 ■ 10'
зр]
dZ.
(пер]
'14 /$а2 = -3.5462- 10 % (3)
Аналогично, для фантома рис. 2 («маленькая» и «средняя» неоднородности имеют поверхностные проводимости <7 = 100)
Zfep]= 0.42274, Zfep] = 0.51382, Z'fSp'i= 0.55304, Zf&P']= 0.58477, Z:nep]= 0.61348, ^ПЄр:)= 0.64512 , Z'fepi= 0.71527, Z'^P]= 0.70581,
Z'fsp} = 0.62219, Zj4™1^ 0.5844, Zj%ep3= 0.57388, Zj™p3= 0.55679,
Т-І^пер] T-1 пер I
Z13 - 0.51494, Zj44 - 0.42276;
(neP]_
5ZfepJ/5<71=-2.7344 ■ 10“% aZ;JIiepJ/5ff1= -4 (пер] /л^ _ r ,^-5 я 7(neP]
5Z%lctjJ/5<7: =-5.881 ■ 10“% 5Z4nepJ/3^i= -6,7723 ■ 10
= -4.9438 ■ 10
-9
-9
dZ.
(пер]
/3<Т1=-7.5773 ■ 10“% 5Z^nep]/5<7i= -8.0129 ■ 10
dzf^/da^-7.6976 . 10“% azf^/^i = -4.5936 . 10“\ dZf^/da^- 6.7293 ■ io“% dZ^^/да^= -3.9477 ■ io“% 5ZJ™PJ/3(7!= 5.5676 ■ 10“\ 3Zj^eP V5orl= 6.2508 ■ 10“% 5Zj3ep^/5<j1= 6.1067 ■ io“% dZf^ / даг= 3.0936 ■ ю“%
dZ.
(пер]
/da2= -1.6975 ■ 10“% dZ'fsp^f da2= -2.9141 ■ 10
1-9
dZ'lnep]/3<T2= -3.3527 ■ 10'
r(nep]
3Z4nep]/3<T2= -3.8699 -10“%
3Z21cijJ/3<72= -4,4562 -10“% дҐ^РІ/da2= -4.8276 -10“% 3Z^nep]/З<72=—4.6718 - 10“% 3Zgnep-)/3<J2 = -2.6817- 10 dz\nep]/3<72= -4.4256 ■ 10“% dZj'^/d^ -2.7921 ■ 10 3Zj“epV3(72=з.5892 ■ 10“% 3Zj%ep3/5^2=4'0381 ■10
-9
^-5
dZ'f*P}/да2=3.953' 10“% 3Zj4epj/5t72=1.9947 ■ 10
— 5
(4)
170
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
Передаточные сопротивления и их производные по поверхностной проводимости для тех же конечных элементов для фантома рис. 3 («средняя» неоднородность с поверхностной проводимостью а = 0) приобретают следующие значения
Z'fsp}= 0.50917, Z'fSp} = 0.68079, Z'fep}= 0.7649, Z'f*P}= 0.8273, Z'^ep}= 0.86768, Z^ep]= 0.89494, Z'f^= 0.93088, Z'^ep}= 0.65597,
Z^ =0.3026, Z^ep]= 0.037939, Z^ep]=-0.11344, Zjy6P:)=-0.035354, Zl^ep:)=0.16983, Z:n/P:)=0.30559;
13 ’14 ’
'(пеР]/д^_ Д7'СпеР]
IdZ'fSpi/д^=-0.015025, dZ'fSpJ/3c71=-0.029793,
5Z^nep]/5^1=-0.038952, dZ'fep]/3(Т1=-0.04651, 5Z:nep-)/5c71=-0.05i443, 5 Z^-’/Wj =-0.054209,
5ZrP]/3^=—0.055083, 5Z^nepY5(Ti =-0.03451, ,
5Z^nep V5t7i= 0.0056016, 5Zj^0p V5iJi= 0.060454,
^Zjnep]/5^= 0 1076^ З7'^ер^5і71= 0.10208,
3Z^ep V5ori= 0.06039, 5Zj"ep Vй^1= 0.020249; gZ',nep]/3(J2=_0.0050922, SZ^ Y3^2=-0.010642, dZ'fepyda2=-0.014834, 5Z^nepY5cj2=-o.019211, dZ'fep} fda2=-0.023183, £^пер:)/д^2=_0'026863’
5Z^nep3/3(72=-0-032573, dZ'^Pi/да3=-0.049803,
5ZgCiep')/5^2=—0.046096, 5Z^ep“)/5c72=—0.0 ІЗ 965,
3z;7p]/0£T2= 0.026692, 3Z;^epVд£Т2= 0.036697,
5Zjnep]/0 022583j 52j-"ep:)/0£72= 0.0071607. (5)
Теперь рассмотрим фантом рис. 4, на котором произведено разбиение на 14 зон проводимостей [15 - 18, 32, 33]. При этом полагаем, что каждая такая k-я зона имеет постоянную поверхностную проводимость ак. Пусть теперь зона, примыкающая сверху к «экватору» фантома рис.4, имеет проводимость сг2, а зона, находящаяся на самом верху на рис.4 - а2. Тогда передаточные сопротивления для равномерного фантома (£72= £Т2 =...= £714), вычисленные на тех же, что и рис.1...рис.3 измерительных электродах, имеют значения
.(пер}_
Z'fvy=0.42526, Z™» =0.52703, Z'^V) =0.57956, Z^4 =0.62754, Z^.nep]= 0 67047, Z^nep:)=0.71303, Z'fep}=0.79461, Z'^p]= 0.79461,
(пері
(пер].
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
171
Радіоелектроніка біомедичних технологій
=0.71303,
Z^O.67047, Z=_rj=0.62754,
'(пєр]_г '11
Z'f*91=0.57956,
1 z 5
Zj3 н—0.52703, Zj4 0.42526,
а производные ПО поверхностным проводимостям ЗОН tJj, <72 значения
1 ^ ,=-0.13998. 0Z2 ^/ОаЛ =-0.14394,
следующие
5z;:nep:)/3^
!= -0.14694, 3Z^nep^/3^1= -0.14889, 5Z;;nepJ/5^= -0.15124, 5ZgnepY5(Ti =-0.15384, 5Z;nep:)/5*Ti= - 0.161, dzf*p}/даг= -0.21398, dZf^/да1= -0.18833, 3Z^0pV5ori= -0.17719, 3Zjyep^/даг= —0.16738, 3Zj'™P75^= -0.15784, ^Zjnep]/5c7i=_0 14863^ 3Z^ep^/5^= -0.132; 5Zjnepj/5^2= 0.51143, dZ'f^]/da2= 0.71444, dzfsp}/da2= 1.1318, dZ^ep}/da2= 0.77858,
dZ:nep7d<T;
dZ^/da,
dzf^/da,
dZ'f/^fda,
dzf*p]fda2
= 0.70978, 2= 0.66234, = 0.60935, = 0.58005, = 0.52415,
(6)
dZfpi/да2= 0.68691, dZ(nep ^ jQa = 0.62536,
dZ^0 p-}/da2= 0.59666, dZ'/Tp] f da2= 0.55643,
dZf*p} /da2= 0.4442.
Производные от передаточных сопротивлений в (3) имеют величину порядка 6-10-4 %... 1-10-4 % от абсолютных значений передаточных сопротивлений. Аналогично малые значения производных получаются для фантома рис.1 и в случае, когда дифференцирование проводится по поверхностным проводимостям конечных элементов, расположенных на обводе контура в крайней верхней части фантома. Обозначим эти поверхностные проводимости как <т3=1 и ст4= 1. Тогда производные
dZ'fep'^/да3= 3.0853-10 6, dZ^iepJ/du3= 1,1705 ■ 10"
(пер]
3Z.
Дер]
/да3= 4.0073 ■ ю 5 ит.д.;
(пер]
щ
/да±= 3.6959- ю_6, dZfspi/да4= 1,3501- ю"
dZ'^/d
аА = 4.1117 ■ 10 3 и т.д.
(пер]
При этом разница между производными 3Z^nepJ/да1 и dZ'Z1^/да2 в (3)
не превышает 10_б, а для 3Z>nepJ/да3 и 3Z£nepj/3:т5 — 10_J.
Ещё худший результат получен в (4). И только при очень больших изме-
172
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
нениях поверхностных проводимостей (<т=0) в (5) значения производных становятся соизмеримыми со значениями передаточных сопротивлений.
Для результатов (3) и (4) значения элементов обратной матрицы в (1)
[dZm ПЄр ЗЬ[„ /йсгц.] будут порядка 10б для (3) и 109 для (4), что не будет
способствовать сходимости итерационной процедуры реконструкции образа по проекциям.
Отметим также, что фантом рис.1 содержит 840 узлов и 776 конечных элементов. И если решение прямой задачи методом модификаций (учитывая разреженность матрицы проводимости системы уравнений узловых напряжений) не представляет сложности, то как составлении, так и обращение матрицы производных в (1) представляет собой громоздкую задачу, решение которой весьма проблематично, учитывая плохую обусловленность матрицы производных и её порядок. Тем более, что структура матрицы производных \dZm П0р ЗЬ[1.; /даk ] неизвестна и матрицу следует считать полностью заполненной ненулевыми элементами.
Приведенные в (6) результаты показывают, что, учитывая плохую обусловленность матрицы производных от передаточных сопротивлений по поверхностным проводимостям конечных элементов, значительно более привлекательным является метод зон проводимости, в котором такими конечными элементами являются сами зоны, состоящие из большого числа элементарных конечных элементов.
а б
Рис.5
При введении зон проводимостей производится дополнительная структуризация схемы фантома, в результате чего матрица производных 776-го порядка преобразуется в 16 (для разных узлов включения источника тока по кросс-методу) матриц производных 14-го порядка.
При этом растр для результирующего изображение для зон проводимости (в отличие от растра рис. 1 классического метода конечных элементов) будет, например, иметь вид рис.5, а,б.
Вісник Національного технічного університету України "КПІ" 173
Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
Производные fdGj в этом случае соизмеримы со значениями пе-
редаточных сопротивлений, а матрицы таких производных достаточно хорошо обусловлены.
Выводы
Проведенный сравнительный анализ классического метода реконструкции образа по проекциям и метода зон проводимостей показал:
1. Громоздкость процедуры формирования матрицы производных для классического подхода, при котором такое формирование следует выполнять на каждой итерации. Эта громоздкость определяется большими порядками матрицы производных.
2. Отсутствие структуризации и плохую обусловленность матрицы производных (при классическом подходе), вызванную весьма малой (в общем случае) чувствительностью передаточных сопротивлений к изменению поверхностной проводимости конечных элементов.
3. Простоту формирования матрицы производных при использовании метода зон проводимости, поскольку её порядок в рассмотренном случае рамен 14, а при использовании 32 электродов — 30.
4. Значительно лучшую обусловленность матрицы производных при методе зон проводимости в сравнении с классическим методом.
5. Хорошую геометрическую разрешающую способность метода зон проводимости, обеспечиваемую растром, образованным наложением поверхностных проводимостей зон, полученных для каждой отдельно схемы подключения источника тока.
Литература
1. Brown B.H. Electrical Impedance Tomography/ B.H. Brown, D.C. Barber // Clinical Physics and Physiological Measurement.—1992.—v.13.—Sappl. A, 207p.
2. Электроимпедансная томография / Я.С. Пеккер, К.С. Бразовский, В.Ю. Усов, М.П. Плотников, О.С. Уманский.— Томск: ООО «Издательство научно-технической литературы»,2004.—190с.
3. Корженевский А. В. Квазистатическая электромагнитная томография для биомедицины. Автореферат диссертации на соискание ученой степени д.ф.-м.н. Москва, 2009 —32c.
4. Корженевский А. В. Способ получения томографического изображения тела и электроимпедансный томограф / А.В. Корженевский, Ю.С. Культисов, В.А. Черепнин. Патент РФ на изобретение № 2127075, 1996 .
5. Дорожковець М.М. Фізичні принципи побудови вимірювальних томографічних систем промислового застосування / М.М. Дорожковець // Вимірювальна техніка та метрологія.— 1996.— №52.— С.113 —121.
6. Cormack A.M. Early two-dimensional reconstruction and recent topics stemming from it/ A.M. Cormack // Nobel Lectures in Physiology or Medicine 1971—1980. World Scientific Publishing Co - 1992. - P.551 - 563.
7. Murray T. Electrical Impedance Computed Tomography Based on a Finite Elements Model / T. Murray, Y. Kagawa // IEEE Trans. on Biomed. Eng. - 1985. - v. 32. -P.177 - 184.
174
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
8. Yorkey T. A comparison of impedance tomographic reconstruction algorithms / T.J. Yorkey, J.G. Webster // Clin. Phys. Physiol. Meas . - 1987. - v. 8. - Suppl. A. - P. 55.
9. Cheney M. Electrical Impedance Tomography / M. Cheney, D. Isaakson, J.C. Newell // SIAM Reviev. - 1999. - v .41 — № 1 . - P. 85 - 101.
10. Hua P. et al. Finite Element Modeling of Electrode - Skin Contact Impedance Tomography // IEEE Trans. on Biomed. Eng.— 1993. - v. 40. - P. 335 - 343.
11. Рыбин А.И. Решение обратной задачи импедансной томографии методом модификаций / А.И. Рыбин, И.А. Манюк // Радиоэлектроника.- 1998.- № 4. - С. 36 - 44 (Изв. вузов).
12. Рибіна І.О. Моделювання кінцевого елемента для об’ємного фантома в імпе-дансній томографії / І.О. Рибіна // Вісник НТУУ «КПІ». Сер. Радіотехніка. Радіоапара-тобудування.—2011.— №46.—С.72 —84.
13. Физика визуализации изображений в медицине. Под ред. С. Уэбба.—М.: Мир,1991, т.1.—408с., т.2— 408с.
14. Рибіна І.О. Моделювання кінцевого елемента в імпедансній томографії / О.І. Рибіна, Є.В. Гайдаєнко // Вісник НТУУ «КПІ». Сер. Радіотехніка. Радіоапаратобуду-вання.—2010.—№41.—С .19 —24.
15. Рыбина И.А. Решение прямой задачи импедансной томографии методами теории цепей / И.А. Рыбина // Вісник НТУУ «КПІ». Сер. Радіотехніка. Радіоапаратобуду-вання.—2010.—№43.—С 4.—13.
16. Рибіна І.О. Метод променів провідностей та моделювання фантома в імпедансній томографії / І.О. Рибіна // Вісник ЖДТУ.—2010.—т.8.—4.—С.21 — 28.
17. Рибіна І.О. Обчислення похідних від передаточного опору по поверхневій провідності кінцевих елементів при розв’язанні зворотної задачі імпедансної томографії / І.О. Рибіна, О.І. Рибін, О.Б. Шарпан // Вісник НТУУ «КПІ». Сер. Радіотехніка. Радіоапаратобудування. - 2011.—№44.— С.5—11.
18. Рибіна І.О. Розв’язання зворотної задачі імпедансної томографії методами зон провідностей та зворотної проекції / І.О. Рибіна, О.І. Рибін, О.Б. Шарпан // Вісник НТУУ «КПІ». Сер. Радіотехніка. Радіоапаратобудування. - 2011.—№45.—С.5—18.
19. Hounsfield G.N. Computed Medical Imaging / G.N. Hounsfield // Nobel Lectures in Physiology or Medicine 1971—1980. World Scientific Publishing Co. -1992. - P.568 - 586 .
20. Mansfield P. Snap-short MRI / P. Mansfield // Les Prix Nobel. The Nobel Prizes 2003. Nobel Foundation. - 2004. - P. 245 - 257.
21. Хорнак Дж. П. Основы МРТ.
22. Мэнсфилд П. Быстрая магнитно-резонансная томография/П. Мэнсфилд // Успехи физических наук. - 2005. - т.175.— № 10. - С. 1044 - 1052.
23. Сильвестр П., Феррари З. Метод конечных элементов для радиоинженеров и инженеров-электриков.— М.: Мир, 1986.—229с.
24. Рыбин А.И. Решение задач моделирования обращением матрицы методом взаимных производных/ А.И. Рыбин // Радиоэлектроника.—1978.—№6.—С.35—47. (Изв. вузов).
25. Рыбин А.И. Численно-символьный метод анализа электрических цепей обобщенным методом модификаций / А.И. Рыбин // Праці Інституту електродинаміки НАН України: Сб. наукових праць.—2002.—№1(1).—С.26—30.
26. Основи теорії кіл: Підручник для ВНЗ. Ч.2 /Ю.О. Коваль, Л.В. Гринченко, І.О. Милютченко, О.І. Рибін / За заг. ред. В.М. Шокала та В.І. Правди.—Харків: ХНУРЕ: Колегіум, 2006.—668с.
27. Трохименко ЯК. Статистический анализ линейных электронных цепей постоянного тока / Я.К.Трохименко, И.С. Каширский, А.И. Рыбин // Радиоэлектроника.—
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
175
Радіоелектроніка біомедичних технологій
1974. —№6.—С.69—73. (Изв.вузов).
28. Трохименко Я.К. Вероятностный аналіз линейных электронных цепей переменного тока / Я.К. Трохименко, И.С. Каширский , А.И. Рыбин // Радиоэлектроника.—
1975. —№6.—С.35—40. (Изв. вузов).
29. Трохименко Я.К. Вероятностный аналіз линейных электронных цепей методом разбиения на подсхемы / Я.К. Трохименко, И.С. Каширский, А.И. Рыбин // Радиоэлек-троника.—1976.—№6.—С.70—77. (Изв. вузов).
30. Рыбин А.И., Трохименко Я.К. Символьный аналіз электронных цепей с использованием матрицы взаимных производных / А.И. Рыбин, Я.К. Трохименко // Радиоэлектроника.—1979.—№6.—С.32—37. (Изв. вузов).
31. Rybin A.I., Trochimenko J.K., Rodionova M.V. Analysis of Irregular Networks using Diacoptic Modification of Admittance Matrix / A.I. Rybin, J.K. Trochimenko, M.V. Rodionova // Electronic Engeneering (Czech rep.).—1993.—№4.— P.9 —12.
32. Рибіна І.О. Обчислення похідних від передаточного опору по поверхневій провідності кінцевих елементів при розв’язанні зворотної задачі імпедансної томографії методом зон провідностей / І.О. Рибіна, О.І. Рибін, О.Б. Шарпан // Вісник НТУУ «КПІ». Сер. Радіотехніка. Радіоапаратобудування. — 2011. — №44.— С.5—21.
33. Сушко І.О. Алгоритм розв’язання прямої задачі імпедансної томографії методом модифікацій / І.О. Сушко // Вісник НТУУ «КПІ» Сер. Радіотехніка. Радіоапаратобудування. - 2011. - № 47. - С.165 - 175.
Сушко И.А., Рыбин А.И. Сравнение классического метода решения обратной задачи импедансной томографии с методом «зон» проводимости. Проведено сравнение свойств классического итерационного метода реконструкции образа в импедансной томографии и метода зон проводимости с точки зрения организации вычислительных процедур, ожидаемой точности и трудоёмкости получения образа. На расчётных примерах показана весьма малая чувствительность передаточных сопротивлений к изменению поверхностной проводимости отдельного конечного элемента. При этом малая чувствительность проявляется, как для конечных элементов, расположенных в центре фантома, так и для конечных элементов, расположенных по обводу фантома. Малая чувствительность вызывает плохую обусловленность матрицы производных, усугубляемую большим порядком матрицы. Всё это приводит к численной неустойчивости при решении обратной задачи (задачи реконструкции). Метод зон проводимости позволяет обойти указанные трудности: приводит (за счёт дополнительной структуризации фантома) к радикальному уменьшению порядка матрицы производных, повышению чувствительности и обусловленности такой матрицы. Для оценок приведены значения передаточных сопротивлений и производных для равномерного фантома и фантомов с неоднородностями.
Ключевые слова: импедансная томография, фантом, конечный элемент, зоны проводимости, прямая задача, обратная задача, задача реконструкции изображения, матрица производных.
Сушко І.О., Рибін О.І. Порівняння класичного методу розв’язання зворотної задачі імпедансної томографії з методом «зон» провідності. Проведено порівняння властивостей класичного ітераційного методу реконструкції образу в імпедансній томографії і методу зон провідності з точки зору організації процедур обчислення, очікуваної точності та трудоємності отримання образу. На прикладах розрахунку показана доволі мала чутливість передаточних опорів до зміни поверхневої провідності окремого кінцевого елементу. При цьому мала чутливість проявляється, як для кін-
176
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
Радіоелектроніка біомедичних технологій
цевих елементів, розташованих в центрі фантому, так і для кінцевих елементів, розташованих по обводу фантому. Мала чутливість є причиною поганої обумовленості матриці похідних, що посилюється великим порядком матриці. Все це призводить до чисельної нестійкості при розв ’язанні зворотної задачі (задачі реконструкції). Метод зон провідності дозволяє оминути вказані труднощі: призводить (за рахунок додаткової структуризації фантому) до радикального зменшення порядку матриці похідних, підвищення чутливості та обумовленості такої матриці. Для оцінки приведені значення передаточних опорів та похідних для рівномірного фантому та фантомів з неодно-рідностями.
Ключові слова: імпедансна томографія, фантом, кінцевий елемент, зони провідності, пряма задача, зворотна задача, задача реконструкції зображення, матриця похідних.
Sushko I.O., Rybin O.I. Comparison of classical and conductivity “zones” methods for solving EIT inverse problem. Comparison of classical iteration imagine reconstruction method and conductivity zones method properties in EITfrom the point of view of calculation organization, expecting precision and computational cost of imagine receiving is carried out. Reasonably small transfer resistances sensitivity to surface conductivity changes of separate finite element is shown. Small sensitivity is appeared for finite elements, which are situated in the central part of phantom, and also for finite elements, which are situated in the phantom outline. Small sensitivity is the cause of derivative matrix bad conditionality, that is aggravated with big matrix order. It led to quantity instability by solving of inverse problem (reconstruction problem). Conductivity zones method allows to avoid some difficulties: leads to derivative matrix order reduction, increasing of sensitivity and matrix conditionality. Transfer resistances and derivatives values for uniform phantom and for phantoms with different inhomogeneities for assessment are proposed.
Keywords: Electrical Impedance Tomography, phantom, finite element, conductivity zones, forward problem, inverse problem, imagine reconstruction problem, derivative matrix.
Вісник Національного технічного університету України "КПІ" Серія — Радіотехніка. Радіоапаратобудування.-2012.-№49
177
