将 nearest.nii.gz 中标签值为 1 的像素提取出来

可以使用以下命令将 nearest.nii.gz 中标签值为 1 的像素提取出来，生成一个新的 NIfTI 文件：

fslmaths nearest.nii.gz -thr 1.5 -uthr 1.5 -bin 1_mask.nii.gz

其中，`-thr 1.5` 表示将值大于等于 1.5 的像素赋值为 1；`-uthr 1.5` 表示将值小于等于 1.5 的像素赋值为 1；`-bin` 表示将结果二值化，生成的 `1_mask.nii.gz` 文件中只有标签值为 1 的像素值为 1，其他像素值为 0。

相关问题

nearest.nii.gz 是将BN_Atlas图谱配准到个体b0图像上的结果，怎样将nearest.nii.gz 中图谱编号为1的脑区提取出来？

可以使用以下步骤将 nearest.nii.gz 中图谱编号为 1 的脑区提取出来：

1. 使用 FSLeyes 或类似的软件打开 BN_Atlas 中的 `BN_Atlas_246_1mm.nii.gz` 文件，找到编号为 1 的脑区的标签值。

2. 使用 FSL 的 `fslmaths` 命令将 nearest.nii.gz 中标签值为 1 的像素提取出来，生成一个新的 NIfTI 文件，命令如下：

   fslmaths nearest.nii.gz -thr 1.5 -uthr 1.5 -bin 1_mask.nii.gz

其中，`-thr 1.5` 表示将值大于等于 1.5 的像素赋值为 1；`-uthr 1.5` 表示将值小于等于 1.5 的像素赋值为 1；`-bin` 表示将结果二值化，生成的 `1_mask.nii.gz` 文件中只有标签值为 1 的像素值为 1，其他像素值为 0。

3. 使用 `fslmaths` 命令将原始图像与新生成的脑区掩膜相乘，提取出标签值为 1 的脑区，命令如下：

   fslmaths nearest.nii.gz -mul 1_mask.nii.gz BN_Atlas_1.nii.gz

其中，`-mul` 表示进行像素间的乘法运算，生成的 `BN_Atlas_1.nii.gz` 文件中只有标签值为 1 的脑区像素值保留，其他像素值为 0。

antsRegistrationSyN.sh -m brain.nii.gz -f brain_subjects.nii.gz -n 90 -t sr -o S1 antsRegistrationSyN.sh -m brain_subjects.nii.gz -f b0.nii.gz -n 90 -o S2 antsApplyTransforms -i BN.nii.gz -r brain_subjects.nii.gz -t S10GenericAffine.mat S11Warp.nii.gz -o inter.nii.gz -n NearestNeighbor antsApplyTransforms -i inter.nii.gz -r b0.nii.gz -t S20GenericAffine.mat S21Warp.nii.gz -o nearest.nii.gz -n NearestNeighbor

这是一个基于 ANTs 工具箱的脑部图像配准流程，可以将一个参考脑图像（brain.nii.gz）与多个脑图像（brain_subjects.nii.gz）进行配准。这个流程主要步骤包括：首先用 antsRegistrationSyN.sh 命令将参考脑图像与多个脑图像进行非线性配准，生成变换矩阵 S1。然后再用 antsRegistrationSyN.sh 命令将已配准的脑图像与 b0.nii.gz 进行线性配准，生成变换矩阵 S2。最后用 antsApplyTransforms 命令将 BN.nii.gz 按照变换矩阵 S1 和 S2 进行变换，最终得到经过配准后的 nearest.nii.gz 图像。其中，-n 90 表示使用 90 个正交步长进行配准，-t sr 表示使用 SyN 正则化方法进行非线性配准，-o 表示输出文件名。-n NearestNeighbor 表示插值方法采用最近邻插值。