U f{jf@s:dZddlZddlZddlZddZedkZer6edS)aYCompute SSP/PCA projections for EOG artifacts. Examples -------- .. code-block:: console $ mne compute_proj_eog -i sample_audvis_raw.fif -a \ --l-freq 1 --h-freq 35 \ --rej-grad 3000 --rej-mag 4000 --rej-eeg 100 or .. code-block:: console $ mne compute_proj_eog -i sample_audvis_raw.fif -a \ --l-freq 1 --h-freq 35 \ --rej-grad 3000 --rej-mag 4000 --rej-eeg 100 \ --proj sample_audvis_ecg-proj.fif to exclude ECG artifacts from projection computation. Nc%Csddlm}|t}|jdddddd|jd d d d d d|jddd ddd|jddddddd|jddddddd|jddddd dd|jd!d"d d#d$d|jd%d&d d'd(d|jd)d*d d+d$d|jd,d-d d.d/d|jd0d1d2d3d4d5|jd6d7d8d9d:d;d<|jd=d>d?d@d5|jdAdBddCdDd|jdEdFdGddHd$d|jdIdJd dKdLd|jdMdNd dOdPd|jdQdRd dSdTd|jdUdVd dWdXd|jdYdZd9d[d;d<|jd\d]d9d^d;d<|jd_d`dad@d5|jdbdcdddded|jdfdgdhd@d5|jdidjd dkdld|jdmdndodpdqd@d|\}}|j}|d@krP|td$|j }|j }|j }|j }|j } |j} |j} |j} |j} |j}|j}|j}|j}tdrt|jdst|jdtt|jdtt|jdu}|j}|j}|j}|j}|j }|j!}|j"}|j#}|d@k rFt$|dv}dwdx|%D}W5QRXt&dy|ng}|'dzsb|'d{rp|d@d|}n |d@d}}|d~}|r|d}n|d}t(j)j*||d} |d@k rt(j)*|}!n| }!d@}"t(j+j,| |!||||| | | |||||"||||| | ||d;d\}#}$| -|d@k r$|!-|d@k rHt&d|t(.||#}#t/|t0rlt1j23|rlt14|t&d|t(5||#t&d|t(6||$d@S)z Run command.r) get_optparserz-iz--inraw_inzInput raw FIF fileFILE)desthelpmetavarz--tmintminfloatzTime before event in secondsgɿ)rtyperdefaultz--tmaxtmaxzTime after event in secondsg?z-gz--n-gradn_gradintz&Number of SSP vectors for gradiometersz-mz--n-magn_magz'Number of SSP vectors for magnetometersz-ez--n-eegn_eegzNumber of SSP vectors for EEGz--l-freql_freqz"Filter low cut-off frequency in Hzz--h-freqh_freqz#Filter high cut-off frequency in Hz#z --eog-l-freq eog_l_freqz?Filter low cut-off frequency in Hz used for EOG event detectionz --eog-h-freq eog_h_freqz@Filter high cut-off frequency in Hz used for EOG event detection z-pz --preloadpreloadz7Temporary file used during computation (to save memory)T)rrr z-az --averageaverage store_truezCompute SSP after averagingF)ractionrr z--projprojz$Use SSP projections from a fif file.Nz --filtersize filter_lengthz#Number of taps to use for filteringiz-jz--n-jobsn_jobsz!Number of jobs to run in parallelz --rej-gradrej_gradzBGradiometers rejection parameter in fT/cm (peak to peak amplitude)iz --rej-magrej_magz@Magnetometers rejection parameter in fT (peak to peak amplitude)i z --rej-eegrej_eegu7EEG rejection parameter in µV (peak to peak amplitude)2z --rej-eogrej_eogu7EOG rejection parameter in µV (peak to peak amplitude)geAz --avg-refavg_refzAdd EEG average reference projz --no-projno_projz5Exclude the SSP projectors currently in the fiff filez--bad bad_fnamez5Text file containing bad channels list (one per line)z --event-idevent_idzID to use for eventsiz --event-rawraw_event_fnamez#raw file to use for event detectionz--tstarttstartz-Start artifact detection after tstart secondsgz-cz --channelch_namestringz&Custom EOG channel(s), comma separatedgvIh%<=gV瞯)ZgradmagZeegZeogrcSsg|] }|qS)rstrip).0wr/r/`/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mne/commands/mne_compute_proj_eog.py szrun..zBad channels read : %sz_raw.fifz-raw.fifiz _eog-eve.fifz_eog_avg-proj.fifz _eog-proj.fif)r)raw raw_eventrr r rrrrrrrrejectflatbadsr%r&r(rrr*r+copyz#Including SSP projections from : %szWriting EOG projections in %szWriting EOG events in %s)7Zmne.commands.utilsr__file__ add_option parse_argsr print_helpsysexitrr r rrrrrrrrrrdictr r r!r"r$r%r&r'r(rr)r*r+open readlinesprintendswithmneioZ read_raw_fifZ preprocessingZcompute_proj_eogcloseZ read_proj isinstancestrospathexistsremoveZ write_projZ write_events)%rparseroptionsargsrrr r rrrrrrrrrrr8r%r&r'r(Z proj_fnamer)r*r+Zfidr:prefixZeog_event_fnameZeog_proj_fnamer6r7r9Zprojseventsr/r/r3runs~                                           rU__main__)__doc__rLr@rGrU__name__Zis_mainr/r/r/r3s)