for frameID in xrange(startID+1, stopID): # Get current frame frame1 = loadFrame(frameManager, frameID, matrixID) active_cells1 = frameManager.get_num_active_cells_matrix(frameID, matrixID) # Compute slip vector if (active_cells0 > thresh_active_cells_translation and active_cells1 > thresh_active_cells_translation): slipvector = NCC.normalized_cross_correlation(frame0, frame1) slipvectors_delta = np.vstack((slipvectors_delta, slipvector)) slipvectors = np.vstack((slipvectors, slipvectors[-1]+slipvector)) slipvector_ncc_1 = NCC.normalized_cross_correlation2(frame0, frame1) slipvectors_ncc_1_delta = np.vstack((slipvectors_ncc_1_delta, slipvector_ncc_1)) slipvectors_ncc_1 = np.vstack((slipvectors_ncc_1, slipvectors_ncc_1[-1]+slipvector_ncc_1)) slipvector_ncc_2 = NCC.normalized_cross_correlation3(frame0, frame1) slipvectors_ncc_2_delta = np.vstack((slipvectors_ncc_2_delta, slipvector_ncc_2)) slipvectors_ncc_2 = np.vstack((slipvectors_ncc_2, slipvectors_ncc_2[-1]+slipvector_ncc_2)) slipvector_pc = NCC.normalized_cross_correlation4(frame0, frame1) slipvectors_pc_delta = np.vstack((slipvectors_pc_delta, slipvector_pc)) slipvectors_pc = np.vstack((slipvectors_pc, slipvectors_pc[-1]+slipvector_pc)) frame0 = frame1 active_cells0 = active_cells1 else: slipvectors_delta = np.vstack((slipvectors_delta, np.zeros(2))) slipvectors = np.vstack((slipvectors, slipvectors[-1])) slipvectors_ncc_1_delta = np.vstack((slipvectors_ncc_1_delta, np.zeros(2))) slipvectors_ncc_1 = np.vstack((slipvectors_ncc_1, slipvectors_ncc_1[-1]))
elapsed[0] = (time.clock() - start) # NCC 1 start = time.clock() for frameID in xrange(0, len_sequence): frame1 = tsframes3D[:, :, rand_frame_idx[frameID]] slipvectors_ncc_1[frameID] = NCC.normalized_cross_correlation2( frame0, frame1) frame0 = frame1 elapsed[1] = (time.clock() - start) # NCC 2 start = time.clock() for frameID in xrange(0, len_sequence): frame1 = tsframes3D[:, :, rand_frame_idx[frameID]] slipvectors_ncc_2[frameID] = NCC.normalized_cross_correlation3( frame0, frame1) frame0 = frame1 elapsed[2] = (time.clock() - start) # PC start = time.clock() for frameID in xrange(0, len_sequence): frame1 = tsframes3D[:, :, rand_frame_idx[frameID]] slipvectors_pc[frameID] = NCC.normalized_cross_correlation4(frame0, frame1) frame0 = frame1 elapsed[3] = (time.clock() - start) # Centroid start = time.clock() for frameID in xrange(0, len_sequence): frame1 = tsframes3D[:, :, rand_frame_idx[frameID]]