[Insight-users] difference in MultiResImageRegistration2 example output
Jagjeet Nain
jagjeet.nain at gmail.com
Sat Oct 29 18:50:54 EDT 2011
Hi Matt,
Thanks for your response.
This shows test as "passed".
with regards
Nain
On Sat, Oct 29, 2011 at 6:10 PM, Matt McCormick
<matt.mccormick at kitware.com>wrote:
> 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
> >
> >
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20111029/bb1be476/attachment.htm>
More information about the Insight-users
mailing list