Skip to content
Snippets Groups Projects
Commit 8d5fb8b3 authored by Alán Muñoz's avatar Alán Muñoz
Browse files

bugfixes

Former-commit-id: 14c7b6bdfe4a278930f9a58d6345f08af17c18ad
parent 0856b9dd
No related branches found
No related tags found
No related merge requests found
...@@ -229,13 +229,13 @@ def get_joinable(tracks, smooth=False, tol=0.1, window=5, degree=3) -> dict: ...@@ -229,13 +229,13 @@ def get_joinable(tracks, smooth=False, tol=0.1, window=5, degree=3) -> dict:
savgol_on_srs = lambda x: non_uniform_savgol(x.index, x.values, window, degree) savgol_on_srs = lambda x: non_uniform_savgol(x.index, x.values, window, degree)
contig = clean.groupby(["trap"]).apply(get_contiguous_pairs) contig = clean.groupby(["trap"]).apply(get_contiguous_pairs)
contig = contig.loc[contig.apply(len) > 0] contig = contig.loc[contig.apply(len) > 0]
linear = set([k for v in contig.values for i in v for j in i for k in j]) flat = set([k for v in contig.values for i in v for j in i for k in j])
smoothed_tracks = clean.loc[linear].apply(savgol_on_srs, 1) smoothed_tracks = clean.loc[flat].apply(savgol_on_srs, 1)
else: else:
contig = tracks.groupby(["trap"]).apply(get_contiguous_pairs) contig = tracks.groupby(["trap"]).apply(get_contiguous_pairs)
contig = contig.loc[contig.apply(len) > 0] contig = contig.loc[contig.apply(len) > 0]
linear = set([k for v in contig.values for i in v for j in i for k in j]) flat = set([k for v in contig.values for i in v for j in i for k in j])
smoothed_tracks = tracks.loc[linear].apply(lambda x: np.array(x.values), axis=1) smoothed_tracks = tracks.loc[flat].apply(lambda x: np.array(x.values), axis=1)
# fetch edges from ids TODO (IF necessary, here we can compare growth rates) # fetch edges from ids TODO (IF necessary, here we can compare growth rates)
idx_to_edge = lambda preposts: [ idx_to_edge = lambda preposts: [
......
...@@ -114,14 +114,15 @@ class PostProcessor: ...@@ -114,14 +114,15 @@ class PostProcessor:
with h5py.File(self._filename, "r") as f: with h5py.File(self._filename, "r") as f:
prev_idchanges = self._signal.get_id_changes() prev_idchanges = self._signal.get_id_changes()
changes_history = prev_idchanges + [merge_events] # + [picks] changes_history = prev_idchanges + merge_events
self._writer.write("/id_changes", data=changes_history) self._writer.write("id_changes", data=changes_history)
# self._writer.write( # self._writer.write(
# "/postprocessing/merge_events/", # "/postprocessing/merge_events/",
# data=merge_events, # data=merge_events,
# meta={"source": "/cell_info/"}, # meta={"source": "/cell_info/"},
# ) # )
# changes_history += picks # changes_history += picks
# TODO add picks dataset
# picks = self.picker.run(self._signal[self.targets["prepost"]["picker"][0]]) # picks = self.picker.run(self._signal[self.targets["prepost"]["picker"][0]])
# self._writer.write() # self._writer.write()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment