Colour Upper Stereographic Projection Using IPF colouration

[angus-hammond]

Hi, I have a project where it would be really valuable to display an upper stereo-graphic projection where data is coloured based on the standard IPF colour key and data-points are sized according to grain size. I see in the MTEX documentation the code to achieve this for an inverse pole figure is:

plotIPDF(big_grains('olivine').meanOrientation,colors,vector3d.Z,...
  'MarkerSize',0.05*big_grains('olivine').area,'markerEdgeColor','k')

However after some fiddling, I have only been able to plot grain orientations scaled to grain area and I am unable to apply the colour key:

My current code is as follows:

run("EBSD9.m");

% Define the crystal shape for visualization
cS = crystalShape.hex(ebsd.CS);
close all;
plot(cS, 'faceAlpha', 0.2);

% 1. Reconstruct the grains
[grains, ebsd.grainId] = calcGrains(ebsd('indexed'));

% 2. Remove very small grains
ebsd(grains(grains.grainSize < 5)) = [];

% 3. Redo grain reconstruction
[grains, ebsd.grainId] = calcGrains(ebsd('indexed'));

% 4. Plot the grain boundaries
plot(grains.boundary, 'lineWidth', 1.5, 'micronbar', 'off');

% 5. Select only very large grains
big_grains = grains(grains.grainSize > 150);

% 6. Plot the crystals
hold on;
plot(big_grains('calcite'), 0.8 * cS, 'linewidth', 2, 'colored');
hold off;
legend off;

close all;
ipfKey = ipfHSVKey(ebsd('calcite'));
plot(ipfKey, '3d');

plot(grains.boundary, 'lineWidth', 1.5, 'micronbar', 'off');

hold on;
plot(big_grains('calcite'), ipfKey);
hold off;
legend off;

plot(ipfKey, 'complete', 'upper');

h = Miller({1, 0, -1, 4}, {2, 1, -3, 4}, ebsd('calcite').CS);
annotate(h.symmetrise, 'labeled', 'backgroundColor', 'w');

% Convert grain orientations to colors
colors = ipfKey.orientation2color(big_grains('calcite').meanOrientation);

% Define Miller indices for the crystallographic direction(s)
h = Miller(0, 0, 0, 1, big_grains('calcite').CS); % Example: [100] direction

% Plot the pole figure
figure;
plotPDF(big_grains('calcite').meanOrientation, h, 'stereo', ...
    'MarkerSize', 0.05 * big_grains('calcite').area, ...
    'markerEdgeColor', 'k', 'color', colors);
colorbar;

[kilir]

Hi,
the color argument should go second.
E.g.

mtexdata titanium;
grains = ebsd.calcGrains;
h=Miller(0,0,1,ebsd('t').CS);

gs = 0.03* grains('t').area;

ck=ipfHSVKey(ebsd('t').CS);
ck.inversePoleFigureDirection=vector3d(1,0,0);
rgb=ck.orientation2color(grains('t').meanOrientation);

plotIPDF(grains('t').meanOrientation, rgb, ck.inversePoleFigureDirection, 'MarkerSize', gs);
nextAxis
plotPDF(grains('t').meanOrientation, rgb, h, 'MarkerSize', gs);

Cheers,

Rüdiger