antsRegistrationSyNQuick.sh运行成功但不出结果文件

问答讨论
1

antsRegistrationSyNQuick.sh -d 3 -f corr_T1.nii.gz -m corr_mag.nii.gz -o magtoT1

大家好，我使用这个命令计算配准，配准成功运行但是没有任何文件生成是为什么呢？求解答，谢谢啦

下面是运行过程(删减版)：

--------------------------------------------------------------------------------------
 Mapping parameters
--------------------------------------------------------------------------------------
 Dimensionality:           3
 Output name prefix:       magtoT1
 Fixed images:             /nhkj/run/corr_T1.nii.gz
 Moving images:            /nhkj/run/corr_mag.nii.gz
 Mask images:
 Initial transforms:
 Number of threads:        1
 Spline distance:          26
 SyN gradient step:        0.1
 Transform type:           s
 MI histogram bins:        32
 Precision:                d
 Use histogram matching:   0
 Repro                     0
======================================================================================
 antsRegistration call:
--------------------------------------------------------------------------------------
antsRegistration --verbose 1 --dimensionality 3 --float 0 --collapse-output-transforms 1 --output [ magtoT1,magtoT1Warped.nii.gz,magtoT1InverseWarped.nii.gz ] --interpolation Linear --use-histogram-matching 0 --winsorize-image-intensities [ 0.005,0.995 ] --initial-moving-transform [ /nhkj/run/corr_T1.nii.gz,/nhkj/run/corr_mag.nii.gz,1 ] --transform Rigid[ 0.1 ] --metric MI[ /nhkj/run/corr_T1.nii.gz,/nhkj/run/corr_mag.nii.gz,1,32,Regular,0.25 ] --convergence [ 1000x500x250x0,1e-6,10 ] --shrink-factors 8x4x2x1 --smoothing-sigmas 3x2x1x0vox --transform Affine[ 0.1 ] --metric MI[ /nhkj/run/corr_T1.nii.gz,/nhkj/run/corr_mag.nii.gz,1,32,Regular,0.25 ] --convergence [ 1000x500x250x0,1e-6,10 ] --shrink-factors 8x4x2x1 --smoothing-sigmas 3x2x1x0vox --transform SyN[ 0.1,3,0 ] --metric MI[ /nhkj/run/corr_T1.nii.gz,/nhkj/run/corr_mag.nii.gz,1,32] --convergence [ 100x70x50x0,1e-6,10 ] --shrink-factors 8x4x2x1 --smoothing-sigmas 3x2x1x0vox
--------------------------------------------------------------------------------------
All_Command_lines_OK
Using double precision for computations.
=============================================================================
The composite transform comprises the following transforms (in order):
  1. Center of mass alignment using fixed image: /nhkj/run/corr_T1.nii.gz and moving image: /nhkj/run/corr_mag.nii.gz (type = Euler3DTransform)
=============================================================================
  number of levels = 4
  number of levels = 4
  number of levels = 4
  fixed image: /nhkj/run/corr_T1.nii.gz
  moving image: /nhkj/run/corr_mag.nii.gz
  fixed image: /nhkj/run/corr_T1.nii.gz
  moving image: /nhkj/run/corr_mag.nii.gz
  fixed image: /nhkj/run/corr_T1.nii.gz
  moving image: /nhkj/run/corr_mag.nii.gz
Dimension = 3
Number of stages = 3
Use histogram matching = false
Winsorize image intensities = true
  Lower quantile = 0.005
  Upper quantile = 0.995


Stage 1 State
   Image metric = Mattes
   Fixed image = Image (0x55fcc21be240)
  RTTI typeinfo:   itk::Image<double, 3u>
  Reference Count: 2
  Modified Time: 1724
  Debug: Off
  Object Name:
  Observers:
   none
  Source: (none)
  Source output name: (none)
  Release Data: Off
  Data Released: False
  Global Release Data: Off
  PipelineMTime: 0
  UpdateMTime: 1507
  RealTimeStamp: 0 seconds
  LargestPossibleRegion:
   Dimension: 3
   Index: [0, 0, 0]
   Size: [208, 240, 170]
  BufferedRegion:
   Dimension: 3
   Index: [0, 0, 0]
   Size: [208, 240, 170]
  RequestedRegion:
   Dimension: 3
   Index: [0, 0, 0]
   Size: [208, 240, 170]
  Spacing: [1, 1, 1]
  Origin: [111.016, 113.315, -48.459]
  Direction:
-0.996475 -0.083883 -0.000604155
0.08388 -0.996468 0.00386228
-0.000926 0.00379799 0.999992

  IndexToPointMatrix:
-0.996475 -0.083883 -0.000604155
0.08388 -0.996468 0.00386228
-0.000926 0.00379799 0.999992

  PointToIndexMatrix:
-0.996475 0.08388 -0.000926
-0.083883 -0.996468 0.00379799
-0.000604155 0.00386228 0.999992

  Inverse Direction:
-0.996475 0.08388 -0.000926
-0.083883 -0.996468 0.00379799
-0.000604155 0.00386228 0.999992

  PixelContainer:
   ImportImageContainer (0x55fcc2212050)
   RTTI typeinfo:   itk::ImportImageContainer<unsigned long, double>
   Reference Count: 1
   Modified Time: 1504
   Debug: Off
   Object Name:
   Observers:
   none
   Pointer: 0x7f0e630fc010
   Container manages memory: true
   Size: 8486400
   Capacity: 8486400

   Moving image = Image (0x55fcc21c33e0)
 

*** Running Euler3DTransform registration ***


Total elapsed time: 142.2
 Registration finished. The antsRegistration call was:
--------------------------------------------------------------------------------------
2

看起来好像没什么问题，也许你可以试试加上-t 'a’的选项，看看在这种简单情况下，是否能够生成结果？

antsRegistrationSyNQuick.sh -d 3 -f corr_T1.nii.gz -m corr_mag.nii.gz -o magtoT1 -t 'a'
3

老师，还是显示能跑成功，但是仍然是没有任何文件生成:'(

4

那确实不太清楚了，我没有遇到过这种情况。我开始猜测也许是内存不足的问题，所以让你试试简单一些的选项。另外，你也许可以试试其他模态的数据，比如用T1配准MNI T1模板，虽然我觉得可能跟模态关系不大。

我也去ANTs的Github仓库的Issues（https://github.com/ANTsX/ANTs/issues）里搜索了一下，没有发现类似的情况，说明这很可能不是一个bug。也许你可以自己再去问一下开发者。