May I ask about these few lines:

# getting two points with maximum distance to norm angle by them
# this is needed for rotation invariant feature
am = r_array.argmax()
max_points = [am / t_points, am % t_points]

argmax() returns the position of the largest value (distance) in r_array. So what does [am/t_points, am%t_points] do? How does this give two points with the max dist? Isn’t this supposed to just return you two elements (I got a floating point value for max_points[0])? Or did I misunderstand? (My “points” variable is a 150 by 150 2D matrix (I have selected 150 points).

Any help or advice would be appreciated, thanks!