手把手教你画漂亮的功能连接矩阵图（附代码）

《本文同步发布于“脑之说”微信公众号，欢迎搜索关注~~》

科研图片的重要性是非常重要的，审稿人看到一个低质量的或者欠美观的图片，很可能影响文章的审稿。本期我将介绍一种画功能连接矩阵热力图的程序。该程序使用简单，有了改代码再也不用担心类似下图的矩阵热力图。
以下几幅图都是来自一些神经、精神、心理学的*期刊，包括nature medicine，nature communications、brain等。这些图是以矩热力图的形式展示的功能连接或者统计学结果图。我们可以看到属于不同的网络的节点的连接被一些网格线分割开，很容易让别人识别脑网络。




那么这种图是怎么做出来的呢？本期我们将给大家介绍一个专门用于做这种矩阵热力图的程序，并教大家如何使用改程序做出高清的用于直接投稿用的图片。
以下图片是使用笔者改程序做出来的矩阵热力图，这个图反映的是sub001的功能连接。

这个功能连接是一个114114的矩阵，表示有114个节点。另外我们有一个叫做netIndex.mat的文件，里面储存的是114个节点的网络标签，即每一个节点属于第几个网络。netIndex.mat是一个1114的向量，比如第1个位置的元素值表示第1个节点属于哪个脑网络。

进入具体操作正题：
第1步，下载代码和例子数据：
代码下载地址（进入github项目后可以下载我的整个项目，也可以只下载该代码，最简单的方式就是直接把代码复制到matlab新建的函数里面）：https://github.com/lichao312214129/lc_rsfmri_tools_matlab/blob/master/Workstation/code_workstation2018_dynamicFC/visualization/lc_netplot.m
数据下载地址（将其中的netIndex.mat和sub001.mat下载下来）: https://github.com/lichao312214129/lc_rsfmri_tools_matlab/tree/master/Workstation/code_workstation2018_dynamicFC/visualization/toy_data
第2步，整理下载的文件：
将下载的两个.mat文件（netIndex.mat 和sub001.mat）放在一个叫做toy_data的文件夹下。并将下载的代码（lc_netplot.m）放在与toy_data同级的文件下。
第3步，运行代码：
代码setpath到matlab路径后，在matlab命令窗口输入”lc_netplot”后，程序会在窗口打印出两个使用的例子，分别是EXAMPLE1和EXAMPLE 2.复制其中一个例子，到窗口，回车即可运行代码。运行结束后，图片会保存到toy_data文件夹下。

注：代码使用非常简单，另外代码里面有详细的使用信息，大家可以自行修改。该代码使用brewermap来生成网络legends的colormap。
代码作者黎超，如有使用问题欢迎咨询。
邮箱：[email protected]或者[email protected]