JQFonseca / met_CDT_EBSD

A programme for plotting orientations from EBSD
3 stars 3 forks source link

EBSD Progression Discussion #5

Open Pstavroulakis opened 5 years ago

Pstavroulakis commented 5 years ago

Hello everyone. Maybe we could use this as a place where we could discuss some issues we have and need resolution.

Pstavroulakis commented 5 years ago

In Marks part:

Mark

XYco=[] # Creates an empty list to hold XY co-ordinates XY=[0,0] # Creates array to put values into within the list

for q in range(len(coordinate)): Co=coordinate[q] # Indexes array within list Co = Co[0] # Indexes value withing array if Co[2]>=0: XY[0]=Co[0]/(Co[2]+1) # Calculates the x and y co-ordinates for the stereographic plot and creates a matrix to hold them (XY). XY[1]=Co[1]/(Co[2]+1) XYco.append(XY)

The line: Co = Co[0] # Indexes value withing array, will produce an error in the if section. I propose that the code is changed with:

Mark

XYco=[] # Creates an empty list to hold XY co-ordinates XY=[0,0] # Creates array to put values into within the list

for q in range(len(coordinate)): if coordinate[q][0]>=0: XY[0]=coordinate[q][0]/(coordinate[q][2]+1) # Calculates the x and y co-ordinates for the stereographic plot and creates a matrix to hold them (XY). XY[1]=coordinate[q][1]/(coordinate[q][2]+1) XYco.append(XY)

Pstavroulakis commented 5 years ago

In Elliots part:

Elliot

A is the postion that is being measureb, B is the postion to the right, C is the postion below xdim is the number of points along the x axis thus the point below is xdim away in data_list x = 1 A_angs = data_list[x] B_angs = data_list[x+1] C_angs = data_list[x+(xdim-1)]

uses the G_matrix function on the assoaciated Eular angles GA = G_Matrix(A_angs) GB = G_Matrix(B_angs) GC = G_Matrix(C_angs)

want to calculate a matrix which will describe a rotation from oriehtnation b and c to orientation a. To go from b to a first must inverse b then apply a np.matmul = matrix multiplication np.linalg.inv = inverse of of matrix GA_B = np.matmul(np.linalg.inv(GB),GA) GA_C = np.matmul(np.linalg.inv(GC),GA)

need to multiply the rotation from a to b&c by the family of symmetry rotations for the system

FamilyPlaneCubic()

GA_B_PF = [] GA_C_PF = [] for n in range(len(familyplane)): GA_B_PF.append(np.matmul(GA_B,familyplane[n])) GA_C_PF.append(np.matmul(GA_B,familyplane[n]))

The line: FamilyPlaneCubic() is not required and produces an error. familyplane[] has already been generated further up the code so simply calling it will be fine.

ElliotCN commented 5 years ago

In Elliots part:

Elliot

A is the postion that is being measureb, B is the postion to the right, C is the postion below xdim is the number of points along the x axis thus the point below is xdim away in data_list x = 1 A_angs = data_list[x] B_angs = data_list[x+1] C_angs = data_list[x+(xdim-1)]

uses the G_matrix function on the assoaciated Eular angles GA = G_Matrix(A_angs) GB = G_Matrix(B_angs) GC = G_Matrix(C_angs)

want to calculate a matrix which will describe a rotation from oriehtnation b and c to orientation a. To go from b to a first must inverse b then apply a np.matmul = matrix multiplication np.linalg.inv = inverse of of matrix GA_B = np.matmul(np.linalg.inv(GB),GA) GA_C = np.matmul(np.linalg.inv(GC),GA)

need to multiply the rotation from a to b&c by the family of symmetry rotations for the system

FamilyPlaneCubic()

GA_B_PF = [] GA_C_PF = [] for n in range(len(familyplane)): GA_B_PF.append(np.matmul(GA_B,familyplane[n])) GA_C_PF.append(np.matmul(GA_B,familyplane[n]))

The line: FamilyPlaneCubic() is not required and produces an error. familyplane[] has already been generated further up the code so simply calling it will be fine.

I forgot we had swapped the stereograph and misorientation sections so I ended up doing some of the misorientation which isn't actually my assigned section. Was gonna delete it but thought it might be useful.