Объявление

Свернуть
Пока нет объявлений.

Картирование ЭЭГ

Свернуть
X
  • Фильтр
  • Время
  • Показать
Очистить всё
новые сообщения

  • Картирование ЭЭГ

    Всем привет, бегло посмотрев форум не нашел толком такой темы. Вообщем поставили в институте задачу, по написанию программы картирования сигнала ээг. Сами данные анализа сняты с энцефалографа Нейрон-Спектр 5 и сохранены в формате ASCII. Теперь пользуясь ими надо написать программу на delphi для картирования.
    Порывшись в литературе нашел алгоритм картирования в книжке: Сахаров В.Л., Андреенко А.С. "МЕТОДЫ МАТЕМАТИЧЕСКОЙ ОБРАБОТКИ
    ЭЛЕКТРОЭНЦЕФАЛОГРАММ" Вроде все подробно и просто распиисано, но на 21-ой странице приведена формула для подсчета значений в узлах сетки U = Σ( A[Ch]*μ[Ch ]) где Ch - номер текущего канала;
    N - число каналов(Σ от 1 до N); A - средняя амплитуда каждого отведения; μ - коэффициент, получаемый в результате расчетов с использованием уравнения Пуассона, характеризует затухание сигнала в диэлектрической среде, каковой является мозг человека, а также зависит от расстояния между точкой расположения датчика и точкой в узле, где рассчитывается значение потенциала. При этом считалось, что влияние каждого отведения не должно перекрываться более чем на 20%, а также не должно быть участков не охваченных полем какого-либо канала.

    Вот с коэффициентом μ возникли непонятки, как его считать??? Книжка Сахарова ссылается на книгу Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1977, 736 с. Но там я толком ничего не нашел.

    Вообщем суть вопроса состоит в том: может кто-то сталкивался с подобной проблемой и может помочь с нахождением этого коэффициента, а возможно кто-то знает другой алгоритм картирования?

    Книжку Сахарова в формате pdf можно слить здесь, возможно понадобится кому-нибудь в будущем: http://filesurf.ru/76321 (файл eeg_saharov.pdf размером 584,756 байт)

  • #2
    Re: Картирование ЭЭГ

    1. В нейрософтовских прогах картирование уже встроено, но если хочется лично написать что-то своё то тут надо 2 варианта а). для дела б). для "понта"(не сочтите за слэнг)
    а). Нужды в визуализации особой нет т.к. полученные значения и так можно загонять в обработчики более высокого уровня(не считая эксперта :))))
    б). Означенный коэффициент вряд ли эффективен потому как
    1). плотность расположения датчиков не одинакова и следовательно в самом картировании возникают области большего или меньшего наложения оссюда - коррекция формы и площади "закрашиваемой области": явная сложность
    2). Для каждого электрода свой коэффициент хотябы из-за анатомических особенностей опять возникают вопосы расчёта
    отюда рекомендую следующие подходы
    а). Поставить коэффицинт от фонаря и поиграться с качеством визуализации (ценность подхода не высока но занимательно)
    б). допустить что между точками амплитуда меняется линейно (или по любому другому закону) и так прорисовать (получается немного резко но читаемо)
    в). представть значения амплитуд как матрицу и аппроксимировать сплайнами , вейвлетами и т.д. (внешне красиво, но содержательность страдает особенно если
    криво выбрать фунукцию).
    еще можно добавть нормиронание на максимальную амплитуду. для наглядности

    Впрочем всё же желательно что б Вы четко сформулировали цель тогда возможно что-то конкретное

    Комментарий


    • #3
      Re: Картирование ЭЭГ

      Кстати в сахаровеэти алгоритмы более-менее явно изложены

      Комментарий


      • #4
        Re: Картирование ЭЭГ

        лучше не заморачиваться и использовать обычную двухмерную интерполяцию на неравномерной сетке. Если хочешь более изощренно, то ищи алгоритмы сферической интреполяции. Привожу код на матлаб, некоторые умельцы переделывают под С++


        % вектор со значениями амплитуд на отведениях
        Am=[132.45 117.9 74.25 107.2 95.75 95.3 76.1 65.7 73.75 94.6 110.2 79.95 58.8 60.85 76.65 89.9 101.3 84.15 124.2]
        % трехмерные коррдинаты электродов ( из Гнездицкий В.В. Обратная задача ЭЭГ и клиническая электроэнцефалография)
        x=[94 94 54 54 66 54 54 0 0 0 0 0 -54 -54 -68 -54 -54 -85 -85];
        y=[30 -30 66 50 0 -50 -66 80 66 0 -66 -80 66 46 0 -46 -66 30 -30];
        z=[34 34 36 74 97 74 36 36 96 120 96 36 36 77 98 77 36 34 34];

        x1=[-25 25 -50 -25 0 25 50 -50 -25 0 25 50 -50 -25 0 25 50 -25 25];
        y1=[50 50 25 25 25 25 25 0 0 0 0 0 -25 -25 -25 -25 -25 -50 -50]; % относительные координаты на плоскости для целей интерполяции

        % ОБЫЧНАЯ ИНТЕРПОЛЯЦИЯ (алгоритм матлаб v4)
        N=.5;
        [X1,Y1] = meshgrid(-50:N:50, -50:N:50);
        AmC = griddata(x1,y1, Am, X1,Y1, 'v4');
        Rz=sqrt(X1.^2+Y1.^2); AmC(Rz>50)=NaN; % убираем лишнюю информацию
        surf (AmC)
        colorbar
        SHADING interp, grid minor, axis off
        view(2)

        %figure
        % ПСЕВДОСФЕРИЧЕСКАЯ ИНТЕРПОЛЯЦИЯ
        N=3;
        [X1,Y1] = meshgrid(-50:N:50, -50:N:50);
        XC = griddata(x1,y1, x, X1,Y1, 'v4');
        YC = griddata(x1,y1, y, X1,Y1, 'v4');
        AmC = griddata(x1,y1, Am, X1,Y1, 'cubic');
        XC=XC+AmC-AmC;
        YC=YC+AmC-AmC;
        surf (XC, YC, (AmC))
        SHADING FACETED;
        grid minor, axis off,
        view(-90, 90)
        colorbar

        вот результат:

        Комментарий


        • #5
          Re: Картирование ЭЭГ

          По-мойму в Сахарове этом инчего толком не написано, буду искать еще.
          asd, спасибо за пример в матлабе, гляну на выходных.
          Последний раз редактировалось Neuro-78; 14.11.2008, 09:02 PM.

          Комментарий


          • #6
            Re: Картирование ЭЭГ

            Толком ничего но идей там пруд пруди. Однако всё же разговор был бы предметен если бы Вы четко сформулировали цель визуализаци.. У нас студент на курсовой по такй же тематике OPEN GL прикрутил и шарик переливался у него всеми цветами. медики аж рыдали от восторга, хотя содержательная часть была никакой.

            Комментарий


            • #7
              Re: Картирование ЭЭГ

              Возник также вопрос, кто как ищет основные ритмы? А то я подозреваю что у меня не правильно выходит, получаются плавные изменения без скачков киких-либо.

              Комментарий


              • #8
                Re: Картирование ЭЭГ

                Не совсем понятно, что вам требуется.
                Если вас интресуют дельта, тета альфа и бета частотные диапазоны (так говорить более правильно, потому, как ритмы - волны (или гармоники, что опять же более точно, после фурье преобразования) , которые находятся примерно в этих диапазонах, что нестрого, так как например альфа ритм может быть и 7 Гц и 15 Гц (а в среднем 10 Гц)- индивидуальные особенности никто не отменял) - то ищут их классически - фурье преобразование. К дельта обычно относят диапазон .5-3 Гц, тета - 4-6 Гц, альфа 8-13 Гц, бета 14-40 Гц (в разных руководствах могут отличаться границы частотных диапазонов - и деление это честно говоря весьма условное). Амплитуду диапазонов указывать вообще по моему личному мнению неправильно, так как вариации значительны, и часто зависят от аппаратуры. То что нет скачков - а они должны быть? в природе это встречается очень редко.

                Комментарий


                • #9
                  Re: Картирование ЭЭГ

                  небольшое дополнение к коду (украшательство)
                  %figure
                  % ПСЕВДОСФЕРИЧЕСКАЯ ИНТЕРПОЛЯЦИЯ
                  N=3;
                  [X1,Y1] = meshgrid(-50:N:50, -50:N:50);
                  XC = griddata(x1,y1, x, X1,Y1, 'v4');
                  YC = griddata(x1,y1, y, X1,Y1, 'v4');
                  ZC = griddata(x1,y1, z, X1,Y1, 'cubic');
                  AmC = griddata(x1,y1, Am, X1,Y1, 'cubic');

                  XC=XC+AmC-AmC;
                  YC=YC+AmC-AmC;

                  I=mat2gray(AmC);
                  warp(XC, YC, ZC, I)
                  colormap jet
                  axis off
                  view([130 70]);
                  Последний раз редактировалось asd; 01.09.2009, 06:34 PM.

                  Комментарий


                  • #10
                    Re: Картирование ЭЭГ

                    Можно я тоже вставлю 5 копеек :) ?

                    1. Как заметили ув. форумчане, единственным правильным способом в данном случаи является эллиптическая (в упрощении сферическая) сплайн интерполяция на неравномерной сетке. Задача это не такая уж и простая, поэтому большинство используют упрощения (на подобии мю коэффициентов у Сахарова).

                    2. Поэтому, чтобы получить хороший результат крайне рекомендую разобраться в матлабовской ф-ии gdatav4(), которая выполняет сплайн интерполяцию в плоскости на неравномерной сетке (для сплайна как указывалось используется V4). Точное название алгоритма - "Biharmonic spline interpolation of GEOS-3 and SEASAT altimeter data".
                    Но этот алгоритм нуждается в доработке, в серьезных работах просто так его использовать нельзя, т.к. он дает сильные искажения на граничных участках. Как и любой сплайн он не знает о поведении сигнала за границей и использует для этого априорную ин-ю о производной сигнала в этой точке, что для задач ЭЭГ неверно. Поэтому без сферических интерполяций не обойтись.

                    Комментарий


                    • #11
                      Re: Картирование ЭЭГ

                      Вопрос опять: "ЗАЧЕМ?"

                      Комментарий

                      Обработка...
                      X
                      yobit wp super cache yobit.net официальный сайт yobit.net freewallet bgogo