mcDETECT 1.0.3__py3-none-any.whl → 1.0.6__py3-none-any.whl

mcDETECT/model.py

@@ -304,7 +304,7 @@ class mcDETECT:
         return sphere
     
     
-    # [MAIN] dataframe, metadata of identified synapses
+    # [MAIN] dataframe, synapse metadata
     def detect(self):
         
         _, data_low, data_high = self.dbscan()
@@ -319,7 +319,7 @@ class mcDETECT:
             return self.nc_filter(sphere_low, sphere_high)
     
     
-    # [MAIN] anndata, spatial transcriptome profile of identified synapses
+    # [MAIN] anndata, synapse spatial transcriptome profile
     def profile(self, synapse, genes = None, print_itr = False):
         
         if genes is None:
@@ -348,4 +348,86 @@ class mcDETECT:
         adata.var['genes'] = genes
         adata.var_names = genes
         adata.var_keys = genes
-        return adata
+        return adata
+    
+    
+    # [MAIN] anndata, spot-level gene expression
+    def spot_expression(self, grid_len, genes = None):
+        
+        if genes is None:
+            genes = list(np.unique(self.transcripts['target']))
+            transcripts = self.transcripts
+        else:
+            transcripts = self.transcripts[self.transcripts['target'].isin(genes)]
+        
+        # construct bins
+        x_bins, y_bins = self.construct_grid(grid_len = None)
+        
+        # initialize data
+        X = np.zeros((len(genes), len(x_bins) * len(y_bins)))
+        global_x, global_y = [], []
+        
+        # coordinates
+        for i in list(x_bins):
+            center_x = i + 0.5 * grid_len
+            for j in list(y_bins):
+                center_y = j + 0.5 * grid_len
+                global_x.append(center_x)
+                global_y.append(center_y)
+        
+        # count matrix
+        for k_idx, k in enumerate(genes):
+            target_gene = transcripts[transcripts['target'] == k]
+            count_gene, _, _ = np.histogram2d(target_gene['global_x'], target_gene['global_y'], bins = [x_bins, y_bins])
+            X[k_idx, :] = count_gene.flatten()
+            if k_idx % 50 == 0:
+                print("{} out of {} genes profiled!".format(k_idx + 1, len(genes)))
+        
+        # spot id
+        spot_id = []
+        for i in range(len(global_x)):
+            id = 'spot_' + str(i)
+            spot_id.append(id)
+        
+        # assemble data
+        adata = anndata.AnnData(X = np.transpose(X))
+        adata.obs['spot_id'] = spot_id
+        adata.obs['global_x'] = global_x
+        adata.obs['global_y'] = global_y
+        adata.var['genes'] = genes
+        adata.var_names = genes
+        adata.var_keys = genes
+        return adata
+    
+    
+    # [MAIN] anndata, spot-level synapse metadata
+    def spot_synapse(self, synapse, spot):
+        
+        x_grid, y_grid = list(np.unique(spot.obs['global_x'])), list(np.unique(spot.obs['global_y']))
+        diameter = x_grid[1] - x_grid[0]
+        
+        indicator, synapse_count, synapse_radius, synapse_size, synapse_score = [], [], [], [], []
+        for i in x_grid:
+            x_min_temp = i
+            x_max_temp = i + diameter
+            for j in y_grid:
+                y_min_temp = j
+                y_max_temp = j + diameter
+                syn_temp = synapse[(synapse['sphere_x'] > x_min_temp) & (synapse['sphere_x'] < x_max_temp) & (synapse['sphere_y'] > y_min_temp) & (synapse['sphere_y'] < y_max_temp)]
+                indicator.append(int(syn_temp.shape[0] > 0))
+                synapse_count.append(syn_temp.shape[0])
+                if syn_temp.shape[0] == 0:
+                    synapse_radius.append(0)
+                    synapse_size.append(0)
+                    synapse_score.append(0)
+                else:
+                    synapse_radius.append(np.nanmean(syn_temp['sphere_r']))
+                    synapse_size.append(np.nanmean(syn_temp['size']))
+                    synapse_score.append(np.nanmean(syn_temp['in_nucleus']))
+        
+        spot.obs['indicator'] = indicator
+        spot.obs['syn_count'] = synapse_count
+        spot.obs['syn_radius'] = synapse_radius
+        spot.obs['syn_size'] = synapse_size
+        spot.obs['syn_score'] = synapse_score
+        return spot

{mcdetect-1.0.3.dist-info → mcdetect-1.0.6.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: mcDETECT
-Version: 1.0.3
+Version: 1.0.6
 Summary: mcDETECT: Decoding 3D Spatial Synaptic Transcriptomes with Subcellular-Resolution Spatial Transcriptomics
 Home-page: https://github.com/chen-yang-yuan/mcDETECT
 Author: Chenyang Yuan

mcdetect-1.0.6.dist-info/RECORD

@@ -0,0 +1,7 @@
+mcDETECT/__init__.py,sha256=smjxcWMgqIGgbQ47Hx0TRXtiJLHcuXizBk9DTfAZp3s,27
+mcDETECT/model.py,sha256=RpZ-x3T2E-NVynHc_pAitYxR-fT-FfugcZaSie_xXoM,21178
+mcdetect-1.0.6.dist-info/LICENSE,sha256=uxq-shEWOGTIGVnQLmpElILmfCkuUhFZRAMnZUiKvtg,1070
+mcdetect-1.0.6.dist-info/METADATA,sha256=aUMpnnSnn95yNOfaHtLmF3fPI9I1lJMWXqB18xwZsjM,2819
+mcdetect-1.0.6.dist-info/WHEEL,sha256=jB7zZ3N9hIM9adW7qlTAyycLYW9npaWKLRzaoVcLKcM,91
+mcdetect-1.0.6.dist-info/top_level.txt,sha256=WwzBojt5U-T2hZ8llO6XgpM9OFIBkWQQldQKu19O8EY,9
+mcdetect-1.0.6.dist-info/RECORD,,

mcdetect-1.0.3.dist-info/RECORD

@@ -1,7 +0,0 @@
-mcDETECT/__init__.py,sha256=smjxcWMgqIGgbQ47Hx0TRXtiJLHcuXizBk9DTfAZp3s,27
-mcDETECT/model.py,sha256=d10VQzt2nycqMlfhC5tvGThDSYAbzKdZR_kmtXBo_MM,17926
-mcdetect-1.0.3.dist-info/LICENSE,sha256=uxq-shEWOGTIGVnQLmpElILmfCkuUhFZRAMnZUiKvtg,1070
-mcdetect-1.0.3.dist-info/METADATA,sha256=6v5umeimWywaDJROtXpVvjL0EiuWIIVraxP4xXEaKK0,2819
-mcdetect-1.0.3.dist-info/WHEEL,sha256=jB7zZ3N9hIM9adW7qlTAyycLYW9npaWKLRzaoVcLKcM,91
-mcdetect-1.0.3.dist-info/top_level.txt,sha256=WwzBojt5U-T2hZ8llO6XgpM9OFIBkWQQldQKu19O8EY,9
-mcdetect-1.0.3.dist-info/RECORD,,