ANTsX / ANTsR

R interface to the ANTs biomedical image processing library
https://antsx.github.io/ANTsR
Apache License 2.0
127 stars 35 forks source link

fitTransformToPairedPoints diffeo example doesn't work #363

Closed SlicerMorph closed 2 years ago

SlicerMorph commented 2 years ago

I am trying to learn the new fitTransformToPairedPoints function introduced to ANTsR. I am following the examples and all works fine except the diffeo one. Below is the error message and the sessionInfo().

>  domainImage <- antsImageRead( getANTsRData( "r16" ) )
>      xfrm <- fitTransformToPairedPoints( moving, fixed, transformType = "Diffeo", domainImage = domainImage, numberOfFittingLevels = 6 )
Error in sigma > 0 : 
  comparison (6) is possible only for atomic and list types

> antsVersions
  Dependency                                   GitTag
1       ANTs 2251ca1dfca9d27a58623eb36784704a032d97a0
2  ANTsRCore caaa6c54b83f6166789296d8170a7caf9d8c8cc5
3    ANTsURL        https://github.com/ANTsX/ANTs.git
4        ITK e21a56d1227c5433066237060368cb4532b8a9d2

> sessionInfo()
R version 4.1.2 (2021-11-01)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 20.04.4 LTS

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.9.0
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.9.0

locale:
 [1] LC_CTYPE=en_CA.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_CA.UTF-8        LC_COLLATE=en_CA.UTF-8    
 [5] LC_MONETARY=en_CA.UTF-8    LC_MESSAGES=en_CA.UTF-8   
 [7] LC_PAPER=en_CA.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_CA.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] ANTsR_0.5.7.5   ANTsRCore_0.7.5

loaded via a namespace (and not attached):
[1] compiler_4.1.2      RcppEigen_0.3.3.9.1 magrittr_2.0.1     
[4] Matrix_1.4-0        tools_4.1.2         ITKR_0.5.3.3       
[7] Rcpp_1.0.7          grid_4.1.2          lattice_0.20-45    
ntustison commented 2 years ago

I've got a meeting in 10 minutes but I took a quick look and fixed an obvious bug. Try the example again. If the error persists, I'll look after my meeting.

SlicerMorph commented 2 years ago

Fixed, thank you!