Giter VIP home page Giter VIP logo

ortools4mtex / ortools Goto Github PK

View Code? Open in Web Editor NEW
32.0 5.0 7.0 111.21 MB

A function library for OR discovery, OR analysis and the plotting of publication-ready figures of martensitic transformations or OR-related phase transitions in MTEX.

Home Page: https://github.com/ORTools4MTEX/ORTools

License: MIT License

MATLAB 99.95% Shell 0.05%
martensite austenite ebsd parent child grain reconstruction texture crystal orientation-relationship

ortools's People

Contributors

azdiargazder avatar frankniessen avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar

ortools's Issues

variant analysis and doubt

Discussed in #17

Hi, frank,
Thank you for sharing the code of variant analysis in lath martensite. One of the problems that has puzzled me is that if I draw the ratio of boundary length of V1-V2... V1-V24 variant pairs and the score of variant pairs that follow the K-S relationship.
Only V1-V2 is given in plotMap_variantPairs. m in ORtools; V1-V3 (5); V1-V4; V1-V6 information is not given for other variant pairs.
I wrote a part of it myself, but it didn't work out the same way as other organizations, so I'd like you to look at it or better to provide the relevant code. I would appreciate it if you could provide the variantpairs boundary length density and the variantpairs fraction code.

%% Determine the variant pairs
% Derive variant types 1 to 6 from variants 1 to 24
% variantGrains.prop.variantType = variantGrains.variantId - (variantGrains.packetId-1) * 24/4;
variantGrains.prop.variantType = variantGrains.variantId;
% Get the boundary segments and neighbouring variants
variantBoundary = variantGrains.boundary(variantGrains.CS,variantGrains.CS);
variantBoundaryIds = variantBoundary.grainId;

[c,~] = ismember(variantBoundaryIds,variantGrains.id);
variantBoundaryIds(any(~c,2),:) = [];
variantBoundary(any(~c,2)) = [];

%Identify any special boundaries, here V1-V2 : V1:V6
varTypes = variantGrains(variantGrains.id2ind(variantBoundaryIds)).variantType;
% cond(1,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,2),2) | ...
% any(ismember(varTypes,3),2) & any(ismember(varTypes,4),2) | ...
% any(ismember(varTypes,5),2) & any(ismember(varTypes,6),2); %V1-V2
% cond(2,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,3),2) | ...
% any(ismember(varTypes,1),2) & any(ismember(varTypes,5),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,4),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,6),2) | ...
% any(ismember(varTypes,3),2) & any(ismember(varTypes,5),2) | ...
% any(ismember(varTypes,4),2) & any(ismember(varTypes,6),2); %V1-V3(V5)
% cond(3,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,6),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,3),2) | ...
% any(ismember(varTypes,4),2) & any(ismember(varTypes,5),2); %V1-V6
% cond(4,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,4),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,5),2) | ...
% any(ismember(varTypes,3),2) & any(ismember(varTypes,6),2); %V1-V4

%-------------------------------------------------------------------
vec1=any(ismember(varTypes,1),2);
vec2=any(ismember(varTypes,2),2);
vec3=any(ismember(varTypes,3),2);
vec4=any(ismember(varTypes,4),2);
vec5=any(ismember(varTypes,5),2);
vec6=any(ismember(varTypes,6),2);
vec7=any(ismember(varTypes,7),2);
vec8=any(ismember(varTypes,8),2);
vec9=any(ismember(varTypes,9),2);
vec10=any(ismember(varTypes,10),2);
vec11=any(ismember(varTypes,11),2);
vec12=any(ismember(varTypes,12),2);
vec13=any(ismember(varTypes,13),2);
vec14=any(ismember(varTypes,14),2);
vec15=any(ismember(varTypes,15),2);
vec16=any(ismember(varTypes,16),2);
vec17=any(ismember(varTypes,17),2);
vec18=any(ismember(varTypes,18),2);
vec19=any(ismember(varTypes,19),2);
vec20=any(ismember(varTypes,20),2);
vec21=any(ismember(varTypes,21),2);
vec22=any(ismember(varTypes,22),2);
vec23=any(ismember(varTypes,23),2);
vec24=any(ismember(varTypes,24),2);
%-------------------------------------------------------------
%======================================
cond(1,:) = vec1&vec2|vec3&vec4|vec5&vec6|vec7&vec8|vec9&vec10|vec11&vec12|vec13&vec14|vec15&vec16|vec17&vec18|vec19&vec20|vec21&vec22|vec23&vec24;

cond(2,:) = vec1&vec3|vec2&vec4|vec3&vec5|vec4&vec6|vec5&vec7|vec6&vec8|vec7&vec9|vec8&vec10|vec9&vec11|vec10&vec12|vec11&vec13|vec12&vec14|vec13&vec15|vec14&vec16|vec15&vec17|vec16&vec18|vec17&vec19|vec18&vec20|vec19&vec21|vec20&vec22|vec21&vec23|vec22&vec24|vec1&vec5|vec2&vec6|vec3&vec7|vec4&vec8|vec5&vec9|vec6&vec10|vec7&vec11|vec8&vec12|vec9&vec13|vec10&vec14|vec11&vec15|vec12&vec16|vec13&vec17|vec14&vec18|vec15&vec19|vec16&vec20|vec17&vec21|vec18&vec22|vec19&vec23|vec20&vec24;

cond(3,:) = vec1&vec4|vec2&vec5|vec3&vec6|vec4&vec7|vec5&vec8|vec6&vec9|vec7&vec10|vec8&vec11|vec9&vec12|vec10&vec13|vec11&vec14|vec12&vec15|vec13&vec16|vec14&vec17|vec15&vec18|vec16&vec19|vec17&vec20|vec18&vec21|vec19&vec22|vec20&vec23|vec21&vec24;

cond(4,:) = vec1&vec6|vec2&vec7|vec3&vec8|vec4&vec9|vec5&vec10|vec6&vec11|vec7&vec12|vec8&vec13|vec9&vec14|vec10&vec15|vec11&vec16|vec12&vec17|vec13&vec18|vec14&vec19|vec15&vec20|vec16&vec21|vec17&vec22|vec18&vec23|vec19&vec24;

%======================================

cond(5,:) = vec1&vec7|vec2&vec8|vec3&vec9|vec4&vec10|vec5&vec11|vec6&vec12|vec7&vec13|vec8&vec14|vec9&vec15|vec10&vec16|vec11&vec17|vec12&vec18|vec13&vec19|vec14&vec20|vec15&vec21|vec16&vec22|vec17&vec23|vec18&vec24;

cond(6,:) = vec1&vec8|vec2&vec9|vec3&vec10|vec4&vec11|vec5&vec12|vec6&vec13|vec7&vec14|vec8&vec15|vec9&vec16|vec10&vec17|vec11&vec18|vec12&vec19|vec13&vec20|vec14&vec21|vec15&vec22|vec16&vec23|vec17&vec24;

cond(7,:) = vec1&vec9|vec2&vec10|vec3&vec11|vec4&vec12|vec5&vec13|vec6&vec14|vec7&vec15|vec8&vec16|vec9&vec17|vec10&vec18|vec11&vec19|vec12&vec20|vec13&vec21|vec14&vec22|vec15&vec23|vec16&vec24|vec1&vec19|vec2&vec20|vec3&vec21|vec4&vec22|vec5&vec23|vec6&vec24;

cond(8,:) = vec1&vec10|vec2&vec11|vec3&vec12|vec4&vec13|vec5&vec14|vec6&vec15|vec7&vec16|vec8&vec17|vec9&vec18|vec10&vec19|vec11&vec20|vec12&vec21|vec13&vec22|vec14&vec23|vec15&vec24|vec1&vec14|vec2&vec15|vec3&vec16|vec4&vec17|vec5&vec18|vec6&vec19|vec7&vec20|vec8&vec21|vec9&vec22|vec10&vec23|vec11&vec24;

%========================================

cond(9,:) = vec1&vec11|vec2&vec12|vec3&vec13|vec4&vec14|vec5&vec15|vec6&vec16|vec7&vec17|vec8&vec18|vec9&vec19|vec10&vec20|vec11&vec21|vec12&vec22|vec13&vec23|vec14&vec24|vec1&vec13|vec2&vec14|vec3&vec15|vec4&vec16|vec5&vec17|vec6&vec18|vec7&vec19|vec8&vec20|vec9&vec21|vec10&vec22|vec11&vec23|vec12&vec24;

cond(10,:) = vec1&vec12|vec2&vec13|vec3&vec14|vec4&vec15|vec5&vec16|vec6&vec17|vec7&vec18|vec8&vec19|vec9&vec20|vec10&vec21|vec11&vec22|vec12&vec23|vec13&vec24|vec1&vec20|vec2&vec21|vec3&vec22|vec4&vec23|vec5&vec24;

cond(11,:) = vec1&vec15|vec2&vec16|vec3&vec17|vec4&vec18|vec5&vec19|vec6&vec20|vec7&vec21|vec8&vec22|vec9&vec23|vec10&vec24|vec1&vec23|vec2&vec24;

cond(12,:) = vec1&vec16|vec2&vec17|vec3&vec18|vec4&vec19|vec5&vec20|vec6&vec21|vec7&vec22|vec8&vec23|vec9&vec24;

%=========================================

cond(13,:) = vec1&vec17|vec2&vec18|vec3&vec19|vec4&vec20|vec5&vec21|vec6&vec22|vec7&vec23|vec8&vec24;

cond(14,:) = vec1&vec18|vec2&vec19|vec3&vec20|vec4&vec21|vec5&vec22|vec6&vec23|vec7&vec24|vec1&vec22|vec2&vec23|vec3&vec24;

cond(15,:) = vec1&vec21|vec2&vec22|vec3&vec23|vec4&vec24;

cond(16,:) = vec1&vec24|vec2&vec3|vec4&vec5|vec6&vec7|vec8&vec9|vec10&vec11|vec12&vec13|vec14&vec15|vec16&vec17|vec18&vec19|vec20&vec21|vec22&vec23;

for ii = 1:size(cond,1)
variantPairs_boundary{ii} = variantBoundary(cond(ii,:));
end

Variant pairing question

Discussed in #26

Hi, my name is Long Wang. Thank you very much for your reply.
Because I am not familiar with variant pairing, I am based on your code, do not know whether it is correct or not. I don't understand the definition of variant pairs in your code(plotMap_variantPairs.m,ORTools(src)), so can you please explain the derivation of 16 variant pairs(doi:10.1016/j.actamat.2011.12.018)?The density and fraction of boundary length of martensitic variant pairs,sunch as V1-V2; V1-V3 (V5); V1-V6; V1-V4;V1-V7;V1-V8;V1-V9(V19) and so on; Can all the remaining variant pairs be supplemented? Sorry to bother you, but I have launched the code myself, but I don't know if it is right. Can you have a look?
image
%Identify any special boundaries, here V1-V2 : V1:V6
varTypes = variantGrains(variantGrains.id2ind(variantBoundaryIds)).variantType;
% cond(1,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,2),2) | ...
% any(ismember(varTypes,3),2) & any(ismember(varTypes,4),2) | ...
% any(ismember(varTypes,5),2) & any(ismember(varTypes,6),2); %V1-V2
% cond(2,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,3),2) | ...
% any(ismember(varTypes,1),2) & any(ismember(varTypes,5),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,4),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,6),2) | ...
% any(ismember(varTypes,3),2) & any(ismember(varTypes,5),2) | ...
% any(ismember(varTypes,4),2) & any(ismember(varTypes,6),2); %V1-V3(V5)
% cond(3,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,6),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,3),2) | ...
% any(ismember(varTypes,4),2) & any(ismember(varTypes,5),2); %V1-V6
% cond(4,:) = any(ismember(varTypes,1),2) & any(ismember(varTypes,4),2) | ...
% any(ismember(varTypes,2),2) & any(ismember(varTypes,5),2) | ...
% any(ismember(varTypes,3),2) & any(ismember(varTypes,6),2); %V1-V4

%-------------------------------------------------------------------
vec1=any(ismember(varTypes,1),2);
vec2=any(ismember(varTypes,2),2);
vec3=any(ismember(varTypes,3),2);
vec4=any(ismember(varTypes,4),2);
vec5=any(ismember(varTypes,5),2);
vec6=any(ismember(varTypes,6),2);
vec7=any(ismember(varTypes,7),2);
vec8=any(ismember(varTypes,8),2);
vec9=any(ismember(varTypes,9),2);
vec10=any(ismember(varTypes,10),2);
vec11=any(ismember(varTypes,11),2);
vec12=any(ismember(varTypes,12),2);
vec13=any(ismember(varTypes,13),2);
vec14=any(ismember(varTypes,14),2);
vec15=any(ismember(varTypes,15),2);
vec16=any(ismember(varTypes,16),2);
vec17=any(ismember(varTypes,17),2);
vec18=any(ismember(varTypes,18),2);
vec19=any(ismember(varTypes,19),2);
vec20=any(ismember(varTypes,20),2);
vec21=any(ismember(varTypes,21),2);
vec22=any(ismember(varTypes,22),2);
vec23=any(ismember(varTypes,23),2);
vec24=any(ismember(varTypes,24),2);
%-------------------------------------------------------------
%======================================
cond(1,:) = vec1&vec2|vec3&vec4|vec5&vec6|vec7&vec8|vec9&vec10|vec11&vec12|vec13&vec14|vec15&vec16|vec17&vec18|vec19&vec20|vec21&vec22|vec23&vec24;

cond(2,:) = vec1&vec3|vec2&vec4|vec3&vec5|vec4&vec6|vec5&vec7|vec6&vec8|vec7&vec9|vec8&vec10|vec9&vec11|vec10&vec12|vec11&vec13|vec12&vec14|vec13&vec15|vec14&vec16|vec15&vec17|vec16&vec18|vec17&vec19|vec18&vec20|vec19&vec21|vec20&vec22|vec21&vec23|vec22&vec24|vec1&vec5|vec2&vec6|vec3&vec7|vec4&vec8|vec5&vec9|vec6&vec10|vec7&vec11|vec8&vec12|vec9&vec13|vec10&vec14|vec11&vec15|vec12&vec16|vec13&vec17|vec14&vec18|vec15&vec19|vec16&vec20|vec17&vec21|vec18&vec22|vec19&vec23|vec20&vec24;

cond(3,:) = vec1&vec4|vec2&vec5|vec3&vec6|vec4&vec7|vec5&vec8|vec6&vec9|vec7&vec10|vec8&vec11|vec9&vec12|vec10&vec13|vec11&vec14|vec12&vec15|vec13&vec16|vec14&vec17|vec15&vec18|vec16&vec19|vec17&vec20|vec18&vec21|vec19&vec22|vec20&vec23|vec21&vec24;

cond(4,:) = vec1&vec6|vec2&vec7|vec3&vec8|vec4&vec9|vec5&vec10|vec6&vec11|vec7&vec12|vec8&vec13|vec9&vec14|vec10&vec15|vec11&vec16|vec12&vec17|vec13&vec18|vec14&vec19|vec15&vec20|vec16&vec21|vec17&vec22|vec18&vec23|vec19&vec24;

%======================================

cond(5,:) = vec1&vec7|vec2&vec8|vec3&vec9|vec4&vec10|vec5&vec11|vec6&vec12|vec7&vec13|vec8&vec14|vec9&vec15|vec10&vec16|vec11&vec17|vec12&vec18|vec13&vec19|vec14&vec20|vec15&vec21|vec16&vec22|vec17&vec23|vec18&vec24;

cond(6,:) = vec1&vec8|vec2&vec9|vec3&vec10|vec4&vec11|vec5&vec12|vec6&vec13|vec7&vec14|vec8&vec15|vec9&vec16|vec10&vec17|vec11&vec18|vec12&vec19|vec13&vec20|vec14&vec21|vec15&vec22|vec16&vec23|vec17&vec24;

cond(7,:) = vec1&vec9|vec2&vec10|vec3&vec11|vec4&vec12|vec5&vec13|vec6&vec14|vec7&vec15|vec8&vec16|vec9&vec17|vec10&vec18|vec11&vec19|vec12&vec20|vec13&vec21|vec14&vec22|vec15&vec23|vec16&vec24|vec1&vec19|vec2&vec20|vec3&vec21|vec4&vec22|vec5&vec23|vec6&vec24;

cond(8,:) = vec1&vec10|vec2&vec11|vec3&vec12|vec4&vec13|vec5&vec14|vec6&vec15|vec7&vec16|vec8&vec17|vec9&vec18|vec10&vec19|vec11&vec20|vec12&vec21|vec13&vec22|vec14&vec23|vec15&vec24|vec1&vec14|vec2&vec15|vec3&vec16|vec4&vec17|vec5&vec18|vec6&vec19|vec7&vec20|vec8&vec21|vec9&vec22|vec10&vec23|vec11&vec24;

%========================================

cond(9,:) = vec1&vec11|vec2&vec12|vec3&vec13|vec4&vec14|vec5&vec15|vec6&vec16|vec7&vec17|vec8&vec18|vec9&vec19|vec10&vec20|vec11&vec21|vec12&vec22|vec13&vec23|vec14&vec24|vec1&vec13|vec2&vec14|vec3&vec15|vec4&vec16|vec5&vec17|vec6&vec18|vec7&vec19|vec8&vec20|vec9&vec21|vec10&vec22|vec11&vec23|vec12&vec24;

cond(10,:) = vec1&vec12|vec2&vec13|vec3&vec14|vec4&vec15|vec5&vec16|vec6&vec17|vec7&vec18|vec8&vec19|vec9&vec20|vec10&vec21|vec11&vec22|vec12&vec23|vec13&vec24|vec1&vec20|vec2&vec21|vec3&vec22|vec4&vec23|vec5&vec24;

cond(11,:) = vec1&vec15|vec2&vec16|vec3&vec17|vec4&vec18|vec5&vec19|vec6&vec20|vec7&vec21|vec8&vec22|vec9&vec23|vec10&vec24|vec1&vec23|vec2&vec24;

cond(12,:) = vec1&vec16|vec2&vec17|vec3&vec18|vec4&vec19|vec5&vec20|vec6&vec21|vec7&vec22|vec8&vec23|vec9&vec24;

%=========================================

cond(13,:) = vec1&vec17|vec2&vec18|vec3&vec19|vec4&vec20|vec5&vec21|vec6&vec22|vec7&vec23|vec8&vec24;

cond(14,:) = vec1&vec18|vec2&vec19|vec3&vec20|vec4&vec21|vec5&vec22|vec6&vec23|vec7&vec24|vec1&vec22|vec2&vec23|vec3&vec24;

cond(15,:) = vec1&vec21|vec2&vec22|vec3&vec23|vec4&vec24;

cond(16,:) = vec1&vec24|vec2&vec3|vec4&vec5|vec6&vec7|vec8&vec9|vec10&vec11|vec12&vec13|vec14&vec15|vec16&vec17|vec18&vec19|vec20&vec21|vec22&vec23;

Problems in Parent Austenite Reconstruction

Dear Frank,
Hi
I have some questions. I really appreciate if you could help me.

  • If i want to plot only specific variants (for example only V1 and V2) in my orientation map, Could you please guide me what command i need to use ?
  • Also, If i want to show that which parent austenite grain boundaries are not KS (misorientation from theoretical ones is more than 5 degree), could you please help me with that too?
  • Last question, after reconstruction how can i export the orientation (Euler angles) of parent austenite grains as an EBSD data?
    Thank you for your time.

How to export .cpr file after phase reconstruction

Hello.
I want to export .cpr file to re-analysis in Aztec program.
I tried .ctf file format and it was possible by this job.ebsd.export('filename.ctf')
However, there is a problem when I export .cpr format
Actually, I can not understand meaning of the Pfname in this Syntax export_crc(ebsd,Pfname)
I am waiting your reply and have a nice day.

Thanks,

ORTools-Example 2-Parent Grain Reconstruction in Titanium alloys error

Hello.

I have trouble when I practice with ORTools-Example2-Paraent grain reconstruction in titanium alloys

Actually, the it is the same problem with Naga Chaitanya in comment.
I am also when I doing peak fitting(3:36 of video)

Undefined function 'findpeaks' for input arguments of type
'double'.

Error in gaussFit (line 46)
[peakValue,peakPositionIdx,peakWidth,peakHeight] = findpeaks(y,
'MinPeakDist',classInterval.fine,
'MinPeakHeight',ypos,'Annotate','extents','WidthReference','halfheight');

Error in defineORs (line 45)
[misoRange] =
gaussFit(classRange,classInterval,counts);

Error in ORTools_example2 (line 54)
job = defineORs(job);

this error is occurred.
How can I solve this problem?

thank you and have a nice day

defineORs(job)

Dear Frank,

I am trying to use this example on my own EBSD analysis however, there is no Beta in my materials which leads the defineORs(job) to go to isempty condition. So I get this window asking me to define an OR using Miller-Bravais indices:
image
What would be the proper values?
Regards,
Mohsen

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.