
New invocation of recon-all 



Wed Jun 20 09:44:29 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri
/usr/local/freesurfer/bin//recon-all
-s sub089mri -all -wsatlas
subjid sub089mri
setenv SUBJECTS_DIR /home/richfrye/Desktop/MRI/sub089
FREESURFER_HOME /usr/local/freesurfer
Actual FREESURFER_HOME /usr/local/freesurfer
build-stamp.txt: freesurfer-Linux-centos4_x86_64-stable-pub-v3.0.5
Linux localhost.localdomain 2.6.9-55.ELsmp #1 SMP Fri Apr 20 16:36:54 EDT 2007 x86_64 x86_64 x86_64 GNU/Linux
cputime      unlimited
filesize     unlimited
datasize     unlimited
stacksize    10240 kbytes
coredumpsize 0 kbytes
memoryuse    unlimited
vmemoryuse   unlimited
descriptors  1024 
memorylocked 32 kbytes
maxproc      16372 
########################################
program versions used
$Id: recon-all,v 1.17.2.16 2006/12/12 19:44:21 nicks Exp $
$Id: mri_motion_correct.fsl,v 1.4 2006/01/01 21:56:44 greve Exp $
mri_convert -all-info 
ProgramName: mri_convert ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_convert.c,v 1.121.2.1 2006/08/28 18:00:04 greve Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
$Id: talairach,v 1.2 2005/12/12 23:04:35 greve Exp $
mri_convert --version 
stable3
# $Id: mritotal.in,v 1.12 2004/02/12 05:55:18 rotor Exp $
Program nu_correct, built from:
Package MNI N3, version 1.10, compiled by nicks@minerva (x86_64-unknown-linux-gnu) on 2005-11-15 at 20:37:35
ProgramName: mri_normalize ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_normalize.c,v 1.44.2.2 2006/06/28 23:38:13 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_watershed ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_watershed.cpp,v 1.48 2006/02/27 16:57:57 segonne Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_segment ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_segment.c,v 1.28.2.1 2006/04/12 02:02:58 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_fill ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_fill.c,v 1.101 2006/02/22 23:05:06 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_tessellate ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_tessellate.c,v 1.27.2.1 2006/05/02 21:41:01 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_smooth ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_smooth.c,v 1.13.2.1 2006/05/02 21:41:01 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_inflate ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_inflate.c,v 1.30.2.1 2006/05/02 21:41:01 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_sphere ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_sphere.c,v 1.38 2006/02/08 00:17:40 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_fix_topology ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_fix_topology.c,v 1.41.2.1 2007/01/05 16:58:33 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_euler_number ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_euler_number.c,v 1.4 2003/09/05 04:45:41 kteich Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_make_surfaces ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_make_surfaces.c,v 1.70.2.6 2007/01/09 19:15:32 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mris_register ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mris_register.c,v 1.32 2005/12/07 14:22:57 fischl Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mrisp_paint ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mrisp_paint.c,v 1.6 2005/02/14 04:30:19 segonne Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_label2label ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_label2label.c,v 1.23.2.1 2006/09/20 20:42:55 greve Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_em_register ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_em_register.c,v 1.52 2006/01/11 21:01:40 greve Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_ca_normalize ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_ca_normalize.c,v 1.31 2005/08/15 14:14:22 fischl Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_ca_register ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_ca_register.c,v 1.48.2.7 2006/10/24 01:48:52 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
ProgramName: mri_ca_label ProgramArguments: -all-info ProgramVersion: $Name: stable3 $ TimeStamp: 07/06/20-13:44:29-GMT CVS: $Id: mri_ca_label.c,v 1.68.2.2 2006/08/16 00:18:42 nicks Exp $ User: richfrye Machine: localhost.localdomain Platform: Linux PlatformVersion: 2.6.9-55.ELsmp CompilerName: GCC CompilerVersion: 30400
#######################################
GCA RB_all_2006-02-15.gca
GCASkull RB_all_withskull_2006-02-15.gca
AvgCurvTif average.curvature.filled.buckner40.tif
GCS curvature.buckner40.filled.desikan_killiany.gcs
#######################################
#--------------------------------------------
#@# MotionCor Wed Jun 20 09:44:29 EDT 2007
Found 2 runs
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/001.mgz
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/002.mgz
/home/richfrye/Desktop/MRI/sub089/sub089mri

 mri_motion_correct.fsl -o /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz -wild /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/001.mgz /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/002.mgz 

tmpdir is /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449
--------------------------------------------------------------
mri_motion_correct.fsl logfile is /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz.mri_motion_correct.fsl.log
--------------------------------------------------------------
-----------------------------------------
Converting /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/001.mgz 
Wed Jun 20 09:44:45 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri
mri_convert /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/001.mgz /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-1.img -odt float 
reading from /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/001.mgz...
TR=8.55, TE=4.03, TI=0.00, flip angle=8.00
i_ras = (-0, -1, 0)
j_ras = (-0, -0, -1)
k_ras = (1, 0, 0)
changing data type from 4 to 3 (noscale = 0)...
writing to /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-1.img...
Analyze Output Matrix
-0.000  -0.000   0.940  -89.096;
-0.938  -0.000   0.000   138.803;
 0.000  -0.938   0.000   85.533;
 0.000   0.000   0.000   1.000;
--------------------
INFO: set hdr.hist.orient to -1
-----------------------------------------
Converting /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/002.mgz 
Wed Jun 20 09:44:49 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri
mri_convert /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/002.mgz /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.img -odt float 
reading from /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig/002.mgz...
TR=8.55, TE=4.03, TI=0.00, flip angle=8.00
i_ras = (-0, -1, 0)
j_ras = (-0, -0, -1)
k_ras = (1, 0, 0)
changing data type from 4 to 3 (noscale = 0)...
writing to /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.img...
Analyze Output Matrix
-0.000  -0.000   0.940  -89.096;
-0.938  -0.000   0.000   138.803;
 0.000  -0.938   0.000   85.533;
 0.000   0.000   0.000   1.000;
--------------------
INFO: set hdr.hist.orient to -1
-----------------------------------------
Motion Correcting /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.img
Wed Jun 20 09:44:54 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri
fsl_rigid_register -maxangle 90 -r /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-1.img -i /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.img -o /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/resampled.img -tmpdir /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449 -cleanup
\n\n
$Id: fsl_rigid_register,v 1.14 2006/02/13 20:13:19 greve Exp $
Wed Jun 20 09:44:54 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/resampled.img
Reference is analyze, not converting
Input is analyze, not converting
Output is analyze, not converting
--------------------------------------
/home/richfrye/Desktop/MRI/sub089/sub089mri
calling eval flirt -in /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.hdr -out /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/resampled.hdr -bins 256 -cost corratio -searchrx -90 90 -searchry -90 90 -searchrz -90 90 -interp trilinear -dof 6 -ref /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-1.hdr -omat /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/resampled.fslmat
WARNING: Both reference and input images have an sform matrix set
  The output image will use the sform from the reference image
 
Started at Wed Jun 20 09:44:54 EDT 2007 
Ended   at Wed Jun 20 09:49:17 EDT 2007
 
fsl_rigid_register Done
 
To check results, run:
tkmedit -f /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-1.img -aux /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/resampled.img
 
mv /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/resampled.img /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.img
-----------------------------------------
Averaging 
Wed Jun 20 09:49:17 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri
inhibiting isotropic volume interpolation
1 of 2: reading /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-1.img...
2 of 2: reading /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/cor-2.img...
writing to /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/avgvol.img...
alignment and averaging took 0 minutes and 18 seconds.
Analyze Output Matrix
 0.000   0.000   0.940  -89.096;
-0.938   0.000   0.000   138.803;
 0.000  -0.938   0.000   85.533;
 0.000   0.000   0.000   1.000;
--------------------
INFO: set hdr.hist.orient to -1
-----------------------------------------
Converting Average to output 
Wed Jun 20 09:49:35 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri
mri_convert /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/avgvol.img /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz -tr 8.5512 -te 4.0256 -TI 0 -flip_angle 0.139626
mri_convert /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/avgvol.img /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz -tr 8.5512 -te 4.0256 -TI 0 -flip_angle 0.139626 
reading from /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/tmp-mri_motion_correct.fsl-20449/avgvol.img...
TR=8.55, TE=0.00, TI=0.00, flip angle=0.00
i_ras = (0, -1, 0)
j_ras = (0, 0, -1)
k_ras = (1, 0, 0)
writing to /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz...
-----------------------------------------
Started at: Wed Jun 20 09:44:33 EDT 2007
Ended at:   Wed Jun 20 09:49:46 EDT 2007
mri_motion_correct.fsl: done
/home/richfrye/Desktop/MRI/sub089/sub089mri

 mri_convert /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig.mgz --conform 

mri_convert /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig.mgz --conform 
reading from /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/rawavg.mgz...
TR=8.55, TE=4.03, TI=0.00, flip angle=8.00
i_ras = (0, -1, 0)
j_ras = (0, 0, -1)
k_ras = (1, 0, 0)
Original Data has (0.9375, 0.9375, 0.94) mm size and (256, 256, 181) voxels.
Data is conformed to 1 mm size and 256 voxels for all directions
changing data type from 3 to 0 (noscale = 0)...
MRIchangeType: Building histogram 
Reslicing using trilinear interpolation 
writing to /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig.mgz...

 mri_add_xform_to_header -c /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/transforms/talairach.xfm /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig.mgz /home/richfrye/Desktop/MRI/sub089/sub089mri/mri/orig.mgz 

INFO: extension is mgz
#--------------------------------------------
#@# Nu Intensity Correction Wed Jun 20 09:50:09 EDT 2007

 mri_nu_correct.mni --i orig.mgz --o nu.mgz --n 2 

/home/richfrye/Desktop/MRI/sub089/sub089mri/mri
/usr/local/freesurfer/bin//mri_nu_correct.mni
--i orig.mgz --o nu.mgz --n 2
nIters 2
$Id: mri_nu_correct.mni,v 1.3.2.1 2006/11/24 20:16:38 nicks Exp $
Linux localhost.localdomain 2.6.9-55.ELsmp #1 SMP Fri Apr 20 16:36:54 EDT 2007 x86_64 x86_64 x86_64 GNU/Linux
Wed Jun 20 09:50:09 EDT 2007
Program nu_correct, built from:
Package MNI N3, version 1.10, compiled by nicks@minerva (x86_64-unknown-linux-gnu) on 2005-11-15 at 20:37:35
tmpdir is ./tmp.mri_nu_correct.mni.20645
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri
mri_convert orig.mgz ./tmp.mri_nu_correct.mni.20645/nu0.mnc
mri_convert orig.mgz ./tmp.mri_nu_correct.mni.20645/nu0.mnc 
reading from orig.mgz...
TR=8.55, TE=4.03, TI=0.00, flip angle=8.00
i_ras = (-1, 0, 0)
j_ras = (0, 0, -1)
k_ras = (0, 1, 0)
writing to ./tmp.mri_nu_correct.mni.20645/nu0.mnc...
 
--------------------------------------------------------
Iteration 1 Wed Jun 20 09:50:18 EDT 2007
nu_correct -clobber ./tmp.mri_nu_correct.mni.20645/nu0.mnc ./tmp.mri_nu_correct.mni.20645/nu1.mnc
[richfrye@localhost.localdomain:/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/] [2007-06-20 09:50:18] running:
  /usr/local/freesurfer/mni/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir /usr/tmp/nu_correct_20711/ ./tmp.mri_nu_correct.mni.20645/nu0.mnc ./tmp.mri_nu_correct.mni.20645/nu1.imp

Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Number of iterations: 50 
CV of field change: 0.00148634
[richfrye@localhost.localdomain:/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/] [2007-06-20 09:52:08] running:
  /usr/local/freesurfer/mni/bin/make_template -quiet -shrink 3 ./tmp.mri_nu_correct.mni.20645/nu0.mnc /usr/tmp/nu_correct_20711//template.mnc

Transforming slices:......................................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Transforming slices:................................................................................................................................................................................................................................................................Done
 
 
--------------------------------------------------------
Iteration 2 Wed Jun 20 09:52:15 EDT 2007
nu_correct -clobber ./tmp.mri_nu_correct.mni.20645/nu1.mnc ./tmp.mri_nu_correct.mni.20645/nu2.mnc
[richfrye@localhost.localdomain:/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/] [2007-06-20 09:52:15] running:
  /usr/local/freesurfer/mni/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir /usr/tmp/nu_correct_21358/ ./tmp.mri_nu_correct.mni.20645/nu1.mnc ./tmp.mri_nu_correct.mni.20645/nu2.imp

Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Processing:.................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Number of iterations: 28 
CV of field change: 0.000998715
[richfrye@localhost.localdomain:/home/richfrye/Desktop/MRI/sub089/sub089mri/mri/] [2007-06-20 09:53:19] running:
  /usr/local/freesurfer/mni/bin/make_template -quiet -shrink 3 ./tmp.mri_nu_correct.mni.20645/nu1.mnc /usr/tmp/nu_correct_21358//template.mnc

Transforming slices:......................................................................................Done
Not implemented yet in cache_volume_range_has_changed()
Not implemented yet in cache_volume_range_has_changed()
Transforming slices:................................................................................................................................................................................................................................................................Done
 
 
 
mri_convert ./tmp.mri_nu_correct.mni.20645/nu2.mnc nu.mgz --like orig.mgz
mri_convert ./tmp.mri_nu_correct.mni.20645/nu2.mnc nu.mgz --like orig.mgz 
reading from ./tmp.mri_nu_correct.mni.20645/nu2.mnc...
TR=0.00, TE=0.00, TI=0.00, flip angle=0.00
i_ras = (-1, 0, 0)
j_ras = (0, 0, -1)
k_ras = (0, 1, 0)
INFO: transform src into the like-volume: orig.mgz
writing to nu.mgz...
 
 
Wed Jun 20 09:53:40 EDT 2007
mri_nu_correct.mni done
#--------------------------------------------
#@# Talairach Wed Jun 20 09:53:40 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 talairach --i nu.mgz --xfm transforms/talairach.auto.xfm 

/home/richfrye/Desktop/MRI/sub089/sub089mri/mri
/usr/local/freesurfer/bin//talairach
--i nu.mgz --xfm transforms/talairach.auto.xfm
$Id: talairach,v 1.2 2005/12/12 23:04:35 greve Exp $
Linux localhost.localdomain 2.6.9-55.ELsmp #1 SMP Fri Apr 20 16:36:54 EDT 2007 x86_64 x86_64 x86_64 GNU/Linux
Wed Jun 20 09:53:41 EDT 2007
tmpdir is transforms/tmp.talairach.21766
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri
mri_convert nu.mgz transforms/tmp.talairach.21766/src.mnc
mri_convert nu.mgz transforms/tmp.talairach.21766/src.mnc 
reading from nu.mgz...
TR=8.55, TE=4.03, TI=0.00, flip angle=8.00
i_ras = (-1, 0, 0)
j_ras = (0, 0, -1)
k_ras = (0, 1, 0)
writing to transforms/tmp.talairach.21766/src.mnc...
--------------------------------------------
mritotal -verbose -debug -clobber -modeldir /usr/local/freesurfer/mni/bin/../share/mni_autoreg -protocol icbm transforms/tmp.talairach.21766/src.mnc transforms/talairach.auto.xfm
 
 
Wed Jun 20 09:54:34 EDT 2007
talairach done
#--------------------------------------------
#@# Intensity Normalization Wed Jun 20 09:54:34 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_normalize -g 1 nu.mgz T1.mgz 

using max gradient = 1.000
reading from nu.mgz...
normalizing image...
talairach transform
 0.855   0.023  -0.021  -0.707;
-0.028   0.778  -0.305  -59.446;
 0.013   0.371   0.946   90.036;
 0.000   0.000   0.000   1.000;
INFO: Modifying talairach volume c_(r,a,s) based on average_305
building Voronoi diagram...
performing soap bubble smoothing...
3d normalization pass 1 of 2
building Voronoi diagram...
performing soap bubble smoothing...
3d normalization pass 2 of 2
building Voronoi diagram...
performing soap bubble smoothing...
writing output to T1.mgz
3D bias adjustment took 7 minutes and 17 seconds.
white matter peak found at 110
gm peak at 62 (62), valley at 25 (25)
csf peak at 31, setting threshold to 51
white matter peak found at 110
gm peak at 56 (56), valley at 24 (24)
csf peak at 28, setting threshold to 46
#--------------------------------------------
#@# Skull Stripping Wed Jun 20 10:01:54 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_watershed -keep brainmask.auto.mgz brainmask.mgz brainmask.mgz -atlas T1.mgz brainmask.auto.mgz 


Mode:          Atlas analysis

*********************************************************
The input file is T1.mgz
The output file is brainmask.auto.mgz
If this is incorrect, please exit with CTL-C


*************************WATERSHED**************************
preflooding height equal to 25 percent
Sorting...
      T1-weighted MRI image
      modification of the preflooding height to 15 percent
      Count how many 110 voxels are present :  216849

      Find the largest 110-component...
             heap usage = 219432 Kbytes.
             removing small segments (less than 1 percent of maxarea).done
      And identify it as the main brain basin...done
      Main component: 105660 voxels
      first estimation of the COG coord: x=122 y=159 z=125 r=88
      first estimation of the main basin volume: 2890158 voxels
      global maximum in x=153, y=120, z=168, Imax=255
      CSF=14, WM_intensity=110, WM_VARIANCE=5
      WM_MIN=110, WM_HALF_MIN=110, WM_HALF_MAX=110, WM_MAX=110 
done
Analyze

      main basin size= 1156184 voxels, voxel volume =1.000 
                     = 1156184 mmm3 = 1156.184 cm3
done
PostAnalyze...
      ambiguous basin, merged: at least 24 ambiguous voxels; size: 12289 voxels
      ambiguous basin, merged: at least 41 ambiguous voxels; size: 66004 voxels
      ambiguous basin, merged: at least 28 ambiguous voxels; size: 20303 voxels
      ambiguous basin, merged: at least 10 ambiguous voxels; size: 844 voxels
      ambiguous basin, merged: at least 13 ambiguous voxels; size: 2133 voxels
      ambiguous basin, merged: at least 5 ambiguous voxels; size: 96 voxels
      ambiguous basin, merged: at least 11 ambiguous voxels; size: 1058 voxels
      ambiguous basin, merged: at least 22 ambiguous voxels; size: 10177 voxels
      ambiguous basin, merged: at least 28 ambiguous voxels; size: 20164 voxels
      ***** 9 basin(s) merged in 5 iteration(s)
      ***** 133068 voxel(s) added to the main basin
      ambiguous basin, non merged: 34 ambiguous voxels; size: 53548 voxels
      ambiguous basin, non merged: 8 ambiguous voxels; size: 6559 voxels
done

****************TEMPLATE DEFORMATION****************

      second estimation of the COG coord: x=120,y=166, z=186, r=108108 iterations
      before analyzing : CSF_MAX=3, TRANSITION=2, GM_MIN=0, GM=102
      after  analyzing : CSF_MAX=2, TRANSITION=18, GM_MIN=27, GM=39
      mri_strip_skull: done peeling brain
      highly tesselated surface with 10242 vertices
      matching...81 iterations

*********************VALIDATION*********************

Rigid alignment...
      scanning 32.00 degree nbhd, min sse = 27.11 at (32.00, -8.00, -32.00)
      scanning 16.00 degree nbhd, min sse = 20.11 at (16.00, -12.00, -16.00)
      scanning  8.00 degree nbhd, min sse = 16.96 at ( 8.00, -6.00, -8.00)
      scanning  4.00 degree nbhd, min sse = 15.76 at ( 4.00,  0.00, -4.00)
done
      before rotation: sse = 40.26, sigma = 96.41
      after  rotation: sse = 15.76, sigma = 30.57
Localization of inacurate regions: Erosion-Dilation steps
      the sse mean is 24.70, its var is 38.40   
      before Erosion-Dilatation 42.17% of inacurate vertices
      after  Erosion-Dilatation 67.14% of inacurate vertices
            48.37% of 'positive' inacurate vertices
            51.63% of 'negative' inacurate vertices

      THE SEGMENTATION IS PROBABLY NOT CORRECT

Scaling of atlas fields onto current surface fields
Correcting the shape of the surface...50 iterations
Scaling of atlas fields onto current surface fields

********FINAL ITERATIVE TEMPLATE DEFORMATION********
Compute Local values csf/gray
Fine Segmentation...
      mri_strip_skull: done peeling brainComputing COG for corrected surface
Initial Configuration: NRG=151845.467177, ( 4.817486 , 85.185130 , -81.401351 )
Final Minimum found: NRG=143772.253688, ( 2.755039 , 83.172085 , -82.976267 )done

      The shape of the surface was incorrect,
      hence we rigidly realign the surface with the template
Rigid alignment...
      scanning 32.00 degree nbhd, min sse = 19.44 at (32.00,  8.00, -32.00)
      scanning 16.00 degree nbhd, min sse = 12.14 at (16.00, -16.00, -16.00)
      scanning  8.00 degree nbhd, min sse =  9.34 at ( 8.00, -8.00, -2.00)
      scanning  4.00 degree nbhd, min sse =  8.24 at ( 4.00, -4.00, -1.00)
done
      before rotation: sse = 36.34, sigma = 107.02
      after  rotation: sse = 8.24, sigma = 14.27
Localization of inacurate regions: Erosion-Dilation steps
      the sse mean is 15.15, its var is 27.06   
      before Erosion-Dilatation 25.46% of inacurate vertices
      after  Erosion-Dilatation 39.80% of inacurate vertices
            58.51% of 'positive' inacurate vertices
            41.49% of 'negative' inacurate vertices

      THE SEGMENTATION IS PROBABLY NOT CORRECT
102 iterations

      mri_strip_skull: done peeling brain

Brain Size = 6329924 voxels, voxel volume = 1.000 mm3
           = 6329924 mmm3 = 6329.924 cm3


******************************
Save...done
Keeping brain edits brainmask.auto.mgz brainmask.mgz
Keeping edits ...
Saving kept edits to brainmask.mgz .....
#-------------------------------------
#@# EM Registration Wed Jun 20 10:11:18 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_em_register -mask brainmask.mgz nu.mgz /usr/local/freesurfer/average/RB_all_2006-02-15.gca transforms/talairach.lta 

using MR volume brainmask.mgz to mask input volume...
reading 1 input volumes...
logging results to talairach.log
reading '/usr/local/freesurfer/average/RB_all_2006-02-15.gca'...
setting gca type = Normal gca type
average std = 7.4   using min determinant for regularization = 5.5
0 singular and 1662 ill-conditioned covariance matrices regularized
reading 'nu.mgz'...
total sample mean = 87.9 (483 zeros)
register_mri: find_optimal_transform
find_optimal_transform: nsamples 2242, passno 0, spacing 8
************************************************
First Search limited to translation only.
************************************************
Found translation: (-18.8, -50.6, -56.2): log p = -89985.1
Result so far: scale 1.000: max_log_p=-68168.0, old_max_log_p =-89985.1 (thresh=-89895.2)
 1.104   0.404   0.214  -152.650;
-0.317   1.079  -0.404   23.781;
-0.291   0.279   0.975  -57.119;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-62684.6, old_max_log_p =-68168.0 (thresh=-68099.8)
 1.161   0.510   0.076  -167.494;
-0.527   1.142  -0.651   75.095;
-0.221   0.509   0.829  -81.178;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-62289.0, old_max_log_p =-62684.6 (thresh=-62621.9)
 1.161   0.510   0.076  -167.494;
-0.587   1.277  -0.563   37.792;
-0.155   0.366   0.968  -84.472;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-57470.8, old_max_log_p =-62289.0 (thresh=-62226.7)
 1.276   0.422   0.278  -174.186;
-0.473   1.524  -0.545  -48.902;
-0.356   0.266   1.081  -62.067;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-54228.7, old_max_log_p =-57470.8 (thresh=-57413.4)
 1.661  -0.174  -0.023  -49.622;
-0.107   1.861  -0.636  -153.954;
 0.243   0.375   1.193  -198.210;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-51018.7, old_max_log_p =-54228.7 (thresh=-54174.5)
 1.804   0.744   0.019  -255.835;
-0.887   2.152  -0.542  -151.473;
-0.036   0.069   1.523  -163.836;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-51018.7, old_max_log_p =-51018.7 (thresh=-50967.6)
 1.804   0.744   0.019  -255.835;
-0.887   2.152  -0.542  -151.473;
-0.036   0.069   1.523  -163.836;
 0.000   0.000   0.000   1.000;
dt = 5.00e-06, momentum=0.80, tol=1.00e-05
l_intensity = 1.0000
Aligning input volume to GCA...
Transform matrix
 1.80413   0.74438   0.01863  -255.83484;
-0.88666   2.15236  -0.54215  -151.47304;
-0.03615   0.06936   1.52315  -163.83559;
 0.00000   0.00000   0.00000   1.00000;
nsamples 2242
Quasinewton: input matrix
 1.80413   0.74438   0.01863  -255.83484;
-0.88666   2.15236  -0.54215  -151.47304;
-0.03615   0.06936   1.52315  -163.83559;
 0.00000   0.00000   0.00000   1.00000;
pass 2 through quasi-newton minimization...
pass 3 through quasi-newton minimization...
pass 4 through quasi-newton minimization...
pass 5 through quasi-newton minimization...
fp = 50514.679688, fret = 50514.679688, ftol = 0.000010, EPS = 0.000000
dt = 5.00e-06, momentum=0.80, tol=1.00e-07
l_intensity = 1.0000
Aligning input volume to GCA...
Transform matrix
 1.80413   0.74438   0.01863  -255.83484;
-0.88666   2.15236  -0.54215  -151.47304;
-0.03615   0.06936   1.52315  -163.83559;
 0.00000   0.00000   0.00000   1.00000;
nsamples 242283
Quasinewton: input matrix
 1.80413   0.74438   0.01863  -255.83484;
-0.88666   2.15236  -0.54215  -151.47304;
-0.03615   0.06936   1.52315  -163.83559;
 0.00000   0.00000   0.00000   1.00000;
pass 2 through quasi-newton minimization...
pass 3 through quasi-newton minimization...
pass 4 through quasi-newton minimization...
freeing gibbs priors...done.
bounding unknown intensity as < 16.0 or > 799.3 
************************************************
spacing=8, using 2242 sample points, tol=1.00e-05...
************************************************
resetting wm mean[0]: 109 --> 116
resetting gm mean[0]: 68 --> 68
input volume #1 is the most T1-like
using real data threshold=17.0
using (99, 132, 166) as brain centroid...
mean wm in atlas = 116, using box (80,109,150) --> (117, 154,182) to find MRI wm
before smoothing, mri peak at 116
after smoothing, mri peak at 76, scaling input intensities by 1.526
scaling channel 0 by 1.52632
initial log_p = -173951.6

****************************************
Nine parameter search.  iteration 0 nscales = 0 ...
****************************************
****************************************
Nine parameter search.  iteration 1 nscales = 1 ...
****************************************
****************************************
Nine parameter search.  iteration 2 nscales = 2 ...
****************************************
****************************************
Nine parameter search.  iteration 3 nscales = 3 ...
****************************************
****************************************
Nine parameter search.  iteration 4 nscales = 4 ...
****************************************
****************************************
Nine parameter search.  iteration 5 nscales = 5 ...
****************************************
****************************************
Nine parameter search.  iteration 6 nscales = 6 ...
****************************************
reducing scale to 0.2500
***********************************************
Computing MAP estimate using 2242 samples...
***********************************************
dfp_em_step_func: 007: -log(p) = 51018.7
after pass:transform: ( 1.80, 0.74, 0.02, -255.83)
                      ( -0.89, 2.15, -0.54, -151.47)
                      ( -0.04, 0.07, 1.52, -163.84)
dfp_em_step_func: 008: -log(p) = 50830.7
after pass:transform: ( 1.80, 0.74, 0.02, -255.83)
                      ( -0.89, 2.15, -0.54, -151.47)
                      ( -0.04, 0.06, 1.52, -163.84)
dfp_em_step_func: 009: -log(p) = 50738.9
after pass:transform: ( 1.80, 0.74, 0.02, -255.83)
                      ( -0.89, 2.15, -0.54, -151.47)
                      ( -0.04, 0.06, 1.52, -163.84)
dfp_em_step_func: 010: -log(p) = 50723.4
after pass:transform: ( 1.80, 0.74, 0.02, -255.83)
                      ( -0.89, 2.15, -0.54, -151.47)
                      ( -0.04, 0.06, 1.52, -163.84)
dfp_em_step_func: 011: -log(p) = 50514.7
after pass:transform: ( 1.80, 0.74, 0.02, -255.83)
                      ( -0.89, 2.15, -0.54, -151.47)
                      ( -0.04, 0.06, 1.52, -163.84)
outof QuasiNewtonEMA: 013: -log(p) = 50514.7  tol 0.000010
Resulting transform:
 1.804   0.744   0.019  -255.835;
-0.887   2.152  -0.542  -151.473;
-0.036   0.069   1.523  -163.836;
 0.000   0.000   0.000   1.000;

pass 1, spacing 8: log(p) = -51018.7 (old=-173951.6)
transform before final EM align:
 1.804   0.744   0.019  -255.835;
-0.887   2.152  -0.542  -151.473;
-0.036   0.069   1.523  -163.836;
 0.000   0.000   0.000   1.000;

**************************************************
 EM alignment process ...
 Computing final MAP estimate using 242283 samples. 
**************************************************
dfp_em_step_func: 012: -log(p) = 4515679.5
after pass:transform: ( 1.80, 0.74, 0.02, -255.83)
                      ( -0.89, 2.15, -0.54, -151.47)
                      ( -0.04, 0.07, 1.52, -163.84)
dfp_em_step_func: 013: -log(p) = 4470533.5
after pass:transform: ( 1.80, 0.74, 0.01, -255.83)
                      ( -0.88, 2.17, -0.53, -151.47)
                      ( -0.03, 0.07, 1.53, -163.84)
dfp_em_step_func: 014: -log(p) = 4460801.0
after pass:transform: ( 1.81, 0.75, 0.02, -255.83)
                      ( -0.88, 2.17, -0.53, -151.47)
                      ( -0.04, 0.07, 1.53, -163.84)
dfp_em_step_func: 015: -log(p) = 4457007.0
after pass:transform: ( 1.81, 0.75, 0.02, -255.83)
                      ( -0.88, 2.16, -0.53, -151.47)
   pass 5 through quasi-newton minimization...
pass 6 through quasi-newton minimization...
                   ( -0.04, 0.07, 1.52, -163.84)
dfp_em_step_func: 016: -log(p) = 4456456.5
after pass:transform: ( 1.81, 0.75, 0.02, -255.83)
                      ( -0.88, 2.17, -0.53, -151.47)
                      ( -0.04, 0.07, 1.52, -163.84)
dfp_em_step_func: 017: -log(p) = 4453918.0
after pass:transform: ( 1.81, 0.75, 0.02, -255.83)
                      ( -0.88, 2.17, -0.53, -151.47)
                      ( -0.04, 0.07, 1.52, -163.84)
outof QuasiNewtonEMA: 019: -log(p) = 4453918.0  tol 0.000000
final transform:
 1.809   0.747   0.019  -255.835;
-0.880   2.166  -0.530  -151.473;
-0.037   0.071   1.524  -163.836;
 0.000   0.000   0.000   1.000;

writing output transformation to transforms/talairach.lta...
registration took 23 minutes and 18 seconds.
#--------------------------------------
#@# CA Normalize Wed Jun 20 10:34:36 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_ca_normalize -mask brainmask.mgz nu.mgz /usr/local/freesurfer/average/RB_all_2006-02-15.gca transforms/talairach.lta norm.mgz 

using MR volume brainmask.mgz to mask input volume...
reading 1 input volumes
reading atlas from '/usr/local/freesurfer/average/RB_all_2006-02-15.gca'...
setting gca type = Normal gca type
reading transform from 'transforms/talairach.lta'...
INFO: compute sample coordinates transform
 1.809   0.747   0.019  -255.835;
-0.880   2.166  -0.530  -151.473;
-0.037   0.071   1.524  -163.836;
 0.000   0.000   0.000   1.000;
INFO: transform used
reading input volume from nu.mgz...
resetting wm mean[0]: 109 --> 116
resetting gm mean[0]: 68 --> 68
input volume #1 is the most T1-like
using real data threshold=17.0
using (99, 132, 166) as brain centroid...
mean wm in atlas = 116, using box (80,109,150) --> (117, 154,182) to find MRI wm
before smoothing, mri peak at 116
after smoothing, mri peak at 76, scaling input intensities by 1.526
scaling channel 0 by 1.52632
using 242283 sample points...
finding control points in Left_Cerebral_White_Matter....
found 38929 control points for structure...
bounding box (113, 189, 118) --> (155, 245, 226)
Left_Cerebral_White_Matter: limiting intensities to 128.0 --> 133.0
55 of 63 (87.3%) samples deleted
finding control points in Right_Cerebral_White_Matter....
found 38531 control points for structure...
bounding box (86, 183, 118) --> (130, 236, 226)
Right_Cerebral_White_Matter: limiting intensities to 122.0 --> 137.0
66 of 78 (84.6%) samples deleted
finding control points in Left_Cerebellum_White_Matter....
found 3087 control points for structure...
bounding box (117, 218, 132) --> (138, 242, 167)
finding control points in Right_Cerebellum_White_Matter....
found 2900 control points for structure...
bounding box (92, 217, 130) --> (119, 234, 167)
finding control points in Brain_Stem....
found 3414 control points for structure...
bounding box (105, 217, 153) --> (125, 245, 174)
using 141 total control points for intensity normalization...
bias field = 0.876 +- 0.035
0 of 18 control points discarded
writing normalized volume to norm.mgz...
freeing GCA...done.
normalization took 2 minutes and 4 seconds.
#--------------------------------------
#@# CA Reg Wed Jun 20 10:36:40 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_ca_register -align -nobigventricles -mask brainmask.mgz -T transforms/talairach.lta norm.mgz /usr/local/freesurfer/average/RB_all_2006-02-15.gca transforms/talairach.m3z 

renormalizing sequences with structure alignment, equivalent to:
	-renormalize
	-regularize_mean 0.500
	-regularize 0.500
not handling expanded ventricles...
using MR volume brainmask.mgz to mask input volume...
using previously computed transform transforms/talairach.lta
reading 1 input volumes...
logging results to talairach.log
reading input volume 'norm.mgz'...
reading GCA '/usr/local/freesurfer/average/RB_all_2006-02-15.gca'...
setting gca type = Normal gca type
freeing gibbs priors...done.
average std[0] = 5.0
label assignment complete, 0 changed (0.00%)
renormalizing by structure alignment....
renormalizing input #0
aligning Left_Pallidum...
gca peak = 0.11705 (100)
mri peak = 0.05768 (128)
aligning Right_Pallidum...
invalid transform detected (det=4.1451) 
gca peak = 0.13501 (106)
mri peak = 0.12826 (78)
aligning Right_Cerebral_White_Matter...
gca peak = 0.05652 (116)
mri peak = 0.04126 (80)
unreasonable value (87.0), not in range [90, 130] - rejecting
aligning Left_Cerebral_White_Matter...
gca peak = 0.05682 (116)
mri peak = 0.03031 (84)
Left_Cerebral_White_Matter (2): linear fit = 0.80 x + 0.0 (6947 voxels, overlap=0.018)
aligning Left_Lateral_Ventricle...
gca peak = 0.05153 (25)
mri peak = 0.05274 (89)
CSF peak too bright - rejecting
aligning Right_Lateral_Ventricle...
gca peak = 0.07557 (30)
mri peak = 0.03091 (86)
CSF peak too bright - rejecting
aligning Right_Hippocampus...
gca peak = 0.14805 (66)
mri peak = 0.03655 (120)
overlap = 0, overlap_threshold = 0.001
insufficient overlap 0.0000 in histograms - rejecting
aligning Left_Hippocampus...
gca peak = 0.17297 (66)
mri peak = 0.05224 (62)
aligning Left_Cerebral_Cortex...
gca peak = 0.05671 (68)
mri peak = 0.02007 (95)
uniform distribution in MR - rejecting arbitrary fit
aligning Right_Cerebral_Cortex...
gca peak = 0.06581 (70)
mri peak = 0.02161 (78)
uniform distribution in MR - rejecting arbitrary fit
aligning Right_Caudate...
gca peak = 0.07786 (81)
mri peak = 0.03143 (102)
unreasonable value (113.4), not in range [50, 100] - rejecting
aligning Left_Caudate...
gca peak = 0.07814 (81)
mri peak = 0.04876 (116)
unreasonable value (101.2), not in range [50, 100] - rejecting
aligning Left_Cerebellum_Cortex...
gca peak = 0.06033 (63)
mri peak = 0.02635 (93)
uniform distribution in MR - rejecting arbitrary fit
aligning Right_Cerebellum_Cortex...
gca peak = 0.05990 (64)
mri peak = 0.02902 (131)
uniform distribution in MR - rejecting arbitrary fit
aligning Left_Cerebellum_White_Matter...
gca peak = 0.09299 (92)
mri peak = 0.03012 (34)
Left_Cerebellum_White_Matter (7): linear fit = 1.02 x + 0.0 (1551 voxels, overlap=0.054)
aligning Right_Cerebellum_White_Matter...
invalid transform detected (det=0.1502) 
gca peak = 0.08899 (93)
mri peak = 0.01653 (59)
uniform distribution in MR - rejecting arbitrary fit
aligning Left_Amygdala...
invalid transform detected (det=0.1464) 
gca peak = 0.13911 (71)
mri peak = 0.07991 (113)
aligning Right_Amygdala...
gca peak = 0.15703 (65)
mri peak = 0.07349 (116)
overlap = 0, overlap_threshold = 0.001
insufficient overlap 0.0000 in histograms - rejecting
aligning Left_Thalamus_Proper...
gca peak = 0.06537 (100)
mri peak = 0.02137 (89)
uniform distribution in MR - rejecting arbitrary fit
aligning Right_Thalamus_Proper...
gca peak = 0.07019 (99)
mri peak = 0.05442 (94)
Right_Thalamus_Proper (49): linear fit = 0.98 x + 0.0 (490 voxels, overlap=0.430)
aligning Left_Putamen...
gca peak = 0.09468 (90)
mri peak = 0.03149 (21)
unreasonable value (119.3), not in range [60, 100] - rejecting
aligning Right_Putamen...
gca peak = 0.07184 (82)
mri peak = 0.02569 (149)
uniform distribution in MR - rejecting arbitrary fit
aligning Brain_Stem...
gca peak = 0.05597 (80)
mri peak = 0.02557 (123)
uniform distribution in MR - rejecting arbitrary fit
aligning Right_VentralDC...
gca peak = 0.05879 (97)
mri peak = 0.03021 (84)
Right_VentralDC (60): linear fit = 0.85 x + 0.0 (206 voxels, overlap=0.040)
aligning Left_VentralDC...
gca peak = 0.06932 (98)
mri peak = 0.03769 (118)
unreasonable value (112.7), not in range [80, 110] - rejecting
aligning Third_Ventricle...
invalid transform detected (det=0.0197) 
gca peak = 0.10105 (38)
mri peak = 0.10023 (34)
aligning Fourth_Ventricle...
gca peak = 0.10426 (32)
not using caudate to estimate GM means
setting label Left_Cerebral_White_Matter based on Right_Cerebral_White_Matter = 0.80 x +  0
setting label Left_Cerebellum_White_Matter based on Right_Cerebellum_White_Matter = 1.02 x +  0
setting label Right_Thalamus_Proper based on Left_Thalamus_Proper = 0.98 x +  0
setting label Right_VentralDC based on Left_VentralDC = 0.85 x +  0
uniform distribution in MR - rejecting arbitrary fit
estimating mean gm scale to be 1.00 x + 0.0
estimating mean wm scale to be 0.80 x + 0.0
estimating mean csf scale to be 1.00 x + 0.0
setting left cbm cortex = 1.00 x + 0.00
setting right cbm cortex = 1.00 x + 0.00
4 transforms computed
rms increased - undoing step...
0000: dt=0.000000, rms=1.286 (-0.029%), neg=0, invalid=1511, LIN
rms increased - undoing step...
0002: dt=0.000000, rms=1.286 (-0.030%), neg=0, invalid=1511, LIN
rms increased - undoing step...
0004: dt=0.000000, rms=1.286 (-0.030%), neg=0, invalid=1511, LIN
rms increased - undoing step...
0006: dt=0.000000, rms=1.286 (-0.030%), neg=0, invalid=1511, LIN
rms increased - undoing step...
0008: dt=0.000000, rms=1.286 (-0.039%), neg=0, invalid=1511, LIN
rms increased - undoing step...
0010: dt=0.000000, rms=1.286 (-0.047%), neg=0, invalid=1511, LIN
label assignment complete, 0 changed (0.00%)
**************** pass 1 of 1 ************************
setting smoothness coefficient to 0.01
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=1.272, neg=0, invalid=1511
0013: dt=417.371608, rms=1.216 (4.349%), neg=0, invalid=1511
0014: dt=260.279635, rms=1.186 (2.483%), neg=0, invalid=1511
0015: dt=86.256410, rms=1.170 (1.371%), neg=0, invalid=1511
0016: dt=129.472002, rms=1.165 (0.381%), neg=0, invalid=1511
0017: dt=129.472002, rms=1.150 (1.348%), neg=0, invalid=1511
0018: dt=122.352941, rms=1.146 (0.299%), neg=0, invalid=1511
0019: dt=129.472002, rms=1.130 (1.396%), neg=0, invalid=1511
0020: dt=32.368000, rms=1.129 (0.114%), neg=0, invalid=1511
0021: dt=110.976002, rms=1.126 (0.242%), neg=0, invalid=1511
0022: dt=155.286822, rms=1.115 (1.039%), neg=0, invalid=1511
0023: dt=26.666667, rms=1.113 (0.183%), neg=0, invalid=1511
0024: dt=23.120000, rms=1.112 (0.034%), neg=0, invalid=1511
0025: dt=23.120000, rms=1.112 (0.053%), neg=0, invalid=1511
0026: dt=23.120000, rms=1.110 (0.102%), neg=0, invalid=1511
0027: dt=23.120000, rms=1.108 (0.205%), neg=0, invalid=1511
0028: dt=23.120000, rms=1.104 (0.336%), neg=0, invalid=1511
0029: dt=23.120000, rms=1.100 (0.418%), neg=0, invalid=1511
0030: dt=23.120000, rms=1.095 (0.441%), neg=0, invalid=1511
0031: dt=23.120000, rms=1.090 (0.429%), neg=0, invalid=1511
0032: dt=23.120000, rms=1.086 (0.368%), neg=0, invalid=1511
0033: dt=23.120000, rms=1.083 (0.339%), neg=0, invalid=1511
0034: dt=23.120000, rms=1.078 (0.379%), neg=0, invalid=1511
0035: dt=23.120000, rms=1.077 (0.134%), neg=0, invalid=1511
0036: dt=23.120000, rms=1.075 (0.196%), neg=0, invalid=1511
0037: dt=23.120000, rms=1.072 (0.228%), neg=0, invalid=1511
0038: dt=23.120000, rms=1.070 (0.241%), neg=0, invalid=1511
0039: dt=23.120000, rms=1.067 (0.258%), neg=0, invalid=1511
0040: dt=23.120000, rms=1.064 (0.281%), neg=0, invalid=1511
0041: dt=23.120000, rms=1.061 (0.323%), neg=0, invalid=1511
0042: dt=23.120000, rms=1.057 (0.382%), neg=0, invalid=1511
0043: dt=23.120000, rms=1.053 (0.381%), neg=0, invalid=1511
0044: dt=23.120000, rms=1.052 (0.046%), neg=0, invalid=1511
0045: dt=11.560000, rms=1.052 (0.022%), neg=0, invalid=1511
0046: dt=1.445000, rms=1.052 (0.003%), neg=0, invalid=1511
0047: dt=0.505750, rms=1.052 (0.001%), neg=0, invalid=1511
0048: dt=0.505750, rms=1.052 (0.001%), neg=0, invalid=1511
0049: dt=0.126438, rms=1.052 (0.000%), neg=0, invalid=1511
0050: dt=0.007902, rms=1.052 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=1.052, neg=0, invalid=1511
0051: dt=117.042802, rms=1.045 (0.654%), neg=0, invalid=1511
0052: dt=51.636364, rms=1.042 (0.294%), neg=0, invalid=1511
0053: dt=295.936004, rms=1.031 (1.045%), neg=0, invalid=1511
0054: dt=32.213333, rms=1.027 (0.423%), neg=0, invalid=1511
0055: dt=32.368000, rms=1.026 (0.038%), neg=0, invalid=1511
0056: dt=32.368000, rms=1.026 (0.058%), neg=0, invalid=1511
0057: dt=32.368000, rms=1.024 (0.123%), neg=0, invalid=1511
0058: dt=32.368000, rms=1.022 (0.184%), neg=0, invalid=1511
0059: dt=32.368000, rms=1.020 (0.241%), neg=0, invalid=1511
0060: dt=32.368000, rms=1.017 (0.279%), neg=0, invalid=1511
0061: dt=32.368000, rms=1.014 (0.302%), neg=0, invalid=1511
0062: dt=32.368000, rms=1.011 (0.315%), neg=0, invalid=1511
0063: dt=32.368000, rms=1.008 (0.308%), neg=0, invalid=1511
0064: dt=32.368000, rms=1.005 (0.304%), neg=0, invalid=1511
0065: dt=32.368000, rms=1.002 (0.300%), neg=0, invalid=1511
0066: dt=32.368000, rms=0.999 (0.294%), neg=0, invalid=1511
0067: dt=32.368000, rms=0.996 (0.278%), neg=0, invalid=1511
0068: dt=32.368000, rms=0.993 (0.268%), neg=0, invalid=1511
0069: dt=32.368000, rms=0.991 (0.258%), neg=0, invalid=1511
0070: dt=32.368000, rms=0.988 (0.245%), neg=0, invalid=1511
0071: dt=32.368000, rms=0.986 (0.235%), neg=0, invalid=1511
0072: dt=32.368000, rms=0.984 (0.248%), neg=0, invalid=1511
0073: dt=32.368000, rms=0.981 (0.253%), neg=0, invalid=1511
0074: dt=32.368000, rms=0.979 (0.235%), neg=0, invalid=1511
0075: dt=32.368000, rms=0.977 (0.219%), neg=0, invalid=1511
0076: dt=32.368000, rms=0.975 (0.208%), neg=0, invalid=1511
0077: dt=32.368000, rms=0.973 (0.212%), neg=0, invalid=1511
0078: dt=32.368000, rms=0.970 (0.221%), neg=0, invalid=1511
0079: dt=32.368000, rms=0.968 (0.203%), neg=0, invalid=1511
0080: dt=32.368000, rms=0.967 (0.191%), neg=0, invalid=1511
0081: dt=32.368000, rms=0.965 (0.169%), neg=0, invalid=1511
0082: dt=32.368000, rms=0.963 (0.161%), neg=0, invalid=1511
0083: dt=32.368000, rms=0.962 (0.161%), neg=0, invalid=1511
0084: dt=32.368000, rms=0.960 (0.160%), neg=0, invalid=1511
0085: dt=32.368000, rms=0.959 (0.141%), neg=0, invalid=1511
0086: dt=32.368000, rms=0.958 (0.118%), neg=0, invalid=1511
0087: dt=32.368000, rms=0.957 (0.095%), neg=0, invalid=1511
0088: dt=32.368000, rms=0.956 (0.084%), neg=0, invalid=1511
0089: dt=32.368000, rms=0.955 (0.079%), neg=0, invalid=1511
0090: dt=32.368000, rms=0.955 (0.084%), neg=0, invalid=1511
0091: dt=32.368000, rms=0.954 (0.096%), neg=0, invalid=1511
0092: dt=32.368000, rms=0.953 (0.104%), neg=0, invalid=1511
0093: dt=32.368000, rms=0.952 (0.104%), neg=0, invalid=1511
0094: dt=32.368000, rms=0.951 (0.113%), neg=0, invalid=1511
0095: dt=32.368000, rms=0.950 (0.101%), neg=0, invalid=1511
0096: dt=32.368000, rms=0.949 (0.105%), neg=0, invalid=1511
0097: dt=32.368000, rms=0.948 (0.108%), neg=0, invalid=1511
0098: dt=32.368000, rms=0.947 (0.103%), neg=0, invalid=1511
0099: dt=32.368000, rms=0.946 (0.099%), neg=0, invalid=1511
0100: dt=32.368000, rms=0.945 (0.094%), neg=0, invalid=1511
0101: dt=32.368000, rms=0.944 (0.089%), neg=0, invalid=1511
0102: dt=32.368000, rms=0.943 (0.090%), neg=0, invalid=1511
0103: dt=32.368000, rms=0.942 (0.094%), neg=0, invalid=1511
0104: dt=32.368000, rms=0.941 (0.093%), neg=0, invalid=1511
0105: dt=32.368000, rms=0.940 (0.089%), neg=0, invalid=1511
0106: dt=32.368000, rms=0.940 (0.082%), neg=0, invalid=1511
0107: dt=32.368000, rms=0.939 (0.074%), neg=0, invalid=1511
0108: dt=32.368000, rms=0.938 (0.070%), neg=0, invalid=1511
0109: dt=32.368000, rms=0.938 (0.060%), neg=0, invalid=1511
0110: dt=32.368000, rms=0.937 (0.058%), neg=0, invalid=1511
0111: dt=32.368000, rms=0.937 (0.050%), neg=0, invalid=1511
0112: dt=32.368000, rms=0.936 (0.052%), neg=0, invalid=1511
0113: dt=32.368000, rms=0.936 (0.055%), neg=0, invalid=1511
0114: dt=32.368000, rms=0.935 (0.057%), neg=0, invalid=1511
0115: dt=32.368000, rms=0.935 (0.075%), neg=0, invalid=1511
0116: dt=32.368000, rms=0.934 (0.083%), neg=0, invalid=1511
0117: dt=32.368000, rms=0.933 (0.090%), neg=0, invalid=1511
0118: dt=32.368000, rms=0.932 (0.092%), neg=0, invalid=1511
0119: dt=32.368000, rms=0.931 (0.090%), neg=0, invalid=1511
0120: dt=32.368000, rms=0.930 (0.085%), neg=0, invalid=1511
0121: dt=32.368000, rms=0.930 (0.079%), neg=0, invalid=1511
0122: dt=32.368000, rms=0.929 (0.078%), neg=0, invalid=1511
0123: dt=32.368000, rms=0.928 (0.072%), neg=0, invalid=1511
0124: dt=32.368000, rms=0.928 (0.067%), neg=0, invalid=1511
0125: dt=32.368000, rms=0.927 (0.064%), neg=0, invalid=1511
0126: dt=32.368000, rms=0.926 (0.066%), neg=0, invalid=1511
0127: dt=32.368000, rms=0.926 (0.073%), neg=0, invalid=1511
0128: dt=32.368000, rms=0.925 (0.074%), neg=0, invalid=1511
0129: dt=32.368000, rms=0.924 (0.074%), neg=0, invalid=1511
0130: dt=32.368000, rms=0.924 (0.075%), neg=0, invalid=1511
0131: dt=32.368000, rms=0.923 (0.076%), neg=0, invalid=1511
0132: dt=32.368000, rms=0.922 (0.080%), neg=0, invalid=1511
0133: dt=32.368000, rms=0.922 (0.082%), neg=0, invalid=1511
0134: dt=32.368000, rms=0.921 (0.081%), neg=0, invalid=1511
0135: dt=32.368000, rms=0.920 (0.076%), neg=0, invalid=1511
0136: dt=32.368000, rms=0.919 (0.073%), neg=0, invalid=1511
0137: dt=32.368000, rms=0.919 (0.067%), neg=0, invalid=1511
0138: dt=32.368000, rms=0.918 (0.062%), neg=0, invalid=1511
0139: dt=32.368000, rms=0.918 (0.057%), neg=0, invalid=1511
0140: dt=32.368000, rms=0.917 (0.056%), neg=0, invalid=1511
0141: dt=32.368000, rms=0.917 (0.056%), neg=0, invalid=1511
0142: dt=32.368000, rms=0.916 (0.052%), neg=0, invalid=1511
0143: dt=32.368000, rms=0.916 (0.045%), neg=0, invalid=1511
0144: dt=32.368000, rms=0.915 (0.042%), neg=0, invalid=1511
0145: dt=32.368000, rms=0.915 (0.049%), neg=0, invalid=1511
0146: dt=32.368000, rms=0.914 (0.051%), neg=0, invalid=1511
0147: dt=32.368000, rms=0.914 (0.054%), neg=0, invalid=1511
0148: dt=32.368000, rms=0.914 (0.051%), neg=0, invalid=1511
0149: dt=32.368000, rms=0.913 (0.051%), neg=0, invalid=1511
0150: dt=32.368000, rms=0.913 (0.048%), neg=0, invalid=1511
0151: dt=32.368000, rms=0.912 (0.046%), neg=0, invalid=1511
0152: dt=32.368000, rms=0.912 (0.045%), neg=0, invalid=1511
0153: dt=32.368000, rms=0.911 (0.049%), neg=0, invalid=1511
0154: dt=32.368000, rms=0.911 (0.045%), neg=0, invalid=1511
0155: dt=32.368000, rms=0.910 (0.050%), neg=0, invalid=1511
0156: dt=32.368000, rms=0.910 (0.042%), neg=0, invalid=1511
0157: dt=32.368000, rms=0.910 (0.038%), neg=0, invalid=1511
0158: dt=32.368000, rms=0.909 (0.040%), neg=0, invalid=1511
0159: dt=32.368000, rms=0.909 (0.039%), neg=0, invalid=1511
0160: dt=32.368000, rms=0.909 (0.037%), neg=0, invalid=1511
0161: dt=32.368000, rms=0.908 (0.040%), neg=0, invalid=1511
0162: dt=32.368000, rms=0.908 (0.037%), neg=0, invalid=1511
0163: dt=32.368000, rms=0.908 (0.034%), neg=0, invalid=1511
0164: dt=32.368000, rms=0.907 (0.033%), neg=0, invalid=1511
0165: dt=32.368000, rms=0.907 (0.034%), neg=0, invalid=1511
0166: dt=32.368000, rms=0.907 (0.032%), neg=0, invalid=1511
0167: dt=32.368000, rms=0.907 (0.030%), neg=0, invalid=1511
0168: dt=32.368000, rms=0.906 (0.027%), neg=0, invalid=1511
0169: dt=32.368000, rms=0.906 (0.030%), neg=0, invalid=1511
0170: dt=32.368000, rms=0.906 (0.029%), neg=0, invalid=1511
0171: dt=32.368000, rms=0.905 (0.030%), neg=0, invalid=1511
0172: dt=32.368000, rms=0.905 (0.037%), neg=0, invalid=1511
0173: dt=32.368000, rms=0.905 (0.037%), neg=0, invalid=1511
0174: dt=32.368000, rms=0.905 (0.033%), neg=0, invalid=1511
0175: dt=32.368000, rms=0.904 (0.036%), neg=0, invalid=1511
0176: dt=32.368000, rms=0.904 (0.032%), neg=0, invalid=1511
0177: dt=32.368000, rms=0.904 (0.035%), neg=0, invalid=1511
0178: dt=32.368000, rms=0.903 (0.037%), neg=0, invalid=1511
0179: dt=32.368000, rms=0.903 (0.030%), neg=0, invalid=1511
0180: dt=32.368000, rms=0.903 (0.027%), neg=0, invalid=1511
0181: dt=32.368000, rms=0.903 (0.020%), neg=0, invalid=1511
0182: dt=32.368000, rms=0.902 (0.024%), neg=0, invalid=1511
0183: dt=32.368000, rms=0.902 (0.024%), neg=0, invalid=1511
0184: dt=32.368000, rms=0.902 (0.027%), neg=0, invalid=1511
0185: dt=32.368000, rms=0.902 (0.024%), neg=0, invalid=1511
0186: dt=32.368000, rms=0.901 (0.024%), neg=0, invalid=1511
0187: dt=129.472002, rms=0.901 (-0.000%), neg=0, invalid=1511
setting smoothness coefficient to 0.03
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.906, neg=0, invalid=1511
0188: dt=29.931034, rms=0.903 (0.345%), neg=0, invalid=1511
0189: dt=8.769231, rms=0.902 (0.105%), neg=0, invalid=1511
0190: dt=9.072000, rms=0.899 (0.290%), neg=0, invalid=1511
0191: dt=25.920000, rms=0.896 (0.371%), neg=0, invalid=1511
0192: dt=1.296000, rms=0.896 (-0.002%), neg=0, invalid=1511
0193: dt=1.296000, rms=0.896 (0.001%), neg=0, invalid=1511
0194: dt=1.296000, rms=0.896 (0.007%), neg=0, invalid=1511
0195: dt=1.296000, rms=0.896 (0.020%), neg=0, invalid=1511
0196: dt=1.296000, rms=0.895 (0.033%), neg=0, invalid=1511
0197: dt=1.296000, rms=0.895 (0.051%), neg=0, invalid=1511
0198: dt=1.296000, rms=0.894 (0.068%), neg=0, invalid=1511
0199: dt=1.296000, rms=0.894 (0.084%), neg=0, invalid=1511
0200: dt=1.296000, rms=0.893 (0.096%), neg=0, invalid=1511
0201: dt=1.296000, rms=0.892 (0.100%), neg=0, invalid=1511
0202: dt=1.296000, rms=0.891 (0.098%), neg=0, invalid=1511
0203: dt=1.296000, rms=0.890 (0.097%), neg=0, invalid=1511
0204: dt=1.296000, rms=0.889 (0.083%), neg=0, invalid=1511
0205: dt=1.296000, rms=0.889 (0.074%), neg=0, invalid=1511
0206: dt=1.296000, rms=0.888 (0.060%), neg=0, invalid=1511
0207: dt=1.296000, rms=0.888 (0.055%), neg=0, invalid=1511
0208: dt=1.296000, rms=0.887 (0.056%), neg=0, invalid=1511
0209: dt=1.296000, rms=0.887 (0.068%), neg=0, invalid=1511
0210: dt=1.296000, rms=0.886 (0.078%), neg=0, invalid=1511
0211: dt=1.296000, rms=0.885 (0.088%), neg=0, invalid=1511
0212: dt=1.296000, rms=0.884 (0.102%), neg=0, invalid=1511
0213: dt=1.296000, rms=0.883 (0.105%), neg=0, invalid=1511
0214: dt=1.296000, rms=0.882 (0.103%), neg=0, invalid=1511
0215: dt=1.296000, rms=0.882 (0.094%), neg=0, invalid=1511
0216: dt=1.296000, rms=0.881 (0.085%), neg=0, invalid=1511
0217: dt=1.296000, rms=0.880 (0.073%), neg=0, invalid=1511
0218: dt=1.296000, rms=0.880 (0.061%), neg=0, invalid=1511
0219: dt=1.296000, rms=0.879 (0.055%), neg=0, invalid=1511
0220: dt=1.296000, rms=0.879 (0.054%), neg=0, invalid=1511
0221: dt=1.296000, rms=0.878 (0.061%), neg=0, invalid=1511
0222: dt=1.296000, rms=0.878 (0.069%), neg=0, invalid=1511
0223: dt=1.296000, rms=0.877 (0.076%), neg=0, invalid=1511
0224: dt=1.296000, rms=0.876 (0.078%), neg=0, invalid=1511
0225: dt=1.296000, rms=0.875 (0.084%), neg=0, invalid=1511
0226: dt=1.296000, rms=0.875 (0.077%), neg=0, invalid=1511
0227: dt=1.296000, rms=0.874 (0.074%), neg=0, invalid=1511
0228: dt=1.296000, rms=0.874 (0.069%), neg=0, invalid=1511
0229: dt=1.296000, rms=0.873 (0.067%), neg=0, invalid=1511
0230: dt=1.296000, rms=0.872 (0.063%), neg=0, invalid=1511
0231: dt=1.296000, rms=0.872 (0.065%), neg=0, invalid=1511
0232: dt=1.296000, rms=0.871 (0.072%), neg=0, invalid=1511
0233: dt=1.296000, rms=0.871 (0.071%), neg=0, invalid=1511
0234: dt=1.296000, rms=0.870 (0.075%), neg=0, invalid=1511
0235: dt=1.296000, rms=0.869 (0.074%), neg=0, invalid=1511
0236: dt=1.296000, rms=0.869 (0.077%), neg=0, invalid=1511
0237: dt=1.296000, rms=0.868 (0.077%), neg=0, invalid=1511
0238: dt=1.296000, rms=0.867 (0.072%), neg=0, invalid=1511
0239: dt=1.296000, rms=0.867 (0.069%), neg=0, invalid=1511
0240: dt=1.296000, rms=0.866 (0.068%), neg=0, invalid=1511
0241: dt=1.296000, rms=0.866 (0.065%), neg=0, invalid=1511
0242: dt=1.296000, rms=0.865 (0.062%), neg=0, invalid=1511
0243: dt=1.296000, rms=0.865 (0.059%), neg=0, invalid=1511
0244: dt=1.296000, rms=0.864 (0.062%), neg=0, invalid=1511
0245: dt=1.296000, rms=0.863 (0.064%), neg=0, invalid=1511
0246: dt=1.296000, rms=0.863 (0.067%), neg=0, invalid=1511
0247: dt=1.296000, rms=0.862 (0.065%), neg=0, invalid=1511
0248: dt=1.296000, rms=0.862 (0.058%), neg=0, invalid=1511
0249: dt=1.296000, rms=0.861 (0.057%), neg=0, invalid=1511
0250: dt=1.296000, rms=0.861 (0.053%), neg=0, invalid=1511
0251: dt=1.296000, rms=0.860 (0.051%), neg=0, invalid=1511
0252: dt=1.296000, rms=0.860 (0.046%), neg=0, invalid=1511
0253: dt=1.296000, rms=0.860 (0.048%), neg=0, invalid=1511
0254: dt=1.296000, rms=0.859 (0.048%), neg=0, invalid=1511
0255: dt=1.296000, rms=0.859 (0.045%), neg=0, invalid=1511
0256: dt=1.296000, rms=0.858 (0.046%), neg=0, invalid=1511
0257: dt=1.296000, rms=0.858 (0.049%), neg=0, invalid=1511
0258: dt=1.296000, rms=0.858 (0.046%), neg=0, invalid=1511
0259: dt=1.296000, rms=0.857 (0.043%), neg=0, invalid=1511
0260: dt=1.296000, rms=0.857 (0.048%), neg=0, invalid=1511
0261: dt=1.296000, rms=0.856 (0.045%), neg=0, invalid=1511
0262: dt=1.296000, rms=0.856 (0.045%), neg=0, invalid=1511
0263: dt=1.296000, rms=0.856 (0.038%), neg=0, invalid=1511
0264: dt=1.296000, rms=0.856 (0.007%), neg=0, invalid=1511
0265: dt=1.296000, rms=0.856 (0.007%), neg=0, invalid=1511
0266: dt=1.296000, rms=0.856 (0.013%), neg=0, invalid=1511
0267: dt=1.296000, rms=0.855 (0.017%), neg=0, invalid=1511
0268: dt=1.296000, rms=0.855 (0.022%), neg=0, invalid=1511
0269: dt=1.296000, rms=0.855 (0.025%), neg=0, invalid=1511
0270: dt=1.296000, rms=0.855 (0.028%), neg=0, invalid=1511
0271: dt=1.296000, rms=0.854 (0.035%), neg=0, invalid=1511
0272: dt=1.296000, rms=0.854 (0.036%), neg=0, invalid=1511
0273: dt=1.296000, rms=0.854 (0.034%), neg=0, invalid=1511
0274: dt=1.296000, rms=0.853 (0.037%), neg=0, invalid=1511
0275: dt=1.296000, rms=0.853 (0.038%), neg=0, invalid=1511
0276: dt=1.296000, rms=0.853 (0.040%), neg=0, invalid=1511
0277: dt=1.296000, rms=0.852 (0.041%), neg=0, invalid=1511
0278: dt=1.296000, rms=0.852 (0.041%), neg=0, invalid=1511
0279: dt=1.296000, rms=0.852 (0.039%), neg=0, invalid=1511
0280: dt=1.296000, rms=0.851 (0.040%), neg=0, invalid=1511
0281: dt=1.296000, rms=0.851 (0.040%), neg=0, invalid=1511
0282: dt=1.296000, rms=0.851 (0.043%), neg=0, invalid=1511
0283: dt=1.296000, rms=0.850 (0.041%), neg=0, invalid=1511
0284: dt=1.296000, rms=0.850 (0.040%), neg=0, invalid=1511
0285: dt=1.296000, rms=0.850 (0.003%), neg=0, invalid=1511
0286: dt=1.296000, rms=0.850 (0.008%), neg=0, invalid=1511
0287: dt=1.296000, rms=0.850 (0.011%), neg=0, invalid=1511
0288: dt=1.296000, rms=0.850 (0.014%), neg=0, invalid=1511
0289: dt=1.296000, rms=0.850 (0.016%), neg=0, invalid=1511
0290: dt=1.296000, rms=0.849 (0.016%), neg=0, invalid=1511
0291: dt=9.072000, rms=0.849 (0.017%), neg=0, invalid=1511
0292: dt=9.072000, rms=0.849 (0.028%), neg=0, invalid=1511
0293: dt=6.480000, rms=0.849 (0.018%), neg=0, invalid=1511
0294: dt=0.810000, rms=0.849 (0.003%), neg=0, invalid=1511
0295: dt=0.405000, rms=0.849 (0.001%), neg=0, invalid=1511
0296: dt=0.101250, rms=0.849 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.849, neg=0, invalid=1511
0297: dt=41.185185, rms=0.846 (0.324%), neg=0, invalid=1511
0298: dt=25.920000, rms=0.840 (0.755%), neg=0, invalid=1511
0299: dt=2.268000, rms=0.840 (0.007%), neg=0, invalid=1511
0300: dt=2.268000, rms=0.840 (0.005%), neg=0, invalid=1511
0301: dt=2.268000, rms=0.840 (0.010%), neg=0, invalid=1511
0302: dt=2.268000, rms=0.839 (0.020%), neg=0, invalid=1511
0303: dt=2.268000, rms=0.839 (0.042%), neg=0, invalid=1511
0304: dt=2.268000, rms=0.838 (0.081%), neg=0, invalid=1511
0305: dt=2.268000, rms=0.837 (0.122%), neg=0, invalid=1511
0306: dt=2.268000, rms=0.836 (0.158%), neg=0, invalid=1511
0307: dt=2.268000, rms=0.835 (0.175%), neg=0, invalid=1511
0308: dt=2.268000, rms=0.833 (0.176%), neg=0, invalid=1511
0309: dt=2.268000, rms=0.832 (0.159%), neg=0, invalid=1511
0310: dt=2.268000, rms=0.831 (0.138%), neg=0, invalid=1511
0311: dt=2.268000, rms=0.830 (0.116%), neg=0, invalid=1511
0312: dt=2.268000, rms=0.829 (0.106%), neg=0, invalid=1511
0313: dt=2.268000, rms=0.828 (0.092%), neg=0, invalid=1511
0314: dt=2.268000, rms=0.827 (0.096%), neg=0, invalid=1511
0315: dt=2.268000, rms=0.826 (0.099%), neg=0, invalid=1511
0316: dt=2.268000, rms=0.826 (0.108%), neg=0, invalid=1511
0317: dt=2.268000, rms=0.825 (0.118%), neg=0, invalid=1511
0318: dt=2.268000, rms=0.824 (0.115%), neg=0, invalid=1511
0319: dt=2.268000, rms=0.823 (0.118%), neg=0, invalid=1511
0320: dt=2.268000, rms=0.822 (0.109%), neg=0, invalid=1511
0321: dt=2.268000, rms=0.821 (0.105%), neg=0, invalid=1511
0322: dt=2.268000, rms=0.820 (0.096%), neg=0, invalid=1511
0323: dt=2.268000, rms=0.819 (0.094%), neg=0, invalid=1511
0324: dt=2.268000, rms=0.819 (0.091%), neg=0, invalid=1511
0325: dt=2.268000, rms=0.818 (0.093%), neg=0, invalid=1511
0326: dt=2.268000, rms=0.817 (0.092%), neg=0, invalid=1511
0327: dt=2.268000, rms=0.816 (0.094%), neg=0, invalid=1511
0328: dt=2.268000, rms=0.816 (0.091%), neg=0, invalid=1511
0329: dt=2.268000, rms=0.815 (0.086%), neg=0, invalid=1511
0330: dt=2.268000, rms=0.814 (0.079%), neg=0, invalid=1511
0331: dt=2.268000, rms=0.814 (0.073%), neg=0, invalid=1511
0332: dt=2.268000, rms=0.813 (0.070%), neg=0, invalid=1511
0333: dt=2.268000, rms=0.813 (0.067%), neg=0, invalid=1511
0334: dt=2.268000, rms=0.812 (0.062%), neg=0, invalid=1511
0335: dt=2.268000, rms=0.811 (0.063%), neg=0, invalid=1511
0336: dt=2.268000, rms=0.811 (0.062%), neg=0, invalid=1511
0337: dt=2.268000, rms=0.810 (0.067%), neg=0, invalid=1511
0338: dt=2.268000, rms=0.810 (0.067%), neg=0, invalid=1511
0339: dt=2.268000, rms=0.809 (0.063%), neg=0, invalid=1511
0340: dt=2.268000, rms=0.809 (0.058%), neg=0, invalid=1511
0341: dt=2.268000, rms=0.808 (0.060%), neg=0, invalid=1511
0342: dt=2.268000, rms=0.808 (0.053%), neg=0, invalid=1511
0343: dt=2.268000, rms=0.808 (0.049%), neg=0, invalid=1511
0344: dt=2.268000, rms=0.807 (0.046%), neg=0, invalid=1511
0345: dt=2.268000, rms=0.807 (0.047%), neg=0, invalid=1511
0346: dt=2.268000, rms=0.806 (0.047%), neg=0, invalid=1511
0347: dt=2.268000, rms=0.806 (0.047%), neg=0, invalid=1511
0348: dt=2.268000, rms=0.806 (0.051%), neg=0, invalid=1511
0349: dt=2.268000, rms=0.805 (0.049%), neg=0, invalid=1511
0350: dt=2.268000, rms=0.805 (0.053%), neg=0, invalid=1511
0351: dt=2.268000, rms=0.804 (0.052%), neg=0, invalid=1511
0352: dt=2.268000, rms=0.804 (0.051%), neg=0, invalid=1511
0353: dt=2.268000, rms=0.804 (0.051%), neg=0, invalid=1511
0354: dt=2.268000, rms=0.803 (0.045%), neg=0, invalid=1511
0355: dt=2.268000, rms=0.803 (0.044%), neg=0, invalid=1511
0356: dt=2.268000, rms=0.803 (0.049%), neg=0, invalid=1511
0357: dt=2.268000, rms=0.802 (0.043%), neg=0, invalid=1511
0358: dt=2.268000, rms=0.802 (0.046%), neg=0, invalid=1511
0359: dt=2.268000, rms=0.801 (0.045%), neg=0, invalid=1511
0360: dt=2.268000, rms=0.801 (0.042%), neg=0, invalid=1511
0361: dt=2.268000, rms=0.801 (0.044%), neg=0, invalid=1511
0362: dt=2.268000, rms=0.800 (0.042%), neg=0, invalid=1511
0363: dt=2.268000, rms=0.800 (0.041%), neg=0, invalid=1511
0364: dt=2.268000, rms=0.800 (0.044%), neg=0, invalid=1511
0365: dt=2.268000, rms=0.799 (0.043%), neg=0, invalid=1511
0366: dt=2.268000, rms=0.799 (0.041%), neg=0, invalid=1511
0367: dt=2.268000, rms=0.799 (0.044%), neg=0, invalid=1511
0368: dt=2.268000, rms=0.798 (0.042%), neg=0, invalid=1511
0369: dt=2.268000, rms=0.798 (0.045%), neg=0, invalid=1511
0370: dt=2.268000, rms=0.798 (0.044%), neg=0, invalid=1511
0371: dt=2.268000, rms=0.797 (0.042%), neg=0, invalid=1511
0372: dt=2.268000, rms=0.797 (0.040%), neg=0, invalid=1511
0373: dt=2.268000, rms=0.797 (0.039%), neg=0, invalid=1511
0374: dt=2.268000, rms=0.796 (0.039%), neg=0, invalid=1511
0375: dt=2.268000, rms=0.796 (0.037%), neg=0, invalid=1511
0376: dt=2.268000, rms=0.796 (0.039%), neg=0, invalid=1511
0377: dt=2.268000, rms=0.795 (0.038%), neg=0, invalid=1511
0378: dt=2.268000, rms=0.795 (0.035%), neg=0, invalid=1511
0379: dt=2.268000, rms=0.795 (0.039%), neg=0, invalid=1511
0380: dt=2.268000, rms=0.795 (0.039%), neg=0, invalid=1511
0381: dt=2.268000, rms=0.794 (0.037%), neg=0, invalid=1511
0382: dt=2.268000, rms=0.794 (0.035%), neg=0, invalid=1511
0383: dt=2.268000, rms=0.794 (0.036%), neg=0, invalid=1511
0384: dt=2.268000, rms=0.793 (0.034%), neg=0, invalid=1511
0385: dt=2.268000, rms=0.793 (0.031%), neg=0, invalid=1511
0386: dt=2.268000, rms=0.793 (0.033%), neg=0, invalid=1511
0387: dt=2.268000, rms=0.793 (0.029%), neg=0, invalid=1511
0388: dt=2.268000, rms=0.792 (0.034%), neg=0, invalid=1511
0389: dt=2.268000, rms=0.792 (0.032%), neg=0, invalid=1511
0390: dt=2.268000, rms=0.792 (0.030%), neg=0, invalid=1511
0391: dt=2.268000, rms=0.792 (0.031%), neg=0, invalid=1511
0392: dt=2.268000, rms=0.791 (0.031%), neg=0, invalid=1511
0393: dt=2.268000, rms=0.791 (0.031%), neg=0, invalid=1511
0394: dt=2.268000, rms=0.791 (0.026%), neg=0, invalid=1511
0395: dt=2.268000, rms=0.791 (0.031%), neg=0, invalid=1511
0396: dt=2.268000, rms=0.790 (0.035%), neg=0, invalid=1511
0397: dt=2.268000, rms=0.790 (0.029%), neg=0, invalid=1511
0398: dt=2.268000, rms=0.790 (0.027%), neg=0, invalid=1511
0399: dt=2.268000, rms=0.790 (0.030%), neg=0, invalid=1511
0400: dt=2.268000, rms=0.790 (0.024%), neg=0, invalid=1511
0401: dt=2.268000, rms=0.789 (0.029%), neg=0, invalid=1511
0402: dt=2.268000, rms=0.789 (0.029%), neg=0, invalid=1511
0403: dt=2.268000, rms=0.789 (0.026%), neg=0, invalid=1511
0404: dt=2.268000, rms=0.789 (0.026%), neg=0, invalid=1511
0405: dt=2.268000, rms=0.788 (0.034%), neg=0, invalid=1511
0406: dt=2.268000, rms=0.788 (0.030%), neg=0, invalid=1511
0407: dt=2.268000, rms=0.788 (0.028%), neg=0, invalid=1511
0408: dt=2.268000, rms=0.788 (0.030%), neg=0, invalid=1511
0409: dt=2.268000, rms=0.788 (0.031%), neg=0, invalid=1511
0410: dt=2.268000, rms=0.787 (0.027%), neg=0, invalid=1511
0411: dt=2.268000, rms=0.787 (0.029%), neg=0, invalid=1511
0412: dt=2.268000, rms=0.787 (0.026%), neg=0, invalid=1511
0413: dt=2.268000, rms=0.787 (0.029%), neg=0, invalid=1511
0414: dt=2.268000, rms=0.786 (0.023%), neg=0, invalid=1511
0415: dt=2.268000, rms=0.786 (0.028%), neg=0, invalid=1511
0416: dt=2.268000, rms=0.786 (0.023%), neg=0, invalid=1511
0417: dt=2.268000, rms=0.786 (0.022%), neg=0, invalid=1511
0418: dt=145.152002, rms=0.785 (0.131%), neg=0, invalid=1511
0419: dt=36.288001, rms=0.784 (0.066%), neg=0, invalid=1511
0420: dt=9.072000, rms=0.784 (0.007%), neg=0, invalid=1511
0421: dt=9.072000, rms=0.784 (0.017%), neg=0, invalid=1511
0422: dt=9.072000, rms=0.784 (0.029%), neg=0, invalid=1511
0423: dt=9.072000, rms=0.784 (0.027%), neg=0, invalid=1511
0424: dt=9.072000, rms=0.783 (0.034%), neg=0, invalid=1511
0425: dt=9.072000, rms=0.783 (0.042%), neg=0, invalid=1511
0426: dt=9.072000, rms=0.783 (0.050%), neg=0, invalid=1511
0427: dt=9.072000, rms=0.782 (0.053%), neg=0, invalid=1511
0428: dt=9.072000, rms=0.782 (0.054%), neg=0, invalid=1511
0429: dt=9.072000, rms=0.781 (0.054%), neg=0, invalid=1511
0430: dt=9.072000, rms=0.781 (0.064%), neg=0, invalid=1511
0431: dt=9.072000, rms=0.780 (0.067%), neg=0, invalid=1511
0432: dt=9.072000, rms=0.780 (0.068%), neg=0, invalid=1511
0433: dt=9.072000, rms=0.780 (0.054%), neg=0, invalid=1511
0434: dt=9.072000, rms=0.779 (0.052%), neg=0, invalid=1511
0435: dt=9.072000, rms=0.779 (0.058%), neg=0, invalid=1511
0436: dt=9.072000, rms=0.778 (0.057%), neg=0, invalid=1511
0437: dt=9.072000, rms=0.778 (0.066%), neg=0, invalid=1511
0438: dt=9.072000, rms=0.777 (0.069%), neg=0, invalid=1511
0439: dt=9.072000, rms=0.777 (0.068%), neg=0, invalid=1511
0440: dt=9.072000, rms=0.776 (0.065%), neg=0, invalid=1511
0441: dt=9.072000, rms=0.776 (0.060%), neg=0, invalid=1511
0442: dt=9.072000, rms=0.775 (0.059%), neg=0, invalid=1511
0443: dt=9.072000, rms=0.775 (0.063%), neg=0, invalid=1511
0444: dt=9.072000, rms=0.774 (0.067%), neg=0, invalid=1511
0445: dt=9.072000, rms=0.774 (0.067%), neg=0, invalid=1511
0446: dt=9.072000, rms=0.773 (0.058%), neg=0, invalid=1511
0447: dt=9.072000, rms=0.773 (0.051%), neg=0, invalid=1511
0448: dt=9.072000, rms=0.772 (0.048%), neg=0, invalid=1511
0449: dt=9.072000, rms=0.772 (0.050%), neg=0, invalid=1511
0450: dt=9.072000, rms=0.772 (0.055%), neg=0, invalid=1511
0451: dt=9.072000, rms=0.771 (0.069%), neg=0, invalid=1511
0452: dt=9.072000, rms=0.771 (0.006%), neg=0, invalid=1511
0453: dt=9.072000, rms=0.771 (0.005%), neg=0, invalid=1511
0454: dt=25.920000, rms=0.771 (0.005%), neg=0, invalid=1511
0455: dt=12.960000, rms=0.771 (0.026%), neg=0, invalid=1511
0456: dt=12.960000, rms=0.771 (0.008%), neg=0, invalid=1511
0457: dt=12.960000, rms=0.771 (0.001%), neg=0, invalid=1511
0458: dt=25.920000, rms=0.771 (0.021%), neg=0, invalid=1511
0459: dt=9.072000, rms=0.771 (0.012%), neg=0, invalid=1511
0460: dt=9.072000, rms=0.771 (-0.004%), neg=0, invalid=1511
setting smoothness coefficient to 0.12
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.793, neg=0, invalid=1511
0461: dt=2.800000, rms=0.791 (0.173%), neg=0, invalid=1511
0462: dt=0.700000, rms=0.791 (0.045%), neg=0, invalid=1511
0463: dt=0.700000, rms=0.790 (0.058%), neg=0, invalid=1511
0464: dt=0.700000, rms=0.790 (0.106%), neg=0, invalid=1511
0465: dt=0.700000, rms=0.788 (0.136%), neg=0, invalid=1511
0466: dt=0.700000, rms=0.787 (0.147%), neg=0, invalid=1511
0467: dt=0.700000, rms=0.786 (0.154%), neg=0, invalid=1511
0468: dt=0.700000, rms=0.785 (0.159%), neg=0, invalid=1511
0469: dt=0.700000, rms=0.784 (0.155%), neg=0, invalid=1511
0470: dt=0.700000, rms=0.783 (0.139%), neg=0, invalid=1511
0471: dt=0.700000, rms=0.782 (0.116%), neg=0, invalid=1511
0472: dt=0.700000, rms=0.781 (0.089%), neg=0, invalid=1511
0473: dt=0.700000, rms=0.780 (0.079%), neg=0, invalid=1511
0474: dt=0.700000, rms=0.780 (0.023%), neg=0, invalid=1511
0475: dt=0.700000, rms=0.780 (0.017%), neg=0, invalid=1511
0476: dt=2.800000, rms=0.780 (0.008%), neg=0, invalid=1511
0477: dt=0.700000, rms=0.780 (0.045%), neg=0, invalid=1511
0478: dt=0.700000, rms=0.779 (0.021%), neg=0, invalid=1511
0479: dt=0.700000, rms=0.779 (-0.018%), neg=0, invalid=1511
0480: dt=0.000000, rms=0.779 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.779, neg=0, invalid=1511
0481: dt=0.500000, rms=0.779 (0.008%), neg=0, invalid=1511
0482: dt=2.400000, rms=0.778 (0.205%), neg=0, invalid=1511
0483: dt=2.000000, rms=0.777 (0.139%), neg=0, invalid=1511
0484: dt=0.700000, rms=0.776 (0.026%), neg=0, invalid=1511
0485: dt=0.700000, rms=0.776 (0.074%), neg=0, invalid=1511
0486: dt=0.700000, rms=0.775 (0.122%), neg=0, invalid=1511
0487: dt=0.700000, rms=0.774 (0.143%), neg=0, invalid=1511
0488: dt=0.700000, rms=0.773 (0.150%), neg=0, invalid=1511
0489: dt=0.700000, rms=0.771 (0.161%), neg=0, invalid=1511
0490: dt=0.700000, rms=0.770 (0.167%), neg=0, invalid=1511
0491: dt=0.700000, rms=0.769 (0.179%), neg=0, invalid=1511
0492: dt=0.700000, rms=0.767 (0.183%), neg=0, invalid=1511
0493: dt=0.700000, rms=0.766 (0.185%), neg=0, invalid=1511
0494: dt=0.700000, rms=0.765 (0.183%), neg=0, invalid=1511
0495: dt=0.700000, rms=0.763 (0.178%), neg=0, invalid=1511
0496: dt=0.700000, rms=0.762 (0.174%), neg=0, invalid=1511
0497: dt=0.700000, rms=0.761 (0.161%), neg=0, invalid=1511
0498: dt=0.700000, rms=0.759 (0.160%), neg=0, invalid=1511
0499: dt=0.700000, rms=0.758 (0.164%), neg=0, invalid=1511
0500: dt=0.700000, rms=0.757 (0.169%), neg=0, invalid=1511
0501: dt=0.700000, rms=0.756 (0.166%), neg=0, invalid=1511
0502: dt=0.700000, rms=0.754 (0.158%), neg=0, invalid=1511
0503: dt=0.700000, rms=0.753 (0.148%), neg=0, invalid=1511
0504: dt=0.700000, rms=0.752 (0.136%), neg=0, invalid=1511
0505: dt=0.700000, rms=0.751 (0.136%), neg=0, invalid=1511
0506: dt=0.700000, rms=0.750 (0.134%), neg=0, invalid=1511
0507: dt=0.700000, rms=0.749 (0.133%), neg=0, invalid=1511
0508: dt=0.700000, rms=0.748 (0.134%), neg=0, invalid=1511
0509: dt=0.700000, rms=0.747 (0.129%), neg=0, invalid=1511
0510: dt=0.700000, rms=0.746 (0.121%), neg=0, invalid=1511
0511: dt=0.700000, rms=0.746 (0.111%), neg=0, invalid=1511
0512: dt=0.700000, rms=0.745 (0.108%), neg=0, invalid=1511
0513: dt=0.700000, rms=0.744 (0.108%), neg=0, invalid=1511
0514: dt=0.700000, rms=0.743 (0.100%), neg=0, invalid=1511
0515: dt=0.700000, rms=0.742 (0.099%), neg=0, invalid=1511
0516: dt=0.700000, rms=0.742 (0.096%), neg=0, invalid=1511
0517: dt=0.700000, rms=0.741 (0.092%), neg=0, invalid=1511
0518: dt=0.700000, rms=0.740 (0.092%), neg=0, invalid=1511
0519: dt=0.700000, rms=0.740 (0.086%), neg=0, invalid=1511
0520: dt=0.700000, rms=0.739 (0.081%), neg=0, invalid=1511
0521: dt=0.700000, rms=0.739 (0.006%), neg=0, invalid=1511
0522: dt=0.700000, rms=0.739 (0.012%), neg=0, invalid=1511
0523: dt=0.700000, rms=0.739 (0.017%), neg=0, invalid=1511
0524: dt=0.700000, rms=0.739 (0.024%), neg=0, invalid=1511
0525: dt=0.700000, rms=0.738 (0.029%), neg=0, invalid=1511
0526: dt=0.700000, rms=0.738 (0.034%), neg=0, invalid=1511
0527: dt=0.700000, rms=0.738 (0.037%), neg=0, invalid=1511
0528: dt=0.700000, rms=0.738 (0.040%), neg=0, invalid=1511
0529: dt=0.700000, rms=0.737 (0.047%), neg=0, invalid=1511
0530: dt=0.700000, rms=0.737 (0.044%), neg=0, invalid=1511
0531: dt=0.700000, rms=0.737 (0.046%), neg=0, invalid=1511
0532: dt=0.700000, rms=0.736 (0.048%), neg=0, invalid=1511
0533: dt=0.700000, rms=0.736 (0.049%), neg=0, invalid=1511
0534: dt=0.700000, rms=0.736 (0.048%), neg=0, invalid=1511
0535: dt=0.700000, rms=0.735 (0.050%), neg=0, invalid=1511
0536: dt=0.700000, rms=0.735 (0.045%), neg=0, invalid=1511
0537: dt=0.700000, rms=0.735 (0.051%), neg=0, invalid=1511
0538: dt=0.700000, rms=0.734 (0.048%), neg=0, invalid=1511
0539: dt=0.700000, rms=0.734 (0.043%), neg=0, invalid=1511
0540: dt=0.700000, rms=0.734 (0.043%), neg=0, invalid=1511
0541: dt=0.700000, rms=0.734 (0.003%), neg=0, invalid=1511
0542: dt=0.700000, rms=0.733 (0.003%), neg=0, invalid=1511
0543: dt=0.700000, rms=0.733 (0.005%), neg=0, invalid=1511
0544: dt=3.333333, rms=0.733 (0.017%), neg=0, invalid=1511
0545: dt=2.800000, rms=0.733 (0.022%), neg=0, invalid=1511
0546: dt=2.400000, rms=0.733 (0.019%), neg=0, invalid=1511
0547: dt=2.400000, rms=0.733 (0.016%), neg=0, invalid=1511
0548: dt=2.400000, rms=0.733 (0.036%), neg=0, invalid=1511
0549: dt=2.400000, rms=0.732 (0.056%), neg=0, invalid=1511
0550: dt=2.400000, rms=0.732 (0.070%), neg=0, invalid=1511
0551: dt=2.400000, rms=0.732 (0.016%), neg=0, invalid=1511
0552: dt=2.400000, rms=0.732 (0.011%), neg=0, invalid=1511
0553: dt=0.400000, rms=0.732 (0.002%), neg=0, invalid=1511
0554: dt=0.400000, rms=0.731 (0.003%), neg=0, invalid=1511
0555: dt=0.400000, rms=0.731 (0.006%), neg=0, invalid=1511
0556: dt=0.400000, rms=0.731 (0.010%), neg=0, invalid=1511
0557: dt=0.400000, rms=0.731 (0.012%), neg=0, invalid=1511
0558: dt=0.400000, rms=0.731 (0.014%), neg=0, invalid=1511
0559: dt=0.400000, rms=0.731 (0.017%), neg=0, invalid=1511
0560: dt=0.400000, rms=0.731 (0.016%), neg=0, invalid=1511
0561: dt=0.400000, rms=0.731 (0.015%), neg=0, invalid=1511
setting smoothness coefficient to 0.40
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.791, neg=0, invalid=1511
0562: dt=0.000000, rms=0.791 (-0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.791, neg=0, invalid=1511
0563: dt=0.000000, rms=0.791 (0.000%), neg=0, invalid=1511
setting smoothness coefficient to 1.00
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.905, neg=0, invalid=1511
0564: dt=0.028000, rms=0.905 (0.015%), neg=0, invalid=1511
0565: dt=0.028000, rms=0.905 (0.016%), neg=0, invalid=1511
0566: dt=0.028000, rms=0.905 (0.017%), neg=0, invalid=1511
0567: dt=0.031250, rms=0.905 (0.019%), neg=0, invalid=1511
0568: dt=0.031250, rms=0.904 (0.017%), neg=0, invalid=1511
0569: dt=0.031250, rms=0.904 (0.019%), neg=0, invalid=1511
0570: dt=0.031250, rms=0.904 (0.034%), neg=0, invalid=1511
0571: dt=0.031250, rms=0.904 (0.048%), neg=0, invalid=1511
0572: dt=0.031250, rms=0.903 (0.058%), neg=0, invalid=1511
0573: dt=0.031250, rms=0.902 (0.065%), neg=0, invalid=1511
0574: dt=0.031250, rms=0.902 (0.015%), neg=0, invalid=1511
0575: dt=0.031250, rms=0.902 (0.028%), neg=0, invalid=1511
0576: dt=0.031250, rms=0.902 (0.039%), neg=0, invalid=1511
0577: dt=0.031250, rms=0.901 (0.048%), neg=0, invalid=1511
0578: dt=0.031250, rms=0.901 (0.050%), neg=0, invalid=1511
0579: dt=0.031250, rms=0.900 (0.055%), neg=0, invalid=1511
0580: dt=0.031250, rms=0.900 (0.056%), neg=0, invalid=1511
0581: dt=0.031250, rms=0.899 (0.055%), neg=0, invalid=1511
0582: dt=0.031250, rms=0.899 (0.052%), neg=0, invalid=1511
0583: dt=0.031250, rms=0.898 (0.050%), neg=0, invalid=1511
0584: dt=0.031250, rms=0.898 (0.044%), neg=0, invalid=1511
0585: dt=0.031250, rms=0.898 (0.039%), neg=0, invalid=1511
0586: dt=0.031250, rms=0.897 (0.033%), neg=0, invalid=1511
0587: dt=0.031250, rms=0.897 (0.023%), neg=0, invalid=1511
0588: dt=0.031250, rms=0.897 (0.018%), neg=0, invalid=1511
0589: dt=0.000000, rms=0.897 (0.000%), neg=0, invalid=1511
0590: dt=0.100000, rms=0.897 (-0.003%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.897, neg=0, invalid=1511
0591: dt=0.448000, rms=0.895 (0.257%), neg=0, invalid=1511
0592: dt=0.062500, rms=0.894 (0.031%), neg=0, invalid=1511
0593: dt=0.031250, rms=0.894 (0.015%), neg=0, invalid=1511
0594: dt=0.031250, rms=0.894 (0.014%), neg=0, invalid=1511
0595: dt=0.031250, rms=0.894 (0.014%), neg=0, invalid=1511
0596: dt=0.031250, rms=0.894 (0.027%), neg=0, invalid=1511
0597: dt=0.031250, rms=0.893 (0.036%), neg=0, invalid=1511
0598: dt=0.031250, rms=0.893 (0.045%), neg=0, invalid=1511
0599: dt=0.001953, rms=0.893 (0.001%), neg=0, invalid=1511
0600: dt=0.001953, rms=0.893 (0.001%), neg=0, invalid=1511
0601: dt=0.000977, rms=0.893 (0.000%), neg=0, invalid=1511
0602: dt=0.000438, rms=0.893 (0.000%), neg=0, invalid=1511
0603: dt=0.000219, rms=0.893 (0.000%), neg=0, invalid=1511
0604: dt=0.000219, rms=0.893 (-0.000%), neg=0, invalid=1511
resetting metric properties...
setting smoothness coefficient to 2.00
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.715, neg=0, invalid=1511
0605: dt=0.000078, rms=0.715 (0.001%), neg=0, invalid=1511
0606: dt=0.007000, rms=0.714 (0.070%), neg=0, invalid=1511
0607: dt=0.028000, rms=0.712 (0.276%), neg=0, invalid=1511
0608: dt=0.028000, rms=0.711 (0.265%), neg=0, invalid=1511
0609: dt=0.028000, rms=0.709 (0.254%), neg=0, invalid=1511
0610: dt=0.024000, rms=0.707 (0.206%), neg=0, invalid=1511
0611: dt=0.006000, rms=0.707 (0.050%), neg=0, invalid=1511
0612: dt=0.001750, rms=0.707 (0.014%), neg=0, invalid=1511
0613: dt=0.001750, rms=0.707 (0.015%), neg=0, invalid=1511
0614: dt=0.001750, rms=0.707 (0.027%), neg=0, invalid=1511
0615: dt=0.001750, rms=0.706 (0.014%), neg=0, invalid=1511
0616: dt=0.000219, rms=0.706 (0.002%), neg=0, invalid=1511
0617: dt=0.000109, rms=0.706 (0.001%), neg=0, invalid=1511
0618: dt=0.000055, rms=0.706 (0.000%), neg=0, invalid=1511
0619: dt=0.000027, rms=0.706 (0.000%), neg=0, invalid=1511
0620: dt=0.000027, rms=0.706 (0.000%), neg=0, invalid=1511
0621: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0622: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0623: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0624: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0625: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0626: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0627: dt=0.000027, rms=0.706 (0.001%), neg=0, invalid=1511
0628: dt=0.000027, rms=0.706 (0.002%), neg=0, invalid=1511
0629: dt=0.000027, rms=0.706 (0.002%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.706, neg=0, invalid=1511
0630: dt=0.020000, rms=0.705 (0.253%), neg=0, invalid=1511
0631: dt=0.001750, rms=0.704 (0.021%), neg=0, invalid=1511
0632: dt=0.000219, rms=0.704 (0.003%), neg=0, invalid=1511
0633: dt=0.000055, rms=0.704 (0.001%), neg=0, invalid=1511
0634: dt=0.000055, rms=0.704 (0.001%), neg=0, invalid=1511
0635: dt=0.000027, rms=0.704 (0.000%), neg=0, invalid=1511
label assignment complete, 0 changed (0.00%)
***************** morphing with label term set to 0 *******************************
label assignment complete, 0 changed (0.00%)
**************** pass 1 of 1 ************************
setting smoothness coefficient to 0.01
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0636: dt=0.000000, rms=0.704 (-0.046%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0637: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
setting smoothness coefficient to 0.03
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0638: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0639: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
setting smoothness coefficient to 0.12
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0640: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0641: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
setting smoothness coefficient to 0.40
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0642: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.704, neg=0, invalid=1511
0643: dt=0.000000, rms=0.704 (0.000%), neg=0, invalid=1511
setting smoothness coefficient to 1.00
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.705, neg=0, invalid=1511
0644: dt=0.112000, rms=0.701 (0.488%), neg=0, invalid=1511
0645: dt=0.028000, rms=0.700 (0.115%), neg=0, invalid=1511
0646: dt=0.020000, rms=0.700 (0.082%), neg=0, invalid=1511
0647: dt=0.080000, rms=0.697 (0.322%), neg=0, invalid=1511
0648: dt=0.007000, rms=0.697 (0.027%), neg=0, invalid=1511
0649: dt=0.007000, rms=0.697 (0.027%), neg=0, invalid=1511
0650: dt=0.007000, rms=0.697 (0.051%), neg=0, invalid=1511
0651: dt=0.007000, rms=0.696 (0.073%), neg=0, invalid=1511
0652: dt=0.007000, rms=0.696 (0.027%), neg=0, invalid=1511
0653: dt=0.007000, rms=0.696 (0.050%), neg=0, invalid=1511
0654: dt=0.007000, rms=0.696 (0.026%), neg=0, invalid=1511
0655: dt=0.007000, rms=0.695 (0.050%), neg=0, invalid=1511
0656: dt=0.007000, rms=0.695 (0.071%), neg=0, invalid=1511
0657: dt=0.007000, rms=0.694 (0.089%), neg=0, invalid=1511
0658: dt=0.007000, rms=0.694 (0.025%), neg=0, invalid=1511
0659: dt=0.007000, rms=0.694 (0.048%), neg=0, invalid=1511
0660: dt=0.007000, rms=0.693 (0.068%), neg=0, invalid=1511
0661: dt=0.007000, rms=0.692 (0.086%), neg=0, invalid=1511
0662: dt=0.007000, rms=0.692 (0.101%), neg=0, invalid=1511
0663: dt=0.007000, rms=0.692 (0.024%), neg=0, invalid=1511
0664: dt=0.007000, rms=0.691 (0.045%), neg=0, invalid=1511
0665: dt=0.007000, rms=0.691 (0.065%), neg=0, invalid=1511
0666: dt=0.007000, rms=0.690 (0.082%), neg=0, invalid=1511
0667: dt=0.007000, rms=0.690 (0.023%), neg=0, invalid=1511
0668: dt=0.007000, rms=0.690 (0.023%), neg=0, invalid=1511
0669: dt=0.001750, rms=0.690 (0.006%), neg=0, invalid=1511
0670: dt=0.000109, rms=0.690 (0.000%), neg=0, invalid=1511
0671: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0672: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0673: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0674: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0675: dt=0.000109, rms=0.690 (0.002%), neg=0, invalid=1511
0676: dt=0.000109, rms=0.690 (0.002%), neg=0, invalid=1511
0677: dt=0.000109, rms=0.690 (0.000%), neg=0, invalid=1511
0678: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0679: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0680: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0681: dt=0.000109, rms=0.690 (0.001%), neg=0, invalid=1511
0682: dt=0.000109, rms=0.690 (0.002%), neg=0, invalid=1511
0683: dt=0.000109, rms=0.690 (0.002%), neg=0, invalid=1511
0684: dt=0.000109, rms=0.690 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.690, neg=0, invalid=1511
0685: dt=0.112000, rms=0.687 (0.456%), neg=0, invalid=1511
0686: dt=0.028000, rms=0.686 (0.108%), neg=0, invalid=1511
0687: dt=0.024000, rms=0.685 (0.091%), neg=0, invalid=1511
0688: dt=0.024000, rms=0.685 (0.090%), neg=0, invalid=1511
0689: dt=0.007000, rms=0.684 (0.027%), neg=0, invalid=1511
0690: dt=0.007000, rms=0.684 (0.026%), neg=0, invalid=1511
0691: dt=0.007000, rms=0.684 (0.050%), neg=0, invalid=1511
0692: dt=0.007000, rms=0.683 (0.071%), neg=0, invalid=1511
0693: dt=0.007000, rms=0.683 (0.089%), neg=0, invalid=1511
0694: dt=0.007000, rms=0.682 (0.106%), neg=0, invalid=1511
0695: dt=0.007000, rms=0.681 (0.121%), neg=0, invalid=1511
0696: dt=0.007000, rms=0.680 (0.133%), neg=0, invalid=1511
0697: dt=0.007000, rms=0.679 (0.144%), neg=0, invalid=1511
0698: dt=0.007000, rms=0.678 (0.154%), neg=0, invalid=1511
0699: dt=0.007000, rms=0.677 (0.161%), neg=0, invalid=1511
0700: dt=0.007000, rms=0.676 (0.168%), neg=0, invalid=1511
0701: dt=0.007000, rms=0.675 (0.174%), neg=0, invalid=1511
0702: dt=0.007000, rms=0.674 (0.178%), neg=0, invalid=1511
0703: dt=0.007000, rms=0.673 (0.181%), neg=0, invalid=1511
0704: dt=0.007000, rms=0.671 (0.183%), neg=0, invalid=1511
0705: dt=0.007000, rms=0.670 (0.185%), neg=0, invalid=1511
0706: dt=0.007000, rms=0.669 (0.185%), neg=0, invalid=1511
0707: dt=0.007000, rms=0.668 (0.185%), neg=0, invalid=1511
0708: dt=0.007000, rms=0.667 (0.018%), neg=0, invalid=1511
0709: dt=0.007000, rms=0.667 (0.035%), neg=0, invalid=1511
0710: dt=0.007000, rms=0.667 (0.050%), neg=0, invalid=1511
0711: dt=0.007000, rms=0.666 (0.064%), neg=0, invalid=1511
0712: dt=0.007000, rms=0.666 (0.076%), neg=0, invalid=1511
0713: dt=0.007000, rms=0.666 (0.018%), neg=0, invalid=1511
0714: dt=0.007000, rms=0.666 (0.035%), neg=0, invalid=1511
0715: dt=0.001750, rms=0.666 (0.005%), neg=0, invalid=1511
0716: dt=0.000219, rms=0.666 (0.001%), neg=0, invalid=1511
0717: dt=0.028000, rms=0.665 (0.073%), neg=0, invalid=1511
0718: dt=0.028000, rms=0.665 (0.072%), neg=0, invalid=1511
0719: dt=0.024000, rms=0.664 (0.061%), neg=0, invalid=1511
0720: dt=0.024000, rms=0.664 (0.060%), neg=0, invalid=1511
0721: dt=0.007812, rms=0.664 (0.020%), neg=0, invalid=1511
0722: dt=0.001953, rms=0.664 (0.005%), neg=0, invalid=1511
0723: dt=0.000244, rms=0.664 (0.001%), neg=0, invalid=1511
0724: dt=0.000244, rms=0.664 (0.001%), neg=0, invalid=1511
0725: dt=0.000244, rms=0.664 (0.002%), neg=0, invalid=1511
0726: dt=0.000244, rms=0.664 (0.002%), neg=0, invalid=1511
0727: dt=0.000244, rms=0.664 (0.003%), neg=0, invalid=1511
0728: dt=0.000244, rms=0.664 (0.003%), neg=0, invalid=1511
0729: dt=0.000244, rms=0.664 (0.003%), neg=0, invalid=1511
0730: dt=0.000244, rms=0.664 (0.004%), neg=0, invalid=1511
0731: dt=0.000244, rms=0.664 (0.004%), neg=0, invalid=1511
0732: dt=0.000122, rms=0.664 (0.000%), neg=0, invalid=1511
0733: dt=0.028000, rms=0.663 (0.071%), neg=0, invalid=1511
0734: dt=0.028000, rms=0.663 (0.069%), neg=0, invalid=1511
0735: dt=0.028000, rms=0.662 (0.069%), neg=0, invalid=1511
0736: dt=0.007000, rms=0.662 (0.016%), neg=0, invalid=1511
0737: dt=0.007000, rms=0.662 (0.017%), neg=0, invalid=1511
0738: dt=0.007000, rms=0.662 (0.032%), neg=0, invalid=1511
0739: dt=0.007000, rms=0.661 (0.046%), neg=0, invalid=1511
0740: dt=0.007000, rms=0.661 (0.059%), neg=0, invalid=1511
0741: dt=0.007000, rms=0.661 (0.069%), neg=0, invalid=1511
0742: dt=0.007000, rms=0.660 (0.078%), neg=0, invalid=1511
0743: dt=0.007000, rms=0.659 (0.087%), neg=0, invalid=1511
0744: dt=0.007000, rms=0.659 (0.094%), neg=0, invalid=1511
0745: dt=0.007000, rms=0.658 (0.100%), neg=0, invalid=1511
0746: dt=0.007000, rms=0.657 (0.106%), neg=0, invalid=1511
0747: dt=0.007000, rms=0.657 (0.015%), neg=0, invalid=1511
0748: dt=0.007000, rms=0.657 (0.029%), neg=0, invalid=1511
0749: dt=0.007000, rms=0.657 (0.041%), neg=0, invalid=1511
0750: dt=0.007000, rms=0.657 (0.052%), neg=0, invalid=1511
0751: dt=0.007000, rms=0.656 (0.062%), neg=0, invalid=1511
0752: dt=0.007000, rms=0.656 (0.070%), neg=0, invalid=1511
0753: dt=0.001750, rms=0.656 (0.004%), neg=0, invalid=1511
0754: dt=0.001750, rms=0.656 (0.007%), neg=0, invalid=1511
0755: dt=0.001750, rms=0.656 (0.010%), neg=0, invalid=1511
0756: dt=0.001750, rms=0.655 (0.013%), neg=0, invalid=1511
0757: dt=0.001750, rms=0.655 (0.015%), neg=0, invalid=1511
0758: dt=0.001750, rms=0.655 (0.017%), neg=0, invalid=1511
0759: dt=0.001750, rms=0.655 (0.019%), neg=0, invalid=1511
0760: dt=0.001750, rms=0.655 (0.021%), neg=0, invalid=1511
0761: dt=0.001750, rms=0.655 (0.022%), neg=0, invalid=1511
0762: dt=0.001750, rms=0.655 (0.024%), neg=0, invalid=1511
0763: dt=0.001750, rms=0.655 (0.025%), neg=0, invalid=1511
0764: dt=0.001750, rms=0.654 (0.026%), neg=0, invalid=1511
0765: dt=0.001750, rms=0.654 (0.027%), neg=0, invalid=1511
0766: dt=0.001750, rms=0.654 (0.028%), neg=0, invalid=1511
0767: dt=0.001750, rms=0.654 (0.028%), neg=0, invalid=1511
0768: dt=0.001750, rms=0.654 (0.029%), neg=0, invalid=1511
0769: dt=0.001750, rms=0.653 (0.030%), neg=0, invalid=1511
0770: dt=0.001750, rms=0.653 (0.030%), neg=0, invalid=1511
0771: dt=0.001750, rms=0.653 (0.030%), neg=0, invalid=1511
0772: dt=0.001750, rms=0.653 (0.003%), neg=0, invalid=1511
0773: dt=0.001750, rms=0.653 (0.006%), neg=0, invalid=1511
0774: dt=0.001750, rms=0.653 (0.009%), neg=0, invalid=1511
0775: dt=0.001750, rms=0.653 (0.012%), neg=0, invalid=1511
0776: dt=0.000875, rms=0.653 (0.002%), neg=0, invalid=1511
0777: dt=0.320000, rms=0.649 (0.617%), neg=0, invalid=1511
0778: dt=0.024000, rms=0.649 (0.041%), neg=0, invalid=1511
0779: dt=0.012000, rms=0.648 (0.021%), neg=0, invalid=1511
0780: dt=0.012000, rms=0.648 (0.039%), neg=0, invalid=1511
0781: dt=0.006000, rms=0.648 (0.010%), neg=0, invalid=1511
0782: dt=0.006000, rms=0.648 (0.019%), neg=0, invalid=1511
0783: dt=0.006000, rms=0.648 (0.010%), neg=0, invalid=1511
0784: dt=0.011719, rms=0.648 (0.019%), neg=0, invalid=1511
0785: dt=0.007000, rms=0.648 (0.011%), neg=0, invalid=1511
0786: dt=0.007000, rms=0.648 (0.012%), neg=0, invalid=1511
0787: dt=0.007000, rms=0.647 (0.023%), neg=0, invalid=1511
0788: dt=0.007000, rms=0.647 (0.032%), neg=0, invalid=1511
0789: dt=0.007000, rms=0.647 (0.040%), neg=0, invalid=1511
0790: dt=0.007000, rms=0.647 (0.047%), neg=0, invalid=1511
0791: dt=0.007000, rms=0.646 (0.054%), neg=0, invalid=1511
0792: dt=0.007000, rms=0.646 (0.060%), neg=0, invalid=1511
0793: dt=0.007000, rms=0.646 (0.065%), neg=0, invalid=1511
0794: dt=0.007000, rms=0.645 (0.070%), neg=0, invalid=1511
0795: dt=0.007000, rms=0.645 (0.074%), neg=0, invalid=1511
0796: dt=0.007000, rms=0.644 (0.077%), neg=0, invalid=1511
0797: dt=0.007000, rms=0.644 (0.011%), neg=0, invalid=1511
0798: dt=0.007000, rms=0.644 (0.020%), neg=0, invalid=1511
0799: dt=0.007000, rms=0.644 (0.029%), neg=0, invalid=1511
0800: dt=0.007000, rms=0.644 (0.036%), neg=0, invalid=1511
0801: dt=0.007000, rms=0.643 (0.043%), neg=0, invalid=1511
0802: dt=0.007000, rms=0.643 (0.049%), neg=0, invalid=1511
0803: dt=0.007000, rms=0.643 (0.010%), neg=0, invalid=1511
0804: dt=0.007000, rms=0.643 (0.019%), neg=0, invalid=1511
0805: dt=0.007000, rms=0.643 (0.028%), neg=0, invalid=1511
0806: dt=0.007000, rms=0.642 (0.035%), neg=0, invalid=1511
0807: dt=0.007000, rms=0.642 (0.041%), neg=0, invalid=1511
0808: dt=0.007000, rms=0.642 (0.047%), neg=0, invalid=1511
0809: dt=0.007000, rms=0.641 (0.052%), neg=0, invalid=1511
0810: dt=0.007000, rms=0.641 (0.010%), neg=0, invalid=1511
0811: dt=0.000875, rms=0.641 (0.001%), neg=0, invalid=1511
0812: dt=0.000438, rms=0.641 (0.001%), neg=0, invalid=1511
0813: dt=0.000219, rms=0.641 (0.000%), neg=0, invalid=1511
0814: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0815: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0816: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0817: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0818: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0819: dt=0.000219, rms=0.641 (0.002%), neg=0, invalid=1511
0820: dt=0.000219, rms=0.641 (0.002%), neg=0, invalid=1511
0821: dt=0.000219, rms=0.641 (0.000%), neg=0, invalid=1511
0822: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0823: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0824: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0825: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0826: dt=0.000219, rms=0.641 (0.001%), neg=0, invalid=1511
0827: dt=0.000219, rms=0.641 (0.002%), neg=0, invalid=1511
0828: dt=0.000219, rms=0.641 (0.002%), neg=0, invalid=1511
0829: dt=0.000219, rms=0.641 (0.000%), neg=0, invalid=1511
resetting metric properties...
setting smoothness coefficient to 2.00
blurring input image with Gaussian with sigma=2.000...
0000: dt=0.000, rms=0.639, neg=0, invalid=1511
0830: dt=0.000027, rms=0.639 (0.000%), neg=0, invalid=1511
0831: dt=0.007000, rms=0.639 (0.058%), neg=0, invalid=1511
0832: dt=0.007000, rms=0.639 (0.057%), neg=0, invalid=1511
0833: dt=0.007000, rms=0.638 (0.056%), neg=0, invalid=1511
0834: dt=0.001750, rms=0.638 (0.014%), neg=0, invalid=1511
0835: dt=0.001750, rms=0.638 (0.014%), neg=0, invalid=1511
0836: dt=0.001750, rms=0.638 (0.026%), neg=0, invalid=1511
0837: dt=0.000875, rms=0.638 (0.007%), neg=0, invalid=1511
0838: dt=0.000875, rms=0.638 (0.007%), neg=0, invalid=1511
0839: dt=0.000875, rms=0.638 (0.013%), neg=0, invalid=1511
0840: dt=0.000875, rms=0.638 (0.018%), neg=0, invalid=1511
0841: dt=0.000875, rms=0.638 (0.023%), neg=0, invalid=1511
0842: dt=0.000875, rms=0.637 (0.028%), neg=0, invalid=1511
0843: dt=0.000875, rms=0.637 (0.032%), neg=0, invalid=1511
0844: dt=0.000875, rms=0.637 (0.035%), neg=0, invalid=1511
0845: dt=0.000875, rms=0.637 (0.007%), neg=0, invalid=1511
0846: dt=0.000875, rms=0.637 (0.012%), neg=0, invalid=1511
0847: dt=0.000438, rms=0.637 (0.003%), neg=0, invalid=1511
0848: dt=0.005000, rms=0.637 (0.037%), neg=0, invalid=1511
0849: dt=0.001750, rms=0.637 (0.013%), neg=0, invalid=1511
0850: dt=0.000875, rms=0.637 (0.006%), neg=0, invalid=1511
0851: dt=0.000875, rms=0.636 (0.006%), neg=0, invalid=1511
0852: dt=0.000875, rms=0.636 (0.006%), neg=0, invalid=1511
0853: dt=0.000875, rms=0.636 (0.012%), neg=0, invalid=1511
0854: dt=0.000875, rms=0.636 (0.017%), neg=0, invalid=1511
0855: dt=0.000875, rms=0.636 (0.022%), neg=0, invalid=1511
0856: dt=0.000875, rms=0.636 (0.026%), neg=0, invalid=1511
0857: dt=0.000875, rms=0.636 (0.006%), neg=0, invalid=1511
0858: dt=0.000875, rms=0.636 (0.006%), neg=0, invalid=1511
0859: dt=0.000875, rms=0.636 (0.006%), neg=0, invalid=1511
0860: dt=0.000438, rms=0.636 (0.003%), neg=0, invalid=1511
0861: dt=0.000313, rms=0.636 (0.002%), neg=0, invalid=1511
0862: dt=0.000039, rms=0.636 (0.000%), neg=0, invalid=1511
0863: dt=0.000020, rms=0.636 (0.000%), neg=0, invalid=1511
0864: dt=0.000020, rms=0.636 (0.000%), neg=0, invalid=1511
0865: dt=0.000020, rms=0.636 (0.000%), neg=0, invalid=1511
0866: dt=0.000020, rms=0.636 (0.000%), neg=0, invalid=1511
0867: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0868: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0869: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0870: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0871: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0872: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0873: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0874: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0875: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0876: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0877: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0878: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0879: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0880: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0881: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0882: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0883: dt=0.000020, rms=0.636 (0.001%), neg=0, invalid=1511
0884: dt=0.000020, rms=0.636 (0.000%), neg=0, invalid=1511
blurring input image with Gaussian with sigma=0.500...
0000: dt=0.000, rms=0.636, neg=0, invalid=1511
0885: dt=0.001750, rms=0.636 (0.020%), neg=0, invalid=1511
0886: dt=0.006000, rms=0.635 (0.070%), neg=0, invalid=1511
0887: dt=0.001500, rms=0.635 (0.017%), neg=0, invalid=1511
0888: dt=0.000094, rms=0.635 (0.001%), neg=0, invalid=1511
0889: dt=0.000023, rms=0.635 (0.000%), neg=0, invalid=1511
0890: dt=0.000023, rms=0.635 (0.001%), neg=0, invalid=1511
0891: dt=0.000023, rms=0.635 (0.001%), neg=0, invalid=1511
0892: dt=0.000023, rms=0.635 (0.001%), neg=0, invalid=1511
0893: dt=0.000023, rms=0.635 (0.001%), neg=0, invalid=1511
0894: dt=0.000023, rms=0.635 (0.001%), neg=0, invalid=1511
0895: dt=0.000023, rms=0.635 (0.001%), neg=0, invalid=1511
0896: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0897: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0898: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0899: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0900: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0901: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0902: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0903: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0904: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0905: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0906: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0907: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0908: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0909: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0910: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0911: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0912: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0913: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0914: dt=0.000023, rms=0.635 (0.002%), neg=0, invalid=1511
0915: dt=0.000023, rms=0.635 (0.003%), neg=0, invalid=1511
0916: dt=0.000023, rms=0.635 (0.000%), neg=0, invalid=1511
0917: dt=0.000050, rms=0.635 (0.001%), neg=0, invalid=1511
0918: dt=0.007000, rms=0.634 (0.079%), neg=0, invalid=1511
0919: dt=0.000438, rms=0.634 (0.005%), neg=0, invalid=1511
0920: dt=0.000438, rms=0.634 (0.005%), neg=0, invalid=1511
0921: dt=0.000438, rms=0.634 (0.009%), neg=0, invalid=1511
0922: dt=0.000438, rms=0.634 (0.013%), neg=0, invalid=1511
0923: dt=0.000438, rms=0.634 (0.017%), neg=0, invalid=1511
0924: dt=0.000438, rms=0.634 (0.020%), neg=0, invalid=1511
0925: dt=0.000438, rms=0.634 (0.023%), neg=0, invalid=1511
0926: dt=0.000438, rms=0.633 (0.025%), neg=0, invalid=1511
0927: dt=0.000438, rms=0.633 (0.027%), neg=0, invalid=1511
0928: dt=0.000438, rms=0.633 (0.029%), neg=0, invalid=1511
0929: dt=0.000438, rms=0.633 (0.031%), neg=0, invalid=1511
0930: dt=0.000438, rms=0.633 (0.033%), neg=0, invalid=1511
0931: dt=0.000438, rms=0.632 (0.034%), neg=0, invalid=1511
0932: dt=0.000438, rms=0.632 (0.035%), neg=0, invalid=1511
0933: dt=0.000438, rms=0.632 (0.036%), neg=0, invalid=1511
0934: dt=0.000438, rms=0.632 (0.005%), neg=0, invalid=1511
0935: dt=0.000438, rms=0.632 (0.009%), neg=0, invalid=1511
0936: dt=0.000438, rms=0.632 (0.012%), neg=0, invalid=1511
0937: dt=0.000219, rms=0.632 (0.002%), neg=0, invalid=1511
0938: dt=0.007000, rms=0.631 (0.072%), neg=0, invalid=1511
0939: dt=0.000438, rms=0.631 (0.004%), neg=0, invalid=1511
0940: dt=0.000438, rms=0.631 (0.004%), neg=0, invalid=1511
0941: dt=0.000438, rms=0.631 (0.008%), neg=0, invalid=1511
0942: dt=0.000438, rms=0.631 (0.012%), neg=0, invalid=1511
0943: dt=0.000438, rms=0.631 (0.004%), neg=0, invalid=1511
0944: dt=0.000438, rms=0.631 (0.008%), neg=0, invalid=1511
0945: dt=0.000438, rms=0.631 (0.012%), neg=0, invalid=1511
0946: dt=0.000438, rms=0.631 (0.015%), neg=0, invalid=1511
0947: dt=0.000438, rms=0.631 (0.004%), neg=0, invalid=1511
0948: dt=0.000438, rms=0.631 (0.004%), neg=0, invalid=1511
0949: dt=0.001500, rms=0.631 (0.015%), neg=0, invalid=1511
0950: dt=0.001500, rms=0.631 (0.015%), neg=0, invalid=1511
writing output transformation to transforms/talairach.m3z...

CA registration took 13 hours, 46 minutes and 11 seconds.
#--------------------------------------
#@# CA Reg Inv Thu Jun 21 00:22:52 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_ca_register -invert-and-save transforms/talairach.m3z 

Loading, Inverting, Saving, Exiting ...
Reading transforms/talairach.m3z 
gcam->type = vox
reading labels out of gcam file...
Inverting GCAM
Saving inverse 
#--------------------------------------
#@# Remove Neck Thu Jun 21 00:29:50 EDT 2007

 mri_remove_neck -radius 25 nu.mgz transforms/talairach.m3z /usr/local/freesurfer/average/RB_all_2006-02-15.gca nu_noneck.mgz 

erasing everything more than 25 mm from possible brain
reading atlas '/usr/local/freesurfer/average/RB_all_2006-02-15.gca'...
setting gca type = Normal gca type
reading input volume 'nu.mgz'...
reading transform 'transforms/talairach.m3z'...
gcam->type = vox
reading labels out of gcam file...
removing structures at least 25 mm from brain...
15239279 nonbrain voxels erased
writing output to nu_noneck.mgz...
nonbrain removal took 7 minutes and 16 seconds.
#--------------------------------------
#@# SkullLTA Thu Jun 21 00:37:06 EDT 2007

 mri_em_register -skull -t transforms/talairach.lta nu_noneck.mgz /usr/local/freesurfer/average/RB_all_withskull_2006-02-15.gca transforms/talairach_with_skull.lta 

aligning to atlas containing skull, setting unknown_nbr_spacing = 5
using previously computed transform transforms/talairach.lta
reading 1 input volumes...
logging results to talairach_with_skull.log
reading '/usr/local/freesurfer/average/RB_all_withskull_2006-02-15.gca'...
setting gca type = Normal gca type
average std = 23.5   using min determinant for regularization = 55.2
0 singular and 6473 ill-conditioned covariance matrices regularized
reading 'nu_noneck.mgz'...
total sample mean = 91.1 (1456 zeros)
register_mri: find_optimal_transform
find_optimal_transform: nsamples 3516, passno 0, spacing 8
Result so far: scale 1.000: max_log_p=-29219.3, old_max_log_p =-2022144.0 (thresh=-2020121.9)
 1.323   1.989   0.141  -485.651;
-1.171   1.193  -1.070   183.653;
-1.188   0.435   1.169  -21.439;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-26576.7, old_max_log_p =-29219.3 (thresh=-29190.1)
 0.699   2.727   0.174  -553.854;
-1.213   0.514  -1.011   324.110;
-1.800  -0.155   1.297   160.903;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-24002.4, old_max_log_p =-26576.7 (thresh=-26550.1)
-0.257   2.641   0.523  -469.196;
-1.086   0.345  -1.493   426.564;
-2.483  -1.400   1.151   529.027;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-23359.3, old_max_log_p =-24002.4 (thresh=-23978.4)
-0.289   2.971   0.589  -543.562;
-1.598   0.165  -1.483   534.785;
-2.759  -1.715   1.611   558.150;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-22405.8, old_max_log_p =-23359.3 (thresh=-23336.0)
 0.914   3.415   1.087  -894.601;
-1.822   1.660  -1.372   243.436;
-3.239  -1.562   2.069   505.402;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-21265.5, old_max_log_p =-22405.8 (thresh=-22383.4)
 0.400   3.974   0.708  -888.302;
-1.833   1.034  -2.260   526.747;
-4.167  -1.740   2.180   644.850;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-21145.7, old_max_log_p =-21265.5 (thresh=-21244.2)
-0.195   3.923   1.068  -869.113;
-1.833   1.034  -2.260   526.747;
-4.182  -2.276   2.060   780.738;
 0.000   0.000   0.000   1.000;
Result so far: scale 1.000: max_log_p=-21145.7, old_max_log_p =-21145.7 (thresh=-21124.5)
-0.195   3.923   1.068  -869.113;
-1.833   1.034  -2.260   526.747;
-4.182  -2.276   2.060   780.738;
 0.000   0.000   0.000   1.000;
dt = 5.00e-06, momentum=0.80, tol=1.00e-05
l_intensity = 1.0000
Aligning input volume to GCA...
Transform matrix
-0.19499   3.92350   1.06768  -869.11322;
-1.83280   1.03418  -2.25956   526.74725;
-4.18182  -2.27635   2.06046   780.73822;
 0.00000   0.00000   0.00000   1.00000;
nsamples 3516
Quasinewton: input matrix
-0.19499   3.92350   1.06768  -869.11322;
-1.83280   1.03418  -2.25956   526.74725;
-4.18182  -2.27635   2.06046   780.73822;
 0.00000   0.00000   0.00000   1.00000;
fp = 21145.679688, fret = 21145.679688, ftol = 0.000010, EPS = 0.000000
dt = 5.00e-06, momentum=0.80, tol=1.00e-07
l_intensity = 1.0000
Aligning input volume to GCA...
Transform matrix
-0.19499   3.92350   1.06768  -869.11322;
-1.83280   1.03418  -2.25956   526.74725;
-4.18182  -2.27635   2.06046   780.73822;
 0.00000   0.00000   0.00000   1.00000;
nsamples 385836
Quasinewton: input matrix
-0.19499   3.92350   1.06768  -869.11322;
-1.83280   1.03418  -2.25956   526.74725;
-4.18182  -2.27635   2.06046   780.73822;
 0.00000   0.00000   0.00000   1.00000;
pass 2 through quasi-newton minimization...
pass 3 through quasi-newton minimization...
fp = 2757868.250000, fret = 2757868.250000, ftol = 0.000000, EPS = 0.000000
freeing gibbs priors...done.
bounding unknown intensity as < 20.7 or > 1320.6 
************************************************
spacing=8, using 3516 sample points, tol=1.00e-05...
************************************************
resetting wm mean[0]: 117 --> 126
resetting gm mean[0]: 73 --> 73
input volume #1 is the most T1-like
using real data threshold=10.0
using (98, 192, 170) as brain centroid...
mean wm in atlas = 126, using box (82,180,152) --> (113, 203,188) to find MRI wm
before smoothing, mri peak at 66
after smoothing, mri peak at 67, scaling input intensities by 1.881
scaling channel 0 by 1.8806
****************************************
Nine parameter search.  iteration 0 nscales = 0 ...
****************************************
****************************************
Nine parameter search.  iteration 1 nscales = 1 ...
****************************************
****************************************
Nine parameter search.  iteration 2 nscales = 2 ...
****************************************
****************************************
Nine parameter search.  iteration 3 nscales = 3 ...
****************************************
****************************************
Nine parameter search.  iteration 4 nscales = 4 ...
****************************************
****************************************
Nine parameter search.  iteration 5 nscales = 5 ...
****************************************
****************************************
Nine parameter search.  iteration 6 nscales = 6 ...
****************************************
****************************************
Nine parameter search.  iteration 7 nscales = 7 ...
****************************************
reducing scale to 0.2500
***********************************************
Computing MAP estimate using 3516 samples...
***********************************************
dfp_em_step_func: 008: -log(p) = 21145.7
after pass:transform: ( -0.19, 3.92, 1.07, -869.11)
                      ( -1.83, 1.03, -2.26, 526.75)
                      ( -4.18, -2.28, 2.06, 780.74)
outof QuasiNewtonEMA: 010: -log(p) = 21145.7  tol 0.000010
Resulting transform:
-0.195   3.923   1.068  -869.113;
-1.833   1.034  -2.260   526.747;
-4.182  -2.276   2.060   780.738;
 0.000   0.000   0.000   1.000;

pass 1, spacing 8: log(p) = -21145.7 (old=-2022144.0)
transform before final EM align:
-0.195   3.923   1.068  -869.113;
-1.833   1.034  -2.260   526.747;
-4.182  -2.276   2.060   780.738;
 0.000   0.000   0.000   1.000;

**************************************************
 EM alignment process ...
 Computing final MAP estimate using 385836 samples. 
**************************************************
dfp_em_step_func: 009: -log(p) = 2757924.8
after pass:transform: ( -0.19, 3.92, 1.07, -869.11)
                      ( -1.83, 1.03, -2.26, 526.75)
                      ( -4.18, -2.28, 2.06, 780.74)
dfp_em_step_func: 010: -log(p) = 2757880.2
after pass:transform: ( -0.19, 3.92, 1.07, -869.11)
                      ( -1.83, 1.03, -2.26, 526.75)
                      ( -4.18, -2.28, 2.06, 780.74)
dfp_em_step_func: 011: -log(p) = 2757868.2
after pass:transform: ( -0.19, 3.92, 1.07, -869.11)
                      ( -1.83, 1.03, -2.26, 526.75)
                      ( -4.18, -2.28, 2.06, 780.74)
outof QuasiNewtonEMA: 013: -log(p) = 2757868.2  tol 0.000000
final transform:
-0.195   3.924   1.068  -869.113;
-1.833   1.034  -2.260   526.747;
-4.182  -2.276   2.060   780.738;
 0.000   0.000   0.000   1.000;

writing output transformation to transforms/talairach_with_skull.lta...
registration took 29 minutes and 46 seconds.
#--------------------------------------
#@# SubCort Seg Thu Jun 21 01:06:52 EDT 2007

 mri_ca_label -align -nobigventricles norm.mgz transforms/talairach.m3z /usr/local/freesurfer/average/RB_all_2006-02-15.gca aseg.auto.mgz 

renormalizing sequences with structure alignment, equivalent to:
	-renormalize
	-renormalize_mean 0.500
	-regularize 0.500
not handling expanded ventricles...
reading 1 input volumes...
reading classifier array from /usr/local/freesurfer/average/RB_all_2006-02-15.gca...
setting gca type = Normal gca type
reading input volume from norm.mgz...
average std[0] = 7.4
reading transform from transforms/talairach.m3z...
gcam->type = vox
reading labels out of gcam file...
Atlas used for the 3D morph was /usr/local/freesurfer/average/RB_all_2006-02-15.gca
average std = 7.4   using min determinant for regularization = 5.5
0 singular and 0 ill-conditioned covariance matrices regularized
labeling volume...
renormalizing by structure alignment....
renormalizing input #0
gca peak = 0.13482 (104)
mri peak = 0.05088 (113)
Left_Pallidum (13): linear fit = 1.08 x + 0.0 (715 voxels, overlap=0.510)
gca peak = 0.19307 (105)
uniform distribution in MR - rejecting arbitrary fit
gca peak = 0.07702 (114)
mri peak = 0.03220 (96)
Right_Cerebral_White_Matter (41): linear fit = 0.95 x + 0.0 (1782 voxels, overlap=0.017)
gca peak = 0.07016 (114)
mri peak = 0.03119 (95)
Left_Cerebral_White_Matter (2): linear fit = 0.90 x + 0.0 (9356 voxels, overlap=0.046)
gca peak = 0.10443 (32)
mri peak = 0.05504 (107)
CSF peak too bright - rejecting
gca peak = 0.11491 (25)
mri peak = 0.07306 (51)
gca peak = 0.27319 (66)
uniform distribution in MR - rejecting arbitrary fit
gca peak = 0.29477 (66)
mri peak = 0.04022 (93)
gca peak = 0.06920 (68)
mri peak = 0.04607 (34)
unreasonable value (27.2), not in range [40, 95] - rejecting
gca peak = 0.08287 (69)
mri peak = 0.02440 (34)
unreasonable value (29.3), not in range [40, 95] - rejecting
gca peak = 0.22476 (80)
mri peak = 0.13063 (61)
gca peak = 0.18589 (81)
mri peak = 0.04550 (107)
unreasonable value (105.3), not in range [50, 100] - rejecting
gca peak = 0.08091 (61)
mri peak = 0.02814 (44)
Left_Cerebellum_Cortex (8): linear fit = 0.68 x + 0.0 (844 voxels, overlap=0.011)
gca peak = 0.10638 (61)
mri peak = 0.03403 (221)
Right_Cerebellum_Cortex (47): linear fit = 0.98 x + 0.0 (1446 voxels, overlap=0.000)
gca peak = 0.15522 (88)
mri peak = 0.07889 (91)
Left_Cerebellum_White_Matter (7): linear fit = 1.05 x + 0.0 (307 voxels, overlap=0.746)
gca peak = 0.12507 (87)
mri peak = 0.06386 (222)
Right_Cerebellum_White_Matter (46): linear fit = 2.50 x + 0.0 (831 voxels, overlap=0.000)
gca peak = 0.25800 (69)
mri peak = 0.05479 (103)
gca peak = 0.29193 (67)
uniform distribution in MR - rejecting arbitrary fit
gca peak = 0.11200 (101)
mri peak = 0.05132 (94)
Left_Thalamus_Proper (10): linear fit = 0.93 x + 0.0 (2038 voxels, overlap=0.390)
gca peak = 0.10010 (92)
mri peak = 0.04324 (94)
Right_Thalamus_Proper (49): linear fit = 1.02 x + 0.0 (755 voxels, overlap=0.579)
gca peak = 0.07167 (89)
mri peak = 0.03957 (119)
unreasonable value (117.9), not in range [60, 100] - rejecting
gca peak = 0.09514 (89)
mri peak = 0.40000 (95)
gca peak = 0.07535 (89)
mri peak = 0.03944 (82)
Brain_Stem (16): linear fit = 0.95 x + 0.0 (1746 voxels, overlap=0.174)
gca peak = 0.11225 (101)
mri peak = 0.10820 (69)
gca peak = 0.15721 (97)
mri peak = 0.04031 (62)
unreasonable value (67.9), not in range [80, 110] - rejecting
gca peak = 0.13106 (39)
uniform distribution in MR - rejecting arbitrary fit
gca peak = 0.19587 (25)
mri peak = 0.13559 (32)
not using caudate to estimate GM means
setting label Left_Cerebellum_White_Matter based on Right_Cerebellum_White_Matter = 1.05 x +  0
setting label Left_Pallidum based on Right_Pallidum = 1.08 x +  0
estimating mean gm scale to be 1.00 x + 0.0
estimating mean wm scale to be 0.92 x + 0.0
estimating mean csf scale to be 1.00 x + 0.0
4879 gm and wm labels changed (%27 to gray, %73 to white out of all changed labels)
151 hippocampal voxels changed.
0 amygdala voxels changed.
pass 1: 21996 changed. image ll: -20.834, PF=1.000
pass 2: 4522 changed.
4315 border labels changed to MLE ...
writing labeled volume to aseg.auto.mgz...
auto-labeling took 29 minutes and 15 seconds.

 cp aseg.auto.mgz aseg.mgz 

#--------------------------------------------
#@# ASeg Stats Thu Jun 21 01:36:09 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri

 mri_segstats --seg mri/aseg.mgz --sum stats/aseg.stats --pv mri/norm.mgz --ctab-default --excludeid 0 --brain-vol-from-seg --brainmask mri/brainmask.mgz --in mri/norm.mgz --in-intensity-name norm --in-intensity-units MR --etiv --subject sub089mri 

Using defalt ctab /usr/local/freesurfer/FreeSurferColorLUT.txt
atlas_icv = 34931.5
Loading mri/aseg.mgz
Loading mri/norm.mgz
Loading mri/norm.mgz
Loading mri/brainmask.mgz
# nbrainmaskvoxels 2981094
# brainmaskvolume  2981094.0
# nbrainsegvoxels 173772
# brainsegvolume   173772.0
Voxel Volume is 1 mm^3
Generating list of segmentation ids
Found 686 segmentations
Computing statistics for each segmentation
  1     1  Left-Cerebral-Exterior  0
  2     2  Left-Cerebral-White-Matter 42991  42991
  3     3  Left-Cerebral-Cortex 29681  29681
  4     4  Left-Lateral-Ventricle 1255  1255
  5     5  Left-Inf-Lat-Vent   85  85
  6     6  Left-Cerebellum-Exterior  0
  7     7  Left-Cerebellum-White-Matter 2796  2796
  8     8  Left-Cerebellum-Cortex 3708  3708
  9     9  Left-Thalamus  0
 10    10  Left-Thalamus-Proper 4596  4596
 11    11  Left-Caudate 1730  1730
 12    12  Left-Putamen 3657  3657
 13    13  Left-Pallidum 3093  3093
 14    14  3rd-Ventricle  289  289
 15    15  4th-Ventricle  398  398
 16    16  Brain-Stem 4401  4401
 17    17  Left-Hippocampus  715  715
 18    18  Left-Amygdala  147  147
 19    19  Left-Insula  0
 20    20  Left-Operculum  0
 21    21  Line-1  0
 22    22  Line-2  0
 23    23  Line-3  0
 24    24  CSF  289  289
 25    25  Left-Lesion  0
 26    26  Left-Accumbens-area   85  85
 27    27  Left-Substancia-Nigra  0
 28    28  Left-VentralDC 2266  2266
 29    29  Left-undetermined  0
 30    30  Left-vessel  143  143
 31    31  Left-choroid-plexus  0
 32    32  Left-F3orb  0
 33    33  Left-lOg  0
 34    34  Left-aOg  0
 35    35  Left-mOg  0
 36    36  Left-pOg  0
 37    37  Left-Stellate  0
 38    38  Left-Porg  0
 39    39  Left-Aorg  0
 40    40  Right-Cerebral-Exterior  0
 41    41  Right-Cerebral-White-Matter 19922  19922
 42    42  Right-Cerebral-Cortex 25978  25978
 43    43  Right-Lateral-Ventricle  805  805
 44    44  Right-Inf-Lat-Vent   13  13
 45    45  Right-Cerebellum-Exterior  0
 46    46  Right-Cerebellum-White-Matter 2822  2822
 47    47  Right-Cerebellum-Cortex 7157  7157
 48    48  Right-Thalamus  0
 49    49  Right-Thalamus-Proper 2639  2639
 50    50  Right-Caudate  463  463
 51    51  Right-Putamen  681  681
 52    52  Right-Pallidum   94  94
 53    53  Right-Hippocampus  402  402
 54    54  Right-Amygdala   40  40
 55    55  Right-Insula  0
 56    56  Right-Operculum  0
 57    57  Right-Lesion  0
 58    58  Right-Accumbens-area   15  15
 59    59  Right-Substancia-Nigra  0
 60    60  Right-VentralDC  472  472
 61    61  Right-undetermined  0
 62    62  Right-vessel   16  16
 63    63  Right-choroid-plexus  0
 64    64  Right-F3orb  0
 65    65  Right-lOg  0
 66    66  Right-aOg  0
 67    67  Right-mOg  0
 68    68  Right-pOg  0
 69    69  Right-Stellate  0
 70    70  Right-Porg  0
 71    71  Right-Aorg  0
 72    72  5th-Ventricle    8  8
 73    73  Left-Interior  0
 74    74  Right-Interior    7  7
 75    75  Left-Lateral-Ventricles  0
 76    76  Right-Lateral-Ventricles  0
 77    77  WM-hypointensities  0
 78    78  Left-WM-hypointensities 4331  4331
 79    79  Right-WM-hypointensities 2107  2107
 80    80  non-WM-hypointensities  0
 81    81  Left-non-WM-hypointensities  151  151
 82    82  Right-non-WM-hypointensities   27  27
 83    83  Left-F1  0
 84    84  Right-F1  0
 85    85  Optic-Chiasm  512  512
 86    86  Corpus_Callosum  0
 87    96  Left-Amygdala-Anterior  0
 88    97  Right-Amygdala-Anterior  0
 89    98  Dura  0
 90   100  Left-wm-intensity-abnormality  0
 91   101  Left-caudate-intensity-abnormality  0
 92   102  Left-putamen-intensity-abnormality  0
 93   103  Left-accumbens-intensity-abnormality  0
 94   104  Left-pallidum-intensity-abnormality  0
 95   105  Left-amygdala-intensity-abnormality  0
 96   106  Left-hippocampus-intensity-abnormality  0
 97   107  Left-thalamus-intensity-abnormality  0
 98   108  Left-VDC-intensity-abnormality  0
 99   109  Right-wm-intensity-abnormality  0
100   110  Right-caudate-intensity-abnormality  0
101   111  Right-putamen-intensity-abnormality  0
102   112  Right-accumbens-intensity-abnormality  0
103   113  Right-pallidum-intensity-abnormality  0
104   114  Right-amygdala-intensity-abnormality  0
105   115  Right-hippocampus-intensity-abnormality  0
106   116  Right-thalamus-intensity-abnormality  0
107   117  Right-VDC-intensity-abnormality  0
108   118  Epidermis  0
109   119  Conn-Tissue  0
110   120  SC-Fat/Muscle  0
111   121  Cranium  0
112   122  CSF-SA  0
113   123  Muscle  0
114   124  Ear  0
115   125  Adipose  0
116   126  Spinal-Cord  0
117   127  Soft-Tissue  0
118   128  Nerve  0
119   129  Bone  0
120   130  Air  0
121   131  Orbital-Fat  0
122   132  Tongue  0
123   133  Nasal-Structures  0
124   134  Globe  0
125   135  Teeth  0
126   136  Left-Caudate/Putamen  0
127   137  Right-Caudate/Putamen  0
128   138  Left-Claustrum  0
129   139  Right-Claustrum  0
130   140  Cornea  0
131   142  Diploe  0
132   143  Vitreous-Humor  0
133   144  Lens  0
134   145  Aqueous-Humor  0
135   146  Outer-Table  0
136   147  Inner-Table  0
137   148  Periosteum  0
138   149  Endosteum  0
139   150  R/C/S  0
140   151  Iris  0
141   152  SC-Adipose/Muscle  0
142   153  SC-Tissue  0
143   154  Orbital-Adipose  0
144   193  Left-hippocampal_fissure  0
145   194  Left-CADG-head  0
146   195  Left-subiculum  0
147   196  Left-fimbria  0
148   197  Right-hippocampal_fissure  0
149   198  Right-CADG-head  0
150   199  Right-subiculum  0
151   200  Right-fimbria  0
152   201  alveus  0
153   202  perforant_pathway  0
154   203  parasubiculum  0
155   204  presubiculum  0
156   205  subiculum  0
157   206  CA1  0
158   207  CA2  0
159   208  CA3  0
160   209  CA4  0
161   210  GC-DG  0
162   211  HATA  0
163   212  fimbria  0
164   213  lateral_ventricle  0
165   214  molecular_layer_HP  0
166   215  hippocampal_fissure  0
167   216  entorhinal_cortex  0
168   217  molecular_layer_subiculum  0
169   218  Amygdala  0
170   219  Cerebral_White_Matter  0
171   220  Cerebral_Cortex  0
172   221  Inf_Lat_Vent  0
173   222  Ectorhinal  0
174   222  Perirhinal  0
175   223  Cerebral_White_Matter_Edge  0
176   224  fMRI_Background  0
177   331  Aorta  0
178   332  Left-Common-IliacA  0
179   333  Right-Common-IliacA  0
180   334  Left-External-IliacA  0
181   335  Right-External-IliacA  0
182   336  Left-Internal-IliacA  0
183   337  Right-Internal-IliacA  0
184   338  Left-Lateral-SacralA  0
185   339  Right-Lateral-SacralA  0
186   340  Left-ObturatorA  0
187   341  Right-ObturatorA  0
188   342  Left-Internal-PudendalA  0
189   343  Right-Internal-PudendalA  0
190   344  Left-UmbilicalA  0
191   345  Right-UmbilicalA  0
192   346  Left-Inf-RectalA  0
193   347  Right-Inf-RectalA  0
194   348  Left-Common-IliacV  0
195   349  Right-Common-IliacV  0
196   350  Left-External-IliacV  0
197   351  Right-External-IliacV  0
198   352  Left-Internal-IliacV  0
199   353  Right-Internal-IliacV  0
200   354  Left-ObturatorV  0
201   355  Right-ObturatorV  0
202   356  Left-Internal-PudendalV  0
203   357  Right-Internal-PudendalV  0
204   358  Pos-Lymph  0
205   359  Neg-Lymph  0
206   400  V1  0
207   401  V2  0
208   402  BA44  0
209   403  BA45  0
210   404  BA4a  0
211   405  BA4p  0
212   406  BA6  0
213   407  BA2  0
214   408  BAun1  0
215   409  BAun2  0
216   1000  ctx-lh-unknown  0
217   1001  ctx-lh-bankssts  0
218   1002  ctx-lh-caudalanteriorcingulate  0
219   1003  ctx-lh-caudalmiddlefrontal  0
220   1004  ctx-lh-corpuscallosum  0
221   1005  ctx-lh-cuneus  0
222   1006  ctx-lh-entorhinal  0
223   1007  ctx-lh-fusiform  0
224   1008  ctx-lh-inferiorparietal  0
225   1009  ctx-lh-inferiortemporal  0
226   1010  ctx-lh-isthmuscingulate  0
227   1011  ctx-lh-lateraloccipital  0
228   1012  ctx-lh-lateralorbitofrontal  0
229   1013  ctx-lh-lingual  0
230   1014  ctx-lh-medialorbitofrontal  0
231   1015  ctx-lh-middletemporal  0
232   1016  ctx-lh-parahippocampal  0
233   1017  ctx-lh-paracentral  0
234   1018  ctx-lh-parsopercularis  0
235   1019  ctx-lh-parsorbitalis  0
236   1020  ctx-lh-parstriangularis  0
237   1021  ctx-lh-pericalcarine  0
238   1022  ctx-lh-postcentral  0
239   1023  ctx-lh-posteriorcingulate  0
240   1024  ctx-lh-precentral  0
241   1025  ctx-lh-precuneus  0
242   1026  ctx-lh-rostralanteriorcingulate  0
243   1027  ctx-lh-rostralmiddlefrontal  0
244   1028  ctx-lh-superiorfrontal  0
245   1029  ctx-lh-superiorparietal  0
246   1030  ctx-lh-superiortemporal  0
247   1031  ctx-lh-supramarginal  0
248   1032  ctx-lh-frontalpole  0
249   1033  ctx-lh-temporalpole  0
250   1034  ctx-lh-transversetemporal  0
251   2000  ctx-rh-unknown  0
252   2001  ctx-rh-bankssts  0
253   2002  ctx-rh-caudalanteriorcingulate  0
254   2003  ctx-rh-caudalmiddlefrontal  0
255   2004  ctx-rh-corpuscallosum  0
256   2005  ctx-rh-cuneus  0
257   2006  ctx-rh-entorhinal  0
258   2007  ctx-rh-fusiform  0
259   2008  ctx-rh-inferiorparietal  0
260   2009  ctx-rh-inferiortemporal  0
261   2010  ctx-rh-isthmuscingulate  0
262   2011  ctx-rh-lateraloccipital  0
263   2012  ctx-rh-lateralorbitofrontal  0
264   2013  ctx-rh-lingual  0
265   2014  ctx-rh-medialorbitofrontal  0
266   2015  ctx-rh-middletemporal  0
267   2016  ctx-rh-parahippocampal  0
268   2017  ctx-rh-paracentral  0
269   2018  ctx-rh-parsopercularis  0
270   2019  ctx-rh-parsorbitalis  0
271   2020  ctx-rh-parstriangularis  0
272   2021  ctx-rh-pericalcarine  0
273   2022  ctx-rh-postcentral  0
274   2023  ctx-rh-posteriorcingulate  0
275   2024  ctx-rh-precentral  0
276   2025  ctx-rh-precuneus  0
277   2026  ctx-rh-rostralanteriorcingulate  0
278   2027  ctx-rh-rostralmiddlefrontal  0
279   2028  ctx-rh-superiorfrontal  0
280   2029  ctx-rh-superiorparietal  0
281   2030  ctx-rh-superiortemporal  0
282   2031  ctx-rh-supramarginal  0
283   2032  ctx-rh-frontalpole  0
284   2033  ctx-rh-temporalpole  0
285   2034  ctx-rh-transversetemporal  0
286   3000  wm-lh-unknown  0
287   3001  wm-lh-bankssts  0
288   3002  wm-lh-caudalanteriorcingulate  0
289   3003  wm-lh-caudalmiddlefrontal  0
290   3004  wm-lh-corpuscallosum  0
291   3005  wm-lh-cuneus  0
292   3006  wm-lh-entorhinal  0
293   3007  wm-lh-fusiform  0
294   3008  wm-lh-inferiorparietal  0
295   3009  wm-lh-inferiortemporal  0
296   3010  wm-lh-isthmuscingulate  0
297   3011  wm-lh-lateraloccipital  0
298   3012  wm-lh-lateralorbitofrontal  0
299   3013  wm-lh-lingual  0
300   3014  wm-lh-medialorbitofrontal  0
301   3015  wm-lh-middletemporal  0
302   3016  wm-lh-parahippocampal  0
303   3017  wm-lh-paracentral  0
304   3018  wm-lh-parsopercularis  0
305   3019  wm-lh-parsorbitalis  0
306   3020  wm-lh-parstriangularis  0
307   3021  wm-lh-pericalcarine  0
308   3022  wm-lh-postcentral  0
309   3023  wm-lh-posteriorcingulate  0
310   3024  wm-lh-precentral  0
311   3025  wm-lh-precuneus  0
312   3026  wm-lh-rostralanteriorcingulate  0
313   3027  wm-lh-rostralmiddlefrontal  0
314   3028  wm-lh-superiorfrontal  0
315   3029  wm-lh-superiorparietal  0
316   3030  wm-lh-superiortemporal  0
317   3031  wm-lh-supramarginal  0
318   3032  wm-lh-frontalpole  0
319   3033  wm-lh-temporalpole  0
320   3034  wm-lh-transversetemporal  0
321   4000  wm-rh-unknown  0
322   4001  wm-rh-bankssts  0
323   4002  wm-rh-caudalanteriorcingulate  0
324   4003  wm-rh-caudalmiddlefrontal  0
325   4004  wm-rh-corpuscallosum  0
326   4005  wm-rh-cuneus  0
327   4006  wm-rh-entorhinal  0
328   4007  wm-rh-fusiform  0
329   4008  wm-rh-inferiorparietal  0
330   4009  wm-rh-inferiortemporal  0
331   4010  wm-rh-isthmuscingulate  0
332   4011  wm-rh-lateraloccipital  0
333   4012  wm-rh-lateralorbitofrontal  0
334   4013  wm-rh-lingual  0
335   4014  wm-rh-medialorbitofrontal  0
336   4015  wm-rh-middletemporal  0
337   4016  wm-rh-parahippocampal  0
338   4017  wm-rh-paracentral  0
339   4018  wm-rh-parsopercularis  0
340   4019  wm-rh-parsorbitalis  0
341   4020  wm-rh-parstriangularis  0
342   4021  wm-rh-pericalcarine  0
343   4022  wm-rh-postcentral  0
344   4023  wm-rh-posteriorcingulate  0
345   4024  wm-rh-precentral  0
346   4025  wm-rh-precuneus  0
347   4026  wm-rh-rostralanteriorcingulate  0
348   4027  wm-rh-rostralmiddlefrontal  0
349   4028  wm-rh-superiorfrontal  0
350   4029  wm-rh-superiorparietal  0
351   4030  wm-rh-superiortemporal  0
352   4031  wm-rh-supramarginal  0
353   4032  wm-rh-frontalpole  0
354   4033  wm-rh-temporalpole  0
355   4034  wm-rh-transversetemporal  0
356   5001  wm-lh-centrum-semiovale  0
357   5002  wm-rh-centrum-semiovale  0
358   1100  ctx-lh-Unknown  0
359   1101  ctx-lh-Corpus_callosum  0
360   1102  ctx-lh-G_and_S_Insula_ONLY_AVERAGE  0
361   1103  ctx-lh-G_cingulate-Isthmus  0
362   1104  ctx-lh-G_cingulate-Main_part  0
363   1105  ctx-lh-G_cuneus  0
364   1106  ctx-lh-G_frontal_inf-Opercular_part  0
365   1107  ctx-lh-G_frontal_inf-Orbital_part  0
366   1108  ctx-lh-G_frontal_inf-Triangular_part  0
367   1109  ctx-lh-G_frontal_middle  0
368   1110  ctx-lh-G_frontal_superior  0
369   1111  ctx-lh-G_frontomarginal  0
370   1112  ctx-lh-G_insular_long  0
371   1113  ctx-lh-G_insular_short  0
372   1114  ctx-lh-G_and_S_occipital_inferior  0
373   1115  ctx-lh-G_occipital_middle  0
374   1116  ctx-lh-G_occipital_superior  0
375   1117  ctx-lh-G_occipit-temp_lat-Or_fusiform  0
376   1118  ctx-lh-G_occipit-temp_med-Lingual_part  0
377   1119  ctx-lh-G_occipit-temp_med-Parahippocampal_part  0
378   1120  ctx-lh-G_orbital  0
379   1121  ctx-lh-G_paracentral  0
380   1122  ctx-lh-G_parietal_inferior-Angular_part  0
381   1123  ctx-lh-G_parietal_inferior-Supramarginal_part  0
382   1124  ctx-lh-G_parietal_superior  0
383   1125  ctx-lh-G_postcentral  0
384   1126  ctx-lh-G_precentral  0
385   1127  ctx-lh-G_precuneus  0
386   1128  ctx-lh-G_rectus  0
387   1129  ctx-lh-G_subcallosal  0
388   1130  ctx-lh-G_subcentral  0
389   1131  ctx-lh-G_temporal_inferior  0
390   1132  ctx-lh-G_temporal_middle  0
391   1133  ctx-lh-G_temp_sup-G_temp_transv_and_interm_S  0
392   1134  ctx-lh-G_temp_sup-Lateral_aspect  0
393   1135  ctx-lh-G_temp_sup-Planum_polare  0
394   1136  ctx-lh-G_temp_sup-Planum_tempolare  0
395   1137  ctx-lh-G_and_S_transverse_frontopolar  0
396   1138  ctx-lh-Lat_Fissure-ant_sgt-ramus_horizontal  0
397   1139  ctx-lh-Lat_Fissure-ant_sgt-ramus_vertical  0
398   1140  ctx-lh-Lat_Fissure-post_sgt  0
399   1141  ctx-lh-Medial_wall  0
400   1142  ctx-lh-Pole_occipital  0
401   1143  ctx-lh-Pole_temporal  0
402   1144  ctx-lh-S_calcarine  0
403   1145  ctx-lh-S_central  0
404   1146  ctx-lh-S_central_insula  0
405   1147  ctx-lh-S_cingulate-Main_part_and_Intracingulate  0
406   1148  ctx-lh-S_cingulate-Marginalis_part  0
407   1149  ctx-lh-S_circular_insula_anterior  0
408   1150  ctx-lh-S_circular_insula_inferior  0
409   1151  ctx-lh-S_circular_insula_superior  0
410   1152  ctx-lh-S_collateral_transverse_ant  0
411   1153  ctx-lh-S_collateral_transverse_post  0
412   1154  ctx-lh-S_frontal_inferior  0
413   1155  ctx-lh-S_frontal_middle  0
414   1156  ctx-lh-S_frontal_superior  0
415   1157  ctx-lh-S_frontomarginal  0
416   1158  ctx-lh-S_intermedius_primus-Jensen  0
417   1159  ctx-lh-S_intraparietal-and_Parietal_transverse  0
418   1160  ctx-lh-S_occipital_anterior  0
419   1161  ctx-lh-S_occipital_middle_and_Lunatus  0
420   1162  ctx-lh-S_occipital_superior_and_transversalis  0
421   1163  ctx-lh-S_occipito-temporal_lateral  0
422   1164  ctx-lh-S_occipito-temporal_medial_and_S_Lingual  0
423   1165  ctx-lh-S_orbital-H_shapped  0
424   1166  ctx-lh-S_orbital_lateral  0
425   1167  ctx-lh-S_orbital_medial-Or_olfactory  0
426   1168  ctx-lh-S_paracentral  0
427   1169  ctx-lh-S_parieto_occipital  0
428   1170  ctx-lh-S_pericallosal  0
429   1171  ctx-lh-S_postcentral  0
430   1172  ctx-lh-S_precentral-Inferior-part  0
431   1173  ctx-lh-S_precentral-Superior-part  0
432   1174  ctx-lh-S_subcentral_ant  0
433   1175  ctx-lh-S_subcentral_post  0
434   1176  ctx-lh-S_suborbital  0
435   1177  ctx-lh-S_subparietal  0
436   1178  ctx-lh-S_supracingulate  0
437   1179  ctx-lh-S_temporal_inferior  0
438   1180  ctx-lh-S_temporal_superior  0
439   1181  ctx-lh-S_temporal_transverse  0
440   2100  ctx-rh-Unknown  0
441   2101  ctx-rh-Corpus_callosum  0
442   2102  ctx-rh-G_and_S_Insula_ONLY_AVERAGE  0
443   2103  ctx-rh-G_cingulate-Isthmus  0
444   2104  ctx-rh-G_cingulate-Main_part  0
445   2105  ctx-rh-G_cuneus  0
446   2106  ctx-rh-G_frontal_inf-Opercular_part  0
447   2107  ctx-rh-G_frontal_inf-Orbital_part  0
448   2108  ctx-rh-G_frontal_inf-Triangular_part  0
449   2109  ctx-rh-G_frontal_middle  0
450   2110  ctx-rh-G_frontal_superior  0
451   2111  ctx-rh-G_frontomarginal  0
452   2112  ctx-rh-G_insular_long  0
453   2113  ctx-rh-G_insular_short  0
454   2114  ctx-rh-G_and_S_occipital_inferior  0
455   2115  ctx-rh-G_occipital_middle  0
456   2116  ctx-rh-G_occipital_superior  0
457   2117  ctx-rh-G_occipit-temp_lat-Or_fusiform  0
458   2118  ctx-rh-G_occipit-temp_med-Lingual_part  0
459   2119  ctx-rh-G_occipit-temp_med-Parahippocampal_part  0
460   2120  ctx-rh-G_orbital  0
461   2121  ctx-rh-G_paracentral  0
462   2122  ctx-rh-G_parietal_inferior-Angular_part  0
463   2123  ctx-rh-G_parietal_inferior-Supramarginal_part  0
464   2124  ctx-rh-G_parietal_superior  0
465   2125  ctx-rh-G_postcentral  0
466   2126  ctx-rh-G_precentral  0
467   2127  ctx-rh-G_precuneus  0
468   2128  ctx-rh-G_rectus  0
469   2129  ctx-rh-G_subcallosal  0
470   2130  ctx-rh-G_subcentral  0
471   2131  ctx-rh-G_temporal_inferior  0
472   2132  ctx-rh-G_temporal_middle  0
473   2133  ctx-rh-G_temp_sup-G_temp_transv_and_interm_S  0
474   2134  ctx-rh-G_temp_sup-Lateral_aspect  0
475   2135  ctx-rh-G_temp_sup-Planum_polare  0
476   2136  ctx-rh-G_temp_sup-Planum_tempolare  0
477   2137  ctx-rh-G_and_S_transverse_frontopolar  0
478   2138  ctx-rh-Lat_Fissure-ant_sgt-ramus_horizontal  0
479   2139  ctx-rh-Lat_Fissure-ant_sgt-ramus_vertical  0
480   2140  ctx-rh-Lat_Fissure-post_sgt  0
481   2141  ctx-rh-Medial_wall  0
482   2142  ctx-rh-Pole_occipital  0
483   2143  ctx-rh-Pole_temporal  0
484   2144  ctx-rh-S_calcarine  0
485   2145  ctx-rh-S_central  0
486   2146  ctx-rh-S_central_insula  0
487   2147  ctx-rh-S_cingulate-Main_part_and_Intracingulate  0
488   2148  ctx-rh-S_cingulate-Marginalis_part  0
489   2149  ctx-rh-S_circular_insula_anterior  0
490   2150  ctx-rh-S_circular_insula_inferior  0
491   2151  ctx-rh-S_circular_insula_superior  0
492   2152  ctx-rh-S_collateral_transverse_ant  0
493   2153  ctx-rh-S_collateral_transverse_post  0
494   2154  ctx-rh-S_frontal_inferior  0
495   2155  ctx-rh-S_frontal_middle  0
496   2156  ctx-rh-S_frontal_superior  0
497   2157  ctx-rh-S_frontomarginal  0
498   2158  ctx-rh-S_intermedius_primus-Jensen  0
499   2159  ctx-rh-S_intraparietal-and_Parietal_transverse  0
500   2160  ctx-rh-S_occipital_anterior  0
501   2161  ctx-rh-S_occipital_middle_and_Lunatus  0
502   2162  ctx-rh-S_occipital_superior_and_transversalis  0
503   2163  ctx-rh-S_occipito-temporal_lateral  0
504   2164  ctx-rh-S_occipito-temporal_medial_and_S_Lingual  0
505   2165  ctx-rh-S_orbital-H_shapped  0
506   2166  ctx-rh-S_orbital_lateral  0
507   2167  ctx-rh-S_orbital_medial-Or_olfactory  0
508   2168  ctx-rh-S_paracentral  0
509   2169  ctx-rh-S_parieto_occipital  0
510   2170  ctx-rh-S_pericallosal  0
511   2171  ctx-rh-S_postcentral  0
512   2172  ctx-rh-S_precentral-Inferior-part  0
513   2173  ctx-rh-S_precentral-Superior-part  0
514   2174  ctx-rh-S_subcentral_ant  0
515   2175  ctx-rh-S_subcentral_post  0
516   2176  ctx-rh-S_suborbital  0
517   2177  ctx-rh-S_subparietal  0
518   2178  ctx-rh-S_supracingulate  0
519   2179  ctx-rh-S_temporal_inferior  0
520   2180  ctx-rh-S_temporal_superior  0
521   2181  ctx-rh-S_temporal_transverse  0
522   3100  wm-lh-Unknown  0
523   3101  wm-lh-Corpus_callosum  0
524   3102  wm-lh-G_and_S_Insula_ONLY_AVERAGE  0
525   3103  wm-lh-G_cingulate-Isthmus  0
526   3104  wm-lh-G_cingulate-Main_part  0
527   3105  wm-lh-G_cuneus  0
528   3106  wm-lh-G_frontal_inf-Opercular_part  0
529   3107  wm-lh-G_frontal_inf-Orbital_part  0
530   3108  wm-lh-G_frontal_inf-Triangular_part  0
531   3109  wm-lh-G_frontal_middle  0
532   3110  wm-lh-G_frontal_superior  0
533   3111  wm-lh-G_frontomarginal  0
534   3112  wm-lh-G_insular_long  0
535   3113  wm-lh-G_insular_short  0
536   3114  wm-lh-G_and_S_occipital_inferior  0
537   3115  wm-lh-G_occipital_middle  0
538   3116  wm-lh-G_occipital_superior  0
539   3117  wm-lh-G_occipit-temp_lat-Or_fusiform  0
540   3118  wm-lh-G_occipit-temp_med-Lingual_part  0
541   3119  wm-lh-G_occipit-temp_med-Parahippocampal_part  0
542   3120  wm-lh-G_orbital  0
543   3121  wm-lh-G_paracentral  0
544   3122  wm-lh-G_parietal_inferior-Angular_part  0
545   3123  wm-lh-G_parietal_inferior-Supramarginal_part  0
546   3124  wm-lh-G_parietal_superior  0
547   3125  wm-lh-G_postcentral  0
548   3126  wm-lh-G_precentral  0
549   3127  wm-lh-G_precuneus  0
550   3128  wm-lh-G_rectus  0
551   3129  wm-lh-G_subcallosal  0
552   3130  wm-lh-G_subcentral  0
553   3131  wm-lh-G_temporal_inferior  0
554   3132  wm-lh-G_temporal_middle  0
555   3133  wm-lh-G_temp_sup-G_temp_transv_and_interm_S  0
556   3134  wm-lh-G_temp_sup-Lateral_aspect  0
557   3135  wm-lh-G_temp_sup-Planum_polare  0
558   3136  wm-lh-G_temp_sup-Planum_tempolare  0
559   3137  wm-lh-G_and_S_transverse_frontopolar  0
560   3138  wm-lh-Lat_Fissure-ant_sgt-ramus_horizontal  0
561   3139  wm-lh-Lat_Fissure-ant_sgt-ramus_vertical  0
562   3140  wm-lh-Lat_Fissure-post_sgt  0
563   3141  wm-lh-Medial_wall  0
564   3142  wm-lh-Pole_occipital  0
565   3143  wm-lh-Pole_temporal  0
566   3144  wm-lh-S_calcarine  0
567   3145  wm-lh-S_central  0
568   3146  wm-lh-S_central_insula  0
569   3147  wm-lh-S_cingulate-Main_part_and_Intracingulate  0
570   3148  wm-lh-S_cingulate-Marginalis_part  0
571   3149  wm-lh-S_circular_insula_anterior  0
572   3150  wm-lh-S_circular_insula_inferior  0
573   3151  wm-lh-S_circular_insula_superior  0
574   3152  wm-lh-S_collateral_transverse_ant  0
575   3153  wm-lh-S_collateral_transverse_post  0
576   3154  wm-lh-S_frontal_inferior  0
577   3155  wm-lh-S_frontal_middle  0
578   3156  wm-lh-S_frontal_superior  0
579   3157  wm-lh-S_frontomarginal  0
580   3158  wm-lh-S_intermedius_primus-Jensen  0
581   3159  wm-lh-S_intraparietal-and_Parietal_transverse  0
582   3160  wm-lh-S_occipital_anterior  0
583   3161  wm-lh-S_occipital_middle_and_Lunatus  0
584   3162  wm-lh-S_occipital_superior_and_transversalis  0
585   3163  wm-lh-S_occipito-temporal_lateral  0
586   3164  wm-lh-S_occipito-temporal_medial_and_S_Lingual  0
587   3165  wm-lh-S_orbital-H_shapped  0
588   3166  wm-lh-S_orbital_lateral  0
589   3167  wm-lh-S_orbital_medial-Or_olfactory  0
590   3168  wm-lh-S_paracentral  0
591   3169  wm-lh-S_parieto_occipital  0
592   3170  wm-lh-S_pericallosal  0
593   3171  wm-lh-S_postcentral  0
594   3172  wm-lh-S_precentral-Inferior-part  0
595   3173  wm-lh-S_precentral-Superior-part  0
596   3174  wm-lh-S_subcentral_ant  0
597   3175  wm-lh-S_subcentral_post  0
598   3176  wm-lh-S_suborbital  0
599   3177  wm-lh-S_subparietal  0
600   3178  wm-lh-S_supracingulate  0
601   3179  wm-lh-S_temporal_inferior  0
602   3180  wm-lh-S_temporal_superior  0
603   3181  wm-lh-S_temporal_transverse  0
604   4100  wm-rh-Unknown  0
605   4101  wm-rh-Corpus_callosum  0
606   4102  wm-rh-G_and_S_Insula_ONLY_AVERAGE  0
607   4103  wm-rh-G_cingulate-Isthmus  0
608   4104  wm-rh-G_cingulate-Main_part  0
609   4105  wm-rh-G_cuneus  0
610   4106  wm-rh-G_frontal_inf-Opercular_part  0
611   4107  wm-rh-G_frontal_inf-Orbital_part  0
612   4108  wm-rh-G_frontal_inf-Triangular_part  0
613   4109  wm-rh-G_frontal_middle  0
614   4110  wm-rh-G_frontal_superior  0
615   4111  wm-rh-G_frontomarginal  0
616   4112  wm-rh-G_insular_long  0
617   4113  wm-rh-G_insular_short  0
618   4114  wm-rh-G_and_S_occipital_inferior  0
619   4115  wm-rh-G_occipital_middle  0
620   4116  wm-rh-G_occipital_superior  0
621   4117  wm-rh-G_occipit-temp_lat-Or_fusiform  0
622   4118  wm-rh-G_occipit-temp_med-Lingual_part  0
623   4119  wm-rh-G_occipit-temp_med-Parahippocampal_part  0
624   4120  wm-rh-G_orbital  0
625   4121  wm-rh-G_paracentral  0
626   4122  wm-rh-G_parietal_inferior-Angular_part  0
627   4123  wm-rh-G_parietal_inferior-Supramarginal_part  0
628   4124  wm-rh-G_parietal_superior  0
629   4125  wm-rh-G_postcentral  0
630   4126  wm-rh-G_precentral  0
631   4127  wm-rh-G_precuneus  0
632   4128  wm-rh-G_rectus  0
633   4129  wm-rh-G_subcallosal  0
634   4130  wm-rh-G_subcentral  0
635   4131  wm-rh-G_temporal_inferior  0
636   4132  wm-rh-G_temporal_middle  0
637   4133  wm-rh-G_temp_sup-G_temp_transv_and_interm_S  0
638   4134  wm-rh-G_temp_sup-Lateral_aspect  0
639   4135  wm-rh-G_temp_sup-Planum_polare  0
640   4136  wm-rh-G_temp_sup-Planum_tempolare  0
641   4137  wm-rh-G_and_S_transverse_frontopolar  0
642   4138  wm-rh-Lat_Fissure-ant_sgt-ramus_horizontal  0
643   4139  wm-rh-Lat_Fissure-ant_sgt-ramus_vertical  0
644   4140  wm-rh-Lat_Fissure-post_sgt  0
645   4141  wm-rh-Medial_wall  0
646   4142  wm-rh-Pole_occipital  0
647   4143  wm-rh-Pole_temporal  0
648   4144  wm-rh-S_calcarine  0
649   4145  wm-rh-S_central  0
650   4146  wm-rh-S_central_insula  0
651   4147  wm-rh-S_cingulate-Main_part_and_Intracingulate  0
652   4148  wm-rh-S_cingulate-Marginalis_part  0
653   4149  wm-rh-S_circular_insula_anterior  0
654   4150  wm-rh-S_circular_insula_inferior  0
655   4151  wm-rh-S_circular_insula_superior  0
656   4152  wm-rh-S_collateral_transverse_ant  0
657   4153  wm-rh-S_collateral_transverse_post  0
658   4154  wm-rh-S_frontal_inferior  0
659   4155  wm-rh-S_frontal_middle  0
660   4156  wm-rh-S_frontal_superior  0
661   4157  wm-rh-S_frontomarginal  0
662   4158  wm-rh-S_intermedius_primus-Jensen  0
663   4159  wm-rh-S_intraparietal-and_Parietal_transverse  0
664   4160  wm-rh-S_occipital_anterior  0
665   4161  wm-rh-S_occipital_middle_and_Lunatus  0
666   4162  wm-rh-S_occipital_superior_and_transversalis  0
667   4163  wm-rh-S_occipito-temporal_lateral  0
668   4164  wm-rh-S_occipito-temporal_medial_and_S_Lingual  0
669   4165  wm-rh-S_orbital-H_shapped  0
670   4166  wm-rh-S_orbital_lateral  0
671   4167  wm-rh-S_orbital_medial-Or_olfactory  0
672   4168  wm-rh-S_paracentral  0
673   4169  wm-rh-S_parieto_occipital  0
674   4170  wm-rh-S_pericallosal  0
675   4171  wm-rh-S_postcentral  0
676   4172  wm-rh-S_precentral-Inferior-part  0
677   4173  wm-rh-S_precentral-Superior-part  0
678   4174  wm-rh-S_subcentral_ant  0
679   4175  wm-rh-S_subcentral_post  0
680   4176  wm-rh-S_suborbital  0
681   4177  wm-rh-S_subparietal  0
682   4178  wm-rh-S_supracingulate  0
683   4179  wm-rh-S_temporal_inferior  0
684   4180  wm-rh-S_temporal_superior  0
685   4181  wm-rh-S_temporal_transverse  0

Reporting on 685 segmentations
#--------------------------------------------
#@# Intensity Normalization2 Thu Jun 21 01:46:26 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_normalize -mask brainmask.mgz -aseg aseg.mgz norm.mgz brain.mgz 

using segmentation for initial intensity normalization
reading from norm.mgz...
normalizing image...
removing outliers in the aseg WM...
20664 control points removed
building Voronoi diagram...
performing soap bubble smoothing...
building Voronoi diagram...
performing soap bubble smoothing...
3d normalization pass 1 of 2
building Voronoi diagram...
performing soap bubble smoothing...
3d normalization pass 2 of 2
building Voronoi diagram...
performing soap bubble smoothing...
writing output to brain.mgz
3D bias adjustment took 5 minutes and 21 seconds.
using MR volume brainmask.mgz to mask input volume...
using wm threshold 103.0 for removing exterior voxels
white matter peak found at 110
gm peak at 71 (71), valley at 58 (58)
csf peak at 36, setting threshold to 59
white matter peak found at 110
gm peak at 55 (55), valley at 14 (14)
csf peak at 28, setting threshold to 46
#--------------------------------------------
#@# Mask BFS Thu Jun 21 01:51:55 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_mask -T 5 brain.mgz brainmask.mgz brain.finalsurfs.mgz 

threshold mask volume at 5
writing masked volume to brain.finalsurfs.mgz...
#--------------------------------------------
#@# WM Segmentation Thu Jun 21 01:52:01 EDT 2007

 mri_segment brain.mgz wm.seg.mgz 

doing initial intensity segmentation...
using local statistics to label ambiguous voxels...
computing class statistics for intensity windows...
WM: 110.0 +- 7.1 [80.0 --> 125.0]
GM: 62.1 +- 13.3 [30.0 --> 96.0]
setting bottom of white matter range to 75.5
setting top of gray matter range to 88.8
doing initial intensity segmentation...
using local statistics to label ambiguous voxels...
using local geometry to label remaining ambiguous voxels...

reclassifying voxels using Gaussian border classifier...

removing voxels with positive offset direction...
smoothing T1 volume with sigma = 0.250
removing 1-dimensional structures...
20532 sparsely connected voxels removed...
thickening thin strands....
20 segments, 8904 filled
41771 bright non-wm voxels segmented.
17482 diagonally connected voxels added...
white matter segmentation took 6.1 minutes
writing output to wm.seg.mgz...

 mri_edit_wm_with_aseg -keep-in wm.seg.mgz brain.mgz aseg.mgz wm.asegedit.mgz 

preserving editing changes in input volume...
auto filling took 0.66 minutes
reading wm segmentation from wm.seg.mgz...
0 voxels added to wm to prevent paths from MTL structures to cortex
14106 additional wm voxels added
0 additional wm voxels added
SEG EDIT: 41230 voxels turned on, 707910 voxels turned off.
propagating editing to output volume from wm.seg.mgz
115,126,128 old 0   new 0
115,126,128 old 0   new 0
writing edited volume to wm.asegedit.mgz....

 mri_pretess wm.asegedit.mgz wm norm.mgz wm.mgz 


Iteration Number : 1
pass   1 (xy+):  92 found -  92 modified     |    TOTAL:  92
pass   2 (xy+):   0 found -  92 modified     |    TOTAL:  92
pass   1 (xy-):  59 found -  59 modified     |    TOTAL: 151
pass   2 (xy-):   0 found -  59 modified     |    TOTAL: 151
pass   1 (yz+):  99 found -  99 modified     |    TOTAL: 250
pass   2 (yz+):   0 found -  99 modified     |    TOTAL: 250
pass   1 (yz-):  62 found -  62 modified     |    TOTAL: 312
pass   2 (yz-):   0 found -  62 modified     |    TOTAL: 312
pass   1 (xz+):  78 found -  78 modified     |    TOTAL: 390
pass   2 (xz+):   0 found -  78 modified     |    TOTAL: 390
pass   1 (xz-):  94 found -  94 modified     |    TOTAL: 484
pass   2 (xz-):   0 found -  94 modified     |    TOTAL: 484
Iteration Number : 1
pass   1 (+++):  33 found -  33 modified     |    TOTAL:  33
pass   2 (+++):   0 found -  33 modified     |    TOTAL:  33
pass   1 (+++):  27 found -  27 modified     |    TOTAL:  60
pass   2 (+++):   0 found -  27 modified     |    TOTAL:  60
pass   1 (+++):  25 found -  25 modified     |    TOTAL:  85
pass   2 (+++):   0 found -  25 modified     |    TOTAL:  85
pass   1 (+++):  33 found -  33 modified     |    TOTAL: 118
pass   2 (+++):   0 found -  33 modified     |    TOTAL: 118
Iteration Number : 1
pass   1 (++):  58 found -  58 modified     |    TOTAL:  58
pass   2 (++):   0 found -  58 modified     |    TOTAL:  58
pass   1 (+-):  56 found -  56 modified     |    TOTAL: 114
pass   2 (+-):   0 found -  56 modified     |    TOTAL: 114
pass   1 (--):  55 found -  55 modified     |    TOTAL: 169
pass   2 (--):   0 found -  55 modified     |    TOTAL: 169
pass   1 (-+):  75 found -  75 modified     |    TOTAL: 244
pass   2 (-+):   0 found -  75 modified     |    TOTAL: 244
Iteration Number : 2
pass   1 (xy+):   9 found -   9 modified     |    TOTAL:   9
pass   2 (xy+):   0 found -   9 modified     |    TOTAL:   9
pass   1 (xy-):   8 found -   8 modified     |    TOTAL:  17
pass   2 (xy-):   0 found -   8 modified     |    TOTAL:  17
pass   1 (yz+):   9 found -   9 modified     |    TOTAL:  26
pass   2 (yz+):   0 found -   9 modified     |    TOTAL:  26
pass   1 (yz-):   9 found -   9 modified     |    TOTAL:  35
pass   2 (yz-):   0 found -   9 modified     |    TOTAL:  35
pass   1 (xz+):   5 found -   5 modified     |    TOTAL:  40
pass   2 (xz+):   0 found -   5 modified     |    TOTAL:  40
pass   1 (xz-):   8 found -   8 modified     |    TOTAL:  48
pass   2 (xz-):   0 found -   8 modified     |    TOTAL:  48
Iteration Number : 2
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   4 found -   4 modified     |    TOTAL:   4
pass   2 (+++):   0 found -   4 modified     |    TOTAL:   4
pass   1 (+++):   2 found -   2 modified     |    TOTAL:   6
pass   2 (+++):   0 found -   2 modified     |    TOTAL:   6
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   6
Iteration Number : 2
pass   1 (++):   9 found -   9 modified     |    TOTAL:   9
pass   2 (++):   0 found -   9 modified     |    TOTAL:   9
pass   1 (+-):   2 found -   2 modified     |    TOTAL:  11
pass   2 (+-):   0 found -   2 modified     |    TOTAL:  11
pass   1 (--):   4 found -   4 modified     |    TOTAL:  15
pass   2 (--):   0 found -   4 modified     |    TOTAL:  15
pass   1 (-+):   4 found -   4 modified     |    TOTAL:  19
pass   2 (-+):   0 found -   4 modified     |    TOTAL:  19
Iteration Number : 3
pass   1 (xy+):   0 found -   0 modified     |    TOTAL:   0
pass   1 (xy-):   1 found -   1 modified     |    TOTAL:   1
pass   2 (xy-):   0 found -   1 modified     |    TOTAL:   1
pass   1 (yz+):   0 found -   0 modified     |    TOTAL:   1
pass   1 (yz-):   1 found -   1 modified     |    TOTAL:   2
pass   2 (yz-):   0 found -   1 modified     |    TOTAL:   2
pass   1 (xz+):   0 found -   0 modified     |    TOTAL:   2
pass   1 (xz-):   1 found -   1 modified     |    TOTAL:   3
pass   2 (xz-):   0 found -   1 modified     |    TOTAL:   3
Iteration Number : 3
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
Iteration Number : 3
pass   1 (++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+-):   0 found -   0 modified     |    TOTAL:   0
pass   1 (--):   0 found -   0 modified     |    TOTAL:   0
pass   1 (-+):   0 found -   0 modified     |    TOTAL:   0
Iteration Number : 4
pass   1 (xy+):   0 found -   0 modified     |    TOTAL:   0
pass   1 (xy-):   0 found -   0 modified     |    TOTAL:   0
pass   1 (yz+):   0 found -   0 modified     |    TOTAL:   0
pass   1 (yz-):   0 found -   0 modified     |    TOTAL:   0
pass   1 (xz+):   0 found -   0 modified     |    TOTAL:   0
pass   1 (xz-):   0 found -   0 modified     |    TOTAL:   0
Iteration Number : 4
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+++):   0 found -   0 modified     |    TOTAL:   0
Iteration Number : 4
pass   1 (++):   0 found -   0 modified     |    TOTAL:   0
pass   1 (+-):   0 found -   0 modified     |    TOTAL:   0
pass   1 (--):   0 found -   0 modified     |    TOTAL:   0
pass   1 (-+):   0 found -   0 modified     |    TOTAL:   0

Total Number of Modified Voxels = 922 (out of 76571: 1.204111)
binarizing input wm segmentation...
Ambiguous edge configurations... 

mri_pretess done

#--------------------------------------------
#@# Fill Thu Jun 21 01:58:57 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/mri

 mri_fill -a ../scripts/ponscc.cut.log -xform transforms/talairach.lta -segmentation aseg.mgz wm.mgz filled.mgz 

logging cutting plane coordinates to ../scripts/ponscc.cut.log...
INFO: Using transforms/talairach.lta and its offset for Talairach volume ...
using segmentation aseg.mgz...
reading input volume...done.
searching for cutting planes...voxel to talairach voxel transform
 1.809   0.747   0.019  -255.835;
-0.880   2.166  -0.530  -151.473;
-0.037   0.071   1.524  -163.836;
 0.000   0.000   0.000   1.000;
voxel to talairach voxel transform
 1.809   0.747   0.019  -255.835;
-0.880   2.166  -0.530  -151.473;
-0.037   0.071   1.524  -163.836;
 0.000   0.000   0.000   1.000;
Looking for area (min, max) = (350, 1400)
area[0] = 4335 (min = 350, max = 1400), aspect = 0.77 (min = 0.10, max = 0.75)
need search nearby
using seed (115, 138, 154), TAL = (13.0, 26.0, -10.0)
talairach voxel to voxel transform
 0.473  -0.161  -0.062   86.519;
 0.193   0.391   0.133   130.449;
 0.003  -0.022   0.648   103.543;
 0.000   0.000   0.000   1.000;
done.
writing output to filled.mgz...
filling took 1.6 minutes
INFO: Modifying dst c_(r,a,s), using the transform dst
INFO: Modifying dst c_(r,a,s), using the transform dst
reading segmented volume aseg.mgz...
INFO: Modifying dst c_(r,a,s), using the transform dst
segmentation indicates cc at (115,  138,  154) --> (13.0, 26.0, -10.0)
INFO: Modifying dst c_(r,a,s), using the transform dst
talairach cc position changed to (13.00, 26.00, -10.00)
ERASING BRAINSTEMseed_search_size = 9, min_neighbors = 5
search rh wm seed point around talairach space:(31.00, 26.00, -10.00) SRC: (100.65, 223.66, 200.59)
search lh wm seed point around talairach space (-5.00, 26.00, -10.00), SRC: (117.69, 230.61, 200.68)
compute mri_fill using aseg
Erasing Brain Stem and Cerebellum ...
Define left and right masks using aseg:
Building Voronoi diagram ...
Using the Voronoi diagram to separate WM into two hemispheres ...
Find the largest connected component for each hemisphere ...
#--------------------------------------------
#@# Tessellate lh Thu Jun 21 02:00:35 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/scripts

 mri_tessellate ../mri/filled.mgz 255 ../surf/lh.orig.nofix 

$Id: mri_tessellate.c,v 1.27.2.1 2006/05/02 21:41:01 nicks Exp $
  $Id: mrisurf.c,v 1.441.2.5 2006/11/29 20:28:16 nicks Exp $
slice 130: 697 vertices, 762 faces
slice 140: 2659 vertices, 2721 faces
slice 150: 4357 vertices, 4453 faces
slice 160: 7210 vertices, 7386 faces
slice 170: 12462 vertices, 12719 faces
slice 180: 17570 vertices, 17881 faces
slice 190: 22286 vertices, 22641 faces
slice 200: 27066 vertices, 27516 faces
slice 210: 31897 vertices, 32376 faces
slice 220: 36004 vertices, 36483 faces
slice 230: 37379 vertices, 37768 faces
slice 240: 37379 vertices, 37768 faces
slice 250: 37379 vertices, 37768 faces
using the conformed surface RAS to save vertex points...
writing ../surf/lh.orig.nofix
using ras2vox matrix:
-1.000   0.000   0.000   128.000;
 0.000   0.000   1.000  -128.000;
 0.000  -1.000   0.000   128.000;
 0.000   0.000   0.000   1.000;
#--------------------------------------------
#@# Smooth1 lh Thu Jun 21 02:01:58 EDT 2007

 mris_smooth -nw ../surf/lh.orig.nofix ../surf/lh.smoothwm.nofix 

/home/richfrye/Desktop/MRI/sub089/sub089mri/scripts
smoothing surface tessellation for 10 iterations...
smoothing complete - recomputing first and second fundamental forms...
#--------------------------------------------
#@# Inflation1 lh Thu Jun 21 02:02:00 EDT 2007

 mris_inflate -no-save-sulc ../surf/lh.smoothwm.nofix ../surf/lh.inflated.nofix 

/home/richfrye/Desktop/MRI/sub089/sub089mri/scripts
avg radius = 28.2 mm, total surface area = 18328 mm^2
writing inflated surface to ../surf/lh.inflated.nofix
inflation took 0.4 minutes
Not saving sulc
step 000: RMS=0.262 (target=0.015)   step 005: RMS=0.217 (target=0.015)   step 010: RMS=0.197 (target=0.015)   step 015: RMS=0.192 (target=0.015)   step 020: RMS=0.191 (target=0.015)   step 025: RMS=0.190 (target=0.015)   step 030: RMS=0.189 (target=0.015)   step 035: RMS=0.189 (target=0.015)   step 040: RMS=0.188 (target=0.015)   step 045: RMS=0.189 (target=0.015)   step 050: RMS=0.188 (target=0.015)   step 055: RMS=0.189 (target=0.015)   step 060: RMS=0.189 (target=0.015)   
inflation complete.
Not saving sulc
#--------------------------------------------
#@# QSphere lh Thu Jun 21 02:02:22 EDT 2007
/home/richfrye/Desktop/MRI/sub089/sub089mri/scripts

 mris_sphere -q ../surf/lh.inflated.nofix ../surf/lh.qsphere.nofix 

doing quick spherical unfolding.
$Id: mris_sphere.c,v 1.38 2006/02/08 00:17:40 nicks Exp $
  $Id: mrisurf.c,v 1.441.2.5 2006/11/29 20:28:16 nicks Exp $
reading original vertex positions...
unfolding cortex into spherical form...
tol=1.0e+00, host=local, nav=0, nbrs=1, l_spring_norm=1.000, l_sphere=0.025, l_convex=1.000
mom=0.90, dt=0.90

surface projected - minimizing metric distortion...
vertex spacing 1.66 +- 1.37 (0.00-->21.11) (max @ vno 6955 --> 6956)
face area 0.23 +- 0.43 (-4.23-->5.63)
tol=1.0e-01, host=local, nav=32, nbrs=1, l_nlarea=1.000
using quadratic fit line minimization

tol=1.8e+01, host=local, nav=32, nbrs=1, l_nlarea=1.000
using quadratic fit line minimization

301: dt: 5.61, sse: 1774.8 (0.265, 46.6, 0.738), neg: 15668 (%3.850:%7.66), avgs: 32
302: dt: 16.40, sse: 1627.4 (0.268, 46.6, 0.741), neg: 15922 (%2.851:%7.36), avgs: 32
303: dt: 9.29, sse: 1568.9 (0.266, 46.5, 0.740), neg: 15931 (%2.452:%7.27), avgs: 32
304: dt: 8.43, sse: 1544.5 (0.267, 46.4, 0.742), neg: 15879 (%2.258:%7.06), avgs: 32
305: dt: 0.29, sse: 1543.4 (0.267, 46.4, 0.741), neg: 15876 (%2.250:%7.05), avgs: 32
306: dt: 0.29, sse: 1541.4 (0.267, 46.4, 0.741), neg: 15874 (%2.236:%7.04), avgs: 32
307: dt: 0.29, sse: 1538.7 (0.267, 46.4, 0.741), neg: 15897 (%2.217:%7.05), avgs: 32
308: dt: 0.29, sse: 1535.4 (0.267, 46.4, 0.741), neg: 15895 (%2.194:%7.06), avgs: 32
309: dt: 0.29, sse: 1532.1 (0.267, 46.4, 0.741), neg: 15905 (%2.168:%7.05), avgs: 32
310: dt: 0.29, sse: 1529.3 (0.267, 46.4, 0.741), neg: 15911 (%2.145:%7.04), avgs: 32
vertex spacing 1.66 +- 1.40 (0.00-->20.57) (max @ vno 6955 --> 6956)
face area 0.24 +- 0.44 (-0.77-->5.81)
311: dt: 0.29, sse: 1526.4 (0.267, 46.4, 0.741), neg: 15936 (%2.121:%7.02), avgs: 32
312: dt: 0.29, sse: 1524.6 (0.267, 46.4, 0.741), neg: 15938 (%2.102:%7.02), avgs: 32
313: dt: 0.29, sse: 1521.9 (0.267, 46.4, 0.742), neg: 15918 (%2.079:%7.00), avgs: 32
314: dt: 0.29, sse: 1519.0 (0.267, 46.4, 0.742), neg: 15897 (%2.053:%6.98), avgs: 32
315: dt: 8.31, sse: 1508.4 (0.268, 46.4, 0.743), neg: 15826 (%1.961:%6.79), avgs: 32
vertex spacing 1.66 +- 1.40 (0.00-->20.62) (max @ vno 6955 --> 6956)
face area 0.24 +- 0.44 (-1.02-->5.83)
316: dt: 0.49, sse: 1499.8 (0.268, 46.4, 0.743), neg: 15869 (%1.899:%6.80), avgs: 8
317: dt: 0.15, sse: 1495.0 (0.268, 46.4, 0.743), neg: 15872 (%1.861:%6.76), avgs: 8
318: dt: 0.15, sse: 1489.0 (0.268, 46.4, 0.743), neg: 15907 (%1.815:%6.75), avgs: 8
319: dt: 0.15, sse: 1483.4 (0.269, 46.4, 0.743), neg: 15929 (%1.768:%6.71), avgs: 8
320: dt: 0.15, sse: 1477.7 (0.269, 46.4, 0.743), neg: 15956 (%1.721:%6.69), avgs: 8
vertex spacing 1.66 +- 1.40 (0.00-->20.61) (max @ vno 6955 --> 6956)
face area 0.24 +- 0.44 (-0.58-->5.85)
321: dt: 0.15, sse: 1471.9 (0.269, 46.4, 0.744), neg: 15974 (%1.673:%6.66), avgs: 8
322: dt: 0.15, sse: 1466.9 (0.269, 46.5, 0.744), neg: 16002 (%1.631:%6.61), avgs: 8
323: dt: 0.15, sse: 1462.4 (0.270, 46.5, 0.744), neg: 16021 (%1.591:%6.54), avgs: 8
324: dt: 0.15, sse: 1458.3 (0.270, 46.5, 0.745), neg: 16084 (%1.555:%6.52), avgs: 8
325: dt: 0.15, sse: 1454.4 (0.270, 46.5, 0.745), neg: 16099 (%1.522:%6.45), avgs: 8
326: dt: 0.15, sse: 1452.2 (0.270, 46.5, 0.745), neg: 16124 (%1.502:%6.41), avgs: 8
327: dt: 0.93, sse: 1449.1 (0.270, 46.5, 0.745), neg: 16139 (%1.474:%6.40), avgs: 8
vertex spacing 1.66 +- 1.41 (0.00-->20.60) (max @ vno 6955 --> 6956)
face area 0.24 +- 0.44 (-0.45-->5.87)
328: dt: 0.02, sse: 1447.4 (0.270, 46.5, 0.745), neg: 16174 (%1.463:%6.43), avgs: 2
329: dt: 0.09, sse: 1448.0 (0.271, 46.6, 0.746), neg: 16206 (%1.457:%6.36), avgs: 2
330: dt: 0.09, sse: 1454.7 (0.271, 46.6, 0.746), neg: 16241 (%1.481:%6.34), avgs: 2
vertex spacing 1.66 +- 1.41 (0.00-->20.77) (max @ vno 6955 --> 6956)
face area 0.24 +- 0.44 (-2.16-->5.87)
331: dt: 0.09, sse: 1443.2 (0.271, 46.7, 0.746), neg: 16292 (%1.408:%6.28), avgs: 2
332: dt: 0.09, sse: 1438.2 (0.271, 46.7, 0.747), neg: 16302 (%1.370:%6.20), avgs: 2
333: dt: 0.09, sse: 1434.3 (0.271, 46.8, 0.747), neg: 16354 (%1.339:%6.19), avgs: 2
334: dt: 0.09, sse: 1430.9 (0.272, 46.8, 0.747), neg: 16409 (%1.312:%6.14), avgs: 2
335: dt: 0.09, sse: 1428.5 (0.272, 46.9, 0.748), neg: 16459 (%1.293:%6.08), avgs: 2
336: dt: 0.09, sse: 1432.8 (0.272, 46.9, 0.748), neg: 16474 (%1.315:%6.01), avgs: 2
337: dt: 0.09, sse: 1436.3 (0.272, 47.0, 0.749), neg: 16558 (%1.332:%5.96), avgs: 2
338: dt: 0.09, sse: 1428.2 (0.272, 47.0, 0.749), neg: 16574 (%1.288:%5.86), avgs: 2
339: dt: 0.01, sse: 1426.7 (0.272, 47.0, 0.749), neg: 16583 (%1.280:%5.87), avgs: 2
vertex spacing 1.67 +- 1.42 (0.00-->19.80) (max @ vno 6965 --> 6966)
face area 0.24 +- 0.44 (-0.81-->5.89)
340: dt: 0.00, sse: 1425.9 (0.272, 47.0, 0.749), neg: 16610 (%1.274:%5.86), avgs: 0
vertex spacing 1.67 +- 1.42 (0.00-->19.80) (max @ vno 6965 --> 6966)
face area 0.24 +- 0.44 (-0.81-->5.89)
vertex spacing 1.67 +- 1.42 (0.00-->19.80) (max @ vno 6965 --> 6966)
face area 0.24 +- 0.44 (-0.81-->5.89)
tol=1.8e+01, host=local, nav=32, nbrs=1, l_nlarea=1.000
using quadratic fit line minimization

341: dt: 28.93, sse: 399.6 (0.273, 46.9, 0.749), neg: 16451 (%1.163:%5.95), avgs: 32
342: dt: 2.47, sse: 394.4 (0.273, 46.8, 0.749), neg: 16413 (%1.132:%5.84), avgs: 32
343: dt: 0.29, sse: 405.8 (0.274, 46.9, 0.749), neg: 16458 (%1.194:%5.81), avgs: 32
344: dt: 0.29, sse: 420.6 (0.275, 46.9, 0.750), neg: 16455 (%1.275:%5.71), avgs: 32
345: dt: 0.29, sse: 434.8 (0.278, 47.0, 0.750), neg: 16406 (%1.354:%5.61), avgs: 32
346: dt: 0.29, sse: 447.2 (0.280, 47.1, 0.751), neg: 16443 (%1.430:%5.59), avgs: 32
347: dt: 0.29, sse: 457.2 (0.283, 47.1, 0.751), neg: 16471 (%1.497:%5.58), avgs: 32
348: dt: 0.29, sse: 465.1 (0.286, 47.2, 0.752), neg: 16480 (%1.555:%5.58), avgs: 32
349: dt: 0.29, sse: 472.8 (0.288, 47.2, 0.752), neg: 16515 (%1.607:%5.61), avgs: 32
350: dt: 0.29, sse: 479.8 (0.289, 47.3, 0.753), neg: 16537 (%1.651:%5.61), avgs: 32
vertex spacing 1.67 +- 1.43 (0.00-->18.42) (max @ vno 6965 --> 6966)
face area 0.24 +- 0.45 (-8.24-->8.99)
351: dt: 0.29, sse: 486.0 (0.290, 47.4, 0.753), neg: 16574 (%1.685:%5.63), avgs: 32
352: dt: 0.29, sse: 491.6 (0.290, 47.4, 0.754), neg: 16622 (%1.712:%5.66), avgs: 32
353: dt: 6.97, sse: 489.4 (0.281, 47.6, 0.756), neg: 16760 (%1.653:%5.91), avgs: 32
vertex spacing 1.68 +- 1.44 (0.00-->14.61) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-3.62-->5.88)
354: dt: 3.19, sse: 467.9 (0.279, 47.5, 0.758), neg: 16699 (%1.524:%5.82), avgs: 8
355: dt: 2.00, sse: 441.4 (0.278, 47.5, 0.756), neg: 16695 (%1.377:%5.75), avgs: 8
356: dt: 1.88, sse: 433.4 (0.278, 47.4, 0.757), neg: 16683 (%1.329:%5.80), avgs: 8
357: dt: 1.97, sse: 426.3 (0.278, 47.4, 0.757), neg: 16690 (%1.286:%5.76), avgs: 8
358: dt: 0.15, sse: 431.4 (0.278, 47.4, 0.757), neg: 16730 (%1.314:%5.80), avgs: 8
359: dt: 0.15, sse: 436.5 (0.278, 47.5, 0.758), neg: 16730 (%1.341:%5.78), avgs: 8
360: dt: 0.15, sse: 439.9 (0.279, 47.5, 0.759), neg: 16715 (%1.360:%5.78), avgs: 8
vertex spacing 1.68 +- 1.45 (0.00-->14.36) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-1.22-->5.89)
361: dt: 0.15, sse: 442.8 (0.279, 47.5, 0.759), neg: 16734 (%1.375:%5.78), avgs: 8
362: dt: 0.15, sse: 443.5 (0.279, 47.6, 0.760), neg: 16735 (%1.379:%5.82), avgs: 8
363: dt: 0.15, sse: 442.9 (0.280, 47.6, 0.760), neg: 16742 (%1.375:%5.81), avgs: 8
364: dt: 0.15, sse: 441.0 (0.280, 47.6, 0.760), neg: 16766 (%1.365:%5.85), avgs: 8
365: dt: 0.15, sse: 438.1 (0.280, 47.6, 0.760), neg: 16799 (%1.348:%5.88), avgs: 8
366: dt: 0.15, sse: 435.4 (0.280, 47.6, 0.760), neg: 16792 (%1.333:%5.87), avgs: 8
367: dt: 0.15, sse: 433.9 (0.280, 47.6, 0.760), neg: 16778 (%1.323:%5.85), avgs: 8
368: dt: 1.70, sse: 429.0 (0.280, 47.6, 0.760), neg: 16688 (%1.294:%5.77), avgs: 8
vertex spacing 1.68 +- 1.45 (0.00-->14.15) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-1.43-->5.90)
inflating to sphere (rms error < 2.00)
000: dt: 0.0000, rms radial error=160.707, avgs=0
005: dt: 0.9000, rms radial error=160.466, avgs=0
010: dt: 0.9000, rms radial error=159.940, avgs=0
015: dt: 0.9000, rms radial error=159.260, avgs=0
020: dt: 0.9000, rms radial error=158.484, avgs=0
025: dt: 0.9000, rms radial error=157.646, avgs=0
030: dt: 0.9000, rms radial error=156.760, avgs=0
035: dt: 0.9000, rms radial error=155.858, avgs=0
040: dt: 0.9000, rms radial error=154.947, avgs=0
045: dt: 0.9000, rms radial error=154.035, avgs=0
050: dt: 0.9000, rms radial error=153.126, avgs=0
055: dt: 0.9000, rms radial error=152.260, avgs=0
060: dt: 0.9000, rms radial error=151.412, avgs=0
065: dt: 0.9000, rms radial error=150.568, avgs=0
070: dt: 0.9000, rms radial error=149.730, avgs=0
075: dt: 0.9000, rms radial error=148.896, avgs=0
080: dt: 0.9000, rms radial error=148.067, avgs=0
085: dt: 0.9000, rms radial error=147.243, avgs=0
090: dt: 0.9000, rms radial error=146.424, avgs=0
095: dt: 0.9000, rms radial error=145.609, avgs=0
100: dt: 0.9000, rms radial error=144.799, avgs=0
105: dt: 0.9000, rms radial error=143.994, avgs=0
110: dt: 0.9000, rms radial error=143.194, avgs=0
115: dt: 0.9000, rms radial error=142.398, avgs=0
120: dt: 0.9000, rms radial error=141.607, avgs=0
125: dt: 0.9000, rms radial error=140.820, avgs=0
130: dt: 0.9000, rms radial error=140.039, avgs=0
135: dt: 0.9000, rms radial error=139.262, avgs=0
140: dt: 0.9000, rms radial error=138.490, avgs=0
145: dt: 0.9000, rms radial error=137.722, avgs=0
150: dt: 0.9000, rms radial error=136.958, avgs=0
155: dt: 0.9000, rms radial error=136.200, avgs=0
160: dt: 0.9000, rms radial error=135.445, avgs=0
165: dt: 0.9000, rms radial error=134.695, avgs=0
170: dt: 0.9000, rms radial error=133.950, avgs=0
175: dt: 0.9000, rms radial error=133.209, avgs=0
180: dt: 0.9000, rms radial error=132.472, avgs=0
185: dt: 0.9000, rms radial error=131.740, avgs=0
190: dt: 0.9000, rms radial error=131.012, avgs=0
195: dt: 0.9000, rms radial error=130.301, avgs=0
200: dt: 0.9000, rms radial error=129.594, avgs=0
205: dt: 0.9000, rms radial error=128.892, avgs=0
210: dt: 0.9000, rms radial error=128.194, avgs=0
215: dt: 0.9000, rms radial error=127.499, avgs=0
220: dt: 0.9000, rms radial error=126.809, avgs=0
225: dt: 0.9000, rms radial error=126.123, avgs=0
230: dt: 0.9000, rms radial error=125.441, avgs=0
235: dt: 0.9000, rms radial error=124.763, avgs=0
240: dt: 0.9000, rms radial error=124.088, avgs=0
245: dt: 0.9000, rms radial error=123.418, avgs=0
250: dt: 0.9000, rms radial error=122.751, avgs=0
255: dt: 0.9000, rms radial error=122.089, avgs=0
260: dt: 0.9000, rms radial error=121.430, avgs=0
265: dt: 0.9000, rms radial error=120.775, avgs=0
270: dt: 0.9000, rms radial error=120.124, avgs=0
275: dt: 0.9000, rms radial error=119.477, avgs=0
280: dt: 0.9000, rms radial error=118.834, avgs=0
285: dt: 0.9000, rms radial error=118.194, avgs=0
290: dt: 0.9000, rms radial error=117.558, avgs=0
295: dt: 0.9000, rms radial error=116.926, avgs=0
300: dt: 0.9000, rms radial error=116.297, avgs=0

spherical inflation complete.
epoch 1 (K=20.0), pass 1, starting sse = 1894.37
integrating with navgs=32 and tol=3.200e+00
taking momentum steps...
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=8 and tol=1.671e+00
taking momentum steps...
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=2 and tol=9.648e-01
taking momentum steps...
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=0 and tol=5.570e-01
pass 1 complete, delta sse/iter = 0.25/40 = 0.00618
epoch 2 (K=80.0), pass 1, starting sse = 415.15
integrating with navgs=32 and tol=3.200e+00
taking momentum steps...
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=8 and tol=1.671e+00
taking momentum steps...
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=8 and tol=1.671e+00
integrat369: dt: 0.44, sse: 410.2 (0.280, 47.5, 0.760), neg: 16652 (%1.180:%5.60), avgs: 2
370: dt: 0.15, sse: 400.2 (0.280, 47.5, 0.760), neg: 16667 (%1.121:%5.57), avgs: 2
vertex spacing 1.68 +- 1.45 (0.00-->14.10) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-0.65-->5.91)
371: dt: 0.08, sse: 395.2 (0.280, 47.5, 0.760), neg: 16723 (%1.092:%5.64), avgs: 2
372: dt: 0.09, sse: 393.4 (0.280, 47.4, 0.760), neg: 16783 (%1.080:%5.71), avgs: 2
373: dt: 0.09, sse: 390.7 (0.280, 47.4, 0.760), neg: 16761 (%1.062:%5.71), avgs: 2
374: dt: 0.09, sse: 387.6 (0.281, 47.4, 0.760), neg: 16761 (%1.042:%5.76), avgs: 2
375: dt: 0.09, sse: 383.4 (0.281, 47.4, 0.760), neg: 16737 (%1.014:%5.73), avgs: 2
376: dt: 0.09, sse: 378.5 (0.281, 47.4, 0.760), neg: 16721 (%0.981:%5.71), avgs: 2
377: dt: 0.09, sse: 373.6 (0.281, 47.3, 0.760), neg: 16679 (%0.949:%5.71), avgs: 2
378: dt: 0.09, sse: 369.4 (0.281, 47.3, 0.760), neg: 16671 (%0.920:%5.76), avgs: 2
379: dt: 0.09, sse: 365.6 (0.281, 47.2, 0.760), neg: 16605 (%0.893:%5.73), avgs: 2
380: dt: 0.09, sse: 361.2 (0.282, 47.2, 0.760), neg: 16543 (%0.864:%5.66), avgs: 2
vertex spacing 1.68 +- 1.45 (0.00-->13.84) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-0.48-->5.93)
381: dt: 0.09, sse: 357.3 (0.282, 47.2, 0.760), neg: 16488 (%0.838:%5.61), avgs: 2
382: dt: 0.09, sse: 354.0 (0.282, 47.1, 0.760), neg: 16438 (%0.816:%5.61), avgs: 2
383: dt: 0.11, sse: 352.7 (0.282, 47.1, 0.760), neg: 16508 (%0.809:%5.66), avgs: 2
vertex spacing 1.68 +- 1.45 (0.00-->13.77) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-0.50-->5.93)
384: dt: 0.03, sse: 352.2 (0.282, 47.1, 0.760), neg: 16578 (%0.804:%5.66), avgs: 0
vertex spacing 1.68 +- 1.45 (0.00-->13.77) (max @ vno 24318 --> 24759)
face area 0.24 +- 0.45 (-0.42-->5.93)
tol=1.8e+01, host=local, nav=32, nbrs=1, l_nlarea=1.000
using quadratic fit line minimization

385: dt: 11.60, sse: 167.6 (0.281, 47.0, 0.760), neg: 16255 (%0.765:%5.46), avgs: 32
386: dt: 5.07, sse: 165.4 (0.281, 47.0, 0.760), neg: 16284 (%0.753:%5.38), avgs: 32
387: dt: 0.29, sse: 164.5 (0.281, 47.0, 0.760), neg: 16221 (%0.748:%5.38), avgs: 32
388: dt: 0.29, sse: 165.6 (0.282, 46.9, 0.760), neg: 16078 (%0.753:%5.34), avgs: 32
389: dt: 0.29, sse: 167.3 (0.282, 46.9, 0.760), neg: 15982 (%0.762:%5.32), avgs: 32
390: dt: 0.29, sse: 169.5 (0.282, 46.9, 0.760), neg: 15911 (%0.775:%5.33), avgs: 32
vertex spacing 1.68 +- 1.45 (0.00-->13.71) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-1.00-->5.93)
391: dt: 0.29, sse: 171.8 (0.282, 46.9, 0.761), neg: 15853 (%0.787:%5.32), avgs: 32
392: dt: 0.29, sse: 173.9 (0.282, 46.9, 0.761), neg: 15823 (%0.799:%5.34), avgs: 32
393: dt: 0.29, sse: 176.0 (0.282, 46.9, 0.761), neg: 15805 (%0.811:%5.38), avgs: 32
394: dt: 0.29, sse: 178.2 (0.282, 46.9, 0.761), neg: 15782 (%0.823:%5.39), avgs: 32
395: dt: 0.29, sse: 180.0 (0.282, 46.9, 0.761), neg: 15745 (%0.834:%5.39), avgs: 32
396: dt: 0.29, sse: 181.5 (0.282, 46.9, 0.761), neg: 15710 (%0.842:%5.37), avgs: 32
397: dt: 3.79, sse: 180.4 (0.282, 46.8, 0.761), neg: 15732 (%0.835:%5.34), avgs: 32
vertex spacing 1.68 +- 1.46 (0.00-->13.89) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-1.72-->5.93)
398: dt: 1.67, sse: 173.2 (0.282, 46.8, 0.761), neg: 15744 (%0.793:%5.31), avgs: 8
399: dt: 1.35, sse: 169.6 (0.282, 46.8, 0.761), neg: 15739 (%0.773:%5.28), avgs: 8
400: dt: 0.67, sse: 166.4 (0.282, 46.8, 0.761), neg: 15748 (%0.754:%5.25), avgs: 8
vertex spacing 1.68 +- 1.46 (0.00-->14.04) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.77-->5.93)
401: dt: 1.59, sse: 164.3 (0.282, 46.8, 0.761), neg: 15726 (%0.742:%5.33), avgs: 8
402: dt: 0.15, sse: 165.1 (0.282, 46.8, 0.761), neg: 15712 (%0.746:%5.33), avgs: 8
403: dt: 0.15, sse: 165.4 (0.282, 46.8, 0.761), neg: 15705 (%0.748:%5.32), avgs: 8
404: dt: 0.15, sse: 166.1 (0.282, 46.8, 0.761), neg: 15708 (%0.751:%5.33), avgs: 8
405: dt: 0.15, sse: 166.8 (0.282, 46.8, 0.761), neg: 15701 (%0.755:%5.32), avgs: 8
406: dt: 0.15, sse: 167.1 (0.282, 46.8, 0.761), neg: 15712 (%0.756:%5.35), avgs: 8
407: dt: 0.15, sse: 166.8 (0.282, 46.8, 0.761), neg: 15710 (%0.754:%5.36), avgs: 8
408: dt: 0.15, sse: 166.2 (0.283, 46.8, 0.761), neg: 15684 (%0.751:%5.35), avgs: 8
409: dt: 0.15, sse: 165.5 (0.283, 46.8, 0.761), neg: 15683 (%0.747:%5.38), avgs: 8
410: dt: 0.15, sse: 165.1 (0.283, 46.7, 0.761), neg: 15694 (%0.744:%5.40), avgs: 8
vertex spacing 1.68 +- 1.46 (0.00-->14.58) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.84-->5.94)
411: dt: 0.15, sse: 164.7 (0.283, 46.7, 0.761), neg: 15701 (%0.741:%5.42), avgs: 8
412: dt: 0.77, sse: 162.9 (0.283, 46.7, 0.761), neg: 15697 (%0.731:%5.36), avgs: 8
vertex spacing 1.68 +- 1.46 (0.00-->14.71) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.73-->5.94)
413: dt: 0.25, sse: 153.5 (0.283, 46.7, 0.761), neg: 15550 (%0.678:%5.12), avgs: 2
414: dt: 0.07, sse: 150.5 (0.283, 46.7, 0.761), neg: 15601 (%0.661:%5.09), avgs: 2
415: dt: 0.08, sse: 148.5 (0.283, 46.7, 0.761), neg: 15642 (%0.650:%5.10), avgs: 2
416: dt: 0.09, sse: 147.2 (0.283, 46.7, 0.761), neg: 15694 (%0.642:%5.17), avgs: 2
417: dt: 0.09, sse: 145.5 (0.283, 46.7, 0.761), neg: 15677 (%0.632:%5.16), avgs: 2
418: dt: 0.09, sse: 143.2 (0.283, 46.7, 0.761), neg: 15730 (%0.619:%5.20), avgs: 2
419: dt: 0.09, sse: 141.1 (0.283, 46.7, 0.761), neg: 15704 (%0.605:%5.18), avgs: 2
420: dt: 0.09, sse: 138.0 (0.283, 46.6, 0.761), neg: 15710 (%0.587:%5.19), avgs: 2
vertex spacing 1.68 +- 1.46 (0.00-->15.05) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.41-->5.95)
421: dt: 0.09, sse: 134.6 (0.283, 46.6, 0.761), neg: 15641 (%0.567:%5.10), avgs: 2
422: dt: 0.09, sse: 132.0 (0.283, 46.6, 0.761), neg: 15616 (%0.552:%5.09), avgs: 2
423: dt: 0.09, sse: 130.8 (0.283, 46.5, 0.761), neg: 15566 (%0.543:%5.03), avgs: 2
424: dt: 0.09, sse: 129.1 (0.283, 46.5, 0.761), neg: 15542 (%0.533:%4.99), avgs: 2
425: dt: 0.09, sse: 127.5 (0.284, 46.5, 0.761), neg: 15569 (%0.523:%5.06), avgs: 2
426: dt: 0.09, sse: 125.8 (0.284, 46.4, 0.762), neg: 15508 (%0.512:%5.01), avgs: 2
427: dt: 0.09, sse: 124.3 (0.284, 46.4, 0.761), neg: 15569 (%0.504:%5.01), avgs: 2
428: dt: 0.03, sse: 124.1 (0.284, 46.4, 0.762), neg: 15621 (%0.502:%5.00), avgs: 2
vertex spacing 1.68 +- 1.46 (0.00-->15.40) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.28-->5.95)
429: dt: 0.03, sse: 122.5 (0.284, 46.4, 0.762), neg: 15609 (%0.493:%5.03), avgs: 0
430: dt: 0.00, sse: 122.5 (0.284, 46.4, 0.762), neg: 15651 (%0.493:%5.07), avgs: 0
vertex spacing 1.68 +- 1.46 (0.00-->15.41) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.26-->5.95)
vertex spacing 1.68 +- 1.46 (0.00-->15.41) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.26-->5.95)
tol=1.8e+01, host=local, nav=32, nbrs=1, l_nlarea=1.000
using quadratic fit line minimization

431: dt: 5.68, sse: 92.4 (0.284, 46.4, 0.762), neg: 15520 (%0.485:%5.10), avgs: 32
432: dt: 0.29, sse: 92.5 (0.284, 46.4, 0.762), neg: 15450 (%0.485:%5.07), avgs: 32
433: dt: 0.29, sse: 93.7 (0.284, 46.4, 0.762), neg: 15385 (%0.493:%5.07), avgs: 32
434: dt: 0.29, sse: 95.9 (0.284, 46.4, 0.762), neg: 15344 (%0.505:%5.13), avgs: 32
435: dt: 0.29, sse: 98.3 (0.284, 46.4, 0.762), neg: 15305 (%0.519:%5.16), avgs: 32
436: dt: 0.29, sse: 100.6 (0.284, 46.4, 0.762), neg: 15311 (%0.533:%5.22), avgs: 32
437: dt: 0.29, sse: 102.9 (0.284, 46.4, 0.762), neg: 15306 (%0.546:%5.27), avgs: 32
438: dt: 0.29, sse: 104.9 (0.284, 46.4, 0.762), neg: 15274 (%0.559:%5.28), avgs: 32
439: dt: 0.29, sse: 106.6 (0.284, 46.4, 0.762), neg: 15270 (%0.569:%5.29), avgs: 32
440: dt: 0.29, sse: 108.4 (0.284, 46.4, 0.762), neg: 15260 (%0.580:%5.30), avgs: 32
vertex spacing 1.68 +- 1.46 (0.00-->15.72) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-1.25-->5.95)
441: dt: 0.29, sse: 110.2 (0.284, 46.4, 0.762), neg: 15257 (%0.591:%5.33), avgs: 32
442: dt: 4.01, sse: 109.8 (0.284, 46.4, 0.762), neg: 15270 (%0.588:%5.32), avgs: 32
vertex spacing 1.68 +- 1.46 (0.00-->15.71) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-1.09-->5.95)
443: dt: 1.52, sse: 104.1 (0.284, 46.4, 0.762), neg: 15265 (%0.553:%5.20), avgs: 8
444: dt: 0.82, sse: 101.2 (0.284, 46.4, 0.762), neg: 15250 (%0.535:%5.17), avgs: 8
445: dt: 0.20, sse: 100.4 (0.284, 46.4, 0.762), neg: 15266 (%0.530:%5.16), avgs: 8
446: dt: 0.15, sse: 101.5 (0.284, 46.4, 0.762), neg: 15292 (%0.537:%5.20), avgs: 8
447: dt: 0.15, sse: 102.3 (0.284, 46.4, 0.762), neg: 15316 (%0.541:%5.23), avgs: 8
448: dt: 0.15, sse: 102.6 (0.284, 46.4, 0.762), neg: 15298 (%0.543:%5.23), avgs: 8
449: dt: 0.15, sse: 102.8 (0.284, 46.3, 0.762), neg: 15296 (%0.544:%5.27), avgs: 8
450: dt: 0.15, sse: 102.7 (0.284, 46.3, 0.762), neg: 15298 (%0.543:%5.32), avgs: 8
vertex spacing 1.68 +- 1.46 (0.00-->15.81) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.53-->5.95)
451: dt: 0.15, sse: 102.1 (0.284, 46.3, 0.762), neg: 15284 (%0.540:%5.33), avgs: 8
452: dt: 0.15, sse: 101.6 (0.284, 46.3, 0.762), neg: 15259 (%0.537:%5.31), avgs: 8
453: dt: 0.15, sse: 100.9 (0.284, 46.3, 0.762), neg: 15220 (%0.533:%5.30), avgs: 8
454: dt: 0.15, sse: 100.6 (0.284, 46.3, 0.762), neg: 15207 (%0.531:%5.31), avgs: 8
455: dt: 0.15, sse: 100.3 (0.284, 46.3, 0.762), neg: 15179 (%0.529:%5.28), avgs: 8
456: dt: 0.99, sse: 99.6 (0.284, 46.3, 0.762), neg: 15193 (%0.524:%5.23), avgs: 8
vertex spacing 1.68 +- 1.46 (0.00-->15.91) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.43-->5.95)
457: dt: 0.13, sse: 95.1 (0.284, 46.3, 0.762), neg: 15126 (%0.498:%5.07), avgs: 2
458: dt: 0.07, sse: 92.4 (0.284, 46.2, 0.762), neg: 15164 (%0.483:%5.03), avgs: 2
459: dt: 0.08, sse: 90.7 (0.284, 46.2, 0.762), neg: 15185 (%0.473:%5.05), avgs: 2
460: dt: 0.09, sse: 89.3 (0.284, 46.2, 0.762), neg: 15184 (%0.465:%5.06), avgs: 2
vertex spacing 1.68 +- 1.46 (0.00-->15.94) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.33-->5.95)
461: dt: 0.05, sse: 88.3 (0.284, 46.2, 0.762), neg: 15230 (%0.460:%5.06), avgs: 2
462: dt: 0.08, sse: 87.4 (0.284, 46.2, 0.762), neg: 15293 (%0.455:%5.11), avgs: 2
463: dt: 0.05, sse: 86.5 (0.284, 46.2, 0.762), neg: 15294 (%0.449:%5.08), avgs: 2
464: dt: 0.08, sse: 85.6 (0.284, 46.2, 0.762), neg: 15309 (%0.444:%5.12), avgs: 2
465: dt: 0.05, sse: 85.1 (0.284, 46.2, 0.762), neg: 15335 (%0.442:%5.14), avgs: 2
466: dt: 0.09, sse: 83.9 (0.284, 46.2, 0.762), neg: 15254 (%0.435:%5.09), avgs: 2
467: dt: 0.09, sse: 82.4 (0.284, 46.2, 0.762), neg: 15266 (%0.427:%5.05), avgs: 2
468: dt: 0.09, sse: 81.5 (0.285, 46.1, 0.762), neg: 15243 (%0.422:%5.07), avgs: 2
469: dt: 0.09, sse: 79.8 (0.285, 46.1, 0.762), neg: 15198 (%0.412:%5.03), avgs: 2
470: dt: 0.09, sse: 78.3 (0.285, 46.1, 0.762), neg: 15204 (%0.404:%5.06), avgs: 2
vertex spacing 1.68 +- 1.46 (0.00-->16.15) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.25-->5.96)
471: dt: 0.09, sse: 77.5 (0.285, 46.1, 0.762), neg: 15165 (%0.399:%5.06), avgs: 2
472: dt: 0.09, sse: 76.2 (0.285, 46.0, 0.762), neg: 15168 (%0.392:%5.09), avgs: 2
473: dt: 0.09, sse: 74.7 (0.285, 46.0, 0.762), neg: 15159 (%0.383:%5.10), avgs: 2
474: dt: 0.09, sse: 74.9 (0.285, 46.0, 0.762), neg: 15144 (%0.385:%5.08), avgs: 2
475: dt: 0.09, sse: 73.1 (0.285, 46.0, 0.762), neg: 15154 (%0.375:%5.11), avgs: 2
476: dt: 0.08, sse: 72.0 (0.285, 46.0, 0.762), neg: 15188 (%0.369:%5.00), avgs: 2
477: dt: 0.05, sse: 71.2 (0.285, 46.0, 0.762), neg: 15204 (%0.364:%4.92), avgs: 2
478: dt: 0.00, sse: 71.2 (0.285, 46.0, 0.762), neg: 15204 (%0.364:%4.92), avgs: 2
vertex spacing 1.68 +- 1.46 (0.00-->16.48) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.24-->5.96)
479: dt: 0.01, sse: 71.0 (0.285, 46.0, 0.762), neg: 15276 (%0.363:%4.99), avgs: 0
vertex spacing 1.68 +- 1.46 (0.00-->16.48) (max @ vno 7780 --> 7797)
face area 0.24 +- 0.45 (-0.24-->5.96)
479: dt: 0.00, sse: 1384.1 (0.285, 46.0, 0.762), neg: 15276 (%0.363:%4.99), avgs: 32
writing spherical brain to ../surf/lh.qsphere.nofix
spherical transformation took 0.07 hours
ing with navgs=2 and tol=9.648e-01
taking momentum steps...
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=0 and tol=5.570e-01
pass 1 complete, delta sse/iter = 0.15/44 = 0.00344
epoch 3 (K=320.0), pass 1, starting sse = 174.60
integrating with navgs=32 and tol=3.200e+00
taking momentum steps...
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=8 and tol=1.671e+00
taking momentum steps...
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=2 and tol=9.648e-01
taking momentum steps...
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=0 and tol=5.570e-01
pass 1 complete, delta sse/iter = 0.30/46 = 0.00649
epoch 4 (K=1280.0), pass 1, starting sse = 93.82
integrating with navgs=32 and tol=3.200e+00
taking momentum steps...
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=32 and tol=3.200e+00
integrating with navgs=8 and tol=1.671e+00
taking momentum steps...
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=8 and tol=1.671e+00
integrating with navgs=2 and tol=9.648e-01
taking momentum steps...
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=2 and tol=9.648e-01
integrating with navgs=0 and tol=5.570e-01
pass 1 complete, delta sse/iter = 0.24/49 = 0.00495
final distance error %100000.00
optimization complete.
unfolding took 0.06 hours
#--------------------------------------------
#@# Fix Topology lh Thu Jun 21 02:06:48 EDT 2007

 cp ../surf/lh.orig.nofix ../surf/lh.orig 


 cp ../surf/lh.inflated.nofix ../surf/lh.inflated 

/home/richfrye/Desktop/MRI/sub089/sub089mri/scripts

 mris_fix_topology -mgz -sphere qsphere.nofix -ga sub089mri lh 

reading spherical homeomorphism from 'qsphere.nofix'
using genetic algorithm with optimized parameters

*************************************************************
Topology Correction Parameters
retessellation mode:           genetic search
number of patches/generation : 10
number of generations :        10
surface mri loglikelihood coefficient :         1.0
volume mri loglikelihood coefficient :          10.0
normal dot loglikelihood coefficient :          1.0
quadratic curvature loglikelihood coefficient : 1.0
volume resolution :                             2
eliminate vertices during search :              1
initial patch selection :                       1
select all defect vertices :                    0
ordering dependant retessellation:              0
use precomputed edge table :                    0
smooth retessellated patch :                    2
match retessellated patch :                     1
verbose mode :                                  0

*************************************************************
INFO: assuming .mgz format
$Id: mris_fix_topology.c,v 1.41.2.1 2007/01/05 16:58:33 nicks Exp $
  $Id: mrisurf.c,v 1.441.2.5 2006/11/29 20:28:16 nicks Exp $
before topology correction, eno=-305 (nv=37379, nf=75536, ne=113220, g=153)
using quasi-homeomorphic spherical map to tessellate cortical surface...

Correction of the Topology
Finding true center and radius of Spherical Surface...done
Surface centered at (0,0,0) with radius 100.0 in 13 iterations
marking ambiguous vertices...
35584 ambiguous faces found in tessellation
segmenting defects...
35 defects found, arbitrating ambiguous regions...
analyzing neighboring defects...
      -merging segment 4 into 0
      -merging segment 5 into 0
      -merging segment 31 into 0
      -merging segment 22 into 21
31 defects to be corrected 
0 vertices coincident
vertex 25113: degenerate normal
vertex 26622: degenerate normal
vertex 28630: degenerate normal
vertex 31525: degenerate normal
vertex 25113: degenerate normal
vertex 26622: degenerate normal
vertex 28630: degenerate normal
vertex 31525: degenerate normal
vertex 25113: degenerate normal
vertex 26622: degenerate normal
vertex 28630: degenerate normal
vertex 31525: degenerate normal
Computing Initial Surface Statistics
      -face       loglikelihood: -8.8994  (-4.4497)
      -vertex     loglikelihood: -8.2166  (-4.1083)
      -normal dot loglikelihood: -3.4360  (-3.4360)
      -quad curv  loglikelihood: -6.4670  (-3.2335)
      Total Loglikelihood : -27.0190

CORRECTING DEFECT 0 (vertices=16283, convex hull=3920)
normal vector of length zero at vertex 18644 with 0 faces
vertex 18644 has 0 face
reading input surface /home/richfrye/Desktop/MRI/sub089/sub089mri/surf/lh.qsphere.nofix...
reading brain volume from brain...
reading wm segmentation from wm...
Linux localhost.localdomain 2.6.9-55.ELsmp #1 SMP Fri Apr 20 16:36:54 EDT 2007 x86_64 x86_64 x86_64 GNU/Linux

recon-all exited with ERRORS at Thu Jun 21 12:07:42 EDT 2007




