No subject
Thu May 27 11:52:52 EDT 2010
1. Vector corresponding to EValue =3D -9 which is EVector =3D =
(0.11, 0.11, 0.78) show the direction of
The tube structure.=20
2. Vector corresponding to EValue =3D -634 and -1389 which are =
EVector =3D (0.99, 0.13, 0) and (0.13, 0.99, 0)
respectively represent the plane perpendicular to the tube direction. (I =
name the plane as Plane1)
=20
Goal : I want to extract the neighbour points and their gray values of =
point P and lie on plane Plane1.=20
=20
I thought about the below pipeline:
Step 1: Read image=20
Step 2: Rotate image using EVector =3D (0.11, 0.11, 0.78) (using =
itkResampleImageFilter, rotation transform, trilinear interpolator)
Step 3: Rotate image using EVector =3D (0.99, 0.13, 0) (using =
itkResampleImageFilter, rotation transform, trilinear interpolator)
Step 4: Rotate image using EVector =3D (0.13, 0.99, 0) (using =
itkResampleImageFilter, rotation transform, trilinear interpolator)
Step 5: In the output image from step 4 use neighbourhood iterator to =
obtain the neighbour points.
=20
What i need from ITK users:
1. I request the itk community to check whether above pipeline is =
correct to achieve my goal.
2. If you have any other pipeline in mind it would be useful for =
me.
=20
Thank you in advance.
=20
Regards,
Kana Arunachalam Kannappan
Research Associate
FH O=D6 Forschungs & Entwicklungs GmbH
Stelzhamer Strasse 23,
4600 Wels,
Austria.
Phone: +43 (0)7242 72811 -4420
kana.arunachalam at fh-wels.at
www.fh-ooe.at; www.3dct.at
=20
------_=_NextPart_001_01CB0EC2.6B3BB21F
Content-Type: text/html;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<html xmlns:v=3D"urn:schemas-microsoft-com:vml" =
xmlns:o=3D"urn:schemas-microsoft-com:office:office" =
xmlns:w=3D"urn:schemas-microsoft-com:office:word" =
xmlns:m=3D"http://schemas.microsoft.com/office/2004/12/omml" =
xmlns=3D"http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=3DContent-Type content=3D"text/html; =
charset=3Diso-8859-1">
<meta name=3DGenerator content=3D"Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:266086658;
mso-list-type:hybrid;
mso-list-template-ids:-1900888132 201785359 201785369 201785371 =
201785359 201785369 201785371 201785359 201785369 201785371;}
@list l0:level1
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l1
{mso-list-id:1320773117;
mso-list-type:hybrid;
mso-list-template-ids:-1170698772 201785359 201785369 201785371 =
201785359 201785369 201785371 201785359 201785369 201785371;}
@list l1:level1
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext=3D"edit" spidmax=3D"1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext=3D"edit">
<o:idmap v:ext=3D"edit" data=3D"1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=3DDE-AT link=3Dblue vlink=3Dpurple>
<div class=3DWordSection1>
<p class=3DMsoNormal>Hi Users,<o:p></o:p></p>
<p class=3DMsoNormal><o:p> </o:p></p>
<p class=3DMsoNormal><span lang=3DEN-GB>My task :<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB><o:p> </o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>I have a tube structure and i =
calculate the
hessian matrix and do eigen analysis.<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>For example: =
<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>point P =3D =
[7,7,7]<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Eigen values: EValueMatrix =3D =
[-634, -1389,
-9] (approximate values)<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Eigen vector: EVectorMatrix =
=3D=A0 0.99, 0.13,
0<o:p></o:p></span></p>
<p class=3DMsoNormal><span =
lang=3DEN-GB>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =A00.13,
0.99, 0<o:p></o:p></span></p>
<p class=3DMsoNormal><span =
lang=3DEN-GB>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0 =A00.11,
0.11, 0.78 <o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>From above i say: =
<o:p></o:p></span></p>
<p class=3DMsoListParagraph style=3D'text-indent:-18.0pt;mso-list:l1 =
level1 lfo1'><![if !supportLists]><span
lang=3DEN-GB><span style=3D'mso-list:Ignore'>1.<span style=3D'font:7.0pt =
"Times New Roman"'>
</span></span></span><![endif]><span lang=3DEN-GB>Vector corresponding =
to EValue
=3D -9 which is EVector =3D (0.11, 0.11, 0.78) show the direction =
of<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>The tube structure. =
<o:p></o:p></span></p>
<p class=3DMsoListParagraph style=3D'text-indent:-18.0pt;mso-list:l1 =
level1 lfo1'><![if !supportLists]><span
lang=3DEN-GB><span style=3D'mso-list:Ignore'>2.<span style=3D'font:7.0pt =
"Times New Roman"'>
</span></span></span><![endif]><span lang=3DEN-GB>Vector corresponding =
to EValue
=3D -634 and -1389 which are EVector =3D (0.99, 0.13, 0) and (0.13, =
0.99, 0)<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>respectively represent the plane
perpendicular to the tube direction. (I name the plane as =
Plane1)<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB><o:p> </o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Goal : I want to extract the =
neighbour points
and their gray values of point P and lie on plane Plane1. =
<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB><o:p> </o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>I thought about the below =
pipeline:<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Step 1: Read image =
<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Step 2: Rotate image using =
EVector =3D (0.11,
0.11, 0.78) (using itkResampleImageFilter, rotation transform, trilinear
interpolator)<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Step 3: Rotate image using =
EVector =3D (0.99,
0.13, 0) (using itkResampleImageFilter, rotation transform, trilinear
interpolator)<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Step 4: Rotate image using =
EVector =3D (0.13,
0.99, 0) (using itkResampleImageFilter, rotation transform, trilinear
interpolator)<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Step 5: In the output image from =
step 4 =A0use
neighbourhood iterator to obtain the neighbour =
points.<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB><o:p> </o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>What i need from ITK =
users:<o:p></o:p></span></p>
<p class=3DMsoListParagraph style=3D'text-indent:-18.0pt;mso-list:l0 =
level1 lfo2'><![if !supportLists]><span
lang=3DEN-GB><span style=3D'mso-list:Ignore'>1.<span style=3D'font:7.0pt =
"Times New Roman"'>
</span></span></span><![endif]><span lang=3DEN-GB>I request the itk =
community to
check whether above pipeline is correct to achieve my =
goal.<o:p></o:p></span></p>
<p class=3DMsoListParagraph style=3D'text-indent:-18.0pt;mso-list:l0 =
level1 lfo2'><![if !supportLists]><span
lang=3DEN-GB><span style=3D'mso-list:Ignore'>2.<span style=3D'font:7.0pt =
"Times New Roman"'>
</span></span></span><![endif]><span lang=3DEN-GB>If you have any other =
pipeline
in mind it would be useful for me.<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB><o:p> </o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Thank you in =
advance.<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB><o:p> </o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Regards,<o:p></o:p></span></p>
<p class=3DMsoNormal><span lang=3DEN-GB>Kana Aru</span>nachalam =
Kannappan<o:p></o:p></p>
<p class=3DMsoNormal>Research Associate<o:p></o:p></p>
<p class=3DMsoNormal>FH O=D6 Forschungs & Entwicklungs =
GmbH<o:p></o:p></p>
<p class=3DMsoNormal>Stelzhamer Strasse 23,<o:p></o:p></p>
<p class=3DMsoNormal>4600 Wels,<o:p></o:p></p>
<p class=3DMsoNormal>Austria.<o:p></o:p></p>
<p class=3DMsoNormal>Phone: +43 (0)7242 72811 -4420<o:p></o:p></p>
<p class=3DMsoNormal>kana.arunachalam at fh-wels.at<o:p></o:p></p>
<p class=3DMsoNormal>www.fh-ooe.at; www.3dct.at<o:p></o:p></p>
<p class=3DMsoNormal><o:p> </o:p></p>
</div>
</body>
</html>
------_=_NextPart_001_01CB0EC2.6B3BB21F--
More information about the Insight-users
mailing list