text:
|Publication number
|Publication type
|Application number
|Publication date
|7 Dec 1999
|Filing date
|28 Aug 1998
|Priority date
|4 Sep 1997
|Fee status
|Also published as
|Publication number
|Steven K. Rogers, Philip Amburn, Telford S. Berkey, Randy P. Broussard, Martin P. DeSimio, Jeffrey W. Hoffmeister, Edward M. Ochoa, Thomas P. Rathbun, John E. Rosenstengel
|Original Assignee
|Qualia Computing, Inc.
|Export Citation
|BiBTeX, EndNote, RefMan
|Patent Citations (9), Non-Patent Citations (2), Referenced by (90), Classifications (37), Legal Events (9)
|External Links: USPTO, USPTO Assignment, Espacenet
t.sub.lo =k.sub.lo σ.sub.NN (x,y)+μ.sub.NN (x,y)
t.sub.hi =k.sub.hi σ.sub.NN (x,y)+μ.sub.NN (x,y)
t.sub.lo =k.sub.lo σ.sub.NN (x,y)+μ.sub.NN (x,y)
t.sub.hi =k.sub.hi σ.sub.NN (x,y)+μ.sub.NN (x,y)
t.sub.lo =k.sub.lo σ+μ
t.sub.hi =k.sub.hi σ+μ
t.sub.lo =k.sub.lo σ+μ
t.sub.hi =k.sub.hi σ+μ
t.sub.lo =k.sub.lo σ+μ
t.sub.hi =k.sub.hi σ+μ
t.sub.lo =k.sub.lo σ+μ
t.sub.hi =k.sub.hi σ+μ
This application claims the benefit of U.S. Provisional Application No. 60/057,801, filed Aug. 28, 1997, U.S. Provisional Application No. 60/066,996, filed Nov. 28, 1997, and U.S. Provisional Application No. 60/076,760, filed Mar. 3, 1998, all of which are incorporated herein by reference.
1. Field of the Invention
This invention relates to a method and system for automated detection of clustered microcalcifications from digital images without reduction of radiologist sensitivity.
2. Discussion of Background
Mammography, along with physical examination, is the current procedure of choice for breast cancer screening. Screening mammography has been responsible for an estimated 30 to 35 percent reduction in breast cancer mortality rates. However, in 1996 approximately 185,700 new breast cancer cases were diagnosed and 44,300 women died from this disease. Women have about a 1 in 8 chance of being diagnosed with breast cancer, and 1 in 30 will die of this disease in her lifetime.
Although mammography is a well-studied and standardized methodology, for 10 to 30 percent of women diagnosed with breast cancer, their mammograms were interpreted as negative. Additionally, only 10 to 20 percent of patients referred for biopsy based on mammographic findings prove to have cancer. Further, estimates indicate the malignancies missed by radiologists are evident in two-thirds of the mammograms retrospectively. Missed detections may be attributed to several factors including: poor image quality, improper patient positioning, inaccurate interpretation, fibroglandular tissue obscuration, subtle nature of radiographic findings, eye fatigue, or oversight.
To increase sensitivity, a double reading has been suggested. However, the growing increase in the number of screening mammograms makes this option unlikely. Alternatively, a computer-aided diagnosis (CAD or CADx) system may act as a "second reader" to assist the radiologist in detecting and diagnosing lesions. Several investigators have attempted to analyze mammographic abnormalities with digital computers. However, the known studies are believed to have achieved rates of true-positive detections versus false-positive detections that are undesirably low.
Microcalcifications represent an ideal target for automated detection because subtle microcalcifications are often the first and sometimes the only radiographic findings in early, curable breast cancers, yet individual microcalcifications in a suspicious cluster have a fairly limited range of radiographic appearances. Between 30 and 50 percent of breast carcinomas detected radiographically demonstrate microcalcifications on mammograms, and between 60 and 80 percent of breast carcinomas reveal microcalcifications upon microscopic examination. Any increase in the detection rate of microcalcifications by mammography will lead to further improvements in its efficacy in the detection of early breast cancer.
Although the promise of CAD systems is to increase the ability of physicians to diagnose cancer, the problem is that all CAD systems fail to detect some regions of interest that could be found by a human interpreter. However, human interpreters also miss regions of interest that are subsequently shown to be indicators of cancers. Missing a region that is associated with a cancer is termed a false negative error while associating a normal region with a cancer is termed a false positive error.
It is not yet clear how CAD system outputs are to be incorporated by practicing radiologists into their mammographic analyses. No existing CAD system can claim to find all of the suspicious regions detected by an average radiologist, and they tend to have unacceptably high false positive error rates. However, CAD systems are capable of finding some suspicious regions that may be missed by radiologists.
Accordingly, an object of this invention is to provide a method and system for automated detection of clustered microcalcifications from digital mammograms.
These and other objects are achieved according to the invention by providing a novel method and system for automated detection of clustered microcalcifications from digital mammograms in which a digital mammogram is obtained, parameters necessary for cropping the digital mammogram image are optimized, the digital mammogram is cropped based on the optimized cropping parameters to select breast tissue for further analysis, parameters necessary for detecting clustered microcalcifications are optimized, and clustered microcalcifications in the cropped digital mammogram are detected based on the optimized clustered microcalcification detection parameters.
The detected clustered microcalcifications are then stored as a detections image, the detections image is processed for display, and a computer-aided detection image is produced for review by a radiologist.
The radiologist first reviews the original mammograms and reports a set of suspicious regions of interest, S1. A CAD system, or more particularly, the CAD system of the invention, operates on the original mammogram and reports a second set of suspicious detections or regions of interest, S2. The radiologist then examines the set S2, accepts or rejects members of S2 as suspicious, thus forming a third set of suspicious detections, S3, that is a subset of set S2. The radiologist then creates a fourth set of suspicious detections, S4, that is the union of sets S1 and S2, for subsequent diagnostic workups. CAD system outputs are thereby incorporated with the radiologist's mammographic analysis in a way that optimizes the overall sensitivity of detecting true positive regions of interest.
Other objects and advantages of the invention will be apparent from the following description, the accompanying drawings and the appended claims.
FIG. 1 is a flow diagram illustrating the automated system for the detection of clustered microcalcifications in a digital mammogram;
FIGS. 2 and 3 are flow diagrams illustrating the autocropping method and system of the invention;
FIGS. 4-10 are flow diagrams illustrating in more detail the autocropping method and system of the invention;
FIG. 11 is a flow diagram illustrating in greater detail the clustered microcalcification detector of the invention;
FIG. 12 is a schematic diagram illustrating a 3�3 cross-shaped median filter of the invention;
FIG. 13 is a three-dimensional plot of a Difference of Gaussians (DoG) filter kernel;
FIG. 14 is a cross-sectional view through the center of the DoG filter kernel of FIG. 13;
FIG. 15 is a flow diagram illustrating the global thresholding portion of the microcalcification detection system;
FIG. 16 is a flow diagram illustrating the dual local thresholding of the invention;
FIG. 17 is a flow diagram illustrating combining the results of global and dual-local thresholding;
FIG. 18 is a flow diagram illustrating the sloping local thresholding of the invention;
FIG. 19 is a flow diagram illustrating the clustering method of the invention;
FIG. 20 is a schematic diagram illustrating the clustering method of the invention;
FIG. 21 is a flow diagram illustrating the feature computation process of the invention;
FIG. 22 is a flow diagram illustrating a classifier having one discriminant function per class;
FIG. 23 is a schematic diagram illustrating a multi-layer perceptron neural network for a two-class classifier;
FIG. 24 is a histogram of testing results after detection and classification;
FIG. 25 is a flow diagram illustrating the parameter optimization method of the invention;
FIG. 26 is a plot of a free response receiver operating characteristic curve of the invention before classifying detections;
FIG. 27 is a plot of a free response receiver operating characteristic curve of the invention after classifying detections;
FIG. 28 is a plot of probability density functions showing the relationship between the probabilities of false negative and false positive detections;
FIG. 29 is a plot of probability density functions showing the relationship between the probabilities of true negative and true positive detections;
FIG. 30 is a Venn diagram showing the relationship between radiologist and CAD system detections; and
FIG. 31 is a flow diagram illustrating a method for incorporating computer-aided diagnosis detections with those of a human interpreter for optimal sensitivity.
Referring now to the drawings, wherein like reference numerals designate identical or corresponding parts throughout the several views, and more particularly to FIG. 1 thereof, there is shown a flow diagram illustrating a sequence of steps performed in order to detect the locations of clusters of microcalcifications within a digital mammogram.
In a first step 100, a digital mammogram is obtained using hardware such as digital mammography systems, or by digitizing mammography films using laser or charge-coupled device (CCD) digitizers. In an optimized cropping step 200, a rectangular analysis region containing breast tissue is segmented from the digital mammogram image and a binary mask corresponding to the breast tissue is created for use in later processing steps to decrease the time required for processing the mammogram image. The binary mask is also used to limit detections to areas of the image containing breast tissue.
Clustered microcalcifications are detected in a clustered microcalcification detection step 300. After first filtering the cropped image with a median filter to reduce noise, the image is filtered using an optimized difference of Gaussians (DoG) filter to enhance the microcalcifications. The DoG-filtered image is then subjected to optimized threshold tests to detect potential microcalcifications. The detected microcalcifications are shrunk to single-pixel representations and detections outside of the breast area are removed. The remaining microcalcifications are grouped into clusters. Features are then computed for the clusters. Detected clusters are classified as either suspicious or non-suspicious in a classification step 400.
The parameters used by the autocropping, clustered microcalcification detection, and classification steps 200, 300, 400 are optimized in a parameter-optimizing step 500. The parameters are optimized by parameter-optimizing means that uses a genetic algorithm (GA) so as to maximize the true-positive detection rate while minimizing the false-positive detection rate. Of course, other optimization schemes may be used as well.
The detected clustered microcalcifications are stored in a list of image coordinates. The detection results are processed in a processing step 600 by simply adding an offset to each of the microcalcification coordinates to account for translation of the coordinates incurred as a result of the cropping procedure. Detected clustered microcalcifications are indicated on the digital mammogram by means of rectangles drawn around the clustered microcalcifications in a display step 700. Other indicators may be used such as, for example, arrows pointing to suspected microcalcifications, or ellipses around suspected microcalcifications.
One method of obtaining digital mammograms comprises digitizing radiologic films by means of a laser or charge-coupled device (CCD) scanner. Digital images obtained in this manner typically have a sample spacing of about 100 μm per pixel, with a gray-level resolution of 10 to 12 bits per pixel. In one embodiment of the present invention, radiologic films are scanned using a Model CX812T digitizer manufactured by Radiographic Digital Imaging of Compton, Calif., to produce digital images having 50 μm spacing per pixel and 12 bits of gray-level resolution per pixel. Another possible input source for digital images is a digital mammography unit from Trex Medical Corporation of Danbury, Conn., which has a spatial resolution of about 45 μm per pixel and a gray-level resolution of 14 bits per pixel.
The digital images are stored as digital representations of the original mammogram images on computer-readable storage media. In a preferred embodiment, the digital representations or images are stored on a 2 GB hard drive of a general-purpose computer such as a PC having dual Pentium II� microprocessors running at 200 MHZ, 512 MB of RAM memory, a ViewSonic PT813� monitor, a pointing device, and a Lexmark Optra S1625� printer. The system operates within a Windows NT� operating system.
As may be seen in FIGS. 2 and 3, a digital mammogram image 190 is first cropped to segment an analysis region 296 from the image and produce a binary mask 298 corresponding to breast tissue in the analysis region. Preferably, the cropping is performed automatically, although it could be cropped manually. The image is cropped as a preliminary step because the breast tissue does not cover the whole radiographic film. Focusing the processing of the image on only that portion of the image which breast tissue reduces the time required to process the image. Also, other items appearing on the film, such as labels and patient information, are excluded from consideration, and false-positive indications lying outside of the breast tissue area are eliminated.
Referring to FIGS. 4 through 10, the autocropping process will be described in detail. The image is first subsampled from 50 μm to 400 μm to reduce the amount of data to be processed in step 202. Of course, the image may be downsampled to other resolutions as desired. Not all of the original image data is needed to reliably segment the breast tissue from the remainder of the image. Subsampling every eighth pixel in both the horizontal and vertical directions reduces he amount of data by 64 times. For purposes of segmenting the breast tissue from he rest of the image, the consequent loss of resolution is immaterial.
A white border twenty pixels in width is added around all sides of the subsampled image in step 204. White corresponds to the maximum pixel value possible given the number of bits used to represent each pixel. For images having 12 bits of gray-scale resolution, the maximum gray-scale value is 4095. The bordered image is then thresholded in step 206 with a relatively high threshold value such that most of the breast tissue is guaranteed to be less than the threshold to produce a binary image. In one embodiment of the invention, the threshold is set equal to a predetermined percentage of the gray-scale value of a pixel near the top middle portion of the image. The thresholded image is then inverted, that is, ones become zeroes and zeroes become ones, in step 208. The inverted image is then dilated in step 210. Dilation is a morphological operation in which each pixel in a binary image is turned on, that is, set to a value of one, if any of its neighboring pixels are on. If the pixel is already on, it is left on.
In step 212 the dilated image is cropped to the size of the largest blob. Blobs are contiguous groups of pixels having the value one. This step 212 removes bright borders from the subsampled mammogram representation while ensuring that none of the breast area is reduced. Other techniques that threshold to find the border have a very difficult time dealing with bright areas in the breast adjacent to the border such as, for example, when breast implants are visible in the image. Pixels from the original image, resulting from step 202, corresponding to the locations of the pixels in the cropped blob, are selected for subsequent processing. Note that this is a simple subset of pixels from the input image.
The image from step 212 is histogram equalized in step 214. The average brightness of the image will vary widely from mammogram to mammogram. Moreover, different digitizers having different optical density characteristics are an additional source of variability in brightness levels in the digital representation of the mammogram. The breast mask that is the output of the autocropper is mainly defined by means of a region-growing algorithm that requires a single contrast setting to work properly. However, it has been determined experimentally that a single contrast setting will not work for a wide range of image inputs. Therefore, each image is mapped into a normalized image space using an automatic histogram enhancement process, after which a single contrast setting works well.
First, a histogram of the image is obtained. Typically, most of the data in the breast area will be in the lower histogram bins (corresponding to gray-scale values of about 0-1000), with borders and labels being in the higher bins (corresponding to gray-scale values of about 4000-4095) for 12-bit data. The upper and lower bin values that contain the typical breast data are determined. The lower bin value is the first highest peak encountered when going from the lowest gray-scale value toward the highest gray-scale value. The upper bin is the last zero-value bin encountered when going from the highest gray-scale level toward the lowest gray-scale value. Then the data are reduced to an eight-bit representation and linearly stretched over the range of the data type. For example, values in the lower bins are set to zero. Values of data in the upper bins are set to 255. The rest of the data are then linearly mapped between the lower and upper bins.
After the image has been histogram equalized, the equalized image may be considered to be a matrix. The image matrix is divided into left and right halves, of equal size if possible, and the brighter side is selected in a step 216. The sums of all the pixels in the left and right halves are computed. The sum values are then compared and the side having the greater sum is the brighter side.
Prior to region growing the brighter side, algorithm variables are initialized in step 218. The size of the region-grown mask is preliminarily checked in step 220. If it is large enough, then the mask is acceptable. Otherwise, processing continues to find the mask. The side of the image to be region grown is selected in step 222. In step 224 this region is searched to find its maximum gray-scale value. This maximum value is used to find a pixel to start a region-growing algorithm. Region growing is the process of grouping connected pixels sharing some like characteristic. The choice of characteristic influences the resultant region. The input to a region growing function is a gray-scale image and a starting point to begin growing. The output is a binary image with ones indicating pixels within the grown region, i.e., blobs. Region growing will create a single blob, but that blob may have within it internal holes, that is, pixels that are off. To grow a blob, each of the four nearest neighbors of a pixel of interest are looked at. The contrast ratio is computed for each nearest neighbor pixel. If the contrast ratio is less than a contrast ratio threshold, then the neighbor pixel is set to a one in a binary mask image. Otherwise, the neighbor pixel is set to zero. The region growing algorithm spirals outwardly from the starting or seed pixel, progressively looking at nearest neighbor pixels until done. To those skilled in the art, it is clear that other region growing algorithms may also be applied.
In step 226, region growing begins with the pixel identified from the previous step 224 to produce a binary mask. The size of the mask resulting from step 226 is computed in step 228 and checked in step 230. There may be three points of failure for this approach. First, the brightest point in the search region may be an artifact outside the breast. Therefore, if the resulting mask is not large enough (50 pixels), then the search region is moved closer to the side of the image and searched again. This is repeated three times, each time lowering the contrast value threshold. This corresponds to the path taken through steps 232 and 234. Second, the side selection approach may be in error. Therefore, if a valid breast mask is not found in the first side searched, then the other side of the breast is searched. This corresponds to the path taken through steps 236 and 238. Third, if a valid breast mask is not found on either side, then the whole breast is thresholded and the largest object is taken to be the breast mask in step 240.
Since a constant contrast value is used in the region-growing algorithm, some masks will be too large. Typically, there will be "tails" along the edge of the digitized mammogram image where extra light leaked in while the original mammogram film was being digitized. The tails are reduced by applying a series of erodes and then a series of dilates to the image. Erosion is a morphological operation in which each pixel in a binary image is turned off unless all of its neighbors are on. If the pixel is already off, it is left off. But first, the holes in the mask must be filled in or the multiple erodes may break the mask into disjoint sections. Thus, holes in the mask are closed in step 242 by means of a majority operation. The majority operation is a morphological operation in which each pixel in a binary image is turned on if a majority of its neighboring pixels are on. If the pixel is already on, it is left on.
However, another problem is that some smaller breast masks can not undergo as many erodes as can larger breast masks. Therefore, as a fail-safe measure, the sum of the breast mask is taken before and after the erodes and dilates. If the size is reduced too much (i.e., by more than 50%), the original mask before the morphological operators is used. Thus, a duplicate copy of the mask is made in step 244 before the mask is eroded and dilated in steps 246 and 248, respectively. The size of the resultant mask is then computed in step 250 and compared with the size of the mask from step 242 in step 252. If the new size is less than half the old size, then the duplicate mask, from step 244, is selected in step 254 for subsequent processing. Otherwise, the resultant mask from step 248 is used.
The original image (from step 202) is then cropped to the size of the breast mask just found (either from step 242 or step 248) in step 256. In case the resulting mask is too small for subsequent processing, a crop adjustment is always made in step 258. The adjustment comes in the form of increasing the size of the breast mask bounding box by including additional pixels from the original image in the cropped image.
The cropped image is then automatically histogram enhanced in step 260 as previously described above in connection with step 214. This enhanced image is passed through a loose region growing step 262 to produce a generous mask. This means that the image is subjected to a lower threshold to yield more "on" pixels. This mask is then subjected to hole-closing, eroding, and dilating in steps 264, 266, and 268, respectively, as above, but to a lesser degree.
The same steps described above are repeated one final time in steps 270 through 276, but the crop adjustments are less and the contrast value is increased for a tight region growing step 276. This tight region growing step 276 can afford the higher contrast value since it will be region growing in just the cropped image. This results in a parsimonious estimate of breast tissue. The resulting mask is segmented to find the largest object in step 278 and its bounding box shrunk to just enclose the object in step 280. There may still be some holes in the breast mask. Therefore, after crop adjustments in step 282, the mask is inverted in step 284 and the largest object is found in step 286. This largest object is extracted and then inverted in step 288 to obtain the penultimate mask.
The final mask is obtained by closing holes in the penultimate mask with multiple majority operations and dilations in step 290. The image is then cropped to the size of the resulting mask and the autocropping is complete. An important result from the autocropper is the offset of the cropped image. This is the pixel location in the original image that corresponds to the pixel in the upper left pixel of the cropped image. Keeping track of all the cropping and crop adjustments determines this offset value.
The output of the autocropping process is a rectangular array of pixels representing a binary mask wherein the pixels corresponding to breast tissue are assigned a value of one while the remainder of the pixels are assigned a value of zero. Put another way, the binary mask is a silhouette of the breast made up of ones while the background is made up of zeroes.
Parameters of the autocropper may be optimized to obtain better breast masks. The procedure is described below in the optimization section.
Turning now to FIG. 11, there is seen therein a flow diagram illustrating in greater detail the clustered microcalcification detection system 300 of the invention.
That portion of the digital representation of the mammogram corresponding to the analysis region 296, designated a cropped sub-image 302, produced in the cropping step 200, is first processed to reduce noise in a noise reduction step 310 to reduce digitization noise that contributes to false detections of microcalcifications. The noise-reduced image is then filtered using an optimized target-size-dependent difference of Gaussians (DoG) spatial kernel in step 320 to enhance differences between targets and background, thus creating global and local maxima in the filtered image. The optimized DoG-filtered image is then thresholded in step 340 to segment maxima that represent potential detections of microcalcifications.
The detected maxima are converted to single-pixel coordinate representations in a conversion step 350. The coordinate representations of the detected maxima are compared with the binary mask of the analysis area in a first false-positive removal step 360 to remove false detections outside the breast mask area. The remaining coordinate representations in the analysis area are clustered in a clustering step 370. Features are computed for the remaining clusters in a feature computation step 380 and used to remove non-suspicious detections in a classifying step 400 (FIG. 1). The remaining detections are outputted as detected clustered microcalcifications in an outputting step 600 in the form of cluster coordinates.
Turning now to a more detailed discussion of the steps in the clustered microcalcification detection process, the digital mammogram image is first filtered to reduce noise in the image. Although the main limitation in image quality should be the granularity of the film emulsion, noise is introduced from the process of digitization. This noise may later be detected as a pseudocalcification. In this system, a cross-shaped median filter is used because it is well known to be extremely effective at removing single-pixel noise. The median filter is a non-linear spatial filter that replaces each pixel value with the median of the pixel values within a kernel of chosen size and shape centered at a pixel of interest. Referring to FIG. 12, it may be seen that the cross shape is formed by the set of pixels which include the center pixel and its four nearest neighbors. The cross shape preserves lines and corners better than typical block-shaped median filters and limits the possible substitution to the four nearest neighbors, thereby reducing the potential for edge displacement.
After noise has been reduced, the image is filtered with an optimized DoG kernel to enhance microcalcifications. Filtering is accomplished by convolving the noise-reduced image with the DoG kernel. In an alternative embodiment, filtering is accomplished by first obtaining the fast Fourier transforms (FFTs) of the noise-reduced image and the DoG kernel, then multiplying the FFTs together, and taking the inverse FFT of the result.
The DoG kernel was chosen because neurophysiological experiments provide evidence that the human visual pathway includes a set of "channels" that are spatial frequency selective. Essentially, at each point in the visual field, there are size-tuned filters or masks analyzing an image. The operation of these spatial receptive fields can be approximated closely by a DoG.
The 2-D Gaussian mask is given as: ##EQU1## where c normalizes the sum of mask elements to unity, x and y are horizontal and vertical indices, and σ is the standard deviation. Using Equation 1, the difference of two Gaussians with different σ yields: ##EQU2## It has been shown that when σ2 =1.6 σ1, then the DoG filter's response closely matches the response of human spatial receptive filters. Therefore, with motivation from human physiology, let the ratio of the DoG standard deviation constants be 1:1.6. Then, for a target of size (average width) t pixels, use σ2 =t/2 and, from the rule of thumb, σ1 =σ2 /1.6.
Since microcalcifications typically range from 100 to 300 μm in diameter, potential target sizes for the 50 μm digitized mammograms correspond to 2 to 6 pixels. It has been found that a DoG kernel constructed using an optimization technique for selecting the target size parameter, such as the GA detailed below, has an optimized target size of t=6.01 pixels. The targetsize t will vary depending on such factors as the resolution and scale of the image to be processed. The impulse response of a DoG filter having t=6.01 pixels and σ2 =1.6 σ1 is shown in FIGS. 13 and 14.
Once the noised-reduced cropped image has been DoG filtered to enhance differences between targets and background, the DoG-filtered subimage contains differences in gray levels between potential microcalcifications and background. Although microcalcifications tend to be among the brightest objects in DoG-filtered subimages, they may exist within regions of high average gray levels and thus prove difficult to reliably segment. The thresholding process used in one embodiment of the invention that generally addresses these concerns involves pair-wise pixel "ANDing" of the results of global histogram and locally adaptive thresholding. However, the preferred embodiment of the invention uses sloping local thresholding.
Since targets tend to exist within an image's higher gray levels, then the global threshold may be approximated by finding the level which segments a preselected percentage of the corresponding higher pixel levels in the image histogram. An embodiment of a global thresholding method is illustrated in FIG. 15. Locally adaptive thresholding may be implemented by varying the high and low thresholds based on the local pixel value mean and standard deviation. An embodiment of a dual-local thresholding method is illustrated in FIG. 16.
After computing the image histogram, p(rk), the gray level threshold, g, used to segment a preselected upper fraction,f, of the histogram, is found using: ##EQU3## where rk is the kth gray level, 0≦g≦gmax, and gmax is the maximum gray level in the image.
The locally adaptive thresholds, tlo and thi, are found using
t.sub.lo =k.sub.lo σ.sub.NN (x,y)+μ.sub.NN (x,y) (4)
t.sub.hi k.sub.hi σ.sub.NN (x,y)+μ.sub.NN (x,y) (5)
where klo and khi are used to preselect the multiple of σμNN (x,y), the local standard deviation of gray-level intensities, and μNN (x,y) is the local gray-level mean of the N�N neighborhood centered on the pixel at (x,y) of the DoG-filtered image. Other neighborhood shapes, such as rectangular, circular, and ellipsoidal, may also be used. Pixels whose brightness or gray-level value falls within the threshold interval, that is, tlo <brightness<thi, are set equal to one. Optimization of ƒ, klo, khi, and N is discussed below in connection with the parameter-optimizing process. The results of the global thresholding process may be combined with the results of the local thresholding step by logically ANDing them as shown in FIG. 17. Alternatively, either thresholding method may be used alone.
The preferred thresholding means are illustrated in FIG. 18 wherein it may be seen that an N�N window is centered at a pixel x,y in the input image p(x,y). The mean, μ(x,y), and standard deviation, σ(x,y), of the digital mammogram image pixels under the window are computed. A local threshold value, T(x,y), is computed as:
T(x,y)=A+Bμ(x,y)+Cσ(x,y) (6)
where values for N, A, B, and C are computed during a parameter optimization stage, discussed below. Values for T(x,y) are computed for every x,y location in the image.
The digital mammogram has also been DoG filtered, producing an image d(x,y). Each pixel of the DoG-filtered image d(x,y) is compared to the threshold value T(x,y). Pixels in the locally thresholded image ls (x,y) are set to one where values of the DoG-filtered image are greater than the threshold, and set to zero elsewhere.
The advantage of this novel local sloping thresholding method over prior art thresholding methods is that the threshold is computed from the pixels in a pre-DoG-filtered image rather than from a post-DoG-filtered image. This eliminates the need for background trend correction. In conventional local thresholding, the threshold is computed as:
T(x,y)=Bμ(x,y)+Cσ(x,y) (7)
from the mean and standard deviation of the DoG-filtered image. The problem of using a local threshold computed from the DoG-filtered image is that DoG-filtered images typically have mean values close to zero and standard deviations significantly affected by the presence of targets.
Local thresholds computed from the statistics of the DoG-filtered image suffer from the following adverse effects. First, since the mean value is close to zero, a degree of freedom is lost in the computation of the threshold, which becomes essentially a function of the standard deviation. Second, the absolute brightness of the input image is lost. To keep many spurious detections from occurring, it is desirable to have high thresholds in bright regions. However, the information about the local mean of the input image is not available in the DoG-filtered image. Finally, the standard deviations of DoG-filtered images are increased by detections of targets. This is so because when local bright spots of proper size exist in the original image, large gray-scale values result in the DoG-filtered image. Thus, the presence of targets in a region increases the local standard deviation thereby raising the threshold of that region. The higher threshold reduces the probability of passing a bright spot to subsequent processing stages.
The novel local thresholding method just described solves the above problems by computing thresholds from the input image, which are then applied to the DoG-filtered image. Additionally, the threshold computed here includes an offset term A, which is independent of the local image mean.
After thresholding, detections are converted to single-pixel representations by computing the centroid or center of gravity of groups of contiguous pixels found by the thresholding process. Detections are thus represented as single pixels having a value of logical one while the remaining pixels have a value of logical zero.
False-positive detections outside of the breast area are removed by logically ANDing the binary mask from the autocropper with the single-pixel representations of the detections.
Calcifications associated with malignancies usually occur in clusters and can be extensive. The cluster detection module identifies clusters based on a clustering algorithm as depicted in FIG. 19. Specifically, a suspicious cluster is declared when at least μCsmin or more detected signals are separated by less than a nearest neighbor distance, dnn. Optimization of μCsmin and dnn is discussed below in connection with the parameter optimizing process. FIG. 20 illustrates the clustering process for the case wherein μCsmin =5 and dnn =4.
Additional false-positive clustered microcalcifications are removed by means of a classifier, detailed below. Features are extracted for each of the potential clustered microcalcifications as shown in FIG. 21. The eight features computed for each of the potential clustered microcalcifications in a preferred embodiment are:
1. The larger eigenvalue (λ1) of the covariance matrix of the points in a cluster;
2. The smaller eigenvalue (λ2) of the covariance matrix of the points in a cluster;
3. The ratio of the smaller eigenvalue of the covariance matrix to the larger eigenvalue of the covariance matrix of the points in a cluster. Equivalent to the ratio of the minor axis to the major axis of an ellipse fitted to cover the points in a cluster;
4. Linear density calculated as the number of detected microcalcifications divided by the maximum interpoint distance;
5. Standard deviation of the distances between points in a cluster;
6. Mean minus median of the distances between points in a cluster;
7. Range of points in cluster calculated as maximum interpoint distance minus the minimum interpoint distance; and
8. Density of a cluster calculated as the number of detections divided by the area of a box just large enough to enclose the detections.
Of course, other features could be computed for the potential microcalcification clusters, and the invention is not limited to the number or types of features enumerated herein.
The cluster features are provided as inputs to the classifier, which classifies each potential clustered microcalcification as either suspicious or not suspicious. In practice, the clustered microcalcification detector is only able to locate regions of interest in the digital representation of the original mammogram that may be associated with cancer. In any detector, there is a tradeoff between locating as many potentially suspicious regions as possible versus reducing the number of normal regions falsely detected as being potentially suspicious. CAD systems are designed to provide the largest detection rates possible at the expense of detecting potentially significant numbers of regions that are actually normal. Many of these unwanted detections are removed from consideration by applying pattern recognition techniques.
Pattern recognition is the process of making decisions based on measurements. In this system, regions of interest or detections are located by a detector, and then accepted or rejected for display. The first step in the process is to characterize the detected regions. Toward this end, multiple measurements are computed from each of the detected regions. Each measurement is referred to as a feature. A collection of measurements for a detected region is referred to as a feature vector, wherein each element of the vector represents a feature value. The feature vector is input to a discriminant function.
Referring to FIG. 22, there may be seen therein a classifier having a feature vector x applied to a set of discriminant functions g(x). The classifier shown in FIG. 22 is designed with one discriminant function per class. A discriminant function computes a single value as a function of an input feature vector. Discriminant functions may be learned from training data and implemented in a variety of functional forms. The output of a discriminant function is referred to as a test statistic. Classification is selecting a class according to the discriminant function with the greatest output value. The test statistic is compared to a threshold value. For values of the test statistic above the threshold, the region or detection associated with the feature vector is retained and displayed as potentially suspicious. When the test statistic is below the threshold, the region is not displayed.
Many methods are available for designing discriminant functions. One approach considered for this invention is a class of artificial neural networks. Artificial neural networks require training, whereby the discriminant function is formed with the assistance of labeled training data.
In a preferred embodiment, the classification process is implemented by means of a multi-layer perceptron (MLP) neural network(NN). Of course, other classifier means could be used such as, for example, a statistical quadratric classifier. Only potential clustered microcalcifications classified as suspicious are retained for eventual designation for a radiologist. Alternatively, it may be desirable to iteratively loop between MLP NN analysis of the individual microcalcification detections and the microcalcification clusters.
Referring to FIG. 23, a schematic diagram of an MLP NN may be seen therein. The MLP NN includes a first layer of J hidden layer nodes or perceptrons 410, and one output node or perceptron 420 for each class. The preferred embodiment of the invention uses two output nodes, one each for the class of suspicious detections and the class of non-suspicious detections. Of course, more or fewer classes could be used for classifying clusters of microcalcifications. Each computed feature xi is first multiplied by a weight wij, where i is an index representing the ith feature vector element, and j is an index representing the jth first layer node. The output yj of each first layer perceptron 410 is a nonlinear function of the weighted inputs and is given by: ##EQU4## where d represents the total number of features xi and f(•) is typically a saturating nonlinearity. In this embodiment, f(•)=tanh(•). The first layer or hidden layer node outputs yj are then multiplied by a second layer of weights uj,k and applied to the output layer nodes 420. The output of an output layer node 420 is a nonlinear function of the weighted inputs and is given by: ##EQU5## where k is an index representing the kth output node.
The hyperbolic tangent function is used in a preferred embodiment of the system because it allows the MLP NN to be trained relatively faster as compared to other functions. However, functions other than the hyperbolic tangent may be used to provide the outputs from the perceptrons. For example, linear functions may be used, as well as smoothly varying nonlinear functions, such as the sigmoid function.
The weight values are obtained by training the network. Training consists of repeatedly presenting feature vectors of known class membership as inputs to the network. Weight values are adjusted with a back propagation algorithm to reduce the mean squared error between actual and desired network outputs. Desired outputs of z1 and z2 for a suspicious input are +1 and -1, respectively. Desired outputs of z1 and z2 for non-suspicious inputs are -1 and +1, respectively. Other error metrics and output values may also be used.
In this embodiment of the system, the MLP NN is implemented by means of software running on a general-purpose computer. Altematively, the MLP NN could also be implemented in a hardware configuration by means readily apparent to those with ordinary skill in the art.
After training, each detected clustered microcalcification is classified as either suspicious or not suspicious by means forming the difference z1 -z2, then comparing the difference to a threshold, θ. For values of z1 -z2 greater than or equal to the threshold θ, i.e., z1 -z2 ≧θ, the classifier returns a value of +1 for suspicious clustered microcalcifications, and for values of z1 -z2 <θ, the classifier returns a value of -1 for non-suspicious clustered microcalcifications.
In order to arrive at optimum values for the respective weights, and the number of first layer nodes, the MLP NN was trained with a training set of feature vectors derived from a database of 978 mammogram images.
To develop and test the CAD system of the invention, truth data was first generated. Truth data provides a categorization of the tissue in the digital images as a function of position. Truth data was generated by certified radiologists marking truth boxes over image regions associated with cancer. In addition to the mammogram images, the radiologists also had access to patient histories and pathology reports.
The radiologists identified 57 regions of interest, containing biopsy-confirmed cancers associated with clustered microcalcifications, by means of truth boxes. All 978 images were then processed by the microcalcification detector of the invention to produce a plurality of feature vectors, a subset of which were associated with the 57 truth boxes. Half of the subset feature vectors were randomly chosen, along with about three times as many feature vectors not associated with clustered microcalcifications, to comprise the training set of feature vectors. The MLP NN, having a predetermined number of hidden nodes, was then trained using the training set. The remaining feature vectors were used as a test database to evaluate the performance of the MLP NN after training. Training of the MLP NN was carried out by means of the Levenberg-Marquardt back propagation algorithm.
Alternatively, the MLP NN can be trained with other learning algorithms and may have nonlinearities other than the hyperbolic tangent in either or both layers. In an alternative embodiment with sigmoidal output nodes, the Bayes optimal solution of the problem of classifying clustered microcalcification detections as either suspicious or non-suspicious may be obtained.
In one run of the preferred embodiment during testing, before application of the MLP NN classifier to eliminate false-positive clustered microcalcifications, the detection procedure found about 93% of the true-positive clustered microcalcifications in both the training and test databases while indicating about 10 false-positive clustered microcalcifications per image. It was found that after an MLP NN classifier having 25 first layer nodes was used with the respective optimum weights found during training, 93% of the true-positive detections were retained while 57% of the false-positive detections were successfully removed. Referring to FIG. 24, there may be seen a histogram of the results of testing on the testing database after classification by the MLP NN. Of course, the MLP NN of the invention may be operated with more or fewer first layer nodes as desired.
After the locations of clustered microcalcifications have been determined, they are indicated on the original digitized mammogram image, or a copy of the original image, by drawing rectangular boxes around microcalcifications. Other means for indicating the locations of microcalcifications may be used, such as, for example, placing arrows in the image pointing at detections or drawing ellipses around the detections.
The locations of clustered microcalcifications are passed to the display detections procedure as a list of row and column coordinates of the upper left and lower right pixels bounding each of the clusters. The minimum row and column coordinates and maximum row and column coordinates are computed for each cluster. Bounding boxes defined by the minimum and maximum row and column coordinates are added to the original digitized image, by means well known in the art. The resulting image is then stored as a computer-readable file, displayed on a monitor, or printed as a hard-copy image, as desired.
In one embodiment of the system, the resulting image is saved to a hard disk on a general-purpose computer having dual Pentium II� processors and running a Windows NT� operating system. The resulting image may be viewed on a VGA or SVGA monitor, such as a ViewSonic PT813� monitor, or printed as a hard-copy gray-scale image using a laser printer, such as a Lexmark Optra S1625�. Of course, other hardware elements may be used by those with ordinary skill in the art.
Genetic algorithms (GAs) have been successfully applied to many diverse and difficult optimization problems. A preferred embodiment of this invention uses an implementation of a GA developed by Houck, et al. ("A Genetic Algorithm for Function Optimization," Tech. Rep., NCSU-IE TR 95-09, 1995), which is incorporated by reference herein, to find promising parameter settings. The parameter optimization process of the invention is shown in FIG. 25. This is a novel application of optimization techniques as compared to current computer-aided diagnosis systems require hand tuning by experiment.
GAs search the solution space to maximize a fitness (objective) function by use of simulated evolutionary operators such as mutation and sexual recombination. In this embodiment, the fitness function to be maximized reflects the goals of maximizing the number of true-positive detections while minimizing the number of false-positive detections. GA use requires determination of several issues: objective function design, parameter set representation, population initialization, choice of selection function, choice of genetic operators (reproduction mechanisms) for simulated evolution, and identification of termination criteria.
The design of the objective function is a key factor in the performance of any optimization algorithm. The function optimization problem for detecting clustered microcalcifications may be described as follows: given some finite domain, D, a particular set of cluster detection parameters, x={t, ƒ, klo, khi, N, μCsmin, dnn } where x εD, and an objective function ƒobj : D→, where denotes the set of real numbers, find the x in D that maximizes or minimizes ƒobj. When sloping local thresholding is used in the cluster detector, the parameters N, A, B, and C are optimized. Radiologic imaging systems may be optimized to maximize the TP rate subject to the constraint of minimizing the FP rate. This objective may be recast into the functional form shown in the following equation: ##EQU6## where maximization is the goal. For a particular set of cluster detection parameters, if the minimum acceptable TP rate, TPmin, is exceeded, the objective function returns the negative of the FP rate. Otherwise, if the TP rate falls below TPmin, the objective function returns a constant value, FPpenalty =-10. Other objective functions may also be used.
Since a real-valued GA is an order of magnitude more efficient in CPU time than the binary GA, and provides higher precision with more consistent results across replications, this embodiment of the invention uses a floating-point representation of the GA.
This embodiment also seeds the initial population with some members known beforehand to be in an interesting part of the search space so as to iteratively improve existing solutions. Also, the number of members is limited to twenty so as to reduce the computational cost of evaluating objective functions.
In one embodiment of the invention, normalized geometric ranking is used, as discussed in greater detail in Houck, et al., supra, for the probabilistic selection process used to identify candidates for reproduction. Ranking is less prone to premature convergence caused by individuals that are far above average. The basic idea of ranking is to select solutions for the mating pool based on the relative fitness between solutions. This embodiment also uses the default genetic operation schemes of arithmetic crossover and nonuniform mutation included in Houck, et al.'s GA.
This embodiment continues to search for solutions until the objective function converges. Alternatively, the search could be terminated after a predetermined number of generations. Although termination due to loss of population diversity and/or lack of improvement is efficient when crossover is the primary source of variation in a population, homogeneous populations can be succeeded with better (higher) fitness when using mutation. Crossover refers to generating new members of a population by combining elements from several of the most fit members. This corresponds to keeping solutions in the best part of the search space. Mutation refers to randomly altering elements from the most fit members. This allows the algorithm to exit an area of the search space that may be just a local maximum. Since restarting populations that may have converged proves useful, several iterations of the GA are run until a consistent lack of increase in average fitness is recognized.
Once potentially optimum solutions are found by using the GA, the most fit GA solution may be further optimized by local searches. An alternative embodiment of the invention uses the simplex method to further refine the optimized GA solution.
The autocropping system may also benefit from optimization of its parameters including contrast value, number of erodes, and number of dilates. The method for optimizing the autocropper includes the steps of generating breast masks by hand for some training data, selecting an initial population, and producing breast masks for training data. The method further includes the steps of measuring the percent of overlap of the hand-generated and automatically-generated masks as well as the fraction of autocropped breast tissue outside the hand-generated masks. The method further comprises selecting winning members, generating new members, and iterating in a like manner as described above until a predetermined objective function converges.
In FIGS. 26 and 27, there may be seen therein free response receiver operating characteristic curves for the system of the invention for the outputs of the optimized microcalcification detector and the classifier, respectively. FIG. 26 represents the performance of the optimized detector before classifying detections, while FIG. 27 represents the performance of the system after classifying detections.
Although the GA has been described above in connection with the parameter optimization portion of the preferred embodiment, other optimization techniques are suitable such as, for example, response surface methodology. Of course, processing systems other than those described herein may be optimized by the methods disclosed herein, including the GA.
Performance metrics for detection of suspicious regions associated with cancer are often reported in terms of sensitivity and specificity. Sensitivity measures how well a system finds suspicious regions and is defined as the percentage of suspicious regions detected from the total number of suspicious regions in the cases reviewed. Sensitivity is defined as: ##EQU7## where TP is the number of regions reported as suspicious by a CAD system that are associated with cancers, and FN is the number of regions that are known to be cancerous that are not reported as suspicious. Specificity measures how well the system reports normal regions as normal. Specificity is defined as: ##EQU8## where TN represents regions correctly identified as not suspicious and FP represents regions reported as suspicious that are not cancerous.
Current CAD systems increase specificity by reducing FP. However, FP and TP are coupled quantities. That is, a reduction of FP leads to a reduction of TP. This implies that some of the suspicious regions that could have been detected are missed when the objective is to maintain high specificity.
FIGS. 28 and 29 illustrate relationships between the quantities TP, FP, TN, and FN. A measurement from a screening mammography image is represented by test statistic, x. The probability density function of x is represented by p(x) and the decision threshold is represented by θ. If x is greater than θ, a suspicious region is reported. Areas under the probability density functions represent probabilities of events. From FIG. 28 observe that increasing the threshold reduces the probability of FP decisions. However, observe from FIG. 29 that increasing the threshold simultaneously reduces the probability of TP decisions.
Another metric that exists for CAD systems is positive predictive value (PPV), which is defined as the probability that cancer actually exists when a region of interest is labeled as suspicious. PPV can be calculated from the following equation: ##EQU9## Note that increasing TP or reducing FP increases PPV.
Radiologists and computers find different suspicious regions. FIG. 30 is a Venn diagram depicting a possible distribution of suspicious regions for man and machine detections. Some suspicious regions are found solely by a human interpreter or radiologist, some solely by a CAD system, some are found by both, and some are not found by either.
Referring to FIG. 31, there may be seen a preferred method for incorporating the outputs of a CAD system, and more particularly for the CAD system of the invention, with the observations of a human interpreter of a screening mammography image 10 for optimal sensitivity, wherein a radiologist examines the screening mammography image 10 in a step 20 and reports a set of suspicious regions 30 designated as S1. The CAD system then operates on the image 10 in a step 40 and reports a set of suspicious regions 50 designated as S2. The radiologist then examines set S2 and accepts or rejects members of set S2 as suspicious in a step 60, thereby forming a third set of suspicious regions 70 denoted as S3, which is a subset of S2. The radiologist then creates in a step 80 a set of workup regions 90 denoted as S4 which is the union of sets S1 and S3. The workup regions 90 are then recommended for further examination such as taking additional mammograms with greater resolution, examining the areas of the breast tissue corresponding to the workup regions by means of ultrasound, or performing biopsies of the breast tissue.
While the invention has been described in connection with detecting clustered microcalcifications in mammograms, it should be understood that the methods and systems described herein may also be applicable to other medical images such as chest x-rays.
While the methods herein described and the forms of apparatus a for carrying these methods into effect, constitute preferred embodiments of this invention, it is to be understood that the invention is not limited to these precise methods and forms of apparatus, and that changes may be made in either without departing from the scope of the invention which is defined in the appended claims.
|Cited Patent
|Filing date
|Publication date
|31 Jul 1991
|9 Nov 1993
|Mayo Foundation For Medical Education And Research
|Method and apparatus for diagnosis of breast tumors
|11 Feb 1993
|19 Jul 1994
|E. I. Du Pont De Nemours And Company
|Application of neural networks as an aid in medical diagnosis and general anomaly detection
|28 Apr 1993
|31 Oct 1995
|Arch Development Corporation
|Method and system for differential diagnosis based on clinical and radiological information using artificial neural networks
|21 Jul 1992
|16 Jul 1996
|Arch Development Corporation
|Method for computer-aided detection of clustered microcalcifications from digital mammograms
|22 Jun 1993
|29 Apr 1997
|Mitsubishi Denki Kabushiki Kaisha
|Image processing device for processing grey level images
|1 Dec 1994
|6 May 1997
|University Of Pittsburgh
|Computerized detection of masses and microcalcifications in digital mammograms
|8 Aug 1996
|30 Sep 1997
|Arch Development Corporation
|Computer-aided method for image feature analysis and diagnosis in mammography
|9 Apr 1996
|5 Jan 1999
|Eastman Kodak Company
|Automated method and system for digital image processing of radiologic images utilizing artificial neural networks
|20 Nov 1990
|30 May 1991
|Stereometrix Corporation
|Method for assisting an investigation of a human breast for a malignancy and apparatus for practicing such method
|Carman et al., "Detecting Calcifications and Calcification Clusters in Digitized Mammograms," Digital Mammography '96 (edited by Doi et al.), Elsevier Science B.V., 1996, pp. 253-255, 1996.
|Carman et al., Detecting Calcifications and Calcification Clusters in Digitized Mammograms, Digital Mammography 96 (edited by Doi et al.), Elsevier Science B.V., 1996, pp. 253 255, 1996.
|Citing Patent
|Filing date
|Publication date
|26 Oct 1998
|13 Jun 2000
|R2 Technology, Inc.
|Method and system for computer-aided lesion detection using information from multiple images
|14 Oct 1999
|18 Jul 2000
|Qualia Computing, Inc.
|Method and system for segmenting desired regions in digital mammograms
|14 Oct 1999
|5 Sep 2000
|Qualia Computing, Inc.
|Method and system for combining automated detections from digital mammograms with observed detections of a human interpreter
|14 Oct 1999
|26 Dec 2000
|Qualia Computing, Inc.
|Method and system for segmentation and detection of microcalcifications from digital mammograms
|12 Oct 1999
|20 Mar 2001
|Qualia Computing, Inc.
|Method and system for automated detection of clustered microcalcifications from digital mammograms
|11 Jan 2001
|14 May 2002
|Qualia Computing, Inc.
|Joint optimization of parameters for the detection of clustered microcalcifications in digital mammograms
|8 Jun 1998
|14 May 2002
|Lawrence W. Stark
|Intelligent systems and methods for processing image data based upon anticipated regions of visual interest
|21 Jan 2000
|2 Jul 2002
|Art Advanced Research Technologies, Inc.
|Optical imaging of turbid media with depth-related discrimination
|24 Aug 2001
|29 Apr 2003
|Qualia Computing, Inc.
|Method for combining automated detections from medical images with observed detections of a human interpreter
|8 Jul 1999
|7 Oct 2003
|Perkinelmer Las, Inc.
|Interactive system for analyzing scatter plots
|25 Oct 2002
|18 Nov 2003
|Qualia Computing, Inc.
|Method for combining automated detections from medical images with observed detections of a human interpreter
|15 Apr 2002
|16 Mar 2004
|General Electric Company
|Generalized filtered back-projection reconstruction in digital tomosynthesis
|15 Apr 2002
|20 Apr 2004
|General Electric Company
|Reprojection and backprojection methods and algorithms for implementation thereof
|15 May 2002
|8 Jun 2004
|General Electric Company
|Scatter correction method for non-stationary X-ray acquisitions
|13 Jun 2003
|29 Jun 2004
|Qualia Computing, Inc.
|Method for determining features from detections in a digital image using a bauer-fisher ratio
|1 May 2000
|29 Jun 2004
|Canon Kabushiki Kaisha
|Image extraction apparatus and method
|13 Jun 2003
|13 Jul 2004
|Qualia Computing, Inc.
|Method for analyzing detections in a set of digital images using case based normalcy classification
|23 Jun 2000
|5 Oct 2004
|Icad, Inc.
|Computer aided detection of masses and clustered microcalcifications with single and multiple input image context classification strategies
|15 Apr 2002
|19 Apr 2005
|General Electric Company
|Tomosynthesis X-ray mammogram system and method with automatic drive system
|29 Sep 2003
|29 Nov 2005
|Icad, Inc.
|Use of computer-aided detection system outputs in clinical practice
|28 Sep 2000
|14 Feb 2006
|Koninklijke Philips Electronics N.V.
|Method of processing, and corresponding filtering device
|29 Sep 2000
|21 Feb 2006
|Sharp Laboratories Of America, Inc.
|Video contrast enhancement through partial histogram equalization
|6 Nov 2003
|8 Aug 2006
|Siemens Medical Solutions Usa, Inc.
|System and method for real-time feature sensitivity analysis based on contextual information
|27 Sep 2002
|12 Dec 2006
|General Electric Company
|Method and apparatus for enhancing an image
|21 Nov 2001
|12 Dec 2006
|Inb Vision Ag
|Method of detecting flaws in the structure of a surface
|15 Apr 2002
|15 May 2007
|General Electric Company
|Computer aided detection (CAD) for 3D digital mammography
|6 Nov 2003
|17 Jul 2007
|Siemens Medical Solutions Usa, Inc.
|Computer aided diagnostic assistance for medical imaging
|6 Nov 2002
|20 Nov 2007
|Icad, Inc.
|Information fusion with Bayes networks in computer-aided detection systems
|1 Dec 2003
|20 Nov 2007
|University Of Chicago
|Automated method and system for advanced non-parametric classification of medical images and lesions
|14 Oct 2005
|11 Dec 2007
|Icad, Inc.
|Use of computer-aided detection system outputs in clinical practice
|26 Nov 2003
|25 Dec 2007
|General Electric Company
|Method, system and computer program product for multi-modality registration using virtual cursors
|6 Feb 2006
|3 Jun 2008
|Health Discovery Corporation
|Computer-aided image analysis
|22 Nov 2005
|22 Sep 2009
|Carestream Health, Inc.
|Computer-aided detection of microcalcification clusters
|30 Nov 2004
|8 Dec 2009
|Ge Medical Systems Global Technology Company, Llc
|Method for radiological image processing
|12 Feb 2007
|6 Apr 2010
|General Electric Company
|X-ray device and image-processing method
|15 Apr 2002
|24 Aug 2010
|General Electric Company
|Method and apparatus for providing mammographic image metrics to a clinician
|18 Dec 2006
|16 Nov 2010
|Intel Corporation
|Method and apparatus for local standard deviation based histogram equalization for adaptive contrast enhancement
|14 Aug 2009
|7 Dec 2010
|Carestream Health, Inc.
|Computer aided detection of microcalcification clusters
|26 Jan 2006
|1 Mar 2011
|The United States Of America As Represented By The Secretary Of The Army
|Medical image processing methodology for detection and discrimination of objects in tissue
|11 Sep 2008
|24 Apr 2012
|PACSPLUS Co., Ltd.
|Method, apparatus and program for processing mammographic image
|22 May 2008
|8 May 2012
|International Business Machines Corporation
|Upper troposphere and lower stratosphere wind direction, speed, and turbidity monitoring using digital imaging and motion tracking
|19 Nov 2008
|9 Oct 2012
|Medison Co., Ltd.
|Adaptive image filtering in an ultrasound imaging device
|3 Nov 2009
|25 Dec 2012
|Icad, Inc.
|Systems, computer-readable media, methods, and medical imaging apparatus for the automated detection of suspicious regions of interest in noise normalized X-ray medical imagery
|29 Apr 2011
|7 Oct 2014
|Vucomp, Inc.
|Microcalcification detection classification in radiographic images
|15 Nov 2007
|14 Oct 2014
|At&T Intellectual Property I, L.P.
|System and method of organizing images
|29 Apr 2011
|30 Dec 2014
|Vucomp, Inc.
|Spiculated malignant mass detection and classification in radiographic image
|14 Nov 2014
|17 Feb 2015
|Vucomp, Inc.
|Spiculated malignant mass detection and classification in a radiographic image
|6 Jul 2011
|9 Feb 2016
|Vucomp, Inc.
|Marking system for computer-aided detection of breast abnormalities
|6 Oct 2014
|9 Feb 2016
|Vucomp, Inc.
|Microcalcification detection and classification in radiographic images
|29 Apr 2011
|16 Feb 2016
|Vucomp, Inc.
|Malignant mass detection and classification in radiographic images
|5 Dec 2011
|5 Apr 2016
|Indian Institute Of Technology, Kharagpur
|Method and system to detect the microcalcifications in X-ray images using nonlinear energy operator
|7 Feb 2013
|13 Sep 2016
|Paul Scherrer Institut
|System for non-invasive classification of different types of micro-calcifications in human tissue
|30 Nov 2012
|29 Nov 2016
|Hologic, Inc.
|Method for breast screening in fused mammography
|21 Nov 2001
|13 Jun 2002
|Bernd Michaelis
|Method of detecting flaws in the structure of a surface
|1 May 2000
|29 May 2003
|Masakazu Matsugu
|Image extraction apparatus and method
|21 Feb 2002
|4 Sep 2003
|The University Of Chicago
|Detection of calcifications within a medical image
|15 Apr 2002
|16 Oct 2003
|General Electric Company
|Reprojection and backprojection methods and algorithms for implementation thereof
|15 Apr 2002
|16 Oct 2003
|General Electric Company
|Multi modality X-ray and nuclear medicine mammography imaging system and method
|15 Apr 2002
|16 Oct 2003
|General Electric
|Tomosynthesis X-ray mammogram system and method with automatic drive system
|15 Apr 2002
|16 Oct 2003
|General Electric Company
|Method and apparatus for providing mammographic image metrics to a clinician
|15 Apr 2002
|16 Oct 2003
|General Electric Company
|Computer aided detection (CAD) for 3D digital mammography
|15 May 2002
|20 Nov 2003
|General Electric Company
|Scatter correction method for non-stationary X-ray acquisitions
|14 May 2003
|22 Jan 2004
|Yitzchak Hillman
|Method and a system for combining automated psychiatric profiling from combined input images of brain scans with observed expert and automated interpreter using a neural network
|27 Sep 2002
|1 Apr 2004
|Kaufhold John Patrick
|Method and apparatus for enhancing an image
|6 Nov 2003
|8 Jul 2004
|Siemens Corporate Research Inc.
|System and method for real-time feature sensitivity analysis based on contextual information
|6 Nov 2003
|29 Jul 2004
|Bhavani Duggirala
|Computer aided diagnostic assistance for medical imaging
|1 Dec 2003
|30 Sep 2004
|University Of Chicago
|Automated method and system for advanced non-parametric classification of medical images and lesions
|23 Oct 2003
|28 Apr 2005
|Ajay Kapur
|Systems and methods for viewing an abnormality in different kinds of images
|30 Nov 2004
|23 Jun 2005
|Sylvain Bernard
|Method for radiological image processing
|22 Nov 2005
|25 May 2006
|Zhang Daoxian H
|Automatic image contrast in computer aided diagnosis
|22 Nov 2005
|6 Jul 2006
|Zhang Daoxian H
|Computer aided detection of microcalcification clusters
|14 Oct 2005
|3 Aug 2006
|Rogers Steven K
|Use of computer-aided detection system outputs in clinical practice
|6 Feb 2006
|5 Oct 2006
|Hong Zhang
|Computer-aided image analysis
|12 Feb 2007
|16 Aug 2007
|Muller Serge L W
|X-ray device and image-processing method
|1 May 2006
|13 Sep 2007
|Terry Dolwick
|Attribute based image enhancement and display for medical imaging applications
|18 Dec 2006
|19 Jun 2008
|Shengqi Yan
|Method and apparatus for local standard deviation based histogram equalization for adaptive contrast enhancement
|21 Nov 2007
|10 Jul 2008
|Sony United Kingdom Limited
|Image processing apparatus and method
|19 Nov 2008
|21 May 2009
|Jae Keun Lee
|Adaptive Image Filtering In An Ultrasound Imaging Device
|15 Nov 2007
|21 May 2009
|At & T Labs
|System and method of organizing images
|22 May 2008
|26 Nov 2009
|International Business Machines Corporation
|Upper troposphere and lower stratosphere wind direction, speed, and turbidity monitoring using digital imaging and motion tracking
|14 Aug 2009
|3 Dec 2009
|Zhang Daoxian H
|Computer aided detection of microcalcification clusters
|11 Sep 2008
|11 Mar 2010
|PACSPLUS Co., Ltd.
|Method, apparatus and program for processing mammographic image
|3 Nov 2009
|5 May 2011
|Rosenstengel John E
|Systems, computer-readable media, methods, and medical imaging apparatus for the automated detection of suspicious regions of interest in noise normalized x-ray medical imagery
|29 Apr 2011
|8 Aug 2013
|VuCOMP ,Inc.
|Malignant Mass Detection and Classification in Radiographic Images
|5 Dec 2011
|26 Sep 2013
|Indian Institute Of Technology, Kharagpur
|Method and System to Detect the Microcalcifications in X-Ray Images Using Nonlinear Energy Operator
|7 Feb 2013
|29 Jan 2015
|Paul Scherrer Institut
|System for non-invasive classification of different types of micro-calcifications in human tissue
|29 Apr 2011
|13 Feb 2013
|Malignant mass detection and classification in radiographic images
|31 Jan 2006
|9 Aug 2007
|Sartorius Ag
|Fuel cell system for road vehicles, spacecraft and power stations comprises a fuel cell stack, volume flow adjusting units, a device for measuring the dew point temperatures of fuel and oxidant, a temperature sensor and a current sensor
|13 Aug 2001
|21 Feb 2002
|University Of Maryland, Baltimore County
|Mammography screening to detect and classify microcalcifications
|13 Aug 2001
|6 Sep 2002
|Univ Maryland
|Mammography screening to detect and classify microcalcifications
|U.S. Classification
|382/132, 382/156, 382/270, 382/260
|International Classification
|Cooperative Classification
|European Classification
|9 Dec 1998
Owner name: QUALIA COMPUTING, INC., OHIO
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ROGERS, STEVEN K.;HOFFMEISTER, JEFFREY W.;OCHOA, EDWARD M.;AND OTHERS;REEL/FRAME:009641/0515;SIGNING DATES FROM 19981030 TO 19981102
|6 Mar 2001
|Certificate of correction
|9 Jul 2001
Owner name: GOVERNMENT OF THE UNITED STATES OF AMERICA, THE AS
Free format text: CONFIRMATORY INSTRUMENT NOT AFFECTING TITLE PURSUANT TO EXECUTIVE ORDER 9424;ASSIGNOR:OCHOA, EDWARDM.;REEL/FRAME:011967/0176
Effective date: 20010515
|3 Jun 2003
|Fee payment
Year of fee payment: 4
|26 Jun 2003
|Maintenance fee reminder mailed
|4 Nov 2004
Owner name: QUALIA ACQUISITION CORP., NEW HAMPSHIRE
Free format text: MERGER;ASSIGNOR:QUALIA COMPUTING, INC.;REEL/FRAME:015328/0620
Effective date: 20031231
|29 Dec 2004
Owner name: ICAD, INC., NEW HAMPSHIRE
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:QUALIA ACQUISITION CORP.;REEL/FRAME:015509/0411
Effective date: 20041112
|7 Jun 2007
|Fee payment
Year of fee payment: 8
|7 Jun 2011
|Fee payment
Year of fee payment: 12


text:
Firing Of Hitting Coach Lights a Fire Under The Offense
The Phillies returned home today after a very bad road trip hoping to turn things around against the Rockies. For starters, they fired their hitting coach, Milt Thompson and replaced him with Greg Gross. Thompson is the fall guy for this severely underachieving offense. But would this slap to the face be the wake-up call the Phillies are certainly in need of?
At first, it looked like the same old routine of hitting into double plays and leaving the bases loaded. But in the 5th inning, the offense woke up and figured out how to score some runs. Three straight hits, including an RBI double from Raul Ibanez, got things started.
Ryan Howard even stole a base, as his teammates giggled in the dugout after seeing the big man scamper to second. Jayson Werth hit a sac fly and then Ross Gload, getting a rare start in the outfield, smacked a 2-run homer for his second hit of the night.
The only starter who did not contribute tonight was Wilson Valdez who went 0 for 4. Everyone else, even the pitcher, got a hit. Ibanez finally looks like RAUUUUUULLL again! He had 3 hits and 3 RBI tonight. Ibanez likes hitting 3rd, as I mentioned yesterday. Jimmy Rollins, back in the lead-off spot as I also suggested yesterday, seems to have found his groove as well. J-Roll collected 2 hits and a walk.
On the mound, Roy Halladay returned to his usual dominating form after allowing 6 runs (5 earned) in his last outing. Tonight, he pitched extremely well, allowing only 5 hits and no runs over 8 innings. Oddly enough, the hitting was a bigger story today than 8 shut-out innings from Halladay. Because we expect Halladay to be good; but the offense has been so bad, seeing the Phillies score 6 runs was more exciting.
JC Romero pitched an easy 9th inning and looked very sharp. He hit the leadoff batter, but was lights-out after that. The Phillies won by a score of 6-0 and Halladay finally experienced some run support. Hopefully, this win will be a much needed boost for the team.
The starting pitching for the rest of the weekend has finally been announced. Kyle Kendrick will be brought up to pitch Saturday as Jamie Moyer is placed on the DL. And on Sunday, JA Happ will get his first start since 4/15/10 when an injury sent him to the DL.
As for Moyer, the news is not good. An MRI showed that he has both a sprained ulnar collateral ligament and strained flexor pronator muscle; those injuries are to the left forearm and elbow. The two injuries may require surgery and this may also signal the end of Moyer's career. He is going for another exam early next week and should know for sure after that. Either way, the injuries are more serious than what Happ had this year and Happ missed more than 3 months.
Moyer is in the last year of his contract with the Phillies and it is unlikely that a team would sign a 48-year old pitcher fresh off multiple injuries and possibly surgery. So Moyer may have pitched his last inning. Then again, if anyone can bounce back after such damage, it would be Moyer. I wish him luck :O)
Photo by Jenn Zambri Photography
*Read more about the Phillies at my other home page, Phightin' Phils Phorum in the My Team Rivals network*

subset: commoncrawl


text:
The Broad Institute and Rockefeller University disagree over which scientists should be named as inventors on certain patents involving the gene-editing technology.
Speaking of Science
July 1, 2015|
July 2015's selection of notable quotes
July 1, 2015|
July 2015's selection of notable quotes
July 1, 2015|
Retiring as a professor, and even shutting down your own lab, doesn't necessarily mean quitting research.
July 1, 2015|
Conflict between science and religion continues, with effects on health, politics, and the environment.
June 29, 2015|
Science releases new guidelines for research transparency, hoping to stem the tide of retractions and misconduct.
June 1, 2015|
In the prologue, "Lemurs and the Delights of Fieldwork," author Ian Tattersall shares the paleoanthropological lessons he learned from studying non-human primates in Madagascar.
June 1, 2015|
How to Clone a Mammoth, The Upright Thinkers, The Thirteenth Step, and Humankind
June 1, 2015|
As the science of paleoanthropology developed, human evolutionary trees changed as much as the minds that constructed them.
May 19, 2015|
While some in the scientific and religious communities have declared an end to the tensions between faith and fact, the conflict continues to have impacts on health, politics, and the environment.
May 1, 2015|
A close cooperation between science, politics, and economics has helped to control one of history's most destructive epidemics.
May 1, 2015|
The Genealogy of a Gene, On the Move, The Chimp and the River, and Domesticated


text:
Red Rock Climbing > Black Velvet Peak > Black Velvet Wall, 5.9-5.12a > Dream of Wild Turkeys, 5.10a, 10 Pitches > Climber's Log
Dream of Wild Turkeys, 5.10a, 10 Pitches Climber's Log
[ Sign the Climber's Log ]
|Dow Williams
Date Climbed: May 1, 2015
|Took a 13 year old up there in 3 hrs to the raps.
|Posted Oct 25, 2015 6:12 pm
|One move wonder. Mostly .9 os
|Posted Jul 27, 2013 9:25 pm
Date Climbed: Oct 23, 2010
|Another great route in Black Velvet. We did the upper pitches that many people skip, which were also good.
|Posted May 19, 2012 1:30 pm
|Dow Williams
|Gobbler Direct
Date Climbed: May 21, 2011
|Repeat of both routes with Michael B. Great weekend showing him this combo and climbing Epi the next day, he had done neither. Beautiful night camping near Windy Peak...does not get any better.
|Posted May 23, 2011 7:30 pm
|Dave S
Date Climbed: May 8, 2009
|Climbed with Justin at the start of a three day weekend in Red Rocks. Great climb! The final two 5.9 pitches were definitely worth the extra time!
|Posted May 16, 2009 11:57 pm
|Climbed with Katie in perfect weather.
|Posted Apr 30, 2009 10:11 pm
|Dream of doing this route again
Date Climbed: Sep 26, 2008
|Pitch 4 is the crux and it is short. The slabby moves are well protected with bolts but it is balancey, especially with a mid-size pack. Nothing else on the route was 5.10. Hint: don't place too much pro at the bottom of the pitch.
Our biggest piece was a #3 cam. I'd do the route again without it and I'd probably leave the #2 behind as well. As has been mentioned, there are a lot of bolts and various small cracks all the way up the wall.
|Posted Sep 29, 2008 2:29 pm
|Pitch 4 is definitely the crux.
Date Climbed: May 5, 2008
|This was my second trip up Dream of Wild Turkeys. It had been several years since I did it last and the only thing I remembered from my first experience was the horrible rope drag at the end of some strenuous OW climbing leading to the two hardest crimpy slab moves anywhere.
|Posted May 20, 2008 2:04 pm
|I remember it well
Date Climbed: Nov 22, 1982
|What a great route. Miguel Carmona and I climbed this mega classic in 1982 and bivied on the descent after we got lost at the top of Black Velvet Peak. Great climbing and a great company.
|Posted May 12, 2008 6:27 pm
|Dow Williams
|Dream of Wild Turkeys
Date Climbed: May 5, 2008
|Best weather I have ever had back here. Perfect really, not too windy, great temp. Anybody not doing pitch 9-10 is missing out, pitch 9 is one of the better pitches. Pitch 7 was crimpy, particularly if you have a bad finger. Thought both traverses were tame, I led those....but the corner starting up pitch 3 was full on, had me huffing and puffing a little. Supertopo has the traverse on P 3 at 5.10, but it is not, it is the corner. The traverse has great feet all the way over to the base of 4. The traverse at 5 was no worries either. The short slab move at the top of 4 was the crux of the route, no question. Pitch 2 was fantastic climbing. I like Sour Mash as much.
|Posted May 12, 2008 11:22 am


text:
Tyler Perry's Strangest Role Yet
Yes, the melodrama comedian guy is going a very different route for his next movie. Also today: the mother of all apocalypse yarns is on its way to the big screen, while The Kennedys finds a home on the small.
Tyler Perry, the visionary filmmaker behind such hits as Madea's Gun and Madea's Gonna Shoot You, is taking on a role you probably never thought you'd see him in. No, he's not dating a woman. (BURN.) He won't be dressing up like a crazy old lady either. No, in fact, he'll be a regular human man who has a very serious job. Remember Kiss the Girls and Along Came a Spider, those low-rent thrillers based on popular books that starred Morgan Freeman as Det. Alex Cross, murder investigator extraordinaire? Well, they're continuing the series, just not with Morgan Freeman. Idris Elba, of The Wire and Beyoncé & Ali's Stair-Throwing Hothouse, was originally supposed to play the part, but he backed out. So can you see where I'm going with this? Yes, no joke, Tyler Perry will be playing Alex Cross in a new planned franchise. The first picture will be directed by Rob Cohen, who did xXx and The Mummy: Tomb of the Dragon Emperor. (How big of a mantle do you need to fit all those Oscars, I wonder?) Doesn't that sound great?? Aren't you thrilled?? (In all seriousness, I know Tyler Perry's movies are sort of a third rail in terms of speaking negatively about them, so say what you will about them as Movies, but Mr. Perry is NOT a good actor, so this sounds like an unholy mess.) [THR]
Joaquin Phoenix, done with whatever he was doing with that I Am Here business, is considering jumping back into the regular acting game, and his first choice might be a role in Abraham Lincoln: Vampire Hunter, that movie about the Great Emancipator killing bloodsucking fiends. (No, not Reconstruction-era Carpetbaggers! Actual vampires.) Phoenix would play a mentor to Lincoln, who will be played by Broadway sexypex Benjamin Walker. I think this could be a good choice for ol' Joaquin. It's silly and fun and could be a hit. And maybe it's just dark and weird enough to fuel whatever grim impulses make someone do I Am Here. Or instead of the vampire movie, he could just make a mockumentary sequel called I Was There, about him looking at pictures from the 2006 Oscars and realizing "Hey, I was kind of a big deal once." [Deadline]
Hey Veronica Martians! Here's a big bit of news for you. No, they're not bringing the show back. Let it go. It's just that Kristen Bell has been added to the cast of House of Lies, that Showtime pilot about Don Cheadle as a fast-talking management consultant. Bell will play an ambitious, Ivy League employee. So that's good for her. I think television is where she belongs. And on Showtime she'll get to swear! Of course there's no guarantee this show will actually go to series, but its chances look good. Now if we could just get Percy Daggs III some work, we'd be in business. [EW]
Everyone's into the apocalypse, right? I mean, just look how excited Glenn Beck is now that Egypt is bringing about the end times. Plus there's like Hunger Games and Walking Dead and all that. It's just what people want these days. So, enter one of the most famous apocalypse narratives of the modern age, Stephen King's The Stand. Yes they're making it into a feature film, after 33 years of trying. Of course there was the 1994 miniseries, which I remember watching as a kid and thinking was the most incredibly bleak and exciting and mysterious thing ever but if you watch now is suuuper cheesy and chintzy looking, but this is the big time movie theater adaptation. Could be exciting! Or could be very very bad. Just keep Zak Snyder away from it, huh? Although he did well with apocalypse stuff in Dawn of the Dead, his subsequent movies have not been pleasant. [THR]
ABC has ordered three new pilots. From Desperate Housewives glamorpuss Marc Cherry comes a show called Hallelujah, about a town in Tennessee besieged by a war between good and evil. But then a mysterious stranger comes to town and "brings justice, peace, and faith." So it's basically what? Touched by an Outlaw? Great. Next up is a show called Identity, in which two detectives try to figure out the mystery behind why that movie Identity had the weirdest ending ever. ("Whores don't get a second chance!!") No, it's about identity theft. Lastly, there's Once Upon a Time, a show about a girl going to a town in Maine where fairy tales might be real. It's from the two Lost guys who wrote TRON: Legacy and the script based on the Ouija board. So, it's bound to be good. [Deadline]
Oh, this is too bad for everyone involved. The Kennedys, the History Channel miniseries that they abandoned supposedly because there were too many historical licenses taken, will now air on the ReelzChannel. Mmhm. There's apparently something called the ReelzChannel, based out of
Minnesota Albuquerque, and that's where the show is going to air. Supposedly the channel is available in 60 million homes nationwide, including those with TimeWarner like mine, but... I've certainly never heard of it! You know what I'm going to do, mostly because I feel bad for Barry Pepper and Greg Kinnear? I'm going to find the ReelzChannel and I'm going to watch at least one of these Kennedys episodes. (There are eight in total, and I'm not committing to that many.) I think you should do that too. See if you have ReelzChannel and then get seven of your friends to each watch a part and it'll be like one whole person watched the whole thing! This is civic duty in action, folks. Or, y'know, don't watch it at all. What do I care. [THR]


text:
Planning to watch Showtime's new miniseries, starring Jeremy Irons, The Borgias, which starts this coming Sunday, April 3rd? Well then it's time to brush up on the Borgias and the best place to start is with Mario Puzo's final novel that plunges reader into the colorful tumult of the Italian Renaissance and the family at the center of it, The Family.
The crowning achievement of a truly phenomenal career, Mario Puzo's final novel is a remarkable epic of greed, treachery, sin, and power beyond mortal imagination. It is a journey to a different time and place—when the Church held the ultimate authority and ambition was cloaked in robes of the richest velvet. But most of all, it is the spellbinding story of a father and his children, bonded by blood, devotion, and dark purpose, who would descend into hell to rise to challenge the heavens; a family whose name is forever emblazoned in the annals of infamy . . . Borgia!
This trailer has me sold. The Borgias debuts Sunday, April 3rd at 9pm ET. Browse inside The Family now.


text:
not missy, Missy is tattooed on her right buttock.. if was her..
I think it's Victoria Butts, a lot of European porn chicks work until their parents find out then quit.
Hey guys....Her name is Victoria Butts and you can see her in Sara Jay & Victoria Butts - Phat Azz White Girls #20 video....She was deported back to Russia because her Visa expired....
Who is that woman? Can anyone tell me her name please.....She is sexy and hot as hell!! I would love to fuck her big whooty ass!!!
thats NOT Missy Monroe...different ass cheek tattoo
thats not missy monroe. different tattoos on ass cheek
as much as I say wesley pipes should shut his damn mouth this shit here was funny as hell.
she got a big ass
good for her wearing such a gay costume
I use to have this dvd, my fucking friends must have taken it!!! Sum-bitches!!
Is he breaking into a safe with a butterknife? Don't quit the day job man.
I think missy monroe has missy tatoo'd on her ass, or she does in the porn of her i have ^-^
That is NOT missy monroe, she is a complete different person. She hasn't done any major studio work yet, she does low budget stuff like this, kind of like an underground musician.
i think its missy monroe
now if we can only get that to be real crime punishment
Huge Ass
When he's being sat on he sounds like he's itching to shout "I'm the juggernaut, bitch!"
I didn't mean to, but I burst out laughing at the beginning of this.
Related Videos
Recognize someone?
Help improve Spankwire by telling us who
the stars are in this video.
Thank You!
Pornstars are added once we get
enough submissions with the same name.
Privacy Policy | Terms of Use | 18 U.S.C. 2257 | Contact | DMCA | RTA | FAQ | Feedback | Mobile version | Sitemap


text:
Late last year AMD and ATI merged
into one company, forming not just a bigger entity but creating a very specific
roadmap altogether. AMD's grand scheme with ATI is to develop a single chip
handling both general
purpose computing as well as graphics. The Fusion project
as AMD calls it is this very goal.
Today, NVIDIA reveals that it is not behind when it comes to general purpose
GPU, or GPGPU, computing. Earlier this year the company announced its complete unified device
architecture, or CUDA, Technology, which laid the groundwork for GPGPU
programming for NVIDIA GPUs. CUDA Technology directly competes with AMD's Stream Computing
The GPGPU product lineup will be known as Tesla. Tesla is a top to bottom
product lineup consisting of internal PCIe cards and external high-performance
computing, or HPC, systems – Tesla C870, S870 and D870.
The internal PCIe solution consists of an output-less GeForce 8-series
based card on a PCIe x16 card. The Tesla D870 is NVIDIA's only internal GPGPU
card for desktops. The GPGPU still requires two external PCIe power connectors
and consumes up to 170-watts of power at maximum. NVIDIA claims the Tesla D870
delivers 518 Gigaflops of GPGPU processing power.
Last year, the company announced a highly integrated graphics sub-system named QuadroPlex. Using
a number of GPUs in a tightly integrated system, the QuadroPlex family of
machines accelerated 3D rendering and graphics work. QuadroPlex became the
stepping-stone for the new Tesla C870.
The Tesla C870 GPGPU server packs two GeForce 8-series GPUs in an external
system with packaging similar to the QuadroPlex. The GPGPU delivers one Teraflop of GPGPU computing power while consuming up to 550-watts of power.
Finally, the Tesla S870 comes equipped with four GeForce 8-series GPUs and offers
up to two Teraflops of computing power. The Tesla S870 consumes up to 800-watts of power and fits into a stackable 1U chassis.
Tesla C870 and S870 systems connect to workstation systems via an external PCIe
Gen2 x16 interconnect. The machines contain PCIe switches and can be
daisy-chained with more systems. As with the Tesla D870 GPU card, the Tesla C870 and S870 systems lack output capabilities. Theoretically, customers can
purchase multiple Tesla GPGPU systems and chain them up for big increases in
NVIDIA designed the new Tesla family for everything from graphics rending to
medical research and data farming. At the core level, GPUs are far more
efficient at dealing with parallel computing than general-purpose processors.
This makes Tesla very powerful for cluster-type applications.
The Tesla S870, D870 and C870 carry an MSRP of $12,000, $1,499 and $7,500,


text:
|Forum Rules
|Firearms Safety
|Firearms Photos
|Lost Password
|Email Changes
|Today's Posts
|Mark Forums Read
|Thread Tools
|Search this Thread
|September 16, 2012, 04:56 PM
Senior Member
Join Date: October 13, 2008
Location: Northern Virginia
Posts: 407
First test drive with MagnetoSpeed chronograph
I received my MagnetoSpeed this week and was able to take it out to the range today for the first time to try it out.
It's extremely easy to use, strap the bayonet sensor to the barrel, tighten it, look down the bore to make sure the bullet has a clear path and plug it into the unit. Took about 90 seconds. I was shooting 6.5x55 from a CG-63 which tapers from 19mm to 17.5mm at the muzzle and I didn't need to use the extra spacer for the bullet to clear the sensor. The strap held tight for the 50 rounds I fired with the bayonet sensor attached, didn't need to retighten once.
It registered each shot, I was shooting my handloads which were Sierra Match Kings and Remington Core-Lokt both 140gr. I didn't have to adjust the unit's settings, the default setting and sensitivity worked just fine. Using the unit is simple, saving each string is simply a matter of pushing the button a couple of times. Once to get to the menu, scroll down to archive and press again; next numbered string appears a couple seconds after that automatically.
As to how accurate it is compared to other units, I can't answer that, this is my first chronograph. Reviews claim it's within less a % of an Oehler, unless proved otherwise I'll take their word for it. I didn't bring my other rounds that I have chronographed on someone else's unit to compare results with, I've never chronographed my 6.5x55 rounds and wanted to do them first.
The reason I bought this is convenience. I live in Northern VA, so the ranges are crowded, usually indoors, and you frequently don't have the luxury of time. I frequent the NRA range a lot, so I needed something that I can set-up quick and is easy to use and this fits the bill perfectly. I also only have to carry a 14x7x2" box, no tripod, and no worrying about lighting. I've only used it once so far, but it did what the makers said it would do and I was very happy with my experience with it today...and will be using it again on the next trip to the range. If I have any problems, I'll update the post
I've been reloading for about 2 years now, so it's about time I got and used a chronograph.
Edit: I almost forgot to mention the effects of the bayonet on accuracy. The rounds I shot were batches I'd shot before, so I knew what they should do. The only thing I could tell was my POI was about an 1" higher than it would be without the sensor attached. Groupings were about as tight as I normally get. At 50 yards (I was at the NRA range and that's maximum distance) most of my groupings were around 1", some a little tighter at 1/2", about 3 groupings approached 2" and I'm pretty certain it's the shooter fault there...I'm not a competition shooter. But on the CG-63's barrel, it just seemed to move the POI a little higher. The bayonet sensor is very light.
"And remember, Abraham Lincoln didn't die in vain, he died in Washington D.C." - Firesign Theatre
Last edited by bumnote; September 16, 2012 at 05:25 PM.
|Thread Tools
|Search this Thread


text:
Welcome to the Bedford Academy. We purchased the transformed mansion in 1996 from its former owners the Toronto Ladies' Club and over the years we have restored and transformed this gracious property into what has become a favourite destination in the Annex and Bloor and Bedford area . Our goal and mandate over the years has been to "resist the ordinary" and provide our customers with a casual, friendly and eclectic environment. Our goal is to deliver more than we promise and provide service and food with a twist as we strive to be a little different in what we offer in our gastropub fair.
The Academy is not only a pub, a bistro and a bar; it combines a little of all of these so there is something for everyone. The reason why both of our locals, the "Bedford Academy" and the "bar Wellington", have become such popular destinations is our people. They are service inspired, good people that share in our level of commitment to professionalism and they are fundamental to us making your visit a success.
The Bedford Academy "pub-bistro-bar" has been the favourite destination for the famous and not so famous. Movie stars filming in Toronto, musicians and regular folks have all made it home, which is part of what makes the "Academy" such a special place.


text:
Recycling Initiatives on Campus
Story written by Kelsey Barnett '12
Members of the Trinity Campus Ecological Stewardship Advisory Group (CESAG) recently met in the Heritage Science Center to see how well students were recycling. The demonstration was an effort to improve the College's recycling habits.
Members of CESAG asked custodians to refrain from throwing out the trash and recycle bins in the Heritage Science Center for one week. They then opened up all the bags and bins to see what students think should be recycled and what should be thrown away. After observing each pile, the members sorted them properly.
The initial response was that the bins were much better than when they did this experiment at the Bootsma Bookstore Café last semester. Another thought was that the volume of each pile was similar.
The most common mistake was that plastic bottles, cans, and paper are being thrown away and should be recycled while Styrofoam and other biodegradable items should not be recycled but thrown away.
When CESAG Chair Dr. Thomas Roose, associate professor of physics and science education, was asked what he would grade the Science Center, he said a low A to a high B for what is being recycled, but a C for the garbage pile. Recalling the same experiment at the BBC last fall, Roose said he would have given both piles a D or lower.
"Our desire should be to be good stewards of God's creation," said Taylor Bandstra '13 of Pella, Iowa, and member of CESAG. As a member of CESAG, he says the vision statement of this group explains why the recycling habits of the college must improve.
It reads: "A Reformed worldview mandates that humans be good stewards of God's creation. Trinity seeks to teach this worldview and to practice good stewardship, working to care for and redeem the earth."
"God commands us to care for his creation, no buts or ifs about it," added Bandstra. "That is the only motivation we should need."
CESAG challenges Trinity students, faculty, and staff to be more conscious of what they throw away and to follow signs explaining what is recyclable.
CESAG addresses questions and issue concerning recycling and their mission here: tcc.trnty.edu/about/habitat-for-life.pdf.
{loadposition vid-cesag}