Draft: Develop

Improvements:

• Normalization using AIS_Points now raise exception if given a bad dictionary

New features:

• Added choosing fields for generating images for adding channels
• Added centering of image on last point
• Added chosing bounding box of image
• Added selecting range of values for features

skais/process/basic_features_operations.py

 import cmath
-
+from numba import jit
 import numpy as np
 
 
+@jit(nopython=True, fastmath=True)
 def angular_average_vector(angles):
-    n = len(angles)
-    x = np.sum(np.cos(angles)) / n
-    y = np.sum(np.sin(angles)) / n
+	n = len(angles)
+	x = np.sum(np.cos(angles)) / n
+	y = np.sum(np.sin(angles)) / n
 
-    return np.array([x, y])
+	return np.array([x, y])
 
 
+@jit(nopython=True, fastmath=True)
 def angular_dispersion(angles):
-    x, y = angular_average_vector(angles)
-    return np.sqrt(x ** 2 + y ** 2)
+	x, y = angular_average_vector(angles)
+	return np.sqrt(x ** 2 + y ** 2)
 
 
+@jit(nopython=True, fastmath=True)
 def angular_mean(angles):
-    x, y = angular_average_vector(angles)
-    theta = np.arctan(y/x)
+	x, y = angular_average_vector(angles)
+	theta = np.arctan(y / x)
 
-    if y > 0 and x > 0:
-        return theta
-    elif x <= 0:
-        return np.pi + theta
-    else:
-        return 2*np.pi + theta
+	if y > 0 and x > 0:
+		return theta
+	elif x <= 0:
+		return np.pi + theta
+	else:
+		return 2 * np.pi + theta
 
 
+@jit(nopython=True, fastmath=True)
 def angular_std(angles):
-    return 1 - angular_dispersion(angles)
+	return 1 - angular_dispersion(angles)
 
 
+@jit(nopython=True, fastmath=True)
 def angular_time_derivative(angle1, angle2, ts1, ts2):
-    z1 = np.cos(angle1) + np.sin(angle1) * 1j
-    z2 = np.cos(angle2) + np.sin(angle2) * 1j
+	z1 = np.cos(angle1) + np.sin(angle1) * 1j
+	z2 = np.cos(angle2) + np.sin(angle2) * 1j
 
-    z = z2 / z1
-    return cmath.polar(z)[1] / (ts2 - ts1)
+	z = z2 / z1
+	return cmath.polar(z)[1] / (ts2 - ts1)
 
 
+@jit(nopython=True, fastmath=True)
 def time_derivative(f1, f2, ts1, ts2):
-    return (f2 - f1) / (ts2 - ts1)
+	return (f2 - f1) / (ts2 - ts1)