Error report

[PSBHS]

Hello Frank and Azdiar,

I face an error when I reconstructed the parent austenite at the step of removing badly reconstructed clusters, as shown below:

job.calcParentFromVote
Error using any
Dimension argument must be a positive integer scalar, a vector of unique positive integers, or 'all'.

Error in subsind (line 61)
sub = any(subs{i}, find(size(subs{i}')==max(size(ind)),1));

Error in grain2d/subsref (line 25)
ind = subsind(grains,s(1).subs);

Error in parentGrainReconstructor/calcParentFromVote (line 75)
pOri = variants(job.p2c, job.grainsPrior(doTransform).meanOrientation, job.votes.parentId(doTransform,1));

The ebsd data contains about 50% retained austenite. Is this the reason for the present error?

[AzdiarGazder]

Hi @PSBHS,

I'm afraid I have very little information to help diagnose the problem.
In cases like this, please provide/upload ALL the background information related to the sample, map and script.
Your detailed background information assists us to understand the problem firsthand and then find out what the possible resolution could be.

Warm regards,
Azdi

[PSBHS]

Hello Azdiar,

I checked the error and think it's related to the function 'calcParentFromVote'. I don't know how to fix it.
The following link contains the raw data and my script. Hope for your answer!
https://connecthkuhk-my.sharepoint.com/:f:/g/personal/pans_connect_hku_hk/EtN26sDTrzFEuRxQ4YCyEHAB5vXSvWkIkfNMHmA4bUNcpA?e=a6PpAU

Best

[AzdiarGazder]

Hi @PSBHS,

Thank you for uploading your map and script.
Indeed, I can confirm there is an error originating within Mtex.
I will be working on a solution with Frank and let you know of a solution when we have one at hand.

Warm regards,
Azdi

[AzdiarGazder]

Hi @PSBHS,

Ralf implemented a fix to MTEX overnight.

Please download and use these latest developer versions of MTEX and ORTools.
I ran your script (using peak finder for the initial OR guess) and found that it processed your map without any errors.
Please do the above on your system and confirm if your issue has been resolved.

Warm regards,
Azdi

[AzdiarGazder]

For reference, here is the script I used:

% Reconstructing the prior austenite microstructure from martensite
% with the variant graph approach and analyzing variant pairing
% The MTEX description for reconstruction is given here:
% https://mtex-toolbox.github.io/MaParentGrainReconstruction.html
% *********************************************************************
home; close all; clear variables;
currentFolder;
set(0,'DefaultFigureWindowStyle','normal');
screenPrint('StartUp','ORTools - HR-900-10min.ang');

%% Initialize MTEX
% Startup and set some settings
startup_mtex;
setMTEXpref('xAxisDirection','east');
setMTEXpref('zAxisDirection','outOfPlane');
setMTEXpref('FontSize',14);   
setInterp2Tex;

% Default directories - Do not modify
Ini.dataPath = [strrep(pwd,'\','/'),'/data/'];
Ini.cifPath = [Ini.dataPath,'input/cif/'];
Ini.ebsdPath = [Ini.dataPath,'input/ebsd/'];
Ini.texturePath = [Ini.dataPath,'output/texture/'];
Ini.imagePath = [Ini.dataPath,'output/images/'];

%% Load data
% Load an MTEX dataset into 'ebsd'
% create an EBSD variable containing the data
ebsd = EBSD.load([Ini.ebsdPath,'HR-900-10min.ang'],'interface','ang','setting 2','convertEuler2SpatialReferenceFrame');


%% Compute, filter and smooth grains
screenPrint('SegmentStart','Computing, filtering and smoothing grains');
% Grains are calculated with a 3° threshold
[grains,ebsd.grainId] = calcGrains(ebsd('indexed'),'angle',3*degree,'unitCell');
% EBSD data in small grains are removed
ebsd(grains(grains.grainSize < 10)) = [];
% Recalculate the grains from the remaining data ...
[grains,ebsd.grainId] = calcGrains(ebsd('indexed'),'angle',3*degree,'unitCell');
% ... and smooth the grain boundaries
grains = smooth(grains,5);
% F = splineFilter;
% ebsd = smooth(ebsd,F,'fill');

%% Rename and recolor phases 
screenPrint('SegmentStart','Renaming and recoloring phases');
phaseNames = {'AlphaP','Gamma'};
% Rename 'Fe (alpha, ferrite, bcc)' to 'AlphaP' and 'Fe (gamma, austenite, fcc)' to 'Gamma'
ebsd = renamePhases(ebsd,phaseNames);
% Choose your favourite colors
ebsd = recolorPhases(ebsd);

%% Define and refine parent-to-child orientation relationship
screenPrint('SegmentStart','Define and refine parent-to-child OR');
% Define 'Gamma' as the parent and 'AlphaP' as the child phase
job = setParentGrainReconstructor(ebsd,grains,Ini.cifPath);
% Find the optimal OR via peak-fitting...
job = defineORs(job);
% ... and refine it based on the fit with boundary misorientations
job.calcParent2Child;
% Display information about the OR
ORinfo(job.p2c);

KS = orientation.KurdjumovSachs(job.csParent,job.csChild);
NW = orientation.NishiyamaWassermann(job.csParent,job.csChild);
plotHist_OR_misfit(job,[KS,NW],'legend',{'K-S OR','N-W OR'});


%% Plotting (with ORTools functions)
screenPrint('SegmentStart','Plotting some ORTools maps');
% Use some of the ORTools functions to visualize the determined OR
% and its relation to the microstructure

% Phase map
plotMap_phases(job,'linewidth',1);
%    - There is no retained austenite (gamma)

% Parent and child IPF maps
plotMap_IPF_p2c(job,vector3d.Z,'linewidth',1);
%    - It is clear that martensite has formed from multiple prior 
%    - austenite grains and that some surface scratches led to bad indexing

% Child-child grain boundary misorientation map
plotMap_gB_c2c(job,'linewidth',1);
%    - Misorientation angles of ~15-50° are not present within prior
%    austenite grains and thus delinitate prior austenite grain
%    boundaries

% Plot a map of the OR boundary disorientation, or misfit
plotMap_gB_misfit(job,'linewidth',1,'maxColor',5);
%    - By setting a threshold at 5 degrees, the prior austenite grain
%    boundaries are identified by their large misfit with the OR

% Plot a map of the OR boundary probability 
plotMap_gB_prob(job,'linewidth',1);
%   - the same can be visualized by calculating the probability that a
%     boundary belongs to the OR

%% Reconstruct parent microstructure
%   - Reconstruct the microstructure with the variant graph based approach
job.calcVariantGraph('threshold',5*degree,'tolerance',2.5*degree,'mergeSimilar')
job.clusterVariantGraph('includeSimilar');
% ... plot the votes (high values show high certainty)
figure; plot(job.grains,job.votes.prob(:,1))
mtexColorbar
% ... and calculate the parent orientations
job.calcParentFromVote('minProb',0)
% Plot the reconstructed parent microstructure
figure;
plot(job.parentGrains,job.parentGrains.meanOrientation,'linewidth',1);

%% Remove badly reconstructed clusters
% In order to reconstruct the remaining parent grains, we can calculate the
% votes for surrounding parent grains by the already reconstructed parent
% grains

% compute the votes
job.calcGBVotes('p2c','reconsiderAll')
% assign parent orientations according to the votes
job.calcParentFromVote
% plot the result
figure;
plot(job.parentGrains,job.parentGrains.meanOrientation)

%% Clean reconstructed grains
% Now clean the grains by: 
% - merging grains with similar orientation
job.mergeSimilar('threshold',7.5*degree);
% - merging small inclusions
job.mergeInclusions('maxSize',250);
% This is the cleaned reconstructed parent microstructure
figure;
plot(job.parentGrains,job.parentGrains.meanOrientation,'linewidth',1)

%% Variant analysis
% Calculate the variant IDs of all alpha grains
job.calcVariants;
% ... and plot them
plotMap_variants(job,'grains','linewidth',1); %Grain data
plotMap_variants(job,'linewidth',1); %EBSD data

% We see that each child grain contains several variants on EBSD level
% For analyzing variant pairing, we need the variants on the EBSD level
% to be reconstructed as grains
[newGrains,newEBSD] = computeVariantGrains(job);
%Compare variant indexing for old and new grains
figure; %Old grains
plot(grains(job.csChild),grains(job.csChild).meanOrientation);
figure; %New grains
plot(newGrains(job.csChild),newGrains(job.csChild).meanOrientation);

return
%% Display variant information for a PAG by interactively clicking on a gamma grain
grainClick(job,'noScalebar','noFrame'); %Based on EBSD data

%% Do the same at the grain level 
grainClick(job,'grains','noFrame'); %Based on grain data

%% Refine gamma twins for a PAG by interactively clicking on a gamma grain
grainClick(job,'parentTwins');

[PSBHS]

Hello Azdiar,

I have run the script and no error was reported. Thank you!

Best