Create a new phase in the EBSD data for some grain conditions

[Matheus Brozovic Gariglio]

Hi everyone,

I'm trying to reconstruct the parent beta phase of my microstructure, but there are some grains that are been reconstructed that I don't want them to be, as you can see in the black squares in the image below.

I'm able to separate the grains that I want to be reconstructed with some conditions:

ebsd2 = ebsd('indexed');

cs_A = ebsd2('Alpha').CS;
cs_B = ebsd2('Beta').CS;

beta2alpha = orientation.Burger(cs_B,cs_A);
round2Miller(beta2alpha)

% reconstruct grains
[grains_1,ebsd2.grainId] = calcGrains(ebsd2,'threshold',1.5*degree,...
  'removeQuadruplePoints');
grains_1 = smooth(grains_1);

grains_A = grains_1('Alpha');
%retirer les grains trop petits (nettoyer la phase secondaire)
cond = grains_A.grainSize<90;
grains_A(cond)=[];
%retirer les grains avec un facteur de forme élevé
cond = grains_A.shapeFactor>=2 & grains_A.grainSize<750;
grains_A(cond)=[];
%retirer les grains lamellaires moyens
cond = grains_A.aspectRatio>=2.05 & grains_A.grainSize<=315 ;
grains_A(cond)=[];

%Séléctionner les autres grains qui ne sont pas primaires
grains_a = grains_1('Alpha');
grains_1 = grains_a(~ismember(grains_a.id,grains_A.id));

figure;plot(grains_1)

What I want to do now is to create a third phase in the EBSD data (basically a copy of the 'Alpha' phase) for the grains that I don't want to be reconstructed (grains_A in my code above) to be stored as another name ('AlphaP', for example).

It would be like in the example below.

ebsd2 = EBSD (show methods, plot)

 

 Phase  Orientations  Mineral  Color  Symmetry  Crystal reference frame

     1  11493 (1.4%)     Beta   Navy      m-3m                         

     2  807154 (99%)    Alpha    Red     6/mmm        X||a*, Y||b, Z||c

     3  ?????? (??%)    AlphaP  Red    6/mmm        X||a*, Y||b, Z||c

So when I do the other steps of calculation, those grains will not be taken into account as 'Alpha'.

Does anyone know how to do it, please?

Kind regards,

Matheus

[Rüdiger]

Hi,
I still don't know anything about titanium, but it's a beautiful
microstructure.

Does that example help?

mtexdata titanium

CS=ebsd.CSList
CS{end+1}  =crystalSymmetry('1','mineral','something'); % obviously this
can also be defined during import, put here whatever required
ebsd.CSList = CS
ebsd.phaseMap(end+1) = ebsd.phaseMap(end)+1;
%% task: assing ebsd of largest grain to the new symmetry
[grains,ebsd.grainId]=ebsd.calcGrains;

cond = grains.grainSize==max(grains.grainSize);

grains(cond).CS = CS{end};
ebsd(grains(cond)).phase=ebsd.phaseMap(end);

plot(ebsd);


Cheers,
Rüdiger

[Matheus Brozovic Gariglio]

Hi Rüdiger,

Yes, we there are both phases, they are very beautiful!

Yes, it worked, thanks a lot!

Cheers,

Matheus

[Ralf Hielscher]

Actually, with the new MTEX you can even do

ebsd(cond).CS = csNew;

grains(cond).CS = csNew;

Ralf.