interp/code/datfile2em.m at master · sgbaird-5DOF/interp · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
function [TJs,e1,e2,e3,m1,m2,m3,nTJs] = datfile2em(fpath,nheaderlines,deg2radQ,invQ,normrQ)
arguments
    fpath char = '../../TJ2GBE/TJdata/triples_30000.dat'
    nheaderlines(1,1) double = 0
    deg2radQ(1,1) logical = true
    invQ(1,1) logical = true
    normrQ(1,1) logical = true
end
% READ_DAT  extract TJs, EAs, and norms from triples.dat file.
% Mirrors https://github.com/Yufeng-shen/TJ2GBE/blob/master/Src/Python/Reconstruction.py
%--------------------------------------------------------------------------
% Author(s): Sterling Baird
%
% Date: 2020-12-05
%
%     Input: triples.dat, wrote from the fortran program Torq_gen
%                 size=[numTJ*8,]
%                 In each group, the data is [TJ directon, EA1, GB1, EA2, GB2, EA3, GB3]
%     Output: TJs, direction of the triple junctions
%                 size = [numTJ, 3]
%             EAs, the EA angles of the 3 grains at a TJ
%                 size = [numTJ, 3, 3]
%             norms, normal direction of the 3 GB at a TJ
%                 size = [numTJ, 3, 3]
%
% .dat file format (everything in sample reference frame)
% 1 - TJnum         (lines 1, 8, etc.)
%   2 3 4 - TJs     (lines 2, 9, etc.)
%   5 6 7 - e1      (lines 3, 10, etc.)
%   8 9 10 - m1     (lines 4, 11, etc.)
%   11 12 13 - e2   (lines 5, 12, etc.)
%   14 15 16 - m2   (lines 6, 13, etc.)
%   17 18 19 - e3   (lines 7, 14, etc.)
%   20 21 22 - m3   (lines 8, 15, etc.)
%--------------------------------------------------------------------------
data = importdata(fpath,' ',nheaderlines);
TJs = [data(2:22:end),data(3:22:end),data(4:22:end)]; %lines 2, 9, etc.

e1 = [data(5:22:end),data(6:22:end),data(7:22:end)]; %lines 3, 10, etc.
m1 = [data(8:22:end),data(9:22:end),data(10:22:end)]; %lines 4, 11, etc.

e2 = [data(11:22:end),data(12:22:end),data(13:22:end)]; %lines 5, 12, etc.
m2 = [data(14:22:end),data(15:22:end),data(16:22:end)]; %lines 6, 13, etc.

e3 = [data(17:22:end),data(18:22:end),data(19:22:end)]; %lines 7, 14, etc.
m3 = [data(20:22:end),data(21:22:end),data(22:22:end)]; %lines 8, 15, etc.

% convert Euler angles from degrees to radians
if deg2radQ
    e1 = deg2rad(e1);
    e2 = deg2rad(e2);
    e3 = deg2rad(e3);
end

if invQ
    q1 = eu2qu(e1,-1);
    q2 = eu2qu(e2,-1);
    q3 = eu2qu(e3,-1);

    e1 = qu2eu(q1,1);
    e2 = qu2eu(q2,1);
    e3 = qu2eu(q3,1);
end

if normrQ
    % renormalize the BP normals
    m1 = normr(m1);
    m2 = normr(m2);
    m3 = normr(m3);
end

% no output, but useful as a double check:
TJnum = data(1:22:end); %lines 1, 8, etc.
nTJs = length(TJnum); %#ok<NASGU>

end
%% CODE GRAVEYARD
%{
% nTJs = length(TJnum);
% TJ = zeros(nTJs,3);

% EAs = [...
%     data(5:22:end),data(6:22:end),data(7:22:end); %lines 3, 10, etc.
%     data(8:22:end),data(9:22:end),data(10:22:end); %lines 4, 11, etc.
%     data(11:22:end),data(12:22:end),data(13:22:end)]; %lines 5, 12, etc.
% norms = [...
%     data(14:22:end),data(15:22:end),data(16:22:end); %lines 6, 13, etc.
%     data(17:22:end),data(18:22:end),data(19:22:end); %lines 7, 14, etc.
%     data(20:22:end),data(21:22:end),data(22:22:end)]; %lines 8, 15, etc.
%}