需要订阅 JoVE 才能查看此. 登录或开始免费试用。

Summary

提出了一种用于检查小脑灰质形态测量的标准化管道。该管道结合了高分辨率，最先进的方法，用于优化和自动化的小脑包裹和基于体素的小脑配准，以进行体积量化。

Abstract

多项研究为小脑在各种认知和情感功能中的作用提供了令人信服的证据，远远超出了其与运动控制的历史关联。结构和功能性神经影像学研究进一步完善了对小脑功能性神经解剖学的理解，超越了其解剖划分，突出了检查健康变异性和神经系统疾病中单个小脑亚基的必要性。本文提出了一种用于检查小脑灰质形态测量的标准化管道，该管道结合了高分辨率，最先进的方法，用于优化和自动化小脑包裹（使用U-Net局部约束优化的自动小脑解剖包裹;ACAPULCO）和基于体素的小脑注册（空间无偏倚的下层模板;SUIT）用于体积定量。

该管道广泛适用于一系列神经系统疾病，并且是全自动的，仅对输出进行质量控制所需的手动干预。该管道是免费提供的，并附有大量文档，并且可以在Mac，Windows和Linux操作系统上运行。该管道应用于Friedreich共济失调（FRDA）患者的队列，并提供代表性结果以及组水平推断统计分析的建议。该管道可以促进整个领域的可靠性和可重复性，最终为表征和跟踪神经系统疾病中的小脑结构变化提供强大的方法学方法。

Introduction

小脑是大脑的一部分，历史上与运动控制^1，2，3有关，并且被认为仅与一小部分罕见疾病（例如遗传性共济失调⁴）整体相关。然而，来自非人类灵长类动物解剖学追踪研究以及人类病变和神经影像学研究的趋同研究，为小脑在各种认知^5，6，⁷，情感⁸，^9，10，11和其他非运动功能^中的作用提供了令人信服的证据7^，12（见⁶ 以供审查）。此外，小脑异常越来越多地与广泛的神经和精神疾病有关，包括帕金森病¹³，阿尔茨海默病^14，15，癫痫^16，17，精神分裂症¹⁸和自闭症谱系障碍¹⁹.因此，将小脑纳入人类脑部疾病和规范行为变异的功能和结构模型中变得至关重要。

在解剖学上，小脑可以沿其上下轴分为三个叶:前叶，后叶和絮状硬叶。波瓣进一步细分为10个小叶，用罗马数字I-X^20，21表示（图1）。小脑也可以分为中线（蚯蚓）和侧（半球）区域，分别接收来自脊髓和大脑皮层的输入。前叶，包括小叶I-V，传统上与运动过程相关，并且与脑运动皮层²²具有相互连接。后叶，包括小叶VI-IX，主要与非运动过程¹¹相关，并且与前额叶皮层，后壁和颞上脑皮层^8，23具有相互连接。最后，由小叶X组成的絮状硬核与前庭核具有相互连接，前庭核在姿势和步态²¹期间控制眼球运动和身体平衡。

最近越来越多的使用功能性神经影像学的工作进一步完善了对小脑功能性神经解剖学的理解，超越了其解剖结构。例如，静息状态功能性磁共振成像（fMRI）技术已被用于绘制小脑和大脑之间功能相互作用的模式²⁴。此外，使用基于任务的包裹方法，King及其同事⁷证明小脑在其广度上显示出丰富而复杂的功能特化模式，这与各种运动，情感，社交和认知任务相关的不同功能边界证明了这一点。总的来说，这些研究强调了检查单个小脑亚基的重要性，以发展小脑参与健康变异性和以小脑结构和/或功能改变为特征的神经系统疾病的完整生物学特征。

本研究的重点是使用结构MRI在人体中量化小脑体积局部变化的方法。一般来说，使用MRI数据量化区域脑体积有两种基本方法:基于特征 的分割 和 基于体素的配准。基于特征的分割方法使用解剖学特征和标准化地图集来自动确定次区域之间的边界。用于细分的主流软件包包括FreeSurfer²⁵，BrainSuite²⁶和 FSL-FIRST²⁷。然而，这些包仅提供小脑的粗包裹（例如，标记每个半球的整个灰质和整个白质），从而忽略了单个小脑小叶。这些方法也容易出现误断，特别是周围脉管系统的过度整合。

已经开发了新的机器学习和多图集标记算法，它们提供了更准确和更细粒度的小脑包裹，包括使用隐式多边界进化的小脑小脑自动分类算法（ACCLAIM^28，29），小脑分析工具包（CATK³⁰），多个自动生成的模板（MAGeT³¹），人类小脑及其小脑的快速自动分割（RASCAL^32）），图形切割分割³³和小脑分割（CERES³⁴）。在最近的一篇论文中，比较了最先进的全自动小脑包裹方法，发现CERES2相对于小脑叶的黄金标准手动分割³⁵优于其他方法。最近，Han及其同事³⁶开发了一种名为ACAPULCO（使用U-Net进行局部约束优化的自动小脑解剖包裹）的深度学习算法，该算法与CERES2相当，对健康和萎缩的小脑具有广泛的适用性，以开源Docker和Singularity容器格式用于”现成”实现，并且比其他方法更具时效性。ACAPULCO自动将小脑包裹成28个解剖区域。

与基于特征的分割相比，基于体素的配准方法是将MRI精确映射到模板图像。要实现此映射，原始图像中的体素必须在大小和形状上失真。这种失真的幅度有效地提供了相对于黄金标准模板的每个体素的体积度量。这种形式的体积评估被称为”基于体素的形态测量”³⁷。基于全脑体素的配准方法，如FSL-FLIRT³⁸/FNIRT³⁹、SPM统一分割⁴⁰和CAT12⁴¹，通常用于基于体素的形态测量。然而，这些方法不能很好地解释小脑，导致在颅下区域（小脑，脑干⁴²）的可靠性和有效性较差。为了解释这些局限性，开发了SUIT（空间无偏倚下层模板）算法来优化小脑配准并提高基于体素的形态测量^42，43的准确性。

基于特征的分割和基于体素的配准方法用于估计区域小脑体积具有基本的优缺点。分割方法对于量化解剖学上定义的区域（例如，小叶³⁵）的体积要准确得多。然而，小脑不同功能模块之间的边界不会映射到其解剖叶和裂缝上（相当于大脑的回和沟⁷）。由于基于配准的方法不受解剖学特征的约束，因此可以对小脑进行更细粒度的空间推断和高维结构功能映射⁴⁴。总而言之，细分和注册方法是相辅相成的，可用于回答不同的研究问题。

在这里，提出了一个新的标准化管道，它集成了这些现有的，经过验证的方法，以提供优化和自动化的小脑包裹（ACAPULCO）和基于体素的配准（SUIT），用于体积量化（图2）。该管道建立在既定方法的基础上，包括质量控制协议，使用定性可视化和定量异常值检测，以及使用Freesurfer获得颅内体积（ICV）估计的快速方法。管道是完全自动化的，只需检查质量控制输出就需要手动干预，并且可以在Mac，Windows和Linux操作系统上运行。该管道是免费提供的，没有限制其用于非商业目的，并且可以在完成简短的注册表⁴⁵后从ENIGMA联盟成像协议网页（在”ENIGMA小脑体积管道”下）访问。

所有必需的软件都列在 材料表中，除了下面描述的协议外，下载管道时还可以获得详细的教程，包括现场演示。最后，提供了代表性的结果，来自Friedreich共济失调（FRDA）患者队列和年龄和性别匹配的健康对照组的管道实施，以及组级统计推断分析的建议。

Protocol

注:本研究中使用的数据是莫纳什大学人类研究伦理委员会批准的项目（项目7810）的一部分。与会者提供了书面知情同意书。虽然管道可以在Mac，Windows或Linux操作系统上运行，但ACAPULCO，SUIT和QC管道已经在Linux（Ubuntu）和Mac（Catalina，Big Sur v11.0.1）操作系统上进行了明确的测试。

1. 模块 1: ACAPULCO （解剖学包裹）

1. 数据采集
   1. 以1 mm 3或更低的分辨率收集整个大脑的^{3D T1} 加权MRI图像。建议使用各向同性体素尺寸（通常为 1 mm x 1 mm x 1 mm）和 3-Telsa（或更大）扫描仪。请咨询其放射成像中心的成像专家，以设置和获取符合这些规范的数据。
      注意:T2加权图像有时可用于体积分析;但是，此处介绍的管道仅依赖于 T1 加权数据，并且所使用的某些工具是此类数据所独有的。因此，不能使用 T2 加权图像。
   2. 对图像进行视觉质量评估，以排除大额小脑畸形（例如，大病变）或大量运动伪影，这些运动伪像无法识别主要小脑标志物（例如，主要解剖裂）。不要自动排除萎缩的脑膜，即使很大。
   3. 对于小组研究，还要考虑使用免费提供的标准化工具（如MRISTHC⁴⁶）进行定量质量评估，以进一步识别有问题的数据。
   4. 使用 dcm2niix⁴⁷ 等工具将所有数据转换为 NIFTI-GZ 格式。
2. 建议的数据组织
   1. 获取 材料表中列出的所有必要软件。确保在运行管道之前安装了 Docker⁴⁸ 或 Singularity⁴⁹、Matlab⁵⁰ 和^{SPM12 51}。
      注意:还提供了描述管道的大量书面和视频教程（请参阅 材料表）。
   2. 安装所有必要的软件后，在工作目录中创建文件夹，并将它们标记为”acapulco”，”suit”和”freesurfer”。从命令行使用 mkdir 命令执行此操作。
   3. 在”acapulco”目录中，创建一个 输出 文件夹。在 输出 文件夹中，为研究中的每个受试者创建一个目录，其中包含NIFTI-GZ格式的T1加权图像。
      注意:建议在其他位置保留原始数据的副本。
3. 使用ACAPULCO解剖小脑包裹术
   1. 转到 材料表 并下载运行 ACAPULCO 所需的相关脚本和容器（在 acapulco 管道文件下）。在”acapulco”目录中，放置（i）ACAPULCO Docker OR Singularity容器（分别acapulco_0.2.1.tar.gz”或”.sif”），（ii）QC_scripts存档的内容（3个文件:”QC_Master.R”，”QC_Plots.Rmd”和”QC_Image_Merge.Rmd”），以及（iii）”R.sif”（奇点） 或“calculate_icv.tar”（docker）文件。
   2. 打开终端，然后从命令行在单个映像上运行 ACAPULCO 容器（替换<>如下所示）。等待约 5 分钟以完成处理。
      1. 使用 Docker，键入命令:
         docker load –input acapulco_0.2.1.tar.gz
         docker run -v $PWD:$PWD -w $PWD -t –user $（id -u）:$（id -g） –rm acapulco:latest -i output/<>/<>.nii.gz -o output/<>
      2. 使用奇点，键入命令:
         奇点运行 –cleanenv -B $PWD:$PWD acapulco-0.2.1.sif -i output/<>/<>.nii.gz -o output/<>
   3. 循环访问队列中的所有受试者/扫描。请参阅 材料表 ，获取用于下载管道的 ENIGMA 成像协议网站的链接（在 ENIGMA 小脑体积管道下）和教程手册，其中包含如何创建用于连续处理多个主题的 for-loop 的示例。
   4. 处理后，查找在特定于主题的文件夹中生成的以下文件:
      1. 识别”_n4_mni_seg_post_inverse.nii.gz”:原始（主题空间）中的包裹小脑面罩。
      2. 识别”_n4_mni_seg_post_volumes.csv”:阿卡普尔科产生的28个亚基中每个亚基的体积（以^{mm 3}为单位）;
      3. 识别代表性图像（在”pics”目录中）:矢状，轴向和冠状。
4. 统计异常值检测和质量控制 （QC）
   1. 从终端和”acapulco”目录中，确保QC_scripts的内容位于”acapulco”目录中。运行 QC 脚本:
      1. 使用 Docker，键入命令:
         docker load calculate_icv.tar
         docker run -v $PWD:$PWD -w $PWD –rm -it luhancheng/calculate_icv:latest Rscript
         QC_Master.R输出/
      2. 使用奇点，键入命令:
         奇点执行 -B $PWD:$PWD R.sif Rscript /path/to/QC_Master.R /path/to/acapulco/output
5. 检查由ACAPULCO生成的QC图像
   注意:ACAPULCO 包裹图像的质量检查有一个 3 步流程。
   1. 在Web浏览器中打开”QC_Images.html”，然后快速（每个受试者约10秒）滚动图像以识别明显的故障或系统性问题。记下失败或可疑的包裹图像的主题 ID，以便进行后续检查。
      注意:有关小脑小脑神经解剖学的指南，请参阅 图 3 ;有关”良好”包裹、”细微误包裹”和”全局故障”包裹的示例，请参见下面代表性结果部分中的图 4、 图 5 和 图 6 。
   2. 打开”Plots_for_Outliers.html”以选中定量统计异常值的复选框图。查找位于箱形图的胡须上方或下方的异常值（高于或低于平均值的 2.698 s.d）。将鼠标悬停在数据点上以显示”使用者 ID”。在”异常值”文件的相关列中，标识由”1″表示的异常值.csv，并记下”异常值.csv”中最后一列中标识为每个主题的异常值的段总数。
   3. 手动检查具有一个或多个异常值的每个图像。关键:使用标准的NIFTI图像查看器（例如，FSLEyes或MRICron），将ACAPULCO蒙版叠加到原始T1w图像上，以逐片检查包裹的质量。
      1. 要使用FSLEyes从命令行生成详细的QC的覆盖，i）将目录更改为”acapulco”目录，ii）指定要查看的主题（替换 ）:
         subj=
      2. 将以下代码复制/粘贴到终端（无需手动更改 {subj}，因为这是由上一行设置的:
         t1_image=output/${subj}/${subj}.nii.gz
         acapulco_image=output/${subj}/${subj}_n4_mni_seg_post_inverse.nii.gz
         fsleyes ${t1_image} ${acapulco_image} –overlayType label –lut random_big –outline –outlineWidth 3 ${acapulco_image} –overlayType volume –alpha 50 –cmap random
         注意:需要确定是否包括异常片段，即是否存在包裹错误，或者只是个体解剖结构中的正常变异性？每个包裹区域都是单独考虑的，因此可以为图像排除几个区域，而如果正确，可以保留其余区域。
      3. 是否需要从最终数据集中排除一个或多个宗地区域？
         如果是（已确认异常值），则通过在该主题的”Cerebel_vols.csv”文件的相应单元格中将体积估计值替换为 NA ，将此宗地从分析中排除。
      4. 包裹错误是否会导致某些小脑被排除在面罩之外？
         如果是（例如，如果面罩上缺少特定的小脑叶或出现”切断”），请立即将受试者排除在进一步的分析之外（即，不要继续对这些受试者运行SUIT模块）。

2. 模块 2: SUIT 小脑优化体素形态测量

1. 使用 SUIT 进行基于体素的形态测量分析
   关键:该管道要求ACAPULCO模块已经运行，因为它依赖于生成特定于主题的小脑面罩来优化小脑的配准和归一化到SUIT模板。如果ACA普尔科生成的特定受试者面罩不包括整个小脑，则有必要从SUIT模块中排除。有关单独运行 SUIT 的说明，请参阅⁵²。
   1. 获取 材料表中列出的所有必要软件。确保 SPM12 文件夹和所有子文件夹都在 MATLAB 路径中。确保enigma_suit脚本保存在 “spm12/toolbox” 目录中，并添加到 MATLAB 路径中。要检查 MATLAB 路径，请在 MATLAB 命令窗口中键入 pathtool ，然后单击 使用子文件夹添加 以添加相关文件夹。
   2. 为一个或多个主题运行 SUIT 管道。等待 ~15-20 分钟（如果使用图形用户界面 [GUI]）和 ~5-7 分钟（如果从终端（bash/shell）运行）以完成处理。
      1. 要使用 GUI（主题将串行运行），请在 MATLAB 命令窗口中键入命令:
         suit_enigma_all
      2. 在第一个弹出窗口中，从”acapulco/output”目录中选择要包含在分析中的主题文件夹。单击窗口右侧的各个文件夹，或右键单击并 全选。按 完成。在第二个弹出窗口中，选择将写入分析的 SUIT 目录。
      3. 或者，从 MATLAB 命令行为单个主题调用函数，键入以下命令:
         suit_enigma_all（’/path/to/acapulco/output/subjdir’，’/path/to/suitoutputdir’）
      4. 或者，通过键入以下命令，从 MATLAB 外部的终端窗口为单个主题调用函数:
         matlab -nodisplay -nosplash -r “suit_enigma_all（’/path/to/acapulco/output/subjdir’，’/path/to/suitoutputdir’）， exit”
   3. 请参阅 材料表 ，获取用于下载管道的 ENIGMA 成像协议网站的链接（在 ENIGMA 小脑体积管道下）和教程手册，其中包含如何创建用于连续处理多个主题的 for-loop 的示例。
   4. 查找有关脚本的以下几点。
      1. 确保脚本将 N4 偏置校正、MNI 对齐（刚体）T1 图像和 ACAPULCO 小脑掩模复制到输出目录中。
      2. 确保脚本分割小脑的灰质和白质。
      3. 确保脚本使用 ACAPULCO 掩码更正包裹中的过度包含错误。
      4. 确保脚本 DARTEL 使用 Jacobian 调制将数据规范化并重新切片到 SUIT 空间中，以便每个体素的值与其原始体积成正比。
      5. 检查每个主题的文件夹中是否有以下最终输出:”wd_seg1.nii”（灰质）和”wd_seg2.nii”（白质）。
2. 统计异常值检测和质量控制
   1. 目视检查归一化、调制图像 （wd*） 是否存在重大故障。在 MATLAB 中，键入以下命令:
      spm_display_4D
   2. 从套装子文件夹中手动选择”wd*seg1″图像，或导航到”suit”目录;在过滤器框中插入”^wd.*seg1″（无引号），然后按” Rec “按钮。按 完成。
   3. 滚动浏览图像以确保它们都对齐。参见 图7 ，了解来自健康对照组（A，B）和小脑严重萎缩个体（D）的正确归一化图像。
      注意:在这个阶段，受试者之间的解剖结构非常相似（因为它们已被注册到同一模板中），并且体积差异是通过不同的体素强度编码的。主要的失败将是显而易见的，例如，空白图像，大面积的缺失组织，不寻常的强度梯度（即，明亮的体素全部位于顶部，深色体素全部位于底部）。应从后续步骤中排除这些映像。
   4. 检查异常值的空间协方差。在 MATLAB 中，键入以下命令:
      check_spatial_cov
      1. 按照上一步选择”wd*seg1″图像。出现提示时，选择以下选项: 属性缩放: 是;要共变的变量:否;切片（毫米）: – 48 ， 间隙: 1.
      2. 查看显示每个影像相对于样本中所有其他影像的平均空间协方差的箱线图。识别在 MATLAB 命令窗口中均值下方>2s.d. 的数据点。对于这些，请检查 SUIT 文件夹中的”_n4_mni.nii.gz”图像是否存在伪影（运动、解剖异常）、图像质量问题或预处理错误。
      3. 如果图像质量和预处理是可以接受的，并且上一步中对调制图像的目视检查并不表示分割和归一化存在问题，请将这些数据保留在样本中。否则，请排除这些数据。

3. 模块 3（可选）:使用 FreeSurfer 估算颅内体积 （ICV）

注意:此模块将使用 FreeSurfer 管道来计算 ICV。如果同期群（任何版本）存在现有的自由冲浪者输出，则无需重新运行。

1. 设置免费冲浪者
   1. 确保免费冲浪者已下载并安装⁵³.转到材料 表 并下载相关脚本以运行此模块（在 ICV管道文件下）。使用 FreeSurfer 时，请设置以下变量:
      导出FREESURFER_HOME =
      source $FREESURFER_HOME/SetUpFreeSurfer.sh
   2. 将 替换为以下内容:
      导出SUBJECTS_DIR=/谜/自由冲浪者
2. 运行自由冲浪者自动回复1
   1. 对于单个主题，从”freesurfer”目录（处理时间约20分钟）中，键入命令:
      cd /谜/自由冲浪者
      侦察所有 -i ../input/.nii.gz -s -autorecon1
   2. 有关如何创建用于连续处理多个主题的 for 循环的示例，请参阅教程手册。
3. 智能网联汽车的计算
   1. 数据组织
      1. 在”freesurfer”目录中，放置（i）模块1中使用的Docker OR Singularity容器（分别为”calculate_icv.tar”或”R.sif”）和（ii）xfm2det脚本（参见 材料表）。然后，执行 git 克隆以克隆所需的 ICV 脚本:
         git clone https://github.com/Characterisation-Virtual-Laboratory/calculate_icv
   2. 运行ICV提取（处理时间〜5分钟）
      1. 从”freesurfer”目录，使用奇点（’R.sif’）容器，键入:
         singularity exec –cleanenv -B $PWD:$PWD R.sif calculate_icv/calculate_icv.py –freesurfer_dir=/path/to/freesurfer –acapulco_dir=/path/to/acapulco/QC/Cerebelvolsfile –output_csv_name=Cerebel_vols.csv calculate_icv
      2. 从”freesurfer”目录，使用 Docker 容器，键入:
         docker run -v $PWD:$PWD -w $PWD -rm -it luhancheng/calculate_icv:latest
         calculate_icv/calculate_icv.py –freesurfer_dir=/path/to/Freesurfer —
         acapulco_dir=/path/to/acapulco/QC/Cerebelvolsfile –output_csv_name=Cerebel_vols.csv
         calculate_icv
      3. 在没有容器的情况下运行脚本 – 请参阅 材料表 ，了解其他必需的软件和依赖项。从”免费冲浪者”目录中，键入:
         ./calculate_icv/ calculate_icv.py —freesurfer_dir=/path/to/freesurfer —
         acapulco_dir=/path/to/acapulco/QC/Cerebelvolsfile —
         output_csv_name=Cerebel_vols.csv calculate_icv
         注意:这将计算每个受试者的ICV，并在”Cerebel_vols.csv”文件的末尾附加一列ICV。

Representative Results

Discussion

小脑对广泛的人类运动³，认知⁵⁸，情感¹⁰和语言^7，59功能至关重要，并且与许多神经和精神疾病有关。用于量化区域小脑体积的标准化且易于实施的方法的可用性将有助于越来越详细的”全脑”结构功能映射，完整的疾病建模以及改善定义和跟踪小脑对脑疾病的贡献的机会。这里描述的这种标准化管道结合了最先进的方法，用于自动小脑包裹和更细粒度的小脑灰质形态测量的空间分析在健康和疾病中。

这里介绍的使用ACAPULCO的横截面小脑包裹分析结果表明，患有FRDA（与HC）的人的白质体积显着减少。这些发现支持了先前的FRDA研究，这些研究始终显示FRDA中早期，稳健和进行性的白质体积损失，特别是在齿状核中。此外，上小脑和下小脑蒂以及齿状核进行性神经变性的模式和程度已被证明是FRDA⁴⁴发病年龄的一个因素。SUIT分析的结果揭示了其他发现。具体而言，在对应于双侧小叶I-IV和右V的前叶区域的FRDA（与HC相比）中存在显着的体素水平体积损失，延伸至小叶VI。此外，SUIT分析显示，内侧后叶区域（包括小叶IX，X和Vermis）的FRDA（与HC相比）的体积损失显着。这种组间差异的模式与先前在同一队列FRDA患者中发表的工作相当，使用全脑VBM方法⁵⁵。

定义神经系统和精神疾病中的小脑异常是具有转化影响的高度优先研究领域。追踪和治疗神经系统疾病（特别是那些小脑是神经变性主要部位的神经系统疾病）是小脑受累的完整生物学特征的发展。本文介绍的管道允许个体小脑小叶灰质形态测量与临床测量之间的关系，这些测量被用作探索疾病临床终点的”金标准”。这样的研究可以产生重大的转化影响。例如，在罕见的小脑疾病领域，识别一个亚组患者小脑灰质萎缩的特定特征，该特征映射到或预测临床症状，这将对指导临床实践产生影响。SUIT模块的包含进一步允许解决有趣的研究问题，例如小脑的结构 – 功能映射或小脑⁶⁰的功能梯度分析。

组级统计分析的一般建议
ACAPULCO:每个受试者的每个小脑小叶（^以3毫米为单位）的体积记录在 Cerebel_vols.csv中。在对群体水平效应进行统计推断期间，应控制颅内体积（ICV;也记录在 Cerebel_vols.csv中），以考虑头部大小的变异性。应校正 α 显著性阈值，以考虑跨多个小叶的推断。

套装:灰质小脑VBM可以使用标准MRI处理软件（如SPM或FSL）在wdseg1.nii图像上进行。请参阅 CAT12 手册，了解使用 SPM12⁵⁴ 的 VBM 的精彩介绍。应控制ICV，以考虑头部尺寸的变化。

对于小脑中的 VBM，通常建议使用全宽不超过 3 mm 的半最大值 （FWHM） 的高斯空间平滑核。必须应用适当的统计校正来考虑体素之间的多重比较。通常，建议使用非参数方法（例如，SnPM或FSL-Randomise）。

使用ACAPULCO成功包裹小脑的最关键步骤是在处理前和处理后对T1图像进行一般质量检查。强烈建议用户检查对比度不好的图像（例如，图像上的梯度不一致）以及头部的严重倾斜和运动伪像，所有这些都会影响算法的性能。此外，虽然ACAPULCO算法已经在萎缩的小脑上进行了训练，但它还没有在病变数据上进行训练。预计大脑皮层中的病变预计不会影响算法的性能和随后的包裹准确性;然而，小脑的大梗死可能会产生包裹错误。对小脑面罩进行质量检查后处理至关重要。轻微的包裹错误（例如，小脑小脑小脑的轻微下计和过度内含）有时不被检测为统计异常值;相反，尽管存在明显的包裹错误，但如果数据在正常范围内，则可能会发生不正确的非异常值实例。如果受试者被确定为异常值，则必须逐片对小脑面罩进行后续的详细质量检查，以指导是否包括或排除该受试者的小叶的决策。运行 SUIT 管道（模块 2）时的另一个关键步骤是，它需要 ACAPULCO 模块已经运行。具体而言，SUIT要求在ACAPULCO中生产的小脑面罩进行小脑分离和分割。重要的是要对小脑面罩进行质量检查，以确保小脑完全覆盖。

该协议存在一些限制。首先，虽然ACAPURCO在小脑灰质包裹方面实现了最先进的精度，但它并没有针对白质包裹进行优化;髓质体覆盖了白质的主体，但不能提供所有白质的量度。其次，用于在ACAPULCO中定位和分割小脑的卷积神经网络不能很好地推广到具有不同对比度的图像或未在训练中使用的图像。例如，由于在训练中仅使用3T图像，因此使用在1.5 T扫描仪上获取的图像的包裹质量通常不如3T图像;此外，没有关于对这些图像进行的基本事实的统计数字。最后，通过提供ICV的估计值来控制头部大小对小脑体积估计的混杂效应的管道控制，该估计可以作为在组水平统计分析中没有兴趣的回归量。然而，理想的方法是在进行QC之前在个体水平上计算ICV校正的小脑体积，这样检测到的异常值反映了真正的包裹误差，而不是受试者神经解剖学中的自然变异性（例如，具有大头）。

总之，我们提出了一种用于检查小脑灰质形态测量的标准化管道，该管道对一系列神经系统疾病具有广泛的适用性。该管道的建立是为了允许进行大型，多站点研究和”大型分析”，并公开供研究小组使用，以促进整个领域的可靠性和可重复性。最终，该管道提供了一种强大的方法学方法，用于进一步表征和跟踪神经系统疾病进展的小脑结构变化。目前正在开发一条纵向管道。

Materials

ACAPULCO pipeline files 	0.2.1	http://enigma.ini.usc.edu/protocols/imaging-protocols/	Please make sure to use acapulco version 0.2.1
Docker for Mac		https://docs.docker.com/desktop/mac/install/	macOS must be version 10.14 or newer<br/> Docker requires sudo priviledges<br/> Docker imposes a memory (RAM) constraint on Mac OS. To increase the RAM, open Docker Desktop, go to Preferences and click on resources. Increase the Memory to the maximum
Docker for Windows		https://docs.docker.com/docker-for-windows/install/
ENIGMA SUIT scripts		http://enigma.ini.usc.edu/protocols/imaging-protocols/
FreeSurfer	7	https://surfer.nmr.mgh.harvard.edu/fswiki/DownloadAndInstall	Following variables need to be set everytime you work with Freesurfer:<br/> export FREESURFER_HOME=<img alt="equation" src="/files/ftp_upload/63340/lesser.png"/>freesurfer _installation_directory<img alt="equation" src="/files/ftp_upload/63340/greater.png"/><br/> source $FREESURFER_HOME/SetUpFreeSurfer.sh
export SUBJECTS_DIR=<img alt="equation" src="/files/ftp_upload/63340/lesser.png"/>path<img alt="equation" src="/files/ftp_upload/63340/greater.png"/>/enigma/Freesurfer
FSL (for FSLeyes). Optional	6	https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FslInstallation
ICV pipeline files		http://enigma.ini.usc.edu/protocols/imaging-protocols/	ICV pipeline can be run in two ways: 1) with docker/singularity. You will not require additionl software; 2) without docker/singularity- this involves running the ICV script (calculate_icv.py) manually. You will require the following additional software:<br/> Python version <img alt="equation" src="/files/ftp_upload/63340/greater.png"/>=3.5<br/> Python module pandas<br/> Python module fire<br/> Python module tabulate<br/> Python module Colorama
		https://github.com/Characterisation-Virtual-Laboratory/calculate_icv
MATLAB*	2019 or newer	https://au.mathworks.com/	An academic license is required
Singularity	3.7 or newer	https://www.sylabs.io/docs/	Prefered for high performance computing (HPC) clusters
SPM	12	http://www.fil.ion.ucl.ac.uk/spm/software/spm12/	Make sure spm12 and all subfolders are in your MATLAB path
SUIT Toolbox	3.4	http://www.diedrichsenlab.org/imaging/suit_download.htm	Make sure you place SUIT toolbox in spm12/toolbox directory
Troubleshooting manual and segmentation output examples		http://enigma.ini.usc.edu/protocols/imaging-protocols/
Tutorial manual and video		http://enigma.ini.usc.edu/protocols/imaging-protocols/	Manual and accompanying live demonstration provide detailed step-by-step instructions on how to run the pipeline from start to finish.
*Not freely available; an academic license is required