mtex-toolbox / mtex

MTEX is a free Matlab toolbox for quantitative texture analysis. Homepage:
http://mtex-toolbox.github.io/
GNU General Public License v2.0
281 stars 185 forks source link

triple junctions on merged grains #136

Closed kilir closed 8 years ago

kilir commented 8 years ago

Hi Ralf, triple junction on grains work fine, however, when I merge these grains to get rid of twin boundaries, triple junctions on merged grains still returns junctions that include twin boundaries.

% calcgrains
[grainsl ebsdc.grainId] = calcGrains(ebsdc,'angle',10*degree);

% get tp points
tp = grainsl('q').triplePoints('q','q','q');

figure
plot(grainsl('q'))
hold on
plot(tp,'MarkerSize',8,'MarkerFaceColor','k')
hold off
% looks good 

tp_non_merged

% but we don't want to include twin boundaries  - merge dauphine
gb_qtz = grainsl.boundary('Quartz','Quartz');
rot1 = rotation('axis',Miller(0,0,0,1,CS{2}),'angle',60*degree);
ind1 = angle(gb_qtz('Quartz','Quartz').misorientation,rot1)<6*degree; 
twinBoundary = gb_qtz(ind1);
[grains_mergedc,parentID_grains]= merge(grainsl,twinBoundary);
ebsd_mergedc = ebsdc;
ebsd_mergedc.grainId  = parentID_grains(ebsdc.grainId);

% get tp points
clear tp
tp = grains_mergedc('q').triplePoints('q','q','q');

% there are triple junction along boundaries
figure
plot(grains_mergedc('q'))
hold on
plot(tp,'MarkerSize',8,'MarkerFaceColor','k')
hold off

% seems that triple junctions still relay on the original grain boundary?
figure
plot(grains_mergedc('q'))
hold on
plot(twinBoundary,'lineColor','r','linewidth',2)
hold on
plot(tp,'MarkerSize',8,'MarkerFaceColor','k')
hold off

There are triple junctions indicated where none should be. tp_merged

It looks like the old twin boundaries are still considered, despite grains are merged. tp_merged_dau

Did I miss something or do I have an error in the script above? Or do I have to implicitly exclude/delete triple junctions that belong to a twin boundary?

All the best, Rüdiger

ralfHielscher commented 8 years ago

Hi Rüdiger,

thank you for reporting this. It will be fixed in the next minor release.

Ralf.