%==========================================================================
%==========================================================================
%
%              /$$$$$$$$ /$$$$$$$   /$$$$$$  /$$$$$$$                       
%             | $$_____/| $$__  $$ /$$__  $$| $$__  $$                      
%             | $$      | $$  \ $$| $$  \__/| $$  \ $$                      
%             | $$$$$   | $$$$$$$ |  $$$$$$ | $$  | $$                      
%             | $$__/   | $$__  $$ \____  $$| $$  | $$                      
%             | $$      | $$  \ $$ /$$  \ $$| $$  | $$                      
%             | $$$$$$$$| $$$$$$$/|  $$$$$$/| $$$$$$$/                      
%             |________/|_______/  \______/ |_______/                       
%                                                                           
%                                                                           
%                                                                           
%                                                    /$$                    
%                                                   |__/                    
%   /$$$$$$  /$$   /$$  /$$$$$$   /$$$$$$   /$$$$$$$ /$$  /$$$$$$$  /$$$$$$ 
%  /$$__  $$|  $$ /$$/ /$$__  $$ /$$__  $$ /$$_____/| $$ /$$_____/ /$$__  $$
% | $$$$$$$$ \  $$$$/ | $$$$$$$$| $$  \__/| $$      | $$|  $$$$$$ | $$$$$$$$
% | $$_____/  >$$  $$ | $$_____/| $$      | $$      | $$ \____  $$| $$_____/
% |  $$$$$$$ /$$/\  $$|  $$$$$$$| $$      |  $$$$$$$| $$ /$$$$$$$/|  $$$$$$$
%  \_______/|__/  \__/ \_______/|__/       \_______/|__/|_______/  \_______/
%                                                                           
%                                                                           
%                                                                           
%                     /$$$$$$   /$$$$$$   /$$$$$$   /$$$$$$                 
%                    /$$__  $$ /$$$_  $$ /$$__  $$ /$$$_  $$                
%                   |__/  \ $$| $$$$\ $$|__/  \ $$| $$$$\ $$                
%                     /$$$$$$/| $$ $$ $$  /$$$$$$/| $$ $$ $$                
%                    /$$____/ | $$\ $$$$ /$$____/ | $$\ $$$$                
%                   | $$      | $$ \ $$$| $$      | $$ \ $$$                
%                   | $$$$$$$$|  $$$$$$/| $$$$$$$$|  $$$$$$/                
%                   |________/ \______/ |________/ \______/                 
%                                                                           
%                                                                   
%                                                                           
%==========================================================================
%==========================================================================
 
% The aim of this exercise is to get familiar with mtex and look at some ebsd
% data

% Some general Notes:
% - usually you can find help by typing `help` command
%   if a command with an identical name might be available for different
%   data types try e.g. `help orientation/calcMDF` or `help odf/calcMDF` or
%   help orientation/mean` (`help mean` will be Matlabs number-valued mean)
% - using TAB-completion makes many tasks faster, prevents typos and helps
%   to explore functions and properties available to an object
%   e.g. ebsd.[hit the tabulator key] will present a list of features
%  available with the ebsd object
%% install/startup mtex (in case it isn't already)


%% 0) load the data -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
% we can use the import_wizard and generate a script which imports the ebsd
% data. Use the file ebsd1.ctf


import_wizard
% selcte EBSD and click Next ...
% when asked for the reference frame, x should be east and y south when
% using convertEuler2SpatialReferenceFrame


%% 1) explore the ebsd data *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
%  What can you tell about the ebsd data?
% hint: if you use ebsd.[TAB] you will see some properties and functions
% available for this object


% %########################################################################
% a) How many indexed pixels does the map have?
% %########################################################################


% %########################################################################
% b) Which properties does the ebsd map have?
% %########################################################################

% eg.g
%       mad = mean angular deviation:misfit between detected bands and
%       simulation (in degree, high mad = 'bad' solution)
%       bc = band contrast: ratio of actual intensity in hough space vs
%       average intensity of kikuchi bands (high value = high contrast)
%       bs =band slope : higest intensity gradient at the morgin of kikuchi
%       bands (high value = sharp band)


% %########################################################################
% c) How many percent of the map have pixels with a mad < 0.8 degree? 
% %########################################################################


% %########################################################################
% d) How many pixels are indexed within the first 100 points?
% %########################################################################


% %########################################################################
% e) How large is the map? How many pixels does it have in x,y direction?
% %########################################################################


% %########################################################################
% f) What is the stepsize of the map?
% %########################################################################


%% % 2) look at the map -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*

% use the plot command to plot almost anything in mtex
plot(ebsd)
% maybe not so spectacular at first sight, it will just plot the phase

% mtex will always try to honor the convention
% plot(where,what)
% in case of EBSD data:
% "where" are x,y coordiantes of each measurement point
% "what" can be anything as long as it has the same size/length as the first argument.

% %########################################################################
% Question:
% a) Colorcode the band contrast in a grey-scale color-coding
% hint1: when plotting (where,what) values, both need to have the same
% length
% hint2: to overlay two plots, use 'hold on' and 'hold off'
% %########################################################################

% %########################################################################
% b) now add to the plot above the non-indexed pixels in red.
% %########################################################################



%% 2) inspect orientations -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*

% %########################################################################
% a) access the first 100 Euler angles
% Hint: orientations are stored in ebsd(phase).orientations
% %########################################################################



% %########################################################################
% b) plot a polefigure (point plot)
%    - define crystal directions (@Miller)
%    - either use plotPDF or remember that orientiation*Miller -> specimen direction
% %########################################################################




% %########################################################################
% c) plot an inverse pole figure (point plot) with respect to (0,0,1), (-5,1,0)
%     and (2,-3,0)
%    - define a specimen direction (@vector3d)
%    - either use plotIPDF or remember that orientation*vector3d -> crystal direction
% %########################################################################


%%
close all
%% 3) orientation mapping  -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
%  Colorcode the orientations in the map

% %########################################################################
% a) plot three maps, colorcoded by the three Euler angles.  Do you
% understand what you see?
% %########################################################################



% %########################################################################
% b) make an inverse pole figure colorcoding with respect to one of the
% directions in 2c) and plot it as a map 
% useful commands: ipfHSVKey the property .inversePoleFigureDirection
% %########################################################################



% %########################################################################
% c) colorcode points in a polefigure and in an inverse polefigure indentical
%    to the colorcoding in the map from b)
% %########################################################################



% %########################################################################
% d) plot a color mapping for the c-axis directions
% useful commands: HSVDirectionKey
% note: the c-axis is just a direction, use .direction2color
% %########################################################################



% %########################################################################
% e) plot a Euler angle mapping, an axis-angle mapping and compare with the
%    previous ones and compare. Which is best?
% %########################################################################



%% 4) Selecting data  -----------------------------------------------------
% %########################################################################
% a) Question:
% Do pixels indexed with just few bands (e.g. <=6) show a different 
% polefigure compared to those with more bands (e.g. >9) ? Where are they
% locatd in the map?
% %########################################################################



%%
 close all
 
%% manually select data - ebsd area
% %########################################################################
% b) Select a triangular area in the upper left of the map.
% useful commands: EBSD/inpolygon
% %########################################################################



% %########################################################################
% c) Plot the c-axis inclination along a N-S profile in the ebsd map.
% useful commands: EBSD/spatialProfile
% %########################################################################



%% -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
% 5) More questions:
% -------------------------------------------------------------------------
% a) The map is taken with respect to a kinematic referce system:
%  - the flow/shear plane is E-W (1,0,0)
%  - the instantaneous shortening direction is (2,-3,0)
%  - the plane of maximum shear stress is (-5,1,0)
%  - the long axis of the finite strain ellipsoid is at 7 degree with in
%    the map plane (rho =7*degree)




% %########################################################################
% a) Question: What can you say about the a-axis alignement in the sample
%  with respect to the kinemtic reference frame? 
%  Colorcode the a-axis direction as a function of the distance from the
%  shear plane, the plane of maximum shear stress and the X direction of
%  finite strain ellipsoid.
%
%  Hint: rotate the vector into crystal coordinates and compute the angle 
%  with the reference crystal direction. This way, you don't have to worry
%  about symmetrically different solutions.
% %########################################################################



% %########################################################################
% b) Plot a spot color mapping , color all orientations with <a> 
%  (i)roughly and (ii) exactly up to 25 degree away from the a-maximum  
% %########################################################################
% Hint: this will give you the direction of the maximum in the <a> axes
% polefigure, more on that on Friday
cs = ebsd('q').CS
a= Miller(1,1,-2,0,cs)
aF  =calcDensity(ebsd('q').orientations*a,'halfwidth',10*degree)
[~,maxpos] = max(aF);
% %########################################################################

%
% option 1

% %########################################################################
% c) Compare the <1,1,-2,0> distributions of those orientation which contain 
% c-axes at the periphery of the pole figure to those with c-axes in the center
% Maybe we need to wait till Friday to answer this question.
% %########################################################################



%%
 close all
 

% 