Combined Plots edit page

Explains how to combine several plots, e.g. plotting on the top of an inverse pole figure some important crystal directions.

General Principle

In order to tell MATLAB to plot one plot right on the top of an older plot one has to use the commands hold all and hold off. Let's demonstrate this using a simple example.

plot([2 2])

hold all

plot([1 3])

hold off

Combine Different EBSD Data

First, we want to show up two different EBSD data sets in one plot

let's simulate some EBSD data

cs = crystalSymmetry('-3m');
odf = unimodalODF(orientation.byEuler(0,0,0,cs));
ori = calcOrientations(odf,100);
ori_rotated = calcOrientations(rotate(odf,rotation.byEuler(60*degree,60*degree,0*degree)),100);

plot them as a scatter plot in axis / angle parametrized orientation space

hold all
hold off

a second way would be to superpose the pole figures of both EBSD data sets.

h = [Miller(0,0,0,1,cs),Miller(1,0,-1,0,cs)];
hold all % keep plot
hold off % next plot command deletes all plots

Combine countoured pole figures (smooth ODF plots) with EBSD Data Scatter Plots

You can also combine a contour plot of a model ODF with a scatter plot of single orientations.

mtexColorMap white2black

plot(ori,'DisplayName','EBSD 1',...

plot(ori_rotated,'DisplayName','EBSD 2',...


and, of course, you can do the same with ODF plots:

mtexColorMap white2black

Add Miller Indices to an Inverse Pole Figure Plot

Next, we are going to add some Miller indices to an inverse pole figure plot.

mtexColorMap white2black

hold all % keep plot
hold off % next plot command deletes all plots

Combining different plots in one figure

The next example demonstrates how to arrange arbitrary plots into one figure

% let us import some pole figure data
mtexdata dubna
pf = PoleFigure
  crystal symmetry : Quartz (321, X||a*, Y||b, Z||c*)
  specimen symmetry: 1
  h = (02-21), r = 72 x 19 points
  h = (10-10), r = 72 x 19 points
  h = (10-11)(01-11), r = 72 x 19 points
  h = (10-12), r = 72 x 19 points
  h = (11-20), r = 72 x 19 points
  h = (11-21), r = 72 x 19 points
  h = (11-22), r = 72 x 19 points

next, we compute an ODF out of them

odf = calcODF(pf)
odf = ODF (Quartz → xyz)
  Radially symmetric portion:
    kernel: de la Vallee Poussin, halfwidth 5°
    center: 19836 orientations, resolution: 5°
    weight: 1

now we want to plot the original data alongsite with the recalculated pole figures and with a difference plot

figure('position',[50 50 1200 500])

% set position 1 in a 1x3 matrix as the current plotting position
axesPos = subplot(1,3,1);

% plot pole figure 1 at this position

% set position 2 in a 1x3 matrix as the current plotting position
axesPos = subplot(1,3,2);

% plot the recalculated pole figure at this position

% set position 3 in a 1x3 matrix as the current plotting position
axesPos = subplot(1,3,3);

% plot the difference pole figure at this position