spacr 0.3.81__py3-none-any.whl → 0.4.1__py3-none-any.whl

spacr/toxo.py

@@ -24,117 +24,128 @@ from sklearn.metrics import mean_absolute_error
 
 from matplotlib.gridspec import GridSpec
 
+def custom_volcano_plot(data_path, metadata_path, metadata_column='tagm_location',point_size=50, figsize=20, threshold=0,save_path=None, x_lim=[-0.5, 0.5], y_lims=[[0, 6], [9, 20]]):
 
-def custom_volcano_plot(data_path, metadata_path, metadata_column='tagm_location',
-                        point_size=50, figsize=20, threshold=0,
-                        save_path=None, x_lim=[-0.5, 0.5], y_lims=[[0, 6], [9, 15]]):
+    # Dictionary mapping compartment to color
     
-    markers = [
-        'o',  # Circle
-        'X',  # X-shaped marker
-        '^',  # Upward triangle
-        's',  # Square
-        'v',  # Downward triangle
-        'P',  # Plus-filled pentagon
-        '*',  # Star
-        '+',  # Plus
-        'x',  # Cross
-        '.',  # Point
-        ',',  # Pixel
-        'd',  # Diamond
-        'D',  # Thin diamond
-        'h',  # Hexagon 1
-        'H',  # Hexagon 2
-        'p',  # Pentagon
-        '|',  # Vertical line
-        '_',  # Horizontal line
-    ]
+    colors = {'micronemes':'black',
+              'rhoptries 1':'darkviolet',
+              'rhoptries 2':'darkviolet',
+              'nucleus - chromatin':'blue',
+              'nucleus - non-chromatin':'blue',
+              'dense granules':'teal',
+              'ER 1':'pink',
+              'ER 2':'pink',
+              'unknown':'black',
+              'tubulin cytoskeleton':'slategray',
+              'IMC':'slategray',
+              'PM - peripheral 1':'slategray',
+              'PM - peripheral 2':'slategray',
+              'cytosol':'turquoise',
+              'mitochondrion - soluble':'red',
+              'mitochondrion - membranes':'red',
+              'apicoplast':'slategray',
+              'Golgi':'green',
+              'PM - integral':'slategray',
+              'apical 1':'orange',
+              'apical 2':'orange',
+              '19S proteasome':'slategray',
+              '20S proteasome':'slategray',
+              '60S ribosome':'slategray',
+              '40S ribosome':'slategray',
+              }
+
+    # Increase font size for better readability
+    fontsize = 18
+    plt.rcParams.update({'font.size': fontsize})
     
-    plt.rcParams.update({'font.size': 14})
-
-    # Load data
+    # --- Load data ---
     if isinstance(data_path, pd.DataFrame):
         data = data_path
     else:
         data = pd.read_csv(data_path)
-
-    fontsize = 18
-
-    plt.rcParams.update({'font.size': fontsize})
+        
+    # Extract ‘variable’ and ‘gene_nr’ from your feature notation
     data['variable'] = data['feature'].str.extract(r'\[(.*?)\]')
     data['variable'].fillna(data['feature'], inplace=True)
     data['gene_nr'] = data['variable'].str.split('_').str[0]
     data = data[data['variable'] != 'Intercept']
 
-    # Load metadata
+    # --- Load metadata ---
     if isinstance(metadata_path, pd.DataFrame):
         metadata = metadata_path
     else:
         metadata = pd.read_csv(metadata_path)
+
     metadata['gene_nr'] = metadata['gene_nr'].astype(str)
     data['gene_nr'] = data['gene_nr'].astype(str)
 
-    merged_data = pd.merge(data, metadata[['gene_nr', metadata_column]], on='gene_nr', how='left')
+    # Merge data and metadata
+    merged_data = pd.merge(data, metadata[['gene_nr', metadata_column]], 
+                           on='gene_nr', how='left')
     merged_data[metadata_column].fillna('unknown', inplace=True)
 
-    # Define palette and markers
-    palette = {'pc': 'red', 'nc': 'green', 'control': 'white', 'other': 'gray'}
-    marker_dict = {val: marker for val, marker in zip(
-        merged_data[metadata_column].unique(), markers)}
-
-    # Create the figure with custom spacing
-    fig = plt.figure(figsize=(figsize,figsize))
+    # --- Create figure with "upper" and "lower" subplots sharing the x-axis ---
+    fig = plt.figure(figsize=(figsize, figsize))
     gs = GridSpec(2, 1, height_ratios=[1, 3], hspace=0.05)
-
     ax_upper = fig.add_subplot(gs[0])
     ax_lower = fig.add_subplot(gs[1], sharex=ax_upper)
 
     # Hide x-axis labels on the upper plot
     ax_upper.tick_params(axis='x', which='both', bottom=False, labelbottom=False)
 
+    # List to collect the variables (hits) that meet threshold criteria
     hit_list = []
 
-    # Scatter plot on both axes
+    # --- Scatter plot on both axes ---
     for _, row in merged_data.iterrows():
         y_val = -np.log10(row['p_value'])
+        
+        # Decide which axis to draw on based on the p-value
         ax = ax_upper if y_val > y_lims[1][0] else ax_lower
 
+        # Here is the main change: color by the colors dict
         ax.scatter(
-            row['coefficient'], y_val,
-            color=palette.get(row['condition'], 'gray'),
-            marker=marker_dict.get(row[metadata_column], 'o'),
-            s=point_size, edgecolor='black', alpha=0.6
+            row['coefficient'], 
+            y_val,
+            color=colors.get(row[metadata_column], 'gray'),  # <-- Use your color dict
+            marker='o',  # You can fix a single marker if desired
+            s=point_size, 
+            edgecolor='black', 
+            alpha=0.6
         )
 
-        if row['p_value'] <= 0.05 and abs(row['coefficient']) >= abs(threshold):
+        # Check significance thresholds
+        if (row['p_value'] <= 0.05) and (abs(row['coefficient']) >= abs(threshold)):
             hit_list.append(row['variable'])
 
-    # Set axis limits
+    # --- Adjust axis limits ---
     ax_upper.set_ylim(y_lims[1])
     ax_lower.set_ylim(y_lims[0])
     ax_lower.set_xlim(x_lim)
 
+    # Hide top spines
     ax_lower.spines['top'].set_visible(False)
     ax_upper.spines['top'].set_visible(False)
     ax_upper.spines['bottom'].set_visible(False)
 
-    # Set x-axis and y-axis titles
-    ax_lower.set_xlabel('Coefficient')  # X-axis title on the lower graph
-    ax_lower.set_ylabel('-log10(p-value)')  # Y-axis title on the lower graph
-    ax_upper.set_ylabel('-log10(p-value)')  # Y-axis title on the upper graph
-    
+    # Set x-axis and y-axis labels
+    ax_lower.set_xlabel('Coefficient')  
+    ax_lower.set_ylabel('-log10(p-value)')
+    ax_upper.set_ylabel('-log10(p-value)')
+
     for ax in [ax_upper, ax_lower]:
         ax.spines['right'].set_visible(False)
 
-    # Add threshold lines to both axes
+    # --- Add threshold lines to both axes ---
     for ax in [ax_upper, ax_lower]:
         ax.axvline(x=-abs(threshold), linestyle='--', color='black')
         ax.axvline(x=abs(threshold), linestyle='--', color='black')
 
     ax_lower.axhline(y=-np.log10(0.05), linestyle='--', color='black')
 
-    # Annotate significant points
-    texts_upper, texts_lower = [], []  # Collect text annotations separately
+    # --- Annotate significant points ---
+    texts_upper, texts_lower = [], []
 
     for _, row in merged_data.iterrows():
         y_val = -np.log10(row['p_value'])
@@ -142,38 +153,50 @@ def custom_volcano_plot(data_path, metadata_path, metadata_column='tagm_location
             continue
 
         ax = ax_upper if y_val > y_lims[1][0] else ax_lower
-        text = ax.text(row['coefficient'], y_val, row['variable'], 
-                       fontsize=fontsize, ha='center', va='bottom')
+        text = ax.text(
+            row['coefficient'],
+            y_val,
+            row['variable'],
+            fontsize=fontsize,
+            ha='center',
+            va='bottom'
+        )
 
         if ax == ax_upper:
             texts_upper.append(text)
         else:
             texts_lower.append(text)
 
-    # Adjust text positions to avoid overlap
+    # Attempt to keep text labels from overlapping
     adjust_text(texts_upper, ax=ax_upper, arrowprops=dict(arrowstyle='-', color='black'))
     adjust_text(texts_lower, ax=ax_lower, arrowprops=dict(arrowstyle='-', color='black'))
 
-    # Add a single legend on the lower axis
-    handles = [plt.Line2D([0], [0], marker=m, color='w', markerfacecolor='gray', markersize=10)
-               for m in marker_dict.values()]
-    labels = marker_dict.keys()
-    ax_lower.legend(handles, 
-                    labels, 
-                    bbox_to_anchor=(1.05, 1), 
-                    loc='upper left',
-                    borderaxespad=0.25, 
-                    labelspacing=2, 
-                    handletextpad=0.25, 
-                    markerscale=2, 
-                    prop={'size': fontsize})
-
-
-    # Save and show the plot
+    # --- Add a legend keyed by color (optional) ---
+    # If you'd like a legend that shows what each compartment color represents:
+    legend_handles = []
+    for comp, comp_color in colors.items():
+        # Create a “dummy” scatter for legend
+        legend_handles.append(
+            plt.Line2D([0], [0], marker='o', color=comp_color, 
+                       label=comp, linewidth=0, markersize=8)
+        )
+    # You can adjust the location and styling of the legend to taste:
+    ax_lower.legend(
+        handles=legend_handles,
+        bbox_to_anchor=(1.05, 1),
+        loc='upper left',
+        borderaxespad=0.25,
+        labelspacing=2,
+        handletextpad=0.25,
+        markerscale=1.5,
+        prop={'size': fontsize}
+    )
+
+    # --- Save and show ---
     if save_path:
         plt.savefig(save_path, format='pdf', bbox_inches='tight')
     plt.show()
-    
+
     return hit_list
 
 def go_term_enrichment_by_column(significant_df, metadata_path, go_term_columns=['Computed GO Processes', 'Curated GO Components', 'Curated GO Functions', 'Curated GO Processes']):

spacr/utils.py

@@ -1372,40 +1372,6 @@ def annotate_conditions(df, cells=None, cell_loc=None, pathogens=None, pathogen_
 
     return df
 
-def _split_data_v1(df, group_by, object_type):
-    """
-    Splits the input dataframe into numeric and non-numeric parts, groups them by the specified column,
-    and returns the grouped dataframes.
-
-    Parameters:
-    df (pandas.DataFrame): The input dataframe.
-    group_by (str): The column name to group the dataframes by.
-    object_type (str): The column name to concatenate with 'prcf' to create a new column 'prcfo'.
-
-    Returns:
-    grouped_numeric (pandas.DataFrame): The grouped dataframe containing numeric columns.
-    grouped_non_numeric (pandas.DataFrame): The grouped dataframe containing non-numeric columns.
-    """
-    
-    if 'prcf' not in df.columns:
-        try:
-            df['prcf'] = df['plate'].astype(str) + '_' + df['row_name'].astype(str) + '_' + df['column_name'].astype(str) + '_' + df['field'].astype(str)
-        except Exception as e:
-            print(e)    
-    
-    df['prcfo'] = df['prcf'] + '_' + df[object_type]
-    df = df.set_index(group_by, inplace=False)
-
-    df_numeric = df.select_dtypes(include=np.number)
-    df_non_numeric = df.select_dtypes(exclude=np.number)
-    
-    []
-
-    grouped_numeric = df_numeric.groupby(df_numeric.index).mean()
-    grouped_non_numeric = df_non_numeric.groupby(df_non_numeric.index).first()
-
-    return pd.DataFrame(grouped_numeric), pd.DataFrame(grouped_non_numeric)
-
 def _split_data(df, group_by, object_type):
     """
     Splits the input dataframe into numeric and non-numeric parts, groups them by the specified column,
@@ -5045,22 +5011,22 @@ def generate_cytoplasm_mask(nucleus_mask, cell_mask):
     return cytoplasm_mask
 
 def add_column_to_database(settings):
-    #"""
-    #Adds a new column to the database table by matching on a common column from the DataFrame.
-    #If the column already exists in the database, it adds the column with a suffix.
-    #NaN values will remain as NULL in the database.
-
-    #Parameters:
-    #    settings (dict): A dictionary containing the following keys:
-    #    csv_path (str): Path to the CSV file with the data to be added.
-    #    db_path (str): Path to the SQLite database (or connection string for other databases).
-    #    table_name (str): The name of the table in the database.
-    #    update_column (str): The name of the new column in the DataFrame to add to the database.
-    #    match_column (str): The common column used to match rows.
-
-    #Returns:
-    #    None
-    #"""
+    """
+    Adds a new column to the database table by matching on a common column from the DataFrame.
+    If the column already exists in the database, it adds the column with a suffix.
+    NaN values will remain as NULL in the database.
+
+    Parameters:
+        settings (dict): A dictionary containing the following keys:
+            csv_path (str): Path to the CSV file with the data to be added.
+            db_path (str): Path to the SQLite database (or connection string for other databases).
+            table_name (str): The name of the table in the database.
+            update_column (str): The name of the new column in the DataFrame to add to the database.
+            match_column (str): The common column used to match rows.
+
+    Returns:
+        None
+    """
 
     # Read the DataFrame from the provided CSV path
     df = pd.read_csv(settings['csv_path'])
@@ -5150,8 +5116,8 @@ def correct_metadata_column_names(df):
         df = df.rename(columns={'plate_name': 'plate'})
     if 'column_name' in df.columns:
         df = df.rename(columns={'column_name': 'column'})
-    if 'column_name' in df.columns:
-        df = df.rename(columns={'column_name': 'column'})
+    if 'col' in df.columns:
+        df = df.rename(columns={'col': 'column'})
     if 'row_name' in df.columns:
         df = df.rename(columns={'row_name': 'row_name'})
     if 'grna_name' in df.columns:
@@ -5245,3 +5211,129 @@ def group_feature_class(df, feature_groups=['cell', 'cytoplasm', 'nucleus', 'pat
         , ignore_index=True)
     
     return df
+
+def delete_intermedeate_files(settings):
+    
+    path_orig = os.path.join(settings['src'], 'orig')
+    path_stack = os.path.join(settings['src'], 'stack')
+    merged_stack = os.path.join(settings['src'], 'merged')
+    path_norm_chan_stack = os.path.join(settings['src'], 'norm_channel_stack')
+    path_1 = os.path.join(settings['src'], '1')
+    path_2 = os.path.join(settings['src'], '2')
+    path_3 = os.path.join(settings['src'], '3')
+    path_4 = os.path.join(settings['src'], '4')
+    path_5 = os.path.join(settings['src'], '5')
+    path_6 = os.path.join(settings['src'], '6')
+    path_7 = os.path.join(settings['src'], '7')
+    path_8 = os.path.join(settings['src'], '8')
+    path_9 = os.path.join(settings['src'], '9')
+    path_10 = os.path.join(settings['src'], '10')
+    
+    paths = [path_stack, path_norm_chan_stack, path_1, path_2, path_3, path_4, path_5, path_6, path_7, path_8, path_9, path_10]
+    
+    merged_len = len(merged_stack)
+    stack_len = len(path_stack)
+    
+    if merged_len == stack_len and stack_len != 0:
+        if 'src' in settings:
+            if os.path.exists(settings['src']):
+                if os.path.exists(path_orig):
+                    for path in paths:
+                        if os.path.exists(path):
+                            try:
+                                shutil.rmtree(path)
+                                print(f"Deleted {path}")
+                            except OSError as e:
+                                print(f"{path} could not be deleted: {e}. Delete manually.")
+                else:
+                    print(f"{path_orig} does not exist.")
+            else:
+                print(f"{settings['src']} does not exist.")
+        else:
+            print("No 'src' key in settings dictionary.")
+        
+def filter_and_save_csv(input_csv, output_csv, column_name, upper_threshold, lower_threshold):
+    """
+    Reads a CSV into a DataFrame, filters rows based on a column for values > upper_threshold and < lower_threshold,
+    and saves the filtered DataFrame to a new CSV file.
+
+    Parameters:
+        input_csv (str): Path to the input CSV file.
+        output_csv (str): Path to save the filtered CSV file.
+        column_name (str): Column name to apply the filters on.
+        upper_threshold (float): Upper threshold for filtering (values greater than this are retained).
+        lower_threshold (float): Lower threshold for filtering (values less than this are retained).
+
+    Returns:
+        None
+    """
+    # Read the input CSV file into a DataFrame
+    df = pd.read_csv(input_csv)
+
+    # Filter rows based on the thresholds
+    filtered_df = df[(df[column_name] > upper_threshold) | (df[column_name] < lower_threshold)]
+
+    # Save the filtered DataFrame to a new CSV file
+    filtered_df.to_csv(output_csv, index=False)
+    display(filtered_df)
+
+    print(f"Filtered DataFrame saved to {output_csv}")
+    
+def extract_tar_bz2_files(folder_path):
+    """
+    Extracts all .tar.bz2 files in the given folder into subfolders with the same name as the tar file.
+    
+    Parameters:
+        folder_path (str): Path to the folder containing .tar.bz2 files.
+    """
+    if not os.path.isdir(folder_path):
+        raise ValueError(f"The provided path '{folder_path}' is not a valid folder.")
+    
+    # Iterate over files in the folder
+    for file_name in os.listdir(folder_path):
+        if file_name.endswith('.tar.bz2'):
+            file_path = os.path.join(folder_path, file_name)
+            extract_folder = os.path.join(folder_path, os.path.splitext(os.path.splitext(file_name)[0])[0])
+            
+            # Create the subfolder for extraction if it doesn't exist
+            os.makedirs(extract_folder, exist_ok=True)
+            
+            # Extract the tar.bz2 file
+            try:
+                with tarfile.open(file_path, 'r:bz2') as tar:
+                    tar.extractall(path=extract_folder)
+                print(f"Extracted: {file_name} -> {extract_folder}")
+            except Exception as e:
+                print(f"Failed to extract {file_name}: {e}")
+            
+            
+def calculate_shortest_distance(df, object1, object2):
+    """
+    Calculate the shortest edge-to-edge distance between two objects (e.g., pathogen and nucleus).
+    
+    Parameters:
+    - df: Pandas DataFrame containing measurements
+    - object1: String, name of the first object (e.g., "pathogen")
+    - object2: String, name of the second object (e.g., "nucleus")
+
+    Returns:
+    - df: Pandas DataFrame with a new column for shortest edge-to-edge distance.
+    """
+
+    # Compute centroid-to-centroid Euclidean distance
+    centroid_distance = np.sqrt(
+        (df[f'{object1}_channel_0_centroid_weighted-0'] - df[f'{object2}_channel_0_centroid_weighted-0'])**2 +
+        (df[f'{object1}_channel_0_centroid_weighted-1'] - df[f'{object2}_channel_0_centroid_weighted-1'])**2
+    )
+
+    # Estimate object radii using Feret diameters
+    object1_radius = df[f'{object1}_feret_diameter_max'] / 2
+    object2_radius = df[f'{object2}_feret_diameter_max'] / 2
+
+    # Compute shortest edge-to-edge distance
+    shortest_distance = centroid_distance - (object1_radius + object2_radius)
+
+    # Ensure distances are non-negative (overlapping objects should have distance 0)
+    df[f'{object1}_{object2}_shortest_distance'] = np.maximum(shortest_distance, 0)
+
+    return df

{spacr-0.3.81.dist-info → spacr-0.4.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: spacr
-Version: 0.3.81
+Version: 0.4.1
 Summary: Spatial phenotype analysis of crisp screens (SpaCr)
 Home-page: https://github.com/EinarOlafsson/spacr
 Author: Einar Birnir Olafsson
@@ -32,6 +32,7 @@ Requires-Dist: monai>=1.3.0
 Requires-Dist: captum<1.0,>=0.7.0
 Requires-Dist: seaborn<1.0,>=0.13.2
 Requires-Dist: matplotlib<4.0,>=3.8.3
+Requires-Dist: matplotlib-venn<2.0,>=1.1
 Requires-Dist: adjustText<2.0,>=1.2.0
 Requires-Dist: bottleneck<2.0,>=1.3.6
 Requires-Dist: numexpr<3.0,>=2.8.4

{spacr-0.3.81.dist-info → spacr-0.4.1.dist-info}/RECORD

@@ -1,4 +1,4 @@
-spacr/__init__.py,sha256=fvk5JfLpOqUA1W0yPcsVZnS9qbpXFOceFk09LKolVfw,1627
+spacr/__init__.py,sha256=iPlE-WRc1CjSNGPMbVvEKlJsTQlKjq29VjzvNUCpcOM,1401
 spacr/__main__.py,sha256=bkAJJD2kjIqOP-u1kLvct9jQQCeUXzlEjdgitwi1Lm8,75
 spacr/app_annotate.py,sha256=W9eLPa_LZIvXsXx_-0iDFEU938LBDvRy6prXo0qF4KQ,2533
 spacr/app_classify.py,sha256=urTP_wlZ58hSyM5a19slYlBxN0PdC-9-ga0hvq8CGWc,165
@@ -9,27 +9,28 @@ spacr/app_sequencing.py,sha256=DjG26jy4cpddnV8WOOAIiExtOe9MleVMY4MFa5uTo5w,157
 spacr/app_umap.py,sha256=ZWAmf_OsIKbYvolYuWPMYhdlVe-n2CADoJulAizMiEo,153
 spacr/cellpose.py,sha256=RBHMs2vwXcfkj0xqAULpALyzJYXddSRycgZSzmwI7v0,14755
 spacr/chat_bot.py,sha256=n3Fhqg3qofVXHmh3H9sUcmfYy9MmgRnr48663MVdY9E,1244
-spacr/core.py,sha256=3u2qKmPmTlswvE1uKTF4gi7KQ3sJBHV9No_ysgk7JCU,48487
-spacr/deep_spacr.py,sha256=V3diLyxX-0_F5UxhX_b94ROOvL9eoLvnoUmF3nMBqPQ,43250
+spacr/core.py,sha256=lKeqmsVrGQ8cPU_WkoNGNBWrk-gtR1RkRkwDdnJ0u64,48829
+spacr/deep_spacr.py,sha256=WN64EaQqF87JZg3Uan46t5Y28xsAGD2KMjr2ht6CyDs,54563
 spacr/gui.py,sha256=ARyn9Q_g8HoP-cXh1nzMLVFCKqthY4v2u9yORyaQqQE,8230
-spacr/gui_core.py,sha256=3S1S7rOtGfm5Z-Fb-jk1HxjsLZo6IboBxtXtsDMw1ME,46578
-spacr/gui_elements.py,sha256=fc9z-GP_PmtJy2lIsfi3fx5r5puTJ-zz8ntB9pC53Gc,138246
-spacr/gui_utils.py,sha256=u9RoIOWpAXFEOnUlLpMQZrc1pWSg6omZsJMIhJdRv_g,41211
-spacr/io.py,sha256=LF6lpphw7GSeuoHQijPykjKNF56wNTFEWFZuDQp3O6Q,145739
+spacr/gui_core.py,sha256=U0A7waKgWq_Es9fMwcZbXUZYGzCqt2bgfY3HbxiFXnw,47466
+spacr/gui_elements.py,sha256=HmITDncklKwtdFhxLhtYXOwndsRfgwWIPVi83VlXHB4,146419
+spacr/gui_utils.py,sha256=0rDF23BUGcmjSJvfCiLoxhlGJdHkio1jTxyCzrMXr-g,41211
+spacr/io.py,sha256=oqJwDJWksVdWE0bRAwytTOsjlL0o-J9lr_pQaw2cQ4Y,138288
 spacr/logger.py,sha256=lJhTqt-_wfAunCPl93xE65Wr9Y1oIHJWaZMjunHUeIw,1538
-spacr/measure.py,sha256=2lK-ZcTxLM-MpXV1oZnucRD9iz5aprwahRKw9IEqshg,55085
+spacr/measure.py,sha256=jmOnLBudq3TuY723Cfo1EJBn67P6rlEvL6I-2FSkUgI,55315
 spacr/mediar.py,sha256=FwLvbLQW5LQzPgvJZG8Lw7GniA2vbZx6Jv6vIKu7I5c,14743
-spacr/ml.py,sha256=x19S8OsR5omb8e6MU9I99Nz95J_QvM5siyk-zaAU3p8,82866
+spacr/ml.py,sha256=MrIAtUUxMOibWVL1SjCUnYlizawCp3l3SeY4Y9yEsPw,97251
 spacr/openai.py,sha256=5vBZ3Jl2llYcW3oaTEXgdyCB2aJujMUIO5K038z7w_A,1246
-spacr/plot.py,sha256=gXC7y3uT4sx8KRODeSFWQG_A1CylsuJ5B7HYe_un6so,165177
+spacr/plot.py,sha256=Q5TbsR2NUWhA7z4HyF_2_FAEBFSNMU-G3UNDbRzW6mM,169485
 spacr/sequencing.py,sha256=ClUfwPPK6rNUbUuiEkzcwakzVyDKKUMv9ricrxT8qQY,25227
-spacr/settings.py,sha256=3bgBWBotIO7TZx2nh6JoEaqHNmwbChAiz1gW4xmURQs,81788
+spacr/settings.py,sha256=fEk-9LSSvV1wGsn6xTaJWY7wF7_u8Fc-S1DaDHqZU3I,83997
 spacr/sim.py,sha256=1xKhXimNU3ukzIw-3l9cF3Znc_brW8h20yv8fSTzvss,71173
+spacr/sp_stats.py,sha256=mbhwsyIqt5upsSD346qGjdCw7CFBa0tIS7zHU9e0jNI,9536
 spacr/stats.py,sha256=mbhwsyIqt5upsSD346qGjdCw7CFBa0tIS7zHU9e0jNI,9536
-spacr/submodules.py,sha256=SK8YEs850LAx30YAiwap7ecLpp1_p-bci6H-Or0GLoA,55500
+spacr/submodules.py,sha256=mb2g0igUTws7y6xW1zIJw1E7eQyxsjEj5mk2Z-Qd8uw,67629
 spacr/timelapse.py,sha256=KGfG4L4-QnFfgbF7L6C5wL_3gd_rqr05Foje6RsoTBg,39603
-spacr/toxo.py,sha256=z2nT5aAze3NUIlwnBQcnkARihDwoPfqOgQIVoUluyK0,25087
-spacr/utils.py,sha256=jjZIqzJKl9nnWgj2eJiLFT27gED8hO6rAEsJLlimm-E,222298
+spacr/toxo.py,sha256=TmuhejSIPLBvsgeblsUgSvBFCR1gOkApyTKidooJ5Us,26044
+spacr/utils.py,sha256=of2t5Tq_RKdJ1QRDo4nJ3oEVev_6s2Oko3-lBxl4ScU,226293
 spacr/version.py,sha256=axH5tnGwtgSnJHb5IDhiu4Zjk5GhLyAEDRe-rnaoFOA,409
 spacr/resources/MEDIAR/.gitignore,sha256=Ff1q9Nme14JUd-4Q3jZ65aeQ5X4uttptssVDgBVHYo8,152
 spacr/resources/MEDIAR/LICENSE,sha256=yEj_TRDLUfDpHDNM0StALXIt6mLqSgaV2hcCwa6_TcY,1065
@@ -152,9 +153,9 @@ spacr/resources/icons/umap.png,sha256=dOLF3DeLYy9k0nkUybiZMe1wzHQwLJFRmgccppw-8b
 spacr/resources/images/plate1_E01_T0001F001L01A01Z01C02.tif,sha256=Tl0ZUfZ_AYAbu0up_nO0tPRtF1BxXhWQ3T3pURBCCRo,7958528
 spacr/resources/images/plate1_E01_T0001F001L01A02Z01C01.tif,sha256=m8N-V71rA1TT4dFlENNg8s0Q0YEXXs8slIn7yObmZJQ,7958528
 spacr/resources/images/plate1_E01_T0001F001L01A03Z01C03.tif,sha256=Pbhk7xn-KUP6RSIhJsxQcrHFImBm3GEpLkzx7WOc-5M,7958528
-spacr-0.3.81.dist-info/LICENSE,sha256=SR-2MeGc6SCM1UORJYyarSWY_A-JaOMFDj7ReSs9tRM,1083
-spacr-0.3.81.dist-info/METADATA,sha256=TneYaiWWfNHUuqpSQpRAEzWzYia6KIksYDzNVDcsimw,6032
-spacr-0.3.81.dist-info/WHEEL,sha256=HiCZjzuy6Dw0hdX5R3LCFPDmFS4BWl8H-8W39XfmgX4,91
-spacr-0.3.81.dist-info/entry_points.txt,sha256=BMC0ql9aNNpv8lUZ8sgDLQMsqaVnX5L535gEhKUP5ho,296
-spacr-0.3.81.dist-info/top_level.txt,sha256=GJPU8FgwRXGzKeut6JopsSRY2R8T3i9lDgya42tLInY,6
-spacr-0.3.81.dist-info/RECORD,,
+spacr-0.4.1.dist-info/LICENSE,sha256=SR-2MeGc6SCM1UORJYyarSWY_A-JaOMFDj7ReSs9tRM,1083
+spacr-0.4.1.dist-info/METADATA,sha256=ibaavHM7Kuo-dVo-SkOgjKfL0D8mG8ouYCHEzTqsbwg,6072
+spacr-0.4.1.dist-info/WHEEL,sha256=HiCZjzuy6Dw0hdX5R3LCFPDmFS4BWl8H-8W39XfmgX4,91
+spacr-0.4.1.dist-info/entry_points.txt,sha256=BMC0ql9aNNpv8lUZ8sgDLQMsqaVnX5L535gEhKUP5ho,296
+spacr-0.4.1.dist-info/top_level.txt,sha256=GJPU8FgwRXGzKeut6JopsSRY2R8T3i9lDgya42tLInY,6
+spacr-0.4.1.dist-info/RECORD,,