Try ICA (Independent component analysis) to localize the source in 3d

[annitrolla]

I tried fastICA on our sample. If you look at a resulting picture under Images, you will see why I am confused - initial sources look more unmixed than the final ones. Also, I tried ICA both on python and R, the exact same package for R requires to set number of components (in our case those dipoles hidden inside), however, fastICA in python does not require them and comes up automatically with the amount equal to the number of electrodes. In further work I will set it up to some number (any ideas?). According to the paper, setting it to equal amount is not right: "the number of independent signal sources does not exceed the number of electrodes". There are also few ideas, how to proceed further to estimate the concrete point of a signal, but let us discuss this later.

[kuz]

На первой картинке по оси Y идут микровольты -- сила сигнала. Из-за особенности устройства каждый электрод имеет некий offset, поэтому они "на разной высоте" Первым шагом тут должна быть нормализация (каждого электрода в отдельности) чтобы они были вокруг нуля. То есть то, что они "раздельно" это никак не связано с тем, что они unmixed каким-либо осмысленным образом, у них просто bias у каждого свой.

Вторую картинку я не понимаю. По Y что идет? Ужатые до промежутка -0.15 до 0.15 амплитуды? А почему? А как это связано с положением источника сигнала в 3D пространстве?

Для начала можно попробовать попросить алгоритм искать 1 источник.

[annitrolla]

Первым шагом тут должна быть нормализация

• ICA уже нормализует данные. Для начала можно попробовать попросить алгоритм искать 1 источник.
• Мне кажется, один источник в данном случае не имеет смысла. У нас есть предположение, что наши наблюдения это смешанный результат из нескольких источников, менее двух - это уже не разделение таковых, а избавление сигнала от шума. По Y что идет? -Если я правильно понимаю, по Y нормализованные результаты. А как это связано с положением источника сигнала в 3D пространстве? -Про 3д - я говорила когда-то с костей, ica недостаточно. мы разделяем сигналы с помощью иса, теперь надо придумать как найти, откуда они. он предложил попробовать что-то по типу attenuation - сигналы затихают одинаково. Зная силу каждого источника (веса иса), можно посмотреть, где каждый сигнал затухает, их пересечение - наш результат. Но в статьей что-то сложное описывается, я не успела разобраться до конца. к тому же в нашем варианте половина формул не отражается, альтернативный я нашла, но он другой, не идентичный. С интернетом на конференции полная фигня.

[annitrolla]

Аааа..сорри..слетели весь форматтинг текста :(

[kuz]

Окей, пусть будет 2 источника. В статье что-то про это было, там была какая-то матрица C и какие-то расстояния d, я думаю это та же идея о которой я в начале говорил -- считаем "время" от каждого электрода и та же что и Костин attenuation. Нужно копать в общем.