-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathshared_cells_distance_calculations.py
131 lines (119 loc) · 6.55 KB
/
shared_cells_distance_calculations.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
from brainrender.scene import Scene
from brainrender.actors import Points
from bg_atlasapi.bg_atlas import BrainGlobeAtlas
from brainrender import settings
from brainrender.actors import Cylinder
import bg_space as bg
import numpy as np
from rich import print
from myterial import orange
from pathlib import Path
import pandas as pd
import json
from bg_atlasapi import show_atlases
from bg_atlasapi.bg_atlas import BrainGlobeAtlas
import UpdateME
from UpdateME import cellfinder_output_path, mouse_id, brain_regions_to_evalutate, allen_mouse_10um, estim_tip_coordinates, extra_brain_region_acryonm, estim_shank_radius_um, estim_tip_radius_um, estim_propigation_radius_um
import os
import brainrender
from vedo import Spheres, Sphere
from vedo import Points as vPoints
import matplotlib.pyplot as plt
from scipy.spatial import distance
import seaborn as sns
brainrender.SHADER_STYLE = "cartoon"
import cellfinder_backend
from cellfinder_backend import analyze_data_cellfinder
import pickle
def shared_cell_distance_calculations_histograms(brainrender_folder_path,mouse_id,estim_tip_coordinates):
# Create new folder in your cellfinder output folder
mouseid_estim_tip_coordinates_folder_path = brainrender_folder_path + '/' + str(mouse_id) + '_' + str(estim_tip_coordinates) # create the path for the new folder
histogram_folder_path = mouseid_estim_tip_coordinates_folder_path + '/' + str(mouse_id) + '_' + str(estim_tip_coordinates) + '_histograms/overlapping_gfp_tdTomato'
npy_folder_path = histogram_folder_path + '/' + "overlapping_distance_arrays"
if not os.path.exists(histogram_folder_path):
os.makedirs(histogram_folder_path)
os.makedirs(npy_folder_path)
print('Histogram folder has been created at')
print(f'{ histogram_folder_path}')
print('')
print('Making histograms...')
print('')
histogram_save_path = histogram_folder_path
# npy_save_path = histogram_save_path + '/' + 'distance_arrays'
# Calaculate 3d-space distances for cells relative to estim_tip_coordinates. and save out histograms of those distances
# Path to cellfinder_output points.npy file
cells_path = brainrender_folder_path + '/' + str(mouse_id) + '_' + str(estim_tip_coordinates) + '/shared_cells_' + mouse_id + '_' +str(estim_tip_coordinates) + ".npy"
points = np.load(cells_path)
hist_save_path = histogram_save_path +'/'+ str(mouse_id)+ '_' + str(estim_tip_coordinates)
# Subtract the reference point from each cell coordinate
displacement = points - estim_tip_coordinates
# Calculate the magnitude of the displacement vectors
euclidean_distances = np.linalg.norm(displacement, axis=1)
plt.hist(euclidean_distances, bins=500)
# Add labels and a title
plt.xlabel('Euclidean Distance (um)')
plt.ylabel('number of cells')
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato Euclidean distances from Estim tip')
# Display & save the histogram
hist_save_path = histogram_save_path +'/'+ str(mouse_id)+ '_' + str(estim_tip_coordinates)
plt.savefig(hist_save_path + '_euclidean_distances.png')
np.save(npy_folder_path +'/'+ 'euclidean_distances.npy', euclidean_distances)
plt.close()
sns.histplot(euclidean_distances, kde = True)
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato Euclidean distances from Estim tip')
plt.xlabel('Euclidean Distance (um)')
plt.ylabel('number of cells')
plt.savefig(hist_save_path + '_sns_euclidean_distances.png')
plt.close()
# manhattan distances calculations: sum of the absolute differences of their coordinates
manhattan_distances = np.sum(np.abs(points - estim_tip_coordinates), axis=1)
manhattan_distances = distance.cdist(points, [estim_tip_coordinates], metric='cityblock')
# Plot manhattan distances histogram
plt.hist(manhattan_distances, bins = 100)
# Add labels
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato Manhattan distances from Estim tip')
plt.xlabel('Manhattan distance (um)')
plt.ylabel('number of cells')
# Save the histogram
plt.savefig(hist_save_path + '_manhattan_distances.png')
np.save(npy_folder_path+'/' + 'manhattan_distances.npy', manhattan_distances)
plt.close()
# create subplot with smooth line overlay
sns.histplot(manhattan_distances, kde = True)
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato Manhattan distances from Estim tip')
plt.xlabel('Manhattan distance (um)')
plt.ylabel('number of cells')
plt.savefig(hist_save_path + '_sns_manhattan_distances.png')
plt.close()
p = 3
minkowski_distances = distance.cdist(points, [estim_tip_coordinates], 'minkowski', p=p)
plt.hist(minkowski_distances, bins=100)
plt.xlabel('Minkowski Distance (um) (p={})'.format(p))
plt.ylabel('Cell Count')
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato Minkowski Distances from Estim tip')
plt.savefig(hist_save_path + '_minkowski_distances.png')
np.save(npy_folder_path +'/'+ 'minkowski_distances.npy', minkowski_distances)
plt.close()
sns.histplot(minkowski_distances, kde = True)
plt.xlabel('Minkowski Distance (um) (p={})'.format(p))
plt.ylabel('Cell Count')
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato Minkowski Distances from Estim tip')
plt.savefig(hist_save_path + '_sns_minkowski_distances.png')
plt.close()
chebyshev_distances = distance.cdist(points, [estim_tip_coordinates], 'chebyshev')
plt.hist(chebyshev_distances, bins=100)
plt.xlabel('Chebyshev Distance (um)')
plt.ylabel('Cell count')
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato, Chebyshev Distances from Estim tip')
plt.savefig(hist_save_path + '_chebyshev_distances.png')
np.save(npy_folder_path +'/'+ 'chebyshev_distances.npy', chebyshev_distances)
plt.close()
sns.histplot(minkowski_distances, kde = True)
plt.xlabel('Chebyshev Distance (um)')
plt.ylabel('Cell count')
plt.title(str(mouse_id) + ' Overlapping gfp & tdTomato, Chebyshev Distances from Estim tip')
plt.savefig(hist_save_path + '_sns_chebyshev_distances.png')
plt.close()
else:
print('Skipping histogram creation, because ' + 'histograms for ' + str(mouse_id) + " with estim_tip_coordinates " + str(estim_tip_coordinates) + ' have already been created')
print('')