spacr 0.3.0__py3-none-any.whl → 0.3.1__py3-none-any.whl

spacr/core.py

@@ -1672,9 +1672,9 @@ def preprocess_generate_masks(src, settings={}):
             time_ls=[]
             if check_mask_folder(src, 'cell_mask_stack'):
                 start = time.time()
-                if settings['segmantation_model'] == 'cellpose':
+                if settings['segmentation_mode'] == 'cellpose':
                     generate_cellpose_masks(mask_src, settings, 'cell')
-                elif settings['segmantation_model'] == 'mediar':
+                elif settings['segmentation_mode'] == 'mediar':
                     generate_mediar_masks(mask_src, settings, 'cell')
                 stop = time.time()
                 duration = (stop - start)
@@ -1686,9 +1686,9 @@ def preprocess_generate_masks(src, settings={}):
             time_ls=[]
             if check_mask_folder(src, 'nucleus_mask_stack'):
                 start = time.time()
-                if settings['segmantation_model'] == 'cellpose':
+                if settings['segmentation_mode'] == 'cellpose':
                     generate_cellpose_masks(mask_src, settings, 'nucleus')
-                elif settings['segmantation_model'] == 'mediar':
+                elif settings['segmentation_mode'] == 'mediar':
                     generate_mediar_masks(mask_src, settings, 'nucleus')
                 stop = time.time()
                 duration = (stop - start)
@@ -1700,9 +1700,9 @@ def preprocess_generate_masks(src, settings={}):
             time_ls=[]
             if check_mask_folder(src, 'pathogen_mask_stack'):
                 start = time.time()
-                if settings['segmantation_model'] == 'cellpose':
+                if settings['segmentation_mode'] == 'cellpose':
                     generate_cellpose_masks(mask_src, settings, 'pathogen')
-                elif settings['segmantation_model'] == 'mediar':
+                elif settings['segmentation_mode'] == 'mediar':
                     generate_mediar_masks(mask_src, settings, 'pathogen')
                 stop = time.time()
                 duration = (stop - start)
@@ -1757,7 +1757,7 @@ def preprocess_generate_masks(src, settings={}):
                             files_processed = i+1
                             files_to_process = settings['examples_to_plot']
                             print_progress(files_processed, files_to_process, n_jobs=1, time_ls=time_ls, batch_size=None, operation_type="Plot mask outlines")
-                            print("Successfully completed run")
+                            
                 except Exception as e:
                     print(f'Failed to plot image mask overly. Error: {e}')
             else:

spacr/io.py

@@ -653,7 +653,6 @@ def _rename_and_organize_image_files(src, regex, batch_size=100, pick_slice=Fals
     from .utils import _extract_filename_metadata, print_progress
     
     regular_expression = re.compile(regex)
-    images_by_key = defaultdict(list)
     stack_path = os.path.join(src, 'stack')
     files_processed = 0
     if not os.path.exists(stack_path) or (os.path.isdir(stack_path) and len(os.listdir(stack_path)) == 0):
@@ -665,7 +664,7 @@ def _rename_and_organize_image_files(src, regex, batch_size=100, pick_slice=Fals
             start = time.time()
             batch_filenames = all_filenames[idx:idx+batch_size]
             for filename in batch_filenames:
-                images_by_key = _extract_filename_metadata(batch_filenames, src, images_by_key, regular_expression, metadata_type, pick_slice, skip_mode)
+                images_by_key = _extract_filename_metadata(batch_filenames, src, regular_expression, metadata_type, pick_slice, skip_mode)
             
             if pick_slice:
                 for i, key in enumerate(images_by_key):
@@ -683,16 +682,16 @@ def _rename_and_organize_image_files(src, regex, batch_size=100, pick_slice=Fals
                     files_to_process = len(all_filenames)
                     print_progress(files_processed, files_to_process, n_jobs=1, time_ls=time_ls, batch_size=batch_size, operation_type='Preprocessing filenames')
 
-                    if os.path.exists(output_path):                        
-                        print(f'WARNING: A file with the same name already exists at location {output_filename}')
-                    else:
+                    #if os.path.exists(output_path):                        
+                    #    print(f'WARNING: A file with the same name already exists at location {output_filename}')
+                    if not os.path.exists(output_path):
                         mip_image.save(output_path)
             else:
                 for i, (key, images) in enumerate(images_by_key.items()):
-                    mip = np.max(np.stack(images), axis=0)
-                    mip_image = Image.fromarray(mip)
                     plate, well, field, channel = key[:4]
                     output_dir = os.path.join(src, channel)
+                    mip = np.max(np.stack(images), axis=0)
+                    mip_image = Image.fromarray(mip)
                     os.makedirs(output_dir, exist_ok=True)
                     output_filename = f'{plate}_{well}_{field}.tif'
                     output_path = os.path.join(output_dir, output_filename)
@@ -703,9 +702,9 @@ def _rename_and_organize_image_files(src, regex, batch_size=100, pick_slice=Fals
                     files_to_process = len(all_filenames)
                     print_progress(files_processed, files_to_process, n_jobs=1, time_ls=time_ls, batch_size=batch_size, operation_type='Preprocessing filenames')
 
-                    if os.path.exists(output_path):                        
-                        print(f'WARNING: A file with the same name already exists at location {output_filename}')
-                    else:
+                    #if os.path.exists(output_path):                        
+                    #    print(f'WARNING: A file with the same name already exists at location {output_filename}')
+                    if not os.path.exists(output_path):
                         mip_image.save(output_path)
             images_by_key.clear()
 
@@ -972,8 +971,8 @@ def _mip_all(src, include_first_chan=True):
         if filename.endswith('.npy'):
             # Load the array from the file.
             array = np.load(os.path.join(src, filename))
-            # Normalize the array using custom parameters (q1=2, q2=98).
-            array = normalize_to_dtype(array, q1=2, q2=98, percentiles=None)
+            # Normalize the array 
+            #array = normalize_to_dtype(array, q1=0, q2=99, percentiles=None)
 
             if array.ndim != 3: # Check if the array is not 3-dimensional.
                 # Log a message indicating a zero array will be generated due to unexpected dimensions.
@@ -1671,6 +1670,7 @@ def preprocess_img_data(settings):
                 if plot:
                     print(f'plotting {nr} images from {src}/stack')
                     plot_arrays(src+'/stack', figuresize, cmap, nr=nr, normalize=normalize)
+
                 if all_to_mip:
                     _mip_all(src+'/stack')
                     if plot:

spacr/mediar.py

@@ -22,8 +22,6 @@ try:
     # Now import the dependencies from MEDIAR
     from core.MEDIAR import Predictor, EnsemblePredictor
     from train_tools.models import MEDIARFormer
-
-    print("Imports successful.")
 finally:
     # Remove the temporary __init__.py file after the import
     if os.path.exists(init_file):

spacr/plot.py

@@ -123,7 +123,7 @@ def plot_image_mask_overlay(file, channels, cell_channel, nucleus_channel, patho
 
     fig = _plot_merged_plot(image=image, outlines=outlines, outline_colors=outline_colors, figuresize=figuresize, thickness=thickness)
 
-    return
+    return fig
 
 def plot_masks(batch, masks, flows, cmap='inferno', figuresize=10, nr=1, file_type='.npz', print_object_number=True):
     """

spacr/utils.py

@@ -303,7 +303,10 @@ def _get_cellpose_batch_size():
     except Exception as e:
         return 8
 
-def _extract_filename_metadata(filenames, src, images_by_key, regular_expression, metadata_type='cellvoyager', pick_slice=False, skip_mode='01'):
+def _extract_filename_metadata(filenames, src, regular_expression, metadata_type='cellvoyager', pick_slice=False, skip_mode='01'):
+    
+    images_by_key = defaultdict(list)
+    
     for filename in filenames:
         match = regular_expression.match(filename)
         if match:
@@ -328,7 +331,6 @@ def _extract_filename_metadata(filenames, src, images_by_key, regular_expression
                 if metadata_type =='cq1':
                     orig_wellID = wellID
                     wellID = _convert_cq1_well_id(wellID)
-                    #clear_output(wait=True)
                     print(f'Converted Well ID: {orig_wellID} to {wellID}', end='\r', flush=True)
 
                 if pick_slice:
@@ -338,7 +340,7 @@ def _extract_filename_metadata(filenames, src, images_by_key, regular_expression
                         sliceid = '00'
 
                     if mode == skip_mode:
-                        continue      
+                        continue
                         
                 key = (plate, well, field, channel, mode)
                 with Image.open(os.path.join(src, filename)) as img:

{spacr-0.3.0.dist-info → spacr-0.3.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: spacr
-Version: 0.3.0
+Version: 0.3.1
 Summary: Spatial phenotype analysis of crisp screens (SpaCr)
 Home-page: https://github.com/EinarOlafsson/spacr
 Author: Einar Birnir Olafsson

{spacr-0.3.0.dist-info → spacr-0.3.1.dist-info}/RECORD

@@ -8,24 +8,24 @@ spacr/app_measure.py,sha256=_K7APYIeOKpV6e_LcqabBjvEi7mfq9Fch8175x1x0k8,162
 spacr/app_sequencing.py,sha256=DjG26jy4cpddnV8WOOAIiExtOe9MleVMY4MFa5uTo5w,157
 spacr/app_umap.py,sha256=ZWAmf_OsIKbYvolYuWPMYhdlVe-n2CADoJulAizMiEo,153
 spacr/chris.py,sha256=YlBjSgeZaY8HPy6jkrT_ISAnCMAKVfvCxF0I9eAZLFM,2418
-spacr/core.py,sha256=gOr5P58s4ubzDfO0Ar3Zwce0bo-E8qxdN37Ian82RGw,150744
+spacr/core.py,sha256=uyZdJ93ysd8oXgRX9b-6iTCldJ0CM1dq5VxF-xbZyN8,150703
 spacr/deep_spacr.py,sha256=a2YewgkQvLV-95NYJAutnojvJmX4S8z_wv6Tb-XIgUI,34484
 spacr/graph_learning.py,sha256=1tR-ZxvXE3dBz1Saw7BeVFcrsUFu9OlUZeZVifih9eo,13070
 spacr/gui.py,sha256=zUkIyAuOwwoMDoExxtI-QHRfOhE1R2rulXJDNxwSLGc,7947
 spacr/gui_core.py,sha256=ZUIqvK7x6NzgrmuTRbvwCTTSpU3yWUaId6MZjXv16us,40128
 spacr/gui_elements.py,sha256=OA514FUVRKAcdu9CFVOt7UEzn1vztakQ-rDyKqV0b9A,129771
 spacr/gui_utils.py,sha256=DCI--DNoYDWY1q0Aohd0XwFqjdPM3K5kCgRKiJGTnfc,30697
-spacr/io.py,sha256=o2Sxan6_HylbYNUjNdgeQ4iEpcXRwlgV1r59J2HG0xQ,117292
+spacr/io.py,sha256=VjH_1zXmf0yEdtABnsoabCEIpU0S3wB-7Hog7_ntCdE,117267
 spacr/logger.py,sha256=7Zqr3TuuOQLWT32gYr2q1qvv7x0a2JhLANmZcnBXAW8,670
 spacr/measure.py,sha256=ooMOP2OE0BHUNqIkg0ltwV2FiO6hZDIcRC6A0YmGcws,54875
-spacr/mediar.py,sha256=KEkrO0dlAbOr2tN4c8-vDncTy3IWNzUAwWYiPRofFqU,14864
-spacr/plot.py,sha256=9Ty2cMJpICXPCEE4inGTb1FvlUZkVgqB7lqJaggmjZE,73975
+spacr/mediar.py,sha256=5HaCyZYiOff74PCvHwKj-jSRua0QRoIv1mvElPfVKtY,14830
+spacr/plot.py,sha256=yFC5m54lB8xVnC7vQp50-FvRn6lCjCt2mgi2GeiRwSs,73979
 spacr/sequencing.py,sha256=92KmjFa8Ptwmpf-GtyH3-uX6djFOYR5lJjMBHeciqhs,66921
 spacr/settings.py,sha256=PfIPLyMyBAfOodtdgNT8QzbysNDxTnsONXdI-fKtIDQ,68038
 spacr/sim.py,sha256=FveaVgBi3eypO2oVB5Dx-v0CC1Ny7UPfXkJiiRRodAk,71212
 spacr/sim_app.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 spacr/timelapse.py,sha256=KMYCgHzf9LTZe-lWl5mvH2EjbKRE6OhpwdY13wEumGc,39504
-spacr/utils.py,sha256=D13hYJsVuq1ZMWYGyvtEskU_W2P8ddsI9zX9pTtMAFA,189159
+spacr/utils.py,sha256=JGjL_Tg5ec1qmaf9BV3gqlREUpxjK8hutgPEGj5mAEs,189141
 spacr/version.py,sha256=axH5tnGwtgSnJHb5IDhiu4Zjk5GhLyAEDRe-rnaoFOA,409
 spacr/resources/MEDIAR/.git,sha256=nHbNNUgehWnXyS2LbJZitX4kbpd1urzYgE0WZYvdMfc,53
 spacr/resources/MEDIAR/.gitignore,sha256=Ff1q9Nme14JUd-4Q3jZ65aeQ5X4uttptssVDgBVHYo8,152
@@ -151,9 +151,9 @@ spacr/resources/images/plate1_E01_T0001F001L01A03Z01C03.tif,sha256=Pbhk7xn-KUP6R
 spacr/resources/models/cp/toxo_plaque_cyto_e25000_X1120_Y1120.CP_model,sha256=z8BbHWZPRnE9D_BHO0fBREE85c1vkltDs-incs2ytXQ,26566572
 spacr/resources/models/cp/toxo_plaque_cyto_e25000_X1120_Y1120.CP_model_settings.csv,sha256=fBAGuL_B8ERVdVizO3BHozTDSbZUh1yFzsYK3wkQN68,420
 spacr/resources/models/cp/toxo_pv_lumen.CP_model,sha256=2y_CindYhmTvVwBH39SNILF3rI3x9SsRn6qrMxHy3l0,26562451
-spacr-0.3.0.dist-info/LICENSE,sha256=SR-2MeGc6SCM1UORJYyarSWY_A-JaOMFDj7ReSs9tRM,1083
-spacr-0.3.0.dist-info/METADATA,sha256=gPhwJrEAVxewekQKy4X9qiLlDr3Q4JvwXnq941vcpew,5646
-spacr-0.3.0.dist-info/WHEEL,sha256=HiCZjzuy6Dw0hdX5R3LCFPDmFS4BWl8H-8W39XfmgX4,91
-spacr-0.3.0.dist-info/entry_points.txt,sha256=BMC0ql9aNNpv8lUZ8sgDLQMsqaVnX5L535gEhKUP5ho,296
-spacr-0.3.0.dist-info/top_level.txt,sha256=GJPU8FgwRXGzKeut6JopsSRY2R8T3i9lDgya42tLInY,6
-spacr-0.3.0.dist-info/RECORD,,
+spacr-0.3.1.dist-info/LICENSE,sha256=SR-2MeGc6SCM1UORJYyarSWY_A-JaOMFDj7ReSs9tRM,1083
+spacr-0.3.1.dist-info/METADATA,sha256=Fi3B8Vxgz4IzfYodV9zMgRdFLXTtVtel4OCdE06o5cI,5646
+spacr-0.3.1.dist-info/WHEEL,sha256=HiCZjzuy6Dw0hdX5R3LCFPDmFS4BWl8H-8W39XfmgX4,91
+spacr-0.3.1.dist-info/entry_points.txt,sha256=BMC0ql9aNNpv8lUZ8sgDLQMsqaVnX5L535gEhKUP5ho,296
+spacr-0.3.1.dist-info/top_level.txt,sha256=GJPU8FgwRXGzKeut6JopsSRY2R8T3i9lDgya42tLInY,6
+spacr-0.3.1.dist-info/RECORD,,