SYSTEMS, PROCESSES, METHODS AND MACHINES FOR TRANSFORMING IMAGE DATA INTO SIZING AND VOLUME MEASUREMENTS FOR TISSUE

- Prodo Laboratories

Automated islet measurement systems (AIMS) in combination with tissue volume analysis (TVA) software effectively gauges volumetric and size-based data to generate heretofore unavailable information regarding, for example, populations of islet cells, stem cells and related desiderata.

Skip to: Description  ·  Claims  · Patent History  ·  Patent History
Description
CROSS-REFERENCE TO RELATED APPLICATION

This application claims the full Paris Convention benefit to and priority of U.S. Provisional Application Ser. No. 61/720,153, filed Oct. 30, 2012, the contents of which are incorporated by this reference as if fully set forth herein in their entirety.

BACKGROUND OF THE DISCLOSURE

Various tissue elements and other biological matter require precise data to support characterization, culturing and rendering the same effect for important usages and purposes. Optical imagery likewise provides a powerful tool.

Those skilled in the art readily understand that a first step in managing sensitive tissue can be found by capturing image data and arraying the same for various purposes.

According to the instant disclosures, scanning imagery technology, improved algorithms and alternate optical techniques may be used to inventively gauge, measure and validate aspects of tissue measurement, by volume.

OBJECTS AND SUMMARY OF THE DISCLOSURE

Briefly stated, automated islet measurement systems (AIMS) in combination with tissue volume analysis (TVA) software effectively gauges volumetric and size-based data to generate heretofore unavailable information regarding, for example, populations of islet cells, stem cells and related desiderata.

According to embodiments, there is disclosed an automatic tissue scanning camera (ATSC) system and methods to measure sizes and volumes of, for example, populations of islet cells within a culture flask.

According to embodiments, there is disclosed a system comprising tissue volume analysis (TVA) software, using thresholding to identify individual cells, for example islet cells, and to evaluate cross-sectional linear size of cells, for example islet cells, via an algorithmic simulacrum of manual counting under microscopy.

According to embodiments, there is disclosed a system, wherein said TVA software measures the optical transmission of each pixel within each cell, for example islet cell, using this to calculate thickness of tissue at that pixel based upon optical extinction; and, based on constant pixel area, calculates volume represented by respective pixels; whereby voxels, volume elements, are then summed to generate resultory total cell, for example islet cell, volume.

According to embodiments, there is disclosed a system wherein an automated image measurement algorithm, as defined in Appendix A, is employed.

According to embodiments, there is disclosed a process for generating sizing and volume measurement data by transforming optical imagery, comprising, in combination: providing an AIMS/imaging apparatus/scanner; having TVA algorithms and related software disposed therein; imaging at least a cell; generating predetermined fields for sizing and volume; conforming, registering and (or otherwise validating) the resultory data; and, using the same to create information which can be further compared to laser transaction microscopy and diffusing measurement data to confirm resulting data sets.

According to embodiments, there is disclosed a machine for transforming optical data into cellular sizing and volume measurement information, whereby scanned images are assembled further comprising at least one of islets and spheres of stem cells; and, in an independent unit validating the same, whereby confocal microscopy, likewise capable of diffusing measure, enables optical slicing of images of said islets and spheres of stem cells horizontally, using scattered light to document the subject area, which multiplied by thickness generates an independent validation of the volume generated by said machine.

BRIEF DESCRIPTION OF THE FIGURES

FIG. 1 is a schematic block/flow diagram showing operation of systems, methods, processes and products by the processes, according to the instant disclosures; and

FIG. 2 is a table distribution into depth X of permeant in a thin membrane.

DETAILED DESCRIPTION OF THE DISCLOSURE

Expressly incorporated herein by reference are U.S. Pat. Nos. 8,211,697; 8,213,081; 8,228,499; 8,293,223; 8,236,281; and 8,288,339.

The present inventor has discovered how to use new combinations of optical tools to characterize heretofore un-cognizable configurations of cells, constituents thereof and related desiderata having extensive implications.

Referring now to FIG. 1 and FIG. 2, algorithms embedded within software are accessed (including by remote, wireless and other means known or developed) to drive AIMS/imaging system/scanner 103, which generates data 111 regarding sizing and volume of subject cells, for example, islet cells (as described in Appendices A and B, for example).

Said data 111 is useful for characterizing the physical nature of subject cells, in this case, islet cells, whose spatial orientation and related behavior can now be understood as driven by data 111 and confirmed and registered with validated information—such as the more ellipsoid nature of the islet cells. Likewise, tissue transection microscopy 107 and diffusing measurement technologies, can be used to confirm these findings.

The Volume Measurement Algorithm for Individual Islets: This is a brief summary of the Automated Image Measurement Algorithm for measurement of tissue volumes. This system is based on the following considerations.

The volume of a right parallelepiped may be expressed as: V=A×h

The volume of an arbitrary object can be approximated by summing the discrete element volumes of the small parallelepipeds with uniform cross-section area, A, and height, h, that compose it: V=Σ(A×h)

The digital image of an object is a construction of rectangles whose individual areas in the image are from individual pixels each representing a specific area in object space (dependent on, for instance, “magnification”).

For absorption microscopy of biological tissues (i.e. non-phase contrast microscopy) without the use of stains, the 20 image corresponds to the optical “shadow” cross-section. The individual pixel measurements are a function of the thickness and the optical absorption of the tissue along the optical axis (this approximation is valid for thin sections) as described by the Lambert-Beers Law: Ip=10 exp (−a×tp) where Ip is the intensity measured by a pixel imaging tissue, 10 is the intensity of illumination and tp is the tissue thickness. Local optical transmission can be defined as TP=IP/10 so that TP=(−In TP)/α; where α=absorption coefficient.

Since the volume of an object can be represented by: V0=Σ(Ap×tp); where the area and thickness are for individual pixels representing the object; space, the volume may be calculated from the digital optical image of the object: V0=Σ(Ap×(−InTp)/α). V0=Σ(Ap×{−In(Ip/I0)}/α).

Uncertainty Analysis—Individual Volumes: The uncertainty of the area of an individual pixel in object space is dependent on the uncertainty of measurement of the microscope field size. This uncertainty is less than 2% (20 microns in a 3 mm field size). Σ(Ap×tp)

Thresholding: Objects are identified by binary thresholding so the summation uncertainty is dependent on thresholding: only pixels that are below a threshold intensity are identified as part of an object (particle detection) therefore exclusion of light areas in objects undercounts pixels in the summation. A binary closing function fills light areas identified as “holes” in objects so this defect will only affect thin portions surrounding the outer fringe of identified objects. The effect of thresholding has been evaluated functionally by measuring objects with thresholds set significantly above and below optimal (>6 units in approximately 130). The variation of measured volume between threshold extremes and optimal is less than 10%.

Calibration to micro-spheres: We have assessed particle detection and measurement reliability by comparing particle measurements with the image analysis system to NIST-traceable specifications for Coulter Counter calibrating spheres. The measurement of particle size (area) using image analysis compares very closely with the specifications for calibrating spheres.

TP=Ip/I0—measurement of light intensity: Light intensity measurement with CCDs: CCD pixel outputs are individually linear over many decades of illumination intensity. We correct for CCD pixel-to-pixel sensitivity and local variations of illumination and viewing response by flat-field image processing on all images prior to extraction of digital information. The flat field image, reflecting the effects of illumination, optics and the CCD, has a typical standard deviation of approximately 3% (6 parts in 200). After flat field correction this non-uniformity is reduced to virtually zero.

Test of linearity: Linearity has been confirmed (by use of overlapping calibrated density filters) to be on average within +/−2% over the range of intensity measurements of interest.

Algorithm for measuring transmission: The measurement of background intensity, I0, is based on averaging the measured intensity of light in a ring around each identified object. The ring is created by a series of 3-pixel wide object dilations followed by image subtraction. Errors in this measurement can arise if any field darkening, in particular by an adjacent particle, lies close enough to be within the ring. Our object image density is low (<1%) so this has not been observed to be a problem. This technique for measuring background minimizes the effect of light intensity variations locally per image and over the time of image acquisition.

Effect of focus: The accurate measurement of intensity by microscope is affected by focussing. Objects above or below the focal plane lose “contrast”, that is, approach the brightness of the background. We have measured the effect of placing the focus above and below object tissue masses, that is, above and below the nominally perceived best-focus position. The worst-case out-of-focus condition changes measured tissue mass by less than 10%. Since tissue typically settles rapidly to the bottom of the flask, this effect has the beneficial effect of helping to discriminate against out-of-focus non-tissue image artifacts such as floating fragments or outer flask wall defects; the threshold eliminates lightened artifacts.

Measurement of the absorption coefficient: The optical darkening that makes an object tissue visible in non-phase contrast microscopy is due to absorption of chromophores native to the object. For thin sections, where scattering is negligible, the average absorption is describable by the absorption coefficient, a. We have measured an average a for islet tissues and for NIT cells by measuring transmission of tissue at 650 nm captured in known spaces between microscope slides. The coefficient measured varies between 0.0033 um-1 and 0.0051 um-1 (see attached figure “measurement of absorption coefficient”). These numbers compare favorably to values in the literature for other tissues at this wavelength.

Specific value: For calculation of volume we have adopted the average value of 0.0033 um-1. For detection of relative changes in tissue volume over time, this consistent use of a single value is required.

Use of average absorption coefficient: The absorption coefficient varies for specific morphological features of a tissue. To minimize this effect we have chosen monochromatic, long-wavelength, visible light for illumination because generally the absorption coefficient for biological tissues decreases and converges for longer red wavelengths.

To evaluate the use of the average absorption coefficient algorithm, we squeezed NIT cells between microscope slides separated by rigid spheres of known size. The consequently flat, parallel-sided tissue disks were measured for volume by optical imaging. This was compared to calculation of volume as area times thickness. The regression (see figure “test algorithm”) was essentially a straight line, which confirms the use of an average coefficient for thickness calculation.

Test of volume measurement: The ability to accurately measure the volume of a single object was tested by imaging a plate-like porcine islet manipulated to edge-on and face-on positions. The difference in thicknesses presented by the two different views is a test of the validity of our optical transmission thickness measurement algorithm. Linear dimensions measured using the microscope image versus transmission were within 4% of each other and volumes measured were within 6% of each other (see attached figure “camps by view” and corresponding images).

Overall uncertainty of individual islet volume measurement: Based on the above the uncertainty of measurement of relative individual islet volumes is estimated to be less than 15%. The estimate of absolute volume measurement is dependent on further characterization of the absorption coefficient.

A flask containing a population of particles may be characterized by optical sampling. Since all the islets quickly settle to the bottom of our flasks, we can obtain an accurate sample of the population by imaging as much of the flask area as possible. The accuracy of this sampling is governed by sampling statistics where our population is non-homogeneous in particle size. The volume characteristics of the population can be calculated: VK=VKs×SK.

Where VK is the volume (or number) of islets in a histogram size bin K and VKs is the volume (or number) measured by optical sampling and SK is the sample multiplier for a given bin. The sample multiplier is: SK=1/P=AF/(Aim×m).

Where P is the fraction of the bottom of the flask which is sample by images, AF is the area of the flask bottom and Aim is the object-area of one image and m is the number of non-overlapping images taken. For our work P has been between 0.05 and 0.40 (currently about 10% sampling). The total volume in a flask is calculated: VT=ΣVK.

The uncertainty of sample analysis is based on sampling statistics. The uncertainty of a bin sample is inversely proportional to the square root of n, where n is the number of particles in a given bin. We have tested the sampling variations by measuring several times in succession a flask containing microspheres. The average fractional range (maximum to minimum divided by the average) for eight series was 6.3%.

We have compared the measurement results of image analysis particle counting to results with the Coulter Counter using identical samples of calibrated microspheres. The image analysis count was approximately 15% low for 68 micron diameter spheres and was identical for 200 micron spheres (see attached figure “compare Coulter to IA beads”). We have tested the sampling variations for tissue by measuring several times in succession a flask containing human islets. Preliminary results show similar consistency, although, as expected, the sampling uncertainty increases for larger size particles because the population contains ever fewer large volume particles.

Artisans readily understand the nature and extent of the instant system to be useful for stem cells and related desiderata.

Returning to islet cells, it has long been assumed that a generally spherical configuration was in play, impacting numerous empirical aspects of their storage, usage, and other key issues. Using the teachings of the present disclosure, the ellipsoid (or, American “football-like”) shape of islet cells was confirmed, enabling ongoing research heretofore uncontemplated as currently understood. Such teachings comprise progress in science and the useful arts, and likewise present subject matters for U.S. Letters Patents.

While the method and apparatus have been described in terms of what are presently considered to be the most practical and preferred embodiments, it is to be understood that the disclosure need not be limited to the disclosed embodiments. It is intended to cover various modifications and similar arrangements included within the spirit and scope of the claims, the scope of which should be accorded the broadest interpretation so as to encompass all such modifications and similar structures. The present disclosure includes any and all embodiments of the following claims.

It should also be understood that a variety of changes may be made without departing from the essence of the invention. Such changes are also implicitly included in the description. They still fall within the scope of this invention. It should be understood that this disclosure is intended to yield a patent covering numerous aspects of the invention both independently and as an overall system and in both method and apparatus modes.

Further, each of the various elements of the invention and claims may also be achieved in a variety of manners. This disclosure should be understood to encompass each such variation, be it a variation of an embodiment of any apparatus embodiment, a method or process embodiment, or even merely a variation of any element of these.

Particularly, it should be understood that as the disclosure relates to elements of the invention, the words for each element may be expressed by equivalent apparatus terms or method terms—even if only the function or result is the same.

Such equivalent, broader, or even more generic terms should be considered to be encompassed in the description of each element or action. Such terms can be substituted where desired to make explicit the implicitly broad coverage to which this invention is entitled.

It should be understood that all actions may be expressed as a means for taking that action or as an element which causes that action.

Similarly, each physical element disclosed should be understood to encompass a disclosure of the action which that physical element facilitates.

Any patents, publications, or other references mentioned in this application for patent are hereby incorporated by reference. In addition, as to each term used it should be understood that unless its utilization in this application is inconsistent with such interpretation, common dictionary definitions should be understood as incorporated for each term and all definitions, alternative terms, and synonyms such as contained in at least one of a standard technical dictionary recognized by artisans and the Merriam-Webster's Unabridged Dictionary, the latest edition of which is hereby incorporated by reference.

Finally, all references listed in the Information Disclosure Statement or other information statement filed with the application are hereby appended and hereby incorporated by reference; however, as to each of the above, to the extent that such information or statements incorporated by reference might be considered inconsistent with the patenting of this/these invention(s), such statements are expressly not to be considered as made by the applicant.

In this regard it should be understood that for practical reasons and so as to avoid adding potentially hundreds of claims, the applicant has presented claims with initial dependencies only.

Support should be understood to exist to the degree required under new matter laws—including but not limited to United States Patent Law 35 USC 132 or other such laws—to permit the addition of any of the various dependencies or other elements presented under one independent claim or concept as dependencies or elements under any other independent claim or concept.

To the extent that insubstantial substitutes are made, to the extent that the applicant did not in fact draft any claim so as to literally encompass any particular embodiment, and to the extent otherwise applicable, the applicant should not be understood to have in any way intended to or actually relinquished such coverage as the applicant simply may not have been able to anticipate all eventualities; one skilled in the art, should not be reasonably expected to have drafted a claim that would have literally encompassed such alternative embodiments.

Further, the use of the transitional phrase “comprising” is used to maintain the “open-end” claims herein, according to traditional claim interpretation. Thus, unless the context requires otherwise, it should be understood that the term “compromise” or variations such as “comprises” or “comprising”, are intended to imply the inclusion of a stated element or step or group of elements or steps but not the exclusion of any other element or step or group of elements or steps.

Such terms should be interpreted in their most expansive forms so as to afford the applicant the broadest coverage legally permissible.

Claims

1. An automatic tissue scanning camera (ATSC) system, which comprises, in combination:

a monochrome CCD camera;
a horizontally-mounted culture flask;
trans-illumination means, which means exceeds the deep field image size by a factor of at least two;
wherein placement of the flask between the illuminator and camera coaxially and opposed on crossed, ultra precision mechanical stages driven by processing means, generates rectangular arrays of continuous digital images, and;
whereby the density of any one cell, for example an islet cell, is derived as a function of scattering, absorption, and thickness.

2. A system according to claim 1, further comprising, tissue volume analysis (TVA) software, using thresholding to identify individual cells, for example islet cells, and to evaluate cross-sectional linear size of cells, for example islet cells, via an algorithmic simulacrum of manual counting under microscopy.

3. A system according to claim 2, wherein said TVA software measures the optical transmission of each pixel within each cell, for example islet cell, using this to calculate thickness of tissue at that pixel based upon optical extinction; and,

based on constant pixel area, calculates volume represented by respective pixels;
whereby voxels, volume elements, are then summed to generate resultory total cell, for example islet cell, volume.

4. A system according to claim 1, wherein an automated image measurement algorithm is employed.

5. A system according to claim 2, wherein an automated image measurement algorithm is employed.

6. A system according to claim 3, wherein an automated image measurement algorithm is employed.

7. A process for generating sizing and volume measurement data by transforming optical imagery, comprising, in combination:

providing an AIMS/imaging apparatus/scanner;
having TVA algorithms and related software disposed therein;
imaging at least a cell;
generating predetermined fields for sizing and volume;
conforming, registering and (or otherwise validating) the resultory data; and,
using the same to create information which can be further compared to laser transsection microscopy and diffusing measurement data to confirm resulting data sets.

8. A machine for transforming optical data into cellular sizing and volume measurement information, which, whereby a multiplicity of images are created using illumination spectra between at least about 500 and 700 nm to generate a plurality of sliced scanned images which are assembled based upon cross-sectional linear size used to calculate cell volume for a population, further comprising at least one of islets and spheres of stem cells; in combination with optionally, at least an independent unit for validating the same, whereby confocal microscopy, likewise capable of diffusivity measurement, enables optical slicing of images of said islets and spheres of stem cells horizontally, using scattered light to document the subject area, which when multiplied by thickness, generates an independent validation of the volume generated by said machine.

9. The machine of claim 8, the islet cells being at least one population selected from the group of large animals, humans and swine.

Patent History
Publication number: 20140119633
Type: Application
Filed: Oct 11, 2013
Publication Date: May 1, 2014
Applicant: Prodo Laboratories (Irvine, CA)
Inventor: Roy McCord (Irvine, CA)
Application Number: 14/052,003
Classifications
Current U.S. Class: Cell Analysis, Classification, Or Counting (382/133); Including Optical Measuring Or Testing Means (435/288.7); By Particle Light Scattering (356/336)
International Classification: G01N 15/02 (20060101); G06K 9/00 (20060101); G01N 21/51 (20060101);