mtex-toolbox / mtex

MTEX is a free Matlab toolbox for quantitative texture analysis. Homepage:
http://mtex-toolbox.github.io/
GNU General Public License v2.0
281 stars 185 forks source link

Misorientation Axes in Specimen reference frame #115

Closed zmichels closed 8 years ago

zmichels commented 8 years ago

Hello Ralf and others,

This is probably an ignorant question, and hopefully has a simple answer. I am uncertain about how to plot the axes of misorientations relative to the principal directions of the specimen/acquisition reference frame.... rather than the crystal reference frame. I am probably missing something quite basic.

Here is an example:

>> mori = calcMisorientation(ebsd('Olivine'))

mori = misorientation (show methods, plot)
  size: 97907 x 1
  crystal symmetry : Olivine (mmm)
  crystal symmetry : Olivine (mmm)

The following plotting commands all seem to visualize the misorientation data relative to the crystallographic directions:

>> plotAxisDistribution(mori,'smooth')

screen shot 2015-10-28 at 8 12 07 am

>> plot(axis(mori),'antipodal','lower','smooth')

screen shot 2015-10-28 at 8 12 48 am

How can I plot the axes of the misorientations relative to the specimen reference frame (i.e., a lower hemisphere plot of the misorientation axes that is orientated relative to the xvector, yvector, zvector specimen directions)?

Any help is greatly appreciated.

Thank you, Zach


For context... here are some other plots

plotPDF of mori data:

>> plotPDF(mori,ph,'antipodal','lower','complete','smooth')

screen shot 2015-10-28 at 8 24 26 am

ODF for the sample (a mantle xenolith; sample surface cut randomly relative to any mesoscopic fabric): screen shot 2015-10-28 at 8 26 14 am

ralfHielscher commented 8 years ago

Hi Zach,

given a misorientation it is impossible to retrieve the misorientation axis in specimen coordinates. To compute this axis we need a pair of orientations. A list of such pairs of orientations can be computed by the command

[ori1,ori2] = calcMisorientation(ebsd('Olivine'))

Not the misorientation axis in specimen coordinates is computed by

a = axis(ori1,ori2)

You can plot the axes now by

plot(a,'contourf')

In the next version of MTEX 4.2 you can do this directly with plotAxisDistribution.

Ralf.

zmichels commented 8 years ago

Ralf,

Thank you much!

best, Zach