Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

BUG: is float leads to Description: itk::ERROR: ImageToImageFilter(0x7fb3b2307ac0): Inputs do not occupy the same physical space! #31

Closed
stnava opened this issue Jul 1, 2013 · 5 comments
Closed

BUG: is float leads to Description: itk::ERROR: ImageToImageFilter(0x7fb3b2307ac0): Inputs do not occupy the same physical space! #31

stnava opened this issue Jul 1, 2013 · 5 comments

Comments

@stnava
Copy link
Member

stnava commented Jul 1, 2013

I suspect this error:

Description: itk::ERROR: ImageToImageFilter(0x7fb3b2307ac0): Inputs do not occupy the same physical space!

is due to some tolerance issue that is coded as double but when we run in float mode does not activate ---- this is on OSX

will email example later ....
@Chaircrusher
Copy link
Collaborator

Do you know which filter is actually causing this problem? itk::ImageToImageFilter::VerifyInputInformation calls vnl_vector::is_equal and vnl_matrix::is_equal. These do comparisons within a tolerance set by ImageToImageFilter::m_CoordinateTolerance, which by default is 1.0e-6.

That's actually a fairly liberal tolerance; if the filters are throwing this exception, they either need to implement their own VerifyInputInformation, or there is a real problem in the congruence of the physical spaces.

@chrisgorgo
Copy link
Contributor

I'm experiencing similar error without using float (sticking with double):

kalifornien:/scr/adenauer1/ants$ antsRegistration --dimensionality 3 --output transform --use-histogram-matching 1 --initial-moving-transform ../ants_SyN_masked/transform0Warp.nii.gz --transform BSplineSyN[0.1,26,0,3] --metric CC[/usr/share/fsl/data/standard/MNI152_T1_1mm_brain.nii.gz, /scr/adenauer1/ants/linear.nii,1,4] --convergence 100x70x50x20 --smoothing-sigmas 3x2x1x0 --shrink-factors 6x4x2x1
Using double precision for computations.
=============================================================================
The composite transform is comprised of the following transforms (in order): 
  1. ../ants_SyN_masked/transform0Warp.nii.gz (type = DisplacementFieldTransform)
=============================================================================
  number of levels = 4
  fixed image: /usr/share/fsl/data/standard/MNI152_T1_1mm_brain.nii.gz
  moving image: /scr/adenauer1/ants/linear.nii
Dimension = 3
Number of stages = 1
Use Histogram Matching true
Winsorize Image Intensities false
Lower Quantile = 0
Upper Quantile = 1
Stage 1 State
   Metric = CC
     Fixed Image = Image (0x2f72690)
  RTTI typeinfo:   itk::Image<double, 3u>
  Reference Count: 2
  Modified Time: 1544
  Debug: Off
  Object Name: 
  Observers: 
    none
  Source: (none)
  Source output name: (none)
  Release Data: Off
  Data Released: False
  Global Release Data: Off
  PipelineMTime: 0
  UpdateMTime: 1383
  RealTimeStamp: 0 seconds 
  LargestPossibleRegion: 
    Dimension: 3
    Index: [0, 0, 0]
    Size: [182, 218, 182]
  BufferedRegion: 
    Dimension: 3
    Index: [0, 0, 0]
    Size: [182, 218, 182]
  RequestedRegion: 
    Dimension: 3
    Index: [0, 0, 0]
    Size: [182, 218, 182]
  Spacing: [1, 1, 1]
  Origin: [-90, 126, -72]
  Direction: 
1 0 0
0 -1 0
0 0 1

  IndexToPointMatrix: 
1 0 0
0 -1 0
0 0 1

  PointToIndexMatrix: 
1 0 0
0 -1 0
0 0 1

  Inverse Direction: 
1 0 0
0 -1 0
0 0 1

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

     Moving Image = Image (0x2f6b880)
  RTTI typeinfo:   itk::Image<double, 3u>
  Reference Count: 2
  Modified Time: 1545
  Debug: Off
  Object Name: 
  Observers: 
    none
  Source: (none)
  Source output name: (none)
  Release Data: Off
  Data Released: False
  Global Release Data: Off
  PipelineMTime: 0
  UpdateMTime: 1542
  RealTimeStamp: 0 seconds 
  LargestPossibleRegion: 
    Dimension: 3
    Index: [0, 0, 0]
    Size: [182, 218, 182]
  BufferedRegion: 
    Dimension: 3
    Index: [0, 0, 0]
    Size: [182, 218, 182]
  RequestedRegion: 
    Dimension: 3
    Index: [0, 0, 0]
    Size: [182, 218, 182]
  Spacing: [1, 1, 1]
  Origin: [-90, 126, -72]
  Direction: 
1 0 0
0 -1 0
0 0 1

  IndexToPointMatrix: 
1 0 0
0 -1 0
0 0 1

  PointToIndexMatrix: 
1 0 0
0 -1 0
0 0 1

  Inverse Direction: 
1 0 0
0 -1 0
0 0 1

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

     Weighting = 1
     Sampling Strategy = none
     NumberOfBins = 1
     Radius = 4
     Sampling percentage  = 1
   Transform = BSplineSyN
     Gradient Step = 0.1
     Update Field Sigma (voxel space) = 0
     Total Field Sigma (voxel space) = 0
     Update Field Time Sigma = 0
     Total Field Time Sigma  = 0
     Number of Time Indices = 0
     Number of Time Point Samples = 0
Registration using 1 total stages.

Stage 0
  iterations = 100x70x50x20
  convergence threshold = 1e-06
  convergence window size = 10
  number of levels = 4
  Shrink factors (level 1 out of 4): [6, 6, 6]
  Shrink factors (level 2 out of 4): [4, 4, 4]
  Shrink factors (level 3 out of 4): [2, 2, 2]
  Shrink factors (level 4 out of 4): [1, 1, 1]
  smoothing sigmas per level: [3, 2, 1, 0]
  preprocessing:  histogram matching the images
  using the CC metric (radius = 4, weight = 1)
  Using default NONE metricSamplingStrategy 

*** Running B-spline SyN registration (updateMeshSizeAtBaseLevel = [7, 9, 7], totalMeshSizeAtBaseLevel = [0, 0, 0]) ***

  Current level = 1 of 4
    number of iterations = 100
    shrink factors = [6, 6, 6]
    smoothing sigmas = 3 vox
    required fixed parameters = [30, 36, 30, -86.5, 122.5, -68.5, 6, 6, 6, 1, 0, 0, 0, -1, 0, 0, 0, 1]
XDIAGNOSTIC,Iteration,metricValue,convergenceValue,ITERATION_TIME_INDEX,SINCE_LAST
WDIAGNOSTIC,     1, -7.807433469000e-01, 1.797693134862e+308, 1.5243e+00, 1.5242e+00, 
WDIAGNOSTIC,     2, -7.824985387568e-01, 1.797693134862e+308, 1.8216e+00, 2.9727e-01, 
WDIAGNOSTIC,     3, -7.841857720773e-01, 1.797693134862e+308, 2.1036e+00, 2.8203e-01, 
WDIAGNOSTIC,     4, -7.856564459842e-01, 1.797693134862e+308, 2.3934e+00, 2.8973e-01, 
WDIAGNOSTIC,     5, -7.869333696374e-01, 1.797693134862e+308, 2.6794e+00, 2.8609e-01, 
WDIAGNOSTIC,     6, -7.882013679583e-01, 1.797693134862e+308, 2.9641e+00, 2.8463e-01, 
WDIAGNOSTIC,     7, -7.893345491175e-01, 1.797693134862e+308, 3.2617e+00, 2.9763e-01, 
WDIAGNOSTIC,     8, -7.904293713367e-01, 1.797693134862e+308, 3.5489e+00, 2.8721e-01, 
WDIAGNOSTIC,     9, -7.914969838079e-01, 1.797693134862e+308, 3.8484e+00, 2.9945e-01, 
WDIAGNOSTIC,    10, -7.924752733705e-01, 9.355939664866e-04, 4.1486e+00, 3.0029e-01, 
WDIAGNOSTIC,    11, -7.932677091932e-01, 7.852530111584e-04, 4.4488e+00, 3.0016e-01, 
WDIAGNOSTIC,    12, -7.939260235063e-01, 6.613280173870e-04, 4.7370e+00, 2.8824e-01, 
WDIAGNOSTIC,    13, -7.945001158981e-01, 5.588707971683e-04, 5.0343e+00, 2.9727e-01, 
WDIAGNOSTIC,    14, -7.949769189912e-01, 4.698964055225e-04, 5.3279e+00, 2.9360e-01, 
WDIAGNOSTIC,    15, -7.953406409922e-01, 3.887885582880e-04, 5.6113e+00, 2.8342e-01, 
WDIAGNOSTIC,    16, -7.956315513545e-01, 3.169466665355e-04, 5.9004e+00, 2.8912e-01, 
WDIAGNOSTIC,    17, -7.958655253009e-01, 2.530694620223e-04, 6.1777e+00, 2.7724e-01, 
WDIAGNOSTIC,    18, -7.961183190603e-01, 1.995557973640e-04, 6.4790e+00, 3.0136e-01, 
WDIAGNOSTIC,    19, -7.963042287376e-01, 1.565922737426e-04, 6.7681e+00, 2.8905e-01, 
WDIAGNOSTIC,    20, -7.964269020391e-01, 1.229664638404e-04, 7.0637e+00, 2.9564e-01, 
WDIAGNOSTIC,    21, -7.965822142226e-01, 9.718985744912e-05, 7.3508e+00, 2.8702e-01, 
WDIAGNOSTIC,    22, -7.967126438470e-01, 7.741133670639e-05, 7.6594e+00, 3.0860e-01, 
WDIAGNOSTIC,    23, -7.967315991381e-01, 6.123245785656e-05, 7.9413e+00, 2.8192e-01, 
WDIAGNOSTIC,    24, -7.967673652349e-01, 4.828848305207e-05, 8.2198e+00, 2.7856e-01, 
WDIAGNOSTIC,    25, -7.967621822854e-01, 3.696300890311e-05, 8.5056e+00, 2.8581e-01, 
WDIAGNOSTIC,    26, -7.968399381158e-01, 2.809164724515e-05, 8.7936e+00, 2.8799e-01, 
WDIAGNOSTIC,    27, -7.968674175305e-01, 2.062091300102e-05, 9.0828e+00, 2.8916e-01, 
WDIAGNOSTIC,    28, -7.969471255377e-01, 1.598424331081e-05, 9.3831e+00, 3.0030e-01, 
WDIAGNOSTIC,    29, -7.968810302086e-01, 1.164290206259e-05, 9.6752e+00, 2.9207e-01, 
WDIAGNOSTIC,    30, -7.968807804740e-01, 7.971321558629e-06, 9.9543e+00, 2.7914e-01, 
WDIAGNOSTIC,    31, -7.969195385805e-01, 5.994473765272e-06, 1.0240e+01, 2.8600e-01, 
WDIAGNOSTIC,    32, -7.968526781827e-01, 4.258712957159e-06, 1.0529e+01, 2.8895e-01, 
WDIAGNOSTIC,    33, -7.968278830906e-01, 2.187522083766e-06, 1.0822e+01, 2.9261e-01, 
WDIAGNOSTIC,    34, -7.970773077322e-01, 3.076729032575e-06, 1.1109e+01, 2.8732e-01, 
WDIAGNOSTIC,    35, -7.970083541697e-01, 2.782383086496e-06, 1.1398e+01, 2.8851e-01, 
WDIAGNOSTIC,    36, -7.969717798456e-01, 2.692520773277e-06, 1.1679e+01, 2.8162e-01, 
WDIAGNOSTIC,    37, -7.968661692902e-01, 1.627447761283e-06, 1.1976e+01, 2.9648e-01, 
WDIAGNOSTIC,    38, -7.967909625666e-01, 5.510098876477e-07, 1.2260e+01, 2.8421e-01, 
  Current level = 2 of 4
    number of iterations = 70
    shrink factors = [4, 4, 4]
    smoothing sigmas = 2.0000e+00 vox
    required fixed parameters = [45, 54, 45, -87.5, 123.5, -69.5, 4, 4, 4, 1, 0, 0, 0, -1, 0, 0, 0, 1]
XDIAGNOSTIC,Iteration,metricValue,convergenceValue,ITERATION_TIME_INDEX,SINCE_LAST
WDIAGNOSTIC,     1, -7.710102284037e-01, 1.797693134862e+308, 1.4258e+01, 1.9976e+00, 
WDIAGNOSTIC,     2, -7.748677692146e-01, 1.797693134862e+308, 1.5167e+01, 9.0911e-01, 
WDIAGNOSTIC,     3, -7.790675177968e-01, 1.797693134862e+308, 1.6052e+01, 8.8559e-01, 
WDIAGNOSTIC,     4, -7.828436140554e-01, 1.797693134862e+308, 1.6928e+01, 8.7569e-01, 
WDIAGNOSTIC,     5, -7.861370005545e-01, 1.797693134862e+308, 1.7822e+01, 8.9385e-01, 
WDIAGNOSTIC,     6, -7.891419008214e-01, 1.797693134862e+308, 1.8715e+01, 8.9327e-01, 
WDIAGNOSTIC,     7, -7.917679605200e-01, 1.797693134862e+308, 1.9633e+01, 9.1807e-01, 
WDIAGNOSTIC,     8, -7.939367221719e-01, 1.797693134862e+308, 2.0544e+01, 9.1040e-01, 
WDIAGNOSTIC,     9, -7.956841422919e-01, 1.797693134862e+308, 2.1405e+01, 8.6149e-01, 
WDIAGNOSTIC,    10, -7.974221352482e-01, 2.158541482106e-03, 2.2296e+01, 8.9058e-01, 
WDIAGNOSTIC,    11, -7.989145370382e-01, 1.745216598845e-03, 2.3189e+01, 8.9273e-01, 
WDIAGNOSTIC,    12, -8.002209503154e-01, 1.398050511259e-03, 2.4057e+01, 8.6854e-01, 
WDIAGNOSTIC,    13, -8.013537936133e-01, 1.124786932673e-03, 2.4954e+01, 8.9689e-01, 
WDIAGNOSTIC,    14, -8.021942130124e-01, 9.062730611692e-04, 2.5898e+01, 9.4390e-01, 
WDIAGNOSTIC,    15, -8.031869069772e-01, 7.360726443068e-04, 2.6805e+01, 9.0698e-01, 
WDIAGNOSTIC,    16, -8.040909472253e-01, 6.063921312430e-04, 2.7718e+01, 9.1321e-01, 
WDIAGNOSTIC,    17, -8.048837456096e-01, 5.075622385734e-04, 2.8596e+01, 8.7772e-01, 
WDIAGNOSTIC,    18, -8.055492482351e-01, 4.288776494007e-04, 2.9472e+01, 8.7625e-01, 
WDIAGNOSTIC,    19, -8.061735625521e-01, 3.630483432817e-04, 3.0346e+01, 8.7414e-01, 
WDIAGNOSTIC,    20, -8.068841969649e-01, 3.134600255625e-04, 3.1225e+01, 8.7917e-01, 
WDIAGNOSTIC,    21, -8.074851582061e-01, 2.745334682261e-04, 3.2119e+01, 8.9384e-01, 
WDIAGNOSTIC,    22, -8.080671272958e-01, 2.438535843800e-04, 3.3001e+01, 8.8229e-01, 
WDIAGNOSTIC,    23, -8.084029967312e-01, 2.156962731931e-04, 3.3875e+01, 8.7363e-01, 
WDIAGNOSTIC,    24, -8.087539381860e-01, 1.872314616134e-04, 3.4759e+01, 8.8393e-01, 
WDIAGNOSTIC,    25, -8.091845872457e-01, 1.635418239278e-04, 3.5666e+01, 9.0706e-01, 
WDIAGNOSTIC,    26, -8.095739608328e-01, 1.437545021658e-04, 3.6539e+01, 8.7296e-01, 
WDIAGNOSTIC,    27, -8.097974694316e-01, 1.251539448682e-04, 3.7418e+01, 8.7917e-01, 
WDIAGNOSTIC,    28, -8.099678403527e-01, 1.066927639434e-04, 3.8305e+01, 8.8722e-01, 
WDIAGNOSTIC,    29, -8.101935762963e-01, 8.978529453776e-05, 3.9220e+01, 9.1483e-01, 
WDIAGNOSTIC,    30, -8.103602769553e-01, 7.575727462081e-05, 4.0140e+01, 9.1991e-01, 
WDIAGNOSTIC,    31, -8.104553431146e-01, 6.329303514731e-05, 4.1038e+01, 8.9754e-01, 
WDIAGNOSTIC,    32, -8.106493964653e-01, 5.423259377901e-05, 4.1916e+01, 8.7833e-01, 
WDIAGNOSTIC,    33, -8.106557141115e-01, 4.422748606051e-05, 4.2793e+01, 8.7665e-01, 
WDIAGNOSTIC,    34, -8.106511401274e-01, 3.419112563673e-05, 4.3704e+01, 9.1118e-01, 
WDIAGNOSTIC,    35, -8.107653344432e-01, 2.695630713010e-05, 4.4590e+01, 8.8636e-01, 
WDIAGNOSTIC,    36, -8.108188968806e-01, 2.187176057098e-05, 4.5505e+01, 9.1464e-01, 
WDIAGNOSTIC,    37, -8.108232996195e-01, 1.714638997002e-05, 4.6404e+01, 8.9933e-01, 
WDIAGNOSTIC,    38, -8.108057499749e-01, 1.246265634034e-05, 4.7282e+01, 8.7734e-01, 
WDIAGNOSTIC,    39, -8.107511603521e-01, 8.404178680942e-06, 4.8156e+01, 8.7409e-01, 
WDIAGNOSTIC,    40, -8.107411304122e-01, 5.137766230769e-06, 4.9038e+01, 8.8245e-01, 
WDIAGNOSTIC,    41, -8.106956188368e-01, 1.914029310726e-06, 4.9925e+01, 8.8734e-01, 
WDIAGNOSTIC,    42, -8.106148688378e-01, -3.588547017597e-07, 5.0828e+01, 9.0291e-01, 
  Current level = 3 of 4
    number of iterations = 50
    shrink factors = [2, 2, 2]
    smoothing sigmas = 1.0000e+00 vox
    required fixed parameters = [91, 109, 91, -89.5, 125.5, -71.5, 2, 2, 2, 1, 0, 0, 0, -1, 0, 0, 0, 1]
XDIAGNOSTIC,Iteration,metricValue,convergenceValue,ITERATION_TIME_INDEX,SINCE_LAST
WDIAGNOSTIC,     1, -8.254524930297e-01, 1.797693134862e+308, 5.8316e+01, 7.4881e+00, 
WDIAGNOSTIC,     2, -8.271440202857e-01, 1.797693134862e+308, 6.4644e+01, 6.3276e+00, 
WDIAGNOSTIC,     3, -8.291818601772e-01, 1.797693134862e+308, 7.1016e+01, 6.3716e+00, 
WDIAGNOSTIC,     4, -8.310663360106e-01, 1.797693134862e+308, 7.7363e+01, 6.3472e+00, 
WDIAGNOSTIC,     5, -8.328685915056e-01, 1.797693134862e+308, 8.3689e+01, 6.3264e+00, 
WDIAGNOSTIC,     6, -8.345501120066e-01, 1.797693134862e+308, 8.9965e+01, 6.2758e+00, 
WDIAGNOSTIC,     7, -8.360489253918e-01, 1.797693134862e+308, 9.6313e+01, 6.3482e+00, 
WDIAGNOSTIC,     8, -8.374133556764e-01, 1.797693134862e+308, 1.0266e+02, 6.3453e+00, 
WDIAGNOSTIC,     9, -8.386802132123e-01, 1.797693134862e+308, 1.0908e+02, 6.4217e+00, 
WDIAGNOSTIC,    10, -8.398920856294e-01, 1.122611395519e-03, 1.1548e+02, 6.3990e+00, 
WDIAGNOSTIC,    11, -8.410597081234e-01, 9.603956566840e-04, 1.2183e+02, 6.3484e+00, 
WDIAGNOSTIC,    12, -8.421831593850e-01, 8.180347431815e-04, 1.2812e+02, 6.2893e+00, 
WDIAGNOSTIC,    13, -8.432580487105e-01, 7.044274123090e-04, 1.3450e+02, 6.3879e+00, 
WDIAGNOSTIC,    14, -8.442026206026e-01, 6.110625730097e-04, 1.4080e+02, 6.2949e+00, 
WDIAGNOSTIC,    15, -8.450388868869e-01, 5.335137917391e-04, 1.4729e+02, 6.4919e+00, 
WDIAGNOSTIC,    16, -8.457367170852e-01, 4.667253037059e-04, 1.5366e+02, 6.3714e+00, 
WDIAGNOSTIC,    17, -8.463768167524e-01, 4.071060328244e-04, 1.6015e+02, 6.4838e+00, 
WDIAGNOSTIC,    18, -8.468993108669e-01, 3.516768789485e-04, 1.6642e+02, 6.2771e+00, 
WDIAGNOSTIC,    19, -8.473639735076e-01, 3.000843953057e-04, 1.7279e+02, 6.3613e+00, 
WDIAGNOSTIC,    20, -8.477838777484e-01, 2.531531040796e-04, 1.7917e+02, 6.3886e+00, 
WDIAGNOSTIC,    21, -8.481648435135e-01, 2.118345444264e-04, 1.8556e+02, 6.3835e+00, 
WDIAGNOSTIC,    22, -8.484945180769e-01, 1.766043046207e-04, 1.9198e+02, 6.4234e+00, 
WDIAGNOSTIC,    23, -8.488174270892e-01, 1.482759495936e-04, 1.9834e+02, 6.3564e+00, 
WDIAGNOSTIC,    24, -8.491111129889e-01, 1.256746494572e-04, 2.0488e+02, 6.5409e+00, 
WDIAGNOSTIC,    25, -8.493920570818e-01, 1.080394746984e-04, 2.1160e+02, 6.7217e+00, 
WDIAGNOSTIC,    26, -8.496632422145e-01, 9.405079338617e-05, 2.1795e+02, 6.3456e+00, 
WDIAGNOSTIC,    27, -8.499005703894e-01, 8.310976104373e-05, 2.2431e+02, 6.3630e+00, 
WDIAGNOSTIC,    28, -8.501213245293e-01, 7.398820744416e-05, 2.3072e+02, 6.4068e+00, 
WDIAGNOSTIC,    29, -8.503369038027e-01, 6.634679611450e-05, 2.3712e+02, 6.4059e+00, 
WDIAGNOSTIC,    30, -8.505702428371e-01, 6.021319301225e-05, 2.4344e+02, 6.3232e+00, 
WDIAGNOSTIC,    31, -8.508229788584e-01, 5.560131351919e-05, 2.4979e+02, 6.3487e+00, 
WDIAGNOSTIC,    32, -8.510332075172e-01, 5.162910861196e-05, 2.5618e+02, 6.3837e+00, 
WDIAGNOSTIC,    33, -8.512591540825e-01, 4.856723619606e-05, 2.6262e+02, 6.4452e+00, 
WDIAGNOSTIC,    34, -8.514919234180e-01, 4.624484709519e-05, 2.6913e+02, 6.5041e+00, 
WDIAGNOSTIC,    35, -8.516728539798e-01, 4.404631006546e-05, 2.7549e+02, 6.3596e+00, 
WDIAGNOSTIC,    36, -8.518927314653e-01, 4.233774731277e-05, 2.8187e+02, 6.3865e+00, 
WDIAGNOSTIC,    37, -8.521240988392e-01, 4.091753180960e-05, 2.8826e+02, 6.3918e+00, 
WDIAGNOSTIC,    38, -8.523469427743e-01, 3.956431783670e-05, 2.9460e+02, 6.3397e+00, 
WDIAGNOSTIC,    39, -8.525581328075e-01, 3.815374212134e-05, 3.0102e+02, 6.4195e+00, 
WDIAGNOSTIC,    40, -8.527624507681e-01, 3.682120160997e-05, 3.0736e+02, 6.3335e+00, 
WDIAGNOSTIC,    41, -8.530047797554e-01, 3.607475902256e-05, 3.1383e+02, 6.4719e+00, 
WDIAGNOSTIC,    42, -8.532446348761e-01, 3.547370417845e-05, 3.2023e+02, 6.4057e+00, 
WDIAGNOSTIC,    43, -8.534643431411e-01, 3.495074344054e-05, 3.2672e+02, 6.4811e+00, 
WDIAGNOSTIC,    44, -8.536713657969e-01, 3.445337158786e-05, 3.3311e+02, 6.3994e+00, 
WDIAGNOSTIC,    45, -8.538736381128e-01, 3.352059258878e-05, 3.3943e+02, 6.3186e+00, 
WDIAGNOSTIC,    46, -8.540668990877e-01, 3.240351867933e-05, 3.4581e+02, 6.3753e+00, 
WDIAGNOSTIC,    47, -8.542299340117e-01, 3.101639732365e-05, 3.5219e+02, 6.3805e+00, 
WDIAGNOSTIC,    48, -8.544112326109e-01, 2.952422260145e-05, 3.5850e+02, 6.3124e+00, 
WDIAGNOSTIC,    49, -8.545995655242e-01, 2.798153485686e-05, 3.6488e+02, 6.3825e+00, 
WDIAGNOSTIC,    50, -8.547830074586e-01, 2.639017336151e-05, 3.7128e+02, 6.3995e+00, 
  Current level = 4 of 4
    number of iterations = 20
    shrink factors = [1, 1, 1]
    smoothing sigmas = 0.0000e+00 vox
    required fixed parameters = [182, 218, 182, -90, 126, -72, 1, 1, 1, 1, 0, 0, 0, -1, 0, 0, 0, 1]
XDIAGNOSTIC,Iteration,metricValue,convergenceValue,ITERATION_TIME_INDEX,SINCE_LAST
WDIAGNOSTIC,     1, -8.629684084058e-01, 1.797693134862e+308, 4.2460e+02, 5.3306e+01, 
WDIAGNOSTIC,     2, -8.643027462962e-01, 1.797693134862e+308, 4.7860e+02, 5.4000e+01, 
WDIAGNOSTIC,     3, -8.657988355672e-01, 1.797693134862e+308, 5.3107e+02, 5.2470e+01, 
WDIAGNOSTIC,     4, -8.671380259551e-01, 1.797693134862e+308, 5.8455e+02, 5.3486e+01, 
WDIAGNOSTIC,     5, -8.683649332121e-01, 1.797693134862e+308, 6.3768e+02, 5.3130e+01, 
WDIAGNOSTIC,     6, -8.694751760729e-01, 1.797693134862e+308, 6.9047e+02, 5.2782e+01, 
WDIAGNOSTIC,     7, -8.704964622880e-01, 1.797693134862e+308, 7.4325e+02, 5.2783e+01, 
WDIAGNOSTIC,     8, -8.714456104943e-01, 1.797693134862e+308, 7.9637e+02, 5.3119e+01, 
WDIAGNOSTIC,     9, -8.723368589923e-01, 1.797693134862e+308, 8.4966e+02, 5.3294e+01, 
WDIAGNOSTIC,    10, -8.732077762028e-01, 7.521381350736e-04, 9.0310e+02, 5.3439e+01, 
WDIAGNOSTIC,    11, -8.740660652078e-01, 6.393356191088e-04, 9.5570e+02, 5.2599e+01, 
WDIAGNOSTIC,    12, -8.749568552338e-01, 5.478799309718e-04, 1.0094e+03, 5.3670e+01, 
WDIAGNOSTIC,    13, -8.758422844382e-01, 4.802240059706e-04, 1.0620e+03, 5.2674e+01, 
WDIAGNOSTIC,    14, -8.767362375139e-01, 4.301284816257e-04, 1.1160e+03, 5.3936e+01, 
WDIAGNOSTIC,    15, -8.776336798073e-01, 3.931621778037e-04, 1.1696e+03, 5.3655e+01, 
WDIAGNOSTIC,    16, -8.785285883767e-01, 3.651800232320e-04, 1.2232e+03, 5.3550e+01, 
WDIAGNOSTIC,    17, -8.793879497702e-01, 3.425142021480e-04, 1.2760e+03, 5.2851e+01, 
WDIAGNOSTIC,    18, -8.802080111855e-01, 3.223912193662e-04, 1.3286e+03, 5.2536e+01, 
WDIAGNOSTIC,    19, -8.809854655319e-01, 3.027984535805e-04, 1.3812e+03, 5.2665e+01, 
WDIAGNOSTIC,    20, -8.817130934347e-01, 2.827587968901e-04, 1.4345e+03, 5.3290e+01, 
  Elapsed time (stage 0): 1.4373e+03


Total elapsed time: 1.4373e+03
Exception Object caught: 

itk::ExceptionObject (0x32721e0)
Location: "unknown" 
File: /SCR/software/ants-build/ITKv4-install/include/ITK-4.6/itkImageToImageFilter.hxx
Line: 248
Description: itk::ERROR: ImageToImageFilter(0x2f724a0): Inputs do not occupy the same physical space! 
InputImage Origin: [-8.7500000e+01, 1.2350000e+02, -6.9500000e+01], InputImage_1 Origin: [-9.0000000e+01, 1.2600000e+02, -7.2000000e+01]
    Tolerance: 4.0000000e-06
InputImage Spacing: [4.0000000e+00, 4.0000000e+00, 4.0000000e+00], InputImage_1 Spacing: [1.0000000e+00, 1.0000000e+00, 1.0000000e+00]
    Tolerance: 4.0000000e-06

@ntustison
Copy link
Member

This was a common occurrence for me both inside and outside the registration framework (using B-spline SyN and original SyN) on my university's linux cluster. The same operations on the same data would not reproduce the exception on my Mac Pro. The exception thrown is in the ImageToImageFilter but there is also a related check in DisplacementFieldTransform which was also a constant problem for me. Finally, I simply changed the tolerance values in both places from 1e-6 to 1e-2 and haven't had an issue since. I couldn't think of an alternative solution (outside of unifying such tolerance values into a CMake variable) since these tolerance values are buried pretty deep.

@chrisgorgo
Copy link
Contributor

For me the problem was the initial transformation. When I estimated it using floats and used it as init without floats I got the error

@cookpa cookpa mentioned this issue Mar 7, 2014
Closed
@stnava
Copy link
Member Author

stnava commented Oct 8, 2014

Moving comments to: #74

@stnava stnava closed this as completed Oct 8, 2014
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
4 participants
@chrisgorgo @stnava @ntustison @Chaircrusher