From e7863f0b0ba56e1ce7318489f76a2520179c83e1 Mon Sep 17 00:00:00 2001 From: EricThomson Date: Fri, 22 Mar 2024 11:46:01 -0400 Subject: [PATCH] fixed bug in correlation image calculation --- caiman/summary_images.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/caiman/summary_images.py b/caiman/summary_images.py index f9a841d35..c40af03da 100644 --- a/caiman/summary_images.py +++ b/caiman/summary_images.py @@ -201,7 +201,7 @@ def local_correlations(Y, eight_neighbours: bool = True, swap_dim: bool = True, Returns: rho: d1 x d2 [x d3] matrix, cross-correlation with adjacent pixels """ - + print("Testing corr func") if swap_dim: Y = np.transpose(Y, tuple(np.hstack((Y.ndim - 1, list(range(Y.ndim))[:-1])))) @@ -248,14 +248,14 @@ def local_correlations(Y, eight_neighbours: bool = True, swap_dim: bool = True, rho_d1 = rho_d1 rho_d2 = rho_d2 rho[:-1, :-1] = rho[:-1, :-1] * rho_d2 - rho[1:, 1:] = rho[1:, 1:] * rho_d1 + rho[1:, 1:] = rho[1:, 1:] * rho_d2 rho[1:, :-1] = rho[1:, :-1] * rho_d1 - rho[:-1, 1:] = rho[:-1, 1:] * rho_d2 + rho[:-1, 1:] = rho[:-1, 1:] * rho_d1 else: rho[:-1, :-1] = rho[:-1, :-1] + rho_d2**(order_mean) - rho[1:, 1:] = rho[1:, 1:] + rho_d1**(order_mean) + rho[1:, 1:] = rho[1:, 1:] + rho_d2**(order_mean) rho[1:, :-1] = rho[1:, :-1] + rho_d1**(order_mean) - rho[:-1, 1:] = rho[:-1, 1:] + rho_d2**(order_mean) + rho[:-1, 1:] = rho[:-1, 1:] + rho_d1**(order_mean) neighbors = 8 * np.ones(np.shape(Y)[1:3]) neighbors[0, :] = neighbors[0, :] - 3