[Insight-users] difference in MultiResImageRegistration2 example output
Matt McCormick
matt.mccormick at kitware.com
Sat Oct 29 18:10:57 EDT 2011
Hi Nain,
The documentation in the source code may be inaccurate.
There are unit tests to determine if the code is working correctly. Try running
ctest -R MultiResImageRegistration2 -V
In you build directory. Those results show what is expected.
Thanks,
Matt
On Sat, Oct 29, 2011 at 5:28 PM, Jagjeet Nain <jagjeet.nain at gmail.com> wrote:
> Hi,
> I am bit new to ITK.
> After compiling InsightToolkit-3.20.1 on my Linux suse 11.3(64 bit) box, i
> tried running example provided along with ITK.
> Some example run fine, but MultiResImageRegistration* examples are giving
> slightly different output then given in source file.
> MultiResImageRegistration2 should converge as per the below detail from
> MultiResImageRegistration2.cxx.
> The registration converges after $5$ iterations in the first
> // level, $7$ in the second level and $4$ in the third level. The final
> // results when printed as an array of parameters are
> //
> // \begin{verbatim}
> // [1.00164, 0.00147688, 0.00168372, 1.0027, 12.6296, 16.4768]
> -------------------------------
> but my executable takes lot of iterations for converging. I want to
> emphasize here, that no change has been made to the source file.
> Run was made like below:
>
> ./MultiResImageRegistration2
> ../../InsightToolkit-3.20.1/Examples/Data/BrainT1SliceBorder20.png
> ../../InsightToolkit-3.20.1/Examples/Data/BrainProtonDensitySliceShifted13x17y.png
> out.png
> and its output was
> -------------------------------------
> MultiResolution Level : 0
>
> 0 -1.87377 [0.999988, 5.26965e-05, 2.92263e-05, 1.00007, 2.82271,
> 29.23] 0
> 1 -1.13243 [1, 3.55154e-05, 2.09008e-05, 1.00008, 15.45, 27.1347] 1
> 2 -1.38838 [1.00004, 8.32104e-05, 5.34165e-05, 1.00011, 5.32325,
> 28.6531] 2
> 3 -1.19929 [1.00006, 4.77606e-07, -1.93134e-05, 1.00008, 2.79924,
> 18.7291] 3
> 4 -1.23267 [1.00004, 1.87947e-05, -8.89719e-06, 1.00009, 10.8252,
> 17.0881] 4
> 5 -1.7143 [1, 3.64921e-05, 8.89255e-06, 1.00009, 18.9275, 18.2971] 5
> 6 -1.44958 [1.00002, 2.56425e-05, -1.65816e-06, 1.0001, 12.9226,
> 20.922] 6
> 7 -1.63302 [1.00003, 1.55103e-05, -8.84177e-06, 1.00011, 13.8223,
> 15.7569] 7
> 8 -1.78001 [1.00004, 2.05547e-05, -9.7243e-06, 1.00011, 13.7179,
> 19.9499] 8
> 9 -1.67048 [1.00003, 3.21431e-05, -4.5119e-06, 1.0001, 12.1907,
> 16.9622] 9
> 10 -1.8265 [1.00005, 3.86097e-05, -9.16404e-07, 1.00011, 14.6889,
> 15.9801] 10
> 11 -1.74105 [1.00006, 3.81516e-05, 1.12975e-06, 1.00011, 13.1032,
> 17.4283] 11
> 12 -1.86948 [1.00006, 3.79044e-05, 2.25677e-06, 1.00012, 14.7698,
> 17.0109] 12
> 13 -1.76855 [1.00007, 3.59673e-05, 1.24242e-06, 1.00012, 13.4169,
> 16.7686] 13
> 14 -1.85574 [1.00008, 3.95096e-05, 3.39406e-06, 1.00013, 12.1107,
> 17.196] 14
> 15 -1.82394 [1.00008, 4.14684e-05, 3.53581e-06, 1.00013, 13.1784,
> 16.9337] 15
> 16 -1.87515 [1.00008, 4.48528e-05, 5.76028e-06, 1.00013, 12.4321,
> 17.3992] 16
> 17 -1.84227 [1.00009, 4.56554e-05, 5.19095e-06, 1.00014, 13.1355,
> 17.4183] 17
> 18 -1.87113 [1.0001, 4.4566e-05, 4.72429e-06, 1.00015, 13.7229,
> 17.0309] 18
> 19 -1.84301 [1.0001, 4.53401e-05, 4.7149e-06, 1.00015, 13.2549,
> 17.3438] 19
> 20 -1.87119 [1.0001, 4.90078e-05, 6.5465e-06, 1.00015, 12.7891,
> 17.0278] 20
> 21 -1.86387 [1.0001, 4.93968e-05, 6.5028e-06, 1.00015, 13.2101,
> 17.1875] 21
> 22 -1.87904 [1.00011, 5.11913e-05, 6.90862e-06, 1.00016, 13.1586,
> 17.6349] 22
> 23 -1.85838 [1.00011, 4.99116e-05, 6.32473e-06, 1.00016, 13.2756,
> 17.2941] 23
> 24 -1.87234 [1.00011, 5.06405e-05, 6.38674e-06, 1.00016, 13.142,
> 16.9595] 24
> 25 -1.87851 [1.00011, 5.60623e-05, 9.50702e-06, 1.00016, 12.9261,
> 17.1505] 25
> 26 -1.87278 [1.00011, 5.7016e-05, 9.81952e-06, 1.00017, 13.1566,
> 17.1573] 26
> 27 -1.87788 [1.00011, 5.79397e-05, 1.03623e-05, 1.00016, 13.1869,
> 17.3859] 27
> 28 -1.87119 [1.00012, 5.78678e-05, 1.02708e-05, 1.00017, 13.1865,
> 17.2014] 28
> 29 -1.87867 [1.00012, 5.93291e-05, 1.09259e-05, 1.00017, 13.2404,
> 17.3388] 29
> 30 -1.87139 [1.00012, 6.02041e-05, 1.1396e-05, 1.00017, 13.1466,
> 17.2671] 30
> 31 -1.87646 [1.00012, 6.09047e-05, 1.16142e-05, 1.00017, 13.228,
> 17.2193] 31
> 32 -1.87726 [1.00012, 6.10943e-05, 1.15513e-05, 1.00017, 13.185,
> 17.1572] 32
> 33 -1.87794 [1.00012, 6.12425e-05, 1.16305e-05, 1.00017, 13.1883,
> 17.2175] 33
> 34 -1.87838 [1.00012, 6.23812e-05, 1.21072e-05, 1.00017, 13.2234,
> 17.1842] 34
> 35 -1.87868 [1.00012, 6.25183e-05, 1.21145e-05, 1.00017, 13.2095,
> 17.2203] 35
> 36 -1.87824 [1.00012, 6.2734e-05, 1.21642e-05, 1.00017, 13.1988,
> 17.1913] 36
> 37 -1.87921 [1.00012, 6.28541e-05, 1.22088e-05, 1.00017, 13.2007,
> 17.216] 37
> 38 -1.8782 [1.00012, 6.32819e-05, 1.23541e-05, 1.00017, 13.1978,
> 17.1964] 38
> 39 -1.87855 [1.00012, 6.33837e-05, 1.23923e-05, 1.00017, 13.1993,
> 17.2121] 39
> 40 -1.87925 [1.00013, 6.58186e-05, 1.32416e-05, 1.00018, 13.2, 17.1995]
> 40
> 41 -1.87823 [1.00013, 6.5898e-05, 1.32687e-05, 1.00018, 13.2005,
> 17.2096] 41
> 42 -1.87794 [1.00013, 6.64167e-05, 1.34403e-05, 1.00018, 13.1994,
> 17.2197] 42
> -------------------------------------
> MultiResolution Level : 1
>
> 0 -1.62165 [1.00015, 6.03801e-05, 1.17789e-05, 1.00019, 9.95152,
> 14.8849] 43
> 1 -1.1875 [1.00016, 6.18228e-05, 1.49073e-05, 1.0002, 12.9899, 15.8891]
> 44
> 2 -1.51398 [1.00017, 6.69607e-05, 1.93186e-05, 1.0002, 12.9537,
> 19.0889] 45
> 3 -1.42245 [1.00018, 7.01343e-05, 2.17115e-05, 1.00021, 12.7735,
> 16.5353] 46
> 4 -1.57939 [1.00018, 7.17956e-05, 2.31559e-05, 1.00021, 13.3021,
> 18.5139] 47
> 5 -1.47109 [1.00019, 7.0268e-05, 2.19e-05, 1.00021, 12.7688, 16.9647]
> 48
> 6 -1.61299 [1.00019, 7.05664e-05, 2.1608e-05, 1.00021, 13.97, 17.4892]
> 49
> 7 -1.4959 [1.00019, 7.03122e-05, 2.16426e-05, 1.00021, 12.9971,
> 17.0981] 50
> 8 -1.63338 [1.0002, 6.86924e-05, 2.1011e-05, 1.00022, 13.3451, 16.109]
> 51
> 9 -1.5237 [1.0002, 6.81505e-05, 2.07289e-05, 1.00022, 13.1068, 16.9133]
> 52
> 10 -1.62846 [1.0002, 6.82681e-05, 2.02992e-05, 1.00022, 12.4269,
> 17.4046] 53
> 11 -1.55341 [1.0002, 6.76953e-05, 1.98522e-05, 1.00022, 13.0085,
> 17.0697] 54
> 12 -1.63467 [1.00021, 6.639e-05, 1.92588e-05, 1.00023, 12.8479,
> 16.4181] 55
> 13 -1.57099 [1.00021, 6.704e-05, 1.98256e-05, 1.00023, 12.9785,
> 16.9389] 56
> 14 -1.6313 [1.00021, 6.89057e-05, 2.10809e-05, 1.00023, 13.1016,
> 17.4615] 57
> 15 -1.60176 [1.00021, 6.8814e-05, 2.06596e-05, 1.00023, 12.8971,
> 17.0838] 58
> 16 -1.62815 [1.00021, 6.93416e-05, 2.10816e-05, 1.00023, 13.2334,
> 17.013] 59
> 17 -1.62391 [1.00021, 6.90842e-05, 2.10131e-05, 1.00024, 12.9709,
> 16.9313] 60
> 18 -1.63013 [1.00022, 6.98766e-05, 2.14936e-05, 1.00024, 13.0301,
> 17.1431] 61
> 19 -1.63213 [1.00022, 6.96417e-05, 2.14305e-05, 1.00024, 13.0242,
> 16.9672] 62
> 20 -1.63369 [1.00022, 6.97123e-05, 2.14339e-05, 1.00024, 12.9752,
> 17.0992] 63
> 21 -1.63229 [1.00022, 6.98076e-05, 2.1516e-05, 1.00024, 13.0534,
> 17.0181] 64
> 22 -1.63476 [1.00022, 6.9642e-05, 2.1329e-05, 1.00024, 12.9642,
> 17.0051] 65
> 23 -1.63365 [1.00022, 6.99387e-05, 2.15251e-05, 1.00024, 13.0228,
> 16.9632] 66
> 24 -1.63359 [1.00022, 6.9969e-05, 2.15261e-05, 1.00024, 13.004,
> 17.0176] 67
> 25 -1.63559 [1.00022, 6.98748e-05, 2.14739e-05, 1.00024, 12.9975,
> 16.972] 68
> 26 -1.63425 [1.00022, 6.99632e-05, 2.15374e-05, 1.00024, 13.0102,
> 17.0066] 69
> 27 -1.63572 [1.00022, 6.99392e-05, 2.15124e-05, 1.00024, 12.9926,
> 16.983] 70
> 28 -1.63557 [1.00022, 6.99708e-05, 2.16214e-05, 1.00024, 13.0006,
> 16.9546] 71
> 29 -1.6335 [1.00022, 7.00353e-05, 2.16619e-05, 1.00024, 12.9988,
> 16.9781] 72
> 30 -1.6355 [1.00022, 7.0194e-05, 2.17912e-05, 1.00024, 13.0155,
> 16.9948] 73
> 31 -1.63563 [1.00022, 7.01802e-05, 2.17694e-05, 1.00024, 13.0007,
> 16.9831] 74
> 32 -1.63564 [1.00022, 7.01415e-05, 2.17686e-05, 1.00024, 12.9848,
> 16.9728] 75
> 33 -1.6343 [1.00022, 7.02017e-05, 2.18139e-05, 1.00024, 12.9946,
> 16.9844] 76
> 34 -1.63558 [1.00022, 7.01902e-05, 2.18266e-05, 1.00024, 12.9949,
> 16.9723] 77
> 35 -1.6344 [1.00022, 7.02139e-05, 2.18457e-05, 1.00024, 12.9998,
> 16.9807] 78
> 36 -1.63558 [1.00022, 7.02012e-05, 2.18512e-05, 1.00024, 12.9932,
> 16.9766] 79
> 37 -1.63544 [1.00022, 7.02338e-05, 2.18765e-05, 1.00024, 12.9986,
> 16.9796] 80
> 38 -1.63556 [1.00022, 7.02559e-05, 2.1918e-05, 1.00024, 12.9996,
> 16.9748] 81
> 39 -1.63515 [1.00022, 7.02707e-05, 2.19276e-05, 1.00024, 13, 16.9787]
> 82
> 40 -1.63555 [1.00022, 7.02872e-05, 2.19443e-05, 1.00024, 12.9967,
> 16.981] 83
> 41 -1.63561 [1.00022, 7.02763e-05, 2.19467e-05, 1.00024, 12.9962,
> 16.9779] 84
> 42 -1.63553 [1.00022, 7.0289e-05, 2.19585e-05, 1.00024, 12.9988,
> 16.9782] 85
> 43 -1.63555 [1.00023, 7.03123e-05, 2.19736e-05, 1.00024, 13.001,
> 16.9793] 86
> 44 -1.63556 [1.00023, 7.03133e-05, 2.19762e-05, 1.00024, 12.999,
> 16.9793] 87
> 45 -1.63558 [1.00023, 7.03392e-05, 2.2003e-05, 1.00024, 12.997,
> 16.9797] 88
> 46 -1.63557 [1.00023, 7.0332e-05, 2.20053e-05, 1.00024, 12.997,
> 16.9781] 89
> 47 -1.63555 [1.00023, 7.03397e-05, 2.20109e-05, 1.00024, 12.9983,
> 16.9781] 90
> 48 -1.6355 [1.00023, 7.03503e-05, 2.20177e-05, 1.00024, 12.9995,
> 16.9785] 91
> 49 -1.63558 [1.00023, 7.03612e-05, 2.20263e-05, 1.00025, 12.9992,
> 16.9798] 92
> 50 -1.63564 [1.00023, 7.03605e-05, 2.2028e-05, 1.00025, 12.9983,
> 16.9792] 93
> -------------------------------------
> MultiResolution Level : 2
>
> 0 -1.28555 [1.00023, 7.17596e-05, 2.34728e-05, 1.00025, 12.3864,
> 17.7702] 94
> 1 -1.20689 [1.00023, 7.17569e-05, 2.33515e-05, 1.00025, 13.0162,
> 17.2769] 95
> 2 -1.2882 [1.00023, 7.20532e-05, 2.33524e-05, 1.00026, 12.8846,
> 16.4878] 96
> 3 -1.2738 [1.00023, 7.32712e-05, 2.40743e-05, 1.00026, 13.2028,
> 17.0431] 97
> 4 -1.29413 [1.00023, 7.2987e-05, 2.41593e-05, 1.00026, 12.7043,
> 16.9264] 98
> 5 -1.29154 [1.00023, 7.28283e-05, 2.39671e-05, 1.00026, 13.0967,
> 17.0438] 99
> 6 -1.29226 [1.00023, 7.32004e-05, 2.42714e-05, 1.00026, 12.8308,
> 16.8523] 100
> 7 -1.30052 [1.00023, 7.35853e-05, 2.44347e-05, 1.00026, 13.0505,
> 16.9953] 101
> 8 -1.28733 [1.00023, 7.39501e-05, 2.50317e-05, 1.00026, 12.8635,
> 17.0903] 102
> 9 -1.29795 [1.00023, 7.41911e-05, 2.51215e-05, 1.00026, 12.9954,
> 16.9866] 103
> 10 -1.28539 [1.00023, 7.44061e-05, 2.54409e-05, 1.00026, 12.88,
> 17.0551] 104
> 11 -1.29558 [1.00023, 7.46609e-05, 2.55977e-05, 1.00026, 12.9805,
> 17.0172] 105
> 12 -1.28655 [1.00023, 7.59312e-05, 2.66034e-05, 1.00026, 13.0043,
> 16.9125] 106
> 13 -1.29005 [1.00023, 7.58524e-05, 2.66422e-05, 1.00026, 12.9799,
> 16.9949] 107
> 14 -1.28649 [1.00023, 7.65503e-05, 2.72752e-05, 1.00026, 12.9877,
> 17.0805] 108
> 15 -1.2893 [1.00023, 7.67555e-05, 2.74021e-05, 1.00026, 13.0005,
> 17.013] 109
> 16 -1.28504 [1.00023, 7.71766e-05, 2.77733e-05, 1.00026, 12.9456,
> 17.0158] 110
> 17 -1.28871 [1.00023, 7.72349e-05, 2.78222e-05, 1.00026, 12.9872,
> 17.0017] 111
> 18 -1.28581 [1.00023, 7.74398e-05, 2.80487e-05, 1.00026, 12.9526,
> 17.008] 112
> 19 -1.28831 [1.00023, 7.75173e-05, 2.81092e-05, 1.00026, 12.9804,
> 17.0038] 113
> 20 -1.28645 [1.00023, 7.8997e-05, 2.94688e-05, 1.00026, 12.972,
> 17.0246] 114
> 21 -1.2875 [1.00023, 7.91652e-05, 2.9606e-05, 1.00027, 12.9829,
> 17.0103] 115
> 22 -1.28633 [1.00023, 7.94498e-05, 2.98614e-05, 1.00027, 12.9738,
> 16.9947] 116
> 23 -1.28705 [1.00023, 7.95571e-05, 2.99512e-05, 1.00027, 12.9833,
> 17.0056] 117
> 24 -1.28626 [1.00023, 7.97569e-05, 3.01446e-05, 1.00027, 12.972,
> 17.0035] 118
> 25 -1.28721 [1.00023, 7.9854e-05, 3.0225e-05, 1.00027, 12.981, 17.0054]
> 119
> 26 -1.28646 [1.00023, 8.03125e-05, 3.06464e-05, 1.00027, 12.9739,
> 17.0035] 120
> 27 -1.28705 [1.00023, 8.03901e-05, 3.07116e-05, 1.00027, 12.9796,
> 17.0051] 121
> 28 -1.28658 [1.00022, 8.14113e-05, 3.1626e-05, 1.00027, 12.9851,
> 17.0072] 122
> 29 -1.28612 [1.00022, 8.14614e-05, 3.16754e-05, 1.00027, 12.9806,
> 17.0057] 123
> 30 -1.28651 [1.00022, 8.18571e-05, 3.20354e-05, 1.00027, 12.9765,
> 17.0034] 124
> 31 -1.28683 [1.00022, 8.19314e-05, 3.20996e-05, 1.00027, 12.9799,
> 17.005] 125
> 32 -1.28657 [1.00022, 8.33454e-05, 3.33747e-05, 1.00027, 12.9794,
> 17.0088] 126
> 33 -1.28666 [1.00022, 8.34364e-05, 3.34534e-05, 1.00027, 12.9801,
> 17.0059] 127
> 34 -1.28657 [1.00022, 8.38365e-05, 3.38138e-05, 1.00027, 12.9782,
> 17.0035] 128
> 35 -1.28671 [1.00022, 8.39123e-05, 3.3881e-05, 1.00027, 12.98, 17.0052]
> 129
> 36 -1.28659 [1.00022, 8.44641e-05, 3.43794e-05, 1.00027, 12.9783,
> 17.0062] 130
> 37 -1.28675 [1.00022, 8.45361e-05, 3.4442e-05, 1.00027, 12.9798,
> 17.0058] 131
> 38 -1.28663 [1.00022, 8.49391e-05, 3.48032e-05, 1.00027, 12.9796,
> 17.0042] 132
> 39 -1.28664 [1.00022, 8.5028e-05, 3.48835e-05, 1.00027, 12.9797,
> 17.0055] 133
> 40 -1.28663 [1.00022, 8.778e-05, 3.73549e-05, 1.00027, 12.9789, 17.005]
> 134
> 41 -1.28671 [1.00022, 8.784e-05, 3.74083e-05, 1.00027, 12.9795,
> 17.0054] 135
> 42 -1.28666 [1.00022, 8.82254e-05, 3.77538e-05, 1.00027, 12.9798,
> 17.0062] 136
> 43 -1.28667 [1.00022, 8.8309e-05, 3.78288e-05, 1.00027, 12.9795,
> 17.0056] 137
> 44 -1.28667 [1.00022, 8.87527e-05, 3.82256e-05, 1.00027, 12.98,
> 17.0056] 138
> 45 -1.28664 [1.00022, 8.88064e-05, 3.82742e-05, 1.00027, 12.9796,
> 17.0056] 139
> 46 -1.28666 [1.00022, 8.91826e-05, 3.86117e-05, 1.00027, 12.9792,
> 17.0058] 140
> 47 -1.28671 [1.00022, 8.92549e-05, 3.86761e-05, 1.00027, 12.9795,
> 17.0057] 141
> 48 -1.28669 [1.00022, 8.95278e-05, 3.89205e-05, 1.00027, 12.9793,
> 17.0056] 142
> 49 -1.2867 [1.00022, 8.96052e-05, 3.89897e-05, 1.00027, 12.9795,
> 17.0057] 143
> 50 -1.28669 [1.00022, 8.9917e-05, 3.9269e-05, 1.00027, 12.9793,
> 17.0057] 144
> 51 -1.28671 [1.00022, 8.99883e-05, 3.93327e-05, 1.00027, 12.9795,
> 17.0057] 145
> 52 -1.2867 [1.00022, 9.02924e-05, 3.96051e-05, 1.00027, 12.9794,
> 17.0057] 146
> Optimizer stop condition: RegularStepGradientDescentOptimizer: Step too
> small after 53 iterations. Current step (8.50706e-05) is less than minimum
> step (0.0001).
> Optimizer Stopping Condition = 2
> Result =
> Translation X = 12.9794
> Translation Y = 17.0057
> Iterations = 54
> Metric value = -1.2867
>
>
> Any idea, why my run is taking too many iteration to complete ???
> Cheers
> Nain
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>
More information about the Insight-users
mailing list