<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:st1="urn:schemas-microsoft-com:office:smarttags" xmlns="http://www.w3.org/TR/REC-html40">

<head>
<meta http-equiv=Content-Type content="text/html; charset=Windows-1254">
<meta name=Generator content="Microsoft Word 11 (filtered medium)">
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="State"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="PlaceName"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="PlaceType"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="place"/>
<!--[if !mso]>
<style>
st1\:*{behavior:url(#default#ieooui) }
</style>
<![endif]-->
<style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:11;
        panose-1:0 0 0 0 0 0 0 0 0 0;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman";}
h4
        {margin-top:12.0pt;
        margin-right:0cm;
        margin-bottom:3.0pt;
        margin-left:43.2pt;
        text-indent:-43.2pt;
        page-break-after:avoid;
        mso-list:l0 level4 lfo1;
        font-size:14.0pt;
        font-family:"Times New Roman";}
a:link, span.MsoHyperlink
        {color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {color:purple;
        text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman";}
p.Style2, li.Style2, div.Style2
        {margin-top:12.0pt;
        margin-right:0cm;
        margin-bottom:3.0pt;
        margin-left:43.2pt;
        text-indent:-43.2pt;
        page-break-after:avoid;
        mso-list:l0 level4 lfo1;
        font-size:12.0pt;
        font-family:Arial;
        font-weight:bold;}
@page Section1
        {size:595.3pt 841.9pt;
        margin:72.0pt 69.65pt 72.0pt 69.65pt;}
div.Section1
        {page:Section1;}
 /* List Definitions */
 @list l0
        {mso-list-id:497774475;
        mso-list-template-ids:1719176970;}
@list l0:level1
        {mso-level-text:%1;
        mso-level-tab-stop:21.6pt;
        mso-level-number-position:left;
        margin-left:21.6pt;
        text-indent:-21.6pt;}
@list l0:level2
        {mso-level-text:"%1\.%2";
        mso-level-tab-stop:28.8pt;
        mso-level-number-position:left;
        margin-left:28.8pt;
        text-indent:-28.8pt;}
@list l0:level3
        {mso-level-text:"%1\.%2\.%3";
        mso-level-tab-stop:36.0pt;
        mso-level-number-position:left;
        margin-left:36.0pt;
        text-indent:-36.0pt;}
@list l0:level4
        {mso-level-style-link:"Heading 4";
        mso-level-text:"%1\.%2\.%3\.%4";
        mso-level-tab-stop:43.2pt;
        mso-level-number-position:left;
        margin-left:43.2pt;
        text-indent:-43.2pt;}
@list l0:level5
        {mso-level-text:"%1\.%2\.%3\.%4\.%5";
        mso-level-tab-stop:50.4pt;
        mso-level-number-position:left;
        margin-left:50.4pt;
        text-indent:-50.4pt;}
@list l0:level6
        {mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6";
        mso-level-tab-stop:57.6pt;
        mso-level-number-position:left;
        margin-left:57.6pt;
        text-indent:-57.6pt;}
@list l0:level7
        {mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6\.%7";
        mso-level-tab-stop:64.8pt;
        mso-level-number-position:left;
        margin-left:64.8pt;
        text-indent:-64.8pt;}
@list l0:level8
        {mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6\.%7\.%8";
        mso-level-tab-stop:72.0pt;
        mso-level-number-position:left;
        margin-left:72.0pt;
        text-indent:-72.0pt;}
@list l0:level9
        {mso-level-text:"%1\.%2\.%3\.%4\.%5\.%6\.%7\.%8\.%9";
        mso-level-tab-stop:79.2pt;
        mso-level-number-position:left;
        margin-left:79.2pt;
        text-indent:-79.2pt;}
ol
        {margin-bottom:0cm;}
ul
        {margin-bottom:0cm;}
-->
</style>

</head>

<body lang=EN-GB link=blue vlink=purple>

<div class=Section1>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>hi, here's a follow up on the colour ellipsoids.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>i thought of the following procedure after the two rounds of
discussions. next step is to identify all the objects to draw the colour
ellipsoid. <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><b><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11;font-weight:bold'># 1. read STRUCTURED_POINTS
data file<o:p></o:p></span></font></b></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'>vtkStructuredPointsReader
reader<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'>            reader
SetFileName &quot;tial.vtk&quot;<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'>            reader Update;<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><b><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11;font-weight:bold'># 2. extract
point tensor data<o:p></o:p></span></font></b></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'>vtkStructuredPoints tialdata<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-indent:36.0pt;text-autospace:none'><font size=2
face=11><span lang=EN-ZA style='font-size:11.0pt;font-family:11'>set tialdata
[reader GetOutput] <o:p></o:p></span></font></p>

<p class=MsoNormal style='text-indent:36.0pt;text-autospace:none'><font size=2
face=11><span lang=EN-ZA style='font-size:11.0pt;font-family:11'>set myTensors
[tialdata GetPointData]<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><b><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11;font-weight:bold'># 3. create
an array to store (40*40) unit normals <o:p></o:p></span></font></b></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
style='font-size:11.0pt;font-family:11'>vtkDoubleArray</span></font><font
size=2 face=11><span style='font-size:11.0pt;font-family:11'> <span lang=EN-ZA>normals<o:p></o:p></span></span></font></p>

<p class=MsoNormal style='text-indent:36.0pt;text-autospace:none'><font size=2
face=11><span lang=EN-ZA style='font-size:11.0pt;font-family:11'>normals
SetNumberOfComponents 3<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><b><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11;font-weight:bold'># 4. create
a stress vector array by (tensor . unit vectors)<o:p></o:p></span></font></b></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
style='font-size:11.0pt;font-family:11'>vtkDoubleArray</span></font><font
size=2 face=11><span style='font-size:11.0pt;font-family:11'> <span lang=EN-ZA>stressVectors<o:p></o:p></span></span></font></p>

<p class=MsoNormal style='text-indent:36.0pt;text-autospace:none'><font size=2
face=11><span lang=EN-ZA style='font-size:11.0pt;font-family:11'>stressVectors
SetNumberOfComponents 3<o:p></o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal style='text-autospace:none'><b><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11;font-weight:bold'># 5. create
a scalar array from the arcos[(stressVectors . unit vectors)/|stressVectors|]<o:p></o:p></span></font></b></p>

<p class=MsoNormal style='text-autospace:none'><font size=2 face=11><span
style='font-size:11.0pt;font-family:11'>vtkDoubleArray</span></font><font
size=2 face=11><span style='font-size:11.0pt;font-family:11'> <span lang=EN-ZA>stress
angle<o:p></o:p></span></span></font></p>

<p class=MsoNormal style='text-indent:36.0pt;text-autospace:none'><font size=2
face=11><span lang=EN-ZA style='font-size:11.0pt;font-family:11'>angle SetNumberOfComponents
1<o:p></o:p></span></font></p>

<p class=MsoPlainText><b><font size=2 face=11><span style='font-size:11.0pt;
font-family:11;font-weight:bold'><o:p>&nbsp;</o:p></span></font></b></p>

<p class=MsoPlainText><b><font size=2 face=11><span style='font-size:11.0pt;
font-family:11;font-weight:bold'># 6. use angle array to create a scalar lookuptable<o:p></o:p></span></font></b></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><b><font size=2 face=11><span style='font-size:11.0pt;
font-family:11;font-weight:bold'># 7. use sphere geometry to map tensors<o:p></o:p></span></font></b></p>

<p class=MsoPlainText><b><font size=2 face=11><span style='font-size:11.0pt;
font-family:11;font-weight:bold'><o:p>&nbsp;</o:p></span></font></b></p>

<p class=MsoPlainText><b><font size=2 face=11><span style='font-size:11.0pt;
font-family:11;font-weight:bold'># 8. use anglelookuptable to color the sphere.</span></font></b><font
size=2 face=11><span style='font-size:11.0pt;font-family:11'><o:p></o:p></span></font></p>

<p class=MsoPlainText><b><font size=2 face=11><span style='font-size:11.0pt;
font-family:11;font-weight:bold'><o:p>&nbsp;</o:p></span></font></b></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>-------------<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>The data file “</span></font><font size=2 face=11><span
lang=EN-ZA style='font-size:11.0pt;font-family:11'>tial.vtk” is</span></font><font
size=2 face=11><span lang=EN-US style='font-size:11.0pt;font-family:11'><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'># vtk DataFile Version 3.0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>Depth Resolved Stress State of <st1:State w:st="on"><st1:place
 w:st="on">Ti-6Al-4V</st1:place></st1:State><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>ASCII<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>DATASET STRUCTURED_POINTS<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>DIMENSIONS 1 1 3<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>ORIGIN 0 0 0 <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>SPACING 10 10 10 <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>POINT_DATA 3<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>TENSORS stress double<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>57        0          0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>0          146      0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>0          0          89<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>57        0          0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>0          46        0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>0          0          89<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>57        0          0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>0          146      0<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>0          0          -89<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'> <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>as suggested i’ll be working with the vtkSphereSource,  vtkDoubleArray,
vtkMath classes.</span></font><font size=2 face=11><span lang=EN-US
style='font-size:11.0pt;font-family:11'><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>regards,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>mecit yaman<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><st1:place w:st="on"><st1:PlaceType w:st="on"><font
  size=2 face=11><span lang=EN-US style='font-size:11.0pt;font-family:11'>university</span></font></st1:PlaceType><font
 size=2 face=11><span lang=EN-US style='font-size:11.0pt;font-family:11'> of <st1:PlaceName
 w:st="on">cape town</st1:PlaceName></span></font></st1:place><font size=2
face=11><span lang=EN-US style='font-size:11.0pt;font-family:11'><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span lang=EN-US style='font-size:
11.0pt;font-family:11'>-----Original Message-----<br>
From: Ugur BOZKAYA [mailto:u.bozkaya@fmcon.com] <br>
Sent: 02 August 2006 10:47 PM<br>
To: Mecit Yaman<br>
Cc: vtkusers@vtk.org<br>
Subject: Re: [vtkusers] FW: colour mapping stress ellipsoids</span></font><font
size=2 face=11><span style='font-size:11.0pt;font-family:11'><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>Hi,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>I suggest you not to use vtkTensorGlyph. Here is what is
confusing, you <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>need 400 values for each tensor but vtkTensorGlyph accepts only
one for <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>each tensor.(I might be wrong since I did not use it for
sometime) May <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>be it can be manipulated but this will be more complicated.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>You will need to create each tensor by your own. Normally in
vtk, you <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>can just take a vtkSphereSource and define its axes and
resolution to <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>create a polydata that results a tensor elipsoid. One way is
you can <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>modify this class and add a dataarray to this class that will
hold a <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>value for each point in the elipsoid. Anyway if you want to
create your <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>own class vtkSphereSource should guide you.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>About holding the data values; the geometry and data values of
this <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>object is two independent things. Simply if you have a gray
level image, <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>the geometry is a rectangular grid MxN (say 100x200). So you
have MxN <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>points with 3D coordinates such as pnt0(0,0,0), pnt1(1,0,0), <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>pnt2(2,0,0)...pntM-1(M-1,0,0) in the first row.(assuming unit
spacing <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>with 0,0,0 origin). Of course how this points are connected is
important <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>this is done by inserting cells to data. But this is not an
image yet, <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>we have to define data values(grey values). So you create a
dataarray <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>such as vtkDoubleArray. Since this is a grey level image only
one <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>component is necessary for each point. The vtkDoubleArray
should be MxN <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>in size. If this was an RGB image the dataarray would have 3
components. <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>The point and the data values are matched according to their
order. So <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>pnt0 will have the value of first value in vtkDoubleArray and
so on... <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>You do this by setscalars() function.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>In your case you will create a tensor in a specific coordinate.
You have <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>to calculate the location of 400 pnts for this tensor. For each
of these <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>pnts just calculate the color value and add it to your
dataarray. Then <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>you have to define cells to construct the surface.(briefly
create <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>triangles from 3 closest pnts on the surface and add it as a
cell)<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>This whole process will be done for each tensor of course.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>For the calculation of color values, you can use vnl which is
included <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>in itk to allow matlab like operations. Otherwise just write or
find a <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>function that will handle tensor product. dot product and norm
can be <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>found in vtkMath class.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>good luck<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>Mecit Yaman wrote:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; hi, thanks for the replies.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; i need to clarify one thing straight. the coloring i am
trying is not <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; restricted to three orthogonal dimension as ugur assumed.
for each <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; tensor, i need a set of colors calculated for each theta
and phi <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; direction (in the spherical coordinates.) assuming a
thetaResolution <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; of 20 and phiresolution of 20, color scalar array will be
400 scalars. <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; i tried to explain how i get the colors in the first
message. <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; basically a colour is calculated from the surface normal
and vector at <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; that location (the vector is calculated from the tensor by
a scalar <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; product of the tensor and the surface normal) (see my
first message <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; below.) although the visualization idea is complicated, i
know that <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; vtk script should not necessarily be as complex.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; a visualization task with vtk follows<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; source &gt; filters&gt; mapper &gt; actor &gt;rendering<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; i understand i need to scale and aling a vtkTensorGlyph using
my <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; tensor. then i need to create a vtkDoubleArray with 400
components <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; using the same tensor and vtkTensorGlpyh surface normals.
what is not <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; clear to me is: how to calculate the color array and store
them in a <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; vtkDoubleArray from the tensor and surface normals. i
suppose i need <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; to getSurfaceNormals and multiply them component by
component with my <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; tensor to obtain another set of vectors. and with this new
set i need <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; to do scalar multiplication with the surface normals to
get the scalar <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; color array.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; once i have a scalar array i assign them to the
vtkTensorGlyph, as <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; luca wrote, with
vtkPolyData::GetPointData::SetScalars(array).<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; i am familiar with vtk theoretically and used it for
straightforward <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; tasks. i will appreciate any help regarding how to
practically create <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; the color array. btw i’m using java to write the scripts.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; cheers,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; mecit yaman<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; <st1:place w:st="on"><st1:PlaceType w:st="on">university</st1:PlaceType>
 of <st1:PlaceName w:st="on">cape town</st1:PlaceName></st1:place><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;
------------------------------------------------------------------------<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; *From:* Luca Pallozzi Lavorante
[mailto:lplavorante@gmail.com]<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; *Sent:* 25 July 2006 01:55 PM<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; *To:* yaman@myhorizon.org<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; *Subject:* Re: [vtkusers] colour mapping stress ellipsoids<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Yaman, I do agree with Ugur.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; I have used the vtkTensorGlyph class in order to get
strain <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; ellipsoids. I also have tried to color the deformed object
and, as far <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; as I know, the only way is to associate scalar values to
the object´s <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; points through the vtkFloatArray or vtkDoubleArray. Once
you created <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; the array, the association is made via a call like <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; vtkPolyData::GetPointData::SetScalars(array).<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; I hope this helps<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Luca<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; On 7/24/06, *Ugur BOZKAYA* &lt;u.bozkaya@fmcon.com <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &lt;mailto:u.bozkaya@fmcon.com&gt;&gt; wrote:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Selam Mecit,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; As far as I know, there is no class for this purpose in
VTK currently.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; vtkTensorGlyph and vtkGlyph3D can be used for typical
tensor<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; visualization. But the tensors are shown only in one
color(RGB is also<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; possible I guess, but it will not solve your problem). If
you are<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; familiar enough with VTK, you can change these classes for
your purpose.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Those glyph shapes such as ellipsoids (deformed spheres)
are already<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; defined as vtkPolyData, so you can use other shapes.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; I assume you are planning to use 3 orthogonal directions
for coloring.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; You only need a vtkDataArray (e.g. vtkDoubleArray) with 3
components. So<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; for each tensor location you will have 3 values.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; I guess it will be easier to write your own class which
looks simpler.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; You will have your output as vtkPolyData, which is
actually composed of<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; your tensors. May be you can have a function which
calculates a<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; tensorshape(multicolored vtkPolyData) for a given location
then you can<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; add the output of this function to your volume output for
each location.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Good Luck,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Ugur<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Mecit Yaman wrote:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; hi,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; i am trying to visualize second order tensors using
stress ellipsoids <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; with a<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; different color mapping. the usual colouring is to
map a scalar <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; value, e.g.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; components, effective stress, trace or other scalars
on the surface. <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; i want<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; to do a colouring where each point/polygon on the
ellipsoid will be <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; coloured<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; differently. the color scheme is defined as follows:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; the color along direction n(i) of the ellipsoid is
obtained from the <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; angle<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; between the n(i), direction vector, and the stress
vector along n(i). <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; (the<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; stress vector along a direction can be obtained by
multiplying the stress<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; tensor, T(ij) with the direction vector, n(i)). Thus:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; T(ij) tensor<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; n(i) direction vector<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; stress vector(i) = T(ij) . n(i)<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; angle between stress vector and direction vector<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; cos(angle) = stress vector (i) . n(i) / ( |stress
vector(i)| | n(i)|)<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; how does one colour map these angles, which are
different for every<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; direction on an ellipsoid? do I need to change stress
ellipsoid into a<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; polygonal data type? how do I calculate and store the
angles in VTK?<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; i am able to render the icons in Mathematica, and
will be happy to share<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; information with anyone interested.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; thanks,<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; M Yaman<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; <st1:place w:st="on"><st1:PlaceType w:st="on">University</st1:PlaceType>
 of <st1:PlaceName w:st="on">Cape Town</st1:PlaceName></st1:place><o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; _______________________________________________<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; This is the private VTK discussion list.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; Please keep messages on-topic. Check the FAQ at: <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; http://www.vtk.org/Wiki/VTK_FAQ<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; Follow this link to subscribe/unsubscribe:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt; http://www.vtk.org/mailman/listinfo/vtkusers<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; &gt;<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; _______________________________________________<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; This is the private VTK discussion list.<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Please keep messages on-topic. Check the FAQ at: <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; http://www.vtk.org/Wiki/VTK_FAQ<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Follow this link to subscribe/unsubscribe:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; http://www.vtk.org/mailman/listinfo/vtkusers<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; ------------------------------------------------------------------------<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;<o:p>&nbsp;</o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; _______________________________________________<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; This is the private VTK discussion list. <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Please keep messages on-topic. Check the FAQ at:
http://www.vtk.org/Wiki/VTK_FAQ<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; Follow this link to subscribe/unsubscribe:<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt; http://www.vtk.org/mailman/listinfo/vtkusers<o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'>&gt;   <o:p></o:p></span></font></p>

<p class=MsoPlainText><font size=2 face=11><span style='font-size:11.0pt;
font-family:11'><o:p>&nbsp;</o:p></span></font></p>

</div>

</body>

</html>