A complete brief intro to EBSD processing using SAnTex

This jupyter notebook will demonstrate SAnTex library and the modules available within it. This notebook starts from loading an EBSD data, to getting phases information, to getting the data, and plotting it, to calculation of ODF, PDF and IPF, rotating the EBSD data to match convention, to cleaning EBSD data, and then to the calculation of seismic anisotropy by loading material parameters

Load the required modules

Firstly, we have to load the required modules from the SAnTex package. This ensures that the methods which are written in the package are avaiable throughout the notebook

import numpy as np
from santex.ebsd import EBSD
from santex.tensor import Tensor
from santex.anisotropy import Anisotropy
from santex.material import Material

In this code snippet, we are importing necessary modules for various functionalities. numpy is imported as np, a widely used library for numerical operations in Python. The santex module is imported four times, each time importing a different sub-module: EBSD, Tensor, Anisotropy, and Material.

EBSD contains codes for analyzing EBSD (Electron Backscatter Diffraction) files, which are used to study crystallographic properties of materials. Tensor holds codes for tensor operations, which deals with the elastic stiffness tensors of the materials. Anisotropy module is responsible for calculating seismic anisotropy, a property related to the variation of seismic wave speed in different directions within a material. Material refers to a catalogue of minerals or materials present in a library, possibly used for referencing or analyzing material properties. By importing these modules, we gain access to pre-written functions and tools that facilitate the respective analyses and operations mentioned above.

Load the ebsd ctf file

Let’s load an EBSD file fdmrn02x01.ctf (Chatzaras et al., 2016). For loading the EBSD file, we need to use EBSD module to load an EBSD file written in ctf format. For this notebook, we are using “ebsd.ctf” which is stored in the same directory as this notebook. But you may place the file in another location, for e.g. data/ebsd.ctf and for loading that ebsd file, we need to use ebsd = EBSD(“data/ebsd.ctf”)

ebsd = EBSD("fdmrn02x01.ctf")

Get the available phases information

Let’s check what are the phases information present in the ebsd file and the phase abundences

x = ebsd.phases()
phases_data = ebsd.phases_data

phases_data

[(0, 'NaN', 10.891769157994323),
 (1, 'Forsterite', 54.43386944181646),
 (2, 'Chromite', 1.7733207190160831),
 (3, 'Diopside   CaMgSi2O6', 19.77256385998108),
 (4, 'Enstatite  Opx AV77', 13.128476821192054)]

Load the ebsd file in a dataframe

A brief introduction to what an EBSD file is: An EBSD file is a data file used in materials science to store information about the microscopic structure of materials, particularly focusing on crystallographic properties. Within this file, Euler Bunge Angles (Phi1, Phi, Phi2) describe the orientation of crystals within the material, akin to longitude, latitude, and rotation for Earth but at a microscopic scale. The Mean Angular Deviation (MAD) measures how much individual crystal orientations deviate from the average, providing insights into the material’s structural consistency. Coordinates X and Y pinpoint specific positions within the material, while ‘Phase’ distinguishes different regions with unique properties, such as varying crystal types. Additionally, terms like ‘Bands,’ ‘Error,’ ‘BC (Band Contrast),’ and ‘BS (Band Slope)’ detail aspects like diffraction patterns, measurement quality, and intensity characteristics, all crucial for understanding the material’s crystalline structure and behavior.

All these information are stored in different columns with every row representing a crystal within an EBSD map.

Let’s load the ebsd data in a dataframe using the following directive df = ebsd.get_ebsd_data() and check the information

df = ebsd.get_ebsd_data()
df

	Phase	X	Y	Bands	Error	Euler1	Euler2	Euler3	MAD	BC	BS
0	0	0.0000	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
1	0	30.0094	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
2	0	60.0188	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
3	0	90.0282	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
4	0	120.0375	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
...	...	...	...	...	...	...	...	...	...	...	...
528495	0	31569.8744	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528496	0	31599.8838	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528497	0	31629.8931	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528498	0	31659.9025	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528499	0	31689.9119	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0

528500 rows × 11 columns

Let’s now remove the unwanted phases at index = [0, 2]

df = ebsd.filter_by_phase_number_list(df = df, phase_list = [0, 2])
df

	Phase	X	Y	Bands	Error	Euler1	Euler2	Euler3	MAD	BC	BS
245	4	7352.2996	0.0000	9	0	115.1774	85.8245	46.3400	0.6917	82	37
246	4	7382.3090	0.0000	8	0	115.0523	85.6425	46.1002	0.7481	71	36
247	4	7412.3184	0.0000	8	0	115.4108	85.5154	46.4322	0.4028	70	37
248	4	7442.3278	0.0000	10	0	115.5122	85.7559	46.5266	0.3767	68	34
249	4	7472.3372	0.0000	7	0	115.6362	85.7542	46.7128	0.2837	81	36
...	...	...	...	...	...	...	...	...	...	...	...
528031	1	17645.5191	14974.6838	10	0	81.8986	122.0925	107.8227	0.2657	141	56
528032	1	17675.5285	14974.6838	11	0	81.8970	122.6484	107.4921	0.4498	143	56
528033	1	17705.5379	14974.6838	11	0	81.7746	122.4250	107.4315	0.2815	142	56
528034	1	17735.5473	14974.6838	11	0	81.7065	122.2363	107.5466	0.3332	143	57
528035	1	17765.5567	14974.6838	11	0	81.9176	122.1538	107.6018	0.3617	142	55

461565 rows × 11 columns

Here you can see that the dataframe suggests at row 0, with a coordinate of X = 0, y = 0, or at X = 50, Y = 0, we have Phase = 1, and looking at the previous cell, we can say that the Phase 1 corresponds to Forsterite. At row 298068, we can check that it suggests Phase = 0, which suggests that its “NaN” from phases data which means thats a non indexed phase, which happens because the ebsd instrument isn’t able to give a measurement at that location.

Let’s look at the ebsd map and interpret further

But before that, we need to rotate the ebsd file to match the SEM orientation

rotated = ebsd.rotate_ebsd(ebsd_df = df, angles=[180, 0, 0], keepXY=True)
filtered_df = ebsd.filter_MAD(rotated, 0.7)

100%|█████████████████████████████████████████████████████████| 461565/461565 [00:41<00:00, 11012.88it/s]
/Users/usin8611/anaconda3/envs/santex/lib/python3.9/site-packages/santex/ebsd/rotateEBSD.py:81: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame.
Try using .loc[row_indexer,col_indexer] = value instead

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy
  euler_df[['Euler1', 'Euler2', 'Euler3']] = results

ebsd.plot(cmap = "viridis")

In the above map, we can see the coordinates of x and y and the corresponding pixels with each representing a pixels, we can be able to look which phases each pixels correspond to. Let’s save the above plot with a dpi of 300, with legend plotted lower right for seeing upper right clearly, and color the phases as Forsterite as Green, Diopside as yellow and enstatite as red

phase_colors = ['green', 'yellow', 'red']
ebsd.plot(data = filtered_df, phase_colors=phase_colors, save_image=True, image_filename="fdmrn02x01.svg", dpi = 300)

Pole Figures, Pole Density function, Inverse pole figures

Now, the above plot doesnt account for the crystal symmetry. Let’s look at the pole figures and pole density function

Pole Figures: Pole figures are a way to visualize the orientation of crystallographic axes relative to a fixed sample reference frame. They help in understanding the preferred orientation of grains in a polycrystalline material by projecting crystallographic directions onto a stereographic plane.

Pole Density Function: This is like zooming in on a specific part of the pile of books and counting how many books are stacked in a particular orientation. It gives a more detailed view of how the books are aligned in a specific area, which can be useful for studying local patterns or properties.

In simple words ODF and PDF consider crystal symmetries and a fixed sample reference

WARNING:
    MTEX and SAnTex (via Orix) do not necessarily use the same point group symmetry
    operation conventions. MTEX has different default crystal axes alignments for
    different crystal systems, while Orix always assumes X || a, Z || c*.
    Users should verify conventions and expectations before trusting the results.

df = filtered_df
ebsd.pf(df)

ebsd.pdf(df)

Lets look at inverse pole figure. Inverse pole figure fixes the crystal symmetry and looks at orientation by changing the specimen reference

ebsd.ipf(df)

ebsd.ipf_colorkey(df, phase = 1, crystal_symmetry='D2')

Plot the ebsd file, default: x2east, zoutofplane

Lets plot the ebsd map relative to the instrument reference system. We can check the map corresponding to the odf and pdf too.

The default plotting system ebsd.plot(df) plots x to the east and z axis out of the plane

ebsd.plot(df)

## Placeholder for defining densities and stiffness tensors the traditional way


# rho_Fo = 3355
# rho_diop = 3286
# rho_ens = 3306

# cij_Fo = np.array([[320.5, 68.15, 71.6, 0, 0, 0],
#                 [68.15, 196.5, 76.8, 0, 0, 0],
#                 [71.6, 76.8, 233.5, 0, 0, 0],
#                 [0, 0, 0, 64, 0, 0],
#                 [0, 0, 0, 0, 77, 0],
#                 [0, 0, 0, 0, 0, 78.7]])
# cij_ens = np.array([[236.90, 79.60, 63.20, 0.00, 0.00, 0.00],
#                 [79.60, 180.50, 56.80, 0.00, 0.00, 0.00],
#                 [63.20, 56.80, 230.40, 0.00, 0.00, 0.00],
#                 [0.00, 0.00, 0.00, 84.30, 0.00, 0.00],
#                 [0.00, 0.00, 0.00, 0.00, 79.40, 0.00],
#                 [0.00, 0.00, 0.00, 0.00, 0.00, 80.10]])
# cij_diop = np.array([[228.10, 78.80, 70.20, 0.00, 7.90, 0.00],
#                 [78.80, 181.10, 61.10, 0.00, 5.90, 0.00],
#                 [70.20, 61.10, 245.40, 0.00, 39.70, 0.00],
#                 [0.00, 0.00, 0.00, 78.90, 0.00, 6.40],
#                 [7.90, 5.90, 39.70, 0.00, 68.20, 0.00],
#                 [0.00, 0.00, 0.00, 6.40, 0.00, 78.10]])
# cij = [cij_Fo, cij_ens, cij_diop]
# density = [rho_Fo, rho_ens, rho_diop]

Calculate the stiffness tensors and the density at Pressure = 1.4GPa and Temperature = 1100 K (Chatzaras et al. 2016)

Now, lets import the material stiffness tensors and density at different pressure and temperature. We do this by creating a material instance of the Material class, and then calling the load_density function and voigthighPT function of the Material class in the following way. The resultant is stored in a variable, which in this case is rho_Fo and cij_Fo, and consequently we create an array to keep these matrix and density in a different variable, here it is called cij, density. The voigt stiffness matrix and the density of material is called at pressure = 1.4GPa and Temperature = 1100K

material_instance = Material()
rho_Fo = material_instance.load_density("Forsterite", 1.4, 1100)
rho_diop = material_instance.load_density("Diopside", 1.4, 1100)
rho_ens = material_instance.load_density("Enstatite", 1.4, 1100)
cij_Fo = material_instance.voigt_high_PT('Forsterite', PRESSURE = 1.4, TEMP = 1100)
cij_ens = material_instance.voigt_high_PT('Enstatite', PRESSURE = 1.4, TEMP = 1100)
cij_diop = material_instance.voigt_high_PT('Diopside', PRESSURE = 1.4, TEMP = 1100)

cij = [cij_Fo, cij_ens, cij_diop]
density = [rho_Fo, rho_ens, rho_diop]

Get the euler angles from ebsd file

Now, we also obtain the euler angles corresponding to the phases whose matrix and density information we extracted above, and then store them in an array format in the variable called euler_angles

forsterite = ebsd.get_euler_angles(phase = 1, data=df)
enstatite = ebsd.get_euler_angles(phase = 4, data=df)
diopside = ebsd.get_euler_angles(phase = 3, data=df)
euler_angles = [forsterite, enstatite, diopside]

Lets calculate anisotropy from the ebsd file. we first calculate the average tensor and average density by giving the directive average_tensor, average_density = ebsd.getAnisotropyForEBSD(cij, euler_angles, density), and then calculate anisotropy after creating an instance of the Anisotropy class, which can be called using anis = Anisotropy(average_tensor*10**9, average_density). Afterwards we plot the anisotropy using the directive anis.plot()

average_tensor, average_density = ebsd.get_anisotropy_for_ebsd(cij, euler_angles, density)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot()

Now, lets calculate some values related to anisotropy measures

anis.anisotropy_values()

Max Vp:  8823.853748923915
Min Vp:  7998.386279432881
Max Vs1:  4934.374218670314
Min Vs1:  4791.7034984600705
Max Vs2:  4814.058835429751
Min Vs2:  4571.539467545672
Max vs anisotropy percent:  6.682765226570884
Min vs anisotropy percent:  0.026571727723363536
P wave anisotropy percent:  9.81400179880404
S1 Wave anisotropy percent:  2.9337770961661094
S2 Wave anisotropy percent:  5.1679042732352025
Velocity difference, maxdvs:  317.89041227343387
Vp/Vs1 ratio:  9.045070493104753
Mean vp:  8279.614460029607
Mean Vs1:  4865.249680161301
Mean Vs2:  4638.154680841542

{'maxvp': 8823.853748923915,
 'minvp': 7998.386279432881,
 'maxvs1': 4934.374218670314,
 'minvs1': 4791.7034984600705,
 'maxvs2': 4814.058835429751,
 'minvs2': 4571.539467545672,
 'max_vs_anisotropy_percent': 6.682765226570884,
 'min_vs_anisotropy_percent': 0.026571727723363536,
 'p_wave_anisotropy_percent': 9.81400179880404,
 's1_wave_anisotropy_percent': 2.9337770961661094,
 's2_wave_anisotropy_percent': 5.1679042732352025,
 'maxdvs': 317.89041227343387,
 'AVpVs1': 9.045070493104753}

Lets look at the average tensors of all the phases

cij

[array([[305.878,  66.318,  71.816,   0.   ,   0.   ,   0.   ],
        [ 66.318, 183.204,  78.874,   0.   ,   0.   ,   0.   ],
        [ 71.816,  78.874, 219.538,   0.   ,   0.   ,   0.   ],
        [  0.   ,   0.   ,   0.   ,  56.728,   0.   ,   0.   ],
        [  0.   ,   0.   ,   0.   ,   0.   ,  68.924,   0.   ],
        [  0.   ,   0.   ,   0.   ,   0.   ,   0.   ,  69.938]]),
 array([[224.196,  72.398,  50.486,   0.   ,   0.   ,   0.   ],
        [ 72.398, 167.126,  60.462,   0.   ,   0.   ,   0.   ],
        [ 50.486,  60.462, 212.108,   0.   ,   0.   ,   0.   ],
        [  0.   ,   0.   ,   0.   ,  77.152,   0.   ,   0.   ],
        [  0.   ,   0.   ,   0.   ,   0.   ,  72.448,   0.   ],
        [  0.   ,   0.   ,   0.   ,   0.   ,   0.   ,  72.35 ]]),
 array([[225.754,  80.726,  77.638,   0.   ,   8.938,   0.   ],
        [ 80.726, 174.992,  57.664,   0.   ,  10.006,   0.   ],
        [ 77.638,  57.664, 218.672,   0.   ,  51.88 ,   0.   ],
        [  0.   ,   0.   ,   0.   ,  71.988,   0.   ,  11.522],
        [  8.938,  10.006,  51.88 ,   0.   ,  67.364,   0.   ],
        [  0.   ,   0.   ,   0.   ,  11.522,   0.   ,  75.504]])]

Calculate the stiffness tensors and the density at ambient pressure and temperature

Now, lets import the material stiffness tensors and density at ambient pressure and temperature, and calculate seismic anisotropy

material_instance = Material()
rho_Fo = material_instance.load_density("Forsterite")
rho_diop = material_instance.load_density("Diopside")
rho_ens = material_instance.load_density("Enstatite")
cij_Fo = material_instance.voigt_high_PT('Forsterite')
cij_ens = material_instance.voigt_high_PT('Enstatite')
cij_diop = material_instance.voigt_high_PT('Diopside')

cij = [cij_Fo, cij_ens, cij_diop]
density = [rho_Fo, rho_ens, rho_diop]

forsterite = ebsd.get_euler_angles(phase = 1, data=df)
enstatite = ebsd.get_euler_angles(phase = 4, data=df)
diopside = ebsd.get_euler_angles(phase = 3, data=df)
euler_angles = [forsterite, enstatite, diopside]

average_tensor, average_density = ebsd.get_anisotropy_for_ebsd(cij, euler_angles, density)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot(savefig = True, figname= "fdmrn02x01_amb_bulk", save_format = 'svg')

Calculate the stiffness tensors and the density at pressure = 1.4GPa and temperature = 1100K, and assumed melt of 7%

Now, lets import the material stiffness tensors and density at ambient pressure and temperature, and calculate seismic anisotropy

cij = [cij_Fo, cij_ens, cij_diop]
density = [rho_Fo, rho_ens, rho_diop]

forsterite = ebsd.get_euler_angles(phase = 1, data=df)
enstatite = ebsd.get_euler_angles(phase = 4, data=df)
diopside = ebsd.get_euler_angles(phase = 3, data=df)
euler_angles = [forsterite, enstatite, diopside]

average_tensor, average_density = ebsd.get_anisotropy_for_ebsd(cij, euler_angles, density, melt = 7)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot(savefig = True, figname= "fdmrn02x01_ptmelt_bulk", save_format = 'svg')

Rotating the dataset by 90 degree around the z-axis

Lets rotate the ebsd by 90 degree around the z-axis, and not keeping x and y fixed by evoking the rotateEBSD nmethod. That means that the X and Y will change after we make the rotation. There are 3 angles through which we can make rotations. The following directive stores the rotated ebsd file in another dataframe called rotated here

rotated = ebsd.rotate_ebsd(df, [180, 0, 90], keepXY = False)

/Users/usin8611/anaconda3/envs/santex/lib/python3.9/site-packages/joblib/externals/loky/process_executor.py:752: UserWarning: A worker stopped while some jobs were given to the executor. This can be caused by a too short worker timeout or by a memory leak.
  warnings.warn(
100%|█████████████████████████████████████████████████████████| 444045/444045 [00:36<00:00, 12002.13it/s]
/Users/usin8611/anaconda3/envs/santex/lib/python3.9/site-packages/santex/ebsd/rotateEBSD.py:129: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame.
Try using .loc[row_indexer,col_indexer] = value instead

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy
  df[['Euler1', 'Euler2', 'Euler3', 'X', 'Y']] = results

Lets look at the resultant dataframe rotated and compare it with the previous dataframe df

In previous datafame, X and Y started from 0, 0 and then X and Y incremented.

This is a sample of previous unrotated dataframe

Phase   X   Y   Bands   Error   Euler1  Euler2  Euler3  MAD BC  BS 0   0   0.0 0.0 0   5   360.0   180.0   0.0 0.0 0   0 1   0   50.0    0.0 0   5   360.0   180.0   0.0 0.0 0   0 2   0   100.0   0.0 0   5   360.0   180.0   0.0 0.0 0   0 3   0   150.0   0.0 0   5   360.0   180.0   0.0 0.0 0   0 4   0   200.0   0.0 0   5   360.0   180.0   0.0 0.0 0   0

rotated

	Phase	X	Y	Bands	Error	Euler1	Euler2	Euler3	MAD	BC	BS
245	4	-1.632538e-12	-7352.2996	9	0	205.1774	85.8245	46.3400	0.6917	82	37
247	4	-1.645865e-12	-7412.3184	8	0	205.4108	85.5154	46.4322	0.4028	70	37
248	4	-1.652529e-12	-7442.3278	10	0	205.5122	85.7559	46.5266	0.3767	68	34
249	4	-1.659192e-12	-7472.3372	7	0	205.6362	85.7542	46.7128	0.2837	81	36
250	4	-1.665856e-12	-7502.3466	10	0	205.6750	85.6075	46.6743	0.4016	73	37
...	...	...	...	...	...	...	...	...	...	...	...
528031	1	1.497468e+04	-17645.5191	10	0	171.8986	122.0925	107.8227	0.2657	141	56
528032	1	1.497468e+04	-17675.5285	11	0	171.8970	122.6484	107.4921	0.4498	143	56
528033	1	1.497468e+04	-17705.5379	11	0	171.7746	122.4250	107.4315	0.2815	142	56
528034	1	1.497468e+04	-17735.5473	11	0	171.7065	122.2363	107.5466	0.3332	143	57
528035	1	1.497468e+04	-17765.5567	11	0	171.9176	122.1538	107.6018	0.3617	142	55

444045 rows × 11 columns

Lets look at the rotated dataframe visually. Now here in the bottom plot we can see that with the same plot command, the x and y values now correspond to different pixels

ebsd.plot(rotated)

Let’s now again calculate and plot the anisotropy maps using the procedure explained before

forsterite = ebsd.get_euler_angles(phase = 1, data=rotated)
enstatite = ebsd.get_euler_angles(phase = 2, data=rotated)
diopside = ebsd.get_euler_angles(phase = 3, data=rotated)
euler_angles = [forsterite, enstatite, diopside]

average_tensor, average_density = ebsd.get_anisotropy_for_ebsd(cij, euler_angles, density)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot()

anis.anisotropy_values()

Max Vp:  8887.111183165302
Min Vp:  8024.4623724207395
Max Vs1:  4991.568023273352
Min Vs1:  4855.654517513256
Max Vs2:  4882.735691095535
Min Vs2:  4642.7902826485015
Max vs anisotropy percent:  6.348485085847311
Min vs anisotropy percent:  0.027701355657550178
P wave anisotropy percent:  10.201875158561121
S1 Wave anisotropy percent:  2.760443469153875
S2 Wave anisotropy percent:  5.037945602340787
Velocity difference, maxdvs:  305.8997558255105
Vp/Vs1 ratio:  9.262001273802962
Mean vp:  8313.309412398487
Mean Vs1:  4929.986539946562
Mean Vs2:  4709.178011889468

{'maxvp': 8887.111183165302,
 'minvp': 8024.4623724207395,
 'maxvs1': 4991.568023273352,
 'minvs1': 4855.654517513256,
 'maxvs2': 4882.735691095535,
 'minvs2': 4642.7902826485015,
 'max_vs_anisotropy_percent': 6.348485085847311,
 'min_vs_anisotropy_percent': 0.027701355657550178,
 'p_wave_anisotropy_percent': 10.201875158561121,
 's1_wave_anisotropy_percent': 2.760443469153875,
 's2_wave_anisotropy_percent': 5.037945602340787,
 'maxdvs': 305.8997558255105,
 'AVpVs1': 9.262001273802962}

Rotating ebsd for matching SEM orientation

ebsd = EBSD("ebsd.ctf")
df = ebsd.get_ebsd_data()
df = ebsd.filter_by_phase_number_list(df = df, phase_list = [0, 2])

Let’s now rotate ebsd which plots x to south and z into plane, without keeping keeping x and y same

rotated_df = ebsd.plot_rotate_ebsd(sample_ref = ["x2west", "zIntoPlane"], ebsd_df = df, keepXY = False)

100%|█████████████████████████████████████████████████████████| 184447/184447 [00:13<00:00, 14186.45it/s]

Lets check how the rotated dataframe looks like. We have rotated the dataframe without keeping X and Y intact, essentially we have changed the plotting conventions

rotated_df

	Phase	X	Y	Bands	Error	Euler1	Euler2	Euler3	MAD	BC	BS
690	1	-34500.0	4.225031e-12	11	0	327.8425	144.2913	153.1576	0.7369	112	74
691	1	-34550.0	4.231155e-12	11	0	327.8451	144.0361	153.6517	0.6896	134	80
692	1	-34600.0	4.237278e-12	11	0	328.1161	144.1026	153.8256	0.6810	123	77
693	1	-34650.0	4.243401e-12	11	0	328.0592	143.9367	153.7876	0.8815	122	80
694	1	-34700.0	4.249524e-12	11	0	328.6093	143.8863	154.4789	0.8696	127	80
...	...	...	...	...	...	...	...	...	...	...	...
297374	1	-700.0	2.100000e+04	10	0	261.9566	64.0177	125.3487	0.7277	107	63
297375	1	-750.0	2.100000e+04	12	0	262.2138	64.1060	125.0667	0.6766	150	78
297376	1	-800.0	2.100000e+04	12	0	262.2349	64.0659	125.1368	0.6592	148	85
297377	1	-850.0	2.100000e+04	12	0	262.3693	64.1572	125.0572	0.6526	158	83
297378	1	-900.0	2.100000e+04	12	0	262.2358	64.2606	125.0286	0.6438	156	83

184447 rows × 11 columns

Lets plot the new ebsd map with different plotting conventions

ebsd.plot(rotated_df)

Lets check the anisotropy now, and compare it with previous one

forsterite = ebsd.get_euler_angles(phase = 1, data=rotated_df)
enstatite = ebsd.get_euler_angles(phase = 2, data=rotated_df)
diopside = ebsd.get_euler_angles(phase = 3, data=rotated_df)
euler_angles = [forsterite, enstatite, diopside]

average_tensor, average_density = ebsd.get_anisotropy_for_ebsd(cij, euler_angles, density)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot()

anis.anisotropy_values()

Max Vp:  8657.982082576274
Min Vp:  7915.162814563828
Max Vs1:  5018.475703191029
Min Vs1:  4752.420087233196
Max Vs2:  4813.023974764757
Min Vs2:  4716.876384343068
Max vs anisotropy percent:  5.828421540277186
Min vs anisotropy percent:  0.0116412173187977
P wave anisotropy percent:  8.96413170370132
S1 Wave anisotropy percent:  5.445879715933032
S2 Wave anisotropy percent:  2.017808933958061
Velocity difference, maxdvs:  284.0918760746672
Vp/Vs1 ratio:  4.674999301124261
Mean vp:  8424.382342662422
Mean Vs1:  4935.801692347352
Mean Vs2:  4759.749654899821

{'maxvp': 8657.982082576274,
 'minvp': 7915.162814563828,
 'maxvs1': 5018.475703191029,
 'minvs1': 4752.420087233196,
 'maxvs2': 4813.023974764757,
 'minvs2': 4716.876384343068,
 'max_vs_anisotropy_percent': 5.828421540277186,
 'min_vs_anisotropy_percent': 0.0116412173187977,
 'p_wave_anisotropy_percent': 8.96413170370132,
 's1_wave_anisotropy_percent': 5.445879715933032,
 's2_wave_anisotropy_percent': 2.017808933958061,
 'maxdvs': 284.0918760746672,
 'AVpVs1': 4.674999301124261}

Lets now rotate the ebsd file with some rotation angles

ebsd = EBSD("fdmrn02x01.ctf")
df = ebsd.get_ebsd_data()
df = ebsd.filter_by_phase_number_list(df = df, phase_list = [4, 5, 6, 7])

ebsd.plot(data = df, rotation_angle=0)

ebsd.plot(rotation_angle=90)

ebsd.plot(rotation_angle=180)

ebsd.plot(rotation_angle=180, mirror=True)

ebsdfile = EBSD("fdmrn02x01.ctf")

df = ebsdfile.get_ebsd_data()

rotated = ebsdfile.rotate_ebsd(df, [70, 160, 50], keepXY = False)

100%|█████████████████████████████████████████████████████████| 528500/528500 [00:30<00:00, 17173.08it/s]

ebsdfile.plot(rotated)

Now, for the above rotation of EBSD data for custom angles, we can calculate anisotropy in the same way described above

rotated_df = rotated

rho_Fo = 3355
rho_diop = 3286
rho_ens = 3306

cij_Fo = np.array([[320.5, 68.15, 71.6, 0, 0, 0],
                [68.15, 196.5, 76.8, 0, 0, 0],
                [71.6, 76.8, 233.5, 0, 0, 0],
                [0, 0, 0, 64, 0, 0],
                [0, 0, 0, 0, 77, 0],
                [0, 0, 0, 0, 0, 78.7]])
cij_ens = np.array([[236.90, 79.60, 63.20, 0.00, 0.00, 0.00],
                [79.60, 180.50, 56.80, 0.00, 0.00, 0.00],
                [63.20, 56.80, 230.40, 0.00, 0.00, 0.00],
                [0.00, 0.00, 0.00, 84.30, 0.00, 0.00],
                [0.00, 0.00, 0.00, 0.00, 79.40, 0.00],
                [0.00, 0.00, 0.00, 0.00, 0.00, 80.10]])
cij_diop = np.array([[228.10, 78.80, 70.20, 0.00, 7.90, 0.00],
                [78.80, 181.10, 61.10, 0.00, 5.90, 0.00],
                [70.20, 61.10, 245.40, 0.00, 39.70, 0.00],
                [0.00, 0.00, 0.00, 78.90, 0.00, 6.40],
                [7.90, 5.90, 39.70, 0.00, 68.20, 0.00],
                [0.00, 0.00, 0.00, 6.40, 0.00, 78.10]])
cij = [cij_Fo, cij_ens, cij_diop]
density = [rho_Fo, rho_ens, rho_diop]

forsterite = ebsdfile.get_euler_angles(phase = 1, data=rotated_df)
enstatite = ebsdfile.get_euler_angles(phase = 2, data=rotated_df)
diopside = ebsdfile.get_euler_angles(phase = 3, data=rotated_df)
euler_angles = [forsterite, enstatite, diopside]

average_tensor, average_density = ebsdfile.get_anisotropy_for_ebsd(cij, euler_angles, density)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot()

Cleaning EBSD data

The steps of cleaning EBSD data includes removing unwanted phases and removing pixels whose mean angular deviation is more than some specified threshold, lets say 0.7 and to remove grains with mean misorientation more than some specified threshold, lets say 10 degrees. In the above steps, we aim to clean EBSD data on these principles.

ebsdfile = EBSD("fdmrn02x01.ctf")

phases_available = ebsdfile.phases()
print(phases_available)

[(0, 'NaN', 10.891769157994323), (1, 'Forsterite', 54.43386944181646), (2, 'Chromite', 1.7733207190160831), (3, 'Diopside   CaMgSi2O6', 19.77256385998108), (4, 'Enstatite  Opx AV77', 13.128476821192054)]

df = ebsdfile.get_ebsd_data()

df

	Phase	X	Y	Bands	Error	Euler1	Euler2	Euler3	MAD	BC	BS
0	0	0.0000	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
1	0	30.0094	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
2	0	60.0188	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
3	0	90.0282	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
4	0	120.0375	0.0000	0	5	0.0	0.0	0.0	0.0	0	0
...	...	...	...	...	...	...	...	...	...	...	...
528495	0	31569.8744	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528496	0	31599.8838	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528497	0	31629.8931	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528498	0	31659.9025	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0
528499	0	31689.9119	14974.6838	0	5	0.0	0.0	0.0	0.0	0	0

528500 rows × 11 columns

ebsdfile.plot(df)

Removing data with MAD > 0.7

We now aim to remove pixels with mean angular deviation greater than 0.7, which can be done by invoking the filterMAD() method of the ebsd object

filtered_df = ebsd.filter_MAD(df, 0.7)
print(filtered_df)

        Phase           X           Y  Bands  Error  Euler1  Euler2  Euler3  \
0           0      0.0000      0.0000      0      5     0.0     0.0     0.0
1           0     30.0094      0.0000      0      5     0.0     0.0     0.0
2           0     60.0188      0.0000      0      5     0.0     0.0     0.0
3           0     90.0282      0.0000      0      5     0.0     0.0     0.0
4           0    120.0375      0.0000      0      5     0.0     0.0     0.0
...       ...         ...         ...    ...    ...     ...     ...     ...
528495      0  31569.8744  14974.6838      0      5     0.0     0.0     0.0
528496      0  31599.8838  14974.6838      0      5     0.0     0.0     0.0
528497      0  31629.8931  14974.6838      0      5     0.0     0.0     0.0
528498      0  31659.9025  14974.6838      0      5     0.0     0.0     0.0
528499      0  31689.9119  14974.6838      0      5     0.0     0.0     0.0

        MAD  BC  BS
0       0.0   0   0
1       0.0   0   0
2       0.0   0   0
3       0.0   0   0
4       0.0   0   0
...     ...  ..  ..
528495  0.0   0   0
528496  0.0   0   0
528497  0.0   0   0
528498  0.0   0   0
528499  0.0   0   0

[510161 rows x 11 columns]

Here we again check the ebsd plot after removing pixels with MAD>0.7

ebsdfile.plot(filtered_df)

Calculating grain boundaries with threshold of 10 degrees misorientation

Let’s now calculate grain boundaries with misorientation of 10 degrees. we calculate that by invoking the calcGrains() method. Here, we can modify the angles to classify the grains by modifying the parameter threshold

phases_names = ebsdfile.phases_names()
phases_names = phases_names['phase'].tolist()
phases_names.insert(0, "na")

df_grain_boundary = ebsdfile.calc_grains(df = filtered_df, threshold = 10, phase_names=phases_names, downsampling_factor=10)

Grain 0: Dominant Phase - Enstatite  Opx AV77, Size - 694
Grain 1: Dominant Phase - Forsterite, Size - 72
Grain 2: Dominant Phase - Diopside   CaMgSi2O6, Size - 47
Grain 3: Dominant Phase - Forsterite, Size - 6
Grain 4: Dominant Phase - Diopside   CaMgSi2O6, Size - 65
Grain 5: Dominant Phase - Forsterite, Size - 334
Grain 6: Dominant Phase - Diopside   CaMgSi2O6, Size - 12
Grain 7: Dominant Phase - Forsterite, Size - 1298
Grain 8: Dominant Phase - Forsterite, Size - 408
Grain 9: Dominant Phase - Forsterite, Size - 31
Grain 10: Dominant Phase - Forsterite, Size - 156
Grain 11: Dominant Phase - Forsterite, Size - 11
Grain 12: Dominant Phase - Forsterite, Size - 202
Grain 13: Dominant Phase - Enstatite  Opx AV77, Size - 58
Grain 14: Dominant Phase - Forsterite, Size - 838
Grain 15: Dominant Phase - Enstatite  Opx AV77, Size - 28
Grain 16: Dominant Phase - Forsterite, Size - 1010
Grain 17: Dominant Phase - Forsterite, Size - 17
Grain 18: Dominant Phase - Diopside   CaMgSi2O6, Size - 20
Grain 19: Dominant Phase - Forsterite, Size - 61
Grain 20: Dominant Phase - Forsterite, Size - 431
Grain 21: Dominant Phase - Chromite, Size - 40
Grain 22: Dominant Phase - Diopside   CaMgSi2O6, Size - 36
Grain 23: Dominant Phase - Forsterite, Size - 842
Grain 24: Dominant Phase - Diopside   CaMgSi2O6, Size - 9
Grain 25: Dominant Phase - Forsterite, Size - 1
Grain 26: Dominant Phase - Enstatite  Opx AV77, Size - 389
Grain 27: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 28: Dominant Phase - Forsterite, Size - 485
Grain 29: Dominant Phase - Forsterite, Size - 196
Grain 30: Dominant Phase - Forsterite, Size - 162
Grain 31: Dominant Phase - Forsterite, Size - 7
Grain 32: Dominant Phase - Forsterite, Size - 519
Grain 33: Dominant Phase - Forsterite, Size - 603
Grain 34: Dominant Phase - Forsterite, Size - 536
Grain 35: Dominant Phase - Enstatite  Opx AV77, Size - 406
Grain 36: Dominant Phase - Enstatite  Opx AV77, Size - 253
Grain 37: Dominant Phase - Forsterite, Size - 46
Grain 38: Dominant Phase - Diopside   CaMgSi2O6, Size - 237
Grain 39: Dominant Phase - Forsterite, Size - 59
Grain 40: Dominant Phase - Forsterite, Size - 3
Grain 41: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 42: Dominant Phase - Enstatite  Opx AV77, Size - 200
Grain 43: Dominant Phase - Diopside   CaMgSi2O6, Size - 65
Grain 44: Dominant Phase - Forsterite, Size - 19
Grain 45: Dominant Phase - Forsterite, Size - 6
Grain 46: Dominant Phase - Forsterite, Size - 21
Grain 47: Dominant Phase - Forsterite, Size - 9
Grain 48: Dominant Phase - Forsterite, Size - 13
Grain 49: Dominant Phase - Enstatite  Opx AV77, Size - 239
Grain 50: Dominant Phase - Forsterite, Size - 47
Grain 51: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 52: Dominant Phase - Forsterite, Size - 172
Grain 53: Dominant Phase - Diopside   CaMgSi2O6, Size - 111
Grain 54: Dominant Phase - Forsterite, Size - 58
Grain 55: Dominant Phase - Forsterite, Size - 258
Grain 56: Dominant Phase - Forsterite, Size - 31
Grain 57: Dominant Phase - Forsterite, Size - 12
Grain 58: Dominant Phase - Forsterite, Size - 5
Grain 59: Dominant Phase - Chromite, Size - 13
Grain 60: Dominant Phase - Diopside   CaMgSi2O6, Size - 56
Grain 61: Dominant Phase - Forsterite, Size - 6
Grain 62: Dominant Phase - Diopside   CaMgSi2O6, Size - 52
Grain 63: Dominant Phase - Diopside   CaMgSi2O6, Size - 198
Grain 64: Dominant Phase - Diopside   CaMgSi2O6, Size - 14
Grain 65: Dominant Phase - Forsterite, Size - 268
Grain 66: Dominant Phase - Forsterite, Size - 12
Grain 67: Dominant Phase - Diopside   CaMgSi2O6, Size - 852
Grain 68: Dominant Phase - Chromite, Size - 1
Grain 69: Dominant Phase - Forsterite, Size - 378
Grain 70: Dominant Phase - Forsterite, Size - 590
Grain 71: Dominant Phase - Enstatite  Opx AV77, Size - 186
Grain 72: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 73: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 74: Dominant Phase - Forsterite, Size - 116
Grain 75: Dominant Phase - Forsterite, Size - 2
Grain 76: Dominant Phase - Forsterite, Size - 2
Grain 77: Dominant Phase - Forsterite, Size - 6
Grain 78: Dominant Phase - Forsterite, Size - 2
Grain 79: Dominant Phase - Forsterite, Size - 2
Grain 80: Dominant Phase - Forsterite, Size - 2
Grain 81: Dominant Phase - Forsterite, Size - 61
Grain 82: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 83: Dominant Phase - Forsterite, Size - 94
Grain 84: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 85: Dominant Phase - Forsterite, Size - 1
Grain 86: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 87: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 88: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 89: Dominant Phase - Forsterite, Size - 10
Grain 90: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 91: Dominant Phase - Diopside   CaMgSi2O6, Size - 240
Grain 92: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 93: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 94: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 95: Dominant Phase - Forsterite, Size - 121
Grain 96: Dominant Phase - Forsterite, Size - 8
Grain 97: Dominant Phase - Forsterite, Size - 10
Grain 98: Dominant Phase - Enstatite  Opx AV77, Size - 11
Grain 99: Dominant Phase - Forsterite, Size - 16
Grain 100: Dominant Phase - Forsterite, Size - 129
Grain 101: Dominant Phase - Diopside   CaMgSi2O6, Size - 84
Grain 102: Dominant Phase - Forsterite, Size - 11
Grain 103: Dominant Phase - Diopside   CaMgSi2O6, Size - 17
Grain 104: Dominant Phase - Diopside   CaMgSi2O6, Size - 29
Grain 105: Dominant Phase - Chromite, Size - 82
Grain 106: Dominant Phase - Forsterite, Size - 24
Grain 107: Dominant Phase - Forsterite, Size - 2
Grain 108: Dominant Phase - Diopside   CaMgSi2O6, Size - 46
Grain 109: Dominant Phase - Diopside   CaMgSi2O6, Size - 9
Grain 110: Dominant Phase - Forsterite, Size - 12
Grain 111: Dominant Phase - Forsterite, Size - 139
Grain 112: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 113: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 114: Dominant Phase - Forsterite, Size - 69
Grain 115: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 116: Dominant Phase - Diopside   CaMgSi2O6, Size - 24
Grain 117: Dominant Phase - Forsterite, Size - 1174
Grain 118: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 119: Dominant Phase - Forsterite, Size - 2
Grain 120: Dominant Phase - Forsterite, Size - 106
Grain 121: Dominant Phase - Enstatite  Opx AV77, Size - 291
Grain 122: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 123: Dominant Phase - Forsterite, Size - 9
Grain 124: Dominant Phase - Diopside   CaMgSi2O6, Size - 42
Grain 125: Dominant Phase - Enstatite  Opx AV77, Size - 40
Grain 126: Dominant Phase - Forsterite, Size - 34
Grain 127: Dominant Phase - Forsterite, Size - 350
Grain 128: Dominant Phase - Forsterite, Size - 15
Grain 129: Dominant Phase - Forsterite, Size - 4
Grain 130: Dominant Phase - Forsterite, Size - 45
Grain 131: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 132: Dominant Phase - Chromite, Size - 1
Grain 133: Dominant Phase - Forsterite, Size - 162
Grain 134: Dominant Phase - Forsterite, Size - 195
Grain 135: Dominant Phase - Enstatite  Opx AV77, Size - 157
Grain 136: Dominant Phase - Forsterite, Size - 11
Grain 137: Dominant Phase - Chromite, Size - 1
Grain 138: Dominant Phase - Diopside   CaMgSi2O6, Size - 185
Grain 139: Dominant Phase - Chromite, Size - 17
Grain 140: Dominant Phase - Diopside   CaMgSi2O6, Size - 54
Grain 141: Dominant Phase - Diopside   CaMgSi2O6, Size - 81
Grain 142: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 143: Dominant Phase - Forsterite, Size - 231
Grain 144: Dominant Phase - Enstatite  Opx AV77, Size - 199
Grain 145: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 146: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 147: Dominant Phase - Diopside   CaMgSi2O6, Size - 163
Grain 148: Dominant Phase - Forsterite, Size - 56
Grain 149: Dominant Phase - Forsterite, Size - 1050
Grain 150: Dominant Phase - Forsterite, Size - 46
Grain 151: Dominant Phase - Diopside   CaMgSi2O6, Size - 11
Grain 152: Dominant Phase - Chromite, Size - 1
Grain 153: Dominant Phase - Forsterite, Size - 125
Grain 154: Dominant Phase - Forsterite, Size - 35
Grain 155: Dominant Phase - Diopside   CaMgSi2O6, Size - 199
Grain 156: Dominant Phase - Forsterite, Size - 74
Grain 157: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 158: Dominant Phase - Chromite, Size - 9
Grain 159: Dominant Phase - Forsterite, Size - 31
Grain 160: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 161: Dominant Phase - Diopside   CaMgSi2O6, Size - 13
Grain 162: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 163: Dominant Phase - Diopside   CaMgSi2O6, Size - 79
Grain 164: Dominant Phase - Forsterite, Size - 340
Grain 165: Dominant Phase - Forsterite, Size - 84
Grain 166: Dominant Phase - Forsterite, Size - 237
Grain 167: Dominant Phase - Forsterite, Size - 326
Grain 168: Dominant Phase - Diopside   CaMgSi2O6, Size - 72
Grain 169: Dominant Phase - Enstatite  Opx AV77, Size - 7
Grain 170: Dominant Phase - Forsterite, Size - 1
Grain 171: Dominant Phase - Forsterite, Size - 14
Grain 172: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 173: Dominant Phase - Forsterite, Size - 328
Grain 174: Dominant Phase - Forsterite, Size - 98
Grain 175: Dominant Phase - Forsterite, Size - 30
Grain 176: Dominant Phase - Diopside   CaMgSi2O6, Size - 22
Grain 177: Dominant Phase - Diopside   CaMgSi2O6, Size - 87
Grain 178: Dominant Phase - Forsterite, Size - 242
Grain 179: Dominant Phase - Forsterite, Size - 278
Grain 180: Dominant Phase - Enstatite  Opx AV77, Size - 4
Grain 181: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 182: Dominant Phase - Forsterite, Size - 2
Grain 183: Dominant Phase - Forsterite, Size - 799
Grain 184: Dominant Phase - Diopside   CaMgSi2O6, Size - 64
Grain 185: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 186: Dominant Phase - Forsterite, Size - 42
Grain 187: Dominant Phase - Forsterite, Size - 41
Grain 188: Dominant Phase - Diopside   CaMgSi2O6, Size - 53
Grain 189: Dominant Phase - Forsterite, Size - 531
Grain 190: Dominant Phase - Diopside   CaMgSi2O6, Size - 63
Grain 191: Dominant Phase - Enstatite  Opx AV77, Size - 6
Grain 192: Dominant Phase - Forsterite, Size - 1
Grain 193: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 194: Dominant Phase - Forsterite, Size - 13
Grain 195: Dominant Phase - Forsterite, Size - 24
Grain 196: Dominant Phase - Forsterite, Size - 121
Grain 197: Dominant Phase - Chromite, Size - 4
Grain 198: Dominant Phase - Diopside   CaMgSi2O6, Size - 25
Grain 199: Dominant Phase - Diopside   CaMgSi2O6, Size - 82
Grain 200: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 201: Dominant Phase - Forsterite, Size - 21
Grain 202: Dominant Phase - Forsterite, Size - 80
Grain 203: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 204: Dominant Phase - Forsterite, Size - 4
Grain 205: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 206: Dominant Phase - Forsterite, Size - 70
Grain 207: Dominant Phase - Forsterite, Size - 230
Grain 208: Dominant Phase - Diopside   CaMgSi2O6, Size - 13
Grain 209: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 210: Dominant Phase - Forsterite, Size - 340
Grain 211: Dominant Phase - Forsterite, Size - 3
Grain 212: Dominant Phase - Diopside   CaMgSi2O6, Size - 37
Grain 213: Dominant Phase - Forsterite, Size - 218
Grain 214: Dominant Phase - Forsterite, Size - 8
Grain 215: Dominant Phase - Forsterite, Size - 28
Grain 216: Dominant Phase - Enstatite  Opx AV77, Size - 3
Grain 217: Dominant Phase - Forsterite, Size - 3
Grain 218: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 219: Dominant Phase - Diopside   CaMgSi2O6, Size - 8
Grain 220: Dominant Phase - Forsterite, Size - 1
Grain 221: Dominant Phase - Forsterite, Size - 2
Grain 222: Dominant Phase - Forsterite, Size - 6
Grain 223: Dominant Phase - Forsterite, Size - 1
Grain 224: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 225: Dominant Phase - Forsterite, Size - 7
Grain 226: Dominant Phase - Forsterite, Size - 4
Grain 227: Dominant Phase - Diopside   CaMgSi2O6, Size - 55
Grain 228: Dominant Phase - Forsterite, Size - 787
Grain 229: Dominant Phase - Forsterite, Size - 98
Grain 230: Dominant Phase - Forsterite, Size - 3
Grain 231: Dominant Phase - Forsterite, Size - 157
Grain 232: Dominant Phase - Forsterite, Size - 6
Grain 233: Dominant Phase - Chromite, Size - 5
Grain 234: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 235: Dominant Phase - Diopside   CaMgSi2O6, Size - 20
Grain 236: Dominant Phase - Forsterite, Size - 9
Grain 237: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 238: Dominant Phase - Forsterite, Size - 17
Grain 239: Dominant Phase - Forsterite, Size - 2
Grain 240: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 241: Dominant Phase - Forsterite, Size - 1
Grain 242: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 243: Dominant Phase - Forsterite, Size - 3
Grain 244: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 245: Dominant Phase - Enstatite  Opx AV77, Size - 131
Grain 246: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 247: Dominant Phase - Forsterite, Size - 126
Grain 248: Dominant Phase - Forsterite, Size - 39
Grain 249: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 250: Dominant Phase - Diopside   CaMgSi2O6, Size - 160
Grain 251: Dominant Phase - Enstatite  Opx AV77, Size - 629
Grain 252: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 253: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 254: Dominant Phase - Forsterite, Size - 1
Grain 255: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 256: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 257: Dominant Phase - Forsterite, Size - 10
Grain 258: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 259: Dominant Phase - Chromite, Size - 1
Grain 260: Dominant Phase - Diopside   CaMgSi2O6, Size - 10
Grain 261: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 262: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 263: Dominant Phase - Diopside   CaMgSi2O6, Size - 94
Grain 264: Dominant Phase - Forsterite, Size - 52
Grain 265: Dominant Phase - Forsterite, Size - 24
Grain 266: Dominant Phase - Forsterite, Size - 349
Grain 267: Dominant Phase - Forsterite, Size - 175
Grain 268: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 269: Dominant Phase - Chromite, Size - 1
Grain 270: Dominant Phase - Forsterite, Size - 2
Grain 271: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 272: Dominant Phase - Forsterite, Size - 54
Grain 273: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 274: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 275: Dominant Phase - Forsterite, Size - 306
Grain 276: Dominant Phase - Forsterite, Size - 1
Grain 277: Dominant Phase - Forsterite, Size - 77
Grain 278: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 279: Dominant Phase - Diopside   CaMgSi2O6, Size - 9
Grain 280: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 281: Dominant Phase - Forsterite, Size - 10
Grain 282: Dominant Phase - Diopside   CaMgSi2O6, Size - 28
Grain 283: Dominant Phase - Chromite, Size - 3
Grain 284: Dominant Phase - Chromite, Size - 57
Grain 285: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 286: Dominant Phase - Forsterite, Size - 375
Grain 287: Dominant Phase - Enstatite  Opx AV77, Size - 9
Grain 288: Dominant Phase - Diopside   CaMgSi2O6, Size - 9
Grain 289: Dominant Phase - Diopside   CaMgSi2O6, Size - 98
Grain 290: Dominant Phase - Forsterite, Size - 1
Grain 291: Dominant Phase - Forsterite, Size - 20
Grain 292: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 293: Dominant Phase - Diopside   CaMgSi2O6, Size - 69
Grain 294: Dominant Phase - Forsterite, Size - 55
Grain 295: Dominant Phase - Forsterite, Size - 5
Grain 296: Dominant Phase - Diopside   CaMgSi2O6, Size - 156
Grain 297: Dominant Phase - Forsterite, Size - 20
Grain 298: Dominant Phase - Enstatite  Opx AV77, Size - 680
Grain 299: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 300: Dominant Phase - Forsterite, Size - 2
Grain 301: Dominant Phase - Diopside   CaMgSi2O6, Size - 63
Grain 302: Dominant Phase - Forsterite, Size - 15
Grain 303: Dominant Phase - Forsterite, Size - 383
Grain 304: Dominant Phase - Enstatite  Opx AV77, Size - 5
Grain 305: Dominant Phase - Forsterite, Size - 34
Grain 306: Dominant Phase - Forsterite, Size - 17
Grain 307: Dominant Phase - Forsterite, Size - 6
Grain 308: Dominant Phase - Forsterite, Size - 2
Grain 309: Dominant Phase - Forsterite, Size - 3
Grain 310: Dominant Phase - Enstatite  Opx AV77, Size - 32
Grain 311: Dominant Phase - Diopside   CaMgSi2O6, Size - 150
Grain 312: Dominant Phase - Forsterite, Size - 4
Grain 313: Dominant Phase - Forsterite, Size - 30
Grain 314: Dominant Phase - Forsterite, Size - 2
Grain 315: Dominant Phase - Forsterite, Size - 1
Grain 316: Dominant Phase - Forsterite, Size - 2
Grain 317: Dominant Phase - Forsterite, Size - 123
Grain 318: Dominant Phase - Forsterite, Size - 5
Grain 319: Dominant Phase - Forsterite, Size - 216
Grain 320: Dominant Phase - Forsterite, Size - 28
Grain 321: Dominant Phase - Forsterite, Size - 12
Grain 322: Dominant Phase - Forsterite, Size - 1
Grain 323: Dominant Phase - Forsterite, Size - 179
Grain 324: Dominant Phase - Forsterite, Size - 1
Grain 325: Dominant Phase - Forsterite, Size - 7
Grain 326: Dominant Phase - Enstatite  Opx AV77, Size - 178
Grain 327: Dominant Phase - Forsterite, Size - 210
Grain 328: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 329: Dominant Phase - Diopside   CaMgSi2O6, Size - 8
Grain 330: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 331: Dominant Phase - Forsterite, Size - 332
Grain 332: Dominant Phase - Forsterite, Size - 492
Grain 333: Dominant Phase - Enstatite  Opx AV77, Size - 141
Grain 334: Dominant Phase - Chromite, Size - 20
Grain 335: Dominant Phase - Forsterite, Size - 79
Grain 336: Dominant Phase - Forsterite, Size - 85
Grain 337: Dominant Phase - Enstatite  Opx AV77, Size - 39
Grain 338: Dominant Phase - Forsterite, Size - 90
Grain 339: Dominant Phase - Forsterite, Size - 1
Grain 340: Dominant Phase - Forsterite, Size - 551
Grain 341: Dominant Phase - Forsterite, Size - 104
Grain 342: Dominant Phase - Forsterite, Size - 188
Grain 343: Dominant Phase - Forsterite, Size - 7
Grain 344: Dominant Phase - Forsterite, Size - 27
Grain 345: Dominant Phase - Enstatite  Opx AV77, Size - 146
Grain 346: Dominant Phase - Chromite, Size - 53
Grain 347: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 348: Dominant Phase - Diopside   CaMgSi2O6, Size - 36
Grain 349: Dominant Phase - Forsterite, Size - 12
Grain 350: Dominant Phase - Chromite, Size - 27
Grain 351: Dominant Phase - Forsterite, Size - 5
Grain 352: Dominant Phase - Forsterite, Size - 9
Grain 353: Dominant Phase - Forsterite, Size - 11
Grain 354: Dominant Phase - Diopside   CaMgSi2O6, Size - 52
Grain 355: Dominant Phase - Forsterite, Size - 1
Grain 356: Dominant Phase - Forsterite, Size - 181
Grain 357: Dominant Phase - Diopside   CaMgSi2O6, Size - 16
Grain 358: Dominant Phase - Enstatite  Opx AV77, Size - 4
Grain 359: Dominant Phase - Enstatite  Opx AV77, Size - 11
Grain 360: Dominant Phase - Chromite, Size - 56
Grain 361: Dominant Phase - Diopside   CaMgSi2O6, Size - 39
Grain 362: Dominant Phase - Diopside   CaMgSi2O6, Size - 8
Grain 363: Dominant Phase - Chromite, Size - 1
Grain 364: Dominant Phase - Forsterite, Size - 18
Grain 365: Dominant Phase - Forsterite, Size - 36
Grain 366: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 367: Dominant Phase - Forsterite, Size - 11
Grain 368: Dominant Phase - Forsterite, Size - 203
Grain 369: Dominant Phase - Forsterite, Size - 69
Grain 370: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 371: Dominant Phase - Forsterite, Size - 46
Grain 372: Dominant Phase - Chromite, Size - 21
Grain 373: Dominant Phase - Diopside   CaMgSi2O6, Size - 19
Grain 374: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 375: Dominant Phase - Chromite, Size - 1
Grain 376: Dominant Phase - Forsterite, Size - 88
Grain 377: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 378: Dominant Phase - Forsterite, Size - 1
Grain 379: Dominant Phase - Forsterite, Size - 12
Grain 380: Dominant Phase - Diopside   CaMgSi2O6, Size - 14
Grain 381: Dominant Phase - Diopside   CaMgSi2O6, Size - 166
Grain 382: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 383: Dominant Phase - Forsterite, Size - 37
Grain 384: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 385: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 386: Dominant Phase - Forsterite, Size - 272
Grain 387: Dominant Phase - Forsterite, Size - 32
Grain 388: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 389: Dominant Phase - Forsterite, Size - 5
Grain 390: Dominant Phase - Diopside   CaMgSi2O6, Size - 47
Grain 391: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 392: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 393: Dominant Phase - Forsterite, Size - 10
Grain 394: Dominant Phase - Diopside   CaMgSi2O6, Size - 60
Grain 395: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 396: Dominant Phase - Forsterite, Size - 251
Grain 397: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 398: Dominant Phase - Forsterite, Size - 53
Grain 399: Dominant Phase - Forsterite, Size - 373
Grain 400: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 401: Dominant Phase - Forsterite, Size - 1
Grain 402: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 403: Dominant Phase - Diopside   CaMgSi2O6, Size - 701
Grain 404: Dominant Phase - Forsterite, Size - 22
Grain 405: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 406: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 407: Dominant Phase - Forsterite, Size - 87
Grain 408: Dominant Phase - Enstatite  Opx AV77, Size - 29
Grain 409: Dominant Phase - Diopside   CaMgSi2O6, Size - 28
Grain 410: Dominant Phase - Forsterite, Size - 16
Grain 411: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 412: Dominant Phase - Diopside   CaMgSi2O6, Size - 17
Grain 413: Dominant Phase - Diopside   CaMgSi2O6, Size - 41
Grain 414: Dominant Phase - Diopside   CaMgSi2O6, Size - 83
Grain 415: Dominant Phase - Forsterite, Size - 9
Grain 416: Dominant Phase - Diopside   CaMgSi2O6, Size - 91
Grain 417: Dominant Phase - Forsterite, Size - 40
Grain 418: Dominant Phase - Forsterite, Size - 1
Grain 419: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 420: Dominant Phase - Enstatite  Opx AV77, Size - 3
Grain 421: Dominant Phase - Forsterite, Size - 1
Grain 422: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 423: Dominant Phase - Diopside   CaMgSi2O6, Size - 82
Grain 424: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 425: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 426: Dominant Phase - Diopside   CaMgSi2O6, Size - 16
Grain 427: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 428: Dominant Phase - Diopside   CaMgSi2O6, Size - 91
Grain 429: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 430: Dominant Phase - Forsterite, Size - 44
Grain 431: Dominant Phase - Enstatite  Opx AV77, Size - 28
Grain 432: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 433: Dominant Phase - Diopside   CaMgSi2O6, Size - 32
Grain 434: Dominant Phase - Diopside   CaMgSi2O6, Size - 8
Grain 435: Dominant Phase - Forsterite, Size - 138
Grain 436: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 437: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 438: Dominant Phase - Diopside   CaMgSi2O6, Size - 185
Grain 439: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 440: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 441: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 442: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 443: Dominant Phase - Forsterite, Size - 22
Grain 444: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 445: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 446: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 447: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 448: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 449: Dominant Phase - Enstatite  Opx AV77, Size - 230
Grain 450: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 451: Dominant Phase - Diopside   CaMgSi2O6, Size - 35
Grain 452: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 453: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 454: Dominant Phase - Chromite, Size - 1
Grain 455: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 456: Dominant Phase - Forsterite, Size - 17
Grain 457: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 458: Dominant Phase - Forsterite, Size - 192
Grain 459: Dominant Phase - Diopside   CaMgSi2O6, Size - 8
Grain 460: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 461: Dominant Phase - Forsterite, Size - 36
Grain 462: Dominant Phase - Forsterite, Size - 16
Grain 463: Dominant Phase - Forsterite, Size - 23
Grain 464: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 465: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 466: Dominant Phase - Diopside   CaMgSi2O6, Size - 77
Grain 467: Dominant Phase - Chromite, Size - 140
Grain 468: Dominant Phase - Chromite, Size - 77
Grain 469: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 470: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 471: Dominant Phase - Chromite, Size - 18
Grain 472: Dominant Phase - Diopside   CaMgSi2O6, Size - 62
Grain 473: Dominant Phase - Diopside   CaMgSi2O6, Size - 35
Grain 474: Dominant Phase - Diopside   CaMgSi2O6, Size - 52
Grain 475: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 476: Dominant Phase - Diopside   CaMgSi2O6, Size - 4
Grain 477: Dominant Phase - Forsterite, Size - 1
Grain 478: Dominant Phase - Forsterite, Size - 149
Grain 479: Dominant Phase - Chromite, Size - 7
Grain 480: Dominant Phase - Diopside   CaMgSi2O6, Size - 141
Grain 481: Dominant Phase - Forsterite, Size - 1
Grain 482: Dominant Phase - Diopside   CaMgSi2O6, Size - 75
Grain 483: Dominant Phase - Forsterite, Size - 1
Grain 484: Dominant Phase - Chromite, Size - 1
Grain 485: Dominant Phase - Diopside   CaMgSi2O6, Size - 104
Grain 486: Dominant Phase - Chromite, Size - 2
Grain 487: Dominant Phase - Forsterite, Size - 136
Grain 488: Dominant Phase - Chromite, Size - 1
Grain 489: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 490: Dominant Phase - Chromite, Size - 1
Grain 491: Dominant Phase - Forsterite, Size - 133
Grain 492: Dominant Phase - Diopside   CaMgSi2O6, Size - 12
Grain 493: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 494: Dominant Phase - Enstatite  Opx AV77, Size - 16
Grain 495: Dominant Phase - Enstatite  Opx AV77, Size - 38
Grain 496: Dominant Phase - Chromite, Size - 5
Grain 497: Dominant Phase - Forsterite, Size - 21
Grain 498: Dominant Phase - Forsterite, Size - 5
Grain 499: Dominant Phase - Forsterite, Size - 24
Grain 500: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 501: Dominant Phase - Forsterite, Size - 3
Grain 502: Dominant Phase - Forsterite, Size - 2
Grain 503: Dominant Phase - Diopside   CaMgSi2O6, Size - 73
Grain 504: Dominant Phase - Chromite, Size - 1
Grain 505: Dominant Phase - Forsterite, Size - 1
Grain 506: Dominant Phase - Enstatite  Opx AV77, Size - 91
Grain 507: Dominant Phase - Forsterite, Size - 37
Grain 508: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 509: Dominant Phase - Diopside   CaMgSi2O6, Size - 9
Grain 510: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 511: Dominant Phase - Diopside   CaMgSi2O6, Size - 44
Grain 512: Dominant Phase - Chromite, Size - 4
Grain 513: Dominant Phase - Diopside   CaMgSi2O6, Size - 137
Grain 514: Dominant Phase - Diopside   CaMgSi2O6, Size - 56
Grain 515: Dominant Phase - Diopside   CaMgSi2O6, Size - 40
Grain 516: Dominant Phase - Chromite, Size - 12
Grain 517: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 518: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 519: Dominant Phase - Chromite, Size - 4
Grain 520: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 521: Dominant Phase - Diopside   CaMgSi2O6, Size - 34
Grain 522: Dominant Phase - Forsterite, Size - 17
Grain 523: Dominant Phase - Enstatite  Opx AV77, Size - 383
Grain 524: Dominant Phase - Forsterite, Size - 3
Grain 525: Dominant Phase - Diopside   CaMgSi2O6, Size - 93
Grain 526: Dominant Phase - Forsterite, Size - 545
Grain 527: Dominant Phase - Diopside   CaMgSi2O6, Size - 6
Grain 528: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 529: Dominant Phase - Forsterite, Size - 7
Grain 530: Dominant Phase - Diopside   CaMgSi2O6, Size - 63
Grain 531: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 532: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 533: Dominant Phase - Diopside   CaMgSi2O6, Size - 14
Grain 534: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 535: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 536: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 537: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 538: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 539: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 540: Dominant Phase - Diopside   CaMgSi2O6, Size - 33
Grain 541: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 542: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 543: Dominant Phase - Diopside   CaMgSi2O6, Size - 111
Grain 544: Dominant Phase - Diopside   CaMgSi2O6, Size - 34
Grain 545: Dominant Phase - Forsterite, Size - 1
Grain 546: Dominant Phase - Forsterite, Size - 4
Grain 547: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 548: Dominant Phase - Chromite, Size - 12
Grain 549: Dominant Phase - Diopside   CaMgSi2O6, Size - 29
Grain 550: Dominant Phase - Forsterite, Size - 26
Grain 551: Dominant Phase - Diopside   CaMgSi2O6, Size - 29
Grain 552: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 553: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 554: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 555: Dominant Phase - Forsterite, Size - 11
Grain 556: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 557: Dominant Phase - Diopside   CaMgSi2O6, Size - 22
Grain 558: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 559: Dominant Phase - Forsterite, Size - 16
Grain 560: Dominant Phase - Forsterite, Size - 38
Grain 561: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 562: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 563: Dominant Phase - Diopside   CaMgSi2O6, Size - 3
Grain 564: Dominant Phase - Forsterite, Size - 2
Grain 565: Dominant Phase - Enstatite  Opx AV77, Size - 66
Grain 566: Dominant Phase - Forsterite, Size - 2
Grain 567: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 568: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 569: Dominant Phase - Chromite, Size - 10
Grain 570: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 571: Dominant Phase - Diopside   CaMgSi2O6, Size - 53
Grain 572: Dominant Phase - Forsterite, Size - 3
Grain 573: Dominant Phase - Chromite, Size - 5
Grain 574: Dominant Phase - Forsterite, Size - 30
Grain 575: Dominant Phase - Diopside   CaMgSi2O6, Size - 42
Grain 576: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 577: Dominant Phase - Diopside   CaMgSi2O6, Size - 30
Grain 578: Dominant Phase - Enstatite  Opx AV77, Size - 6
Grain 579: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 580: Dominant Phase - Diopside   CaMgSi2O6, Size - 34
Grain 581: Dominant Phase - Chromite, Size - 1
Grain 582: Dominant Phase - Chromite, Size - 11
Grain 583: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 584: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 585: Dominant Phase - Diopside   CaMgSi2O6, Size - 80
Grain 586: Dominant Phase - Enstatite  Opx AV77, Size - 47
Grain 587: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 588: Dominant Phase - Chromite, Size - 5
Grain 589: Dominant Phase - Chromite, Size - 10
Grain 590: Dominant Phase - Chromite, Size - 1
Grain 591: Dominant Phase - Diopside   CaMgSi2O6, Size - 35
Grain 592: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 593: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 594: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 595: Dominant Phase - Forsterite, Size - 39
Grain 596: Dominant Phase - Enstatite  Opx AV77, Size - 1
Grain 597: Dominant Phase - Forsterite, Size - 1
Grain 598: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 599: Dominant Phase - Forsterite, Size - 1
Grain 600: Dominant Phase - Chromite, Size - 23
Grain 601: Dominant Phase - Chromite, Size - 10
Grain 602: Dominant Phase - Diopside   CaMgSi2O6, Size - 2
Grain 603: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 604: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 605: Dominant Phase - Forsterite, Size - 4
Grain 606: Dominant Phase - Diopside   CaMgSi2O6, Size - 10
Grain 607: Dominant Phase - Chromite, Size - 28
Grain 608: Dominant Phase - Diopside   CaMgSi2O6, Size - 7
Grain 609: Dominant Phase - Forsterite, Size - 11
Grain 610: Dominant Phase - Diopside   CaMgSi2O6, Size - 1
Grain 611: Dominant Phase - Diopside   CaMgSi2O6, Size - 5
Grain 612: Dominant Phase - Enstatite  Opx AV77, Size - 4
Grain 613: Dominant Phase - Forsterite, Size - 4
Grain 614: Dominant Phase - Forsterite, Size - 1

Removing grains with grain size less than 10 pixels

Now, once the grains are calculated above, the next step is to remove grains whose size is less than 10 pixels

filtered_df_grain_boundary = ebsdfile.filter_by_grain_size(df_grain_boundary, phases_names, min_grain_size=10)

Grain 0: Dominant Phase - Enstatite  Opx AV77, Size - 694
Grain 1: Dominant Phase - Forsterite, Size - 72
Grain 2: Dominant Phase - Diopside   CaMgSi2O6, Size - 47
Grain 4: Dominant Phase - Diopside   CaMgSi2O6, Size - 65
Grain 5: Dominant Phase - Forsterite, Size - 334
Grain 6: Dominant Phase - Diopside   CaMgSi2O6, Size - 12
Grain 7: Dominant Phase - Forsterite, Size - 1298
Grain 8: Dominant Phase - Forsterite, Size - 408
Grain 9: Dominant Phase - Forsterite, Size - 31
Grain 10: Dominant Phase - Forsterite, Size - 156
Grain 11: Dominant Phase - Forsterite, Size - 11
Grain 12: Dominant Phase - Forsterite, Size - 202
Grain 13: Dominant Phase - Enstatite  Opx AV77, Size - 58
Grain 14: Dominant Phase - Forsterite, Size - 838
Grain 15: Dominant Phase - Enstatite  Opx AV77, Size - 28
Grain 16: Dominant Phase - Forsterite, Size - 1010
Grain 17: Dominant Phase - Forsterite, Size - 17
Grain 18: Dominant Phase - Diopside   CaMgSi2O6, Size - 20
Grain 19: Dominant Phase - Forsterite, Size - 61
Grain 20: Dominant Phase - Forsterite, Size - 431
Grain 21: Dominant Phase - Chromite, Size - 40
Grain 22: Dominant Phase - Diopside   CaMgSi2O6, Size - 36
Grain 23: Dominant Phase - Forsterite, Size - 842
Grain 26: Dominant Phase - Enstatite  Opx AV77, Size - 389
Grain 28: Dominant Phase - Forsterite, Size - 485
Grain 29: Dominant Phase - Forsterite, Size - 196
Grain 30: Dominant Phase - Forsterite, Size - 162
Grain 32: Dominant Phase - Forsterite, Size - 519
Grain 33: Dominant Phase - Forsterite, Size - 603
Grain 34: Dominant Phase - Forsterite, Size - 536
Grain 35: Dominant Phase - Enstatite  Opx AV77, Size - 406
Grain 36: Dominant Phase - Enstatite  Opx AV77, Size - 253
Grain 37: Dominant Phase - Forsterite, Size - 46
Grain 38: Dominant Phase - Diopside   CaMgSi2O6, Size - 237
Grain 39: Dominant Phase - Forsterite, Size - 59
Grain 42: Dominant Phase - Enstatite  Opx AV77, Size - 200
Grain 43: Dominant Phase - Diopside   CaMgSi2O6, Size - 65
Grain 44: Dominant Phase - Forsterite, Size - 19
Grain 46: Dominant Phase - Forsterite, Size - 21
Grain 48: Dominant Phase - Forsterite, Size - 13
Grain 49: Dominant Phase - Enstatite  Opx AV77, Size - 239
Grain 50: Dominant Phase - Forsterite, Size - 47
Grain 52: Dominant Phase - Forsterite, Size - 172
Grain 53: Dominant Phase - Diopside   CaMgSi2O6, Size - 111
Grain 54: Dominant Phase - Forsterite, Size - 58
Grain 55: Dominant Phase - Forsterite, Size - 258
Grain 56: Dominant Phase - Forsterite, Size - 31
Grain 57: Dominant Phase - Forsterite, Size - 12
Grain 59: Dominant Phase - Chromite, Size - 13
Grain 60: Dominant Phase - Diopside   CaMgSi2O6, Size - 56
Grain 62: Dominant Phase - Diopside   CaMgSi2O6, Size - 52
Grain 63: Dominant Phase - Diopside   CaMgSi2O6, Size - 198
Grain 64: Dominant Phase - Diopside   CaMgSi2O6, Size - 14
Grain 65: Dominant Phase - Forsterite, Size - 268
Grain 66: Dominant Phase - Forsterite, Size - 12
Grain 67: Dominant Phase - Diopside   CaMgSi2O6, Size - 852
Grain 69: Dominant Phase - Forsterite, Size - 378
Grain 70: Dominant Phase - Forsterite, Size - 590
Grain 71: Dominant Phase - Enstatite  Opx AV77, Size - 186
Grain 73: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 74: Dominant Phase - Forsterite, Size - 116
Grain 81: Dominant Phase - Forsterite, Size - 61
Grain 83: Dominant Phase - Forsterite, Size - 94
Grain 87: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 89: Dominant Phase - Forsterite, Size - 10
Grain 91: Dominant Phase - Diopside   CaMgSi2O6, Size - 240
Grain 95: Dominant Phase - Forsterite, Size - 121
Grain 97: Dominant Phase - Forsterite, Size - 10
Grain 98: Dominant Phase - Enstatite  Opx AV77, Size - 11
Grain 99: Dominant Phase - Forsterite, Size - 16
Grain 100: Dominant Phase - Forsterite, Size - 129
Grain 101: Dominant Phase - Diopside   CaMgSi2O6, Size - 84
Grain 102: Dominant Phase - Forsterite, Size - 11
Grain 103: Dominant Phase - Diopside   CaMgSi2O6, Size - 17
Grain 104: Dominant Phase - Diopside   CaMgSi2O6, Size - 29
Grain 105: Dominant Phase - Chromite, Size - 82
Grain 106: Dominant Phase - Forsterite, Size - 24
Grain 108: Dominant Phase - Diopside   CaMgSi2O6, Size - 46
Grain 110: Dominant Phase - Forsterite, Size - 12
Grain 111: Dominant Phase - Forsterite, Size - 139
Grain 114: Dominant Phase - Forsterite, Size - 69
Grain 116: Dominant Phase - Diopside   CaMgSi2O6, Size - 24
Grain 117: Dominant Phase - Forsterite, Size - 1174
Grain 120: Dominant Phase - Forsterite, Size - 106
Grain 121: Dominant Phase - Enstatite  Opx AV77, Size - 291
Grain 124: Dominant Phase - Diopside   CaMgSi2O6, Size - 42
Grain 125: Dominant Phase - Enstatite  Opx AV77, Size - 40
Grain 126: Dominant Phase - Forsterite, Size - 34
Grain 127: Dominant Phase - Forsterite, Size - 350
Grain 128: Dominant Phase - Forsterite, Size - 15
Grain 130: Dominant Phase - Forsterite, Size - 45
Grain 133: Dominant Phase - Forsterite, Size - 162
Grain 134: Dominant Phase - Forsterite, Size - 195
Grain 135: Dominant Phase - Enstatite  Opx AV77, Size - 157
Grain 136: Dominant Phase - Forsterite, Size - 11
Grain 138: Dominant Phase - Diopside   CaMgSi2O6, Size - 185
Grain 139: Dominant Phase - Chromite, Size - 17
Grain 140: Dominant Phase - Diopside   CaMgSi2O6, Size - 54
Grain 141: Dominant Phase - Diopside   CaMgSi2O6, Size - 81
Grain 142: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 143: Dominant Phase - Forsterite, Size - 231
Grain 144: Dominant Phase - Enstatite  Opx AV77, Size - 199
Grain 147: Dominant Phase - Diopside   CaMgSi2O6, Size - 163
Grain 148: Dominant Phase - Forsterite, Size - 56
Grain 149: Dominant Phase - Forsterite, Size - 1050
Grain 150: Dominant Phase - Forsterite, Size - 46
Grain 151: Dominant Phase - Diopside   CaMgSi2O6, Size - 11
Grain 153: Dominant Phase - Forsterite, Size - 125
Grain 154: Dominant Phase - Forsterite, Size - 35
Grain 155: Dominant Phase - Diopside   CaMgSi2O6, Size - 199
Grain 156: Dominant Phase - Forsterite, Size - 74
Grain 157: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 159: Dominant Phase - Forsterite, Size - 31
Grain 160: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 161: Dominant Phase - Diopside   CaMgSi2O6, Size - 13
Grain 163: Dominant Phase - Diopside   CaMgSi2O6, Size - 79
Grain 164: Dominant Phase - Forsterite, Size - 340
Grain 165: Dominant Phase - Forsterite, Size - 84
Grain 166: Dominant Phase - Forsterite, Size - 237
Grain 167: Dominant Phase - Forsterite, Size - 326
Grain 168: Dominant Phase - Diopside   CaMgSi2O6, Size - 72
Grain 171: Dominant Phase - Forsterite, Size - 14
Grain 173: Dominant Phase - Forsterite, Size - 328
Grain 174: Dominant Phase - Forsterite, Size - 98
Grain 175: Dominant Phase - Forsterite, Size - 30
Grain 176: Dominant Phase - Diopside   CaMgSi2O6, Size - 22
Grain 177: Dominant Phase - Diopside   CaMgSi2O6, Size - 87
Grain 178: Dominant Phase - Forsterite, Size - 242
Grain 179: Dominant Phase - Forsterite, Size - 278
Grain 183: Dominant Phase - Forsterite, Size - 799
Grain 184: Dominant Phase - Diopside   CaMgSi2O6, Size - 64
Grain 186: Dominant Phase - Forsterite, Size - 42
Grain 187: Dominant Phase - Forsterite, Size - 41
Grain 188: Dominant Phase - Diopside   CaMgSi2O6, Size - 53
Grain 189: Dominant Phase - Forsterite, Size - 531
Grain 190: Dominant Phase - Diopside   CaMgSi2O6, Size - 63
Grain 193: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 194: Dominant Phase - Forsterite, Size - 13
Grain 195: Dominant Phase - Forsterite, Size - 24
Grain 196: Dominant Phase - Forsterite, Size - 121
Grain 198: Dominant Phase - Diopside   CaMgSi2O6, Size - 25
Grain 199: Dominant Phase - Diopside   CaMgSi2O6, Size - 82
Grain 201: Dominant Phase - Forsterite, Size - 21
Grain 202: Dominant Phase - Forsterite, Size - 80
Grain 203: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 205: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 206: Dominant Phase - Forsterite, Size - 70
Grain 207: Dominant Phase - Forsterite, Size - 230
Grain 208: Dominant Phase - Diopside   CaMgSi2O6, Size - 13
Grain 210: Dominant Phase - Forsterite, Size - 340
Grain 212: Dominant Phase - Diopside   CaMgSi2O6, Size - 37
Grain 213: Dominant Phase - Forsterite, Size - 218
Grain 215: Dominant Phase - Forsterite, Size - 28
Grain 218: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 227: Dominant Phase - Diopside   CaMgSi2O6, Size - 55
Grain 228: Dominant Phase - Forsterite, Size - 787
Grain 229: Dominant Phase - Forsterite, Size - 98
Grain 231: Dominant Phase - Forsterite, Size - 157
Grain 235: Dominant Phase - Diopside   CaMgSi2O6, Size - 20
Grain 237: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 238: Dominant Phase - Forsterite, Size - 17
Grain 240: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 242: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 245: Dominant Phase - Enstatite  Opx AV77, Size - 131
Grain 247: Dominant Phase - Forsterite, Size - 126
Grain 248: Dominant Phase - Forsterite, Size - 39
Grain 250: Dominant Phase - Diopside   CaMgSi2O6, Size - 160
Grain 251: Dominant Phase - Enstatite  Opx AV77, Size - 629
Grain 252: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 257: Dominant Phase - Forsterite, Size - 10
Grain 258: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 260: Dominant Phase - Diopside   CaMgSi2O6, Size - 10
Grain 261: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 263: Dominant Phase - Diopside   CaMgSi2O6, Size - 94
Grain 264: Dominant Phase - Forsterite, Size - 52
Grain 265: Dominant Phase - Forsterite, Size - 24
Grain 266: Dominant Phase - Forsterite, Size - 349
Grain 267: Dominant Phase - Forsterite, Size - 175
Grain 271: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 272: Dominant Phase - Forsterite, Size - 54
Grain 275: Dominant Phase - Forsterite, Size - 306
Grain 277: Dominant Phase - Forsterite, Size - 77
Grain 281: Dominant Phase - Forsterite, Size - 10
Grain 282: Dominant Phase - Diopside   CaMgSi2O6, Size - 28
Grain 284: Dominant Phase - Chromite, Size - 57
Grain 286: Dominant Phase - Forsterite, Size - 375
Grain 289: Dominant Phase - Diopside   CaMgSi2O6, Size - 98
Grain 291: Dominant Phase - Forsterite, Size - 20
Grain 293: Dominant Phase - Diopside   CaMgSi2O6, Size - 69
Grain 294: Dominant Phase - Forsterite, Size - 55
Grain 296: Dominant Phase - Diopside   CaMgSi2O6, Size - 156
Grain 297: Dominant Phase - Forsterite, Size - 20
Grain 298: Dominant Phase - Enstatite  Opx AV77, Size - 680
Grain 299: Dominant Phase - Diopside   CaMgSi2O6, Size - 27
Grain 301: Dominant Phase - Diopside   CaMgSi2O6, Size - 63
Grain 302: Dominant Phase - Forsterite, Size - 15
Grain 303: Dominant Phase - Forsterite, Size - 383
Grain 305: Dominant Phase - Forsterite, Size - 34
Grain 306: Dominant Phase - Forsterite, Size - 17
Grain 310: Dominant Phase - Enstatite  Opx AV77, Size - 32
Grain 311: Dominant Phase - Diopside   CaMgSi2O6, Size - 150
Grain 313: Dominant Phase - Forsterite, Size - 30
Grain 317: Dominant Phase - Forsterite, Size - 123
Grain 319: Dominant Phase - Forsterite, Size - 216
Grain 320: Dominant Phase - Forsterite, Size - 28
Grain 321: Dominant Phase - Forsterite, Size - 12
Grain 323: Dominant Phase - Forsterite, Size - 179
Grain 326: Dominant Phase - Enstatite  Opx AV77, Size - 178
Grain 327: Dominant Phase - Forsterite, Size - 210
Grain 331: Dominant Phase - Forsterite, Size - 332
Grain 332: Dominant Phase - Forsterite, Size - 492
Grain 333: Dominant Phase - Enstatite  Opx AV77, Size - 141
Grain 334: Dominant Phase - Chromite, Size - 20
Grain 335: Dominant Phase - Forsterite, Size - 79
Grain 336: Dominant Phase - Forsterite, Size - 85
Grain 337: Dominant Phase - Enstatite  Opx AV77, Size - 39
Grain 338: Dominant Phase - Forsterite, Size - 90
Grain 340: Dominant Phase - Forsterite, Size - 551
Grain 341: Dominant Phase - Forsterite, Size - 104
Grain 342: Dominant Phase - Forsterite, Size - 188
Grain 344: Dominant Phase - Forsterite, Size - 27
Grain 345: Dominant Phase - Enstatite  Opx AV77, Size - 146
Grain 346: Dominant Phase - Chromite, Size - 53
Grain 348: Dominant Phase - Diopside   CaMgSi2O6, Size - 36
Grain 349: Dominant Phase - Forsterite, Size - 12
Grain 350: Dominant Phase - Chromite, Size - 27
Grain 353: Dominant Phase - Forsterite, Size - 11
Grain 354: Dominant Phase - Diopside   CaMgSi2O6, Size - 52
Grain 356: Dominant Phase - Forsterite, Size - 181
Grain 357: Dominant Phase - Diopside   CaMgSi2O6, Size - 16
Grain 359: Dominant Phase - Enstatite  Opx AV77, Size - 11
Grain 360: Dominant Phase - Chromite, Size - 56
Grain 361: Dominant Phase - Diopside   CaMgSi2O6, Size - 39
Grain 364: Dominant Phase - Forsterite, Size - 18
Grain 365: Dominant Phase - Forsterite, Size - 36
Grain 367: Dominant Phase - Forsterite, Size - 11
Grain 368: Dominant Phase - Forsterite, Size - 203
Grain 369: Dominant Phase - Forsterite, Size - 69
Grain 371: Dominant Phase - Forsterite, Size - 46
Grain 372: Dominant Phase - Chromite, Size - 21
Grain 373: Dominant Phase - Diopside   CaMgSi2O6, Size - 19
Grain 376: Dominant Phase - Forsterite, Size - 88
Grain 377: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 379: Dominant Phase - Forsterite, Size - 12
Grain 380: Dominant Phase - Diopside   CaMgSi2O6, Size - 14
Grain 381: Dominant Phase - Diopside   CaMgSi2O6, Size - 166
Grain 382: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 383: Dominant Phase - Forsterite, Size - 37
Grain 384: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 386: Dominant Phase - Forsterite, Size - 272
Grain 387: Dominant Phase - Forsterite, Size - 32
Grain 390: Dominant Phase - Diopside   CaMgSi2O6, Size - 47
Grain 393: Dominant Phase - Forsterite, Size - 10
Grain 394: Dominant Phase - Diopside   CaMgSi2O6, Size - 60
Grain 395: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 396: Dominant Phase - Forsterite, Size - 251
Grain 398: Dominant Phase - Forsterite, Size - 53
Grain 399: Dominant Phase - Forsterite, Size - 373
Grain 403: Dominant Phase - Diopside   CaMgSi2O6, Size - 701
Grain 404: Dominant Phase - Forsterite, Size - 22
Grain 405: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 407: Dominant Phase - Forsterite, Size - 87
Grain 408: Dominant Phase - Enstatite  Opx AV77, Size - 29
Grain 409: Dominant Phase - Diopside   CaMgSi2O6, Size - 28
Grain 410: Dominant Phase - Forsterite, Size - 16
Grain 412: Dominant Phase - Diopside   CaMgSi2O6, Size - 17
Grain 413: Dominant Phase - Diopside   CaMgSi2O6, Size - 41
Grain 414: Dominant Phase - Diopside   CaMgSi2O6, Size - 83
Grain 416: Dominant Phase - Diopside   CaMgSi2O6, Size - 91
Grain 417: Dominant Phase - Forsterite, Size - 40
Grain 419: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 422: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 423: Dominant Phase - Diopside   CaMgSi2O6, Size - 82
Grain 425: Dominant Phase - Diopside   CaMgSi2O6, Size - 15
Grain 426: Dominant Phase - Diopside   CaMgSi2O6, Size - 16
Grain 428: Dominant Phase - Diopside   CaMgSi2O6, Size - 91
Grain 430: Dominant Phase - Forsterite, Size - 44
Grain 431: Dominant Phase - Enstatite  Opx AV77, Size - 28
Grain 433: Dominant Phase - Diopside   CaMgSi2O6, Size - 32
Grain 435: Dominant Phase - Forsterite, Size - 138
Grain 438: Dominant Phase - Diopside   CaMgSi2O6, Size - 185
Grain 443: Dominant Phase - Forsterite, Size - 22
Grain 449: Dominant Phase - Enstatite  Opx AV77, Size - 230
Grain 451: Dominant Phase - Diopside   CaMgSi2O6, Size - 35
Grain 456: Dominant Phase - Forsterite, Size - 17
Grain 458: Dominant Phase - Forsterite, Size - 192
Grain 461: Dominant Phase - Forsterite, Size - 36
Grain 462: Dominant Phase - Forsterite, Size - 16
Grain 463: Dominant Phase - Forsterite, Size - 23
Grain 466: Dominant Phase - Diopside   CaMgSi2O6, Size - 77
Grain 467: Dominant Phase - Chromite, Size - 140
Grain 468: Dominant Phase - Chromite, Size - 77
Grain 471: Dominant Phase - Chromite, Size - 18
Grain 472: Dominant Phase - Diopside   CaMgSi2O6, Size - 62
Grain 473: Dominant Phase - Diopside   CaMgSi2O6, Size - 35
Grain 474: Dominant Phase - Diopside   CaMgSi2O6, Size - 52
Grain 478: Dominant Phase - Forsterite, Size - 149
Grain 480: Dominant Phase - Diopside   CaMgSi2O6, Size - 141
Grain 482: Dominant Phase - Diopside   CaMgSi2O6, Size - 75
Grain 485: Dominant Phase - Diopside   CaMgSi2O6, Size - 104
Grain 487: Dominant Phase - Forsterite, Size - 136
Grain 491: Dominant Phase - Forsterite, Size - 133
Grain 492: Dominant Phase - Diopside   CaMgSi2O6, Size - 12
Grain 494: Dominant Phase - Enstatite  Opx AV77, Size - 16
Grain 495: Dominant Phase - Enstatite  Opx AV77, Size - 38
Grain 497: Dominant Phase - Forsterite, Size - 21
Grain 499: Dominant Phase - Forsterite, Size - 24
Grain 503: Dominant Phase - Diopside   CaMgSi2O6, Size - 73
Grain 506: Dominant Phase - Enstatite  Opx AV77, Size - 91
Grain 507: Dominant Phase - Forsterite, Size - 37
Grain 511: Dominant Phase - Diopside   CaMgSi2O6, Size - 44
Grain 513: Dominant Phase - Diopside   CaMgSi2O6, Size - 137
Grain 514: Dominant Phase - Diopside   CaMgSi2O6, Size - 56
Grain 515: Dominant Phase - Diopside   CaMgSi2O6, Size - 40
Grain 516: Dominant Phase - Chromite, Size - 12
Grain 521: Dominant Phase - Diopside   CaMgSi2O6, Size - 34
Grain 522: Dominant Phase - Forsterite, Size - 17
Grain 523: Dominant Phase - Enstatite  Opx AV77, Size - 383
Grain 525: Dominant Phase - Diopside   CaMgSi2O6, Size - 93
Grain 526: Dominant Phase - Forsterite, Size - 545
Grain 530: Dominant Phase - Diopside   CaMgSi2O6, Size - 63
Grain 531: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 533: Dominant Phase - Diopside   CaMgSi2O6, Size - 14
Grain 540: Dominant Phase - Diopside   CaMgSi2O6, Size - 33
Grain 543: Dominant Phase - Diopside   CaMgSi2O6, Size - 111
Grain 544: Dominant Phase - Diopside   CaMgSi2O6, Size - 34
Grain 548: Dominant Phase - Chromite, Size - 12
Grain 549: Dominant Phase - Diopside   CaMgSi2O6, Size - 29
Grain 550: Dominant Phase - Forsterite, Size - 26
Grain 551: Dominant Phase - Diopside   CaMgSi2O6, Size - 29
Grain 552: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 554: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 555: Dominant Phase - Forsterite, Size - 11
Grain 556: Dominant Phase - Diopside   CaMgSi2O6, Size - 18
Grain 557: Dominant Phase - Diopside   CaMgSi2O6, Size - 22
Grain 559: Dominant Phase - Forsterite, Size - 16
Grain 560: Dominant Phase - Forsterite, Size - 38
Grain 561: Dominant Phase - Diopside   CaMgSi2O6, Size - 23
Grain 565: Dominant Phase - Enstatite  Opx AV77, Size - 66
Grain 569: Dominant Phase - Chromite, Size - 10
Grain 571: Dominant Phase - Diopside   CaMgSi2O6, Size - 53
Grain 574: Dominant Phase - Forsterite, Size - 30
Grain 575: Dominant Phase - Diopside   CaMgSi2O6, Size - 42
Grain 577: Dominant Phase - Diopside   CaMgSi2O6, Size - 30
Grain 580: Dominant Phase - Diopside   CaMgSi2O6, Size - 34
Grain 582: Dominant Phase - Chromite, Size - 11
Grain 585: Dominant Phase - Diopside   CaMgSi2O6, Size - 80
Grain 586: Dominant Phase - Enstatite  Opx AV77, Size - 47
Grain 589: Dominant Phase - Chromite, Size - 10
Grain 591: Dominant Phase - Diopside   CaMgSi2O6, Size - 35
Grain 592: Dominant Phase - Diopside   CaMgSi2O6, Size - 45
Grain 595: Dominant Phase - Forsterite, Size - 39
Grain 600: Dominant Phase - Chromite, Size - 23
Grain 601: Dominant Phase - Chromite, Size - 10
Grain 603: Dominant Phase - Diopside   CaMgSi2O6, Size - 21
Grain 606: Dominant Phase - Diopside   CaMgSi2O6, Size - 10
Grain 607: Dominant Phase - Chromite, Size - 28
Grain 609: Dominant Phase - Forsterite, Size - 11

Now, lets check the number of pixels which were originally present in the EBSD data

len(df)

528500

The number of pixels originally in the ebsd data was 298068, now, lets check the length of the ebsd data after filtering with MAD>0.7

len(filtered_df)

510161

The number of pixels with MAD less than 0.7 is 235679. Now, lets check the length of pixels whose grains misorientation angle is less than 10 degrees

len(filtered_df_grain_boundary)

44398

ebsdfile.plot(filtered_df_grain_boundary)

Calculating Anisotropy

Now, after cleaning the data for MAD and grain size, we now calculate the seismic anisotropy of the above data

material_instance = Material()
rho_Fo = material_instance.load_density("Forsterite", 1.4, 1100)
rho_diop = material_instance.load_density("Diopside", 1.4, 1100)
rho_ens = material_instance.load_density("Enstatite", 1.4, 1100)
cij_Fo = material_instance.voigt_high_PT('Forsterite', PRESSURE = 1.4, TEMP = 1100)
cij_ens = material_instance.voigt_high_PT('Enstatite', PRESSURE = 1.4, TEMP = 1100)
cij_diop = material_instance.voigt_high_PT('Diopside', PRESSURE = 1.4, TEMP = 1100)

cij = [cij_Fo, cij_ens, cij_diop]
density = [rho_Fo, rho_ens, rho_diop]

forsterite = ebsdfile.get_euler_angles(phase = 1, data=filtered_df_grain_boundary)
enstatite = ebsdfile.get_euler_angles(phase = 2, data=filtered_df_grain_boundary)
diopside = ebsdfile.get_euler_angles(phase = 3, data=filtered_df_grain_boundary)
euler_angles = [forsterite, enstatite, diopside]

average_tensor, average_density = ebsdfile.get_anisotropy_for_ebsd(cij, euler_angles, density)
anis = Anisotropy(average_tensor*10**9, average_density)
anis.plot()

anis.anisotropy_values()

Max Vp:  8909.2511558578
Min Vp:  8016.424918867894
Max Vs1:  4932.585937221644
Min Vs1:  4777.78027777745
Max Vs2:  4812.697675943905
Min Vs2:  4558.366203280216
Max vs anisotropy percent:  6.873087103896376
Min vs anisotropy percent:  0.07689722204489084
P wave anisotropy percent:  10.549962471787126
S1 Wave anisotropy percent:  3.1884618152727127
S2 Wave anisotropy percent:  5.42801705209903
Velocity difference, maxdvs:  326.3719941181398
Vp/Vs1 ratio:  9.580521312998327
Mean vp:  8311.616111003932
Mean Vs1:  4865.174775177715
Mean Vs2:  4627.852257034686

{'maxvp': 8909.2511558578,
 'minvp': 8016.424918867894,
 'maxvs1': 4932.585937221644,
 'minvs1': 4777.78027777745,
 'maxvs2': 4812.697675943905,
 'minvs2': 4558.366203280216,
 'max_vs_anisotropy_percent': 6.873087103896376,
 'min_vs_anisotropy_percent': 0.07689722204489084,
 'p_wave_anisotropy_percent': 10.549962471787126,
 's1_wave_anisotropy_percent': 3.1884618152727127,
 's2_wave_anisotropy_percent': 5.42801705209903,
 'maxdvs': 326.3719941181398,
 'AVpVs1': 9.580521312998327}

ebsdfile.pf(df)

ebsdfile.pdf(df)

ebsdfile.ipf(df)

ebsdfile.ipf_colorkey(df)