cs = crystalSymmetry('1','mineral','test');
id = [...
0 0 0 0 0 0; ...
0 1 1 1 1 0; ...
0 1 1 1 1 0; ...
0 1 0 0 1 0; ...
0 1 0 0 1 0; ...
0 1 1 1 0 0; ...
0 0 0 0 0 0]==1;
rot = rotation.id(size(id));
rot(id) = rotation.rand;
ebsd = EBSDsquare(rot,2*ones(size(rot)),1:2,{'not indexed',cs},[1 1]);
plot(ebsd,ebsd.orientations)
grains = calcGrains(ebsd,'removeQuadruplePoints')
grains = grain2d
Phase Grains Pixels Mineral Symmetry Crystal reference frame
2 2 42 test 1 X||a*, Y||b*, Z||c*
boundary segments: 61 (54 µm)
inner boundary segments: 0 (0 µm)
triple points: 0
Id Phase Pixels phi1 Phi phi2 GOS
1 2 15 303 18 129 0
2 2 27 0 0 0 0
gB = grains.boundary
grains = merge(grains,gB(end),'calcMeanOrientation')
gB = grainBoundary
Segments length mineral 1 mineral 2
35 28 µm notIndexed test
26 26 µm test test
grains = grain2d
Phase Grains Pixels Mineral Symmetry Crystal reference frame
2 2 42 test 1 X||a*, Y||b*, Z||c*
boundary segments: 61 (54 µm)
inner boundary segments: 0 (0 µm)
triple points: 0
Id Phase Pixels phi1 Phi phi2 GOS
1 2 15 303 18 129 0
2 2 27 0 0 0 0
grains = smooth(grains,1,'moveTriplePoints')
grains = grain2d
Phase Grains Pixels Mineral Symmetry Crystal reference frame
2 2 42 test 1 X||a*, Y||b*, Z||c*
boundary segments: 61 (51 µm)
inner boundary segments: 0 (0 µm)
triple points: 0
Id Phase Pixels phi1 Phi phi2 GOS
1 2 15 303 18 129 0
2 2 27 0 0 0 0
plot(grains.boundary)
id = 2;
gB = grains(id).boundary;
plot(gB,gB.curvature(10),'linewidth',6)
mtexColorMap blue2red
setColorRange(0.5*[-1,1])