cs = crystalSymmetry('6/mmm', [3.2093 3.2093 5.2103], 'X||a*', 'Y||b', 'Z||c*', 'mineral', 'Mg', 'color', 'light blue');
data=[
0 | 0 | 2 | 1.66917 |
1 | 0 | 5 | 0.171019 |
1 | 0 | 3 | 0.896366 |
1 | 0 | 2 | 0.659014 |
1 | 0 | 1 | 3.712668 |
1 | 0 | 0 | 0.584565 |
1 | 1 | 0 | 0.348282 |
1 | 1 | 2 | 0.324243 |
0 | 2 | 1 | 0.194963 |
0 | 2 | 3 | 0.109812 |
1 | 2 | 1 | 0.079389 |
1 | 1 | 4 | 0.084234 |
];
v = Miller(vector3d(data(:,1:3).'),cs)
plot(v.project2FundamentalRegion,data(:,4),'FundamentalRegion')
S2G = plotS2Grid(cs.fundamentalSector)
S2G=Miller(S2G,cs);
v.resolution = 10*degree % this is used as an parameter of interpolation
d = interp(v,data(:,4),S2G);
plotx2south
contourf(S2G,d)