In the section Denoising and Filling Missing Data we have discussed how to work with noisy EBSD data the contained non indexed pixels. Hereby, we made the assumption that the grid before and after the operations is the same.
In this section we explain how to interpolate an EBSD map at positions that do not belong to the grid. Lets us consider a simple example
mtexdata twins;
[grains, ebsd.grainId] = calcGrains(ebsd('indexed'));
% this command here is important :)
ebsd = ebsd.project2FundamentalRegion(grains);
plot(ebsd('indexed'),ebsd('indexed').orientations)
ebsd = EBSD
Phase Orientations Mineral Color Symmetry Crystal reference frame
0 46 (0.2%) notIndexed
1 22833 (100%) Magnesium LightSkyBlue 6/mmm X||a*, Y||b, Z||c*
Properties: bands, bc, bs, error, mad, x, y
Scan unit : um

Now we can use the command interp
to interpolate the orientation at arbitrary coordinates x
and y
.
x = 30.5; y = 5.5;
e1 = interp(ebsd,x,y)
e1 = EBSD
Phase Orientations Mineral Color Symmetry Crystal reference frame
1 1 (100%) Magnesium LightSkyBlue 6/mmm X||a*, Y||b, Z||c*
Id Phase phi1 Phi phi2 x y bands bc bs error mad grainId oldId
1 1 163 112 186 30.5 5.5 10 160 255 0 0.4 36 3109
Scan unit : um
By default the command interp
performs inverse distance interpolation. This is different to
e2 = ebsd('xy',x,y)
e2 = EBSD
Phase Orientations Mineral Color Symmetry Crystal reference frame
1 1 (100%) Magnesium LightSkyBlue 6/mmm X||a*, Y||b, Z||c*
Id Phase phi1 Phi phi2 bands bc bs error mad x y grainId
3109 1 163 112 186 10 160 255 0 0.4 30.6 5.4 36
Scan unit : um
which returns the nearest neighbour EBSD measurement. Lets have a look at the difference
angle(e1.orientations,e2.orientations)./degree
ans =
0.1508
Change of the measurement grid
The command interp
can be used to evaluate the EBSD map on a different grid, which might have higher or lower resolution or might even be rotated. Lets demonstrate this
% define a rotated coarse grid
omega = 5*degree;
[xmin, xmax, ymin, ymax] = ebsd.extent;
x = linspace(xmin-cos(omega)*ymax,xmax,100);
y = linspace(ymin-sin(omega)*xmax,ymax,50);
[x,y] = meshgrid(x,y);
xy = [cos(omega) -sin(omega); sin(omega) cos(omega) ] * [x(:),y(:)].';
% define the EBSD data set on this new grid
ebsdNewGrid = interp(ebsd,xy(1,:),xy(2,:))
% plot the regridded EBSD data set
plot(ebsdNewGrid('indexed'),ebsdNewGrid('indexed').orientations)
ebsdNewGrid = EBSD
Phase Orientations Mineral Color Symmetry Crystal reference frame
0 2624 (52%) notIndexed
1 2376 (48%) Magnesium LightSkyBlue 6/mmm X||a*, Y||b, Z||c*
Properties: x, y, bands, bc, bs, error, mad, grainId, oldId
Scan unit : um

Note, that we have not rotated the EBSD data but only the grid. All orientations as well as the position of all grains remains unchanged.
Another example is the change from a square to an hexagonal grid or vice versa. In this case the command interp
is implicitely called by the command gridify
. In order to demonstrate this functionality we start by EBSD data on a hex grid
mtexdata ferrite silent
plot(ebsd,ebsd.orientations)

and resample the data on a square grid. To do so we first define a smaller square unit cell corresponding to the hexagonal unit cell
% define a square unit cell
hexUnitCell = abs(round(ebsd.unitCell,4));
minUnit = min(hexUnitCell(hexUnitCell>0));
squnitCell = minUnit * [-1 -1;-1 1; 1 1; 1 -1];
% use the square unit cell for gridify
ebsd = ebsd.gridify('unitCell',squnitCell);
plot(ebsd,ebsd.orientations)
