estoy tratando de hacer restas entre filas de 2 columnas. Para esto uso una funcion que resta la fila1-fila2, la eleva al cuadrado y le saca raíz cuadrada para eliminar signo negativo. Después se va a la fila2-fila3 y así hasta el final de los datos. Posteriormente la función calcula lo mismo pero para un separación de 2 filas (fila1-fila3), 3 filas y así hasta filas-1. Los resultados se guardan en dataFrame. Los datos con los que estoy trabajando están acomodados así:
df1
Out[44]:
TRACK_ID POSITION_X POSITION_Y POSITION_T
0 0 1 1 35.36
1 0 2 2 35.52
2 0 3 3 35.68
3 0 4 4 35.84
4 0 1 1 35.36
5 0 4 3 34.88
6 0 2 3 34.40
7 0 6 4 33.92
8 0 4 2 33.44
La función parece funcionar, el problema es que me di cuenta que para algunos valores de separación entre las restas, me repite exactamente los mimos valores de la columna previa. Ejemplo:
rad
Out[28]:
0 1 2 3 4 5 6 \
0 1.414214 2.828427 4.242641 0.000000 3.605551 2.236068 2.236068
1 1.414214 2.828427 1.414214 2.236068 1.000000 4.472136 4.472136
2 1.414214 2.828427 1.000000 1.000000 3.162278 1.414214 1.414214
3 4.242641 1.000000 2.236068 2.000000 2.000000 NaN NaN
4 3.605551 2.236068 5.830952 3.162278 NaN NaN NaN
5 2.000000 2.236068 1.000000 NaN NaN NaN NaN
6 4.123106 2.236068 NaN NaN NaN NaN NaN
7 2.828427 NaN NaN NaN NaN NaN NaN
8 NaN NaN NaN NaN NaN NaN NaN
7 8
0 3.162278 NaN
1 NaN NaN
2 NaN NaN
3 NaN NaN
4 NaN NaN
5 NaN NaN
6 NaN NaN
7 NaN NaN
8 NaN NaN
Columna 5 y 6 es lo mismo.
Este es mi código completo:
df1 = df[['TRACK_ID','POSITION_X','POSITION_Y','POSITION_T']].copy()
#Parameter input
N = df1.groupby('TRACK_ID').size()
max_time = N*(0.160)
frames = max_time/N
t_step=frames.item()
data = pd.DataFrame({'N':N,'max_time':max_time,'frames':frames})
print(data)
t=np.linspace(0.160, max_time.item(), N)
#funcion para calcular las diferencias
def radial(df1, coords=['POSITION_X', 'POSITION_Y']):
tau = t.copy()
shifts = np.divide(tau,t_step).astype(float) #matrix que se ocupa para construir las diferencias entre valores de filas
print(shifts)
radials = list()
for i, shift in enumerate(shifts):
diffs = np.array(df1[coords] - df1[coords].shift(-shift))
sqdist = np.square(diffs).sum(axis=1)
r = np.sqrt(sqdist)
radials.append(r)
radial_disp = pd.DataFrame({'radials':radials})
return radials
radial_d = radial(df1, coords=['POSITION_X', 'POSITION_Y'])
radd = pd.DataFrame.from_records(radial_d) #horizontal
rad = radd.transpose() #vertical
Ya modifique algunas partes de la función y me había dado cuenta que mi variable shifts
que establece la separación entre las restas me arrojaba resultados repetidos porque estaban como int
lo arregle poniendo float
pero el resultado sigue siendo el mismo. Porque se repite el computo para la misma columna?? Gracias por leer mi post
df1
, o al menos un ejemplo con la salida esperada para crear un ejemplo mínimo verificable que podamos reproducir tu código. Por otro lado, siempre es positivo añadir comentarios al código explicando que es cada variable y que se hace en cada momento para facilitar la comprensión del código. Saludos.