@@ -212,13 +212,14 @@ def update_temporal_components(Y, A, b, Cin, fin, bl=None, c1=None, g=None, sn=N
212
212
YrA = YA - AA .T .dot (Cin ).T
213
213
# creating the patch of components to be computed in parallel
214
214
parrllcomp , len_parrllcomp = caiman .source_extraction .cnmf .utilities .update_order_greedy (AA [:nr , :][:, :nr ])
215
- logger .info ("entering the deconvolution " )
215
+ logger .info (f"Entering Deconvolution, { C . shape = } { S . shape = } " )
216
216
C , S , bl , YrA , c1 , sn , g , lam = update_iteration (parrllcomp , len_parrllcomp , nb , C , S , bl , nr ,
217
217
ITER , YrA , c1 , sn , g , Cin , T , nA , dview , debug , AA , kwargs )
218
218
ff = np .where (np .sum (C , axis = 1 ) == 0 ) # remove empty components
219
219
if np .size (ff ) > 0 : # Eliminating empty temporal components
220
220
ff = ff [0 ]
221
- logger .info (f'removing { len (ff )} empty spatial component(s)' )
221
+ logger .info (f'removing { len (ff )} empty temporal component(s)' )
222
+
222
223
keep = list (range (A .shape [1 ]))
223
224
for i in ff :
224
225
keep .remove (i )
@@ -258,9 +259,6 @@ def update_iteration(parrllcomp, len_parrllcomp, nb, C, S, bl, nr,
258
259
Cin: np.ndarray
259
260
current estimate of temporal components (K x T)
260
261
261
- g: np.ndarray
262
- Global time constant (not used)
263
-
264
262
bl: np.ndarray
265
263
baseline for fluorescence trace for each column in A
266
264
@@ -314,18 +312,21 @@ def update_iteration(parrllcomp, len_parrllcomp, nb, C, S, bl, nr,
314
312
bl: float
315
313
same as input
316
314
315
+ YrA: np.ndarray
316
+ matrix of spatial component filtered raw data, after all contributions have been removed.
317
+ YrA corresponds to the residual trace for each component and is used for faster plotting (K x T)
318
+
317
319
c1: float
318
320
same as input
319
321
320
- g: float
322
+ sn: float
321
323
same as input
322
324
323
- sn: float
325
+ g: float
324
326
same as input
325
327
326
- YrA: np.ndarray
327
- matrix of spatial component filtered raw data, after all contributions have been removed.
328
- YrA corresponds to the residual trace for each component and is used for faster plotting (K x T)
328
+ lam: (undocumented)
329
+
329
330
"""
330
331
331
332
logger = logging .getLogger ("caiman" )
0 commit comments