Hi Douglas,

Please find below the mri_coreg terminal output :

$Id: mri_coreg.c,v 1.27 2016/04/30 15:11:49 greve Exp $

cwd /NAS/tupac/matthieu/FS5.3/207118_M0_2014-01-29/pet

cmdline mri_coreg --s 207118_M0_2014-01-29 --mov BS7_PET.lps.nii.gz --reg Pet2T1.BS7.register.dof6.mri_coreg.lta --regdat Pet2T1.BS7.register.dof6.mri_coreg.dat

sysname Linux

hostname yakuza

machine x86_64

user matthieu

dof 6

nsep 2

cras0 1

ftol 0.000000

linmintol 0.001000

bf 1

bflim 30.000000

bfnsamp 30

SmoothRef 0

SatPct 99.990000

MovOOB 0

optschema 1

Reading in mov BS7_PET.lps.nii.gz

Reading in ref /NAS/tupac/matthieu/FS5.3//207118_M0_2014-01-29/mri/brainmask.mgz

Reading in and applying refmask /NAS/tupac/matthieu/FS5.3//207118_M0_2014-01-29/mri/aparc+aseg.mgz

Setting cras translation parameters to align centers

Creating random numbers for coordinate dithering

Performing intensity dithering

Initial parameters 20.9141 11.8161 149.1538 0.0000 0.0000 0.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000

Separation list (2): 4 2 min = 2

DoSmoothing 1

DoCoordDither 1

DoIntensityDither 1

nitersmax 4

ftol 1.000e-07

linmintol 1.000e-03

SatPct 99.990000

Hist FWHM 7.000000 7.000000

nthreads 1

movsat = 10897.7666

mov gstd 0.8459 0.8459 0.8459

Smoothing mov

refsat = 119.0000

ref gstd 0.8459 0.8459 0.8459

Smoothing ref

COREGpreproc() done

Testing if mov and target overlap

Init cost -1.0011578057

nhits = 262144 out of 16777216, Percent Overlap: 100.0

Initial RefRAS-to-MovRAS

1.00000 0.00000 0.00000 20.91408;

0.00000 1.00000 0.00000 11.81612;

0.00000 0.00000 1.00000 149.15384;

0.00000 0.00000 0.00000 1.00000;

Initial RefVox-to-MovVox

1.00000 0.00000 0.00000 0.00000;

0.00000 0.00000 -1.00000 255.00000;

0.00000 -1.00000 0.00000 255.00000;

0.00000 0.00000 0.00000 1.00000;

sep = 4 -----------------------------------

COREGoptBruteForce() 30 1 30

Turning on MovOOB for BruteForce Search

#BF# sep= 4 iter=0 lim=30.0 delta=2.00 4.91408 21.81612 119.15384 2.00000 0.00000 0.00000 -1.0078773

Turning MovOOB back off after brute force search

---------------------------------

Init Powel Params dof = 6

Starting OpenPowel2(), sep = 4

InitialCost -1.0076701641

#@# 4 188 4.91408 21.81612 119.15384 2.00000 0.00000 0.00000 -1.0076702

fs_powell::minimize

nparams 6

maxfev 4

ftol 0.000000

linmin_xtol_ 0.001000

powell nthiter 0: fret = -1.007670

#@# 4 197 4.93866 21.81612 119.15384 2.00000 0.00000 0.00000 -1.0076704

#@# 4 198 4.96905 21.81612 119.15384 2.00000 0.00000 0.00000 -1.0076705

#@# 4 209 4.97536 20.19808 119.15384 2.00000 0.00000 0.00000 -1.0077189

#@# 4 210 4.97536 19.55424 119.15384 2.00000 0.00000 0.00000 -1.0077529

#@# 4 211 4.97536 18.51247 119.15384 2.00000 0.00000 0.00000 -1.0078081

#@# 4 212 4.97536 16.82685 119.15384 2.00000 0.00000 0.00000 -1.0078434

#@# 4 213 4.97536 16.77889 119.15384 2.00000 0.00000 0.00000 -1.0078460

#@# 4 214 4.97536 16.70129 119.15384 2.00000 0.00000 0.00000 -1.0078515

#@# 4 215 4.97536 16.57572 119.15384 2.00000 0.00000 0.00000 -1.0078617

#@# 4 216 4.97536 16.37254 119.15384 2.00000 0.00000 0.00000 -1.0078768

#@# 4 217 4.97536 14.89926 119.15384 2.00000 0.00000 0.00000 -1.0079292

#@# 4 218 4.97536 14.85724 119.15384 2.00000 0.00000 0.00000 -1.0079303

#@# 4 219 4.97536 12.95773 119.15384 2.00000 0.00000 0.00000 -1.0079467

#@# 4 220 4.97536 13.40911 119.15384 2.00000 0.00000 0.00000 -1.0079498

#@# 4 223 4.97536 13.62039 119.15384 2.00000 0.00000 0.00000 -1.0079516

#@# 4 229 4.97536 13.61634 117.53580 2.00000 0.00000 0.00000 -1.0086172

#@# 4 230 4.97536 13.61634 80.23408 2.00000 0.00000 0.00000 -1.0266197

#@# 4 235 4.97536 13.61634 69.09492 2.00000 0.00000 0.00000 -1.0285894

#@# 4 236 4.97536 13.61634 69.00322 2.00000 0.00000 0.00000 -1.0285913

#@# 4 237 4.97536 13.61634 68.30045 2.00000 0.00000 0.00000 -1.0285933

#@# 4 239 4.97536 13.61634 68.58790 2.00000 0.00000 0.00000 -1.0285972

#@# 4 241 4.97536 13.61634 68.56419 2.00000 0.00000 0.00000 -1.0285974

#@# 4 244 4.97536 13.61634 68.54494 2.00000 0.00000 0.00000 -1.0285981

#@# 4 250 4.97536 13.61634 68.54494 3.00000 0.00000 0.00000 -1.0286810

#@# 4 251 4.97536 13.61634 68.54494 4.61803 0.00000 0.00000 -1.0287892

#@# 4 255 4.97536 13.61634 68.54494 5.29042 0.00000 0.00000 -1.0287902

#@# 4 256 4.97536 13.61634 68.54494 5.02220 0.00000 0.00000 -1.0287969

#@# 4 266 4.97536 13.61634 68.54494 5.02421 1.00000 0.00000 -1.0288255

#@# 4 272 4.97536 13.61634 68.54494 5.02421 1.14711 0.00000 -1.0288261

#@# 4 275 4.97536 13.61634 68.54494 5.02421 1.09092 0.00000 -1.0288264

#@# 4 287 4.97536 13.61634 68.54494 5.02421 1.03666 1.00000 -1.0288312

#@# 4 291 4.97536 13.61634 68.54494 5.02421 1.03666 0.61803 -1.0288329

#@# 4 296 4.97536 13.61634 68.54494 5.02421 1.03666 0.72161 -1.0288330

#@# 4 297 4.97536 13.61634 68.54494 5.02421 1.03666 0.75100 -1.0288332

#@# 4 305 4.97536 13.61634 68.54494 5.02421 1.03666 0.74092 -1.0288333

#@# 4 315 4.97542 13.60814 68.49433 5.02724 1.03770 0.74166 -1.0288336

powell nthiter 1: fret = -1.028834

#@# 4 319 3.35739 13.60814 68.49433 5.02724 1.03770 0.74166 -1.0288464

#@# 4 320 4.02089 13.60814 68.49433 5.02724 1.03770 0.74166 -1.0288594

#@# 4 324 3.97462 13.60814 68.49433 5.02724 1.03770 0.74166 -1.0288596

#@# 4 327 4.00121 13.60814 68.49433 5.02724 1.03770 0.74166 -1.0288597

#@# 4 336 3.99013 14.60814 68.49433 5.02724 1.03770 0.74166 -1.0288848

#@# 4 341 3.99013 14.70860 68.49433 5.02724 1.03770 0.74166 -1.0288850

#@# 4 342 3.99013 14.68519 68.49433 5.02724 1.03770 0.74166 -1.0288851

#@# 4 356 3.99013 14.67684 68.49433 5.02724 1.03770 0.70149 -1.0288855

#@# 4 365 3.99013 14.67684 68.49433 6.02724 1.03770 0.68860 -1.0288972

#@# 4 384 3.99013 14.67684 68.49433 6.02724 1.06723 0.68860 -1.0288974

#@# 4 385 3.99013 14.67684 68.49433 6.02724 1.13709 0.68860 -1.0288981

#@# 4 386 3.99013 14.67684 68.49433 6.02724 1.24502 0.68860 -1.0288985

#@# 4 391 3.99013 14.67684 68.49433 6.02724 1.25019 0.68860 -1.0288986

#@# 4 402 3.99037 14.64484 68.29682 6.03904 1.25424 0.69150 -1.0289003

#@# 4 404 3.99030 14.65511 68.36022 6.03525 1.25294 0.69057 -1.0289006

powell nthiter 2: fret = -1.028901

#@# 4 413 3.80284 14.65511 68.36022 6.03525 1.25294 0.69057 -1.0289013

#@# 4 415 3.82815 14.65511 68.36022 6.03525 1.25294 0.69057 -1.0289015

#@# 4 422 3.82147 15.65511 68.36022 6.03525 1.25294 0.69057 -1.0289214

#@# 4 426 3.82147 15.27315 68.36022 6.03525 1.25294 0.69057 -1.0289227

#@# 4 427 3.82147 15.40752 68.36022 6.03525 1.25294 0.69057 -1.0289237

#@# 4 441 3.82147 15.39508 68.36022 6.03525 1.25294 0.44847 -1.0289254

#@# 4 443 3.82147 15.39508 68.36022 6.03525 1.25294 0.30487 -1.0289257

#@# 4 446 3.82147 15.39508 68.36022 6.03525 1.25294 0.32211 -1.0289260

#@# 4 457 3.82147 15.39508 68.36022 6.19153 1.25294 0.32011 -1.0289278

#@# 4 461 3.82147 15.39508 68.36022 6.18538 1.25294 0.32011 -1.0289279

#@# 4 469 3.82147 15.39508 68.36022 6.18638 0.63491 0.32011 -1.0289334

#@# 4 470 3.82147 15.39508 68.36022 6.18638 0.25294 0.32011 -1.0289339

#@# 4 475 3.82147 15.39508 68.36022 6.18638 0.21850 0.32011 -1.0289340

#@# 4 476 3.82147 15.39508 68.36022 6.18638 0.19721 0.32011 -1.0289341

#@# 4 480 3.82147 15.39508 68.36022 6.18638 0.20785 0.32011 -1.0289342

#@# 4 494 3.82131 15.41560 68.48689 6.17881 0.21339 0.31825 -1.0289359

powell nthiter 3: fret = -1.028936

#@# 4 503 3.20328 15.41560 68.48689 6.17881 0.21339 0.31825 -1.0289463

#@# 4 505 3.36275 15.41560 68.48689 6.17881 0.21339 0.31825 -1.0289472

#@# 4 508 3.28979 15.41560 68.48689 6.17881 0.21339 0.31825 -1.0289474

#@# 4 510 3.31021 15.41560 68.48689 6.17881 0.21339 0.31825 -1.0289475

#@# 4 518 3.31196 15.79757 68.48689 6.17881 0.21339 0.31825 -1.0289484

#@# 4 519 3.31196 15.62773 68.48689 6.17881 0.21339 0.31825 -1.0289508

#@# 4 527 3.31196 15.62417 68.48689 6.17881 0.21339 -0.29978 -1.0289558

#@# 4 529 3.31196 15.62417 68.48689 6.17881 0.21339 -0.13832 -1.0289564

#@# 4 534 3.31196 15.62417 68.48689 6.17881 0.21339 -0.13215 -1.0289565

#@# 4 561 3.31196 15.62417 68.48689 6.17881 0.14601 -0.13674 -1.0289570

#@# 4 565 3.31196 15.62417 68.48689 6.17881 0.13447 -0.13674 -1.0289571

Powell done niters total = 4

OptTimeSec 24.2 sec

OptTimeMin 0.40 min

nEvals 580

Final parameters 3.31203771 15.61334038 68.42005920 6.18280220 0.13176514 -0.13576585

Final cost -1.028957125573420

---------------------------------

sep = 2 -----------------------------------

---------------------------------

Init Powel Params dof = 6

Starting OpenPowel2(), sep = 2

InitialCost -1.0287065506

#@# 2 582 3.31204 15.61334 68.42006 6.18280 0.13177 -0.13577 -1.0287066

fs_powell::minimize

nparams 6

maxfev 4

ftol 0.000000

linmin_xtol_ 0.001000

powell nthiter 0: fret = -1.028707

#@# 2 589 3.10392 15.61334 68.42006 6.18280 0.13177 -0.13577 -1.0287081

#@# 2 604 3.09369 14.99531 68.42006 6.18280 0.13177 -0.13577 -1.0287125

#@# 2 606 3.09369 15.04368 68.42006 6.18280 0.13177 -0.13577 -1.0287126

#@# 2 621 3.09369 15.04550 68.23710 6.18280 0.13177 -0.13577 -1.0287137

#@# 2 622 3.09369 15.04550 68.21909 6.18280 0.13177 -0.13577 -1.0287138

#@# 2 630 3.09369 15.04550 68.21119 7.18280 0.13177 -0.13577 -1.0287204

#@# 2 635 3.09369 15.04550 68.21119 7.02697 0.13177 -0.13577 -1.0287205

#@# 2 636 3.09369 15.04550 68.21119 7.06769 0.13177 -0.13577 -1.0287206

#@# 2 650 3.09369 15.04550 68.21119 7.08699 0.03346 -0.13577 -1.0287210

#@# 2 657 3.09369 15.04550 68.21119 7.08699 0.02114 0.86423 -1.0287254

powell nthiter 1: fret = -1.028725

#@# 2 681 3.12473 15.04550 68.21119 7.08699 0.02114 0.95289 -1.0287255

#@# 2 695 3.12664 15.42746 68.21119 7.08699 0.02114 0.95289 -1.0287271

#@# 2 696 3.12664 15.34841 68.21119 7.08699 0.02114 0.95289 -1.0287273

#@# 2 720 3.12664 15.32519 68.20757 7.08297 0.02114 0.95289 -1.0287274

#@# 2 724 3.12664 15.32519 68.20757 7.07547 0.02114 0.95289 -1.0287275

#@# 2 745 3.12664 15.32519 68.20757 7.07447 0.01759 0.81959 -1.0287278

#@# 2 762 3.14700 15.49804 68.20535 7.06672 0.01540 0.72493 -1.0287280

powell nthiter 2: fret = -1.028728

#@# 2 792 3.14765 15.49494 68.20538 7.06686 0.01544 0.59641 -1.0287281

#@# 2 803 3.14765 15.49494 68.11650 7.06686 0.01544 0.59810 -1.0287282

#@# 2 804 3.14765 15.49494 68.14099 7.06686 0.01544 0.59810 -1.0287284

#@# 2 805 3.14765 15.49494 68.15399 7.06686 0.01544 0.59810 -1.0287285

#@# 2 807 3.14765 15.49494 68.16380 7.06686 0.01544 0.59810 -1.0287286

powell nthiter 3: fret = -1.028729

#@# 2 867 3.12259 15.49539 68.16066 7.07157 0.01111 0.59788 -1.0287287

Powell done niters total = 4

OptTimeSec 161.2 sec

OptTimeMin 2.69 min

nEvals 938

Final parameters 3.12265301 15.49595165 68.16065216 7.06536150 0.01209922 0.59759676

Final cost -1.028728664232616

---------------------------------

mri_coreg utimesec 241.181825

mri_coreg stimesec 0.360032

mri_coreg ru_maxrss 504732

mri_coreg ru_ixrss 0

mri_coreg ru_idrss 0

mri_coreg ru_isrss 0

mri_coreg ru_minflt 381121

mri_coreg ru_majflt 0

mri_coreg ru_nswap 0

mri_coreg ru_inblock 71512

mri_coreg ru_oublock 16

mri_coreg ru_msgsnd 0

mri_coreg ru_msgrcv 0

mri_coreg ru_nsignals 0

mri_coreg ru_nvcsw 505

mri_coreg ru_nivcsw 28383

Final RefRAS-to-MovRAS

0.99995 0.01043 0.00021 3.12265;

-0.01038 0.99235 0.12300 15.49595;

0.00107 -0.12300 0.99241 68.16065;

0.00000 0.00000 0.00000 1.00000;

Final RefVox-to-MovVox

0.99995 0.00021 -0.01043 19.00880;

-0.01038 0.12300 -0.99235 243.54027;

-0.00107 -0.99241 -0.12300 188.07765;

0.00000 0.00000 0.00000 1.00000;

Final parameters 3.1227 15.4960 68.1607 7.0654 0.0121 0.5976

nhits = 1287172 out of 16777216, Percent Overlap: 61.4

mri_coreg RunTimeSec 241.7 sec

To check run:

tkregisterfv --mov BS7_PET.lps.nii.gz --targ /NAS/tupac/matthieu/FS5.3//207118_M0_2014-01-29/mri/brainmask.mgz --reg Pet2T1.BS7.register.dof6.mri_coreg.lta --s 207118_M0_2014-01-29 --surfs

mri_coreg done

Best regards,

Matthieu

2016-09-30 21:07 GMT+02:00 Douglas N Greve <greve@nmr.mgh.harvard.edu>:

Can you send the mri_coreg terminal output?

On 09/29/2016 11:19 AM, Matthieu Vanhoutte wrote:
> Dear Freesurfer's experts,
>
> I would like to register PET image onto T1 anatomical. I used
> bbregister with the following command and it failed :
> /bbregister --s SUBJID --init-fsl --t2 --mov PET.nii.gz --reg
> Pet2T1.register.dof6.dat --lta Pet2T1.register.dof6.lta --init-reg-out
> Pet2T1.init.register.dof6.dat --o rPET.nii.gz /
> /
> /
> Having heard about mri_coreg solving some bbregister fails, I tested
> it with the following command :
> /mri_coreg --s SUBJID --mov PET.nii.gz --reg
> Pet2T1.register.dof6.mri_coreg.lta --regdat
> Pet2T1.register.dof6.mri_coreg.dat/
> /
> /
> But it failed in the same way than bbregister. Do you have any advice
> on this problem ?
>
> Best regards,
> Matthieu
>
>
> _______________________________________________
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

--
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422

Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.