[MINC-users] minctracc error after initialization from register-based xfm

[Audette, Michel]

Dear all, 

I am attempting a nonlinear registration between 2 volumes, where the extent of one volume is significantly less than that of the other (a high res MR of the inner ear of individual A, 0.078mm resolution, and a fairly high res CT of individual B, 0.2mm in x and y, 0.3mm in z), after reasonably successful rigid registration from 5 homologous pairs via the register program. I get the following error: 

maudette at icaw164201:~/work/earData> minctracc -transformation HumanS16885_etzoldCT.xfm HumanS16885_Stack_nuc_xyz.mnc etzold_emma_12792115_003002_mri_xyz.mnc HumanS16885_Stack_nuc_xyz_nonlin.mnc
WARNING: (rotmat_to_ang.c:187) step one: rz not in the range -pi/2..pi/2
Cannot convert R to radians!  -0.717591    0.045836   -0.694954
   0.408952   -0.779968   -0.473716
  -0.563756   -0.624138    0.540954
Error in minctracc in file minctracc.c, line 536
Could not initialize transformation parameters

Does this ring a bell with anyone? [etzold_emma_12792115_003002_mri_xyz.mnc is in fact a CT volume converted from Dicom, not mr data...]

Thanks for your kind support. 

Cheers, 

Michel


Michel Audette, Ph.D. 
Innovation Center Computer Assisted Surgery (ICCAS) 
Philipp-Rosenthal-Strasse 55
Leipzig, Germany
Phone: ++49 (0) 341 / 97 - 1 20 13
Fax: ++49 (0) 341 / 97 - 1 20 09


P.S: 
Other info: 
maudette at icaw164201:~/work/earData> minctracc -version
The program <minctracc> was built from:
Package mni_autoreg 0.98u, compiled by root at icaw164201 (x86_64-unknown-linux-gnu) on Thu Mar 16 13:13:40 CET 2006

Also, one volume has cosine info:  
       int zspace ;
                zspace:varid = "MINC standard variable" ;
                zspace:vartype = "dimension____" ;
                zspace:version = "MINC Version    1.0" ;
                zspace:comments = "Z increases from patient inferior to superior" ;
                zspace:spacing = "regular__" ;
                zspace:alignment = "centre" ;
                zspace:step = 0.299075200106658 ;
                zspace:start = -67.2154025304974 ;
                zspace:direction_cosines = 0., -0.0784590962903196, 0.99691733368886 ;
        int yspace ;
                yspace:varid = "MINC standard variable" ;
                yspace:vartype = "dimension____" ;
                yspace:version = "MINC Version    1.0" ;
                yspace:comments = "Y increases from patient posterior to anterior" ;
                yspace:spacing = "regular__" ;
                yspace:alignment = "centre" ;
                yspace:step = -0.1953125 ;
                yspace:start = 193.094103137385 ;
                yspace:direction_cosines = 0., 0.99691733368886, 0.0784590962903196 ;
        int xspace ;
                xspace:varid = "MINC standard variable" ;
                xspace:vartype = "dimension____" ;
                xspace:version = "MINC Version    1.0" ;
                xspace:comments = "X increases from patient left to right" ;
                xspace:spacing = "regular__" ;
                xspace:alignment = "centre" ;
                xspace:step = -0.1953125 ;
                xspace:start = 26.80468775 ;
                xspace:direction_cosines = 1., 0., 0. ;

The other not: 

        double xspace ;
                xspace:varid = "MINC standard variable" ;
                xspace:vartype = "dimension____" ;
                xspace:version = "MINC Version    1.0" ;
                xspace:comments = "X increases from patient left to right" ;
                xspace:spacing = "regular__" ;
                xspace:alignment = "centre" ;
                xspace:step = 0.078125 ;
                xspace:start = 0. ;
        double yspace ;
                yspace:varid = "MINC standard variable" ;
                yspace:vartype = "dimension____" ;
                yspace:version = "MINC Version    1.0" ;
                yspace:comments = "Y increases from patient posterior to anterior" ;
                yspace:spacing = "regular__" ;
                yspace:alignment = "centre" ;
                yspace:step = 0.078125 ;
                yspace:start = 0. ;
        double zspace ;
                zspace:varid = "MINC standard variable" ;
                zspace:vartype = "dimension____" ;
                zspace:version = "MINC Version    1.0" ;
                zspace:comments = "Z increases from patient inferior to superior" ;
                zspace:spacing = "regular__" ;
                zspace:alignment = "centre" ;
                zspace:step = 0.078125 ;
                zspace:start = 0. ;