diff --git a/ISBI2018_Abstract/Slepians.aux b/ISBI2018_Abstract/Slepians.aux index 83ce3c0..a147d72 100644 --- a/ISBI2018_Abstract/Slepians.aux +++ b/ISBI2018_Abstract/Slepians.aux @@ -1,19 +1,21 @@ \relax \@writefile{toc}{\contentsline {section}{\numberline {1} Introduction}{1}} \newlabel{sec:intro}{{1}{1}} \@writefile{toc}{\contentsline {section}{\numberline {2} Methods}{1}} \newlabel{sec:methods}{{2}{1}} +\@writefile{toc}{\contentsline {subsection}{\numberline {2.1} Graph signal processing basics}{1}} \citation{C2} \bibstyle{IEEEbib} \bibdata{strings,refs} \bibcite{Lamp86}{1} \bibcite{C2}{2} +\@writefile{toc}{\contentsline {subsection}{\numberline {2.2} Slepian vectors}{2}} \@writefile{toc}{\contentsline {section}{\numberline {3} Results}{2}} \newlabel{sec:results}{{3}{2}} \@writefile{toc}{\contentsline {subsection}{\numberline {3.1} Simulated data}{2}} \@writefile{toc}{\contentsline {subsection}{\numberline {3.2} Real fMRI data}{2}} \@writefile{toc}{\contentsline {section}{\numberline {4} Discussion}{2}} \newlabel{sec:discussion}{{4}{2}} \@writefile{lof}{\contentsline {figure}{\numberline {1}{\ignorespaces Example of placing a figure with experimental results.}}{2}} \newlabel{fig:res}{{1}{2}} \@writefile{toc}{\contentsline {section}{\numberline {5} References}{2}} diff --git a/ISBI2018_Abstract/Slepians.log b/ISBI2018_Abstract/Slepians.log index 9b9f538..6ea38d4 100644 --- a/ISBI2018_Abstract/Slepians.log +++ b/ISBI2018_Abstract/Slepians.log @@ -1,320 +1,319 @@ -This is pdfTeX, Version 3.14159265-2.6-1.40.18 (TeX Live 2017) (preloaded format=pdflatex 2017.8.24) 12 OCT 2017 17:15 +This is pdfTeX, Version 3.14159265-2.6-1.40.18 (TeX Live 2017) (preloaded format=pdflatex 2017.8.24) 12 OCT 2017 18:50 entering extended mode restricted \write18 enabled. file:line:error style messages enabled. %&-line parsing enabled. **Slepians.tex (./Slepians.tex LaTeX2e <2017-04-15> Babel <3.10> and hyphenation patterns for 84 language(s) loaded. (/usr/local/texlive/2017/texmf-dist/tex/latex/base/article.cls Document Class: article 2014/09/29 v1.4h Standard LaTeX document class (/usr/local/texlive/2017/texmf-dist/tex/latex/base/size10.clo File: size10.clo 2014/09/29 v1.4h Standard LaTeX file (size option) ) \c@part=\count79 \c@section=\count80 \c@subsection=\count81 \c@subsubsection=\count82 \c@paragraph=\count83 \c@subparagraph=\count84 \c@figure=\count85 \c@table=\count86 \abovecaptionskip=\skip41 \belowcaptionskip=\skip42 \bibindent=\dimen102 ) (./spconf.sty) (/usr/local/texlive/2017/texmf-dist/tex/latex/amsmath/amsmath.sty Package: amsmath 2016/11/05 v2.16a AMS math features \@mathmargin=\skip43 For additional information on amsmath, use the `?' option. (/usr/local/texlive/2017/texmf-dist/tex/latex/amsmath/amstext.sty Package: amstext 2000/06/29 v2.01 AMS text (/usr/local/texlive/2017/texmf-dist/tex/latex/amsmath/amsgen.sty File: amsgen.sty 1999/11/30 v2.0 generic functions \@emptytoks=\toks14 \ex@=\dimen103 )) (/usr/local/texlive/2017/texmf-dist/tex/latex/amsmath/amsbsy.sty Package: amsbsy 1999/11/29 v1.2d Bold Symbols \pmbraise@=\dimen104 ) (/usr/local/texlive/2017/texmf-dist/tex/latex/amsmath/amsopn.sty Package: amsopn 2016/03/08 v2.02 operator names ) \inf@bad=\count87 LaTeX Info: Redefining \frac on input line 213. \uproot@=\count88 \leftroot@=\count89 LaTeX Info: Redefining \overline on input line 375. \classnum@=\count90 \DOTSCASE@=\count91 LaTeX Info: Redefining \ldots on input line 472. LaTeX Info: Redefining \dots on input line 475. LaTeX Info: Redefining \cdots on input line 596. \Mathstrutbox@=\box26 \strutbox@=\box27 \big@size=\dimen105 LaTeX Font Info: Redeclaring font encoding OML on input line 712. LaTeX Font Info: Redeclaring font encoding OMS on input line 713. \macc@depth=\count92 \c@MaxMatrixCols=\count93 \dotsspace@=\muskip10 \c@parentequation=\count94 \dspbrk@lvl=\count95 \tag@help=\toks15 \row@=\count96 \column@=\count97 \maxfields@=\count98 \andhelp@=\toks16 \eqnshift@=\dimen106 \alignsep@=\dimen107 \tagshift@=\dimen108 \tagwidth@=\dimen109 \totwidth@=\dimen110 \lineht@=\dimen111 \@envbody=\toks17 \multlinegap=\skip44 \multlinetaggap=\skip45 \mathdisplay@stack=\toks18 LaTeX Info: Redefining \[ on input line 2817. LaTeX Info: Redefining \] on input line 2818. ) (/usr/local/texlive/2017/texmf-dist/tex/latex/graphics/graphicx.sty Package: graphicx 2014/10/28 v1.0g Enhanced LaTeX Graphics (DPC,SPQR) (/usr/local/texlive/2017/texmf-dist/tex/latex/graphics/keyval.sty Package: keyval 2014/10/28 v1.15 key=value parser (DPC) \KV@toks@=\toks19 ) (/usr/local/texlive/2017/texmf-dist/tex/latex/graphics/graphics.sty Package: graphics 2017/04/14 v1.1b Standard LaTeX Graphics (DPC,SPQR) (/usr/local/texlive/2017/texmf-dist/tex/latex/graphics/trig.sty Package: trig 2016/01/03 v1.10 sin cos tan (DPC) ) (/usr/local/texlive/2017/texmf-dist/tex/latex/graphics-cfg/graphics.cfg File: graphics.cfg 2016/06/04 v1.11 sample graphics configuration ) Package graphics Info: Driver file: pdftex.def on input line 99. (/usr/local/texlive/2017/texmf-dist/tex/latex/graphics-def/pdftex.def File: pdftex.def 2017/01/12 v0.06k Graphics/color for pdfTeX (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/infwarerr.sty Package: infwarerr 2016/05/16 v1.4 Providing info/warning/error messages (HO) ) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/ltxcmds.sty Package: ltxcmds 2016/05/16 v1.23 LaTeX kernel commands for general use (HO) ) \Gread@gobject=\count99 )) \Gin@req@height=\dimen112 \Gin@req@width=\dimen113 ) (/usr/local/texlive/2017/texmf-dist/tex/latex/amsfonts/amssymb.sty Package: amssymb 2013/01/14 v3.01 AMS font symbols (/usr/local/texlive/2017/texmf-dist/tex/latex/amsfonts/amsfonts.sty Package: amsfonts 2013/01/14 v3.01 Basic AMSFonts support \symAMSa=\mathgroup4 \symAMSb=\mathgroup5 LaTeX Font Info: Overwriting math alphabet `\mathfrak' in version `bold' (Font) U/euf/m/n --> U/euf/b/n on input line 106. )) (./Slepians.aux) \openout1 = `Slepians.aux'. -LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 24. -LaTeX Font Info: ... okay on input line 24. -LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 24. -LaTeX Font Info: ... okay on input line 24. -LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 24. -LaTeX Font Info: ... okay on input line 24. -LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 24. -LaTeX Font Info: ... okay on input line 24. -LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 24. -LaTeX Font Info: ... okay on input line 24. -LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 24. -LaTeX Font Info: ... okay on input line 24. -LaTeX Font Info: Try loading font information for OT1+ptm on input line 24. +LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 26. +LaTeX Font Info: ... okay on input line 26. +LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 26. +LaTeX Font Info: ... okay on input line 26. +LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 26. +LaTeX Font Info: ... okay on input line 26. +LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 26. +LaTeX Font Info: ... okay on input line 26. +LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 26. +LaTeX Font Info: ... okay on input line 26. +LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 26. +LaTeX Font Info: ... okay on input line 26. +LaTeX Font Info: Try loading font information for OT1+ptm on input line 26. (/usr/local/texlive/2017/texmf-dist/tex/latex/psnfss/ot1ptm.fd File: ot1ptm.fd 2001/06/04 font definitions for OT1/ptm. ) (/usr/local/texlive/2017/texmf-dist/tex/context/base/mkii/supp-pdf.mkii [Loading MPS to PDF converter (version 2006.09.02).] \scratchcounter=\count100 \scratchdimen=\dimen114 \scratchbox=\box28 \nofMPsegments=\count101 \nofMParguments=\count102 \everyMPshowfont=\toks20 \MPscratchCnt=\count103 \MPscratchDim=\dimen115 \MPnumerator=\count104 \makeMPintoPDFobject=\count105 \everyMPtoPDFconversion=\toks21 ) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/pdftexcmds.sty Package: pdftexcmds 2017/03/19 v0.25 Utility functions of pdfTeX for LuaTeX (HO ) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/ifluatex.sty Package: ifluatex 2016/05/16 v1.4 Provides the ifluatex switch (HO) Package ifluatex Info: LuaTeX not detected. ) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/ifpdf.sty Package: ifpdf 2017/03/15 v3.2 Provides the ifpdf switch ) Package pdftexcmds Info: LuaTeX not detected. Package pdftexcmds Info: \pdf@primitive is available. Package pdftexcmds Info: \pdf@ifprimitive is available. Package pdftexcmds Info: \pdfdraftmode found. ) (/usr/local/texlive/2017/texmf-dist/tex/latex/oberdiek/epstopdf-base.sty Package: epstopdf-base 2016/05/15 v2.6 Base part for package epstopdf (/usr/local/texlive/2017/texmf-dist/tex/latex/oberdiek/grfext.sty Package: grfext 2016/05/16 v1.2 Manage graphics extensions (HO) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/kvdefinekeys.sty Package: kvdefinekeys 2016/05/16 v1.4 Define keys (HO) )) (/usr/local/texlive/2017/texmf-dist/tex/latex/oberdiek/kvoptions.sty Package: kvoptions 2016/05/16 v3.12 Key value format for package options (HO) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/kvsetkeys.sty Package: kvsetkeys 2016/05/16 v1.17 Key value parser (HO) (/usr/local/texlive/2017/texmf-dist/tex/generic/oberdiek/etexcmds.sty Package: etexcmds 2016/05/16 v1.6 Avoid name clashes with e-TeX commands (HO) Package etexcmds Info: Could not find \expanded. (etexcmds) That can mean that you are not using pdfTeX 1.50 or (etexcmds) that some package has redefined \expanded. (etexcmds) In the latter case, load this package earlier. ))) Package epstopdf-base Info: Redefining graphics rule for `.eps' on input line 4 38. Package grfext Info: Graphics extension search list: (grfext) [.png,.pdf,.jpg,.mps,.jpeg,.jbig2,.jb2,.PNG,.PDF,.JPG,.JPE G,.JBIG2,.JB2,.eps] (grfext) \AppendGraphicsExtensions on input line 456. (/usr/local/texlive/2017/texmf-dist/tex/latex/latexconfig/epstopdf-sys.cfg File: epstopdf-sys.cfg 2010/07/13 v1.3 Configuration of (r)epstopdf for TeX Liv e )) LaTeX Font Info: Font shape `OT1/ptm/bx/n' in size <12> not available -(Font) Font shape `OT1/ptm/b/n' tried instead on input line 27. -LaTeX Font Info: Try loading font information for U+msa on input line 27. +(Font) Font shape `OT1/ptm/b/n' tried instead on input line 29. +LaTeX Font Info: Try loading font information for U+msa on input line 29. (/usr/local/texlive/2017/texmf-dist/tex/latex/amsfonts/umsa.fd File: umsa.fd 2013/01/14 v3.01 AMS symbols A ) -LaTeX Font Info: Try loading font information for U+msb on input line 27. +LaTeX Font Info: Try loading font information for U+msb on input line 29. (/usr/local/texlive/2017/texmf-dist/tex/latex/amsfonts/umsb.fd File: umsb.fd 2013/01/14 v3.01 AMS symbols B ) LaTeX Font Info: Font shape `OT1/ptm/bx/n' in size <10> not available -(Font) Font shape `OT1/ptm/b/n' tried instead on input line 29. +(Font) Font shape `OT1/ptm/b/n' tried instead on input line 31. LaTeX Font Info: Font shape `OT1/ptm/bx/it' in size <10> not available -(Font) Font shape `OT1/ptm/b/it' tried instead on input line 33. - -Underfull \hbox (badness 1132) in paragraph at lines 62--63 -\OT1/ptm/m/it/10 quency \OT1/ptm/m/n/10 struc-tural pat-terns on the graph (i.e -., for which - [] - -[1{/usr/local/texlive/2017/texmf-var/fonts/map/pdftex/updmap/pdftex.map} +(Font) Font shape `OT1/ptm/b/it' tried instead on input line 35. + [1{/usr/local/texlive/2017/texmf-var/fonts/map/pdftex/updmap/pdftex.map} ] +Overfull \hbox (4.24074pt too wide) in paragraph at lines 80--81 +[]\OT1/ptm/m/n/10 If we de-note the new set by the ma-trix $\OT1/cmr/bx/n/10 S +\OT1/cmr/m/n/10 = [[]\OMS/cmsy/m/n/10 j[]j [] j[]\OT1/cmr/m/n/10 ]$\OT1/ptm/m/n +/10 , + [] + Package epstopdf Info: Source file: (epstopdf) date: 2015-09-21 14:44:34 (epstopdf) size: 381618 bytes (epstopdf) Output file: (epstopdf) date: 2017-10-12 12:15:57 (epstopdf) size: 12642 bytes (epstopdf) Command: -(epstopdf) \includegraphics on input line 93. +(epstopdf) \includegraphics on input line 115. Package epstopdf Info: Output file is already uptodate. - File: image1-eps-converted-to.pdf Graphic file (type pdf) -Package pdftex.def Info: image1-eps-converted-to.pdf used on input line 93. +Package pdftex.def Info: image1-eps-converted-to.pdf used on input line 115. (pdftex.def) Requested size: 241.84842pt x 56.43285pt. Package epstopdf Info: Source file: (epstopdf) date: 2015-09-21 14:44:34 (epstopdf) size: 387374 bytes (epstopdf) Output file: (epstopdf) date: 2017-10-12 12:15:57 (epstopdf) size: 13074 bytes (epstopdf) Command: -(epstopdf) \includegraphics on input line 100. +(epstopdf) \includegraphics on input line 122. Package epstopdf Info: Output file is already uptodate. File: image2-eps-converted-to.pdf Graphic file (type pdf) -Package pdftex.def Info: image2-eps-converted-to.pdf used on input line 100. +Package pdftex.def Info: image2-eps-converted-to.pdf used on input line 122. (pdftex.def) Requested size: 113.81102pt x 70.17142pt. Package epstopdf Info: Source file: (epstopdf) date: 2015-09-21 14:44:34 (epstopdf) size: 229254 bytes (epstopdf) Output file: (epstopdf) date: 2017-10-12 12:15:58 (epstopdf) size: 10135 bytes (epstopdf) Command: -(epstopdf) \includegraphics on input line 107. +(epstopdf) \includegraphics on input line 129. Package epstopdf Info: Output file is already uptodate. File: image3-eps-converted-to.pdf Graphic file (type pdf) -Package pdftex.def Info: image3-eps-converted-to.pdf used on input line 107. +Package pdftex.def Info: image3-eps-converted-to.pdf used on input line 129. (pdftex.def) Requested size: 113.81102pt x 81.7108pt. (./Slepians.bbl) [2 <./image1-eps-converted-to.pdf> <./image2-eps-converted-to .pdf> <./image3-eps-converted-to.pdf>] (./Slepians.aux) ) Here is how much of TeX's memory you used: - 2512 strings out of 492995 - 33366 string characters out of 6132704 - 98878 words of memory out of 5000000 - 6000 multiletter control sequences out of 15000+600000 + 2514 strings out of 492995 + 33379 string characters out of 6132704 + 99892 words of memory out of 5000000 + 6002 multiletter control sequences out of 15000+600000 23915 words of font info for 59 fonts, out of 8000000 for 9000 1141 hyphenation exceptions out of 8191 - 39i,10n,26p,810b,283s stack positions out of 5000i,500n,10000p,200000b,80000s + 39i,10n,26p,1268b,283s stack positions out of 5000i,500n,10000p,200000b,80000s {/usr/local/texlive/2017/texmf-dist/fonts/enc/dvips/base/8r.e nc} < /usr/local/texlive/2017/texmf-dist/fonts/type1/public/amsfonts/cm/cmmi10.pfb> -Output written on Slepians.pdf (2 pages, 181044 bytes). +Output written on Slepians.pdf (2 pages, 193939 bytes). PDF statistics: 104 PDF objects out of 1000 (max. 8388607) 73 compressed objects within 1 object stream 0 named destinations out of 1000 (max. 500000) 16 words of extra memory for PDF output out of 10000 (max. 10000000) diff --git a/ISBI2018_Abstract/Slepians.pdf b/ISBI2018_Abstract/Slepians.pdf index bd5c6b0..c58e287 100644 Binary files a/ISBI2018_Abstract/Slepians.pdf and b/ISBI2018_Abstract/Slepians.pdf differ diff --git a/ISBI2018_Abstract/Slepians.synctex.gz b/ISBI2018_Abstract/Slepians.synctex.gz index bf2ffec..17c695a 100644 Binary files a/ISBI2018_Abstract/Slepians.synctex.gz and b/ISBI2018_Abstract/Slepians.synctex.gz differ diff --git a/ISBI2018_Abstract/Slepians.tex b/ISBI2018_Abstract/Slepians.tex index 65cf5e9..e34faf9 100755 --- a/ISBI2018_Abstract/Slepians.tex +++ b/ISBI2018_Abstract/Slepians.tex @@ -1,133 +1,155 @@ % Template for ISBI-2018 paper; to be used with: % spconf.sty - ICASSP/ICIP LaTeX style file, and % IEEEbib.bst - IEEE bibliography style file. % -------------------------------------------------------------------------- \documentclass{article} \usepackage{spconf,amsmath,graphicx} \usepackage{amssymb} +\DeclareMathOperator*{\argmin}{\arg\!\min} + % Example definitions. % -------------------- \def\x{{\mathbf x}} \def\L{{\cal L}} % Title. % ------ \title{Slepians: a generalized graph Fourier transform to resolve localised functional brain interactions} % % Single address. % --------------- \name{Thomas AW Bolton$^{\star \dagger \ddagger}$, Younes Farouj$^{\star \dagger \ddagger}$, Dimitri Van De Ville$^{\star \dagger}$\thanks{$^{\ddagger}$ Both authors contributed equally.}} \address{$^{\star}$ Institute of Bioengineering, Ecole Polytechnique Federale de Lausanne (EPFL), Switzerland \\ $^{\dagger}$ Department of Radiology and Medical Informatics, University of Geneva (UNIGE), Switzerland} \begin{document} %\ninept % \maketitle % \begin{abstract} Neuroscience research is boring, because it lacks exotic mathematical tricks. Here, we introduce such a convoluted approach in the form of Slepians. Neuroscience research is boring, because it lacks exotic mathematical tricks. Here, we introduce such a convoluted approach in the form of Slepians. Neuroscience research is boring, because it lacks exotic mathematical tricks. Here, we introduce such a convoluted approach in the form of Slepians. Neuroscience research is boring, because it lacks exotic mathematical tricks. Here, we introduce such a convoluted approach in the form of Slepians. Neuroscience research is boring, because it lacks exotic mathematical tricks. Here, we introduce such a convoluted approach in the form of Slepians. \end{abstract} % \begin{keywords} Graph Fourier transform, Slepians, Graph signal processing \end{keywords} % \section{Introduction} \label{sec:intro} % How can we study the brain from the structural point of view ? Magnetic resonance imaging (MRI) has opened a number of avenues for the study of the brain. At the structural level, diffusion MRI enables to resolve the physical connectivity that exists between different regions. Viewed from the prism of graph analysis, this information is embedded in an adjacency matrix, which can then be decomposed into a set of canonical structural elements, commonly termed the \textit{eigenmodes}. % What about the functional side ? We would like to combine both aspects ! In parallel, functional MRI offers a window on brain activity over time, and on its dynamics at rest or upon cognitive challenge. There is an exquisite relationship between brain function and the underlying structural scaffold; for this reason, bimodal analytical approaches that can combine those two pieces of information are particularly tailored, and an emerging topic of interest. % Describe more precisely that in GSP, we express the signal as a combination of those building blocks In particular, graph signal processing (GSP) has recently gained momentum for this purpose. In this framework, functional information is regarded as temporal signals on a graph defined from structural measurements, and viewed as a linear combination of eigenmodes. This has the double advantage of enabling the study of structure/function relationships, and of permitting a wide set of signal processing operations in the spectral domain to improve the quality of functional signals. % Why Slepians can be nice: more accurate local decomposition that is also robust to changes in the chosen subset of nodes In such approaches, a similar importance is typically assigned to each brain region (i.e., each node of the studied network) at the decomposition stage. However, in some settings, it may be desirable to enhance the accuracy of the analysis on a subset of particularly important areas (for example, the brain regions expected to respond to a given paradigm). At the same time, some flexibility is also desired when it comes to defining this subset, so that rough prior knowledge be enough information to provide. % Summary of what is to come Here, we show how this can be achieved through a generalised decomposition into Slepian vectors. On simulated data, we show that if provided with a set of nodes of interest, noisy functional signals can be recovered more accurately within this set compared to a standard decomposition. In addition, we show how error estimates are robust to changes in the selected subset. On a real data example, we exemplify the potential of Slepians to reveal subtler, localised interaction patterns in the context of a visual stimulation task. \section{Methods} \label{sec:methods} -% Eigenmodes (Laplacian eigendecomposition). -Let a graph $\mathcal{G} = (\mathcal{V},\mathbf{W})$ characterised by the set of $N$ nodes $\mathcal{V}$, linked as described by the symmetrical adjacency matrix $\mathbf{W} \in \mathbb{R}^{N \times N}$. In the brain application considered here, $w_{i,j}$ will be large if brain regions $i$ and $j$ are strongly physically connected. The resulting Laplacian matrix $\mathbf{L} = \mathbf{D} - \mathbf{W}$, where $\mathbf{D}$ is the diagonal degree matrix with $d_{i,i} = \sum_j w_{i,j}$, admits an eigendecomposition as $\mathbf{L} = \mathbf{V} \mathbf{\Lambda} \mathbf{V^{\top}}$. The matrix $\mathbf{V} = [\mathbf{v_1}|\mathbf{v_2}|\ldots|\mathbf{v_N}]$ contains the eigenmodes as its columns, arranged in ascending eigenvalue order $\lambda_1 < \lambda_2 < \ldots < \lambda_N$. +\subsection{Graph signal processing basics} + +% Introducing the eigenmodes +Let a graph $\mathcal{G} = (\mathcal{V},\mathbf{W})$ characterised by the set of $N$ nodes $\mathcal{V}$, linked as described by the symmetrical adjacency matrix $\mathbf{W} \in \mathbb{R}^{N\times N}$. In the brain application considered here, $w_{i,j}$ will be large if brain regions $i$ and $j$ are strongly physically connected. The resulting Laplacian matrix $\mathbf{L} = \mathbf{D} - \mathbf{W}$, where $\mathbf{D}$ is the diagonal degree matrix with $d_{i,i} = \sum_j w_{i,j}$, admits an eigendecomposition as $\mathbf{L} = \mathbf{V} \mathbf{\Lambda} \mathbf{V^{\top}}$. The matrix $\mathbf{V} = [\mathbf{v_1}|\mathbf{v_2}|\ldots|\mathbf{v_N}]$ contains the eigenmodes as its columns, arranged in ascending eigenvalue order $\lambda_1 < \lambda_2 < \ldots < \lambda_N$. They satisfy $\mathbf{v_k^{\top}} \mathbf{v_l} = \delta_{k,l}$. Because $\lambda_k = \mathbf{v_k}^{\top} \mathbf{L} \mathbf{v_k} = \sum_{i\neq j} W_{i,j} ([v_k]_i - [v_k]_j)^2$, eigenmodes of smaller eigenvalues will represent \textit{low frequency} structural patterns on the graph (i.e., for which strongly connected nodes show similar values), while eigenmodes of larger eigenvalues will denote less organised patterns with respect to the graph structure. + +% Introducing the GFT +From this description, a signal $\mathbf{x} \in \mathbb{R}^{N\times 1}$ can conveniently be expressed in the graph domain as $\mathbf{\hat{x}} = \mathbf{V^{\top}} \mathbf{x}$, which is known as the \textit{graph Fourier transform}; conversely, we also have $\mathbf{x} = \mathbf{V} \mathbf{\hat{x}}$. Each element of $\mathbf{\hat{x}}$ then represents the strength with which an eigenmode is contributing to the signal at hand. In addition, note that the eigenmodes are also solutions to the Laplacian embedding problem, where the goal is to find a mapping of the graph nodes on a line so that connected ones stay as close as possible, that is, to find $\mathbf{x}$ as: + +\begin{equation} +\mathbf{x^{*}} = \argmin_\mathbf{x}{\mathbf{x^{\top}} \mathbf{L} \mathbf{x}} = \argmin_\mathbf{x}{\mathbf{x^{\top}} \mathbf{V} \mathbf{\Lambda} \mathbf{V^{\top}} \mathbf{x}}, +\end{equation} + +with $\mathbf{x^{\top}} \mathbf{x} = 1$ and $\mathbf{x^{\top}} \mathbf{1} = 0$. -% Link to the classical graph cut problem -Because $\lambda_k = \mathbf{v_k}^{\top} \mathbf{L} \mathbf{v} = \sum_{i\neq j} W_{i,j} ([v_k]_i - [v_k]_j)^2$, eigenmodes of smaller eigenvalues will represent \textit{low frequency} structural patterns on the graph (i.e., for which strongly connected nodes show similar values), while eigenmodes of larger eigenvalues will denote less organised patterns with respect to the graph structure. Note that the eigenmodes are also solutions to the Laplacian embedding problem, where the goal is to find a mapping of the graph nodes on a line so that connected ones stay as close as possible, that is, to minimise $\mathbf{x^{\top}} \mathbf{L} \mathbf{x}$ with $\mathbf{x^{\top}} \mathbf{x} = 1$ and $\mathbf{x^{\top}} \mathbf{1} = 0$. +\subsection{Slepian vectors} -% Slepians -Slepians: a generalised GFT concept, $\mu$ and $\xi$. +% Slepians: (1) concentration, (2) generalised Laplacian embedding +The goal is to derive an alternative set of basis vectors, under the constraint that their energy should be localised within a predefined subset of nodes $\mathcal{S}$. Futher, we want this alternative set to be derived from a bandwidth-limited subset of original eigenvectors, with dimension $K\leq N$. To highlight the selected nodes, we define $\mathbf{M}$ as the diagonal matrix with $m_{i,i} = 1$ if a node is included, and $0$ otherwise. Further, we define $\mathbf{V_T} \in \mathbb{R}^{N \times K}$ as the trimmed set of eigenmodes. + + +If we denote the new set by the matrix $\mathbf{S} = [\mathbf{s_1}|\mathbf{s_2}|\ldots|\mathbf{s_M}]$, and using the equality $\mathbf{\Lambda} = \mathbf{\Lambda}^{1/2} \mathbf{V^{\top}} \mathbf{V} \mathbf{\Lambda}^{1/2}$, Slepian vectors $\mathbf{s_k}$, $k = 1,\ldots,M$ can be seen as the solutions of a generalised Laplacian embedding formulation: + +\begin{equation} +\mathbf{s^{*}} = \argmin_\mathbf{s}{\mathbf{s^{\top}} \mathbf{V} \mathbf{\Lambda}^{1/2} \mathbf{V^{\top}} \mathbf{S} \mathbf{V} \mathbf{\Lambda}^{1/2} \mathbf{V^{\top}} \mathbf{s}}. +\end{equation} + +Indeed, we retrieve Equation (1) if $\mathbf{V_T} = \mathbf{V}$ and $\mathbf{S} = \mathbf{I}$. % Processing functional signals on the structural basis + \section{Results} \label{sec:results} \subsection{Simulated data} Quantification of optimal retrieval error for signals corrupted with added noise. Probed across noise levels, across ????? \subsection{Real fMRI data} HCP example subject on which $W$ and $X$ were computed; visual task, comparison in our subset of interest between GFT and Slepians. \section{Discussion} \label{sec:discussion} % Below is an example of how to insert images. Delete the ``\vspace'' line, % uncomment the preceding line ``\centerline...'' and replace ``imageX.ps'' % with a suitable PostScript file name. % ------------------------------------------------------------------------- \begin{figure}[htb] \begin{minipage}[b]{1.0\linewidth} \centering \centerline{\includegraphics[width=8.5cm]{image1}} % \vspace{2.0cm} \centerline{(a) Result 1}\medskip \end{minipage} % \begin{minipage}[b]{.48\linewidth} \centering \centerline{\includegraphics[width=4.0cm]{image2}} % \vspace{1.5cm} \centerline{(b) Results 3}\medskip \end{minipage} \hfill \begin{minipage}[b]{0.48\linewidth} \centering \centerline{\includegraphics[width=4.0cm]{image3}} % \vspace{1.5cm} \centerline{(c) Result 4}\medskip \end{minipage} % \caption{Example of placing a figure with experimental results.} \label{fig:res} % \end{figure} This is an example reference \cite{C2}. % To start a new column (but not a new page) and help balance the last-page % column length use \vfill\pagebreak. % ------------------------------------------------------------------------- \vfill \pagebreak % References should be produced using the bibtex program from suitable % BiBTeX files (here: strings, refs, manuals). The IEEEbib.bst bibliography % style file from IEEE produces unsorted bibliography list. % ------------------------------------------------------------------------- \bibliographystyle{IEEEbib} \bibliography{strings,refs} \end{document}