generate IPF colormap

628 views
Skip to first unread message

Alec Day

unread,
Sep 2, 2014, 10:26:03 PM9/2/14
to mtex...@googlegroups.com
Hello

Is it possible to use mtex to produce an inverse pole figure colormap (see: http://www.uni-ulm.de/typo3temp/pics/602211b629.jpg) and to automatically assign a colour based upon the miller index of the crystal in the specimen reference frame.

Thank you,

Alec Day

Ralf Hielscher

unread,
Sep 3, 2014, 12:29:34 AM9/3/14
to mtex...@googlegroups.com
This is what MTEX up to 3.5 does automatically if you plot EBSD data. To see the inverse pole figure triangle simple type

plot(ebsd)

colorbar

Starting with MTEX 4.0 this has to be coded more exhaustive, i.e., one has to explicitly specify an orientation mapping:

oM = ipdfHSVOrientationMapping(ebsd('phaseName'))

colors
= oM.orientation2color(ebsd('phaseName').orientations);

plot
(ebsd('phaseName'),colors)


I hope this helps,

Ralf.

Alec Day

unread,
Sep 3, 2014, 10:50:00 AM9/3/14
to mtex...@googlegroups.com
Hi Ralf,

Sorry I don't think that is quite what I want (at least i don't think it is :p), as I'm not using EBSD data.

What I want is the ability to generate the inverse pole figure color map, and say if I had a crystal plane of {111}, the code would output the color associated from the location of {111} on the IPF colormap, i.e. blue from the one linked (red for {001} and green for {101}).

So given any orientation i'd like the position in the IPF along with its color defined by the RGB vector at that point in the colormap

Ralf Hielscher

unread,
Sep 22, 2014, 3:00:34 AM9/22/14
to mtex...@googlegroups.com
Hi Alec Day,

maybe in MTEX 4.0 :)

After defining an orientation to color mapping for a specific crystal symmetry by
cs = crystalSymmetry('your symmetry')
oM
= ipdfHSVOrientationMapping(cs)

you can plot the colorized triangle by 
plot(oM)

to get the rgb values for some orientations ori do
rgb = oM.orientation2color(ori);

to get the position in the inverse pole figure triangle do
h = inv(ori) * oM.inversePoleFigureDirection
h
= h.project2FundamentalRegion

I hope this helps,
Ralf.


On Wednesday, September 3, 2014 4:26:03 AM UTC+2, Alec Day wrote:
Reply all
Reply to author
Forward
0 new messages