Kernel: Python 3
In [ ]:
Quantifying spinal cord regeneration in zebrafish
We want to quantify the amount of fluorescent cells in a wounded zebrafish embryo spinal column:
The key steps are:
estimating the position and width of the cord
estimating the average fluorescence along the length of the cord
In [1]:
SciPy to estimate coordinates
First, we get just the top and bottom rows of pixels, and use a 1D gaussian filter to smooth the signal.
In [2]:
<matplotlib.text.Text at 0x1096a3090>
With smooth curves, we can get the mode (the position of the center) and width of the signal.
In [3]:
(421, 181, 739, 142)
scikit-image to trace the profile
Now, use measure.profile_line
to trace from (0, top_mode
) to (-1, bottom_mode
).
In [8]:
Finally, plot the trace.
In [9]:
<matplotlib.text.Text at 0x10a9fa0d0>
From this trace, we can compute various summary statistics (e.g. min/max, gap width, slope, etc), and plot these over time as the wound recovers.