#5 The timing of ndi should include conversions

jamesavery · Mar 7, 2024 · d324edc · d324edc
1 parent b41402e
commit d324edc
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/processing_steps/0900_generate_gauss_c.py b/src/processing_steps/0900_generate_gauss_c.py
@@ -101,8 +101,8 @@ def gauss_kernel(sigma):
         Image.fromarray(toint((np.max(np.abs(result),axis=2)!=0).astype(float))).save(f'{output_dir}/{sample}-gauss-yz-nonzero.png')
 
     if verify: # generate ndimage comparison
-        control = (implant_mask>0).astype(internal_type)
         start = timeit.default_timer()
+        control = (implant_mask>0).astype(internal_type)
         for _ in tqdm(range(reps), desc='ndimage repititions'):
             control[:] = ndi.gaussian_filter(control, sigma_voxels, mode='constant', cval=0)
             control[implant_mask] = 1 # Illuminate