graph TB

%% structural
subgraph MRI structural
    tone(T1)
    surf(FreeSurfer surfaces)
    src(Source space)
    head(Scalp surface)
    tone   -->  surf
    surf   -->  src
    surf --> head

    %% bem
    flash(Flash5/30)
    bem(Conductor Model)
    surf --> bem
    flash --> bem

    surf --> fslab
    fslab("Cortical Labels<br/>(*.annot, *.label)")
end

%% forward
trans("MRI<->head<br>transformation")
info --> fwd
trans --> fwd
bem --> fwd
src --> fwd

%% importing
acq(Raw data from system)
acq   -->  raw
raw --> info
info("Measurement information")
raw(Preprocessed data)

%% events
events("Events")
raw     -->  events

%% epoching
epo(Epochs)
events --> epo
raw       -->  epo

%% evoked
evk(Evoked)
epo  -->  evk

%% covariance
cov(Covariance)
epo      -->  cov

%% make inverse
head-->trans
info-->trans
fwd(Forward)
inv(Inverse)
info --> inv
cov      -->  inv
fwd      -->  inv

%% apply inverse
stc(SourceEstimate)
evk       -->  stc
epo --> stc
inv       -->  stc
raw --> stc

%% label time course
ltc("Label time courses")
fslab --> ltc
stc --> ltc

subgraph Outcomes
  %% connectivity
  conn("Connectivity")
  decoding("Decoding (ML)")
  stats("Statistics")

end

ltc --> conn

%% ML
epo --> decoding

%% stats
ltc-->stats
stc --> stats
epo -->stats

%% %% %% %% %%
%% classes  %%
%% %% %% %% %%

%% red: freesurfer objects
classDef source font-weight:bold,fill:#fcc,stroke:#633,stroke-width:4px;
class bem,tone,flash,src,surf,head,fslab source;

%% Blue: python objects
classDef obj font-weight:bold,fill:#bce,stroke:#225,stroke-width:4px;
class acq,raw,events,info,trans,epo,evk,cov,inv,arr,fwd,stc,evarr,lab obj;
class ltc,decoding,stats,conn obj