[Insight-users] Mail on Registration.
valli gummadi
valli gummadi" <valli_gummadi at rediffmail.com
9 Apr 2004 05:04:12 -0000
This is a multipart mime message
--Next_1081487052---0-202.54.124.151-8287
Content-type: text/html;
charset=iso-8859-1
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline
<P>=0A <BR>=0AHi All,<BR>=0AI am posting the same message again and a=
gain because I am really struck to Proceed further in my Registration appli=
cation. And I am in need of this registration process.And,I could achieve t=
he registration process when both datasets are in same view.only problem is=
when both are in different view.Any help is favourable to me.Thank you in =
advance.<BR>=0A<BR>=0A I have been working in Registration of MRI and SPECT=
Images of same <BR>=0Apatient Brain.<BR>=0AMR dimensions are 256 X 256 X 1=
20 and 0.859382 \ 0.859375 \ 1.60000.<BR>=0AMy moving image is SPECT =
image.Actual dimensions are 128 X 128 X 47 with voxel dimensions 1.802125\1=
.802125\3.60425.<BR>=0ASPECT image undergoes a ForceCubic operation and Its=
dimensions changes to 128X 128 X 95 and voxel dimesnions to 1.802125 \ 1.8=
02125 \ 1.802125.<BR>=0AAnd it is in Transverse View.<BR>=0AAs MR image is =
in Coronal view i changed the Spect to coronal before <BR>=0Agiving it to t=
he Registration. So,After changing its view to coronal its dimensions becom=
e<BR>=0A128 X 95 X 128 and 1.802125 \ 1.802125 \ 1.802125.SPECT=
with these dimensions is given to the registration process.<BR>=0AI used t=
he following itk classes:<BR>=0AMattesMutualInformationImageToImageMetric,R=
ecursiveMultiResolutionPyramidImageFilter,<BR>=0AQuaternionRigidTransform,R=
egularStepGradientDescentOptimizer,RecursiveMultiResolutionPyramidImageFilt=
er.<BR>=0A<BR>=0AAnd the input parameters are :<BR>=0ATranslate scale value=
is 1/10000, iterations are 150, scales =3D3D 1, max step length =3D 4.0, m=
in step length =3D 0.005, bins =3D 35, samples =3D<BR>=0A1000, levels=
=3D 3.<BR>=0A<BR>=0AI added observer to the process and i am getting the f=
ollowing values :<BR>=0AFirst Level:<BR>=0A 0 -0.320374 =
[0.0100438, -0.0350144, 0.0105636, 0.986362, =3D<BR>=0A-0.553508, -0=
.00411211, 0.831858]<BR>=0A 1 -0.321063 [0.030953=
, -0.0382486, -0.00937501, 0.999558, =3D<BR>=0A-0.00125347, -0.440459, 1.54=
15]<BR>=0A 2 -0.371695 [0.0659527, -0.0315038, -0=
.0302047, 0.974523, =3D<BR>=0A0.0133357, -1.43741, 1.60113]<BR>=0A 3 =
-0.389297 [0.0748733, -0.0149963, -0.0518501, 0.992978,=
=3D<BR>=0A0.962366, -1.64413, 1.8366]<BR>=0A 4 -0.473205 =
; [0.0765801, -0.0308507, -0.0461271, 1.0002, 0.831316, =3D<BR>=0A-1=
.52403, 2.30357]<BR>=0A 5 -0.489634 [0.0867224, -=
0.0318137, -0.0522786, 1.01073, 1.15721, =3D<BR>=0A-1.6065, 2.67335]<BR>=0A=
-------------------------------<BR>=0A20 -0.526699  =
; [0.0869593, -0.0339316, -0.0559707, 1.01556, 1.18742, =3D<BR>=0A-1.66461,=
2.72034]<BR>=0A 21 -0.526875 [0.086819, -0.03405=
73, -0.056085, 1.01558, 1.19139, =3D<BR>=0A-1.65994, 2.72518]<BR>=0A 22&nbs=
p; -0.527021 [0.0867698, -0.0341474, -0.0562826, 1.0156=
, 1.19699, =3D<BR>=0A-1.65843, 2.73041]<BR>=0A 23 -0.527126&nb=
sp; [0.0866777, -0.0342297, -0.0564333, 1.01563, =3D<BR>=0A1.20222, =
-1.65538, 2.73534]<BR>=0A-----------------------------------------------<BR=
>=0A47 -0.528611 [0.0852474, -0.0356566, -0.05851=
8, 1.015, 1.25176, =3D<BR>=0A-1.57546, 2.82544]<BR>=0A 48 -0.5=
28646 [0.0852217, -0.0356863, -0.0585624, 1.015, 1.25251, =3D<=
BR>=0A-1.57258, 2.82797]<BR>=0A 49 -0.528681 [0.0=
852033, -0.0357116, -0.0586277, 1.01498, =3D<BR>=0A1.25373, -1.57004, 2.830=
68]<BR>=0A---------------------------------------------<BR>=0A 144&nb=
sp; -0.525427 [0.0845718, -0.0363607, -0.0603638, 1.014=
22, =3D<BR>=0A1.22735, -1.23997, 2.96168]<BR>=0A 145 -0.525436=
[0.0845504, -0.0363369, -0.06038, 1.01423, 1.22755, =3D<BR>=
=0A-1.2361, 2.96218]<BR>=0A 146 -0.52494 [0.08453=
61, -0.0363315, -0.0603705, 1.01426, =3D<BR>=0A1.22752, -1.23235, 2.96328]<=
BR>=0A 147 -0.52495 [0.0847016, -0.036275, -0.060=
5773, 1.01413, 1.22964, =3D<BR>=0A-1.23218, 2.96654]<BR>=0A 148  =
; -0.525374 [0.0846794, -0.0363122, -0.0605381, 1.01414, =3D<B=
R>=0A1.22845, -1.23069, 2.96692]<BR>=0A 149 -0.525412 &n=
bsp; [0.0846573, -0.036348, -0.0605014, 1.01415, =3D<BR>=0A1.22734, -1.2291=
3, 2.96732]<BR>=0ASecond Level:<BR>=0A0 -0.418512 =
[0.0887451, -0.0189369, -0.077675, 1.00696, 2.20275, =3D<BR>=0A-1.37414, 3=
.13128]<BR>=0A 1 -0.418945 [0.0836208, -0.0236872=
, -0.0630348, 1.01192, 2.22691, =3D<BR>=0A-0.917966, 3.33384]<BR>=0A 2 =
; -0.413912 [0.0652313, -0.0161832, -0.0486405, 1.00361=
, 2.10126, =3D<BR>=0A-0.505873, 3.08141]<BR>=0A 3 -0.371655&nb=
sp; [0.0629191, -0.00723327, -0.0532294, 1.01451, =3D<BR>=0A2.54857,=
-0.291537, 3.14261]<BR>=0A 4 -0.361762 [0.071574=
1, -0.0131997, -0.059176, 1.02901, 2.72654, =3D<BR>=0A-0.369328, 3.60296]<B=
R>=0A 5 -0.367962 [0.0714853, -0.0257059, -0.0513=
781, 1.02656, 2.41281, =3D<BR>=0A-0.219652, 3.96205]<BR>=0A 6 =
-0.388251 [0.0817307, -0.0310034, -0.0601911, 1.01481, 2.18368=
, =3D<BR>=0A-0.662284, 3.99708]<BR>=0A 7 -0.41388 =
[0.0881204, -0.03356, -0.0638268, 0.991765, 1.87106, =3D<BR>=0A-1.04525, 3=
.92623]<BR>=0A----------------------------<BR>=0A53 -0.444609&=
nbsp; [0.0904033, -0.0317655, -0.0649183, 0.986647, =3D<BR>=0A2.2865=
2, -1.12252, 4.21855]<BR>=0A 54 -0.444668 [0.0903=
786, -0.0317895, -0.0648504, 0.986712, =3D<BR>=0A2.29139, -1.12121, 4.22452=
]<BR>=0A 55 -0.444666 [0.0904111, -0.0317457, -0.=
064888, 0.986666, =3D<BR>=0A2.29744, -1.12326, 4.22902]<BR>=0A 56 &nb=
sp; -0.44467 [0.0903941, -0.0317951, -0.0648093, 0.986712, 2.3=
017, =3D<BR>=0A-1.12212, 4.23546]<BR>=0A 57 -0.44468 &nb=
sp; [0.0904166, -0.0317319, -0.064856, 0.986649, 2.30811, =3D<BR>=0A-1.1241=
1, 4.23945]<BR>=0A 58 -0.444699 [0.0903994, -0.03=
17871, -0.0647737, 0.986687, =3D<BR>=0A2.31225, -1.12307, 4.24599]<BR>=0A--=
--------------------------------<BR>=0A 146 -0.445491 &n=
bsp; [0.0900287, -0.0319401, -0.0639129, 0.985568, =3D<BR>=0A2.73939, -1.19=
044, 4.6637]<BR>=0A 147 -0.445511 [0.0900287, -0.=
031949, -0.0639108, 0.985541, =3D<BR>=0A2.74461, -1.19175, 4.66937]<BR>=0A =
148 -0.445525 [0.0900161, -0.0319461, -0.0639007,=
0.985548, =3D<BR>=0A2.75007, -1.19225, 4.67494]<BR>=0A 149 -0=
.445532 [0.090016, -0.0319676, -0.0638845, 0.985523, =3D<BR>=
=0A2.75507, -1.19321, 4.68086]<BR>=0AThird Level:<BR>=0A0 -0.3=
02895 [0.112791, -0.0365895, -0.0844165, 1.00896, 3.33379, =3D=
<BR>=0A-1.60939, 5.38112]<BR>=0A 1 -0.322343 [0.0=
992007, -0.0304807, -0.0782828, 1.00517, 3.25472, =3D<BR>=0A-1.32882, 4.975=
22]<BR>=0A 2 -0.313556 [0.101683, -0.0331957, -0.=
0747196, 1.01122, 3.1463, =3D<BR>=0A-1.1917, 5.15378]<BR>=0A 3 =
-0.328914 [0.0990012, -0.0300732, -0.0743876, 1.0092, 3.13413=
, =3D<BR>=0A-1.19278, 5.02946]<BR>=0A 4 -0.330295 =
[0.0997371, -0.0316464, -0.0740325, 1.00807, 3.0786, =3D<BR>=0A-1.21719, 5=
.0444]<BR>=0A 5 -0.332905 [0.100176, -0.0317067, =
-0.0743968, 1.00936, 3.0995, =3D<BR>=0A-1.2133, 5.06726]<BR>=0A------------=
---------------<BR>=0A35 -0.338694 [0.100619, -0.=
031991, -0.0749333, 1.01022, 3.03826, =3D<BR>=0A-1.2607, 5.02906]<BR>=0A 36=
-0.33863 [0.100617, -0.0320048, -0.0749167, 1.01=
019, 3.03661, =3D<BR>=0A-1.26134, 5.02824]<BR>=0A 37 -0.338678=
[0.100626, -0.0320129, -0.0749159, 1.01021, 3.03511, =3D<BR>=
=0A-1.26227, 5.02739]<BR>=0A 38 -0.338677 [0.1006=
2, -0.0319927, -0.0749215, 1.01021, 3.03419, =3D<BR>=0A-1.26311, 5.02589]<B=
R>=0A-------------------------<BR>=0A59 -0.338709 =
[0.100639, -0.0320017, -0.0748891, 1.01031, 3.00697, =3D<BR>=0A-1.27972, 5=
.00324]<BR>=0A 60 -0.338664 [0.100631, -0.0319735=
, -0.0749054, 1.01029, 3.00644, =3D<BR>=0A-1.28063, 5.0016]<BR>=0A 61 =
-0.338709 [0.100639, -0.0320041, -0.0748886, 1.01031, =
3.00458, =3D<BR>=0A-1.28112, 5.00129]<BR>=0A 62 -0.338667 =
; [0.10063, -0.0319716, -0.074908, 1.0103, 3.00421, =3D<BR>=0A-1.282=
05, 4.9996]<BR>=0A 63 -0.33871 [0.100645, -0.0320=
181, -0.074882, 1.01033, 3.00229, =3D<BR>=0A-1.28225, 4.99989]<BR>=0A<BR>=
=0AFirst and Second Levels took all the 150 iterations and at third levels =
process came out at 65th iteration. I tried to change the tranlate scale va=
lue to less than 1000. Becuase for mr volumes 100 is sufficient.But, I am g=
ettign the itk errorr message "All the sampled points are moved outof =
the image regios.<BR>=0APlease anybody explain me the better way to achieve=
good registration <BR>=0Aoutput.<BR>=0A<BR>=0ARegards,<BR>=0ASrivalli.<BR>=
=0A=0A</P>=0A<br><br>=0A<A target=3D"_blank" HREF=3D"http://clients.rediff.=
com/signature/track_sig.asp"><IMG SRC=3D"http://ads.rediff.com/RealMedia/ad=
s/adstream_nx.cgi/www.rediffmail.com/inbox.htm at Bottom" BORDER=3D0 VSPACE=3D=
0 HSPACE=3D0></a>=0A
--Next_1081487052---0-202.54.124.151-8287
Content-type: text/plain;
charset=iso-8859-1
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline
=A0=0AHi All,=0AI am posting the same message again and again because I am=
really struck to Proceed further in my Registration application. And I am =
in need of this registration process.And,I could achieve the registration p=
rocess when both datasets are in same view.only problem is when both are in=
different view.Any help is favourable to me.Thank you in advance.=0A=0A I =
have been working in Registration of MRI and SPECT Images of same =0Apatien=
t Brain.=0AMR dimensions are 256 X 256 X 120 and 0.859382 \ 0.859375 \ 1.6=
0000.=0AMy moving image is SPECT image.Actual dimensions are 128 X 128 X 47=
with voxel dimensions 1.802125\1.802125\3.60425.=0ASPECT image undergoes a=
ForceCubic operation and Its dimensions changes to 128X 128 X 95 and voxel=
dimesnions to 1.802125 \ 1.802125 \ 1.802125.=0AAnd it is in Transverse Vi=
ew.=0AAs MR image is in Coronal view i changed the Spect to coronal before =
=0Agiving it to the Registration. So,After changing its view to coronal its=
dimensions become=0A128 X 95 X 128 and 1.802125 \ 1.802125 \ 1.802125.SP=
ECT with these dimensions is given to the registration process.=0AI used th=
e following itk classes:=0AMattesMutualInformationImageToImageMetric,Recurs=
iveMultiResolutionPyramidImageFilter,=0AQuaternionRigidTransform,RegularSte=
pGradientDescentOptimizer,RecursiveMultiResolutionPyramidImageFilter.=0A=0A=
And the input parameters are :=0ATranslate scale value is 1/10000, iteratio=
ns are 150, scales =3D3D 1, max step length =3D 4.0, min step length =3D 0.=
005, bins =3D 35, samples =3D=0A1000, levels =3D 3.=0A=0AI added observer =
to the process and i am getting the following values :=0AFirst Level:=0A 0 =
-0.320374 [0.0100438, -0.0350144, 0.0105636, 0.986362, =3D=0A-0.55350=
8, -0.00411211, 0.831858]=0A 1 -0.321063 [0.030953, -0.0382486, -0.00=
937501, 0.999558, =3D=0A-0.00125347, -0.440459, 1.5415]=0A 2 -0.371695 =
[0.0659527, -0.0315038, -0.0302047, 0.974523, =3D=0A0.0133357, -1.43741, =
1.60113]=0A 3 -0.389297 [0.0748733, -0.0149963, -0.0518501, 0.992978,=
=3D=0A0.962366, -1.64413, 1.8366]=0A 4 -0.473205 [0.0765801, -0.0308=
507, -0.0461271, 1.0002, 0.831316, =3D=0A-1.52403, 2.30357]=0A 5 -0.4896=
34 [0.0867224, -0.0318137, -0.0522786, 1.01073, 1.15721, =3D=0A-1.6065, =
2.67335]=0A-------------------------------=0A20 -0.526699 [0.0869593,=
-0.0339316, -0.0559707, 1.01556, 1.18742, =3D=0A-1.66461, 2.72034]=0A 21 =
-0.526875 [0.086819, -0.0340573, -0.056085, 1.01558, 1.19139, =3D=0A-1=
.65994, 2.72518]=0A 22 -0.527021 [0.0867698, -0.0341474, -0.0562826, =
1.0156, 1.19699, =3D=0A-1.65843, 2.73041]=0A 23 -0.527126 [0.0866777,=
-0.0342297, -0.0564333, 1.01563, =3D=0A1.20222, -1.65538, 2.73534]=0A-----=
------------------------------------------=0A47 -0.528611 [0.0852474,=
-0.0356566, -0.058518, 1.015, 1.25176, =3D=0A-1.57546, 2.82544]=0A 48 -=
0.528646 [0.0852217, -0.0356863, -0.0585624, 1.015, 1.25251, =3D=0A-1.57=
258, 2.82797]=0A 49 -0.528681 [0.0852033, -0.0357116, -0.0586277, 1.0=
1498, =3D=0A1.25373, -1.57004, 2.83068]=0A---------------------------------=
------------=0A 144 -0.525427 [0.0845718, -0.0363607, -0.0603638, 1.=
01422, =3D=0A1.22735, -1.23997, 2.96168]=0A 145 -0.525436 [0.0845504,=
-0.0363369, -0.06038, 1.01423, 1.22755, =3D=0A-1.2361, 2.96218]=0A 146 =
-0.52494 [0.0845361, -0.0363315, -0.0603705, 1.01426, =3D=0A1.22752, -1.=
23235, 2.96328]=0A 147 -0.52495 [0.0847016, -0.036275, -0.0605773, 1.=
01413, 1.22964, =3D=0A-1.23218, 2.96654]=0A 148 -0.525374 [0.0846794,=
-0.0363122, -0.0605381, 1.01414, =3D=0A1.22845, -1.23069, 2.96692]=0A 149 =
-0.525412 [0.0846573, -0.036348, -0.0605014, 1.01415, =3D=0A1.22734, =
-1.22913, 2.96732]=0ASecond Level:=0A0 -0.418512 [0.0887451, -0.01893=
69, -0.077675, 1.00696, 2.20275, =3D=0A-1.37414, 3.13128]=0A 1 -0.418945=
[0.0836208, -0.0236872, -0.0630348, 1.01192, 2.22691, =3D=0A-0.917966, =
3.33384]=0A 2 -0.413912 [0.0652313, -0.0161832, -0.0486405, 1.00361, =
2.10126, =3D=0A-0.505873, 3.08141]=0A 3 -0.371655 [0.0629191, -0.0072=
3327, -0.0532294, 1.01451, =3D=0A2.54857, -0.291537, 3.14261]=0A 4 -0.36=
1762 [0.0715741, -0.0131997, -0.059176, 1.02901, 2.72654, =3D=0A-0.36932=
8, 3.60296]=0A 5 -0.367962 [0.0714853, -0.0257059, -0.0513781, 1.0265=
6, 2.41281, =3D=0A-0.219652, 3.96205]=0A 6 -0.388251 [0.0817307, -0.0=
310034, -0.0601911, 1.01481, 2.18368, =3D=0A-0.662284, 3.99708]=0A 7 -0.=
41388 [0.0881204, -0.03356, -0.0638268, 0.991765, 1.87106, =3D=0A-1.0452=
5, 3.92623]=0A----------------------------=0A53 -0.444609 [0.0904033,=
-0.0317655, -0.0649183, 0.986647, =3D=0A2.28652, -1.12252, 4.21855]=0A 54 =
-0.444668 [0.0903786, -0.0317895, -0.0648504, 0.986712, =3D=0A2.29139=
, -1.12121, 4.22452]=0A 55 -0.444666 [0.0904111, -0.0317457, -0.06488=
8, 0.986666, =3D=0A2.29744, -1.12326, 4.22902]=0A 56 -0.44467 [0.0903=
941, -0.0317951, -0.0648093, 0.986712, 2.3017, =3D=0A-1.12212, 4.23546]=0A =
57 -0.44468 [0.0904166, -0.0317319, -0.064856, 0.986649, 2.30811, =3D=
=0A-1.12411, 4.23945]=0A 58 -0.444699 [0.0903994, -0.0317871, -0.0647=
737, 0.986687, =3D=0A2.31225, -1.12307, 4.24599]=0A------------------------=
----------=0A 146 -0.445491 [0.0900287, -0.0319401, -0.0639129, 0.985=
568, =3D=0A2.73939, -1.19044, 4.6637]=0A 147 -0.445511 [0.0900287, -0=
.031949, -0.0639108, 0.985541, =3D=0A2.74461, -1.19175, 4.66937]=0A 148 =
-0.445525 [0.0900161, -0.0319461, -0.0639007, 0.985548, =3D=0A2.75007, -=
1.19225, 4.67494]=0A 149 -0.445532 [0.090016, -0.0319676, -0.0638845,=
0.985523, =3D=0A2.75507, -1.19321, 4.68086]=0AThird Level:=0A0 -0.30289=
5 [0.112791, -0.0365895, -0.0844165, 1.00896, 3.33379, =3D=0A-1.60939, 5=
.38112]=0A 1 -0.322343 [0.0992007, -0.0304807, -0.0782828, 1.00517, 3=
.25472, =3D=0A-1.32882, 4.97522]=0A 2 -0.313556 [0.101683, -0.0331957=
, -0.0747196, 1.01122, 3.1463, =3D=0A-1.1917, 5.15378]=0A 3 -0.328914 =
[0.0990012, -0.0300732, -0.0743876, 1.0092, 3.13413, =3D=0A-1.19278, 5.029=
46]=0A 4 -0.330295 [0.0997371, -0.0316464, -0.0740325, 1.00807, 3.078=
6, =3D=0A-1.21719, 5.0444]=0A 5 -0.332905 [0.100176, -0.0317067, -0.0=
743968, 1.00936, 3.0995, =3D=0A-1.2133, 5.06726]=0A------------------------=
---=0A35 -0.338694 [0.100619, -0.031991, -0.0749333, 1.01022, 3.03826=
, =3D=0A-1.2607, 5.02906]=0A 36 -0.33863 [0.100617, -0.0320048, -0.07=
49167, 1.01019, 3.03661, =3D=0A-1.26134, 5.02824]=0A 37 -0.338678 [0.=
100626, -0.0320129, -0.0749159, 1.01021, 3.03511, =3D=0A-1.26227, 5.02739]=
=0A 38 -0.338677 [0.10062, -0.0319927, -0.0749215, 1.01021, 3.03419, =
=3D=0A-1.26311, 5.02589]=0A-------------------------=0A59 -0.338709 [=
0.100639, -0.0320017, -0.0748891, 1.01031, 3.00697, =3D=0A-1.27972, 5.00324=
]=0A 60 -0.338664 [0.100631, -0.0319735, -0.0749054, 1.01029, 3.00644=
, =3D=0A-1.28063, 5.0016]=0A 61 -0.338709 [0.100639, -0.0320041, -0.0=
748886, 1.01031, 3.00458, =3D=0A-1.28112, 5.00129]=0A 62 -0.338667 [0=
.10063, -0.0319716, -0.074908, 1.0103, 3.00421, =3D=0A-1.28205, 4.9996]=0A =
63 -0.33871 [0.100645, -0.0320181, -0.074882, 1.01033, 3.00229, =3D=
=0A-1.28225, 4.99989]=0A=0AFirst and Second Levels took all the 150 iterati=
ons and at third levels process came out at 65th iteration. I tried to chan=
ge the tranlate scale value to less than 1000. Becuase for mr volumes 100 i=
s sufficient.But, I am gettign the itk errorr message "All the sampled poin=
ts are moved outof the image regios.=0APlease anybody explain me the better=
way to achieve good registration =0Aoutput.=0A=0ARegards,=0ASrivalli.=0A
--Next_1081487052---0-202.54.124.151-8287--