способ неинвазивного электрофизиологического исследования сердца
Классы МПК: | A61B5/0402 электрокардиография, те ЭКГ A61B5/055 с применением электронного или ядерного магнитного резонанса (ЯМР), например получение изображения с помощью магнитного резонанса A61B6/03 томографы с применением вычислительной техники |
Автор(ы): | Ревишвили Амиран Шотаевич (RU), Калинин Виталий Викторович (RU), Калинин Александр Викторович (RU) |
Патентообладатель(и): | Ревишвили Амиран Шотаевич (RU) |
Приоритеты: |
подача заявки:
2008-11-27 публикация патента:
27.04.2011 |
Изобретение относится к кардиологии, сердечно-сосудистой хирургии, функциональной диагностике и клинической электрофизиологии сердца. Способ неинвазивного электрофизиологического исследования сердца включает следующие стадии: закрепление регистрирующих электродов на поверхности грудной клетки; регистрация ЭКГ; обработка ЭКГ-сигналов в режиме реального времени; ретроспективная обработка полученных ЭКГ; компьютерная или магнитно-резонансная томография грудной клетки; построение и редактирование компьютерных воксельных моделей органов грудной клетки и сердца; построение полигональных моделей торса и сердца; автоматическое определение координат регистрирующих электродов на поверхности грудной клетки; интерполяция значений ЭКГ-сигналов в узлы полигональной сетки; реконструкция потенциала электрического поля в заданных точках; визуализация результатов реконструкции электрического поля сердца; клиническая оценка результатов. При этом для построения воксельной модели используют алгоритм факторизации сдвига - деформации для преобразования просмотра. Стадия построения полигональных моделей включает: фильтрацию исходных воксельных моделей, построение триангуляционной поверхности, разреживание и улучшение качества сетки с использованием метода пуассоновской реконструкции. Изобретение позволяет повысить точность неинвазивной диагностики нарушений сердечного ритма и других сердечно-сосудистых заболеваний.7 з.п. ф-лы, 19 ил.
Формула изобретения
1. Способ неинвазивного электрофизиологического исследования сердца, включающий следующие стадии:
закрепление одноразовых регистрирующих электродов на поверхности грудной клетки;
регистрация ЭКГ во множестве однополюсных отведений с поверхности грудной клетки;
обработка ЭКГ-сигналов в режиме реального времени;
ретроспективная обработка полученных ЭКГ;
компьютерная томография (КТ) или магнитно-резонансная томография (МРТ) грудной клетки пациента закрепленными электродами;
построение по томографическим данным и редактирование компьютерных воксельных моделей органов грудной клетки и сердца, при этом для построения воксельной модели используют алгоритм факторизации сдвига - деформации для преобразования просмотра (Shear-Warp Factorization of the Viewing Transformation);
построение при помощи компьютерной программы полигональных моделей торса и сердца, причем стадия построения полигональных моделей включает следующие этапы:
фильтрация исходных вексельных моделей для уменьшения уровня случайного шума;
построение триангуляционной поверхности методом «марширующих кубов» или «методом исчерпывания» («advancing front method»);
разреживание и улучшение качества сетки с использованием метода пуассоновской реконструкции (Poisson Surface Reconstruction);
определение координат регистрирующих электродов на поверхности грудной клетки проводят в автоматическом режиме по данным КТ и МРТ;
интерполяция значений ЭКГ-сигналов в узлы полигональной сетки, которую осуществляют с использованием радиальных базисных функций;
реконструкция потенциала электрического поля в заданных точках грудной клетки, эпикардиальной поверхности сердца, поверхности межжелудочковой и межпредсердной перегородок;
визуализация результатов реконструкции электрического поля сердца в виде эпикардиальных электрограмм, изохронных и изопотенциальных карт, а также динамических карт (propagation maps) на полигональных моделях сердца и его структур;
клиническая оценка результатов.
2. Способ по п.1, в котором для КТ используют наклеиваемые металлические хлор-серебряные электроды, а для МРТ - наклеиваемые графитовые электроды.
3. Способ по п.1, в котором одноразовые электроды закрепляют в виде горизонтальных пяти - восьми поясов, расположенных на одинаковых расстояниях по вертикали, причем первый пояс располагают на уровне грудинно-ключичного сочленения, а последний пояс - на уровне нижнего края реберной поверхности и каждый пояс включает от 16 до 30 электродов, расположенных на одинаковых расстояниях по окружности грудной клетки.
4. Способ по п.1, в котором реконструкцию потенциала электрического поля сердца проводят путем численного решения задачи Коши для уравнения Лапласа методом граничных элементов, включающим решение возникающей в результате применения метода граничных элементов итоговой системы матрично-векторных уравнений
при помощи итерационной процедуры
при этом на каждом шаге итерационного процесса для решения уравнения (13) используют регуляризирующий метод решения, выбранный из группы: метод регуляризации Тихонова, в котором параметр регуляризации определяют по формуле
где - параметр регуляризации, 0 - малый действительный параметр, зависящий от погрешности задания граничных условий обратной задачи электрокардиографии, р - положительный действительный параметр, зависящий от скорости сходимости итерационной процедуры (11)-(13), - положительный действительный параметр, зависящий от точности начального приближения (11) в итерационной процедуре (11)-(13),
k - номер итерации в итерационной процедуре (11)-(13),
или
регуляризирующий алгоритм на основе SVD-разложения матрицы уравнения (13) с заменой нулями сингулярных чисел, меньших заданного положительного числа , причем параметр определяют согласно формуле:
= 0+ ·p-(k/2), где 0 - малый действительный параметр, зависящий от погрешности задания граничных условий обратной задачи электрокардиографии, р - положительный действительный параметр, зависящий от скорости сходимости итерационной процедуры (11)-(13), - положительный действительный параметр, зависящий от точности начального приближения (11) в итерационной процедуре (11)-(13), k - номер итерации в итерационной процедуре (11)-(13); или
регуляризирующий алгоритм решения уравнения (13) на основе итерационного метода обобщенных минимальных невязок (Generalized minimal residual method) с ограничением числа итераций, причем требуемое число итераций, требуемых для решения уравнения (13), определяют по формуле: n=n0+ ·k,
где n - число итераций алгоритма обобщенных минимальных невязок, k - номер итерации в итерационной процедуре (11)-(13), n0 и - положительные целые числа, зависящие от точности начального приближения (11) и скорости сходимости итерационной процедуры (11)-(13),
общее число итераций алгоритма (11)-(13) определяют по принципу невязки (принцип Морозова).
5. Способ по п.4, в котором для решения уравнения
системы матрично-векторных уравнений используется метод регуляризации Тихонова, причем параметр регуляризации определяют по формуле
где - параметр регуляризации, 0 - малый действительный параметр, зависящий от погрешности задания граничных условий обратной задачи электрокардиографии, р - положительный действительный параметр, зависящий от скорости сходимости итерационной процедуры, - положительный действительный параметр, зависящий от точности начального приближения в итерационной процедуре, k - номер итерации.
6. Способ по п.4, в котором для решения уравнения
системы матрично-векторных уравнений используется регуляризирующий алгоритм на основе SVD-разложения матрицы уравнения с заменой нулями сингулярных чисел, меньших заданного положительного числа , причем параметр определяют согласно формуле
= 0+ ·p-(k/2),
где 0 - малый действительный параметр, зависящий от погрешности задания граничных условий обратной задачи электрокардиографии;
р - положительный действительный параметр;
зависящий от скорости сходимости итерационной процедуры;
- положительный действительный параметр, зависящий от точности начального приближения в итерационной процедуре;
k - номер итерации.
7. Способ по п.4, в котором для решения уравнения
системы матрично-векторных уравнений используется регуляризирующий алгоритм на основе итерационного метода обобщенных минимальных невязок (Generalized minimal residual method) с ограничением числа итераций, причем требуемое число итераций определяют по формуле
n=n0+ ·k,
где n - число итераций алгоритма;
k - номер итерации в общей итерационной процедуре;
n0 и - положительные целые числа, зависящие от точности начального приближения и скорости сходимости процедуры (11)-(13)
8. Способ по п.4, в котором уравнения системы матрично-векторных уравнений решаются на основе «быстрого мультипольного метода» (Fast Multipole Method).
Описание изобретения к патенту
Класс A61B5/0402 электрокардиография, те ЭКГ
Класс A61B5/055 с применением электронного или ядерного магнитного резонанса (ЯМР), например получение изображения с помощью магнитного резонанса
Класс A61B6/03 томографы с применением вычислительной техники