| Name | cp_measure JSON |
| Version |
0.1.2
JSON |
| download |
| home_page | None |
| Summary | cp-measure implements CellProfiler Measurements. You can integrate them to your workflow or use them from CellProfiler 5 once they are integrated there. |
| upload_time | 2024-09-27 19:40:27 |
| maintainer | None |
| docs_url | None |
| author | Alan Munoz |
| requires_python | <3.12,>=3.9 |
| license | None |
| keywords |
|
| VCS |
|
| bugtrack_url |
|
| requirements |
No requirements were recorded.
|
| Travis-CI |
No Travis.
|
| coveralls test coverage |
No coveralls.
|
#+TITLE: Cell Profiler measurements
Do you need to use [[https://github.com/CellProfiler][CellProfiler]] features, but you want to do it in a programmatic way? Look no more, this package was developed by and for the click-a-phobic scientists.
* Quick overview
** Installation
#+begin_src bash
pip install cp-measure
#+end_src
*** Poetry
If you want a development environment.
#+begin_src bash
git clone git@github.com:afermg/cp_measure.git
cd cp_measure
poetry install
#+end_src
** Usage
Commonly users want to calculate all the features. There are four type of measurements, based on their inputs:
- Type 1: 1 image + 1 set of masks (e.g., intensity)
- Type 2: 2 images + 1 set of masks (e.g., colocalization)
- Type 3: 2 sets of masks (e.g., number of neighbors)
- Type 4: 1 image + 2 sets of masks (e.g., skeleton)
This shows the simplest way to use the first set (1 image, 1 mask set), which currently follows the style of the scikit-image syule (1 image, 1 matrix with non-overlapping labels).
#+begin_src python
import numpy as np
from cp_measure.bulk import get_fast_measurements
measurements = get_fast_measurements()
print(measurements.keys())
# dict_keys(['radial_distribution', 'radial_zernikes', 'intensity', 'sizeshape', 'zernike', 'ferret', 'texture', 'granularity'])
size = 200
rng = np.random.default_rng(42)
pixels = rng.integers(low=0, high=10, size=(size, size))
masks = np.zeros_like(pixels)
masks[5:-6, 5:-6] = 1
results = {}
for name, v in measurements.items():
results = {**results, **v(masks, pixels)}
"""
{'RadialDistribution_FracAtD_1of4': array([0.0477544]),
'RadialDistribution_MeanFrac_1of4': array([0.98888986]),
'RadialDistribution_RadialCV_1of4': array([0.04705031]),
...
'RadialDistribution_ZernikeMagnitude_0_0': array([4.49027183]),
'RadialDistribution_ZernikePhase_0_0': array([1.57079633]),
..
'RadialDistribution_ZernikeMagnitude_9_9': array([0.00250042]),
'RadialDistribution_ZernikePhase_9_9': array([-2.81002186]),
'Intensity_IntegratedIntensity': array([160425.]),
'Intensity_MeanIntensity': array([4.49105568]),
'Intensity_StdIntensity': array([2.8698785]),
'Intensity_MinIntensity': array([0.]),
'Intensity_MaxIntensity': array([9.]),
'Intensity_IntegratedIntensityEdge': array([3253.]),
'Intensity_MeanIntensityEdge': array([4.32579787]),
'Intensity_StdIntensityEdge': array([2.87255435]),
'Intensity_MinIntensityEdge': array([0.]),
'Intensity_MaxIntensityEdge': array([9.]),
'Intensity_MassDisplacement': array([0.36112335]),
'Intensity_LowerQuartileIntensity': array([2.]),
'Intensity_MedianIntensity': array([4.]),
'Intensity_MADIntensity': array([2.]),
'Intensity_UpperQuartileIntensity': array([7.]),
'Location_CenterMassIntensity_X': array([98.79390369]),
'Location_CenterMassIntensity_Y': array([99.29653732]),
'Location_CenterMassIntensity_Z': array([0.]),
'Location_MaxIntensity_X': array([5.]),
'Location_MaxIntensity_Y': array([5.]),
'Location_MaxIntensity_Z': array([0.]),
'SpatialMoment_0_0': array([35721.]),
'SpatialMoment_0_1': array([3357774.]),
...
'SpatialMoment_2_3': array([7.04675332e+14]),
'CentralMoment_0_0': array([35721.]),
'CentralMoment_0_1': array([0.]),
...
'CentralMoment_2_3': array([0.]),
'NormalizedMoment_0_0': array([nan]),
'NormalizedMoment_0_1': array([nan]),
'NormalizedMoment_0_2': array([0.083331]),
...
'NormalizedMoment_3_3': array([0.]),
'HuMoment_0': array([0.166662]),
...
'InertiaTensor_0_0': array([2976.66666667]),
'InertiaTensor_0_1': array([-0.]),
'InertiaTensor_1_0': array([-0.]),
'InertiaTensor_1_1': array([2976.66666667]),
'InertiaTensorEigenvalues_0': array([2976.66666667]),
'InertiaTensorEigenvalues_1': array([2976.66666667]),
'Zernike_0_0': array([0.64333829]),
...
'Zernike_9_9': array([3.9043712e-19]),
'MinFeretDiameter': array([188.]),
'MaxFeretDiameter': array([265.87214973]),
'AngularSecondMoment_3_00_256': array([0.01237016]),
'Contrast_3_00_256': array([13.32731615]),
'Correlation_3_00_256': array([0.00052763]),
'Variance_3_00_256': array([6.66717588]),
'InverseDifferenceMoment_3_00_256': array([0.27445799]),
'SumAverage_3_00_256': array([9.97394114]),
'SumVariance_3_00_256': array([13.34138738]),
'SumEntropy_3_00_256': array([3.87614681]),
'Entropy_3_00_256': array([6.33842161]),
'DifferenceVariance_3_00_256': array([0.00372448]),
'DifferenceEntropy_3_00_256': array([2.98563863]),
'InfoMeas1_3_00_256': array([-0.00035948]),
'InfoMeas2_3_00_256': array([0.04771106]),
'AngularSecondMoment_3_01_256': array([0.01237633]),
'Contrast_3_01_256': array([13.34467007]),
'Correlation_3_01_256': array([0.00074894]),
'Variance_3_01_256': array([6.67733598]),
'InverseDifferenceMoment_3_01_256': array([0.27467073]),
'SumAverage_3_01_256': array([9.96441954]),
'SumVariance_3_01_256': array([13.36467386]),
'SumEntropy_3_01_256': array([3.87774892]),
'Entropy_3_01_256': array([6.33805833]),
'DifferenceVariance_3_01_256': array([0.00361481]),
'DifferenceEntropy_3_01_256': array([2.9890093]),
'InfoMeas1_3_01_256': array([-0.00041841]),
'InfoMeas2_3_01_256': array([0.05146776]),
'AngularSecondMoment_3_02_256': array([0.01237207]),
'Contrast_3_02_256': array([13.36134306]),
'Correlation_3_02_256': array([-0.0008637]),
'Variance_3_02_256': array([6.67490643]),
'InverseDifferenceMoment_3_02_256': array([0.27233695]),
'SumAverage_3_02_256': array([9.95842397]),
'SumVariance_3_02_256': array([13.33828266]),
'SumEntropy_3_02_256': array([3.87581509]),
'Entropy_3_02_256': array([6.33830874]),
'DifferenceVariance_3_02_256': array([0.00369281]),
'DifferenceEntropy_3_02_256': array([2.98609235]),
'InfoMeas1_3_02_256': array([-0.00037672]),
'InfoMeas2_3_02_256': array([0.04884048]),
'AngularSecondMoment_3_03_256': array([0.0123648]),
'Contrast_3_03_256': array([13.28753834]),
'Correlation_3_03_256': array([0.00346848]),
'Variance_3_03_256': array([6.66689316]),
'InverseDifferenceMoment_3_03_256': array([0.27439131]),
'SumAverage_3_03_256': array([9.97160996]),
'SumVariance_3_03_256': array([13.38003429]),
'SumEntropy_3_03_256': array([3.87790345]),
'Entropy_3_03_256': array([6.33872973]),
'DifferenceVariance_3_03_256': array([0.0036703]),
'DifferenceEntropy_3_03_256': array([2.98626744]),
'InfoMeas1_3_03_256': array([-0.0002397]),
'InfoMeas2_3_03_256': array([0.03896704]),
'Granularity_1': array([23.96446938]),
...
'Granularity_16': array([100.])}
"""
#+end_src
*** Call specific measurements
If you need a specific measurement/feature you can just import it. Note that measurements come in sets, so you have to fetch the one that you specifically require from the resultant dictionary. Any available measurement can be found using code as follows:
#+begin_src python
import numpy as np
from cp_measure.minimal.measureobjectsizeshape import get_sizeshape
mask = np.zeros((50, 50))
mask[5:-6, 5:-6] = 1
get_sizeshape(mask, None) # pixels, the second argument, is not necessary for this measurement
#+end_src
The other available functions are as follows:
- measureobjectintensitydistribution.get_radial_zernikes,
- measureobjectintensity.get_intensity,
- measureobjectsizeshape.get_zernike,
- measureobjectsizeshape.get_ferret,
- measuregranularity.get_granularity,
- measuretexture.get_texture,
And for Type 2 functions:
- measurecolocalization.get_correlation_pearson
- measurecolocalization.get_correlation_manders_fold
- measurecolocalization.get_correlation_rwc
- measurecolocalization.get_correlation_costes
- measurecolocalization.get_correlation_overlap
* Pending measurements
You can follow progress [[https://docs.google.com/spreadsheets/d/1_7jQ8EjPwOr2MUnO5Tw56iu4Y0udAzCJEny-LQMgRGE/edit?usp=sharing][here]].
*** Done
- Type 1 and 2 measurements in sklearn style (process multiple masks per image)
*** Pending
- Type 3 and 4 measurements
*** Additional notes
The Image-wide functions will not be implemented directly, as they were originally implemented as additional functions to the Object (mask) functions. We will adjust the existing functions assume that an image-wide measurement is the same as measuring an object with the same size as the intensity image.
* Additional notes
- This is not optimised for efficiency (yet). We aim to reproduce the 'vanilla' results of CellProfiler with minimal code changes. Optimisations will be implemented once we come up with a standard interface for functionally-focused CellProfiler components.
- The functions exposed perform minimal checks. They will fail if provided with empty masks. Not all functions will fail if provided with masks only.
Raw data
{
"_id": null,
"home_page": null,
"name": "cp_measure",
"maintainer": null,
"docs_url": null,
"requires_python": "<3.12,>=3.9",
"maintainer_email": null,
"keywords": null,
"author": "Alan Munoz",
"author_email": null,
"download_url": "https://files.pythonhosted.org/packages/3e/34/dd31db3434d55da4eaa30184d03ae21b3948a15a8b062f6ee7d19d014df9/cp_measure-0.1.2.tar.gz",
"platform": null,
"description": "#+TITLE: Cell Profiler measurements\n\nDo you need to use [[https://github.com/CellProfiler][CellProfiler]] features, but you want to do it in a programmatic way? Look no more, this package was developed by and for the click-a-phobic scientists.\n\n* Quick overview\n** Installation\n#+begin_src bash\n pip install cp-measure\n#+end_src\n*** Poetry\nIf you want a development environment.\n#+begin_src bash\n git clone git@github.com:afermg/cp_measure.git\n cd cp_measure\n poetry install \n#+end_src\n\n** Usage\nCommonly users want to calculate all the features. There are four type of measurements, based on their inputs:\n- Type 1: 1 image + 1 set of masks (e.g., intensity)\n- Type 2: 2 images + 1 set of masks (e.g., colocalization)\n- Type 3: 2 sets of masks (e.g., number of neighbors)\n- Type 4: 1 image + 2 sets of masks (e.g., skeleton)\n \nThis shows the simplest way to use the first set (1 image, 1 mask set), which currently follows the style of the scikit-image syule (1 image, 1 matrix with non-overlapping labels).\n#+begin_src python\nimport numpy as np\n\nfrom cp_measure.bulk import get_fast_measurements\n\nmeasurements = get_fast_measurements()\nprint(measurements.keys())\n# dict_keys(['radial_distribution', 'radial_zernikes', 'intensity', 'sizeshape', 'zernike', 'ferret', 'texture', 'granularity'])\n\nsize = 200\nrng = np.random.default_rng(42)\npixels = rng.integers(low=0, high=10, size=(size, size))\n\nmasks = np.zeros_like(pixels)\nmasks[5:-6, 5:-6] = 1\n\nresults = {}\nfor name, v in measurements.items():\n results = {**results, **v(masks, pixels)}\n\n\"\"\"\n{'RadialDistribution_FracAtD_1of4': array([0.0477544]),\n 'RadialDistribution_MeanFrac_1of4': array([0.98888986]),\n 'RadialDistribution_RadialCV_1of4': array([0.04705031]),\n ...\n 'RadialDistribution_ZernikeMagnitude_0_0': array([4.49027183]),\n 'RadialDistribution_ZernikePhase_0_0': array([1.57079633]),\n ..\n 'RadialDistribution_ZernikeMagnitude_9_9': array([0.00250042]),\n 'RadialDistribution_ZernikePhase_9_9': array([-2.81002186]),\n 'Intensity_IntegratedIntensity': array([160425.]),\n 'Intensity_MeanIntensity': array([4.49105568]),\n 'Intensity_StdIntensity': array([2.8698785]),\n 'Intensity_MinIntensity': array([0.]),\n 'Intensity_MaxIntensity': array([9.]),\n 'Intensity_IntegratedIntensityEdge': array([3253.]),\n 'Intensity_MeanIntensityEdge': array([4.32579787]),\n 'Intensity_StdIntensityEdge': array([2.87255435]),\n 'Intensity_MinIntensityEdge': array([0.]),\n 'Intensity_MaxIntensityEdge': array([9.]),\n 'Intensity_MassDisplacement': array([0.36112335]),\n 'Intensity_LowerQuartileIntensity': array([2.]),\n 'Intensity_MedianIntensity': array([4.]),\n 'Intensity_MADIntensity': array([2.]),\n 'Intensity_UpperQuartileIntensity': array([7.]),\n 'Location_CenterMassIntensity_X': array([98.79390369]),\n 'Location_CenterMassIntensity_Y': array([99.29653732]),\n 'Location_CenterMassIntensity_Z': array([0.]),\n 'Location_MaxIntensity_X': array([5.]),\n 'Location_MaxIntensity_Y': array([5.]),\n 'Location_MaxIntensity_Z': array([0.]),\n 'SpatialMoment_0_0': array([35721.]),\n 'SpatialMoment_0_1': array([3357774.]),\n ...\n 'SpatialMoment_2_3': array([7.04675332e+14]),\n 'CentralMoment_0_0': array([35721.]),\n 'CentralMoment_0_1': array([0.]),\n ...\n 'CentralMoment_2_3': array([0.]),\n 'NormalizedMoment_0_0': array([nan]),\n 'NormalizedMoment_0_1': array([nan]),\n 'NormalizedMoment_0_2': array([0.083331]),\n ...\n 'NormalizedMoment_3_3': array([0.]),\n 'HuMoment_0': array([0.166662]),\n ...\n 'InertiaTensor_0_0': array([2976.66666667]),\n 'InertiaTensor_0_1': array([-0.]),\n 'InertiaTensor_1_0': array([-0.]),\n 'InertiaTensor_1_1': array([2976.66666667]),\n 'InertiaTensorEigenvalues_0': array([2976.66666667]),\n 'InertiaTensorEigenvalues_1': array([2976.66666667]),\n 'Zernike_0_0': array([0.64333829]),\n ...\n 'Zernike_9_9': array([3.9043712e-19]),\n 'MinFeretDiameter': array([188.]),\n 'MaxFeretDiameter': array([265.87214973]),\n 'AngularSecondMoment_3_00_256': array([0.01237016]),\n 'Contrast_3_00_256': array([13.32731615]),\n 'Correlation_3_00_256': array([0.00052763]),\n 'Variance_3_00_256': array([6.66717588]),\n 'InverseDifferenceMoment_3_00_256': array([0.27445799]),\n 'SumAverage_3_00_256': array([9.97394114]),\n 'SumVariance_3_00_256': array([13.34138738]),\n 'SumEntropy_3_00_256': array([3.87614681]),\n 'Entropy_3_00_256': array([6.33842161]),\n 'DifferenceVariance_3_00_256': array([0.00372448]),\n 'DifferenceEntropy_3_00_256': array([2.98563863]),\n 'InfoMeas1_3_00_256': array([-0.00035948]),\n 'InfoMeas2_3_00_256': array([0.04771106]),\n 'AngularSecondMoment_3_01_256': array([0.01237633]),\n 'Contrast_3_01_256': array([13.34467007]),\n 'Correlation_3_01_256': array([0.00074894]),\n 'Variance_3_01_256': array([6.67733598]),\n 'InverseDifferenceMoment_3_01_256': array([0.27467073]),\n 'SumAverage_3_01_256': array([9.96441954]),\n 'SumVariance_3_01_256': array([13.36467386]),\n 'SumEntropy_3_01_256': array([3.87774892]),\n 'Entropy_3_01_256': array([6.33805833]),\n 'DifferenceVariance_3_01_256': array([0.00361481]),\n 'DifferenceEntropy_3_01_256': array([2.9890093]),\n 'InfoMeas1_3_01_256': array([-0.00041841]),\n 'InfoMeas2_3_01_256': array([0.05146776]),\n 'AngularSecondMoment_3_02_256': array([0.01237207]),\n 'Contrast_3_02_256': array([13.36134306]),\n 'Correlation_3_02_256': array([-0.0008637]),\n 'Variance_3_02_256': array([6.67490643]),\n 'InverseDifferenceMoment_3_02_256': array([0.27233695]),\n 'SumAverage_3_02_256': array([9.95842397]),\n 'SumVariance_3_02_256': array([13.33828266]),\n 'SumEntropy_3_02_256': array([3.87581509]),\n 'Entropy_3_02_256': array([6.33830874]),\n 'DifferenceVariance_3_02_256': array([0.00369281]),\n 'DifferenceEntropy_3_02_256': array([2.98609235]),\n 'InfoMeas1_3_02_256': array([-0.00037672]),\n 'InfoMeas2_3_02_256': array([0.04884048]),\n 'AngularSecondMoment_3_03_256': array([0.0123648]),\n 'Contrast_3_03_256': array([13.28753834]),\n 'Correlation_3_03_256': array([0.00346848]),\n 'Variance_3_03_256': array([6.66689316]),\n 'InverseDifferenceMoment_3_03_256': array([0.27439131]),\n 'SumAverage_3_03_256': array([9.97160996]),\n 'SumVariance_3_03_256': array([13.38003429]),\n 'SumEntropy_3_03_256': array([3.87790345]),\n 'Entropy_3_03_256': array([6.33872973]),\n 'DifferenceVariance_3_03_256': array([0.0036703]),\n 'DifferenceEntropy_3_03_256': array([2.98626744]),\n 'InfoMeas1_3_03_256': array([-0.0002397]),\n 'InfoMeas2_3_03_256': array([0.03896704]),\n 'Granularity_1': array([23.96446938]),\n ...\n 'Granularity_16': array([100.])}\n\"\"\"\n#+end_src\n\n*** Call specific measurements\nIf you need a specific measurement/feature you can just import it. Note that measurements come in sets, so you have to fetch the one that you specifically require from the resultant dictionary. Any available measurement can be found using code as follows:\n#+begin_src python\n import numpy as np\n\n from cp_measure.minimal.measureobjectsizeshape import get_sizeshape\n \n mask = np.zeros((50, 50))\n mask[5:-6, 5:-6] = 1\n get_sizeshape(mask, None) # pixels, the second argument, is not necessary for this measurement\n#+end_src\n\nThe other available functions are as follows:\n- measureobjectintensitydistribution.get_radial_zernikes,\n- measureobjectintensity.get_intensity,\n- measureobjectsizeshape.get_zernike,\n- measureobjectsizeshape.get_ferret,\n- measuregranularity.get_granularity,\n- measuretexture.get_texture,\n\nAnd for Type 2 functions:\n- measurecolocalization.get_correlation_pearson\n- measurecolocalization.get_correlation_manders_fold\n- measurecolocalization.get_correlation_rwc\n- measurecolocalization.get_correlation_costes\n- measurecolocalization.get_correlation_overlap\n \n\n* Pending measurements \nYou can follow progress [[https://docs.google.com/spreadsheets/d/1_7jQ8EjPwOr2MUnO5Tw56iu4Y0udAzCJEny-LQMgRGE/edit?usp=sharing][here]].\n\n*** Done\n- Type 1 and 2 measurements in sklearn style (process multiple masks per image)\n*** Pending\n- Type 3 and 4 measurements\n \n*** Additional notes\nThe Image-wide functions will not be implemented directly, as they were originally implemented as additional functions to the Object (mask) functions. We will adjust the existing functions assume that an image-wide measurement is the same as measuring an object with the same size as the intensity image.\n\n\n* Additional notes\n- This is not optimised for efficiency (yet). We aim to reproduce the 'vanilla' results of CellProfiler with minimal code changes. Optimisations will be implemented once we come up with a standard interface for functionally-focused CellProfiler components.\n- The functions exposed perform minimal checks. They will fail if provided with empty masks. Not all functions will fail if provided with masks only.\n",
"bugtrack_url": null,
"license": null,
"summary": "cp-measure implements CellProfiler Measurements. You can integrate them to your workflow or use them from CellProfiler 5 once they are integrated there.",
"version": "0.1.2",
"project_urls": null,
"split_keywords": [],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "9c81f1f865b918d46303b0ceda0d129a740c07644496789ec3fc383b115872b5",
"md5": "69d6ad3777632bf619435e82880c37e6",
"sha256": "90cc841427cfac51c30ef4b2b960d50fd4765400b5bdf75f2a9fd0cacfb675cf"
},
"downloads": -1,
"filename": "cp_measure-0.1.2-py3-none-any.whl",
"has_sig": false,
"md5_digest": "69d6ad3777632bf619435e82880c37e6",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": "<3.12,>=3.9",
"size": 36495,
"upload_time": "2024-09-27T19:40:26",
"upload_time_iso_8601": "2024-09-27T19:40:26.463379Z",
"url": "https://files.pythonhosted.org/packages/9c/81/f1f865b918d46303b0ceda0d129a740c07644496789ec3fc383b115872b5/cp_measure-0.1.2-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "3e34dd31db3434d55da4eaa30184d03ae21b3948a15a8b062f6ee7d19d014df9",
"md5": "b8ffce0b2dc04027a0eb0c691e4ffd1e",
"sha256": "c10b67da92b034a7ba2d8b010ae0f2590399e9a33cf2ad259bfe462b96272009"
},
"downloads": -1,
"filename": "cp_measure-0.1.2.tar.gz",
"has_sig": false,
"md5_digest": "b8ffce0b2dc04027a0eb0c691e4ffd1e",
"packagetype": "sdist",
"python_version": "source",
"requires_python": "<3.12,>=3.9",
"size": 80095,
"upload_time": "2024-09-27T19:40:27",
"upload_time_iso_8601": "2024-09-27T19:40:27.651842Z",
"url": "https://files.pythonhosted.org/packages/3e/34/dd31db3434d55da4eaa30184d03ae21b3948a15a8b062f6ee7d19d014df9/cp_measure-0.1.2.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-09-27 19:40:27",
"github": false,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"lcname": "cp_measure"
}
