Differences

This shows you the differences between two versions of the page.

--- evaluation:plotting-grain-orientations [2019/12/03 14:52]
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.
-<code>
+<code matlab plotGrainPF.m>
 %%
 % Define crystal and sample symmetry
@@ Line 48: / Line 48: @@
 h = Miller({1,0,0},{0,1,0},{0,0,1},cs_ol)
 plotPDF(ori0,h,'antipodal','MarkerSize',10, 'MarkerFaceColor', 'b', 'MarkerEdgeColor', 'b')
+</code>
+===== Loading and plotting grain orientations in MTeX, with colors! =====
+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.
+In this case grains are assigned a color based on the ODF intensity for their orientations.
+{{ :processing:resultsevalutation:grains-indexing-1.jpg?direct&600 |}}
+<code matlab plotColorGrainPF.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 orientations, colored based on ODF
+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])
+mtexColorbar
 </code>
@@ Line 60: / Line 153: @@
 {{ :processing:resultsevalutation:indexations-1-2-3.jpg?direct&600 |}}
-<code>
+<code matlab compareGrainPF.m>
 %%
 % Define crystal and sample symmetry

TIMEleSS Multigrain Wiki

User Tools

Site Tools

Differences

Page Tools