Area Fraction with a given orientation

[grandr...@gmail.com]

Hello, I have an EBSD map, and I want to determine the area fraction that is within 15 degrees of the most common orientation.  I'm just missing the last step;

cs=ebsd ('Magnesium').CS

%first calculate the odf to generate continuous data.

odf=calcODF(ebsd('Magnesium').orientations,'halfwidth',5*degree);

%select the maximum from the odf, and load the 3d vector into the

% variable 'modes', and skip return of 'values' by using ~ in that place

[modes, ~] = calcModes(odf,1);

%transform the mode output into the appropriate crystallographic vector

dir = modes(1).*Miller(0,0,0,1,cs);

I don't want to calculate fiber volumes, because those are a percentage of the odf.  The end point of this exercise is to be able to quantitatively compare two maps and determine which has more scatter away from the single major texture peak.  If one map has 95% area fraction within 15 degrees of the mode, and the other has only 80% then obviously the second EBSD map had higher scatter.  If anyone can suggest a better way to do this, please let me know.

Jessica

[Rüdiger Kilian]

Hi,
in case you want it wrt to a reference orientation, how about (in my example for phase 'q':
oref=orientation('Euler',[0,0,0]*degree,cs) % some reference orientation
m=oref\ebsd('q').orientations
in=m.angle<30*degree;
frac=sum(in')/length(ebsd('q'))*100

Your "dir" is a vector (axis) which corresponds to a crystal direction in specimen space, if you want the percentage of all those points that share this vector within a cone of some width, you could do:   
cdir=ebsd('q').orientations.*Miller(0,0,0,1,cs)
refvector=vector3d(0,0,1) %some reference vector
in=angle(cdir,refvector,'antipodal')< 15*degree;
frac=sum(in)/length(ebsd('q'))*100

I hope this helps.
Cheers,
Rüdiger

---

From: mtex...@googlegroups.com [mtex...@googlegroups.com] on behalf of grandr...@gmail.com [grandr...@gmail.com]
Sent: Saturday, May 21, 2016 5:18 AM
To: MTEX
Subject: {MTEX} Area Fraction with a given orientation

--
If you want to reduce the number of emails you get through this forum login to https://groups.google.com/forum/?fromgroups=#!forum/mtexmail, click "My membership" and select "Don't send me email updates". You can still get emails on selected topics by staring them.
---
You received this message because you are subscribed to the Google Groups "MTEX" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mtexmail+u...@googlegroups.com.
Visit this group at https://groups.google.com/group/mtexmail.
For more options, visit https://groups.google.com/d/optout.

[grandr...@gmail.com]

Thank you, this is exactly what I needed.

Jessica