Measuring the radial intensity distribution of equally spaced regions of circular-like shapes
Hi I'm currently stuck on a problem that I'm hoping someone may be able to help me with. I've tried looking around for a solution but haven't had too much luck. I am a biologist starting to integrate automated image-analysis into my workflow so forgive me if this isn't the right place to ask this question. I have performed segmentation with skimage to identify objects within a set of images (example below) obtained via microscopy as seen below:
What I am wanting to do now however is take each one of these objects and measure their individual radial profiles by creating an n-number of equally spaced rings within the object. From this I want to extract the mean and total intensity from each of those rings in order to quantify the variation of radial intensity (example below). The biggest issue I'm having is actually finding a way to re-segment these objects into equally spaced rings with skimage. I have tried creating multiple rings from the object-center however as the object-bounds are not perfect circles this does not achieve what we want.
FYI - This is also a function available using the software CellProfiler which I have used in the past to accomplish this (example of function can be found here) however I am attempting to recreate this functionality independently of CellProfiler so we can run our analysis more efficiently. Any guidance or even nudges in the right direction would be greatly appreciated.
Let's say we have the coordinates of the object-center (x0, y0) along with the set of M coordinates of the object-edge H={(hx1, hy1), (hx2, hy2), ..., (hxM, hyM)}.
In this context, our aim is to find the coordinates of N rings or closed curves that equally segment the object. All the rings have the same center as the object located at (x0, y0).
I suggest that we can do it through calculating the coordinates of N evenly spaced sections segmenting all the lines lying between the object-center and the object edge as L={ {(x0, y0), (hx1, hy1)}, {(x0, y0), (hx2, hy2)}, ..., {(x0, y0), (hxM, hyM)} }.
E.g, utilizing this formula https://www.geeksforgeeks.org/section-formula-point-divides-line-given-ratio/
- Activating Anaconda Environment in VsCode
- How to implement the Softmax function in Python?
- Dataframe - Select the minimum set of rows to cover all possible values of each columns
- Non-equi join in polars
- Tkinter and matplotlib - NavigationToolbar2Tk - cursor position precision
- Media Image is not loaded in the output of Django Project
- selenium upload part not working on windows. Path related?
- How to encode UTF8 filename for HTTP headers? (Python, Django)
- Split columns in Csv file
- In the Django admin site, how do I change the display format of time fields?
- Condition in pytest fixture depending on some test-suite results
- How to make a class JSON serializable
- Apache Airflow unit and integration test
- TypeError: 'MultiPolygon' object is not iterable with Plotly Choropleth
- How do I get the current time in milliseconds in Python?
- Python: Decoding base64 encoded strings
- How to fix the error when try to plot data with pandas?
- Python Image - Finding largest branch from image skeleton
- Enumerate is Returning A None Type
- Python finite difference functions?
- How to extract hex color codes from a colormap
- Regex: Get "London, UK" from "Evolution Recruitment (Agency) (London, UK)"
- Visual Studio Code - How to add multiple paths to python path?
- What are the Spark transformations that causes a Shuffle?
- How to find the second lowest lists into a nested list by their second value?
- python-docx: Parse a table to Panda Dataframe
- Deposit and Withdraw on Account
- Why does this code snippet give the error list object is not callable?
- How to find out if argparse argument has been actually specified on command line?
- Print current call stack from a method in code