In [1]:

    

%matplotlib qt
import mne
import os
import numpy as np


    

In [2]:

    

eeg_dat = mne.io.read_raw_edf('./session1/a_1.edf', preload=True)


    

    




Extracting edf Parameters from ./session1/a_1.edf...
Setting channel info structure...
Interpolating stim channel. Events may jitter.
Creating Raw.info structure...
Reading 0 ... 301249  =      0.000 ...  1204.996 secs...






    




<ipython-input-2-cf38b9597152>:1: RuntimeWarning: Interpolating stim channel. Events may jitter.
  eeg_dat = mne.io.read_raw_edf('./session1/a_1.edf', preload=True)






    




Ready.

In [3]:

    

eeg_dat.ch_names


    

    
Out[3]:





[u'EEG FP1-REF',
 u'EEG FP2-REF',
 u'EEG F3-REF',
 u'EEG F4-REF',
 u'EEG C3-REF',
 u'EEG C4-REF',
 u'EEG P3-REF',
 u'EEG P4-REF',
 u'EEG O1-REF',
 u'EEG O2-REF',
 u'EEG F7-REF',
 u'EEG F8-REF',
 u'EEG T3-REF',
 u'EEG T4-REF',
 u'EEG T5-REF',
 u'EEG T6-REF',
 u'EEG A1-REF',
 u'EEG A2-REF',
 u'EEG FZ-REF',
 u'EEG CZ-REF',
 u'EEG PZ-REF',
 u'EEG ROC-REF',
 u'EEG LOC-REF',
 u'EEG EKG1-REF',
 u'EEG T1-REF',
 u'EEG T2-REF',
 u'PHOTIC-REF',
 u'IBI',
 u'BURSTS',
 'STI 014']

In [4]:

    

from collections import defaultdict
ch_types = defaultdict(str)
channel_names = []
channel_types = []

for ch in eeg_dat.ch_names:
    if 'EEG' in ch:
        ch_types[ch] = 'eeg'
    elif 'EMG' in ch:
        ch_types[ch] = 'emg'
    elif 'STI' in ch: 
        ch_types[ch] = 'stim'
    else:
        ch_types[ch] = 'misc'
        
eeg_dat.set_channel_types(ch_types)


# Example of one way to start fixing up the channels
#  Priority again is to keep everything consistent with MNE, so first check for MNE functions that can do 
#   things like this, and that can update this information insode the data.info object. 

from collections import defaultdict
ch_types = defaultdict(str)

channel_names = []
channel_types =[]

# lists to store the array names

for ch in eeg_dat.ch_names:    
    if 'EEG' in ch:
        ch_types[ch] = 'eeg'
        
    elif 'EMG' in ch:
        ch_types[ch] = 'emg'
        # Keeps track that this is an EMG channels
        
    elif 'STI' in ch:
        ch_types[ch] = 'stim'
        # Keeps track that this is an STI channel
    
    else:
        ch_types[ch] = 'misc'
        # Keeps track of MISC channels
        
eeg_dat.set_channel_types(ch_types)
        
for ch in eeg_dat.ch_names:    
    if 'EEG' in ch:
        ch = ch[4:7]
        ch =''.join(e for e in ch if e.isalnum())
        channel_names.append(ch)
        channel_types.append('eeg')
        ch_types[ch] = 'eeg'
        
    elif 'EMG' in ch:
        ch = ch[4:7]
        ch =''.join(e for e in ch if e.isalnum())
        channel_names.append(ch)
        channel_types.append('emg')
        ch_types[ch] = 'emg'
        # Keeps track that this is an EMG channels
        
    elif 'STI' in ch:
        channel_names.append(ch)
        channel_types.append('stim')
        ch_types[ch] = 'stim'
        # Keeps track that this is an STI channel
    
    else:
        channel_names.append(ch)
        channel_types.append('misc')
        ch_types[ch] = 'misc'
        # Keeps track of MISC channels
        
# The EEG channels use the standard naming strategy.
# By supplying the 'montage' parameter, approximate locations
montage = 'standard_1020'

info = mne.create_info(channel_names, eeg_dat.info['sfreq'], channel_types, montage)
print(info)


    

    




The following EEG sensors did not have a position specified in the selected montage: [u'ROC', u'LOC', u'EKG', u'T1', u'T2']. Their position has been left untouched.
<Info | 15 non-empty fields
    bads : list | 0 items
    ch_names : list | FP1, FP2, F3, F4, C3, C4, P3, P4, O1, ...
    chs : list | 30 items (EEG: 26, STIM: 1, MISC: 3)
    comps : list | 0 items
    custom_ref_applied : bool | False
    dev_head_t : 'mne.transforms.Transform | 3 items
    events : list | 0 items
    highpass : float | 0.0 Hz
    hpi_meas : list | 0 items
    hpi_results : list | 0 items
    lowpass : float | 125.0 Hz
    meas_date : numpy.ndarray | 1969-12-31 19:00:00
    nchan : int | 30
    projs : list | 0 items
    sfreq : float | 250.0 Hz
    acq_pars : NoneType
    acq_stim : NoneType
    buffer_size_sec : NoneType
    ctf_head_t : NoneType
    description : NoneType
    dev_ctf_t : NoneType
    dig : NoneType
    experimenter : NoneType
    file_id : NoneType
    hpi_subsystem : NoneType
    kit_system_id : NoneType
    line_freq : NoneType
    meas_id : NoneType
    proj_id : NoneType
    proj_name : NoneType
    subject_info : NoneType
    xplotter_layout : NoneType
>






    




<ipython-input-4-1d9a003652fe>:81: RuntimeWarning: The following EEG sensors did not have a position specified in the selected montage: [u'ROC', u'LOC', u'EKG', u'T1', u'T2']. Their position has been left untouched.
  info = mne.create_info(channel_names, eeg_dat.info['sfreq'], channel_types, montage)

In [5]:

    

eeg_dat.info


    

    
Out[5]:





<Info | 16 non-empty fields
    bads : list | 0 items
    buffer_size_sec : float | 1.0
    ch_names : list | EEG FP1-REF, EEG FP2-REF, EEG F3-REF, EEG F4-REF, ...
    chs : list | 30 items (EEG: 26, STIM: 1, MISC: 3)
    comps : list | 0 items
    custom_ref_applied : bool | False
    dev_head_t : 'mne.transforms.Transform | 3 items
    events : list | 0 items
    highpass : float | 0.0 Hz
    hpi_meas : list | 0 items
    hpi_results : list | 0 items
    lowpass : float | 125.0 Hz
    meas_date : int | 1369747734
    nchan : int | 30
    projs : list | 0 items
    sfreq : float | 250.0 Hz
    acq_pars : NoneType
    acq_stim : NoneType
    ctf_head_t : NoneType
    description : NoneType
    dev_ctf_t : NoneType
    dig : NoneType
    experimenter : NoneType
    file_id : NoneType
    hpi_subsystem : NoneType
    kit_system_id : NoneType
    line_freq : NoneType
    meas_id : NoneType
    proj_id : NoneType
    proj_name : NoneType
    subject_info : NoneType
    xplotter_layout : NoneType
>

In [6]:

    

eeg_dat.plot()


    

    
Out[6]:

In [7]:

    

eeg_dat.plot_psd()


    

    




tmax defaults to 60. in 0.14 but will change to np.inf in 0.15. Set it explicitly to avoid this warning
Effective window size : 8.192 (s)






    




<ipython-input-7-650886116877>:1: DeprecationWarning: tmax defaults to 60. in 0.14 but will change to np.inf in 0.15. Set it explicitly to avoid this warning
  eeg_dat.plot_psd()
/usr/local/lib/python2.7/dist-packages/mne/viz/raw.py:614: DeprecationWarning: In version 0.15 average will default to False and spatial_colors will default to True.
  'spatial_colors will default to True.', DeprecationWarning)






    
Out[7]:

In [8]:

    

notches = np.arange(60, 61, 60)
eeg_dat.notch_filter(notches)
print(eeg_dat.get_data(start=0).shape)
filtered = eeg_dat.copy().filter(1, 70, h_trans_bandwidth=10)
print(filtered.get_data(start=0).shape)


    

    




Setting up band-stop filter from 59 - 61 Hz
Filter length of 3300 samples (13.200 sec) selected
(30, 301250)
Setting up band-pass filter from 1 - 70 Hz
l_trans_bandwidth chosen to be 1.0 Hz
Filter length of 1650 samples (6.600 sec) selected
(30, 301250)

In [9]:

    

filtered.plot_psd()
ch_names = filtered.info['ch_names']
picks = filtered.pick_channels(ch_names=ch_names)


    

    




tmax defaults to 60. in 0.14 but will change to np.inf in 0.15. Set it explicitly to avoid this warning
Effective window size : 8.192 (s)






    




<ipython-input-9-1667abd226ea>:1: DeprecationWarning: tmax defaults to 60. in 0.14 but will change to np.inf in 0.15. Set it explicitly to avoid this warning
  filtered.plot_psd()
/usr/local/lib/python2.7/dist-packages/mne/viz/raw.py:614: DeprecationWarning: In version 0.15 average will default to False and spatial_colors will default to True.
  'spatial_colors will default to True.', DeprecationWarning)

In [10]:

    

filtered.ch_names
data = filtered.get_data(start=0, stop=250)

print(data.shape)


    

    




(30, 250)

In [11]:

    

data2 = mne.time_frequency.stft(data, wsize=140, tstep=2)

print(data2.shape)


    

    




Number of frequencies: 71
Number of time steps: 125
(30, 71, 125)

In [12]:

    

filtered.ch_names


    

    
Out[12]:





[u'EEG FP1-REF',
 u'EEG FP2-REF',
 u'EEG F3-REF',
 u'EEG F4-REF',
 u'EEG C3-REF',
 u'EEG C4-REF',
 u'EEG P3-REF',
 u'EEG P4-REF',
 u'EEG O1-REF',
 u'EEG O2-REF',
 u'EEG F7-REF',
 u'EEG F8-REF',
 u'EEG T3-REF',
 u'EEG T4-REF',
 u'EEG T5-REF',
 u'EEG T6-REF',
 u'EEG A1-REF',
 u'EEG A2-REF',
 u'EEG FZ-REF',
 u'EEG CZ-REF',
 u'EEG PZ-REF',
 u'EEG ROC-REF',
 u'EEG LOC-REF',
 u'EEG EKG1-REF',
 u'EEG T1-REF',
 u'EEG T2-REF',
 u'PHOTIC-REF',
 u'IBI',
 u'BURSTS',
 'STI 014']

In [13]:

    

print(np.real(data2))


    

    




[[[  3.62785821e-06   3.73357744e-06   3.83422639e-06 ...,  -1.82435841e-06
    -1.77196366e-06  -1.71837090e-06]
  [ -1.75781275e-06  -1.98661436e-06  -2.20745162e-06 ...,   1.19708173e-06
     1.06045008e-06   9.16291490e-07]
  [ -2.18998869e-07   8.75869719e-08   3.96174056e-07 ...,   8.47002970e-09
     4.08391233e-07   7.80058066e-07]
  ..., 
  [  1.79845342e-08   1.73459731e-08   1.60485218e-08 ...,   3.23797277e-09
     3.47693870e-09   3.19526013e-09]
  [ -1.80964815e-08  -1.78715607e-08  -1.75409018e-08 ...,  -3.20724167e-09
    -3.22393227e-09  -3.61555011e-09]
  [  1.80461999e-08   1.80480101e-08   1.79402733e-08 ...,   3.41099358e-09
     3.46684819e-09   3.10662502e-09]]

 [[ -1.58345649e-06  -1.62418619e-06  -1.66664426e-06 ...,  -1.79137534e-06
    -1.69434442e-06  -1.60063294e-06]
  [  1.16950680e-06   1.14353293e-06   1.10509528e-06 ...,   4.17074478e-07
     2.98583390e-07   1.81208357e-07]
  [ -1.60268219e-06  -1.51367487e-06  -1.36375206e-06 ...,   6.16267467e-07
     8.42815800e-07   1.03759001e-06]
  ..., 
  [ -6.93639074e-08  -6.60262136e-08  -6.26740568e-08 ...,  -3.88204341e-08
    -4.06472637e-08  -4.22115392e-08]
  [  6.85462632e-08   6.89013756e-08   6.63449871e-08 ...,   4.07141196e-08
     4.18672692e-08   4.14845511e-08]
  [ -6.94524428e-08  -6.83363326e-08  -6.93045964e-08 ...,  -4.17313352e-08
    -4.14314321e-08  -4.21648548e-08]]

 [[ -1.70164622e-07  -1.62354662e-07  -1.49549428e-07 ...,  -2.58304233e-07
    -2.39649876e-07  -2.15724394e-07]
  [  3.09288443e-07   3.09073167e-07   3.09185681e-07 ...,   4.66525129e-08
     1.75986850e-10  -3.96819615e-08]
  [ -2.74336149e-07  -3.08590273e-07  -3.30154700e-07 ...,   3.56526406e-07
     5.41771730e-07   7.08441188e-07]
  ..., 
  [ -8.73172394e-09  -8.27742216e-09  -7.64498913e-09 ...,  -1.94275099e-08
    -2.13613195e-08  -2.13070132e-08]
  [  8.44656995e-09   8.54606264e-09   8.45766483e-09 ...,   2.12540665e-08
     2.11587942e-08   2.21147041e-08]
  [ -8.71953269e-09  -8.53066143e-09  -8.42685448e-09 ...,  -2.13667198e-08
    -2.19824531e-08  -2.12923701e-08]]

 ..., 
 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]]

In [14]:

    

print(np.imag(data2))


    

    




[[[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  2.49696822e-06   2.41234307e-06   2.30195133e-06 ...,   1.37969972e-06
     1.45223053e-06   1.50832214e-06]
  [ -1.88082350e-06  -1.89895368e-06  -1.86343275e-06 ...,  -2.13670269e-06
    -2.06938913e-06  -1.92535692e-06]
  ..., 
  [ -1.99722870e-09  -5.14840126e-09  -8.12845373e-09 ...,   1.32473738e-09
     7.11055083e-10   1.82994268e-10]
  [  9.70920796e-10   2.59536993e-09   4.17425141e-09 ...,  -6.54509191e-10
    -3.72369798e-10  -2.80032884e-11]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [ -1.39418386e-07  -5.65006651e-08   2.03325409e-08 ...,   1.62655554e-06
     1.59867965e-06   1.55921961e-06]
  [ -3.73333657e-07  -6.93960512e-07  -9.92364817e-07 ...,  -1.26634651e-06
    -1.14711932e-06  -9.86651687e-07]
  ..., 
  [  6.24504898e-09   1.83578819e-08   2.99384626e-08 ...,  -1.53964104e-08
    -8.31382886e-09  -8.11760701e-10]
  [ -3.07641091e-09  -9.27067857e-09  -1.52638460e-08 ...,   7.86656278e-09
     4.16712804e-09   4.52730042e-10]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  4.37280826e-08   7.80265613e-08   1.12923956e-07 ...,   4.88941382e-07
     4.81300130e-07   4.68814264e-07]
  [  1.83804273e-07   1.38627849e-07   8.78394924e-08 ...,  -9.77372796e-07
    -8.85481700e-07  -7.56823864e-07]
  ..., 
  [  6.80296860e-10   2.21525632e-09   3.64021083e-09 ...,  -9.13597821e-09
    -5.50097880e-09  -1.68052753e-09]
  [ -3.57975387e-10  -1.11780500e-09  -1.87601714e-09 ...,   4.65325557e-09
     2.76119641e-09   8.14746507e-10]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 ..., 
 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]]

In [15]:

    

print(np.absolute(data2))


    

    




[[[  3.62785821e-06   3.73357744e-06   3.83422639e-06 ...,   1.82435841e-06
     1.77196366e-06   1.71837090e-06]
  [  3.05364634e-06   3.12506571e-06   3.18932949e-06 ...,   1.82662968e-06
     1.79820129e-06   1.76483024e-06]
  [  1.89353044e-06   1.90097253e-06   1.90508149e-06 ...,   2.13671947e-06
     2.10930197e-06   2.07737571e-06]
  ..., 
  [  1.80950931e-08   1.80938889e-08   1.79896306e-08 ...,   3.49848493e-09
     3.54890153e-09   3.20049593e-09]
  [  1.81225089e-08   1.80590317e-08   1.80307407e-08 ...,   3.27334407e-09
     3.24536570e-09   3.61565856e-09]
  [  1.80461999e-08   1.80480101e-08   1.79402733e-08 ...,   3.41099358e-09
     3.46684819e-09   3.10662502e-09]]

 [[  1.58345649e-06   1.62418619e-06   1.66664426e-06 ...,   1.79137534e-06
     1.69434442e-06   1.60063294e-06]
  [  1.17778761e-06   1.14492790e-06   1.10528231e-06 ...,   1.67917660e-06
     1.62632366e-06   1.56971407e-06]
  [  1.64559054e-06   1.66517050e-06   1.68659646e-06 ...,   1.40833912e-06
     1.42345397e-06   1.43180815e-06]
  ..., 
  [  6.96444707e-08   6.85308157e-08   6.94575333e-08 ...,   4.17621307e-08
     4.14887912e-08   4.22193438e-08]
  [  6.86152644e-08   6.95222629e-08   6.80782073e-08 ...,   4.14671236e-08
     4.20741392e-08   4.14870214e-08]
  [  6.94524428e-08   6.83363326e-08   6.93045964e-08 ...,   4.17313352e-08
     4.14314321e-08   4.21648548e-08]]

 [[  1.70164622e-07   1.62354662e-07   1.49549428e-07 ...,   2.58304233e-07
     2.39649876e-07   2.15724394e-07]
  [  3.12364348e-07   3.18770085e-07   3.29161974e-07 ...,   4.91162022e-07
     4.81300162e-07   4.70490672e-07]
  [  3.30218615e-07   3.38298148e-07   3.41640018e-07 ...,   1.04036948e-06
     1.03807247e-06   1.03666353e-06]
  ..., 
  [  8.75818513e-09   8.56872676e-09   8.46740774e-09 ...,   2.14684475e-08
     2.20582578e-08   2.13731838e-08]
  [  8.45415226e-09   8.61885576e-09   8.66322888e-09 ...,   2.17574844e-08
     2.13381999e-08   2.21297074e-08]
  [  8.71953269e-09   8.53066143e-09   8.42685448e-09 ...,   2.13667198e-08
     2.19824531e-08   2.12923701e-08]]

 ..., 
 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]

 [[  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  ..., 
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]
  [  0.00000000e+00   0.00000000e+00   0.00000000e+00 ...,   0.00000000e+00
     0.00000000e+00   0.00000000e+00]]]

In [ ]:

    




    

In [ ]: