obtaining subgrains Id from specific grains

[louisMabon]

Hi,
I am currently working on forged materials. During my analysis, I need to separate 2 populations from one phase. I already serarated my grain in the alpha phase but I don't know how to obtain their subgrain ID.

I have here separated the alpha phase into the blue and red population.

I would like to obtain the subgrain ID from the red grains. And if possible to convert some subgrains into blue grains.

here is my code so far:

% define crystal symmetry variables
cs_A = ebsd2('Alpha').CS;
cs_B = ebsd2('Beta').CS;
%_____________________________________________________________________________
%TEST PASSAGE GRAINS ET NON SUBGRAINS 

[grains,ebsd2.grainId] = calcGrains(ebsd2,'angle',HAGB);
grains = smooth(grains,5);
G_a = grains('Alpha');

condi = G_a.aspectRatio > 2 &  G_a.grainSize < 100;
Alpha_sec_grain = G_a(condi);

Alpha_sec_grain_id = Alpha_sec_grain.id;

Alpha_Sec = grains(Alpha_sec_grain_id);
Alpha_Sec_true_Id = vertcat(to_add,Alpha_sec_grain_id);

Alpha_Prim2_Id = setdiff(G_a.id,Alpha_Sec_true_Id);
Alpha_Prim2 = grains(Alpha_Prim2_Id);
Alpha_Prim2 = Alpha_Prim2('Alpha');

Alpha_Sec3 = grains(Alpha_Sec_true_Id);
Alpha_Sec3 = Alpha_Sec3('Alpha');

figure; hold on 
plot(Alpha_Prim2,'Facecolor','r');hold on 
plot (Alpha_Sec3,'Facecolor','b'); 
hold on 

linewidth = 0.75;
plot(subgrains_A.boundary('Alpha','Alpha'),'displayName','\theta > 1°','linewidth',linewidth)
plot(grains_A.boundary('Alpha','Alpha'),'displayName','\theta > 10°','linecolor','yellow','linewidth',linewidth)
[subgrains,ebsd2.grainId] = calcGrains(ebsd2,'angle',LAGB);
subgrains = smooth(subgrains,5);

subgrains_Alpha_Sec3 = subgrains(Alpha_Sec3);
subgrains_Alpha_Sec3 = subgrains_Alpha_Sec3('Alpha');
plot (subgrains_Alpha_Sec3)

§§§§ but when i call subgrains it doesn't work

[subgrains,ebsd2.grainId] = calcGrains(ebsd2,'angle',LAGB);
subgrains = smooth(subgrains,5);


subgrains_Alpha_Sec3 = subgrains(Alpha_Sec3);
subgrains_Alpha_Sec3 = subgrains_Alpha_Sec3('Alpha');
plot (subgrains_Alpha_Sec3)

thanks for your help

[ralfHielscher]

Could you try to explain a little bit more in detail what you want to achieve? I do not understand what is a subgrain in your context and how you compute it.

Ralf.

[louisMabon]

Yes sure.
I am working on forged titanium. During forging , we will deform the alpha phase into a globular shape ( red phase). After forging we will cool the material, leading us to the apparition of needle like alpha phase alpso named seconday alpha phase (blue phase).
On my case , there are some of the secondary alpha phase who have appeared on the globular alpha and are lowly desoriented with each other.

My objective is to pick manually those secondary alpha phase that were asigned as subgrains of the globular one (red phase) and switch them into the blue phase.

[ralfHielscher]

I try to understand your code.

Alpha_sec_grain = G_a(condi);
Alpha_sec_grain_id = Alpha_sec_grain.id;
Alpha_Sec = grains(Alpha_sec_grain_id);

What is the sense of these three lines of code. Will not Alpha_sec_grain be the same as Alpha_Sec ?
Ralf.

[louisMabon]

Hi , yes it is the case , I want to separate the grains stored in the Alpha_sec_XXXX variables from the ones stored in the Alpha_prim_XXX variables. But my code isn't very clear. Il will change the variables anme for it to be easier to understand.

[louisMabon]

% define crystal symmetry variables
cs_A = ebsd2('Alpha').CS;
cs_B = ebsd2('Beta').CS;
%_____________________________________________________________________________
%TEST PASSAGE GRAINS ET NON SUBGRAINS 

[grains,ebsd2.grainId] = calcGrains(ebsd2,'angle',HAGB);
grains = smooth(grains,5);
G_a = grains('Alpha');

condi = G_a.aspectRatio > 2 &  G_a.grainSize < 100;
condi_2 =  G_a.grainSize < 50;

Alpha_sec_grain_2= G_a(condi_2);
Alpha_sec_grain = G_a(condi);

Alpha_sec_grain_id = Alpha_sec_grain.id;
Alpha_sec_grain_2_id = Alpha_sec_grain_2.id;

% assemble all the secondary alpha phases respecting the condition condi and condi_2 in a variable Alpha_sec
Alpha_sec_grain_id = vertcat(Alpha_sec_grain_id,Alpha_sec_grain_2_id);
Alpha_Sec = grains(Alpha_sec_grain_id);

% primary alpha Alpha_Prim are all the alpha grains that does not respect the conditions
Alpha_Prim_Id = setdiff(G_a.id,Alpha_Sec_true_Id);
Alpha_Prim = grains(Alpha_Prim2_Id);
Alpha_Prim = Alpha_Prim('Alpha');

[subgrains,ebsd2.grainId] = calcGrains(ebsd2,'angle',LAGB);
subgrains = smooth(subgrains,5);

% some of the Alpha prim are touching Alpha sec and are lowly desoriented one another, Need to convert some of the %Alpha prim subgrains into Alpha sec subgrains 

subgrains_Alpha_Sec = subgrains(Alpha_Sec);
subgrains_Alpha_Prim = subgrains(Alpha_Prim);

figure; hold on 
plot(subgrains_Alpha_Sec,'Facecolor','r');hold on 
plot (subgrains_Alpha_Prim,'Facecolor','b'); 

[louisMabon]

but when I plot this , it's as if it doesn't take into account the attributions of the grains between Alpha_Prim and Alpha_sec

[ralfHielscher]

I think I finally figured out what you want to achieve. I cleaned up your code a bit. The problem is that the last two lines are not MTEX syntax. Those you need to replace by

subgrains_Alpha_Sec = subgrains( any( Alpha_Sec.checkInside(subgrains.centroid),2))

Your clean code ...

% define crystal symmetry variables
cs_A = ebsd2('Alpha').CS;
cs_B = ebsd2('Beta').CS;

[grains,ebsd2.grainId] = calcGrains(ebsd2,'angle',HAGB);
grains = smooth(grains,5);
G_a = grains('Alpha');

condi = G_a.aspectRatio > 2 &  G_a.grainSize < 100;
condi_2 =  G_a.grainSize < 50;

Alpha_sec_grain_2= G_a(condi_2);
Alpha_sec_grain = G_a(condi);

% assemble all the secondary alpha phases respecting the condition condi or condi_2 in a variable Alpha_sec
Alpha_Sec = G_a(condi | condi2);

% primary alpha Alpha_Prim are all the alpha grains that does not respect the conditions
Alpha_Prim = G_a(~condi & ~condi_2);

[subgrains,ebsd2.grainId] = calcGrains(ebsd2,'angle',LAGB);
subgrains = smooth(subgrains,5);

% some of the Alpha prim are touching Alpha sec and are lowly desoriented one another, Need to convert some of the %Alpha prim subgrains into Alpha sec subgrains 
subgrains_Alpha_Sec = subgrains(Alpha_Sec);       % this is not MTEX syntax
subgrains_Alpha_Prim = subgrains(Alpha_Prim);  % this is not MTEX syntax

[louisMabon]

Yes it worked , thank you verry much !