diff --git a/debugging/test_session_alignment.py b/debugging/test_session_alignment.py index dade7692a4..385c7c99a9 100644 --- a/debugging/test_session_alignment.py +++ b/debugging/test_session_alignment.py @@ -45,14 +45,24 @@ # list out carefully all notes # handle the case where the passed recordings are not motion correction recordings. -SAVE = False + +# 1) get all commits and PRs in order. Work on the original PR +# 2) investigate why the expected peaks do not drop when recording_amplitude_scalings (rename) is used +# 3) think about and add new neurons that are introduced when shifted + +# 4) add interpolation of the histograms prior to cross correlation +# 5) add robust cross-correlation +# 6) add trimmed methods +# 7) add better way to estimate chunk length. + +SAVE = True PLOT = False BIN_UM = 2 if SAVE: scalings = [np.ones(25), np.r_[np.zeros(10), np.ones(15)]] # TODO: there is something wrong here, because why are the maximum histograms not removed? recordings_list, _ = generate_session_displacement_recordings( - non_rigid_gradient=None, # 0.05, # None, + non_rigid_gradient=0.1, # 0.05, # None, num_units=15, recording_durations=(100, 100, 100, 100), recording_shifts=( @@ -82,7 +92,7 @@ corrected_recordings_list, motion_objects_list, extra_info = session_alignment.run_inter_session_displacement_correction( - recordings_list, peaks_list, peak_locations_list, bin_um=BIN_UM, histogram_estimation_method="entire_session", alignment_method="mean_crosscorr", rigid=False, log_scale=True, num_nonrigid_bins=7 + recordings_list, peaks_list, peak_locations_list, bin_um=BIN_UM, histogram_estimation_method="entire_session", alignment_method="mean_crosscorr", rigid=False, log_scale=True, num_nonrigid_bins=24 ) plotting.SessionAlignmentWidget(