How to get local. max. values on a misorientation axis plot?

[Laszlo Aradi]

Hello everyone,

is there a way to quantify the local max. densities on a misorentation axis distribution plot? In my case I'd like to quantify which rot. axis is more dominant in forsterite, because my samples show smooth transition between [001] and [010] as rotational axis for low angle subgrain boundaries (e.g. max at [001]/max at [010]). I can count the isolines manually (btw 'contour' always draws the multiplies of MUD?), but would be more elegant, precise and less time consuming, if I could implement it in my script. My only idea was to draw a  'circle' at 010 and 001 and somehow get the max. intensity within it, but unfortunately i could not make it work.

The basic script i'm using is coming from David's workshop. Thank you for your help!

cheers

Laszlo

edit: attached a wrong figure

%% Misorientation axes for specific angles 2 to 10 degrees forsterite**************************************************************************

close all

lower_limit = 2

upper_limit = 10

% logical condition

select_angle_range = gb_m_angle_cor >= lower_limit & gb_m_angle_cor < upper_limit;

% select angles

selected_angles = gb_m_angle_cor(select_angle_range);

% select axes using angle range

selected_axes = gb_m_axis_cor(select_angle_range);

N_axes = length(selected_axes)

figure

% plots whole upper hemisphere

plot(selected_axes,'contourf','upper','antipodal','fundamentalRegion')

colorbar

% add crystallographic axes

annotate([Miller(1,0,0,ebsd('Forsterite').CS,'uvw'),...

          Miller(0,1,0,ebsd('Forsterite').CS,'uvw'),...

          Miller(0,0,1,ebsd('Forsterite').CS,'uvw')],...

         'all','labeled','BackgroundColor','w')

% save file

set(gcf,'units','normalized','outerposition',[0.1 0.1 0.5 0.5])

view([90 90])

export_fig GN1411_For_For_MDF_gb_Misorientation_angle_2_10_degrees

[David Mainprice]

Dear Laszlo,

While you have the plot active in Matlab/MTEX you click on any point it will give you azimuth, polar angle and density in multiples of uniform distribution, something like 6.5 at (85,5) where 6.5 x uniform, 85 degrees polar angle from plot centre and 5 azimuth anti-clockwise from X.

all the best David

Le 21 sept. 2015 à 14:32, Laszlo Aradi <aradi.la...@gmail.com> a écrit :

Hello everyone,

is there a way to quantify the local max. densities on a misorentation axis distribution plot? In my case I'd like to quantify which rot. axis is more dominant in forsterite, because my samples show smooth transition between [001] and [010] as rotational axis for low angle subgrain boundaries (e.g. max at [001]/max at [010]). I can count the isolines manually (btw 'contour' always draws the multiplies of MUD?), but would be more elegant, precise and less time consuming, if I could implement it in my script. My only idea was to draw a  'circle' at 010 and 001 and somehow get the max. intensity within it, but unfortunately i could not make it work.

The basic script i'm using is coming from David's workshop. Thank you for your help!

cheers

Laszlo

%% Misorientation axes for specific angles 2 to 10 degrees forsterite**************************************************************************

close all

lower_limit = 2

upper_limit = 10

% logical condition

select_angle_range = gb_m_angle_cor >= lower_limit & gb_m_angle_cor < upper_limit;

% select angles

selected_angles = gb_m_angle_cor(select_angle_range);

% select axes using angle range

selected_axes = gb_m_axis_cor(select_angle_range);

N_axes = length(selected_axes)

figure

% plots whole upper hemisphere

plot(selected_axes,'contourf','upper','antipodal','fundamentalRegion')

colorbar

% add crystallographic axes

annotate([Miller(1,0,0,ebsd('Forsterite').CS,'uvw'),...

          Miller(0,1,0,ebsd('Forsterite').CS,'uvw'),...

          Miller(0,0,1,ebsd('Forsterite').CS,'uvw')],...

         'all','labeled','BackgroundColor','w')

% save file

set(gcf,'units','normalized','outerposition',[0.1 0.1 0.5 0.5])

view([90 90])

export_fig GN1411_For_For_MDF_gb_Misorientation_angle_2_10_degrees

--
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 http://groups.google.com/group/mtexmail.
For more options, visit https://groups.google.com/d/optout.
<GN1411_Dio_Dio_MDF_gb_Misorientation_angle_2_10_degrees.png>

David Mainprice

15 Les Romarins

34270 Les Matelles

Tel. +33 4 67 67 03 30

(In France 04 67 67 03 30)

E-mail: davidma...@icloud.com

[Laszlo Aradi]

Dear David,

thank you for your answer, i totally forgot about the data cursor mode. However, i'm still curious if there are anyway to do this automatically from predetermined positions ( e.g. [0 0]) of the plot (manually it can be done by right clicking the data tip and Export Cursor Data to Workspace), since i plan to use this on a huge amount of data and figure.

cheers

Laszlo

[David Mainprice]

Dear Laszlo,

You can use [modes, values] = calcModes(odf,n).  modes are the Euler angles of the maxima in the ODF or MDF and values of the maxima in x uniform, n is number of maxima you want to search.

I think you can use this on an MDF as well.

You can then you can used modes to annotate pole figures, inverse pole figures etc.

annotate(modes,'marker','s','MarkerFaceColor','black’)

all the best David

[Laszlo Aradi]

Dear David,

thank you again for your help and kindness.The method works on only 'Miller type' arguments, and mathlab tells me selected_axes (which i'd liket to plot) is not that kind:

Undefined function 'calcModes' for input arguments of type 'Miller'.

However according to Mathlab's workspace selected_axes is a 5581x1 big Miller type matrix.

cheers

Laszlo

[Laszlo Aradi]

A bit different problem, but i didn't want to open another topic. Is it possible to move "out" the annotations from the plot (so in this case the [001] and [010] annotations does not cover the plot), unfortunately 'northoutside' and other parameters does not seem to be working.

cheers

Laszlo