This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
evaluation:plotting-grain-orientations [2019/12/04 11:45] smerkel |
evaluation:plotting-grain-orientations [2019/12/05 15:37] (current) smerkel |
||
---|---|---|---|
Line 20: | Line 20: | ||
Here is an example of a script that will load these orientations and plot them in a pole figures. The example is for olivine and can be adjusted to your needs. | Here is an example of a script that will load these orientations and plot them in a pole figures. The example is for olivine and can be adjusted to your needs. | ||
- | <code> | + | <code matlab plotGrainPF.m> |
%% | %% | ||
% Define crystal and sample symmetry | % Define crystal and sample symmetry | ||
Line 60: | Line 60: | ||
{{ :processing:resultsevalutation:grains-indexing-1.jpg?direct&600 |}} | {{ :processing:resultsevalutation:grains-indexing-1.jpg?direct&600 |}} | ||
- | <code> | + | <code matlab plotColorGrainPF.m> |
%% | %% | ||
% Define crystal and sample symmetry | % Define crystal and sample symmetry | ||
Line 92: | Line 92: | ||
% Plot orientations, colored based on ODF | % Plot orientations, colored based on ODF | ||
plotPDF(ori0,eval(odf,ori0),h,'antipodal','MarkerSize',10, 'MarkerEdgeColor', 'b') | plotPDF(ori0,eval(odf,ori0),h,'antipodal','MarkerSize',10, 'MarkerEdgeColor', 'b') | ||
+ | CLim(gcm,[0,10]) | ||
+ | mtexColorbar | ||
+ | |||
+ | |||
+ | </code> | ||
+ | |||
+ | |||
+ | ===== Loading and plotting grain orientations in MTeX, in an inverse pole figure, and with colors! ===== | ||
+ | |||
+ | Here is an example of a script that will load these orientations and plot them in a inverse pole figures. The example is for olivine and can be adjusted to your needs. Grains are assigned a color based on the ODF intensity for their orientations. | ||
+ | |||
+ | {{ :processing:resultsevalutation:ipf-grains-indexing-1.jpg?direct&400 |}} | ||
+ | |||
+ | <code matlab plotColorGrainIPF.m> | ||
+ | %% | ||
+ | % Define crystal and sample symmetry | ||
+ | cs_ol = crystalSymmetry('Pnma',[9.915267 5.886196 4.713400]) | ||
+ | ss = specimenSymmetry() | ||
+ | % Plotting convention | ||
+ | plotx2east | ||
+ | % File with Euler angles | ||
+ | fname = fullfile('/path/to/folder/',... | ||
+ | '','euler.txt'); | ||
+ | |||
+ | %% | ||
+ | % Load file and create the corresponding list of orientations | ||
+ | ori = load(fname) | ||
+ | n = size(ori); | ||
+ | n = n(1); | ||
+ | |||
+ | for i = 1:n | ||
+ | phi1 = ori(i,1)*degree; | ||
+ | Phi = ori(i,2)*degree; | ||
+ | phi2 = ori(i,3)*degree; | ||
+ | ori0(i) = orientation('Euler',phi1,Phi,phi2,cs_ol,ss); | ||
+ | end | ||
+ | |||
+ | %% | ||
+ | % Plots grains using a color scale base on ODF fitted to the grains | ||
+ | |||
+ | % Fit an ODF to the data | ||
+ | odf = calcODF(ori0,15*degree) | ||
+ | |||
+ | |||
+ | % Plot grains as IPF of the compression (if compression is in X), colored based on ODF | ||
+ | plotIPDF(ori0,eval(odf,ori0),[vector3d.X],'antipodal','MarkerSize',10, 'MarkerEdgeColor', 'b') | ||
CLim(gcm,[0,10]) | CLim(gcm,[0,10]) | ||
mtexColorbar | mtexColorbar | ||
Line 107: | Line 153: | ||
{{ :processing:resultsevalutation:indexations-1-2-3.jpg?direct&600 |}} | {{ :processing:resultsevalutation:indexations-1-2-3.jpg?direct&600 |}} | ||
- | <code> | + | <code matlab compareGrainPF.m> |
%% | %% | ||
% Define crystal and sample symmetry | % Define crystal and sample symmetry |