text:
|Publication number
|Publication type
|Application number
|Publication date
|Jun 12, 2007
|Filing date
|May 15, 2003
|Priority date
|May 17, 2002
|Fee status
|Also published as
|Publication number
|David L. Raunig
|Original Assignee
|Pfizer Inc.
|Export Citation
|BiBTeX, EndNote, RefMan
|Patent Citations (7), Non-Patent Citations (11), Referenced by (15), Classifications (25), Legal Events (4)
|External Links: USPTO, USPTO Assignment, Espacenet
This application claims priority to provisional application No. 60/381,064 filed on May 17, 2002 by David L. Raunig, entitled "Apparatus and Method for Statistical Image Analysis", the entire contents of which are incorporated herein by reference, and for which benefit is claimed under 35 U.S.C. �119(e).
This invention relates generally to an apparatus, a system, and method for the analysis of images, and more particularly to a method for the analysis of ultrasonic images. This invention further relates even more particularly to an apparatus and method for the analysis of ultrasonic images using multivariate detection techniques, which, in one exemplary embodiment, may be used to differentiate image pixels representing necrotic material from image pixels representing viable tissue within a specified region of interest, and which, in a further aspect of the invention, may be used in a method for treating cancer.
Changes in tissue types are manifested in ultrasound propagation as changes in acoustic properties, which in turn lead to the bright echoes from the interfaces. In the case of tumors, the interior tissue is generally a homogeneous mixture of viable cancer cells, capillaries, connective tissue and necrotic cells. The amount of necrosis can provide some valuable clues to the viability of the tumor, itself.
Conventionally, since capillaries are the sites of the vital interchange of blood nutrients and cellular metabolites in healthy tissue, a measurement of vascularity may be used as an indicator of tissue viability, often calling for excision of a treated section of tissue for histopathology. The tissue is stained for vascular endothelia, usually the stain CD31, and the microvascular density (MVD) is measured under a microscope. The MVD is a standard measurement of the efficacy of angiogenic and anti-angiogenic drugs. Additionally, optical microscopy can also be used to image and measure the amount of dead, or necrotic tissue within a tumor.
Typical histopathology analysis of a tumor involves removal of the tumor from the animal, fixing the tumor with a fixative agent, usually Neutral Buffered Formalin (NBF) or ethanol (ETOH), and examining the specimen one thin slice at a time. This conventional approach necessarily requires more specimens (i.e. animals) than a longitudinal experiment over time, and is very labor intensive.
Necrosis in tissue occurs as a result of a fatal injury to a cell. Cells, upon death, elicit an inflammatory response and quickly deteriorate. Cell structure disintegrates, cell walls break down (lyse) and the fragmented cell contents, including nuclear material, spill into the interstitial spaces between the cells. Necrotic tissue even appears different from viable cells under the microscope, and it is recognizable by the complete absence of structure. Necrosis is histologically identified by the structure of the tissue using stains differential to nucleic proteins. However, histopathology is necessarily post-mortem, and a longitudinal analysis of the tumor over time is not possible.
Additionally, cancer treatments that target the cells increase the regions of necrosis around the intratumoral blood vessels. The presence of necrosis can therefore be used to analyze such characteristics as tumor morphology, tumor angiogenesis and efficacy of anti-cancer therapy; however, the study of these metrics has traditionally been done using optical microscopy and histopathology biomarkers. In laboratory research, where the number of animals used for research can be quite large, histopathology measurements are time consuming, expensive and require terminating the study of the subject. Additionally, tissue fixing can result in damage to some aspect of the tumor and loss of information.
As indicated above, past study of tumors and development of anti-cancer therapies have primarily relied on terminating the study at specified time points and performing histological tests on the excised tumor. This process is time consuming, requires the use of a large number of subjects, and necessarily precludes the use of longitudinal measurements on the same tumor. In addition to the statistical problems posed by the lack of repeated measurements, the histological analysis of tumors is done using microscopically thin slices of the tumor, providing information on a very limited region, in a non-realtime manner.
Ultrasonic imaging techniques are widely used in imaging portions of the body, including analysis of tumors. Changes in tissue types are manifested in ultrasound propagation as changes in acoustic properties, which in turn lead to the bright echoes from the interfaces.
Clinical ultrasound instruments can be classified into two broad categories: Continuous Wave and Pulsed Wave. Continuous Wave ultrasound is used primarily for its Doppler capabilities as an audio aid to detect flowing blood. Pulsed ultrasound is used for imaging, and nearly all modern pulsed wave ultrasound devices employ color representations of blood flow (CD) to associate movement to a pixel.
Ultrasound imaging equipment has made significant improvements, and now has capabilities that were not possible only a few years ago. With an increase in frequency and resolution, it is now possible to use ultrasound in a number of clinical settings that were previously reserved only for X-ray imaging. Increased resolution and decreased noise can now provide clear images able to detect details that would have, in many situations, required extensive interpretation in the past. In the late 1980's and early 1990's, ultrasonic tissue characterization had mixed results due to the characteristic coherent interference phenomenon known as "speckle". With the improvement in devices and the increase in ultrasound frequency, also came an increase in the bandwidth with concomitant speckle reduction. The engineering advances in clinical ultrasound equipment can now allow tissue characterization, at least in research environments that were previously reserved for special instrumentation.
The ultrasound image is characterized by pixels of differing intensities that may be used to differentiate different regions of interest. Additionally, some tissue types can be subjectively identified by their texture. Intensities of the necrotic core of a tumor are brighter than the viable border tissue, for example. Observing necroses obtained from histopathologic sampling under the microscope and correlating the observed features with the corresponding ultrasound images, shows that there is a coarse texture to the necrotic regions, as illustrated in
In the case of tumors, the interior tissue is generally a homogeneous mixture of viable cancer cells, capillaries, connective tissue and necrotic cells. As mentioned above, the amount of necrosis can provide some valuable clues to the viability of the tumor, itself. Tumor cells lack the normal, programmed cell-death mechanisms that keep normal cells in check. However, tumor cells, like all living cells, can die if exposed to a fatal environment. For a new tumor without medical intervention, the cells farthest away from the peritumoral regions and outside the perfusion range of newly formed microvasculature die from lack of oxygen and nutrients.
The ability of ultrasound to discriminate between viable and necrotic tumor regions has previously been accomplished using the ultrasound backscatter microscope (UBM). The necrotic core of a tumor spheroid is known to be detectable using very high ultrasound frequencies. There are morphological differences between the two regions, as determined by histopathology. In a very high frequency (e.g., 100 MHz) ultrasound image, the necrotic core appears as a bright circular region in the middle of the tumor spheroid, because of the increased density of nuclei. However, such high frequency ultrasound devices are relatively expensive, and are not readily available in clinical settings, in comparison to lower frequency devices, e.g., devices operating in the 4–12 MHz region, for example.
There are three primary characteristics of tissue that can increase backscattered energy, and result in brighter pixels in the resulting image:
However, previous efforts in analyzing image pixel data to place the pixels into one of two groups, e.g., necrotic tissue and viable tissue, have used conventional univariate approaches, as shown in
What is needed then is an apparatus and method that allows more powerful statistical test of differences over time, and which includes the repeated longitudinal measurement of a tumor as it progresses in a single subject.
What is further needed is an apparatus and method for ultrasound imaging, using multiple image characteristics, to provide a means of performing non-histopathologic longitudinal segmentation of an image, including identification and differentiation of tumor necrosis and viable tissue.
What is still further needed is an apparatus and method for multivariate analysis of necroses using tissue differentiation based on pixel gray-scale values, augmented by differences in texture to increase the power of the statistical analysis and lower the incidence of misclassification of the different tissue types.
The present invention solves the aforementioned problems of conventional tumor analysis using histopathology, and the problems of univariate statistical approaches which result in undesirably high rates of misclassification.
The apparatus and method of the invention has, as one object, the provision of a measurement modality for repeated measurements of a tumor that includes differentiation between viable and necrotic tumor tissue, as well as measurements of peritumoral blood flow.
Another object of this invention is to provide a multivariate analysis method to characterize tumor tissue into necrotic and viable regions and provide a surrogate measurement for functional intratumoral microvessel blood flow.
These objectives are accomplished by the use of powerful multivariate statistical analysis methods that take advantage of the joint probability distribution function to increase both specificity and sensitivity of the image segmentation. In addition to the immediate advantages to the specific application to cancer research, this method may be adapted to the segmentation of other tissue regions of interest. The experimental results and the analysis of images obtained by the apparatus and method of the invention show that the multivariate analysis of ultrasound images more accurately characterize tumors by properly segmenting viable from necrotic tissue.
In addition to the immediate advantages of one embodiment of the invention directed to cancer treatment, this method is easily adapted to the segmentation of other tissue regions of interest or, in general, to segmentation of a wide variety of image types into various regions of interest, based upon multivariate statistical analysis of image features corresponding to actual sample features, for example, a surface feature.
Biomedical imaging modalities can provide a quick, non-invasive means of obtaining information on the interior of the tumor. While some imaging modalities are cumbersome and expensive for small animal research, modern ultrasound provides an inexpensive and relatively easy to use clinical imaging modality that can provide detail not previously available, and the enhanced detail allows the analyst to measure the image texture, and use this information to classify regions.
For example, ultrasound is non-invasive, non-ionizing, inexpensive and easily adaptable for use in clinical settings for imaging humans and other animals. Additional Doppler information available within the ultrasound image provides an added dimension for tumor image analysis that includes a surrogate measure for angiogenesis. Further, ultrasound allows for a detailed examination of the tumor over time without excision from the body, saving time and resources, both animal and human.
Ultrasound techniques allow measurement of the actual cross section of the tumor for a particular image slice, and provide the opportunity to use other slices to better estimate the volume of the tumor. Additionally, ultrasound imaging allows the investigator to study tumors in situ, and the tumor growth can be analyzed at intermediate intervals. However, whole-tumor volume calculation alone does not discriminate between viable and necrotic tissue. Measurements taken on the same tumor over time takes advantage of the correlation of those measurements within the same tumor specimen.
One embodiment of the invention includes a method for the analysis of image data represented by multiple pixels, the method characterized by determining a texture value for at least a portion of each of the multiple pixels; initially assigning each of the at least a portion of multiple pixels to one of two or more groups based upon respective determined texture values; calculating an initial statistical distance between at least one pixel and each initial centroid of the two or more groups; comparing each of the calculated initial statistical distances of the at least one pixel with each other; and depending upon the comparison result, reassigning the at least one pixel into a different one of the two or more groups.
In another aspect of this embodiment, a method for the analysis of image data represented by a plurality of pixels is disclosed, the method comprising measuring an intensity level of each of the plurality of pixels; initially assigning each of the plurality of pixels to one of two or more groups; determining an initial centroid for each of the two or more groups; calculating an initial statistical distance between each pixel of the plurality of pixels and each initial centroid for each of the two or more groups; for each pixel of the plurality of pixels, comparing the initial statistical distance for the assigned group with each of the other initial statistical distances associated with respective remaining ones of the two or more groups; reassigning at least one of the plurality of pixels into a different one of the two or more groups based upon the comparison result, wherein, for each pixel of said plurality of pixels, said initial statistical distance is calculated using at least two parameters selected from the group consisting of a pixel intensity, a pixel contrast, and a pixel entropy.
Another embodiment of the invention includes an image analysis system suitable for performing multivariate differentiation between necrotic and viable tissue in a tumor, wherein the system is characterized by an imager outputting image data in a form of pixels; a memory unit storing said image data; a processor operatively responsive to said image data, said processor calculating a vector representation of the image data which includes an intensity, a contrast, and an entropy value for each of the pixels; wherein said processor assigns said each of the pixels into either a necrotic tissue class or a viable tissue classed based upon a texture determined from said vector representation.
In another aspect of this embodiment, an image analysis system suitable for performing multivariate differentiation between necrotic and viable tissue in a tumor is disclosed, wherein the system comprises an imager outputting image data in a form of pixels; a memory unit storing said image data; a processor operatively responsive to said image data, said processor calculating a vector representation of the image data which includes an intensity, a contrast, and an entropy value for each of the pixels; wherein said processor assigns said each of the pixels into either a necrotic tissue class or a viable tissue classed based upon a texture determined from said vector representation, wherein said processor calculates a necrotic Mahalanobis distance between each pixel of the image data and a centroid of the necrotic tissue class and a viable Mahalanobis distance between each pixel of the image data and a centroid of the viable tissue class, wherein, for each pixel of the image data, said processor compares the necrotic and viable Mahalanobis distances and reassigns at least one pixel to a different class based upon the comparison.
In yet another aspect of this embodiment, an ultrasonic imaging apparatus is disclosed, which is characterized by an ultrasonic transceiver which transmits an ultrasonic signal, receives an ultrasonic signal reflected from a sample, and provides an output signal which includes pixel image data representative of the sample; a processor coupled to said ultrasonic transceiver, said processor processing at least a portion of the pixel image data to determine a cooccurrence matrix corresponding to a plurality of image pixels in the pixel image data, wherein said processor classifies each of the plurality of image pixels into one or more classes based upon the cooccurrence matrix.
Yet another embodiment of the invention is directed to a computer readable medium containing computer-executable code for implementing image analysis which is characterized by performing the functions of accepting image data from a plurality of pixels; classifying each of the plurality of pixels in one of a plurality of classes; statistically processing the image data and calculating at least a contrast and an entropy value for each of the plurality of pixels; and reclassifying at least a portion of the plurality of pixels into a different one of the plurality of classes based upon the statistically processed image data
In another aspect of this embodiment, a computer readable medium containing computer-executable code for implementing image analysis is disclosed which includes accepting image data from a plurality of pixels; determining an intensity of each of the plurality of pixels; classifying each of the plurality of pixels in one of a plurality of classes; statistically processing the image data and calculating a contrast and an entropy value for each of the plurality of pixels; and reclassifying at least a portion of the plurality of pixels into a different one of the plurality of classes based upon the statistically processed image data, wherein the computer-executable code for implementing image analysis further includes accepting tumor image data and classifying said at least a portion of the plurality of pixels into one of a necrotic tissue class and a viable tissue class based upon the associated contrast and entropy values.
The present application presents a novel, non-obvious multivariate method of tumor tissue differentiation, using present-day, readily available clinical ultrasound technology. The method and apparatus of the present application makes use of the fact that the multivariate use of ultrasound information, including signal intensity and texture measurements, increases the ability to discriminate between regions of necrotic and viable tumor tissue.
The above capability allows for a more detailed in vivo analysis of tumor growth beyond simple dimensional measurements and may be useful in the analysis of tumor angiogenesis biomarkers from Power Doppler. Consequently, longitudinal studies can be done that take advantage of the within-subjects correlation of repeated measurements, reducing the error variance, expanding the statistical inference to the population of subjects and increasing the power of the analysis.
In one aspect of the invention, the k-means clustering algorithm is modified for use with multiple classification variables. K-means clustering is a nonhierarchical method of pixel clustering, and was chosen over other methods for its relative simplicity and ability to use prior information about the clustering variables. Also, k-means clustering as used here, does not assume a minimum contiguous cluster of pixels, since the necroses may vary significantly in size. Known hierarchical methods, such as Bayesian techniques or region growing/splitting, do not provide for intensity grouping reassignment in the case of error. In one aspect of the invention, multivariate description of each pixel is assigned to each pixel, and this multivariate descriptor is used to cluster the pixels into one of two classes—viable or necrotic tissue.
In a further aspect of the invention, the disclosed and claimed invention provides a novel and non-obvious apparatus and method for multivariate segmentation of ultrasound images, which can be used to measure the necrotic core within a hard-cell carcinoma.
The features and advantages of the invention will be more readily understood upon consideration of the following detailed description of the invention, taken in conjunction with the accompanying drawings in which:
A first aspect of the invention will be explained with reference to
The geometry of a single voxel 310 of scatterers, corresponding to a single pixel within the image, is depicted in
Additionally, there are many regions within the necroses where the dense nuclei are interspersed with regions of cytoplasm. Therefore, necrotic regions are brighter than viable tumor tissue by the higher scatterer density and the increase in reflectivity.
Another characteristic of necrotic tissue that contributes to the increased necrotic pixel brightness is the tendency of the pyknotic nuclei to clump into non-homogeneous distributions within a single pixel. The clumping leads to a shift in the distribution of the pixel amplitudes, and an increase in the median pixel amplitude. The distribution of the received amplitudes for the non-specular regions within a homogenous volume of tissue is typically Rayleigh.
Since the necrotic tissue visibly appeared to be brighter and easily recognizable, and the histopathology images show the necrotic regions to have coarse irregularities associated with the texture, the distribution of necroses has a different type of distribution than viable tissue, due to the texture differences.
These types of differences lend themselves to being measured by ultrasound techniques. Clinical ultrasound instruments can be classified into two broad categories: Continuous Wave (CW)and Pulsed Wave. Continuous wave ultrasound is used primarily for its Doppler capabilities as an audio aid to detect flowing blood. Pulsed ultrasound is used for imaging, and nearly all modern pulsed wave ultrasound devices employ color representations of blood flow (CD) to associate movement to a pixel. Therefore, the ability to determine, in vivo, the extent of the necrotic core ties in well with contemporary studies in tumor angiogenesis.
As generically illustrated in
Equivalently, the power is proportional to the density of scatterers. The theoretical amplitude of the return signal assumes that all of the scatterers within the region of insonification have the same velocity and backscattering cross-section. Deviations from this assumption decrease the amplitude of the maximum Doppler frequency, and distribute the energy to other frequencies.
The reflected signal is preferably received by a receiving transducer 435 that is located near enough to transmitter 415 that the two can be approximated as collocated. The two transducer elements 420, 435 may be tilted slightly toward one another, and placed in a resin that acts as a focal lens. The beams of both transducers preferably intersect at a pre-defined range, and within this intersected range is the optimum axial range for detection. The range may be designed at the factory to be specific to a particular clinical use of the transducer, for example. Range calculations with CW are not possible since the transmitted time cannot be measured. Since the CW transmission is made at a single frequency, any shift in the received frequency can be determined very accurately as a constant Doppler shift within the sample period.
As shown in
Though a CW ultrasound system can accurately determine the Doppler shift due to motion in the line-of-insonification (LOI), the lack of range information does not allow the Doppler information to be tied to a particular area of interest. Typically, CW ultrasound is used in the audio mode, and placed to insonify tissue where the signal is generally known. For example, CW can be effectively used to detect artifacts of blood flow such as fetal heartbeat or carotid pulse.
Pulse Doppler and later-developed Color and Power Doppler, allow using information from many pulses to extract both range and frequency information. In this way, motion may be indicated on an image display (not shown) with a color map that represents different frequencies or velocities for specific tissue of interest. Such devices are also generalized in
Processor 530 preferably calculates various statistics related to image data 515, including, for example, intensity, entropy, and/or contrast for all or a portion of the pixels representing the image of sample 540. In a further aspect, processor 530 may calculate a vector representation of image data which may include the above statistical values and from which a texture of the surface of the sample may be represented. For example, we have verified that a relatively rough texture is correlated with pixels which have a relatively high contrast value, and a fine texture is correlated with pixels having a relatively high entropy value.
Further, in tumor analysis, for example, processor 530 preferably classifies pixels (e.g., either all image pixels or all image pixels in a region of interest) into a class of tissue, for example, either necrotic or viable tissue. Still further in this regard, processor 530 is preferably programmed to calculate the Mahalanobis distance (MD) (defined infra) between each pixel of all or a portion of the pixels and a centroid of either or both the viable tissue and necrotic tissue. In a further aspect, processor 530 may be programmed to determine the cooccurrence matrix corresponding to the pixel image data which may be used to classify the pixels into one or more classes, e.g., necrotic or viable tissue in the case of tumor analysis.
Reclassification of a pixel from the necrotic tissue class to the viable tissue class could occur when, for example, the MD for the pixel with respect to the centroid of the viable tissue group becomes less than the MD for the pixel with respect to the centroid of the necrotic tissue group. Similarly, reclassification of a pixel from the viable tissue class to the necrotic tissue class could occur when the MD for the pixel with respect to the centroid of the necrotic tissue group becomes less than the MD for the pixel with respect to the centroid of the viable tissue group. Although the centroids of these groups may shift as pixels are added or removed from a group, if a pixel is properly assigned to a group, the MD for the pixel with respect to the centroid of its assigned group will remain less than the MD to the centroid of the non-assigned group.
In terms of actual hardware used in one embodiment, a Parks 1051 CW Doppler ultrasound instrument is used to transmit an 8.1 MHz ultrasound wave into various types of biological tissue including necrotic tissue in a tumor. The frequency-shifted signals were collected, and analyzed digitally. In another embodiment, an ATL HDI�5000 ultrasound system was used for imaging done on subcutaneous tumors. The HDI�5000 uses a mean speed in tissue of 1540 cm/s with an error in non-specific tissue, assumed by ATL to be random, with a range of about 2–5%. In an exemplary embodiment, the PRF selected for the subcutaneous tumors was 1500 Hz, allowing resolution of range to approximately 70 μm. The lower limit of frequency detection for the ATL HDI�5000 was approximately 50 Hz, the low setting for the wall filter.
After reception and demodulation of the received ultrasound signal, the image data, either from a CW or pulsed Doppler process, may be stored in a memory, prior to or concurrent with statistical processing in a microprocessor, for example.
The interior of a tumor is not a homogenous mass of cancerous cells but a mixture of different regions that have different functions in the growth of the tumor and the supporting biology. The regions of interest in an ultrasound image of a subcutaneous tumor may include, for example, ultrasound gel, viable tumor cells, necrosis, muscle, fat and skin. The gel, muscle, skin and fat define the edge of the tumor, and are different enough in acoustic impedance from the tumor interior to provide a fairly distinct border in the image. The primary acoustic tissues within the tumor are viable tumor cells, areas of necrosis, blood vessels and fat vacuoles. The viable tumor cells have nearly identical acoustic characteristics to normal tissue, and therefore appear much the same as the surrounding healthy tissue.
Necrotic tumor cells are a common characteristic of nearly every intact tumor that has grown beyond a minimum size. As the cells divide and the tumor grows, the interior cells grow further away from the capillaries that provide oxygen and nutrients from the blood. Tumor cells have the same or greater nutritional needs as normal cells and those that are too far from the servicing capillary will die in the absence of oxygen. Viable cells that have sustained massive injury or, for lack of oxygen, use up their anaerobic energy stores (adenosine triphosphate—ATP), undergo cellular swelling, plasma membrane disruption and leakage of internal contents, including some nuclear cytoplasm, into the extra-cellular space.
Unlike programmed cell death, or apoptosis, the necrotic cell nucleus does not undergo any morphological change and remains relatively intact, except for some of its fluid content. The nucleus becomes denser, or pyknotic and is no longer constrained in its proximity to other nuclei by the cell. The cytoplasm diffuses through the extra-cellular matrix resulting in some areas of a dense mixture of cytoplasm and pyknotic nuclei and other areas where there are few nuclei. The dense regions appear brighter than viable regions and can be detected at clinical ultrasound frequencies. In some instances, the nuclei are dense enough so that pools of cytoplasm form. If the regions of cytoplasm are not filled with cellular debris, then the corresponding image region will be more like viable tissue and should appear as amorphous hypoechoic areas within the hyperechoic necrotic region.
While the present application is concerned in one aspect with the imaging capabilities of an ultrasonic instrument to discriminate between different types of tumor regions, functional intratumoral blood flow is directly related to the degree of viable tissue. Therefore, the ability to determine, in vivo, the extent of the necrotic core ties in well with contemporary techniques in tumor angiogenesis. Further, since perfusion and viability go hand-in-hand, the role of ultrasound imaging and Doppler ultrasound blood flow measurement complement each other in the clinical research environment.
In an example of an ultrasound image of a subcutaneous tumor shown in
For model validation, the a priori knowledge of the distribution of the intensities of each region listed above within a learning set of images can be used to classify the pixels into a particular region. Because there is no guarantee that each image will be statistically identical, the images are preferably stretched and histogram equalized to a uniform 8-bit gray scale map to allow for differences in imaging intensities due to operator, animal, etc. Histogram equalization maps the gray-scale intensity distribution of the images to the desired distribution typically using a maximum likelihood mapping function.
The region centroids and dispersions (standard deviations) provide an initial starting point for clustering. Differences between ultrasound images may distort centroids and dispersion characteristics, but will not change the rank ordering of the groups. Therefore, once initially assigned to a particular class, the pixels can be reassigned depending on their "distance" to each class. In a preferred embodiment, the "distance" is measured as the likelihood function of the multivariate score assigned to the pixel. The algorithm is similar to that shown in
Given normally distributed image clusters, the original pixel groups can be randomly assigned in a n�ive clustering, where no prior information is available. In this case, the first reassignment is usually the largest, and subsequent reassignments depend on the distribution of the data. However, the random initial assignment of pixels to groups can cause some problems in the final clustering. If two or more of the original seed points happen to be too close, then the resulting clusters may actually represent the same type of data and the reassignment of pixels would result in a large error. This problem can be reduced if the seeds are determined prior to the original assignment based on a training set of images. The pixel group assignments then update the centroids, and the minimum squared distances, i.e., the Mahalanobis distances, are then used to reassign the pixels and iterate until no reassignments occur.
Traditional k-means clustering usually starts out with a predefined number of clusters. However, if the number of clusters is fixed, outliers can be transferred into a cluster and seriously bias the true cluster centroid, thus distorting subsequent reassignments. Also, if there were no data points that populate a cluster, then the algorithm for updating cluster centroids would require special rules to account for this. Fortunately, neither of these conditions was observed in extensive investigations of mouse tumor tissue. The use of a Gaussian smoothing mask helped reduce the effect of outliers. The size of the image of the tumor was on the order of 50,000 pixels, and the large number of elements in each grouping made the probability of these two problems negligible. Limiting the clusters to four for the original image segmentation to identify the tumor boundary, and to two for segmentation of the tumor into viable and necrotic tissue may preferably minimize the computational load associated with group assignment and reassignment. Methods exist for using similarity measures to optimally determine the number of clusters, but the regions were defined well enough and limited to a small enough area that this was not necessary.
A relatively easy determination of the centroid is the mean of the group. The mean works well if the data is symmetrically distributed, and is ideal for Gaussian distributions when the standard deviation can be reliably estimated. In the case of ultrasound imaging, pixel intensities within a homogenous media are Rayleigh distributed, and the assumption of Gaussian behavior would result in some misclassification. The actual intensity probability distribution function can be determined using the Gibbs sampler, which is a method of accurately determining the probability distribution without requiring the actual functional form of the distribution.
Though the Gibbs sampling algorithms are relatively fast, the computational load can still be expensive for large images. An alternative to the Gibbs sampler for ultrasound images is low frequency filtering of the image to reduce the effects of speckle. This low-pass filtering approach takes advantage of the Central Limit Theorem (CLT) to transform the image into a suitable form which allows for the use of more mathematically tractable Gaussian statistics. The CLT states that, given
Y=Z 1 +Z 2 +Z 3 + . . . +Z n,
where Z1, Z1, . . . Zn are independent and identically distributed, then the distribution of Y is asymptotically Gaussian. This tendency is true for nearly all distributions (Cauchy is one exception), and the convergence is very rapid for symmetric or nearly symmetric data distributions. If a Gaussian mask with a standard deviation that is not too small is used, then the addition of the weighted pixel intensities results in a filtered intensity distribution that is less right skewed, or more normally distributed.
Therefore, to reduce the error of centroid determination based on the mean of the group, and to help reduce the non-speckle type image noise, the pixel data may be filtered with a 3�3 Gaussian filter with standard deviation of 0.7 pixels, for example. The weighted addition of the pixels result in a smoothed data point that is symmetrically distributed but which does not unacceptably blur the edges. With this approach, the resulting smoothed pixels are Gaussian-like distributed, without sacrificing the differentiation of neighboring pixels.
Once an assignment of the pixels to a group is made, the probability, or likelihood, of a pixel belonging to the jth group, given that it was assigned to that group in the previous iteration, is
Λj(x l ,k)=p(x l |M j(k),x j(k−1),S j(k−1))
where xl is the ith pixel intensity in the kth iteration, Mj is the centroid of the jth group and Sj(k−1) is the variance of the jth group. For pre-determined group transition probabilities, the probability of a pixel being in group j at iteration k if it was in group i at iteration (k−1) is then defined as
The term 1/c is a normalization constant that is the same for all groups, and can therefore be disregarded for comparisons of group probabilities. In a simple, conventional form of k-means clustering, the groups are assumed to be equally dispersed, and the likelihood function is the simple difference between the intensity value of the ith pixel and the centroid of the jth cluster.
D ij =|X i −C j|.
In this case, the normalization constant would be the maximum possible difference between values, or 255 for an 8-bit gray scale. The images analyzed in this research were histogram-equalized to a uniform gray-scale map to stabilize the intensity differences between images. Therefore, each image had intensities that spanned the gray-scale range (0–255). The reassignment of pixels would then be dependent only on the distance, Dij, and so the probability transition values will be 1/r for all pij.
One problem with the conventional univariate distance function Dij as above is that each segment is assumed to have equal dispersion characteristics for all intensities. A pixel with an intensity that was numerically closer to one segment with a narrow distribution of intensities may actually belong to the adjoining segment with a wider dispersion. For a Gaussian distribution, the probability of being in one segment is a function of the standardized distance for each cluster, which measures the distance from the cluster centroid by the intensity standard deviation as shown:
where σj is the standard deviation of the jth cluster. This type of metric allows for different cluster dispersion characteristics. However, the variance of each of the clusters is preferably not updated at each iteration. Movement into a cluster is likely to happen from values near the edges of the cluster of the previous iteration, increasing that cluster variance, and decreasing the variance of the cluster that was exited. As the variance increases in the first cluster, nearby pixel values will be further reassigned to that cluster, and the variance will continue to increase until all pixels are undesirably assigned to a single cluster.
If prior information from a training set of images is sufficient to characterize each cluster, then the use of the prior information for all iterations is sufficient to establish the probability of cluster membership, with the assumption that the prior variance/covariance estimate is the best estimate.
If, however, segment information is multivariate and the variables are reasonably, but not necessarily, uncorrelated, a situation particularly well-suited to use of the apparatus and method of the invention, then k-means clustering may preferably be generalized using a standardized Euclidean distance that takes advantage of both the variance and covariance of the variables. This standardized distance, also known as the Mahalanobis distance (MD), is the generalized Euclidean form of the standardized distance,
D 2 ij=(X i −C j)Σ−1(X i −C j)′,
where Xi is the ith pixel measurements, Σ is the multivariate variance/covariance matrix and Cj is the centroid of the jth group. The multivariate MD, then, is the multivariate analog of the square of the standardized statistical distance previously discussed.
For two segmentation groups, if the groups are multivariate normal, then the probability of being in the jth group is
for p variables. The classification criteria can then be chosen by considering the ratios of probabilities. Using the above equation for j=1,2, the ratio of probabilities is
The criteria of choosing group 1 over group 2 can then be designed into the decision process based on the cost of misclassification. If the cost of misclassification is equal between groups, then the ratio above is compared to 1. Alternatively, if the log of the likelihood ratio is used, then the criteria becomes
If the left side of the above equation is positive, then group 1 is chosen, and if negative, then group 2 is chosen.
If the cost of misclassification has equal weight for both groups, then the probability of misclassification is known to be asymptotically normal. The argument expands to three dimensions for multivariate normality. However, pixel values will nearly always fall between two groups, and the decision for group assignment will be reduced to two primary groups, with the probability of assignment to other groups being much lower than the primary group, and thus able to be ignored.
In the case of variables that are not normally distributed, there are a number of transforms that result in a transformed distribution that closely approximates the normal distribution. These distributions may be determined by the Box-Cox transform, which determines the linear relationship between the log of the binned standard deviation and the log of the binned means. These transforms work for most of the class of exponential functions.
The MD values for the pixels were used in the development of the invention to distinguish between viable and necrotic tissue within the tumor. Equal cost was assigned to misclassification for both groups. Alternatively, if the costs change, then the cost function could be altered, using normal distribution tables, to obtain the desired results.
If a statistical decision is made on the classification of a pixel based on a conventional single variable approach, then the lower and upper bounds are single values. If the decision is based on the individual bounds of two variables, then decision space is a rectangle. However, the joint distribution of two independent variables is an ellipse with major and minor axes parallel to the variable axes. The 2-variable joint decision space, then, is not the same as the individual rectangular space. Furthermore, if the variables are correlated, then the ellipse is slanted. Therefore, if the decision to segment into a group is based on the variables, individually, then the probability of misclassification rises.
The increased probability of misclassification inherent the univariate decision is demonstrated in
It is clear that univariate analysis of multivariate data increases the probability of both Type I and II errors, and increases the probability of a misclassification of pixels. The probability of misclassification could be reduced with multiple comparison models, but these essentially reduce the size of the univariate confidence bounds to fit into the ellipse to protect against Type I errors, and not Type II errors. Additionally, these methods do not generally take the covariance into account. The MD normalizes the ellipse to the variance/covariance of the data, and transforms the ellipse into a circle to take advantage of any correlation between variables, thus minimizing misclassification that occurs with repeated use of univariate clustering. The reasoning above expands similarly for higher dimensions.
Other methods of clustering are known and have been performed and represented by a statistical comparison means, and method of clustering and comparison may also be achieved by other algorithms, such as Hierarchical Clustering, and Mutually Exclusive analyses. One skilled in statistical analysis of data would be familiar with these and alternate means for use in the instant invention.
Application of the multivariate statistical analysis technique above to a particular type of image will now be discussed. A common feature of most image analysis software packages is the use of histograms to characterize the one-dimensional distribution of pixel intensities. Although this is a useful tool in a region where the intensities are randomly distributed, the histogram ignores the structure of any texture, which has a two-dimensional distribution.
Although there are many different ways of defining texture based on anisotropic patterns, there is no apparent defined texture in the tissue of a tumor. The necrotic areas appear as coarse "blobs" of equal intensities, while viable tissue is finer in structure. Coarse texture can be seen in
Independent measurements of texture also increase the information available to the segmentation algorithm. The cooccurrence matrix provides a measurement of the probability that two pixel intensities within a defined neighborhood have the same intensities, regardless of degree of isotropy or actual intensity. This does not assume a particular pattern.
The probability function of the occurrence of two intensities within a connected array of pixels is known as the cooccurrence matrix. The cooccurrence matrix within a window of pixels is defined as the probability that two pixels, separated by a distance vector Dij, have the intensities Ii and Ij. If the distance vector is limited to separation along the rows and the columns, then the cooccurrence function, h, for a pixel separation of δab and intensities (Ii, Ij), is then defined as
h δab(I i ,I j)=p(x l,m =I i |x l+a,m+b =I j),
for all x within the analysis window. A sample from a small region in
Pixel values in a small region shown in FIG. 7
Example of cooccurrence matrix with pixel separation of 1 along the
rows of the image window shown in FIG. 7
The cooccurrence matrix provides a 2-D probability distribution of pixel intensities and, from this distribution, second-order summary statistics can be calculated that represent the dispersion characteristics. Two measures of the coarseness of the texture of a window of intensities are contrast and entropy. If the individual elements in Table II are designated as h(ij), then the contrast and entropy of window X are defined as follows:
Relatively large values of contrast indicate that similar pixel intensities co-occur at large distances, or that the texture is coarse. Entropy is essentially a measure of the homogeneity of the image intensities. Therefore, if the texture is fine, the probabilities of cooccurrence of intensities are close to the same, and the entropy is relatively large (positive). Coarse textures indicate some structure, and so yield small values of entropy. In the image in
The determination of the cooccurrence matrix may be computationally intensive that requires many calculations for each pixel. If the pixel intensities vary approximately linearly and slowly over the extent of the analysis window, then the image can be low-pass filtered and subsampled, and the unsampled pixel values may be further estimated using bilinear interpolation. This method has been shown to be appropriate for ultrasound image analysis using varying connected window sizes.
The profile of pixel intensities through a typical tumor is shown in
Once the intensity, entropy and contrast are determined, the MD is calculated for each of the k segments. The image is then recharacterized into a vector of pixel values. For an image that is m�n pixels, the length of the vector is N. The X matrix is constructed as
where I is the intensity, E is the entropy value and C is the contrast value for each of N pixels. The value for Σ may be determined from a training set of images where the necrotic and viable tissue is known through histopathology images.
If the training set is defined for the jth group as X0j, then the sampled variance/covariance matrix is defined for the vector X as
The result is the unbiased estimate of the variance based on the sampled data from the training set of data. However, N>>100 and so the use of the sampled variance is essentially equal to the estimate of the population variance. The result can be expressed as the matrix of the covariances of each pairing of variables,
The off-diagonal values for σ are all zero if the measurements are uncorrelated. For three variables, the data space is then symmetric about the multivariate mean, and the confidence region for the multivariate normal distribution is spherical. This is a desirable situation, because the information provided by each of the variables is maximized. During one experiment, the values for Σ for the two groupings, viable and necrotic tissue, were determined from a training set using the Hras tumor, since it was well necrotized. The variance/covariance matrix remained constant throughout the clustering process.
An example using a sample set of pixel intensities for necrotic and viable tissue is shown below. The example centroids, C, and variances, Σ, are
A typical pixel vector of values is [82 17 −320], and the deviations from the centroids are calculated as
The Mahalanobis distances are then calculated as
It is clear from the above, since the higher MD value corresponds to the lower probability, that the pixel should be classified as necrotic tissue.
In one embodiment, these matrix and statistical calculations may be carried out in the processor shown in
The operation of an embodiment of the system and method of the invention will now be described, with reference to
Each one-second cine is preferably examined for a single image that does not have any apparent motion that would corrupt the Doppler signal. Additionally, high PD signals would also be seen as bright gray-scale pixels in the conversion from RGB format to 8-bit gray scale. Motion may be detected using the Doppler signal of the ultrasound coupling gel on the skin of the subject. The absence of a Doppler signal in the gel next to the skin is assumed to indicate that the subject is motion-free, so that the image may be extracted.
Additionally, images may be reviewed for systolic events by watching the Doppler signal from a prominent artery. Systole is the peak of the pulse as the blood is pumped throughout the body. Though the detection of flow through the smaller capillaries is generally not possible due to the wall filter of the apparatus, using the systolic blood flow may be desirable to maximize the detection probability in the small peritumoral arteries that feed the tumor.
Next, the tumor preferably is isolated within the image. This may be done by manually tracing the tumor border for every image or, preferably, one may take advantage of the hyperechogenic tumor border to establish the region of interest for further analysis. At the interface between tumor and skin or between tumor and muscle, the reflections are stronger than the backscatter energy from the interior of the tumor itself so that the tumor edge may be detected using known edge detection, interpolation, and edge-connection techniques, for example.
Preferably, analysis of the ultrasound images in gray scale to identify the tumor and to segment the necrotic areas is conducted using the multivariate statistical techniques discussed above. The texture of the necrotic tissue is preferably measured from the cooccurrence matrices using the calculated entropy and contrast values. The intensity, entropy and contrast measurements may be used in a multivariate segmentation of the image to identify the necrosis. To reduce computation time, the image within the detected tumor boundary may be extracted. The number of groupings or clusters may also be reduced to two—viable tissue and necrosis. The multivariate standardized distance, or Mahalanobis distance (MD), is preferably used in the iterative k-means clustering of the pixels of the tumor, in order to segment viable tissue from necrotic tissue.
Once this is done, colorized pixels may be measured and counted on a dilated tumor border expanded to cover pixels within 750 mm of the thinned border, for example. This measurement is preferably used to obtain a surrogate measurement for intratumoral and peritumoral blood flow, respectively. Image analysis may further be accomplished using Matlab version 6.02, for example, and the Matlab Image Processing Toolbox, version 3.0, may also be used for all standard image-processing functions.
Further, in another aspect of the invention, during a treatment protocol for cancer administered over time, the real-time imaging of a tumor described above may be compared in a "before" and "after" treatment timeframe, in order to evaluate the efficacy of the treatment, or to modify the treatment in response to the relative amounts of necrotic or viable tissue seen in the tumor over time.
Algorithms used in the determination of the cooccurrence matrices and the multivariate k-means clustering can be found at the end of this disclosure.
To confirm that the above-described approach provided sound, credible results, experiments using the method and system were set up in three stages. The first stage was to simulate the necrotic areas using a non-central Rayleigh distribution for the necrosis, and a central Rayleigh noise field representing the viable tissue. The coarse texture was also simulated. The next stage was to use a well-studied tumor cell line to develop the multivariate methodology for the detection of necrotic regions and the methods for measuring the blood flow to the tumor. The tumor cell line used was the Hras tumor because of its rapid growth, well-established necroses and moderate angiogenic properties. The third stage was designed to validate the methods of necrosis classification and tumor blood flow with a different tumor type. The tumor cell line used for these data was the B16-F10 tumor, a melanoma.
A simulation was run comparing the univariate segmentation of the intensity image to the multivariate segmentation using intensity, contrast and entropy. Matlab version 6.12 was used to create a region of necrosis, characterized by a non-central Rayleigh distribution with varying noncentrality parameters. The necrotic region was further given texture by creating an image mask of Poisson distributed pixel intensities, dilating the texture regions to correspond to the approximate size noticed in tumor images from previous studies, and then using a Gaussian mask to soften the edges of each texton. The texture mask gave a coarse texture to the necrotic region that was not obvious. An example of the simulated image of a necrotic region and the Poisson distribution image mask shown in
The simulation varied the distribution of the pixel intensities within the necrosis with values of 0.5, 1, 3 and 10. The noise field was Rayleigh-distributed with the parameter set to 1. The results are shown in
Simulated Necrosis Image Segmentation Using Multivariate
and Univariate Segmentation.
RMS Pixel
RMS Pixel
Intensity of
of Tumor
The texture of the necrosis gave the multivariate segmentation algorithm a distinct advantage in identifying the necrotic areas. Even when the necrotic region had a very low signal-to-noise ratio when compared to the surrounding viable region, the contrast and entropy measurements clearly detected the region. The addition of the texture variables provided the multivariate method with good results, while the univariate use of intensity had poor segmentation results, until the noncentrality parameter increased to 10.
Experiments were conducted on actual tumors in test animals, and the multivariate approach of the invention was correlated with training data obtained from histopathologic analysis.
The principal means of tissue segmentation was a k-means cluster method using a standardized distance and 4 clusters. The initial cluster assignment used prior information from a training set of images to determine the distribution characteristics of the desired tissue types. Assigning pixels to a cluster based on a prior model provides insurance against the problems of randomly assigning two clusters very close to each other. Using a standardized distance enables the clustering to account for very tight or very wide distributions. Initial segment summary statistics, i.e., mean and variance, were determined using an image from the first Hras experiment, and the histopathology analysis of tissue type. While the segment assignments could be done with one iteration using the prior segment summary statistics, there was no guarantee that the segments would be exactly as indicated in the learning model. Therefore, the variances for each segment were assumed to remain constant and the cluster centroids (i.e. arithmetic means) were allowed to update to account for variations in image capture.
Prior to the clustering, the image was filtered using a Gaussian 3�3 filter with a standard deviation of 0.5 pixels. The purpose of the filtering was two-fold: (1) to reduce the effects of high-frequency noise and (2) to shift the pixel intensity distribution from Rayleigh to normal, following the Central Limit Theorem (CLT).
Once the gray-scale image was segmented into the four regions as shown in
The interior of the tumor was segmented from the entire tumor by choosing only those pixels that were inside the border by more than 750 μm, or 15 pixels. This prevented the tumor wall from being included in the calculation, and allowed for some error in the detection of the actual border. The border of the tumor includes some fat, muscle, and other connective tissue that looks much like the necrotic areas within the core.
The texture of the necrotic tissue was measured from the cooccurrence matrices using the calculated entropy and contrast values. The intensity, entropy and contrast measurements were used in a multivariate segmentation of the image to identify the necrosis. The image within the detected tumor boundary was extracted to reduce the computation time. The number of clusters was reduced to two: viable tissue and necrosis. The multivariate standardized distance, or Mahalanobis distance, was calculated, and used in the k-means clustering of the pixels of the tumor.
The images for the tumor study were collected and analyzed using the image correlation with available histopathology determination of necrosis as the measure of similarity. The histopathology images were limited due to limited resources. Similarity of the images analyzed for necrosis was measured as the fraction of analyzed necrotic pixels that were truly necrotic, as identified by the histopathology images. Dissimilarity was measured as the percent of pixels that were analyzed as viable but were within the identified necrotic areas. These numbers were obtained using Matlab version 6.02 and the morphological operators in the Image Processing Toolbox version 3.0. The similarity, PS was calculated as
where the I matrices are binary, with "1" being assigned to necrotic or viable tissue for Inecrotic and Iviable, respectively, and "0" assigned to the alternate tissue type.
Image processing of the images for peritumoral blood flow was accomplished using the Image Processing, Signal Processing and Statistics Toolboxes in Matlab Version 6.12. Functional relationships of peritumoral blood flow were determined using the Nonlinear Regression package in Statistica� and Labstats, an Excel� add-in used for curve analysis.
The apparatus, system, and method of the disclosed invention has wide applicability to a variety of industrial and medical processes, as discussed above, including analysis of images, including ultrasound images, using multivariate statistical techniques related to texture features of the sample. Although the use of the invention in medical research and clinical treatment, including the treatment and evaluation of treatments for cancer are apparent, the invention may also find application in manufacturing processes, including processes where failure analysis or quality control functions are important, for example, semiconductor device manufacturing, or other fields where image analysis to determine surface textures of materials are of interest. Further, although the applicability of the method and apparatus of the invention has been demonstrated for ultrasound image analysis, further applicability may be found in the analysis of images related to electron microscope (SEM, TEM), MRI, or other modalities.
Representative examples of algorithms used during development of the novel method, apparatus and computer-readable mediu are provided below:
function [hContrast,hEntropy,vContrast,vEntropy] =co_occurrence(y)
% determine the co-occurrence matrix for 4-bit gray scale.
if p>1
y = rgb2gray(y);
y = floor(y/16)+1;
yh1 = [y(:,2:n),99*ones(m,1)];
yv1 = [y(2:m,:);99*ones(1,n)];
hEntropy = zeros(length([1:16:n−16]),length([1:16:m−16]));
hContrast = zeros(length([1:16:n−16]),length([1:16:m−16]));
ictr = 0;
neoc = 4;
for i = 1:ncoc:m−ncoc
ictr = ictr+1;
jctr = 0;
for j=1:ncoc:n−ncoc
jctr = jctr+1;
ya = y(i:i+ncoc,j:j+ncoc);
yb = yh1(i:i+ncoc,j:j+ncoc);
for ilvl=1:ncoc
for jlvl=1:ncoc
hij = sum(sum(aa.*bb));
tempcncr(ilvl,jlvl) = hij;
E = −E+hij*log(hij+.01);
C = C+(ilvl−jlvl){circumflex over ( )}2*hij;
hEntropy(ictr,jctr) = E;
vEntropy = zeros(length([1:16:n−16]),length([1:16:m−16]));
vContrast = zeros(length([1:16:n−16]),length([1:16:m−16]));
ictr = 0;
for i = 1:ncoc:m−ncoc
ictr = ictr+1;
jctr = 0;
for j=1:ncoc:n−ncoc
jctr = jctr+1;
ya = y(i:i+ncoc,j:j+ncoc);
yb = yvl(i:i+ncoc,j:j+ncoc);
for ilvl=1:ncoc
for jlvl=1:ncoc
hij = sum(sum(aa.*bb));
tempcncr(ilvl,jlvl) = hij;
E = −E+hij*log(hij+.01);
C = C+(ilvl−jlvl){circumflex over ( )}2*hij;
vEntropy(ictr,jctr) = E;
function [c,Yc] = modelk3(Img, C, E,xborder,yborder);
% Determine the k-means clustering given the a priori estimate of the cluster
%centers and dispersion. The input variable is the multivariate input with n variables
% The initial cluster centers for the three variables are
% [Necrosis Viable]
% Centroids
% Intensities
[75 30]
[300 100]
% Contrast
[14 8]
[1e5 1e5]
% Entropy
[−350 −420]
[3 49]
[imrow, imcol]=size(Img);
% Optimum Values
center = [.6 0.1;
.5 0
.5 1];
varcenter = [.01 .002;
1 1;
1 .2];
% These are the Suboptimum for no prior information on the dispersion
% center = [.6 0.1;
.5 0
.5 1];
% varcenter = [1 1;
.8 1;
1 1];
%Extract only the area inside the tumor walls
[J,BW] = ROIFILL(Img,xborder,yborder);
clear J
ITumorindex = find(BW~=0);
ITumor = double(Img).*BW;
= ITumor(ITumorindex);
= imadjust(double(I1)/256,[ ],[0 1],.8);;
% makes up for the changes in intensity
CTumor = double(C).*BW;
= CTumor(ITumorindex);
= imadjust(C1/max(abs(C1)),stretchlim(C1/max(abs(C1))),[0,1],1);
ETumor = double(E).*BW;
= ETumor(ITumorindex);
= imadjust(−E1/max(abs(E1)),stretchlim(−E 1/max(abs(E1))),[0,1],1);
% E1 = double(E(:)); clear E
nrows = length(I1)
%Do [X−C]inv(V)[X−C]′ for Necrosis and for Viable
% This gives an N�1 vector of standardized distances
center1mat = repmat(center(:,1)′,length(I1),1);
center2mat = repmat(center(:,2)′,length(I1),1);
N = zeros(nrows,1);
V = zeros(nrows,1);
% For Entropy Measurements
MVNcntr = ([I1 , C1 , E1]−center1mat);
MVVcntr = ([I1 , C1 , E1]−center2mat);
% Without Entropy Measurements
MVNcntr = ([I1 , C1 ]−center1mat);
MVVcntr = ([I1 , C1 ]−center2mat);
% Partition the vector to reduce problems with memory in the determination of the
Mahalanobis Distance
N1 = MVNcntr*pinv(diag(varcenter(:,1)));
V1 = MVVcntr*pinv(diag(varcenter(:,2)));
for i = 1:nrows
N(i) = sum(N1(i,:).*MVNcntr(i,:));
V(i) = sum(V1(i,:).*MVVcntr(i,:));
[Y,I] = min(D,[ ],2);
moved = nrows;
while moved>nrows/100
Iold = I;
% Recalculate the centroids of the groups
center(1,1) = mean(I1(I==1));
center(1,2) = mean(I1(I==2));
center(2,1) = mean(C1(I==1));
center(2,2) = mean(C1(I==2));
center(3,1) = mean(E1(I==1));
center(3,2) = mean(E1(I==2));
%Recalculate the groups
center1mat = repmat(center(:,1)′,length(I1),1);
center2mat = repmat(center(:,2)′,length(I1),1);
% For Entropy Measurements
MVNcntr = ([I1, C1, E1]−center1mat);
MVVcntr = ([I1, C1, E1]−center2mat);
% Without Entropy Measurements
MVNcntr = ([I1, C1, ]−center1mat);
MVVcntr = ([I1, C1, ]−center2mat);
% Have to break up the vector
N1 = MVNcntr*pinv(diag(varcenter(:,1)));
V1 = MVVcntr*pinv(diag(varcenter(:,2)));
for i = 1:nrows
N(i) = sum(N1(i,:).*MVNcntr(i,:));
V(i) = sum(V1(i,:).*MVVcntr(i,:));
[Y,I] = min(D,[ ],2);
% Calculate the number of values that moved
moved = sum(I~=Iold);
fprintf(2,′moved = %d\n′,moved);
c = zeros(size(Img));
Yc = zeros(size(Img));
Yc(ITumorindex) = Y;
c(ITumorindex) = I;
% figure(4),pcolor(flipud(reshape(c,irnrow,imcol))),shading flat
Additionally, the disclosure shows and describes only the preferred embodiments of the invention, but as aforementioned, it is to be understood that the invention is capable of use in various other combinations, modifications, and environments, and is capable of changes or modifications within the scope of the inventive concept as expressed herein, commensurate with the above teachings, and/or the skill or knowledge of the relevant art.
Such variations are not to be regarded as a departure from the spirit and scope of the invention, and all such modifications as would be obvious to one skilled in the art are intended to be included within the scope of the following claims. The embodiments described hereinabove are further intended to explain best modes known of practicing the invention and to enable others skilled in the art to utilize the invention in such, or other, embodiments and with the various modifications required by the particular applications or uses of the invention. Accordingly, the description is not intended to limit the invention to the form disclosed herein. Also, it is intended that the appended claims be construed to include alternative embodiments.
|Cited Patent
|Filing date
|Publication date
|Nov 25, 1992
|Oct 25, 1994
|Arch Development Corporation
|Method and system for detection of interval change in temporally sequential chest images
|Jun 23, 2000
|Jun 21, 2005
|Litton Systems, Inc.
|Historical comparison of breast tissue by image processing
|Nov 29, 2001
|Aug 8, 2002
|Foran David J.
|Collaborative diagnostic systems
|Feb 20, 2002
|Aug 29, 2002
|Cytokinetics, Inc., A Delaware Corporation
|Method of characterizing potential therapeutics by determining cell-cell interactions
|Apr 12, 2002
|Nov 28, 2002
|Medispectra, Inc.
|Optical methods and systems for rapid screening of the cervix
|Sep 12, 2002
|Apr 17, 2003
|System and method for quantitative assessment of cancers and their change over time
|Oct 2, 2002
|Apr 8, 2004
|Virtualscopics, Llc
|Method and system for assessment of biomarkers by measurement of response to stimulus
|De Backer, et al., "A competitive elliptical clustering algorithm" Pattern Recognition Letters, vol. 20, No. 11-13, Nov. 1999, pp. 1141-1147.
|F. Cavayas, et al., "Chapter 3.8, Pattern Recognition and Computer Vision for Geographic Data Analysis" in The Handbook of Pattern Recognition and Computer vision (2nd Edition), 1998, pp. 646-649; Fig. 7.
|Gordon et al., "Utilization of Experimental Animal Model for Correlative Multispectral MRI and Pathological Analysis of Brain Tumors," Magnetic Resonance Imaging, vol. 17, No. 10, 1999, pp. 1495-1502.
|Jianchao, et al. "Texture Segmentation of Ultrasound B-Scan Image by Sum and Difference Histograms" Proceedings of the Annual International Conference of the Engineering in Medicare and Biology Society, New York, IEEE, vol. 1. No. 2, November 9, 1989, pp. 417-418.
|M Tuceryan, A.K. Jain: "Chapter 2.1, Texture Analysis" in "The Handbook of Pattern Recognition and Computer Vision (2nd Edition)" C.H. Chen, et al., 1998, pp. 207-248; Fig. 7.
|M. Sonka, J.M Fitzpatrick (EDS.): "Handbook of Medical Imaging, vol. 2. Medical Image Processing and Analysis" 2000, SPIE Press, Bellingham, pp. 100-109, 284, 258.
|McNitt-Gray et al., "Contrast Enhancement Maps for Lung Lesions Imaged on CT," Proc. SPIE vol. 3978: Medical Imaging 2000: Physiology and Function from Multidimensional Images, Feb. 2000, pp. 78-83.
|Subramaniam, et al., "Seafloor characterization using texture" Southeastcon 1993, Proceedings, IEEE, Apr. 4, 1993, p. 8.
|Vaidyananthan et al., "Monitoring Brain Tumor Response to Therapy Using MRI Segmentation," Magnetic Resonance Imaging, vol. 15, No. 13, 1997, pp. 323-334.
|Wong, Sh, et al., "Automatic segmentation of ultrasonic image" Proceedings of the Region Ten Conference, Oct. 19-21, 1993, vol. 3, Oct. 19, 1993, pp. 910-913.
|Xiaoguang Wang, et al., "Using Three-Dimensional Features to Improve Terrain Classifications" Computer Vision and Pattern Recognition, 1997, Proceedings., 1997 IEEE Computer Society Conference, Jun. 17, 1997, pp. 915-920.
|Citing Patent
|Filing date
|Publication date
|Nov 27, 2012
|General Electric Company
|Process and system for analyzing the expression of biomarkers in cells
|Feb 5, 2013
|Olympus Medical Systems Corp.
|Image processing device and image processing method in image processing device for identifying features in an image
|Jul 5, 2011
|Apr 29, 2014
|Bvba Dr. K. Coenegrachts
|Method and device for evaluating evolution of tumoural lesions
|Nov 6, 2009
|Jul 1, 2014
|Olympus Corporation
|Image processing apparatus, computer readable recording medium storing therein image processing program, and image processing method
|Feb 27, 2013
|Sep 2, 2014
|General Electric Company
|Process and system for analyzing the expression of biomarkers in a cell
|Nov 28, 2011
|Dec 8, 2015
|Quantitative Insights, Inc.
|Method, system, software and medium for advanced intelligent image analysis and display of medical images and information
|Mar 14, 2013
|May 10, 2016
|Zuci Realty Llc
|Insight and algorithmic clustering for automated synthesis
|Jul 3, 2003
|Apr 15, 2004
|Marcus Pfister
|In vivo small animal image analysis process and apparatus for image evaluation for in vivo small animal imaging
|Nov 6, 2009
|May 13, 2010
|Olympus Corporation
|Image display device, computer readable storage medium storing image processing program, and image processing method
|Oct 16, 2009
|Apr 21, 2011
|General Electric Company
|Method and system for analyzing the expression of biomarkers in cells in situ in their tissue of origin
|Jul 23, 2010
|Apr 21, 2011
|General Electric Company
|Process and system for analyzing the expression of biomarkers in cells
|Oct 20, 2011
|Olympus Medical Systems Corp.
|Image processing device and image processing method in image processing device
|Jul 26, 2012
|Giger Maryellen L
|Method, system, software and medium for advanced intelligent image analysis and display of medical images and information
|Jul 5, 2011
|Apr 18, 2013
|Kenneth Louis Stefan Coenegrachts
|Method and device for evaluating evolution of tumoural lesions
|Oct 24, 2008
|Mar 18, 2010
|Bionet Co., Ltd
|Automatic contour detection method for ultrasonic diagnosis appartus
|U.S. Classification
|382/128, 382/173
|International Classification
|Cooperative Classification
|European Classification
|Aug 12, 2003
Owner name: PFIZER, INC., NEW YORK
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:RAUNIG, DAVID L.;REEL/FRAME:013869/0507
Effective date: 20030725
|Jan 17, 2011
|Maintenance fee reminder mailed
|Jun 12, 2011
|Lapse for failure to pay maintenance fees
|Aug 2, 2011
|Expired due to failure to pay maintenance fee
Effective date: 20110612


text:
Lakeland Mills® Ottoman
Item # WX2-134198SOLD OUT
You May Also Like
CASTLECREEK 3 Piece Padded Outdoor Rocking Chair SetNon-Member $149.99 Buyer's Club $134.99
CASTLECREEK Complete Patio Dining Set, 6 PiecesNon-Member $129.99 Buyer's Club $116.99
Wooden Purple Martin Bird HouseNon-Member $90.99 Buyer's Club $81.89
24 Piece Interlocking Patio Pavers, 12" x 12" setNon-Member $74.99 Buyer's Club $67.49
CASTLECREEK Deluxe Rolling Garden Seat with Turnbar & Tool StorageNon-Member $99.99 Buyer's Club $89.99
CASTLECREEK Bronze Patio Umbrella Base, Celtic KnotNon-Member $29.99 Buyer's Club $26.99
CASTLECREEK Windmill / Weather StationNon-Member $46.74 Buyer's Club $42.07
CASTLECREEK 9' Market Patio UmbrellaNon-Member $29.99 Buyer's Club $26.99
- Details & Specs
- Reviews
- Articles & Tips
Spend your downtime in this comfy, rustic - looking Lakeland Mills® Cedar Furniture!
At home or at the cabin, when it's time to take time out, do it right! Lakeland Mills® Cedar Furniture is not only beautiful, but its contoured design makes for a comfortable, supportive seat.
Whether you're spending some time with that special someone on the Cedar Log Glider, or enjoying some "me" time relaxing on the Lounge Chair, you'll appreciate every stress-reducing moment you spend on this rustic white cedar furniture. And to complete the set, add the matching Ottoman and put your feet up.More:
- Cedar is naturally rot and insect resistant
- Comes partially assembled, and is easy to finish with common tools
- Made in the USA.
- Assembled measurements: 20" d. x 18" h. x 24" w.
- Weighs approx. 22 lbs.
Order ONLINE Today!

subset: commoncrawl


text:
The Timberwolves erased a 22-point deficit and rallied to down the Nets 107-96 last night in Brooklyn. Minnesota's Andrei Kirlienko filled up the box sheet, scoring 16 points with 10 rebounds, six assists and four blocks in the win. Nikola Pekovic paced Minnesota with 21 points on 9-15 shooting. Deron [[ DARE-un ]] Williams posted 18 points and 14 rebounds in the loss for Brooklyn. Joe Johnson added 17 points and Brook Lopez chipped in 15. This was Brooklyn's first loss at the Barclays Center. Minnesota is now 2-and-1 on the season and hosts Orlando on Wednesday.
Home » Sports » Minnesota Sports » Timberwolves » Wolves Come back Beat Nets
Most Popular
- Discussed
- Viewed
- Prevention, Lane Graves, Disney World, America and Alligators – By Glenn Mollette
- Threat to shoot up rural courthouse in next several days
- New Ulm City Council overrides HPC permit denial
- 4 pounds of pot seized from Mankato home
- Morning News July 25
- Morning Sports July 25
- Paynesville man seriously injured in motorcycle crash
- Church in Renville County burns after being stuck by lightning
- New restaurant opening in St. James soon
- Vikings U.S. Bank Stadium holds open house for public
- Winthrop man accused of raping teen girl
- Brown County drowning
- Threat to shoot up rural courthouse in next several days
- New Ulm Police commander gives more details about threat
- Man charged in high speed chase near New Ulm
- Gaylord woman "crossing a line" with Facebook picture
- New Ulm stabbing causes non-life threatening injuries
- Man who police say was main source of Minnesota meth arrested
- Madelia man seriously injured, passenger killed in crash
- Two Arlington residents charged with Welfare fraud


text:
geophysicsnoun [ U ]uk /ˌdʒiː.əˈfɪz.ɪks/ us /ˌdʒiː.əˈfɪz.ɪks/ specialized
adjective uk /ˌdʒiː.əˈfɪz.ɪ.kəl/ us /ˌdʒiː.əˈfɪz.ɪ.kəl/
geophysical and geological phenomena
(Definition of "geophysics" from the Cambridge Advanced Learner's Dictionary & Thesaurus © Cambridge University Press)
Get a quick, free translation!
Create and share your own word lists and quizzes for free!
Word of the Day
to sew two things together, or to repair something by sewing


text: News of the WeekSCHOLARLY PUBLISHING U.S. Eases the Squeeze on 'Sanctioned' Authors Yudhijit Bhattacharjee*With reporting by Jeffrey Mervis. Science 09 Apr 2004:Vol. 304, Issue 5668, pp. 187DOI: 10.1126/science.304.5668.187a Yudhijit BhattacharjeeWith reporting by Jeffrey Mervis.Find this author on Google Scholar Find this author on PubMed Search for this author on this site Article Figures & Data Info & Metrics eLetters PDF Log in to view full text Username Enter your Sciencemag.org username. Password Enter the password that accompanies your username. Forgot your username or password?Log in Join/Subscribe Purchase Article Activate Member Account Renew Subscription Recommend a subscription to your library Help for librarians


text:
Guide Gear® 2 - person Log Swing
Item # WX2-216736SOLD OUT
You May Also Like
CASTLECREEK 4 Feet Log Swing, 2 SeaterNon-Member $159.99 Buyer's Club $143.99
CASTLECREEK Shale Island Patio Dining Set, 10 PieceNon-Member $479.99 Buyer's Club $431.99
CASTLECREEK Backyard Chicken CoopNon-Member $249.99 Buyer's Club $224.99
CASTLECREEK Retractable AwningNon-Member $259.99 Buyer's Club $233.99
Guide Gear Reversible Outdoor Rug, 9' x 12'Non-Member $49.99 Buyer's Club $44.99
CASTLECREEK Solar Deck Post Cap Lights, 2 PackNon-Member $19.99 Buyer's Club $17.99
Ultrasonic Pest Repeller with StrobeNon-Member $24.99 Buyer's Club $22.49
- Details & Specs
- Reviews
- Articles & Tips
Sway to and fro aboard this rustic, all - wood Guide Gear® 2 - person Log Swing. BIG BUCKS OFF !
Don't swing from the trees like Tarzan, swing from logs! The pleasures of swinging motion come naturally when you sit atop this genuine wood Log Swing. It creates a haven in your yard that invites guests to enjoy a swing or two on lazy summer afternoons.
Sit a spell. Relax. Enjoy the quiet moments. Drink in the natural surroundings as you enjoy years of sturdy comfort. Perfect seating for your porch, patio, deck, garden, backyard or lakeside cabin. We made special arrangements to have these units crafted to our exacting specs.
That means no middleman markup. Simply a dynamite discount deal!Details:
- Logs are New Zealand pine. Seat and back are Chinese fir
- All wood has a protective, weather-resistant clear coat
- Includes Swing and frame, plus all hanging hardware. Frame is 67 x 43 x 63"h. Swing is 48 x 30 x 27 1/ 2 "h. Weighs 77 lbs. Ships in 2 cartons. Assembly required.
Please Note : This is a Heavy / Bulky item. $20.00 for heavy / bulky shipping and handling will be charged in addition to regular shipping and handling.
45 x 18 x 2" thick.
State Color, as available in the Shopping Cart.
Built to Sportsman's Guide's exact specifications, Guide Gear ® stands for dependability, long-lasting quality and unmatched value! Guaranteed!


text:
Your Cart
You have no items in your shopping cart.
Status:Ships in 2-4 Weeks
A true blend of functionality and distinctive looks, the 28d Desk Lamp designed by Omer Arbel is anchored by a grey crochet memory cable that can be coiled up to provide a safe resting place for the light. This innovative design allows the 28d Desk Lamp to comfortably sit on any horizontal surface—desk, table, bookshelf, or beyond. The cable can also be extended, making the 28d a versatile fixture that is easy to reposition when necessary. The 20-watt bulb is contained within an opaque glass shade and is surrounded by a blown glass sphere to optimize light dispersion.
Glass: 6" Dia.
Cable: 15' L with plug
1 lb.
blown glass
Lamp Type: XENON.
Bulbs required: 1 X 20W 12V bi-pin xenon lamp (included)
This lamp has a clear glass diffuser.
Enjoy 15% off your first order by signing up for our newsletter.
Thank you, your promo code will arrive in your inbox shortly.
Hotmail users: add email@example.com to your address book for white listing.
Hello again.
You're already on our list, so you'll continue to hear first about our new products and special features.
Promotion code terms and restrictions
Have a question? Email us at email@example.com, or call us at 800-805-7820 Monday - Friday, 9am - 5pm PT.


text:
The situation was great for access to Lake Constance facilities such as boat trips to other areas. There was also lots of good restaurants all within walking distance.
No tea/coffee facilities in the room which would have been nice. I feel that carparking should be free. Otherwise it was a great stay,
Well-positioned, 3 min from the lake and 2 min from the center. Staff friendly and helpful, room looked great and clean, i found breakfast really good! I recommend it to everyone!
It was said that there is internet in the whole hotel but from the 5th floor where our room was we couldn't get it. As others mentioned it before it is true that the bell rings so often and so loud which can be annoying after a while, but it was not a big problem for us, didn't wake us up.
Location is perfect. Easy walk from there to see that town. It's very close to the station and taxis for transport to airport and/or businesses.
Hotel is clean and comfortable and the location is excellent, close to the station and the lakeside. The staff are friendly and helpful. Excellent value for money.
some added extras rediculously expensive ie. internet, breakfast, bathrobes and towels for sauna etc.
Very central location in Friedrichshafen. Easy to access the lake as well as exhibition center via bus. Bus terminal less than 5 min. walking.
Hotel staff tries to be friendly but seems to be stressed. Booking.com sells rooms at rack rate. I believe booking.com should offer better rate as international booking platform.
Rooms were clean and beds were comfortable.
There was no air circulation, a/c or fan available and there was so much street noise that we could not open the windows. The breakfast area was very small and uncomfortable to move around for a buffet.
Splendid view on the lake. Large rooms, very clean. Great breakfast, again with the view.
The central Location, the view over the Lake and the sound of the lapping Water against the pier.
The Bed !!!! creaking and uncomfortable at every turn of the Body. Definately not a Comfort Room for sleeping.
To be highlighted that we arrived at the Hotel at 03.00 p.m. but nobody was present, no messages and calling to the telephone number of the hotel there was no answer. So we had to return later at 06.00 p.m. when the restaurant waiter open us and give us the key of room asking to keep it with us.
Staff had a positive attitude and were helpful. Good food.
Bathrooms a little too small
I am very like the service of hotel. Everything was very good. Personal of the hotel Eden are very friendly and professional. The room is very comfortable and nice. Location of the hotel is excellent – close to lake, shop and center town.
The hotel staff were very friendly and helpful, giving advice for day trips etc. The breakfast buffet was plenty and the food quality was excellent. The room was very comfortable and clean with private patio. The locations was ideal with Meersburg centre only a short walk away. I would recommend it to everyone.
mini Problems in bathroom !
The food was very good, and you felt free to do as you wanted
Clean and warm with friendly staff.
Could be a little noisy and quite difficult to find even with Sat Nav.
The restaurant is great and one has a good variety of delicious dishes to choose. There is also a lounge and once can relax with friends or chat with colleagues during business travels. The decoration is quite interesting and the furniture is modern. The staff is very friendly and helpful.
Unfortunately, the breakfast was poor. There is no buffet and the guest is served at the table by a plate of cheese and a few cold meat (ham, salami...). The content lacks variety: only one type of bread, one type of juice (orange), one type of water (with gas), etc. On Sunday morning the bread was from last day.
usual ibis,clean,convenient.
overpriced,poor breakfast for the money and dinner was dear and not so good. just an overpriced airport hotel.
Spa was missing
everything was excellent....
Prompt efficeient staff, friendly and helpful.
Room too uncomfortable at night too warm, fan or ac would be a bonus.I Closed the window shutters to keep the sun out during the day while I was away but found they had been raised once the room was attended to while I was at my event.This made the room very uncomfortable upon my return.
Perfect location on Lake Constance. Rooms are very attractive, clean and comfortable, Breakfast was outstanding and dinners were very delicious.
Good location; OK breakfast (although stingy with coffee). Staff was friendly and helpful.
Room was so tiny I could barely move around. No A/C and window was so small that opening it did little to cool off the room.
furniture and style of hotel and room; nice garden views (but no balcony and small window)
mediocre, cold reception; difficult access; bad parking lot; lack of food facilities (no restaurant or room service in the hotel, next door restaurant closed; had to drive 8 km to next crowded town for a simple pizza); no information or tasting/sale of local products although hotel in the middle of a vineyard; internet inordinately expensive...
It had a good and varied breakfast. I found it met my expectations.
Here's why:
We're not the seller.
At KAYAK we don't set the prices, so it's not possible for us to guarantee what other companies are selling.
We compile tonnes of data for you
Our service lets you quickly and easily compare results from hundreds of travel sites at once. In this sense, KAYAK is a search engine and we can't guarantee prices we find.
Why aren't prices accurate 100% of the time?
Prices on airline seats, hotel rooms and car rentals can change frequently. Several people may also be trying to buy the same travel option simultaneously. As a result, you may find, on occasion, that certain prices are no longer available.
If you have more questions or if you have a bad price to report to us,
please send us an
and someone from our team will get back to you promptly.
KAYAK is part of The Priceline Group, the world leader in online travel & related services.


text:
Barrier Meditation, 2015
Installation in California residence
Installation at California residence
We work closely with renowned international artist Ricky Swallow on his bronze casting projects.
Made In L.A., Hammer Museum
Whitney Museum of American Art
The Huntington Library and Museum
For many years -- until the artist realized his enduring dream of opening his own foundry -- we cast hundreds of bronze sculptures for Richard MacDonald.
We collaborated with the Norman Rockwell Sculpture Collection to create and cast editions of several bronze sculptures based upon the artist's covers for the Saturday Evening Post.
Celebrating the female form.
Installation in Disney Ambassador Hotel, Tokyo
We can't show you much more than this yet, but here is the underside of a large sculpture we are casting now.
(casting in progress) Installed in Taipei City, Taiwan
Working with artist Armen Balyan, we created and cast this monumental bronze bust of Congressional Medal of Honor recipient Tibor Rubin.
Several times each year, we work with the Art Center College of Design on special student projects. These projects give the students hands-on experience with the process of bronze casting.


text:
MICE : Meetings, Incentives, Conferences, Events and Exhibitions - Event Managers - All you need to Know !
Event Stuff South Africa ! An agglomeration site for the Meeting, Incentive, Conference, Association, Event and Exhibition industry, gathering together news, trends, insights, best practices and a host of other relevant topics. We'll even include the whacky and wonderful – let Event Stuff South Africa do the hard work and we'll bring it to you !
Our aim is to be topical, informed, provocative, broad-based, insightful, amusing and entertaining with news and views picked from around the globe.
Smart Stuff : "The truth is, technology can't make a bad event great. But the opportunities are truly boundless when we adopt it as part of our creative problem-solving." - Julius Solaris, International Speaker and Author!


text:
Window On Your World 25
Listeners' photos at 5pm
Conor Wilson
Hi chaps
I listen to PM every day in my studio, here in Bristol - I enjoy the leavening humour.
I don't usually have two pieces sitting nearly finished on the bench, but today, at 5, I do - they are off to the London Art Fair in January, after glazing and firing of course.
I noticed that the clock on my camera is wrong - it was taken just seconds after 5, honest.


text:
Dead WP experts,
I am using the Iconic One theme now and already found a lot of things on my "to change" list here in the iconic one area. However, three parts regarding the header have not been solved yet.
I want to
1. increase the font size in the header menu
2. change the header text font color and
3. change the background header color
Any tipps?
I allready digged through the style.css and successfully changed the accent color etc...
My site is: http://www.deathmetalmods.de


text:
Hello! We use cookies to give you the best possible experience on JustGiving. You can change your settings here.
Registered Charity Number NIC 100236
Formed to commemorate Mary Peters Olympic Gold Medal win,
the Mary Peters Trust helps young people achieve excellence in sport by financially assisting them through an Awards programme
The Mary Peters Trust Registered charity number n/a
This year's Chairman's Challenge will see IoD Chair, Paul Terrington take on Dame Mary Peters and you too are invited to take aim and take part.
"Great ambassador for Northern Ireland, deserves greater credit all she has done for this wee country"
Anonymous donated
"I was with my RYANI colleagues in Victoria Square and was so busy chatting overlooked making a donation myself in the collection bucket !"
Sheela Lewis donated
£20.00 + £5.00 Gift Aid
"I was at Victoria Square last Saturday with my RYANI colleagues and was so busy chatting I overlooked making a donation"
Thomas McKenna has raised $2,640.05 so far
I'm On the 3rd April I am running the Paris Marathon for Thomas McKenna because we should support their great work.
Aidan McPeake has raised $381.66 so far
Runway Run, Openwave Belfast
About JustGiving
For fundraisers
For event partners and companies
For charities


text:
I am not responding to bug reports in the comments. Please create a ticket here!
I released my first Android game, would you mind giving it a go?
A simple, yet very powerful system for making physical shops using signs and chests.
It's never been that easy to create shops! With features like shop protection and anti-lag protection, you won't have to worry about your server's economy anymore!
"Just works"! Don't worry about your economy plugin, ChestShop supports all of them with Vault).
The plugin by default submits anonymous stats about your server (player count, plugin/server version and type) to a public statistics website - mcstats.org - in order to provide the developers with usage information. If you wish to opt out, you can do so in the PluginMetrics/config.yml file. It also features an auto-update feature, which can be turned off in the config.
- Frequently Asked Questions
- User-made tutorials
- Download the latest version
- Download older versions
- API for plugin developers
- Source code on GitHub
- Thread on Bukkit forums
Additional modules for ChestShop:
- ChestIcon - Icons for Chests, just like the Showcase plugin!
- Notification module - to use with Notifications plugin
- Towny module - to use with Towny plugin
- Libelula Auto Shop - easy creation of admin shops, using Essentials (User-made)
- ChestShop limiter - limits the number of shops per player (User-made)
- ChestShopLogger - store all shop information in a database (User-made)
- ChestShop Notifier - shows the sales made while the user was offline (User-made)
Using shops
It's easy! Just right-click the sign to buy, left-click to sell. (You can configure the buttons in the config file!)
How to create a shop
Player shop
|Place a chest on the ground
|Place a sign 1 block away from the chest (for example, above it)
|On the sign, write:
|First line will be filled by the plugin automatically
|Price line is a combination of buy and sell prices. You need to have a price indicator (like 'B' - price that will apply to items BOUGHT by PLAYERS - or 'S' - the SELL price) near the price - the prices are separated by a colon (:)
|Some of the possible price combinations: "B 5:5 S", "B 5", "S free"
|You can find out the name of the item you're holding by doing /iteminfo
|For example
|A person can buy 64 diamonds for 10 dollars and sell them back for 5 dollars from AcrobotPL
Admin Shop - shop with unlimited items
Firstly, you've got to be sure that you've the ChestShop.admin permission / you are an OP. To make an Admin Shop, simply put "Admin Shop" (configurable in config.yml) on the first line of the sign.
Here is a list of player-made translations: Translation list
Permission nodes
|shop creation and buying/selling - it contains:
|create all types of shops
|buy any item from shop
|sell any item to shop
|open any shop chest, destroy stores
|open any shop chest, destroy stores, create admin shops
A full list of permission nodes
Useful tips by players
|How not to get bankrupt
|A customisable price table
|Alternative price table
Full plugin description
For full plugin description, visit the plugin's thread on Bukkit forums.
Got an error or suggestion?
Simply submit a ticket here!
Has ChestShop helped you or your server?
If you want to thank the team in a more materialistic way, feel free to donate! :)
A shout-out to my friend
I believe that the music he makes is really stunning - you can decide it for yourself!
MinecraftDotOrg Jul 07, 2016 at 15:36 UTC - 0 likes
Would be really helpful if someone could update this D:
HALOween1232332 Jul 07, 2016 at 14:44 UTC - 0 likes
Seems like the plugin got removed from Spigot, i got a version working in 1.10, but it doesnt support the new spawneggs, will this plugin ever be available again ? Any Information ? Thank you already,
Jalulek Jul 06, 2016 at 17:05 UTC - 0 likes
@frizzbee30: Go
Error. Someone have other link?
faisal2030 Jun 17, 2016 at 12:36 UTC - 0 likes
HELLLP When player want to buy from admin shop they cant why?? but iam using bukkit 1.9.4 help me please write back soon
<<code java>>
Phoenix616 Jun 09, 2016 at 14:15 UTC - 0 likes
The uuids and namechanges worked for me since supprot for them were first implemented? Oo
The only thing that doesn't work properly are the metrics due to the getOnlinePlayers method change in 1.9. But I don't really care about themLast edited Jun 09, 2016 by Phoenix616
Infamousbanana Jun 07, 2016 at 22:15 UTC - 0 likes
hi I'm portuguese and I want to know how I can change the B (Buy) and S (Sell) to C (Comprar) and V (Vender) Please answer me its very important!!
Buy = Comprar Sell = Vender
frizzbee30 Jun 06, 2016 at 15:12 UTC - 1 like
use this https://www.spigotmc.org/resources/1-9-chestshop.19511/ works great :)
KB Craft multiworld server
Website http://kbcraftgaming.enjin.com/home
Members blog site and file share http://kbcraftgaming.com
Server play.kbcraftgaming.com
New players are always given a warm welcome!
KevinEssence May 31, 2016 at 00:47 UTC - 1 like
@Phoenix616: Go
This plugin has been broken since uuids came out, it doesn't support whatsoever. You'll have hundreds of players bitching about names not changing on chest after doing one.
Phoenix616 May 29, 2016 at 22:11 UTC - 0 likes
@trollandia: Go
Plugin still works fine with 1.9, no update needed as far as I can tell.
- Date created
- Aug 25, 2011
- Categories
- Last update
- Jul 12, 2015
- Development stage
- Abandoned
- Language
- deDE
- enUS
- esES
- frFR
- itIT
- plPL
- ptBR
- ruRU
- zhCN
- zhTW
- License
- GNU Lesser General Public License version 2.1 (LGPLv2.1)
- Reverse relationships
- 6
- Downloads
- 4,796,188
- Recent files
- Required dependency
- Vault
- Optional dependency
- Lockette
- Towny Advanced
- WorldGuard


text: Share It MainStoriesHarvard Law School Spotlights Bill Alford for Role in the Global Development Summit Harvard Law School Spotlights Bill Alford for Role in the Global Development Summit Meredith Ciaccia | Jun 21, 2013 Harvard Law School recently recognized Special Olympics International Board member, Bill Alford, for his major role in the Global Development Summit. Alford moderated the Summit's closing panel and was also a key player in facilitating the Global Development Summit Research Symposium which continued the dialogue from the Summit with an emphasis on research, monitoring, and evaluation of best practices. About Me:I am the Global Development and Government Relations Coordinator at Special Olympics International.Related:Read our full story here. Read More Stories Making a Difference in India and NepalDavid HircockThank You, Special Olympics LouisianaJude and Paula BernardAgreement in Italy Will Boost Inclusion and RespectMaureen RabbittBelgium Athlete Addresses Experts at European ForumMartha Jo BraycichMy Coaches Were There For MeCarla SimonBond Between Student and Coach Transforms Life of Special Olympics Swimmer at YMCACarolyn GessnerIn Malta, a Proud Mother Shares Her StoryNatasha MicallefAn "Awesome" Secret: The Universal Power of Unified SportsDavid EvangelistaA Unified Gathering Like No OtherNoha Gaballa"We all need to work together and respect each other to be successful"Maureen RabbittMore Schools in the UK want to Play UnifiedMaureen RabbittDo You Have A Story to Share?Share Your Story More Stories


text:
Are you...?
If yes, click here to login.
Sheehy Nissan of Mechanicsville
Phone: (804) 417-0202 Fax: (804) 417-1034 6564 Mechanicsville Turnpk, Mechanicsville, VA 23111 http://sheehymechanicsvillenissan.com/ ! BBB Business Review on Sheehy Nissan of Mechanicsville is being Updated by BBB !
BBB Business Reviews may not be reproduced for sales or promotional purposes.
View Business Review Inquiries
Sheehy Nissan of Mechanicsville has had their Business Review viewed 12 time(s) over the last 40 days. Click to view where these inquiries originated.
Reason for Rating
BBB rating is based on 13 factors. Get the details about the factors considered.
Based on BBB files, this business has a BBB Rating of No Rating. The reason is as follows:
- This business has no rating at this time because BBB file information is being reviewed and/or updated.