Description_Leadfield_Matrix.ipynb

{
 "cells": [
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this project, we will calculate a lead-field matrix that can convert brain activity in the AAL2 atlas to a 64 electrode EEG cap.\n",
    "A lead-field matrix is a matrix `M` s.th.\n",
    "`y_m = M_mn x_n`\n",
    "where `y_m` is the EEG activity at electrode `m` and `x_n` is the neural activity at voxel `n`. Note: I will call the elements of a brain that can generate activity a \"voxel\", but they could also be represented as \"nodes\" in a mesh representation in Field trip.\n",
    "These are the steps:\n",
    "- First, we calculate a high-resolution lead-field matrix using Field Trip. Field trip converts the neural activity that is represented as dipoles on each voxel of the brain into the electric field at the electrodes. These dipoles generate an electric field that is simulated until it reaches the EEG electrodes.\n",
    "- Typically, this problem is done in reverse: you start with EEG data `y` and want to end up with brain activity `x`. We have the inverse problem: we want to know what `y` is, if *ALL* voxels `x` are active.\n",
    "- Electric fields are linear that means we can compute the *overall* lead field matrix that measure the effect of all voxels from the *individual* lead field matrices that measure only the effect of a single voxel.\n",
    "- I can imagine there are two ways of doing this: either there is a way in field trip to \"activate\" the entire brain at once and get the *overall* lead field matrix (probably uses some kind of data structure to hold this), or we would have to go through the brain one voxel `n` (example: `n=1`) at a time, activate that voxel and measure the lead field matrix `M_m1` (it's actually just a column of the overall matrix) for that voxel. We do that `N` times such that we get `N` lead field matrices that we can just sum up to get the overall one: `M_mn = M_m1 + M_m2 + ... + M_mN`\n",
    "- The first way of doing this would be a lot more efficient.\n",
    "- We then get a high-res lead field matrix (the *overall* matrix) with `N` columns.\n",
    "- Export this matrix into a python-usable format (simple numpy array). We can work with this matrix later in Python.\n",
    "- Now we want to create a low-res lead field matrix defined over the AAL2 atlas.\n",
    "- For this, we need to summarize the columns of the high-res matrix such that values that belong to the same AAL2 brain region are reduced into one column.\n",
    "- We can take a column of this matrix that corresponds to voxels that belong to the same brain area and simply take the average EEG activity that this brain area causes by simply taking the average of these columns and using it as one.\n",
    "- Now we should have a low-dimensional lead field matrix with `n` being 120 (or: 80 only cortical areas, excluding subcortical).\n",
    "# Testing results\n",
    "- Now that we have the lead field matrix, we can do cool things with neurolib.\n",
    "- We can simulate brain activity `x` and calculate EEG activity `y` by simply multiplying the activity with the matrix.\n",
    "- We now have EEG activity that we can compare with the original activity `x` and see which effects this introduces. First thing to check would be to simulate a brain that produces slow oscillations and compare the power spectra of `y` and `x`.\n",
    "- Same for a brain that is oscillating at higher frequencies.\n",
    "- The neurolib whole-brain model parameters to simulate slow oscillations or high-frequency oscillations will be shared with you.\n",
    "The end.\n",
    "\n",
    "Matlab Script: https://www.fieldtriptoolbox.org/workshop/ohbm2018/forward/ \n",
    "\n",
    "Forward Model: https://www.fieldtriptoolbox.org/development/module/forward/\n",
    "\n",
    "MNE: https://mne.tools/stable/auto_tutorials/forward/index.html\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'__header__': b'MATLAB 5.0 MAT-file, Platform: GLNXA64, Created on: Mon Jul  3 18:59:45 2023', '__version__': '1.0', '__globals__': [], 'leadfield_fem_eeg': array([[(array([[34, 27, 26]], dtype=uint8), array([[-60, -65, -10],\n",
      "               [-55, -65, -10],\n",
      "               [-50, -65, -10],\n",
      "               ...,\n",
      "               [ 95,  65, 115],\n",
      "               [100,  65, 115],\n",
      "               [105,  65, 115]], dtype=int16), array(['mm'], dtype='<U2'), array([[0],\n",
      "               [0],\n",
      "               [0],\n",
      "               ...,\n",
      "               [0],\n",
      "               [0],\n",
      "               [0]], dtype=uint8), array([[(array(['empty - this was cleared by checkconfig'], dtype='<U39'), array(['empty - this was cleared by checkconfig'], dtype='<U39'), array([[(array([[ 28.15, -71.85,  -1.85],\n",
      "                               [ 16.85,  81.85,   4.15],\n",
      "                               [-73.  ,   0.  ,   5.  ],\n",
      "                               [-72.85, -25.85,  38.  ],\n",
      "                               [-78.  ,   0.  ,  43.  ],\n",
      "                               [-75.85,  27.85,  38.  ],\n",
      "                               [-54.85, -44.85,  -0.85],\n",
      "                               [-57.85, -48.85,  37.15],\n",
      "                               [-63.85, -42.85,  50.85],\n",
      "                               [-65.85, -31.85,  63.85],\n",
      "                               [-66.85, -16.85,  74.85],\n",
      "                               [-67.85,   0.  ,  75.85],\n",
      "                               [-68.85,  20.  ,  73.85],\n",
      "                               [-68.85,  33.85,  65.85],\n",
      "                               [-67.  ,  46.  ,  51.  ],\n",
      "                               [-60.85,  51.85,  36.15],\n",
      "                               [-55.85,  45.85,  -2.85],\n",
      "                               [-24.  , -70.  ,  -5.  ],\n",
      "                               [-35.85, -63.85,  36.15],\n",
      "                               [-41.85, -61.85,  62.  ],\n",
      "                               [-45.85, -48.85,  85.85],\n",
      "                               [-47.85, -26.85, 101.85],\n",
      "                               [-48.85,   0.  , 105.85],\n",
      "                               [-48.85,  27.85, 101.85],\n",
      "                               [-48.85,  51.85,  86.  ],\n",
      "                               [-43.85,  67.85,  61.  ],\n",
      "                               [-37.85,  67.85,  36.15],\n",
      "                               [-23.85,  71.85,  -1.85],\n",
      "                               [-11.15, -74.85,  -7.15],\n",
      "                               [ -8.85, -78.85,  28.85],\n",
      "                               [-12.85, -74.85,  68.85],\n",
      "                               [-18.  , -60.  , 101.  ],\n",
      "                               [-19.85, -32.85, 122.85],\n",
      "                               [-18.85,   0.  , 126.85],\n",
      "                               [-19.85,  33.85, 121.85],\n",
      "                               [-19.  ,  62.  , 101.  ],\n",
      "                               [-13.85,  75.85,  69.  ],\n",
      "                               [-10.  ,  82.  ,  34.  ],\n",
      "                               [-10.  ,  76.  ,  -6.  ],\n",
      "                               [ 11.85, -75.85,  -0.85],\n",
      "                               [ 21.  , -76.85,  37.15],\n",
      "                               [ 19.  , -77.85,  73.85],\n",
      "                               [ 15.85, -61.85, 105.85],\n",
      "                               [ 12.85, -32.85, 127.85],\n",
      "                               [ 14.  ,   0.  , 135.  ],\n",
      "                               [ 15.  ,  33.85, 127.85],\n",
      "                               [ 15.85,  63.85, 105.85],\n",
      "                               [ 17.85,  77.85,  73.85],\n",
      "                               [ 20.  ,  77.85,  37.15],\n",
      "                               [ 15.85,  81.85,   3.15],\n",
      "                               [ 47.85, -74.85,   1.15],\n",
      "                               [ 50.  , -76.  ,  38.  ],\n",
      "                               [ 48.85, -72.85,  72.85],\n",
      "                               [ 48.85, -54.85, 101.85],\n",
      "                               [ 44.85, -27.85, 121.85],\n",
      "                               [ 45.85,   0.15, 127.85],\n",
      "                               [ 46.85,  30.85, 120.85],\n",
      "                               [ 47.85,  56.85, 100.85],\n",
      "                               [ 46.85,  72.85,  71.85],\n",
      "                               [ 46.85,  76.85,  38.15],\n",
      "                               [ 49.85,  73.85,   3.15],\n",
      "                               [ 73.85, -64.85,   8.15],\n",
      "                               [ 74.85, -64.85,  43.  ],\n",
      "                               [ 74.85, -57.85,  70.85],\n",
      "                               [ 76.85, -45.85,  93.85],\n",
      "                               [ 78.85, -22.85, 108.85],\n",
      "                               [ 79.85,   0.  , 112.85],\n",
      "                               [ 77.85,  24.85, 108.85],\n",
      "                               [ 74.85,  46.85,  93.85],\n",
      "                               [ 74.85,  57.85,  68.85],\n",
      "                               [ 73.85,  64.85,  42.  ],\n",
      "                               [ 71.85,  64.85,   5.15],\n",
      "                               [ 99.85, -50.85,   9.15],\n",
      "                               [ 97.85, -48.85,  44.85],\n",
      "                               [ 98.85, -42.85,  63.85],\n",
      "                               [ 99.85, -32.85,  79.85],\n",
      "                               [103.85, -16.85,  88.85],\n",
      "                               [103.85,  -0.15,  92.85],\n",
      "                               [102.85,  18.85,  88.85],\n",
      "                               [ 99.85,  32.85,  77.85],\n",
      "                               [ 96.85,  42.85,  65.85],\n",
      "                               [ 96.85,  47.85,  45.85],\n",
      "                               [ 99.85,  49.85,  10.15],\n",
      "                               [114.85, -27.85,  51.85],\n",
      "                               [120.  ,   0.  ,  60.  ],\n",
      "                               [113.85,  26.85,  52.85],\n",
      "                               [117.85, -28.85,  16.  ],\n",
      "                               [122.  ,   0.  ,  22.  ],\n",
      "                               [116.85,  27.85,  16.  ],\n",
      "                               [ 21.  , -76.85,  37.15],\n",
      "                               [ 74.85, -64.85,  43.  ],\n",
      "                               [ 20.  ,  77.85,  37.15],\n",
      "                               [ 73.85,  64.85,  42.  ],\n",
      "                               [ 47.85, -72.85, -23.  ],\n",
      "                               [ 49.  ,  73.  , -23.  ],\n",
      "                               [ 36.15, -72.85, -20.85],\n",
      "                               [ 36.15,  71.85, -23.15]]), array([[array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')],\n",
      "                               [array(['eeg'], dtype='<U3')]], dtype=object), array([[array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')],\n",
      "                               [array(['V'], dtype='<U1')]], dtype=object), array([[ 28.15, -71.85,  -1.85],\n",
      "                               [ 16.85,  81.85,   4.15],\n",
      "                               [-73.  ,   0.  ,   5.  ],\n",
      "                               [-72.85, -25.85,  38.  ],\n",
      "                               [-78.  ,   0.  ,  43.  ],\n",
      "                               [-75.85,  27.85,  38.  ],\n",
      "                               [-54.85, -44.85,  -0.85],\n",
      "                               [-57.85, -48.85,  37.15],\n",
      "                               [-63.85, -42.85,  50.85],\n",
      "                               [-65.85, -31.85,  63.85],\n",
      "                               [-66.85, -16.85,  74.85],\n",
      "                               [-67.85,   0.  ,  75.85],\n",
      "                               [-68.85,  20.  ,  73.85],\n",
      "                               [-68.85,  33.85,  65.85],\n",
      "                               [-67.  ,  46.  ,  51.  ],\n",
      "                               [-60.85,  51.85,  36.15],\n",
      "                               [-55.85,  45.85,  -2.85],\n",
      "                               [-24.  , -70.  ,  -5.  ],\n",
      "                               [-35.85, -63.85,  36.15],\n",
      "                               [-41.85, -61.85,  62.  ],\n",
      "                               [-45.85, -48.85,  85.85],\n",
      "                               [-47.85, -26.85, 101.85],\n",
      "                               [-48.85,   0.  , 105.85],\n",
      "                               [-48.85,  27.85, 101.85],\n",
      "                               [-48.85,  51.85,  86.  ],\n",
      "                               [-43.85,  67.85,  61.  ],\n",
      "                               [-37.85,  67.85,  36.15],\n",
      "                               [-23.85,  71.85,  -1.85],\n",
      "                               [-11.15, -74.85,  -7.15],\n",
      "                               [ -8.85, -78.85,  28.85],\n",
      "                               [-12.85, -74.85,  68.85],\n",
      "                               [-18.  , -60.  , 101.  ],\n",
      "                               [-19.85, -32.85, 122.85],\n",
      "                               [-18.85,   0.  , 126.85],\n",
      "                               [-19.85,  33.85, 121.85],\n",
      "                               [-19.  ,  62.  , 101.  ],\n",
      "                               [-13.85,  75.85,  69.  ],\n",
      "                               [-10.  ,  82.  ,  34.  ],\n",
      "                               [-10.  ,  76.  ,  -6.  ],\n",
      "                               [ 11.85, -75.85,  -0.85],\n",
      "                               [ 21.  , -76.85,  37.15],\n",
      "                               [ 19.  , -77.85,  73.85],\n",
      "                               [ 15.85, -61.85, 105.85],\n",
      "                               [ 12.85, -32.85, 127.85],\n",
      "                               [ 14.  ,   0.  , 135.  ],\n",
      "                               [ 15.  ,  33.85, 127.85],\n",
      "                               [ 15.85,  63.85, 105.85],\n",
      "                               [ 17.85,  77.85,  73.85],\n",
      "                               [ 20.  ,  77.85,  37.15],\n",
      "                               [ 15.85,  81.85,   3.15],\n",
      "                               [ 47.85, -74.85,   1.15],\n",
      "                               [ 50.  , -76.  ,  38.  ],\n",
      "                               [ 48.85, -72.85,  72.85],\n",
      "                               [ 48.85, -54.85, 101.85],\n",
      "                               [ 44.85, -27.85, 121.85],\n",
      "                               [ 45.85,   0.15, 127.85],\n",
      "                               [ 46.85,  30.85, 120.85],\n",
      "                               [ 47.85,  56.85, 100.85],\n",
      "                               [ 46.85,  72.85,  71.85],\n",
      "                               [ 46.85,  76.85,  38.15],\n",
      "                               [ 49.85,  73.85,   3.15],\n",
      "                               [ 73.85, -64.85,   8.15],\n",
      "                               [ 74.85, -64.85,  43.  ],\n",
      "                               [ 74.85, -57.85,  70.85],\n",
      "                               [ 76.85, -45.85,  93.85],\n",
      "                               [ 78.85, -22.85, 108.85],\n",
      "                               [ 79.85,   0.  , 112.85],\n",
      "                               [ 77.85,  24.85, 108.85],\n",
      "                               [ 74.85,  46.85,  93.85],\n",
      "                               [ 74.85,  57.85,  68.85],\n",
      "                               [ 73.85,  64.85,  42.  ],\n",
      "                               [ 71.85,  64.85,   5.15],\n",
      "                               [ 99.85, -50.85,   9.15],\n",
      "                               [ 97.85, -48.85,  44.85],\n",
      "                               [ 98.85, -42.85,  63.85],\n",
      "                               [ 99.85, -32.85,  79.85],\n",
      "                               [103.85, -16.85,  88.85],\n",
      "                               [103.85,  -0.15,  92.85],\n",
      "                               [102.85,  18.85,  88.85],\n",
      "                               [ 99.85,  32.85,  77.85],\n",
      "                               [ 96.85,  42.85,  65.85],\n",
      "                               [ 96.85,  47.85,  45.85],\n",
      "                               [ 99.85,  49.85,  10.15],\n",
      "                               [114.85, -27.85,  51.85],\n",
      "                               [120.  ,   0.  ,  60.  ],\n",
      "                               [113.85,  26.85,  52.85],\n",
      "                               [117.85, -28.85,  16.  ],\n",
      "                               [122.  ,   0.  ,  22.  ],\n",
      "                               [116.85,  27.85,  16.  ],\n",
      "                               [ 21.  , -76.85,  37.15],\n",
      "                               [ 74.85, -64.85,  43.  ],\n",
      "                               [ 20.  ,  77.85,  37.15],\n",
      "                               [ 73.85,  64.85,  42.  ],\n",
      "                               [ 47.85, -72.85, -23.  ],\n",
      "                               [ 49.  ,  73.  , -23.  ],\n",
      "                               [ 36.15, -72.85, -20.85],\n",
      "                               [ 36.15,  71.85, -23.15]]), array([[array(['LPA'], dtype='<U3')],\n",
      "                               [array(['RPA'], dtype='<U3')],\n",
      "                               [array(['Nz'], dtype='<U2')],\n",
      "                               [array(['Fp1'], dtype='<U3')],\n",
      "                               [array(['Fpz'], dtype='<U3')],\n",
      "                               [array(['Fp2'], dtype='<U3')],\n",
      "                               [array(['AF9'], dtype='<U3')],\n",
      "                               [array(['AF7'], dtype='<U3')],\n",
      "                               [array(['AF5'], dtype='<U3')],\n",
      "                               [array(['AF3'], dtype='<U3')],\n",
      "                               [array(['AF1'], dtype='<U3')],\n",
      "                               [array(['AFz'], dtype='<U3')],\n",
      "                               [array(['AF2'], dtype='<U3')],\n",
      "                               [array(['AF4'], dtype='<U3')],\n",
      "                               [array(['AF6'], dtype='<U3')],\n",
      "                               [array(['AF8'], dtype='<U3')],\n",
      "                               [array(['AF10'], dtype='<U4')],\n",
      "                               [array(['F9'], dtype='<U2')],\n",
      "                               [array(['F7'], dtype='<U2')],\n",
      "                               [array(['F5'], dtype='<U2')],\n",
      "                               [array(['F3'], dtype='<U2')],\n",
      "                               [array(['F1'], dtype='<U2')],\n",
      "                               [array(['Fz'], dtype='<U2')],\n",
      "                               [array(['F2'], dtype='<U2')],\n",
      "                               [array(['F4'], dtype='<U2')],\n",
      "                               [array(['F6'], dtype='<U2')],\n",
      "                               [array(['F8'], dtype='<U2')],\n",
      "                               [array(['F10'], dtype='<U3')],\n",
      "                               [array(['FT9'], dtype='<U3')],\n",
      "                               [array(['FT7'], dtype='<U3')],\n",
      "                               [array(['FC5'], dtype='<U3')],\n",
      "                               [array(['FC3'], dtype='<U3')],\n",
      "                               [array(['FC1'], dtype='<U3')],\n",
      "                               [array(['FCz'], dtype='<U3')],\n",
      "                               [array(['FC2'], dtype='<U3')],\n",
      "                               [array(['FC4'], dtype='<U3')],\n",
      "                               [array(['FC6'], dtype='<U3')],\n",
      "                               [array(['FT8'], dtype='<U3')],\n",
      "                               [array(['FT10'], dtype='<U4')],\n",
      "                               [array(['T9'], dtype='<U2')],\n",
      "                               [array(['T7'], dtype='<U2')],\n",
      "                               [array(['C5'], dtype='<U2')],\n",
      "                               [array(['C3'], dtype='<U2')],\n",
      "                               [array(['C1'], dtype='<U2')],\n",
      "                               [array(['Cz'], dtype='<U2')],\n",
      "                               [array(['C2'], dtype='<U2')],\n",
      "                               [array(['C4'], dtype='<U2')],\n",
      "                               [array(['C6'], dtype='<U2')],\n",
      "                               [array(['T8'], dtype='<U2')],\n",
      "                               [array(['T10'], dtype='<U3')],\n",
      "                               [array(['TP9'], dtype='<U3')],\n",
      "                               [array(['TP7'], dtype='<U3')],\n",
      "                               [array(['CP5'], dtype='<U3')],\n",
      "                               [array(['CP3'], dtype='<U3')],\n",
      "                               [array(['CP1'], dtype='<U3')],\n",
      "                               [array(['CPz'], dtype='<U3')],\n",
      "                               [array(['CP2'], dtype='<U3')],\n",
      "                               [array(['CP4'], dtype='<U3')],\n",
      "                               [array(['CP6'], dtype='<U3')],\n",
      "                               [array(['TP8'], dtype='<U3')],\n",
      "                               [array(['TP10'], dtype='<U4')],\n",
      "                               [array(['P9'], dtype='<U2')],\n",
      "                               [array(['P7'], dtype='<U2')],\n",
      "                               [array(['P5'], dtype='<U2')],\n",
      "                               [array(['P3'], dtype='<U2')],\n",
      "                               [array(['P1'], dtype='<U2')],\n",
      "                               [array(['Pz'], dtype='<U2')],\n",
      "                               [array(['P2'], dtype='<U2')],\n",
      "                               [array(['P4'], dtype='<U2')],\n",
      "                               [array(['P6'], dtype='<U2')],\n",
      "                               [array(['P8'], dtype='<U2')],\n",
      "                               [array(['P10'], dtype='<U3')],\n",
      "                               [array(['PO9'], dtype='<U3')],\n",
      "                               [array(['PO7'], dtype='<U3')],\n",
      "                               [array(['PO5'], dtype='<U3')],\n",
      "                               [array(['PO3'], dtype='<U3')],\n",
      "                               [array(['PO1'], dtype='<U3')],\n",
      "                               [array(['POz'], dtype='<U3')],\n",
      "                               [array(['PO2'], dtype='<U3')],\n",
      "                               [array(['PO4'], dtype='<U3')],\n",
      "                               [array(['PO6'], dtype='<U3')],\n",
      "                               [array(['PO8'], dtype='<U3')],\n",
      "                               [array(['PO10'], dtype='<U4')],\n",
      "                               [array(['O1'], dtype='<U2')],\n",
      "                               [array(['Oz'], dtype='<U2')],\n",
      "                               [array(['O2'], dtype='<U2')],\n",
      "                               [array(['O9'], dtype='<U2')],\n",
      "                               [array(['Iz'], dtype='<U2')],\n",
      "                               [array(['O10'], dtype='<U3')],\n",
      "                               [array(['T3'], dtype='<U2')],\n",
      "                               [array(['T5'], dtype='<U2')],\n",
      "                               [array(['T4'], dtype='<U2')],\n",
      "                               [array(['T6'], dtype='<U2')],\n",
      "                               [array(['M1'], dtype='<U2')],\n",
      "                               [array(['M2'], dtype='<U2')],\n",
      "                               [array(['A1'], dtype='<U2')],\n",
      "                               [array(['A2'], dtype='<U2')]], dtype=object), array(['eeg1010'], dtype='<U7'), array(['mm'], dtype='<U2'))]],\n",
      "                      dtype=[('chanpos', 'O'), ('chantype', 'O'), ('chanunit', 'O'), ('elecpos', 'O'), ('label', 'O'), ('type', 'O'), ('unit', 'O')]), array([[3]], dtype=uint8), array(['pedantic'], dtype='<U8'), array(['yes'], dtype='<U3'), array(['overwrite'], dtype='<U9'), array(['yes'], dtype='<U3'), array(['yes'], dtype='<U3'), array([[(array(['matlab'], dtype='<U6'), array(['matlab'], dtype='<U6'), array(['matlab'], dtype='<U6'), array([], shape=(0, 0), dtype=object))]],\n",
      "                      dtype=[('images', 'O'), ('stats', 'O'), ('signal', 'O'), ('cleanup', 'O')]), array([[(array([[(array(['empty - this was cleared by checkconfig'], dtype='<U39'), array(['empty - this was cleared by checkconfig'], dtype='<U39'), array(['empty - this was cleared by checkconfig'], dtype='<U39'), array([[3]], dtype=uint8), array(['pedantic'], dtype='<U8'), array(['yes'], dtype='<U3'), array(['overwrite'], dtype='<U9'), array(['yes'], dtype='<U3'), array(['yes'], dtype='<U3'), array([[(array(['matlab'], dtype='<U6'), array(['matlab'], dtype='<U6'), array(['matlab'], dtype='<U6'), array([], shape=(0, 0), dtype=object))]],\n",
      "                                      dtype=[('images', 'O'), ('stats', 'O'), ('signal', 'O'), ('cleanup', 'O')]))                                                                                                                                                                                                                                                                                                                                                                                                                                                              ]],\n",
      "                              dtype=[('sourcemodel', 'O'), ('headmodel', 'O'), ('elec', 'O'), ('reducerank', 'O'), ('checkpath', 'O'), ('showlogo', 'O'), ('outputfilepresent', 'O'), ('tracktimeinfo', 'O'), ('trackmeminfo', 'O'), ('toolbox', 'O')]), array(['20230503'], dtype='<U8'), array(['9.14.0.2286388 (R2023a) Update 3'], dtype='<U32'), array(['glnxa64'], dtype='<U7'), array(['ben'], dtype='<U3'), array(['ben'], dtype='<U3'), array(['/home/ben/NIP_Projects'], dtype='<U22'), array([[2023.      ,    7.      ,    3.      ,   18.      ,   28.      ,\n",
      "                                   7.077072]]), array([[982.302233]]), array([[7.86476237e+09]]))                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  ]],\n",
      "                      dtype=[('usercfg', 'O'), ('fieldtrip', 'O'), ('matlab', 'O'), ('computer', 'O'), ('hostname', 'O'), ('user', 'O'), ('pwd', 'O'), ('calltime', 'O'), ('proctime', 'O'), ('procmem', 'O')]), array([[(array(['/home/ben/NIP_Projects/fieldtrip-20230503/ft_prepare_leadfield.m'],\n",
      "                              dtype='<U64'), array(['$Id$'], dtype='<U4'))                          ]],\n",
      "                      dtype=[('name', 'O'), ('id', 'O')]), array(['no'], dtype='<U2'), array(['no'], dtype='<U2'), array(['text'], dtype='<U4'), array(['no'], dtype='<U2'), array(['no'], dtype='<U2'), array([[array(['LPA'], dtype='<U3')],\n",
      "                       [array(['RPA'], dtype='<U3')],\n",
      "                       [array(['Nz'], dtype='<U2')],\n",
      "                       [array(['Fp1'], dtype='<U3')],\n",
      "                       [array(['Fpz'], dtype='<U3')],\n",
      "                       [array(['Fp2'], dtype='<U3')],\n",
      "                       [array(['AF9'], dtype='<U3')],\n",
      "                       [array(['AF7'], dtype='<U3')],\n",
      "                       [array(['AF5'], dtype='<U3')],\n",
      "                       [array(['AF3'], dtype='<U3')],\n",
      "                       [array(['AF1'], dtype='<U3')],\n",
      "                       [array(['AFz'], dtype='<U3')],\n",
      "                       [array(['AF2'], dtype='<U3')],\n",
      "                       [array(['AF4'], dtype='<U3')],\n",
      "                       [array(['AF6'], dtype='<U3')],\n",
      "                       [array(['AF8'], dtype='<U3')],\n",
      "                       [array(['AF10'], dtype='<U4')],\n",
      "                       [array(['F9'], dtype='<U2')],\n",
      "                       [array(['F7'], dtype='<U2')],\n",
      "                       [array(['F5'], dtype='<U2')],\n",
      "                       [array(['F3'], dtype='<U2')],\n",
      "                       [array(['F1'], dtype='<U2')],\n",
      "                       [array(['Fz'], dtype='<U2')],\n",
      "                       [array(['F2'], dtype='<U2')],\n",
      "                       [array(['F4'], dtype='<U2')],\n",
      "                       [array(['F6'], dtype='<U2')],\n",
      "                       [array(['F8'], dtype='<U2')],\n",
      "                       [array(['F10'], dtype='<U3')],\n",
      "                       [array(['FT9'], dtype='<U3')],\n",
      "                       [array(['FT7'], dtype='<U3')],\n",
      "                       [array(['FC5'], dtype='<U3')],\n",
      "                       [array(['FC3'], dtype='<U3')],\n",
      "                       [array(['FC1'], dtype='<U3')],\n",
      "                       [array(['FCz'], dtype='<U3')],\n",
      "                       [array(['FC2'], dtype='<U3')],\n",
      "                       [array(['FC4'], dtype='<U3')],\n",
      "                       [array(['FC6'], dtype='<U3')],\n",
      "                       [array(['FT8'], dtype='<U3')],\n",
      "                       [array(['FT10'], dtype='<U4')],\n",
      "                       [array(['T9'], dtype='<U2')],\n",
      "                       [array(['T7'], dtype='<U2')],\n",
      "                       [array(['C5'], dtype='<U2')],\n",
      "                       [array(['C3'], dtype='<U2')],\n",
      "                       [array(['C1'], dtype='<U2')],\n",
      "                       [array(['Cz'], dtype='<U2')],\n",
      "                       [array(['C2'], dtype='<U2')],\n",
      "                       [array(['C4'], dtype='<U2')],\n",
      "                       [array(['C6'], dtype='<U2')],\n",
      "                       [array(['T8'], dtype='<U2')],\n",
      "                       [array(['T10'], dtype='<U3')],\n",
      "                       [array(['TP9'], dtype='<U3')],\n",
      "                       [array(['TP7'], dtype='<U3')],\n",
      "                       [array(['CP5'], dtype='<U3')],\n",
      "                       [array(['CP3'], dtype='<U3')],\n",
      "                       [array(['CP1'], dtype='<U3')],\n",
      "                       [array(['CPz'], dtype='<U3')],\n",
      "                       [array(['CP2'], dtype='<U3')],\n",
      "                       [array(['CP4'], dtype='<U3')],\n",
      "                       [array(['CP6'], dtype='<U3')],\n",
      "                       [array(['TP8'], dtype='<U3')],\n",
      "                       [array(['TP10'], dtype='<U4')],\n",
      "                       [array(['P9'], dtype='<U2')],\n",
      "                       [array(['P7'], dtype='<U2')],\n",
      "                       [array(['P5'], dtype='<U2')],\n",
      "                       [array(['P3'], dtype='<U2')],\n",
      "                       [array(['P1'], dtype='<U2')],\n",
      "                       [array(['Pz'], dtype='<U2')],\n",
      "                       [array(['P2'], dtype='<U2')],\n",
      "                       [array(['P4'], dtype='<U2')],\n",
      "                       [array(['P6'], dtype='<U2')],\n",
      "                       [array(['P8'], dtype='<U2')],\n",
      "                       [array(['P10'], dtype='<U3')],\n",
      "                       [array(['PO9'], dtype='<U3')],\n",
      "                       [array(['PO7'], dtype='<U3')],\n",
      "                       [array(['PO5'], dtype='<U3')],\n",
      "                       [array(['PO3'], dtype='<U3')],\n",
      "                       [array(['PO1'], dtype='<U3')],\n",
      "                       [array(['POz'], dtype='<U3')],\n",
      "                       [array(['PO2'], dtype='<U3')],\n",
      "                       [array(['PO4'], dtype='<U3')],\n",
      "                       [array(['PO6'], dtype='<U3')],\n",
      "                       [array(['PO8'], dtype='<U3')],\n",
      "                       [array(['PO10'], dtype='<U4')],\n",
      "                       [array(['O1'], dtype='<U2')],\n",
      "                       [array(['Oz'], dtype='<U2')],\n",
      "                       [array(['O2'], dtype='<U2')],\n",
      "                       [array(['O9'], dtype='<U2')],\n",
      "                       [array(['Iz'], dtype='<U2')],\n",
      "                       [array(['O10'], dtype='<U3')],\n",
      "                       [array(['T3'], dtype='<U2')],\n",
      "                       [array(['T5'], dtype='<U2')],\n",
      "                       [array(['T4'], dtype='<U2')],\n",
      "                       [array(['T6'], dtype='<U2')],\n",
      "                       [array(['M1'], dtype='<U2')],\n",
      "                       [array(['M2'], dtype='<U2')],\n",
      "                       [array(['A1'], dtype='<U2')],\n",
      "                       [array(['A2'], dtype='<U2')]], dtype=object), array(['no'], dtype='<U2'), array([], shape=(0, 0), dtype=uint8), array([], shape=(0, 0), dtype=uint8), array([], shape=(0, 0), dtype=uint8), array([], shape=(0, 0), dtype=uint8), array([], shape=(0, 0), dtype=uint8), array([], shape=(0, 0), dtype=object))                                                                                                                                                                                                                                                                                                                               ]],\n",
      "              dtype=[('sourcemodel', 'O'), ('headmodel', 'O'), ('elec', 'O'), ('reducerank', 'O'), ('checkpath', 'O'), ('showlogo', 'O'), ('outputfilepresent', 'O'), ('tracktimeinfo', 'O'), ('trackmeminfo', 'O'), ('toolbox', 'O'), ('callinfo', 'O'), ('version', 'O'), ('lbex', 'O'), ('sel50p', 'O'), ('feedback', 'O'), ('mollify', 'O'), ('patchsvd', 'O'), ('channel', 'O'), ('siunits', 'O'), ('backproject', 'O'), ('normalize', 'O'), ('normalizeparam', 'O'), ('weight', 'O'), ('senstype', 'O'), ('previous', 'O')]), array([[array([], shape=(0, 0), dtype=uint8),\n",
      "                array([], shape=(0, 0), dtype=uint8),\n",
      "                array([], shape=(0, 0), dtype=uint8), ...,\n",
      "                array([], shape=(0, 0), dtype=uint8),\n",
      "                array([], shape=(0, 0), dtype=uint8),\n",
      "                array([], shape=(0, 0), dtype=uint8)]], dtype=object), array([[array(['LPA'], dtype='<U3')],\n",
      "               [array(['RPA'], dtype='<U3')],\n",
      "               [array(['Nz'], dtype='<U2')],\n",
      "               [array(['Fp1'], dtype='<U3')],\n",
      "               [array(['Fpz'], dtype='<U3')],\n",
      "               [array(['Fp2'], dtype='<U3')],\n",
      "               [array(['AF9'], dtype='<U3')],\n",
      "               [array(['AF7'], dtype='<U3')],\n",
      "               [array(['AF5'], dtype='<U3')],\n",
      "               [array(['AF3'], dtype='<U3')],\n",
      "               [array(['AF1'], dtype='<U3')],\n",
      "               [array(['AFz'], dtype='<U3')],\n",
      "               [array(['AF2'], dtype='<U3')],\n",
      "               [array(['AF4'], dtype='<U3')],\n",
      "               [array(['AF6'], dtype='<U3')],\n",
      "               [array(['AF8'], dtype='<U3')],\n",
      "               [array(['AF10'], dtype='<U4')],\n",
      "               [array(['F9'], dtype='<U2')],\n",
      "               [array(['F7'], dtype='<U2')],\n",
      "               [array(['F5'], dtype='<U2')],\n",
      "               [array(['F3'], dtype='<U2')],\n",
      "               [array(['F1'], dtype='<U2')],\n",
      "               [array(['Fz'], dtype='<U2')],\n",
      "               [array(['F2'], dtype='<U2')],\n",
      "               [array(['F4'], dtype='<U2')],\n",
      "               [array(['F6'], dtype='<U2')],\n",
      "               [array(['F8'], dtype='<U2')],\n",
      "               [array(['F10'], dtype='<U3')],\n",
      "               [array(['FT9'], dtype='<U3')],\n",
      "               [array(['FT7'], dtype='<U3')],\n",
      "               [array(['FC5'], dtype='<U3')],\n",
      "               [array(['FC3'], dtype='<U3')],\n",
      "               [array(['FC1'], dtype='<U3')],\n",
      "               [array(['FCz'], dtype='<U3')],\n",
      "               [array(['FC2'], dtype='<U3')],\n",
      "               [array(['FC4'], dtype='<U3')],\n",
      "               [array(['FC6'], dtype='<U3')],\n",
      "               [array(['FT8'], dtype='<U3')],\n",
      "               [array(['FT10'], dtype='<U4')],\n",
      "               [array(['T9'], dtype='<U2')],\n",
      "               [array(['T7'], dtype='<U2')],\n",
      "               [array(['C5'], dtype='<U2')],\n",
      "               [array(['C3'], dtype='<U2')],\n",
      "               [array(['C1'], dtype='<U2')],\n",
      "               [array(['Cz'], dtype='<U2')],\n",
      "               [array(['C2'], dtype='<U2')],\n",
      "               [array(['C4'], dtype='<U2')],\n",
      "               [array(['C6'], dtype='<U2')],\n",
      "               [array(['T8'], dtype='<U2')],\n",
      "               [array(['T10'], dtype='<U3')],\n",
      "               [array(['TP9'], dtype='<U3')],\n",
      "               [array(['TP7'], dtype='<U3')],\n",
      "               [array(['CP5'], dtype='<U3')],\n",
      "               [array(['CP3'], dtype='<U3')],\n",
      "               [array(['CP1'], dtype='<U3')],\n",
      "               [array(['CPz'], dtype='<U3')],\n",
      "               [array(['CP2'], dtype='<U3')],\n",
      "               [array(['CP4'], dtype='<U3')],\n",
      "               [array(['CP6'], dtype='<U3')],\n",
      "               [array(['TP8'], dtype='<U3')],\n",
      "               [array(['TP10'], dtype='<U4')],\n",
      "               [array(['P9'], dtype='<U2')],\n",
      "               [array(['P7'], dtype='<U2')],\n",
      "               [array(['P5'], dtype='<U2')],\n",
      "               [array(['P3'], dtype='<U2')],\n",
      "               [array(['P1'], dtype='<U2')],\n",
      "               [array(['Pz'], dtype='<U2')],\n",
      "               [array(['P2'], dtype='<U2')],\n",
      "               [array(['P4'], dtype='<U2')],\n",
      "               [array(['P6'], dtype='<U2')],\n",
      "               [array(['P8'], dtype='<U2')],\n",
      "               [array(['P10'], dtype='<U3')],\n",
      "               [array(['PO9'], dtype='<U3')],\n",
      "               [array(['PO7'], dtype='<U3')],\n",
      "               [array(['PO5'], dtype='<U3')],\n",
      "               [array(['PO3'], dtype='<U3')],\n",
      "               [array(['PO1'], dtype='<U3')],\n",
      "               [array(['POz'], dtype='<U3')],\n",
      "               [array(['PO2'], dtype='<U3')],\n",
      "               [array(['PO4'], dtype='<U3')],\n",
      "               [array(['PO6'], dtype='<U3')],\n",
      "               [array(['PO8'], dtype='<U3')],\n",
      "               [array(['PO10'], dtype='<U4')],\n",
      "               [array(['O1'], dtype='<U2')],\n",
      "               [array(['Oz'], dtype='<U2')],\n",
      "               [array(['O2'], dtype='<U2')],\n",
      "               [array(['O9'], dtype='<U2')],\n",
      "               [array(['Iz'], dtype='<U2')],\n",
      "               [array(['O10'], dtype='<U3')],\n",
      "               [array(['T3'], dtype='<U2')],\n",
      "               [array(['T5'], dtype='<U2')],\n",
      "               [array(['T4'], dtype='<U2')],\n",
      "               [array(['T6'], dtype='<U2')],\n",
      "               [array(['M1'], dtype='<U2')],\n",
      "               [array(['M2'], dtype='<U2')],\n",
      "               [array(['A1'], dtype='<U2')],\n",
      "               [array(['A2'], dtype='<U2')]], dtype=object), array(['{pos}_chan_ori'], dtype='<U14'))                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  ]],\n",
      "      dtype=[('dim', 'O'), ('pos', 'O'), ('unit', 'O'), ('inside', 'O'), ('cfg', 'O'), ('leadfield', 'O'), ('label', 'O'), ('leadfielddimord', 'O')])}\n"
     ]
    }
   ],
   "source": [
    "import scipy.io\n",
    "import pandas as pd\n",
    "\n",
    "leadfield_mat = scipy.io.loadmat('leadfield_fem_eeg_cfg.sourcemodel.mat')\n",
    "print(leadfield_mat)\n",
    "#df = pd.DataFrame(leadfield_mat)\n",
    "#df.to_csv('leadfield_matrix.csv')\n",
    "#leadfield_mat = {k:v for k, v in leadfield_mat.items() if k[0] != '_'}\n",
    "#data = pd.DataFrame({k: pd.Series(v[0]) for k, v in leadfield_mat.items()}) \n",
    "\n",
    "#data.to_csv(\"eadfield_matrix.csv\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "NI_Project",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.4"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}