US20030129648A1 - Automated DNA array image segmentation and analysis - Google Patents

Automated DNA array image segmentation and analysis Download PDF

Info

Publication number
US20030129648A1
US20030129648A1 US10/375,289 US37528903A US2003129648A1 US 20030129648 A1 US20030129648 A1 US 20030129648A1 US 37528903 A US37528903 A US 37528903A US 2003129648 A1 US2003129648 A1 US 2003129648A1
Authority
US
United States
Prior art keywords
image
dna
grid
frame
grid point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
US10/375,289
Inventor
Soheil Shams
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Biodiscovery Inc
Original Assignee
Biodiscovery Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Biodiscovery Inc filed Critical Biodiscovery Inc
Priority to US10/375,289 priority Critical patent/US20030129648A1/en
Publication of US20030129648A1 publication Critical patent/US20030129648A1/en
Assigned to BIODISCOVERY, INC. reassignment BIODISCOVERY, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: SHAMS, SOHEIL
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10064Fluorescence image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30072Microarray; Biochip, DNA array; Well plate
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S707/00Data processing: database and file management or data structures
    • Y10S707/99931Database or file accessing
    • Y10S707/99933Query processing, i.e. searching
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S707/00Data processing: database and file management or data structures
    • Y10S707/99941Database schema or data structure
    • Y10S707/99944Object-oriented database structure
    • Y10S707/99945Object-oriented database structure processing

Definitions

  • the present invention relates to DNA array image analysis, and, in particular, to automatically segmenting DNA array images into individual DNA spot images for quantification.
  • Cellular behavior is primarily dictated by the selective expression of a subset of genes. Normal growth and differentiation depends on the appropriate genes being expressed in a desired context. Various disease states alter the normal expression of genes as compared to normal tissue. For example, malignant transformation of cancer tissues involves or induces altered gene expression. Through signal transduction cascades and transcriptional networks, alterations of one gene can impact a large number of genes and result in global effects on cell behavior. Regulation of translation and post-transcriptional modification play significant roles, but, invariably, signal transduction pathways lead to the nucleus and changes in gene transcription.
  • One such technique includes use of ordered micro-arrays that allow two color fluorescence detection of hybridization signals. Individual DNA targets are arrayed on a small glass surface and hybridized with fluorescently labeled heterogeneous DNA probes derived from cDNA. The amount of fluorescence at each DNA spot correlates with the abundance of that DNA fragment in the probe mixture.
  • genomic micro-arrays are useful for genomic mapping, genomic ploidy measurements and as hybridization targets for genomic mismatch scanning. Such techniques require rapid quantitative analysis of fluorescent hybridization for hundreds to tens of thousands of DNA spots. As such, there is a severe bottleneck in gene expression data collection due to inadequate methods for processing of individual DNA spot images for determining the quantitative fluorescent hybridization levels.
  • Some existing methods include manual processing of DNA spot images using a generic image processing tool, such as NIH image.
  • a generic image processing tool such as NIH image.
  • NIH image a user visually locates each DNA spot image in a micro-array image, and moves a display pointer to each spot image, and manually defines a small area around the spot image.
  • the image processing tool then reports image intensity values within the small area.
  • the user then manually records the intensity values and continues this process for other visually located DNA spot images in the micro-array image.
  • some existing methods automate the process of locating DNA spot images from micro-array images and quantifying corresponding expression values.
  • Such methods utilize a computer to manually position a cell grid on an area of the micro-array image containing an array of DNA spot images.
  • the grid can be resized and individual columns and rows of the grid can be manually adjusted to better fit the arrayed pattern of DNA spot images.
  • the grid position is then used by the computer to quantify the expression values using the intensity levels at each cell in the grid.
  • the grid used in such methods is either completely fixed in shape, or has limited global flexibility (e.g., resizing and rotating the entire grid).
  • DNA spot image signals derived from the micro-arrays are susceptible to surface noise and laser reflection, due to surface dust. And, nonspecific DNA binding to the silanized surface occurs in a non-uniform pattern creating a varying background of fluorescence over the surface.
  • Existing methods are unable to cope with irregular micro-array pattern, search for DNA image spots, and accurately quantify specific signals while accounting for the local background.
  • Another disadvantage of existing systems is their inability to display micro-array image pixel intensities, corresponding to gene expression values in related DNA spots for example, in an intuitive manner. As such, the user cannot easily determine gene properties in such DNA spots.
  • the present invention satisfies these needs.
  • the present invention provides a method for segmentation of a frame of image information including a plurality of spaced DNA spot images corresponding to a plurality of DNA spots, the image information including image intensity level and intra frame position information corresponding to said DNA spots.
  • the method of the present invention comprises the steps of: (a) transferring the frame of image information into a memory device; (b) selecting a set of image information within said frame including a selected set of the DNA spot images; (c) generating a grid in said memory device, the grid including a plurality of spaced grid points corresponding to said selected DNA spot images, each grid point including position information indicating the position of the grid point within said frame; and (d) modifying a current position of at least one grid point corresponding to a spot image to shift the grid point toward the corresponding spot image. Step (d) can be repeated for said grid point and for all the grid points of the grid.
  • the step of modifying said current position includes: (i) selecting a first bounding area in the frame around the current position of the grid point; (ii) generating a first position update including position information for updating a current position of said grid point to a first new position within the first bounding area, the location of said first new position relative to said current position being a function of intensity level of at least a portion of the image within the first bounding area; (iii) generating a second position update including position information for updating said current position to a second new position in the frame, said second new position being in a geometric arrangement with the position of one or more grid points around said grid point; and (iv) updating said current position with the position information of the first and the second position updates, thereby shifting said grid point toward the corresponding spot image.
  • the DNA spot images can be in a substantially two dimensional array arrangement, and generating the grid can include generating a two dimensional array of grid points spaced according to a predetermined criteria.
  • the method can further include the step of segmenting the selected set of image information by selecting at least one image segment defining a segment area around a grid point and including a spot image with minimum distance from said grid point, said segment area being a function of the spacing between said grid point and one or more neighboring grid points.
  • the selected set of image information can further be segmented into a plurality of image segments corresponding to the plurality of grid points in the grid, each image segment defining a segment area around a corresponding grid point and including a corresponding spot image with minimum distance from said grid point, said segment area being a function of the spacing between said grid point and one or more neighboring grid points, wherein each spot image is contained in a corresponding image segment.
  • the method of the present invention can further include quantifying at least a portion of image information in said image segment to obtain image characteristic values for said image segment.
  • the image characteristic values can include DNA information for a DNA spot corresponding to the DNA spot image in said image segment, said DNA information including gene expression values.
  • the present invention provides a method of displaying image information corresponding to a plurality of DNA spot images of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels.
  • the display method includes the steps of: (a) for each DNA spot image: (1) extracting said background and signal intensity levels from the image characteristic values for the spot image, and (2) determining difference values between the background intensity levels and signal intensity levels; and (b) for each DNA spot: (1) relating the corresponding difference values to a range of graphic values, (2) selecting a graphic value for each difference value, and (3) displaying the selected graphic values.
  • the steps of relating and selecting can include associating each difference value to a segment of a pie chart having multiple segments, and the step of displaying the selected graphic values can include displaying said segments as a pie chart.
  • the area of each segment of each pie chart can be a function of the magnitude of the associated difference value.
  • the present invention provides a software system for configuring a computer system comprising a processor, and a memory device, to perform the steps of the methods of the present invention described above.
  • the present invention also provides a computer system including means for performing the steps of the method of the present invention.
  • the present invention provides a method, software system and computer system for automatically deforming a grid to locate individual DNA spot images and to quantify the spot images for measuring the local signal and background intensity values for the spot images, and to display such values.
  • the method and software system of the present invention automate data quantification and extraction in DNA array image analysis applications.
  • FIG. 1 is a graphic representation of a frame of image information including DNA images spots corresponding to an ideal micro-array of DNA spots;
  • FIG. 2 is a graphic representation of a frame of image information including DNA images spots corresponding to a typical micro-array of DNA spots;
  • FIG. 3 a illustrates the steps of an embodiment of a DNA array image analysis according the present invention
  • FIG. 3 b illustrates the steps of an embodiment of adjusting the position of a grid point in the method of FIG. 3 a;
  • FIG. 4 is a graphic representation of a grid with uniform spacing used to locate DNA spots in a typical example micro-array according to the method of FIG. 3 a;
  • FIG. 5 is a graphic representation of the grid of FIG. 4 deformed according to the method of FIG. 3 a to substantially match placement of DNA spot images in a micro-array;
  • FIG. 6 is a graphic illustration of updating the position of a grid point in a grid according to the method of FIG. 3 b;
  • FIG. 7 is a graphic illustration of segmenting a frame of image information including DNA images spots according to the method of FIG. 3 a;
  • FIG. 8 illustrates an example flow diagram for program instructions for implementing the DNA array image analysis method of FIG. 3 a
  • FIG. 9 illustrates the steps of an embodiment of a display method according the present invention for displaying quantified image information corresponding to DNA spots
  • FIG. 10 illustrates differential gene expression levels for different images of a micro-array displayed as pie charts according to an embodiment of a display method of FIG. 9;
  • FIGS. 11 a - d illustrates differential gene expression levels for different images of a micro-array displayed as bar graphs according to another embodiment of a display method of FIG. 9;
  • FIG. 12 illustrates an example flow diagram for program instructions for implementing the display method of FIG. 9.
  • FIG. 13 is an example block diagram of a computer system for implementing the present invention.
  • the present invention provides a method for automatically locating an array of DNA spot images 10 within a scanned image frame 12 of a DNA micro-array or a DNA macro-array, shown in FIGS. 1 and 2, wherein each spot corresponds to a particular gene or gene fragment.
  • the method of the present invention is applicable to both high-density micro-arrays, where spots are closely packed together on a solid surface, such as glass, with several thousands of spots placed in about 1 cm square area, and to macro-arrays with larger spacing of spots on surfaces such as membrane surfaces.
  • FIG. 2 is a graphic representation of the frame 12 of image information for a micro-array of DNA spots, including the DNA spot images 10 :
  • the image frame 12 is generated by scanning a micro-array with a particular laser frequency.
  • the spot images 10 are not in perfect alignment to each other, and there are large fluctuations in intensity, shape, and size of each spot in the micro-array.
  • the image information includes intensity levels for the spots corresponding to the level of expression of a particular gene.
  • an embodiment of the method of the present invention comprises the steps of: (a) storing the frame 12 of image information in a memory device 14 (step 16 ); (b) selecting a set 18 of image information within said frame 12 including a selected set of the DNA spot images 10 (step 20 ); (c) generating a grid 22 in said memory device 14 , the grid 22 including a plurality of spaced grid points 24 corresponding to said selected DNA spot images 10 , each grid point 24 including position information indicating the position of the grid point within said frame 12 (step 26 ); and (d) modifying a current position of at least one grid point 24 corresponding to a spot image 10 to shift the grid point 24 toward the corresponding spot image 10 (step 28 ). Step 28 can be repeated for said grid point 24 and for all the grid points 24 of the grid 22 .
  • the method of the present invention can be implemented as program instructions for configuring a computer system 34 , further described below, to perform the steps of the method of the present invention described herein.
  • the computer system 34 includes a processor 36 , the memory device 14 , an input device 38 and a display 40 Using the computer 34 , a user selects an image file containing the image frame 12 (control image) for processing, stores the image frame 12 in the memory 14 and displays it on the display 40 as the control image 12 .
  • the control image 12 includes a plurality of pixels each having an intensity level and a position within the control image 12 .
  • the user selects an image region 18 in the control image 12 by defining approximate four corners 42 of the image region 18 using the input device 38 . If not all corners 42 are visible, due to lack of DNA product at a particular location, the user can guess at a rough placement for a missing corner. Anchor spots can be used depending on the experiment to guarantee that all corners are visible.
  • the user then specifies the number of columns, C, and rows, R, of arrayed image spots 10 in the selected region 18 .
  • the computer 34 then automatically generates the grid 22 with equal spacing between each pair of corners having R rows and C columns within the specified region 18 .
  • the grid 22 comprises R ⁇ C grid points 24 , one grid point 24 for each intersection of a row with a column.
  • Each gird point 24 in the grid 22 except for those along the edges of the grid 22 , is displayed as connected to its four neighbors to the right, let, up, and down, through an elastic connection 46 . This placement of the grid points 24 establishes the starting configuration of the dynamic grid 22 as shown in FIG. 4.
  • the grid 22 can be represented in the memory 14 using two matrices: (i) a first matrix comprising an adjacency matrix of size R ⁇ C ⁇ 4 where each row number refers to a particular intersection point in the grid 22 and each column number refers to the neighboring intersection points arranged in a North, West, South, East fashion, and (ii) a second matrix comprising a position matrix of size R ⁇ C ⁇ 2 specifying the absolute location of each grid point 24 in the control image 12 .
  • the computer 34 since it is assumed that the pixel intensity corresponding to DNA spots images 10 in the image region 18 are greater than their surrounding background 50 intensity values, the computer 34 , according to the above steps, automatically shifts each grid point 24 towards local regions with the highest intensity values in subsequent iterations of said steps, wherein each grid point's location in the image frame 12 is modified.
  • FIG. 5 illustrates an example representations of the grid 22 with grid points 24 so shifted. A similar process can be applied to the image frame 12 in reverse video. Referring to FIG.
  • an embodiment of the step of modifying comprises: (i) selecting a first bounding area 52 in the control image 12 around the current position of a grid point 24 (step 54 ); (ii) generating a first position update including position information for updating a current position of said grid point 24 to a first new position 48 within the first bounding area 52 , the location of said first new position 48 relative to said current position being a function of intensity level of at least a portion of the image within the first bounding area 52 (step 56 ); (iii) generating a second position update including position information for updating said current position to a second new position 49 in the control image 12 , said second new position 49 being in a geometric arrangement with the position of one or more grid points 24 around said grid point 24 (step 58 ); and (iv) updating said current position with the position information of the first and the second position updates, thereby shifting said grid point 24 toward the corresponding spot image 10 (step 60 ).
  • FIG. 6 is a graphic illustration of updating the position of a grid point 24 in the grid 22 according to the above steps.
  • the first position update comprises a normalized direction vector d, based on the pixel intensity values within the first bounding area 52 around the grid point 24 .
  • the first bounding area 52 can comprise a bounding box of size r ⁇ r, or a circle of radius r centered on the current position of the grid point 24 .
  • the direction vector d can comprise an average or a weighted sum of vectors defining arrows originating at the center of the first bounding area 52 and ending at a plurality of the pixel locations within the first bounding area 52 .
  • the intensity value at each such pixel location can be used as the weight coefficient for calculating the weighted sum of said vectors.
  • the direction vector d can be based on the direction of the local intensity gradient. Other weighting coefficients can also be utilized.
  • the bounding box is represented as a matrix P in the memory 14 with n columns and m rows, and elements p ij corresponding to image intensity values at a location (i, j) in the bounding box.
  • the second position update is generated to place an additional constraint on the movement of the grid points 24 .
  • the constraint maintains the position of a grid point 24 in a linear geometric arrangement relative to position of one or more of neighboring grid points 24 in vertical and horizontal directions. Other geometric arrangements such as curves can also be selected.
  • the neighboring grid points 24 can be selected by the user, or automatically selected by the computer 34 , to include one or more of first and second order neighbors of the grid point 24 .
  • the second position update comprises another direction vector e, defining an arrow pointing at the mean location of the mid point between the first order neighbors in the horizontal direction to the left and right of the grid point 24 , and the mid point between the first order neighbors in the vertical direction to the top and bottom of the grid point 24 .
  • the direction vector, e attempts to keep the spacing between adjacent grid points 24 equal by using a linear geometric arrangement discussed above.
  • first and second order neighboring grid points around said grid point include spatial vectors: (1) A with elements A x and A y , (2) B with elements B x and B y , (3) C with elements C x and C y , (4) D with elements D x and D y , (5) E with elements E x and E y , (6) F with elements F x and F y , (7) G with elements G x and G y , and (8) H with elements H x and H y , as shown in diagram 1.
  • the computer 34 then linearly combines the direction vectors d and e to obtain a direction vector t for updating the position of the grid point 24 :
  • ⁇ and ⁇ are weighting coefficient parameters, with an example ⁇ or ⁇ range 0 to 10. The larger the value of ⁇ relative to ⁇ , the stiffer are connections 46 between adjacent grid points 24 .
  • is the update rate with an example range of 0 to 1.
  • the upper limit of said range for ⁇ or ⁇ is inversely proportion to an upper limit of ⁇ .
  • the local neighborhood size defined by the first bounding area 52 for each grid point 24 can be gradually reduced after each iteration, or every few iterations, of the modification step 28 described above.
  • the number of iterations is typically around forty and can be increased or decreased by the user to optimize the grid position appropriate for the image spots 10 .
  • the user can instruct the computer 34 to perform further tasks according to the present invention, including: (1) executing more iterations to optimize the location of the grid points 24 , (2) redrawing the grid 22 , (3) canceling out of the grid placement, or (4) accepting the current grid placement and proceed to segmentation and quantification steps 30 , 32 described below All of the above steps can be implemented using program instructions to be executed by a computer.
  • the method of the present invention further includes the step of segmenting the selected region 18 into a plurality of image segments 62 corresponding to the plurality of grid points 24 in the grid (step 32 ).
  • Each image segment 62 defines an area around a corresponding grid point 24 and includes a corresponding spot image 10 with minimum distance from said grid point 24 .
  • the size and shape of each segment 62 for each grid point 24 is a function of the spacing between said grid point 24 and one or more neighboring grid points 24 .
  • the programmed computer automatically segments the image region 18 into the segments 62 each having: (a) a width substantially equal to the smaller of: (i) the distance between the positions of the grid point 24 in the image segment 62 and the midpoint between said grid point 24 and an adjacent grid point to the left of said grid point 24 , and (ii) the distance between the positions of said grid point 24 and the midpoint between said grid point 24 and an adjacent grid point to the right of said grid point 24 ; and (b) a height substantially equal to the smaller of: (i) the distance between the positions of a said grid point 24 and the midpoint between said grid point 24 and an adjacent grid point to the left of said grid point 24 , and (ii) the distance between the positions of said grid point 24 and the midpoint between said grid point 24 an adjacent grid point to the right of said grid point 24 .
  • the method of the present invention can further include the step of quantifying at least a portion of image information in each image segment 62 to obtain image characteristic values for the image segment 62 (step 34 ).
  • Each spot image 1 0 in an image segment 62 can be used to measure the gene expression signal value and local background intensity levels according to a number of different user selected quantification methods. Five example quantification methods are described below.
  • Segmented Intensities This method includes sorting all the pixel intensities within an image segment 62 , selecting a portion, for example the top 10%, of said intensity values, and calculating the mean of the selected intensity values as a signal value. A similar portion, for example the bottom 10%, of intensity values is also selected and its mean value is provided as the local background intensity level.
  • Fixed Circle Total Intensity This method is similar to the Fixed Circle Mean Intensity method described above, except, total sum of all intensity values is provided in place of the mean values.
  • Automatic Circle Detection This method is similar to Fixed Circle Segmented Intensity method except, an automatic spot detection method is used to localize each image spot 10 in the image segment 62 .
  • a detection method can comprise a Hough transform for circle detection, a match-filter approach for optimum match between filters of various sizes to the data, or other detection methods.
  • the present invention provides an automatic method for refining the position of grid points 24 to optimally match the arrayed spot images 10 in micro-array images 12 , using the dynamic elastic grid 22 .
  • the user need only specify the four corners 42 of a region 18 in a micro-array image 12 , and the number of rows and columns of the image spots 10 in the micro-array.
  • Non-rectangular griddling patterns are also contemplated by the present invention.
  • the present invention provides a computer system 34 , described further below, for segmentation of the frame 12 of image information including the plurality of spaced DNA spot images 10 corresponding to the plurality of DNA spots, said image information including image intensity level information corresponding to said DNA spots.
  • the computer system comprises means for performing the above steps of the method of the present invention described herein and shown in FIGS. a and 3 b .
  • Said means include program instructions for configuring a general purpose, or dedicated computer, to perform said steps.
  • the present invention further provides a software system including program instructions for configuring a computer system to perform the above steps described herein and shown in FIGS. a and 3 b.
  • FIG. 8 illustrates an example general flow diagram for the program instructions of the computer system 34 and the software system of the present invention described above.
  • the program instructions include steps for: receiving and storing the image frame 12 in memory 14 (step 64 ); displaying the image frame 12 on the display 40 (step 66 ); obtaining four corners 42 of the image region 18 selected by user via the input device 38 (step 68 ); generating the grid 22 of R rows and C columns in the memory 14 as described above (step 70 ); forming a bounding box 52 of size r ⁇ r around a grid point (step 72 ); calculating the direction vector d as described above (step 74 ); calculating the direction vector e as described above (step 76 ); adjusting position of said grid point 24 with direction vectors d and e (step 78 ); displaying the adjusted grid 22 on the display 40 (step 80 ); determining if user is satisfied with the adjusted grid 22 (step ( 82 ); if not, proceeding to step 72 to adjust other grid point positions
  • the program instructions can be implemented utilizing a program language such as MatLabTM, C, Fortran, C++, and executed by a computer system 34 described below.
  • Mathematical calculations and image display can be implemented utilizing a simulation package or a math library such as MatLab, from MathworksTM, located in Natick.
  • the program instructions and related data are stored in the memory of the computer, to be executed by the processor to interact with the display, input device and storage in performing the steps described above.
  • the program instructions and related data can be used to program a dedicated graphics system to perform the above steps, the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner. In such a system, DNA micro-arrays are scanned into the memory device as image frames for processing as described above.
  • the values provided by one or more of the above quantification methods can be stored, as an ASCII file for example, and also saved in the memory 14 for comparison and display with similar quantified values corresponding to one or more other DNA micro-array images 12 according to another aspect of the present invention.
  • the user can also select one or more non-control images.
  • the grid position determined in steps described above can be applied, with any user defined translation or transformation, to the non-control image to quantify expression values according to the quantification methods described above.
  • the control image 12 is generated by scanning a micro-array with one particular laser frequency
  • the non-control image is generated by scanning the same micro-array with a different frequency laser.
  • each sensitive to one of the two laser frequencies are used to label DNA fragments from two different tissue types, e.g., healthy and diseased tissue. It is one of the main goals of micro-array data analysis to identify those sets of genes that are differentially expressed in different tissues. Extracted signal and background intensity levels for each gene (each DNA spot in the micro-array) can be displayed according to the present invention to visualize the differential gene expression levels between the control and non-control images.
  • an embodiment of the steps of such a display method for displaying image information corresponding to a plurality of DNA spot images 10 of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels comprises the steps of: (a) receiving image characteristic values for DNA spot images 10 in said control and non-control images (step 88 ); (b) for each DNA spot image 10 : (1) extracting said background and signal intensity levels from the image characteristic values for the spot image 10 , and (2) determining difference values between the background intensity levels and signal intensity levels (step 90 ); and (c) for each DNA spot: (1) relating the corresponding difference values to a range of graphic values, (2) selecting a graphic value for each difference value; and (3) displaying the selected graphic values (step 92 ).
  • the graphic values can include graphic objects 93 and color characteristic values as described below.
  • Applying said display method to the control and non-control images described above includes the steps of: (a) determining difference values between the background intensity levels and the signal intensity levels for both the control and non-control images, (b) displaying both difference values for all spots in the micro-array using a plurality of graphic objects 93 such as pie charts 94 or bar graphs 96 , each graphic object corresponding to a DNA spot, (c) probing each graphic object 93 to examine the expression levels, ratios, and other similar information, including displaying corresponding image segments form the control and non-control images for a selected graphic object 93 .
  • graphic objects 93 such as pie charts 94 or bar graphs 96
  • each difference value is associated to a segment 98 of a pie chart 94 having multiple segments, and the segments 98 as displayed a pie chart.
  • the area of each segment 98 of each pie chart 94 can be a function of the magnitude of the associated difference value.
  • color characteristic values can be assigned to the pie segments 98 by: (1) relating the corresponding difference values to a range of color characteristic values; (2) selecting a color characteristic value for each difference value; and (3) displaying the selected color characteristic value in the corresponding pie chart segments 98 .
  • each pie chart segment 98 among a plurality of pie segments can have a different color characteristic value.
  • the color characteristic values can include color, hue, brightness, intensity, and texture.
  • each pie chart 94 includes at least two segments: (1) a segment representing difference values between the background and signal intensity levels of a spot image 10 in the control image, corresponding to a DNA spot, and (2) another segment representing difference values between the background and signal intensity levels of a spot image 10 in the non-control image, corresponding to said DNA spot.
  • Each pie chart 94 can include additional segments for visualizing other differences associated with images 10 of a DNA spot in additional non-control images.
  • the graphic objects 93 can be arranged, for example, in the order in which their corresponding DNA spot images 10 appear in the control and non-control images. The user can also specify a different desired grouping of the graphic objects. The display arrangement can also be different from that of DNA spots.
  • the graphic objects 93 are shown as the bar graphs 96 , wherein each difference value is associated to a segment 98 of a bar graph having multiple segments, and the segments are displayed as a bar graph.
  • the arrangement, segment size and segment attributes of the bar graphs 96 can be identical to those of the pie charts 94 described above. Further, one or more type of the bar graphs 96 can be displayed in the same arrangement as shown for the pie charts 94 in FIG. 10.
  • the present invention provides a computer system 34 for displaying image information corresponding to the plurality of DNA spot images 10 of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels.
  • the computer system 34 comprises means for performing the steps of the display method described above. Said means include program instructions for configuring a general purpose or dedicated computer to perform said steps.
  • the present invention further provides a software system including program instructions for configuring a computer system to perform the steps of said display method.
  • FIG. 12 illustrates a general flow diagram for the program instructions of the display computer system and the display software system of the present invention described above.
  • the program instructions include steps for: receiving and storing said characteristic values in memory 14 (step 100 ); selecting a DNA spot and corresponding DNA spot images' characteristic values (ste 102 ); extracting background and signal intensity values for a image spot 10 (step 104 ); determining difference values between said background and intensity values (step 106 ); determining in step 108 if all spot images corresponding to said DNA spot have been so processed; if not, proceeding to step 104 to other spot images corresponding to said DNA spot, otherwise, relating the difference values to graphic values as described above (step 110 ); selecting a graphic value for each difference value as described above (step 112 ); displaying the selected graphic values as described above (step 114 ); and determining, in step 116 if all images for all DNA spots have been so processed, if not, proceeding to step 102 to process images for other DNA spots.
  • the program instructions can be implemented utilizing a program language such as MatLab, AVS/Expert, and Java, and executed by a computer system described below. Mathematical calculations can be implemented utilizing a simulation package or a math library such as MatLab, from Mathworks.
  • the program instructions and related data are stored in the memory of the computer, to be executed by the processor to interact with the display, input device and storage in performing the steps described above.
  • the program instructions and related data can be used to program a dedicated graphics system to perform the above steps, the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner. In such a system, DNA micro-arrays are scanned into the memory device as image frames for processing as described above.
  • the program instructions and related data can be used to program a dedicated graphics system to perform the above steps, the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner.
  • the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner.
  • image input means such as a scanner.
  • DNA micro-arrays are scanned into the memory device as image frames for processing as described above.
  • a suitable computer system 34 for executing said program instructions can be a dedicated computer such as a computer dedicated to scanning micro-arrays and processing micro-array images, or a general purpose computer system such as a personal computer or a dedicated computer system.
  • FIG. 13 shows a functional block diagram of the computer system 34 embodying the present.
  • a central processing unit (CPU) 36 operates on program instructions in the memory 14 using a processing unit 118 .
  • the CPU 36 also has a clock/calendar logic circuit 120 for maintaining an internal time/date clock.
  • a storage device 122 for storing information pertaining to micro-array images is connected to the CPU 36 over a bus 124 .
  • the micro-array images can be located on a file server 126 over a LAN or local to the CPU.
  • a keyboard 128 or mouse 38 receives instructions from the user concerning the DNA image micro-array analysis as necessary.
  • a scanner 130 allows scanning micro-arrays and obtaining images frames for processing as described above, and a printer 132 allows printing of images and data.
  • the main memory 14 stores the program instructions implementing the method of the present invention.
  • An example of a computer system suitable is a microcomputer equipped with a Pentium IITM microprocessor running at 266 MHZ. Such a system is preferably equipped with at least 64 MB megabytes of random access memory and a 2.0 GB hard drive.
  • the system preferably runs an operating system such as the WindowsTM operating environment. WindowsTM is manufactured by Microsoft Corporation, Redmond, Wash.

Abstract

A segmentation method of a frame of image information including a plurality of spaced DNA spot images corresponding to a plurality of DNA spots. The image information includes image intensity level information corresponding to said DNA spots. The frame is stored in a memory device and a set of image information within said frame including a selected set of the DNA spot images is selected. A grid including a plurality of spaced grid points corresponding to said selected DNA spot images is generated, such that each grid point includes position information indicating the position of the grid point within said frame. The current position of one or more grid points are adjusted by: selecting a first bounding area in the frame around the current position of the grid point; generating a first position update including position information for updating a current position of said grid point to a first new position within the first bounding area, the location of said first new position relative to said current position being a function of intensity level of at least a portion of the image within the first bounding area; generating a second position update including position information for updating said current position to a second new position in the frame, said second new position being in a geometric arrangement with the position of grid points around said grid point; and updating said current position with the position information of the first and the second position updates, thereby shifting said grid point toward the corresponding spot image. A display method diaplays image information corresponding to a plurality of DNA spot images of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels. For each DNA spot image: (1) background and signal intensity levels are extracted from the image characteristic values for the spot image, and (2) difference values between the background intensity levels and signal intensity levels are determined. For each DNA spot: (1) the corresponding difference values are related a range of graphic values, (2) a graphic value for each difference value is selected; and (3) the selected graphic values are displayed.

Description

    RELATED APPLICATION
  • This is a Continuation Application of application Ser. No. 09/020,155, filed Feb. 7, 1998, now pending.[0001]
  • FIELD OF THE INVENTION
  • The present invention relates to DNA array image analysis, and, in particular, to automatically segmenting DNA array images into individual DNA spot images for quantification. [0002]
  • BACKGROUND
  • Cellular behavior is primarily dictated by the selective expression of a subset of genes. Normal growth and differentiation depends on the appropriate genes being expressed in a desired context. Various disease states alter the normal expression of genes as compared to normal tissue. For example, malignant transformation of cancer tissues involves or induces altered gene expression. Through signal transduction cascades and transcriptional networks, alterations of one gene can impact a large number of genes and result in global effects on cell behavior. Regulation of translation and post-transcriptional modification play significant roles, but, invariably, signal transduction pathways lead to the nucleus and changes in gene transcription. [0003]
  • Therefore, there has been enormous interest in the development of techniques that allow the analysis of differential gene expression between different tissues or cell lines. One such technique includes use of ordered micro-arrays that allow two color fluorescence detection of hybridization signals. Individual DNA targets are arrayed on a small glass surface and hybridized with fluorescently labeled heterogeneous DNA probes derived from cDNA. The amount of fluorescence at each DNA spot correlates with the abundance of that DNA fragment in the probe mixture. [0004]
  • Using micro-arrays, gene expression levels can be quantitated at up to thousands of genes simultaneously. As hundreds of the same array can be printed, numerous tissues can be easily analyzed for relative expression levels. As such, the technique provides a powerful new tool for analyzing differential gene expression in numerous biologic problems. In addition to the determination of gene expression differences between tissues, genomic micro-arrays are useful for genomic mapping, genomic ploidy measurements and as hybridization targets for genomic mismatch scanning. Such techniques require rapid quantitative analysis of fluorescent hybridization for hundreds to tens of thousands of DNA spots. As such, there is a severe bottleneck in gene expression data collection due to inadequate methods for processing of individual DNA spot images for determining the quantitative fluorescent hybridization levels. [0005]
  • Some existing methods include manual processing of DNA spot images using a generic image processing tool, such as NIH image. Using such a tool a user visually locates each DNA spot image in a micro-array image, and moves a display pointer to each spot image, and manually defines a small area around the spot image. The image processing tool then reports image intensity values within the small area. The user then manually records the intensity values and continues this process for other visually located DNA spot images in the micro-array image. [0006]
  • However, such manual methods are impractical for micro-arrays with more than a handful of spot images. Further such methods are tedious and repetitive, requiring considerable time and effort. For example, with a micro-array image having about 600. DNA image spots, such manual methods can take about 8 hours of work, and resulting in quantification of only a limited number of image spots which visually seem to have a “good” expression level. As the micro-array density increases and becomes more complex, use of such methods becomes even more prohibitive. For example, current micro-array sizes range from several hundred to 1,200 genes, arrayed in a 1.8×1.8 cm area. As tip fabrication has improved, arrays with greater than 50,000 genes are viable. Such methods are also prone to various errors, including errors in manually recording the intensity values. Further such methods provide inconsistent quantification of intensity values, both for different spot images measured by a single individual, and for multiple individuals making measurements from the same micro-array image. [0007]
  • To alleviate the shortcomings of manual methods, some existing methods automate the process of locating DNA spot images from micro-array images and quantifying corresponding expression values. Such methods utilize a computer to manually position a cell grid on an area of the micro-array image containing an array of DNA spot images. The grid can be resized and individual columns and rows of the grid can be manually adjusted to better fit the arrayed pattern of DNA spot images. The grid position is then used by the computer to quantify the expression values using the intensity levels at each cell in the grid. However, such methods are inflexible since the grid placement requires extensive user interaction to fine-tune the grid. Further, the grid used in such methods is either completely fixed in shape, or has limited global flexibility (e.g., resizing and rotating the entire grid). [0008]
  • Such limitations cause a major handicap in most DNA array image analysis applications since DNA spots are never perfectly formed in a regular grid pattern in a micro-array such as shown in FIG. 1. Although a robot used in spotting DNA fragments on a glass surface has positional accuracy to within ±5 um, larger variations in the precise spacing of the arrayed DNA spots occur due to surface interactions of the solution with the silanized surface and tip variations. Moreover, printing tips are difficult to fabricate and many do not work uniformly. Therefore, as shown in FIG. 2, not only are DNA spots occasionally placed out of the regular grid pattern, but they also vary in size. It is therefore rare to have a fixed grid that can match exactly the pattern in the micro-array. Though in existing methods the grid can be manually resized, rotated, and a column or a row of the grid can be moved, the individual grid cells cannot be manipulated Therefore, such methods are impractical for most DNA array image analysis applications, and specially for high density micro-arrays [0009]
  • Further, DNA spot image signals derived from the micro-arrays are susceptible to surface noise and laser reflection, due to surface dust. And, nonspecific DNA binding to the silanized surface occurs in a non-uniform pattern creating a varying background of fluorescence over the surface. Existing methods are unable to cope with irregular micro-array pattern, search for DNA image spots, and accurately quantify specific signals while accounting for the local background. [0010]
  • Other existing methods do not use a grid at all but apply a “spot” filter to detect locations in the micro-array image which “look-like” DNA spot images. However, using such methods it is difficult to define what a spot should look like. Furthermore, extensive noise and variations in the spot shape, due to the processing and scanning mechanisms, significantly reduce the signal to noise ratio (SNR) of the spot images. Thus, the detection scheme misses many real spots and processes many false patches in the image as real DNA spot images. [0011]
  • Another disadvantage of existing systems is their inability to display micro-array image pixel intensities, corresponding to gene expression values in related DNA spots for example, in an intuitive manner. As such, the user cannot easily determine gene properties in such DNA spots. [0012]
  • There is, therefore, a need for a DNA array image analysis method for automatically segmenting DNA array images into individual DNA spot images for quantification. There is also a need for such method to process irregular micro-array patterns, search for DNA image spots, and accurately quantify, and intuitively display, specific signals while accounting for the local background. [0013]
  • SUMMARY
  • The present invention satisfies these needs. In one embodiment, the present invention provides a method for segmentation of a frame of image information including a plurality of spaced DNA spot images corresponding to a plurality of DNA spots, the image information including image intensity level and intra frame position information corresponding to said DNA spots. The method of the present invention comprises the steps of: (a) transferring the frame of image information into a memory device; (b) selecting a set of image information within said frame including a selected set of the DNA spot images; (c) generating a grid in said memory device, the grid including a plurality of spaced grid points corresponding to said selected DNA spot images, each grid point including position information indicating the position of the grid point within said frame; and (d) modifying a current position of at least one grid point corresponding to a spot image to shift the grid point toward the corresponding spot image. Step (d) can be repeated for said grid point and for all the grid points of the grid. [0014]
  • The step of modifying said current position includes: (i) selecting a first bounding area in the frame around the current position of the grid point; (ii) generating a first position update including position information for updating a current position of said grid point to a first new position within the first bounding area, the location of said first new position relative to said current position being a function of intensity level of at least a portion of the image within the first bounding area; (iii) generating a second position update including position information for updating said current position to a second new position in the frame, said second new position being in a geometric arrangement with the position of one or more grid points around said grid point; and (iv) updating said current position with the position information of the first and the second position updates, thereby shifting said grid point toward the corresponding spot image. The DNA spot images can be in a substantially two dimensional array arrangement, and generating the grid can include generating a two dimensional array of grid points spaced according to a predetermined criteria. [0015]
  • The method can further include the step of segmenting the selected set of image information by selecting at least one image segment defining a segment area around a grid point and including a spot image with minimum distance from said grid point, said segment area being a function of the spacing between said grid point and one or more neighboring grid points. The selected set of image information can further be segmented into a plurality of image segments corresponding to the plurality of grid points in the grid, each image segment defining a segment area around a corresponding grid point and including a corresponding spot image with minimum distance from said grid point, said segment area being a function of the spacing between said grid point and one or more neighboring grid points, wherein each spot image is contained in a corresponding image segment. [0016]
  • The method of the present invention can further include quantifying at least a portion of image information in said image segment to obtain image characteristic values for said image segment. The image characteristic values can include DNA information for a DNA spot corresponding to the DNA spot image in said image segment, said DNA information including gene expression values. [0017]
  • In another aspect, the present invention provides a method of displaying image information corresponding to a plurality of DNA spot images of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels. In one embodiment, the display method includes the steps of: (a) for each DNA spot image: (1) extracting said background and signal intensity levels from the image characteristic values for the spot image, and (2) determining difference values between the background intensity levels and signal intensity levels; and (b) for each DNA spot: (1) relating the corresponding difference values to a range of graphic values, (2) selecting a graphic value for each difference value, and (3) displaying the selected graphic values. The steps of relating and selecting can include associating each difference value to a segment of a pie chart having multiple segments, and the step of displaying the selected graphic values can include displaying said segments as a pie chart. The area of each segment of each pie chart can be a function of the magnitude of the associated difference value. [0018]
  • In another aspect, the present invention provides a software system for configuring a computer system comprising a processor, and a memory device, to perform the steps of the methods of the present invention described above. The present invention also provides a computer system including means for performing the steps of the method of the present invention. [0019]
  • As such, the present invention provides a method, software system and computer system for automatically deforming a grid to locate individual DNA spot images and to quantify the spot images for measuring the local signal and background intensity values for the spot images, and to display such values. The method and software system of the present invention automate data quantification and extraction in DNA array image analysis applications.[0020]
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • These and other features, aspects and advantages of the present invention will become better understood with regard to the following description, appended claims and accompanying drawings where: [0021]
  • FIG. 1 is a graphic representation of a frame of image information including DNA images spots corresponding to an ideal micro-array of DNA spots; [0022]
  • FIG. 2 is a graphic representation of a frame of image information including DNA images spots corresponding to a typical micro-array of DNA spots; [0023]
  • FIG. 3[0024] a illustrates the steps of an embodiment of a DNA array image analysis according the present invention;
  • FIG. 3[0025] b illustrates the steps of an embodiment of adjusting the position of a grid point in the method of FIG. 3a;
  • FIG. 4 is a graphic representation of a grid with uniform spacing used to locate DNA spots in a typical example micro-array according to the method of FIG. 3[0026] a;
  • FIG. 5 is a graphic representation of the grid of FIG. 4 deformed according to the method of FIG. 3[0027] a to substantially match placement of DNA spot images in a micro-array;
  • FIG. 6 is a graphic illustration of updating the position of a grid point in a grid according to the method of FIG. 3[0028] b;
  • FIG. 7 is a graphic illustration of segmenting a frame of image information including DNA images spots according to the method of FIG. 3[0029] a;
  • FIG. 8 illustrates an example flow diagram for program instructions for implementing the DNA array image analysis method of FIG. 3[0030] a;
  • FIG. 9 illustrates the steps of an embodiment of a display method according the present invention for displaying quantified image information corresponding to DNA spots; [0031]
  • FIG. 10 illustrates differential gene expression levels for different images of a micro-array displayed as pie charts according to an embodiment of a display method of FIG. 9; [0032]
  • FIGS. 11[0033] a-d illustrates differential gene expression levels for different images of a micro-array displayed as bar graphs according to another embodiment of a display method of FIG. 9;
  • FIG. 12 illustrates an example flow diagram for program instructions for implementing the display method of FIG. 9; and [0034]
  • FIG. 13 is an example block diagram of a computer system for implementing the present invention.[0035]
  • DESCRIPTION
  • In one embodiment, the present invention provides a method for automatically locating an array of [0036] DNA spot images 10 within a scanned image frame 12 of a DNA micro-array or a DNA macro-array, shown in FIGS. 1 and 2, wherein each spot corresponds to a particular gene or gene fragment. The method of the present invention is applicable to both high-density micro-arrays, where spots are closely packed together on a solid surface, such as glass, with several thousands of spots placed in about 1 cm square area, and to macro-arrays with larger spacing of spots on surfaces such as membrane surfaces.
  • FIG. 2 is a graphic representation of the [0037] frame 12 of image information for a micro-array of DNA spots, including the DNA spot images 10: Typically, the image frame 12 is generated by scanning a micro-array with a particular laser frequency. The spot images 10 are not in perfect alignment to each other, and there are large fluctuations in intensity, shape, and size of each spot in the micro-array. The image information includes intensity levels for the spots corresponding to the level of expression of a particular gene.
  • Referring to FIG. a, an embodiment of the method of the present invention comprises the steps of: (a) storing the [0038] frame 12 of image information in a memory device 14 (step 16); (b) selecting a set 18 of image information within said frame 12 including a selected set of the DNA spot images 10 (step 20); (c) generating a grid 22 in said memory device 14, the grid 22 including a plurality of spaced grid points 24 corresponding to said selected DNA spot images 10, each grid point 24 including position information indicating the position of the grid point within said frame 12 (step 26); and (d) modifying a current position of at least one grid point 24 corresponding to a spot image 10 to shift the grid point 24 toward the corresponding spot image 10 (step 28). Step 28 can be repeated for said grid point 24 and for all the grid points 24 of the grid 22.
  • The method of the present invention can be implemented as program instructions for configuring a [0039] computer system 34, further described below, to perform the steps of the method of the present invention described herein. The computer system 34 includes a processor 36, the memory device 14, an input device 38 and a display 40 Using the computer 34, a user selects an image file containing the image frame 12 (control image) for processing, stores the image frame 12 in the memory 14 and displays it on the display 40 as the control image 12. The control image 12 includes a plurality of pixels each having an intensity level and a position within the control image 12. The user then selects an image region 18 in the control image 12 by defining approximate four corners 42 of the image region 18 using the input device 38. If not all corners 42 are visible, due to lack of DNA product at a particular location, the user can guess at a rough placement for a missing corner. Anchor spots can be used depending on the experiment to guarantee that all corners are visible.
  • The user then specifies the number of columns, C, and rows, R, of arrayed image spots [0040] 10 in the selected region 18. The computer 34 then automatically generates the grid 22 with equal spacing between each pair of corners having R rows and C columns within the specified region 18. The grid 22 comprises R×C grid points 24, one grid point 24 for each intersection of a row with a column. Each gird point 24 in the grid 22, except for those along the edges of the grid 22, is displayed as connected to its four neighbors to the right, let, up, and down, through an elastic connection 46. This placement of the grid points 24 establishes the starting configuration of the dynamic grid 22 as shown in FIG. 4. The grid 22 can be represented in the memory 14 using two matrices: (i) a first matrix comprising an adjacency matrix of size R×C×4 where each row number refers to a particular intersection point in the grid 22 and each column number refers to the neighboring intersection points arranged in a North, West, South, East fashion, and (ii) a second matrix comprising a position matrix of size R×C ×2 specifying the absolute location of each grid point 24 in the control image 12.
  • Since it is assumed that the pixel intensity corresponding to [0041] DNA spots images 10 in the image region 18 are greater than their surrounding background 50 intensity values, the computer 34, according to the above steps, automatically shifts each grid point 24 towards local regions with the highest intensity values in subsequent iterations of said steps, wherein each grid point's location in the image frame 12 is modified. FIG. 5 illustrates an example representations of the grid 22 with grid points 24 so shifted. A similar process can be applied to the image frame 12 in reverse video. Referring to FIG. 3b, an embodiment of the step of modifying (step 28) comprises: (i) selecting a first bounding area 52 in the control image 12 around the current position of a grid point 24 (step 54); (ii) generating a first position update including position information for updating a current position of said grid point 24 to a first new position 48 within the first bounding area 52, the location of said first new position 48 relative to said current position being a function of intensity level of at least a portion of the image within the first bounding area 52 (step 56); (iii) generating a second position update including position information for updating said current position to a second new position 49 in the control image 12, said second new position 49 being in a geometric arrangement with the position of one or more grid points 24 around said grid point 24 (step 58); and (iv) updating said current position with the position information of the first and the second position updates, thereby shifting said grid point 24 toward the corresponding spot image 10 (step 60).
  • The position matrix elements are modified and updated by the [0042] computer 34 during multiple iterations of the above steps. FIG. 6 is a graphic illustration of updating the position of a grid point 24 in the grid 22 according to the above steps. In the embodiment of the invention described herein, the first position update comprises a normalized direction vector d, based on the pixel intensity values within the first bounding area 52 around the grid point 24. The first bounding area 52 can comprise a bounding box of size r×r, or a circle of radius r centered on the current position of the grid point 24. The direction vector d can comprise an average or a weighted sum of vectors defining arrows originating at the center of the first bounding area 52 and ending at a plurality of the pixel locations within the first bounding area 52. The intensity value at each such pixel location can be used as the weight coefficient for calculating the weighted sum of said vectors. The direction vector d can be based on the direction of the local intensity gradient. Other weighting coefficients can also be utilized.
  • An example calculation of the direction vector d for said bounding box of size r×r is described below. The bounding box is represented as a matrix P in the [0043] memory 14 with n columns and m rows, and elements pij corresponding to image intensity values at a location (i, j) in the bounding box. The direction vector d is calculated as: T = i = 1 n j = 1 m p ij s j = i = 1 n P ij T t i = j = 1 m P ij T
    Figure US20030129648A1-20030710-M00001
  • x[0044] L≡Number of pixels from the left edge to the center of P
  • x[0045] R≡Number of pixels from the right edge to the center of P
  • y[0046] t≡Number of pixels from the top edge to the center of P
  • y[0047] b≡Number of pixels from the bottom edge to the center of P X = [ - x L - ( x L + 1 ) - 1012 ( x R + 1 ) x R ] Y = [ - y b - ( y b + 1 ) - 1012 ( y t + 1 ) y t ] dx = j = 1 m s j X j dy = i = 1 n t i Y i
    Figure US20030129648A1-20030710-M00002
  • d=[dx dy]
  • Using a priori information about the location of [0048] DNA spot images 10 in the frame 12, i.e., almost a uniform 2-D array, the second position update is generated to place an additional constraint on the movement of the grid points 24. In the embodiment described herein, the constraint maintains the position of a grid point 24 in a linear geometric arrangement relative to position of one or more of neighboring grid points 24 in vertical and horizontal directions. Other geometric arrangements such as curves can also be selected. The neighboring grid points 24 can be selected by the user, or automatically selected by the computer 34, to include one or more of first and second order neighbors of the grid point 24. In this embodiment, the second position update comprises another direction vector e, defining an arrow pointing at the mean location of the mid point between the first order neighbors in the horizontal direction to the left and right of the grid point 24, and the mid point between the first order neighbors in the vertical direction to the top and bottom of the grid point 24. The direction vector, e, attempts to keep the spacing between adjacent grid points 24 equal by using a linear geometric arrangement discussed above.
  • An example calculation of the direction vector e for a [0049] grid point 24 with a spatial position vector L having elements Lx and Ly is described below. Eight first and second order neighboring grid points around said grid point include spatial vectors: (1) A with elements Ax and Ay, (2) B with elements Bx and By, (3) C with elements Cx and Cy, (4) D with elements Dx and Dy, (5) E with elements Ex and Ey, (6) F with elements Fx and Fy, (7) G with elements Gx and Gy, and (8) H with elements Hx and Hy, as shown in diagram 1.
    Figure US20030129648A1-20030710-P00001
  • The vector e is calculated as: [0050]
  • When said first order neighbors are used: [0051] X ~ = ( A x + B x ) / 2 + ( C x + D x ) / 2 Y ~ = ( A y + B y ) / 2 + ( C y + D y ) / 2 dx = X ~ / 2 - L x dy = Y ~ / 2 - L y
    Figure US20030129648A1-20030710-M00003
  • When first and second order neighbors are used: [0052] X ~ = [ ( A x + B x ) + ( C x + D x ) ] + ( F x + E x ) / 2 + ( G x + H x ) / 2 Y ~ = [ ( A y + B y ) + ( C y + D y ) ] + ( F y + E y ) / 2 + ( G y + H y ) / 2 dx = X ~ / 6 - L x dy = Y ~ / 6 - L y
    Figure US20030129648A1-20030710-M00004
  • e=[dx dy]
  • The [0053] computer 34 then linearly combines the direction vectors d and e to obtain a direction vector t for updating the position of the grid point 24:
  • t=αd+βe
  • Where α and β are weighting coefficient parameters, with an example α or β range 0 to 10. The larger the value of β relative to α, the stiffer are [0054] connections 46 between adjacent grid points 24.
  • The spatial position L of said [0055] grid point 24 is then updated as:
  • L←L+ηt
  • Where η is the update rate with an example range of 0 to 1. The upper limit of said range for α or β is inversely proportion to an upper limit of η. [0056]
  • The local neighborhood size defined by the [0057] first bounding area 52 for each grid point 24 can be gradually reduced after each iteration, or every few iterations, of the modification step 28 described above. The number of iterations is typically around forty and can be increased or decreased by the user to optimize the grid position appropriate for the image spots 10.
  • After a number of iterations, the user can instruct the [0058] computer 34 to perform further tasks according to the present invention, including: (1) executing more iterations to optimize the location of the grid points 24, (2) redrawing the grid 22, (3) canceling out of the grid placement, or (4) accepting the current grid placement and proceed to segmentation and quantification steps 30, 32 described below All of the above steps can be implemented using program instructions to be executed by a computer.
  • Referring to FIG. 7, once the user is satisfied with the grid position, the method of the present invention further includes the step of segmenting the selected [0059] region 18 into a plurality of image segments 62 corresponding to the plurality of grid points 24 in the grid (step 32). Each image segment 62 defines an area around a corresponding grid point 24 and includes a corresponding spot image 10 with minimum distance from said grid point 24. The size and shape of each segment 62 for each grid point 24 is a function of the spacing between said grid point 24 and one or more neighboring grid points 24.
  • As an example, for the two-[0060] dimensional grid 22, the programmed computer automatically segments the image region 18 into the segments 62 each having: (a) a width substantially equal to the smaller of: (i) the distance between the positions of the grid point 24 in the image segment 62 and the midpoint between said grid point 24 and an adjacent grid point to the left of said grid point 24, and (ii) the distance between the positions of said grid point 24 and the midpoint between said grid point 24 and an adjacent grid point to the right of said grid point 24; and (b) a height substantially equal to the smaller of: (i) the distance between the positions of a said grid point 24 and the midpoint between said grid point 24 and an adjacent grid point to the left of said grid point 24, and (ii) the distance between the positions of said grid point 24 and the midpoint between said grid point 24 an adjacent grid point to the right of said grid point 24.
  • The method of the present invention can further include the step of quantifying at least a portion of image information in each [0061] image segment 62 to obtain image characteristic values for the image segment 62 (step 34). Each spot image 1 0 in an image segment 62 can be used to measure the gene expression signal value and local background intensity levels according to a number of different user selected quantification methods. Five example quantification methods are described below.
  • Segmented Intensities: This method includes sorting all the pixel intensities within an [0062] image segment 62, selecting a portion, for example the top 10%, of said intensity values, and calculating the mean of the selected intensity values as a signal value. A similar portion, for example the bottom 10%, of intensity values is also selected and its mean value is provided as the local background intensity level.
  • Fixed Circle Mean Intensity: In this method, a fixed circle of user specified size is centered at each [0063] grid point 24 in the image segment 62 The mean intensity value of all pixels within the circle is provided as the signal value of the image spot 10 in the image segment 62 and the mean intensity value of the surrounding pixels are reported as the background intensity levels.
  • Fixed Circle Total Intensity: This method is similar to the Fixed Circle Mean Intensity method described above, except, total sum of all intensity values is provided in place of the mean values. [0064]
  • Fixed Circle Segmented Intensity: This method is a combination of the above three methods where the mean of certain predefined portion of the intensity values within the circle is provided as the signal value and the mean of another predefined portion of the intensity values outside the circle is provided as the background intensity level. [0065]
  • Automatic Circle Detection: This method is similar to Fixed Circle Segmented Intensity method except, an automatic spot detection method is used to localize each [0066] image spot 10 in the image segment 62. Such a detection method can comprise a Hough transform for circle detection, a match-filter approach for optimum match between filters of various sizes to the data, or other detection methods.
  • Other quantification methods can also be utilized and are contemplated by the present invention. As such, the present invention provides an automatic method for refining the position of grid points [0067] 24 to optimally match the arrayed spot images 10 in micro-array images 12, using the dynamic elastic grid 22. The user need only specify the four corners 42 of a region 18 in a micro-array image 12, and the number of rows and columns of the image spots 10 in the micro-array. Non-rectangular griddling patterns are also contemplated by the present invention.
  • In another aspect, the present invention provides a [0068] computer system 34, described further below, for segmentation of the frame 12 of image information including the plurality of spaced DNA spot images 10 corresponding to the plurality of DNA spots, said image information including image intensity level information corresponding to said DNA spots. In one embodiment, the computer system comprises means for performing the above steps of the method of the present invention described herein and shown in FIGS. a and 3 b. Said means include program instructions for configuring a general purpose, or dedicated computer, to perform said steps. The present invention further provides a software system including program instructions for configuring a computer system to perform the above steps described herein and shown in FIGS. a and 3 b.
  • FIG. 8 illustrates an example general flow diagram for the program instructions of the [0069] computer system 34 and the software system of the present invention described above. Referring to the flow diagram, the program instructions include steps for: receiving and storing the image frame 12 in memory 14 (step 64); displaying the image frame 12 on the display 40 (step 66); obtaining four corners 42 of the image region 18 selected by user via the input device 38 (step 68); generating the grid 22 of R rows and C columns in the memory 14 as described above (step 70); forming a bounding box 52 of size r×r around a grid point (step 72); calculating the direction vector d as described above (step 74); calculating the direction vector e as described above (step 76); adjusting position of said grid point 24 with direction vectors d and e (step 78); displaying the adjusted grid 22 on the display 40 (step 80); determining if user is satisfied with the adjusted grid 22 (step (82); if not, proceeding to step 72 to adjust other grid point positions, otherwise, proceeding to step 84 to segment the selected image 18 into segments 62 as described above; and quantifying DNA spot image information for spot images 10 in the segments 62 as described above (step. 86).
  • The program instructions can be implemented utilizing a program language such as MatLab™, C, Fortran, C++, and executed by a [0070] computer system 34 described below. Mathematical calculations and image display can be implemented utilizing a simulation package or a math library such as MatLab, from Mathworks™, located in Natick. The program instructions and related data are stored in the memory of the computer, to be executed by the processor to interact with the display, input device and storage in performing the steps described above. Alternatively, the program instructions and related data can be used to program a dedicated graphics system to perform the above steps, the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner. In such a system, DNA micro-arrays are scanned into the memory device as image frames for processing as described above.
  • The values provided by one or more of the above quantification methods can be stored, as an ASCII file for example, and also saved in the [0071] memory 14 for comparison and display with similar quantified values corresponding to one or more other DNA micro-array images 12 according to another aspect of the present invention. In addition to the control image 12, the user can also select one or more non-control images. The grid position determined in steps described above can be applied, with any user defined translation or transformation, to the non-control image to quantify expression values according to the quantification methods described above. Typically, the control image 12 is generated by scanning a micro-array with one particular laser frequency, and the non-control image is generated by scanning the same micro-array with a different frequency laser. Different tags, each sensitive to one of the two laser frequencies are used to label DNA fragments from two different tissue types, e.g., healthy and diseased tissue. It is one of the main goals of micro-array data analysis to identify those sets of genes that are differentially expressed in different tissues. Extracted signal and background intensity levels for each gene (each DNA spot in the micro-array) can be displayed according to the present invention to visualize the differential gene expression levels between the control and non-control images.
  • Referring to FIG. 9, an embodiment of the steps of such a display method for displaying image information corresponding to a plurality of [0072] DNA spot images 10 of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels, comprises the steps of: (a) receiving image characteristic values for DNA spot images 10 in said control and non-control images (step 88); (b) for each DNA spot image 10: (1) extracting said background and signal intensity levels from the image characteristic values for the spot image 10, and (2) determining difference values between the background intensity levels and signal intensity levels (step 90); and (c) for each DNA spot: (1) relating the corresponding difference values to a range of graphic values, (2) selecting a graphic value for each difference value; and (3) displaying the selected graphic values (step 92). The graphic values can include graphic objects 93 and color characteristic values as described below.
  • Applying said display method to the control and non-control images described above, includes the steps of: (a) determining difference values between the background intensity levels and the signal intensity levels for both the control and non-control images, (b) displaying both difference values for all spots in the micro-array using a plurality of [0073] graphic objects 93 such as pie charts 94 or bar graphs 96, each graphic object corresponding to a DNA spot, (c) probing each graphic object 93 to examine the expression levels, ratios, and other similar information, including displaying corresponding image segments form the control and non-control images for a selected graphic object 93.
  • Referring to FIG. 10, each difference value is associated to a [0074] segment 98 of a pie chart 94 having multiple segments, and the segments 98 as displayed a pie chart. The area of each segment 98 of each pie chart 94 can be a function of the magnitude of the associated difference value. Further, color characteristic values can be assigned to the pie segments 98 by: (1) relating the corresponding difference values to a range of color characteristic values; (2) selecting a color characteristic value for each difference value; and (3) displaying the selected color characteristic value in the corresponding pie chart segments 98. As such, each pie chart segment 98 among a plurality of pie segments can have a different color characteristic value. The color characteristic values can include color, hue, brightness, intensity, and texture.
  • In the example [0075] pie chart embodiment 94 of the graphic objects 93 shown in FIG. 10, each pie chart 94 includes at least two segments: (1) a segment representing difference values between the background and signal intensity levels of a spot image 10 in the control image, corresponding to a DNA spot, and (2) another segment representing difference values between the background and signal intensity levels of a spot image 10 in the non-control image, corresponding to said DNA spot. Each pie chart 94 can include additional segments for visualizing other differences associated with images 10 of a DNA spot in additional non-control images. As shown in FIG. 10 the graphic objects 93 can be arranged, for example, in the order in which their corresponding DNA spot images 10 appear in the control and non-control images. The user can also specify a different desired grouping of the graphic objects. The display arrangement can also be different from that of DNA spots.
  • Referring to FIGS. 11[0076] a-d, the graphic objects 93 are shown as the bar graphs 96, wherein each difference value is associated to a segment 98 of a bar graph having multiple segments, and the segments are displayed as a bar graph. The arrangement, segment size and segment attributes of the bar graphs 96 can be identical to those of the pie charts 94 described above. Further, one or more type of the bar graphs 96 can be displayed in the same arrangement as shown for the pie charts 94 in FIG. 10.
  • In another aspect, the present invention provides a [0077] computer system 34 for displaying image information corresponding to the plurality of DNA spot images 10 of at least one DNA spot, the image information including image characteristic values including background and signal intensity levels. In one embodiment, the computer system 34 comprises means for performing the steps of the display method described above. Said means include program instructions for configuring a general purpose or dedicated computer to perform said steps. The present invention further provides a software system including program instructions for configuring a computer system to perform the steps of said display method.
  • FIG. 12 illustrates a general flow diagram for the program instructions of the display computer system and the display software system of the present invention described above. Referring to the flow diagram, the program instructions include steps for: receiving and storing said characteristic values in memory [0078] 14 (step 100); selecting a DNA spot and corresponding DNA spot images' characteristic values (ste 102); extracting background and signal intensity values for a image spot 10 (step 104); determining difference values between said background and intensity values (step 106); determining in step 108 if all spot images corresponding to said DNA spot have been so processed; if not, proceeding to step 104 to other spot images corresponding to said DNA spot, otherwise, relating the difference values to graphic values as described above (step 110); selecting a graphic value for each difference value as described above (step 112); displaying the selected graphic values as described above (step 114); and determining, in step 116 if all images for all DNA spots have been so processed, if not, proceeding to step 102 to process images for other DNA spots.
  • The program instructions can be implemented utilizing a program language such as MatLab, AVS/Expert, and Java, and executed by a computer system described below. Mathematical calculations can be implemented utilizing a simulation package or a math library such as MatLab, from Mathworks. The program instructions and related data are stored in the memory of the computer, to be executed by the processor to interact with the display, input device and storage in performing the steps described above. Alternatively, the program instructions and related data can be used to program a dedicated graphics system to perform the above steps, the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner. In such a system, DNA micro-arrays are scanned into the memory device as image frames for processing as described above. Alternatively, the program instructions and related data can be used to program a dedicated graphics system to perform the above steps, the graphics system including a processor, a memory device, display, input device, storage and image input means such as a scanner. In such a system, DNA micro-arrays are scanned into the memory device as image frames for processing as described above. [0079]
  • A [0080] suitable computer system 34 for executing said program instructions can be a dedicated computer such as a computer dedicated to scanning micro-arrays and processing micro-array images, or a general purpose computer system such as a personal computer or a dedicated computer system. FIG. 13 shows a functional block diagram of the computer system 34 embodying the present. A central processing unit (CPU) 36 operates on program instructions in the memory 14 using a processing unit 118. The CPU 36 also has a clock/calendar logic circuit 120 for maintaining an internal time/date clock. A storage device 122 for storing information pertaining to micro-array images is connected to the CPU 36 over a bus 124. The micro-array images can be located on a file server 126 over a LAN or local to the CPU. A keyboard 128 or mouse 38 receives instructions from the user concerning the DNA image micro-array analysis as necessary. A scanner 130 allows scanning micro-arrays and obtaining images frames for processing as described above, and a printer 132 allows printing of images and data. The main memory 14 stores the program instructions implementing the method of the present invention. An example of a computer system suitable is a microcomputer equipped with a Pentium II™ microprocessor running at 266 MHZ. Such a system is preferably equipped with at least 64 MB megabytes of random access memory and a 2.0 GB hard drive. The system preferably runs an operating system such as the Windows™ operating environment. Windows™ is manufactured by Microsoft Corporation, Redmond, Wash.
  • Although the present invention has been described in considerable detail with regard to the preferred versions thereof, other versions are possible. Therefore, the appended claims should not be limited to the descriptions of the preferred versions contained herein. [0081]

Claims (10)

What is claimed is:
1. A method for segmentation of a frame of image information including a plurality of spaced DNA spot images corresponding to a plurality of DNA spots, said image information including image intensity level information corresponding to said DNA spots, the method comprising the steps of:
(a) storing the frame of image information in memory;
(b) selecting a set of image information within said frame including a selected set of the DNA spot images;
(c) generating a grid in memory, the grid including a plurality of spaced grid points corresponding to said selected set of DNA spot images, the grid points having a predefined relationship, each grid point including position information indicating the position of the grid point within said image frame;
(d) segmenting the selected set of image information by selecting at least one image segment defining a segment area around a grid point and including a spot image; and
(e) quantifying at least a portion of image information in said image segment to obtain image characteristic values for said image segment.
2. The method of claim 1, wherein said segment area is a function of the spacing between said grid point and one or more neighboring grid points.
3. The method of claim 1, wherein the image characteristic values include DNA information for a DNA spot corresponding to the DNA spot image in said image segment, said DNA information including gene expression values.
4. The method of claim 1, wherein the frame of image information includes a plurality of pixels each having an intensity level, and wherein the step of quantifying includes: (i) sorting at least a portion of the pixel intensities within said image segment, (ii) selecting a portion of said pixels, and (iii) computing an image characteristic value for the selected pixel values as function of the intensities of at least a portion of the selected pixel values.
5. The method of claim 1, wherein the frame of image information includes a plurality of pixels each having an intensity level, and wherein the step of quantifying includes: (i) selecting a subset of said pixels in said image segment, (ii) computing a first image characteristic value as a function of at least a portion of the intensities of the selected pixel values, and (iii) computing a second image characteristic value as a function of intensities of at least a portion of pixels proximate said subset of pixels.
6. A software system for configuring a computer system comprising a processor, and memory, for segmentation of a frame of image information including a plurality of spaced DNA spot images corresponding to a plurality of DNA spots, said image information including image intensity level and intra frame position information corresponding to said DNA spots, the software system comprising program instructions for:
(a) storing the frame of image information in memory;
(b) selecting a set of image information within said frame including a selected set of the DNA spot images;
(c) generating a grid in memory, the grid including a plurality of spaced grid points corresponding to said selected set of DNA spot images, the grid points having a predefined relationship, each grid point including position information indicating the position of the grid point within said image frame;
(d) segmenting the selected set of image information by selecting at least one image segment defining a segment area around a grid point and including a spot image; and
(e) quantifying at least a portion of image information in said image segment to obtain image characteristic values for said image segment.
7. The software system of claim 6, wherein said segment area is a function of the spacing between said grid point and one or more neighboring grid points.
8. The software system of claim 6, wherein said image characteristic values include DNA information for a DNA spot corresponding to the DNA spot image in said image segment, said DNA information including gene expression values.
9. The software system of claim 6, wherein the frame of image information includes a plurality of pixels each having an intensity level, and wherein the program instructions for quantifying include program instructions for: (i) sorting all the pixel intensities within said image segment, (ii) selecting a portion of said pixels, (iii) computing an image characteristic value for the selected pixel values as function of the intensities of at least a portion of the selected pixel value.
10. The software system of claim 6, wherein the frame image information includes a plurality of pixels each having an intensity level, and wherein the program instructions for quantifying include program instructions for: (i) selecting a subset of said pixels in said image segment, (ii) computing a first image characteristic value as a function of at least a portion of the intensities of the selected pixel values, and (iii) computing a second image characteristic value as a function of intensities of at least a portion of pixels proximate said subset of pixels.
US10/375,289 1998-02-07 2003-02-27 Automated DNA array image segmentation and analysis Pending US20030129648A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US10/375,289 US20030129648A1 (en) 1998-02-07 2003-02-27 Automated DNA array image segmentation and analysis

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US09/020,155 US6349144B1 (en) 1998-02-07 1998-02-07 Automated DNA array segmentation and analysis
US09/591,915 US6577956B1 (en) 1998-02-07 2000-06-09 Automated DNA array image segmentation and analysis
US10/375,289 US20030129648A1 (en) 1998-02-07 2003-02-27 Automated DNA array image segmentation and analysis

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US09/591,915 Continuation US6577956B1 (en) 1998-02-07 2000-06-09 Automated DNA array image segmentation and analysis

Publications (1)

Publication Number Publication Date
US20030129648A1 true US20030129648A1 (en) 2003-07-10

Family

ID=21797056

Family Applications (4)

Application Number Title Priority Date Filing Date
US09/020,155 Expired - Lifetime US6349144B1 (en) 1998-02-07 1998-02-07 Automated DNA array segmentation and analysis
US09/591,915 Expired - Lifetime US6577956B1 (en) 1998-02-07 2000-06-09 Automated DNA array image segmentation and analysis
US09/992,687 Expired - Lifetime US6674882B1 (en) 1998-02-07 2001-11-16 Automated DNA array image segmentation and analysis
US10/375,289 Pending US20030129648A1 (en) 1998-02-07 2003-02-27 Automated DNA array image segmentation and analysis

Family Applications Before (3)

Application Number Title Priority Date Filing Date
US09/020,155 Expired - Lifetime US6349144B1 (en) 1998-02-07 1998-02-07 Automated DNA array segmentation and analysis
US09/591,915 Expired - Lifetime US6577956B1 (en) 1998-02-07 2000-06-09 Automated DNA array image segmentation and analysis
US09/992,687 Expired - Lifetime US6674882B1 (en) 1998-02-07 2001-11-16 Automated DNA array image segmentation and analysis

Country Status (1)

Country Link
US (4) US6349144B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040001623A1 (en) * 2002-06-19 2004-01-01 Commissariat A L'energie Atomique Image analysis process for measuring the signal on biochips
US20100204057A1 (en) * 2009-02-10 2010-08-12 Samsung Electronics Co., Ltd. Substrate for microarray, method of manufacturing microarray using the same and method of obtaining light data from microarray

Families Citing this family (76)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6741344B1 (en) * 1994-02-10 2004-05-25 Affymetrix, Inc. Method and apparatus for detection of fluorescently labeled materials
US5631734A (en) * 1994-02-10 1997-05-20 Affymetrix, Inc. Method and apparatus for detection of fluorescently labeled materials
US6090555A (en) * 1997-12-11 2000-07-18 Affymetrix, Inc. Scanned image alignment systems and methods
DE69737883T2 (en) 1996-04-25 2008-03-06 Bioarray Solutions Ltd. LIGHT-REGULATED, ELECTROKINETIC COMPOSITION OF PARTICLES TO SURFACES
US6990221B2 (en) * 1998-02-07 2006-01-24 Biodiscovery, Inc. Automated DNA array image segmentation and analysis
EP1041514B1 (en) * 1999-03-30 2006-03-01 Fuji Photo Film Co., Ltd. Method and apparatus for selectively displaying measurement result and corresponding images
JP3469504B2 (en) * 1999-06-01 2003-11-25 日立ソフトウエアエンジニアリング株式会社 Microarray chip and indexing method thereof
US6754375B1 (en) * 1999-07-16 2004-06-22 Packard Bioscience Company Method and system for interactively developing at least one grid pattern and computer-readable storage medium having a program for executing the method
US6558623B1 (en) * 2000-07-06 2003-05-06 Robodesign International, Inc. Microarray dispensing with real-time verification and inspection
US7099502B2 (en) * 1999-10-12 2006-08-29 Biodiscovery, Inc. System and method for automatically processing microarrays
US6633669B1 (en) * 1999-10-21 2003-10-14 3M Innovative Properties Company Autogrid analysis
US7027629B2 (en) * 1999-11-05 2006-04-11 Agilent Technologies, Inc. Method of extracting locations of nucleic acid array features
US20050069956A1 (en) * 1999-11-22 2005-03-31 Michael Seul Color-encoding and in-situ interrogation of matrix-coupled chemical compounds
US6834122B2 (en) * 2000-01-22 2004-12-21 Kairos Scientific, Inc. Visualization and processing of multidimensional data using prefiltering and sorting criteria
US6591196B1 (en) * 2000-06-06 2003-07-08 Agilent Technologies Inc. Method and system for extracting data from surface array deposited features
DE60124598T2 (en) * 2000-06-06 2007-09-06 Agilent Technologies, Inc., Palo Alto Method and system for determining the orientation of a microarray
US9709559B2 (en) 2000-06-21 2017-07-18 Bioarray Solutions, Ltd. Multianalyte molecular analysis using application-specific random particle arrays
EP1311839B1 (en) * 2000-06-21 2006-03-01 Bioarray Solutions Ltd Multianalyte molecular analysis using application-specific random particle arrays
US7025933B2 (en) * 2000-07-06 2006-04-11 Robodesign International, Inc. Microarray dispensing with real-time verification and inspection
US7062092B2 (en) * 2000-08-22 2006-06-13 Affymetrix, Inc. System, method, and computer software product for gain adjustment in biological microarray scanner
US6965704B2 (en) * 2000-08-22 2005-11-15 Affymetrix, Inc. System, method, and computer software product for grid alignment of multiple scanned images
US6498863B1 (en) * 2000-09-20 2002-12-24 Media Cybernetics Inc. Method, system, and product for analyzing a digitized image of an array to create an image of a grid overlay
US20030045005A1 (en) * 2000-10-17 2003-03-06 Michael Seul Light-controlled electrokinetic assembly of particles near surfaces
US6829376B2 (en) 2000-10-24 2004-12-07 Affymetrix, Inc. Computer software system, method, and product for scanned image alignment
US7130458B2 (en) * 2000-10-24 2006-10-31 Affymetrix, Inc. Computer software system, method, and product for scanned image alignment
JP3913978B2 (en) * 2000-12-13 2007-05-09 富士フイルム株式会社 Image analysis method and apparatus
US6943036B2 (en) * 2001-04-30 2005-09-13 Agilent Technologies, Inc. Error detection in chemical array fabrication
US7262063B2 (en) 2001-06-21 2007-08-28 Bio Array Solutions, Ltd. Directed assembly of functional heterostructures
US6993172B2 (en) * 2001-06-29 2006-01-31 Agilent Technologies, Inc. Method and system for automated outlying feature and outlying feature background detection during processing of data scanned from a molecular array
US20030003600A1 (en) * 2001-07-02 2003-01-02 Fuji Photo Film Co., Ltd. Biochemical analysis kit and method for exposing stimulable phosphor sheet
ES2661167T3 (en) 2001-10-15 2018-03-27 Bioarray Solutions Ltd. Multiplexed analysis of polymorphic loci by simultaneous consultation and enzyme-mediated detection
US20030156136A1 (en) * 2002-02-15 2003-08-21 Cattell Herbert F. Method and system for visualization of results of feature extraction from molecular array data
US20040006431A1 (en) * 2002-03-21 2004-01-08 Affymetrix, Inc., A Corporation Organized Under The Laws Of Delaware System, method and computer software product for grid placement, alignment and analysis of images of biological probe arrays
US6946285B2 (en) * 2002-04-29 2005-09-20 Agilent Technologies, Inc. Arrays with elongated features
EP1502098B1 (en) * 2002-05-03 2008-09-17 ViaLogy, LLC Method for characterizing microarray output data
US20030220746A1 (en) * 2002-05-21 2003-11-27 Srinka Ghosh Method and system for computing and applying a global, multi-channel background correction to a feature-based data set obtained from scanning a molecular array
US7512496B2 (en) * 2002-09-25 2009-03-31 Soheil Shams Apparatus, method, and computer program product for determining confidence measures and combined confidence measures for assessing the quality of microarrays
AU2003298655A1 (en) 2002-11-15 2004-06-15 Bioarray Solutions, Ltd. Analysis, secure access to, and transmission of array images
US7206438B2 (en) * 2003-04-30 2007-04-17 Agilent Technologies, Inc. Feature locations in array reading
US7269518B2 (en) * 2003-04-30 2007-09-11 Agilent Technologies, Inc. Chemical array reading
US7206439B2 (en) * 2003-04-30 2007-04-17 Agilent Technologies, Inc. Feature locations in array reading
US20040243341A1 (en) * 2003-05-30 2004-12-02 Le Cocq Christian A. Feature extraction methods and systems
US20050038839A1 (en) * 2003-08-11 2005-02-17 Srinka Ghosh Method and system for evaluating a set of normalizing features and for iteratively refining a set of normalizing features
WO2005029705A2 (en) * 2003-09-18 2005-03-31 Bioarray Solutions, Ltd. Number coding for identification of subtypes of coded types of solid phase carriers
WO2005031305A2 (en) * 2003-09-22 2005-04-07 Bioarray Solutions, Ltd. Surface immobilized polyelectrolyte with multiple functional groups capable of covalently bonding to biomolecules
FR2860319B1 (en) * 2003-09-30 2006-02-10 Ct De Lutte Contre Le Cancer F METHOD FOR DETERMINING THE POSITION OF TISSUE SAMPLER CHIPS ON A BLADE OF TRANSPARENT MATERIAL
US7324677B2 (en) * 2003-10-14 2008-01-29 Agilent Technologies, Inc. Feature quantitation methods and system
EP1692298A4 (en) 2003-10-28 2008-08-13 Bioarray Solutions Ltd Optimization of gene expression analysis using immobilized capture probes
US20050089916A1 (en) * 2003-10-28 2005-04-28 Xiongwu Xia Allele assignment and probe selection in multiplexed assays of polymorphic targets
CN1882703B (en) 2003-10-29 2011-07-06 佰尔瑞溶液有限公司 Multiplexed nucleic acid analysis by fragmentation of double-stranded DNA
US7555154B2 (en) * 2004-03-10 2009-06-30 Biomedical Photometrics Inc. Fully automated segmentation of genetic micro-array images
US20050226535A1 (en) * 2004-04-13 2005-10-13 Srinka Ghosh Method and system for rectilinearizing an image of a microarray having a non-rectilinear feature arrangement
WO2005109316A2 (en) * 2004-05-03 2005-11-17 Perkinelmer Las, Inc. Method and apparatus for automatically segmenting a microarray image
US20050276512A1 (en) * 2004-06-01 2005-12-15 3M Innovative Properties Company Selective deconvolution of an image
US7302348B2 (en) 2004-06-02 2007-11-27 Agilent Technologies, Inc. Method and system for quantifying and removing spatial-intensity trends in microarray data
US20060040287A1 (en) * 2004-06-02 2006-02-23 Corson John F Method and system for quantifying random errors and any spatial-intensity trends present in microarray data sets
US20050281462A1 (en) * 2004-06-16 2005-12-22 Jayati Ghosh System and method of automated processing of multiple microarray images
US7363170B2 (en) * 2004-07-09 2008-04-22 Bio Array Solutions Ltd. Transfusion registry network providing real-time interaction between users and providers of genetically characterized blood products
US7848889B2 (en) 2004-08-02 2010-12-07 Bioarray Solutions, Ltd. Automated analysis of multiplexed probe-target interaction patterns: pattern matching and allele identification
US7818129B2 (en) * 2004-08-04 2010-10-19 Agilent Technologies, Inc. Detection of feature boundary pixels during array scanning
US7877213B2 (en) 2004-09-23 2011-01-25 Agilent Technologies, Inc. System and methods for automated processing of multiple chemical arrays
US8486629B2 (en) 2005-06-01 2013-07-16 Bioarray Solutions, Ltd. Creation of functionalized microparticle libraries by oligonucleotide ligation or elongation
US20070154923A1 (en) * 2005-12-29 2007-07-05 Affymetrix, Inc. Method for Gridding and Quality Control of Polymer Arrays
US7664336B2 (en) * 2006-01-26 2010-02-16 Microsoft Corporation Video noise reduction
US9445025B2 (en) 2006-01-27 2016-09-13 Affymetrix, Inc. System, method, and product for imaging probe arrays with small feature sizes
US8055098B2 (en) 2006-01-27 2011-11-08 Affymetrix, Inc. System, method, and product for imaging probe arrays with small feature sizes
US8009889B2 (en) 2006-06-27 2011-08-30 Affymetrix, Inc. Feature intensity reconstruction of biological probe array
JP5011394B2 (en) * 2006-12-01 2012-08-29 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Method for automatically decoding microarray images
WO2009112905A2 (en) * 2008-03-12 2009-09-17 Koninklijke Philips Electronics N.V. Real-time digital image processing architecture
CN101971208A (en) * 2008-03-12 2011-02-09 皇家飞利浦电子股份有限公司 Correction of spot area in measuring brightness of sample in biosensing device
AU2009201228B2 (en) * 2008-03-31 2011-09-22 Starkey Laboratories, Inc. Real ear measurement adaptor with internal sound conduit
EP2369325A1 (en) 2010-03-12 2011-09-28 Eppendorf Ag Array analysis for online detection
US8774494B2 (en) * 2010-04-30 2014-07-08 Complete Genomics, Inc. Method and system for accurate alignment and registration of array for DNA sequencing
JP5728878B2 (en) * 2010-10-13 2015-06-03 オムロン株式会社 Image processing apparatus and image processing method
CN103971589B (en) * 2013-01-28 2017-10-03 腾讯科技(深圳)有限公司 The processing method and processing device that the interest point information of map is made an addition in street view image
US9898820B2 (en) * 2013-03-18 2018-02-20 Life Technologies Corporation Methods and systems for analyzing biological reaction systems

Citations (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4550084A (en) * 1982-01-15 1985-10-29 Allied Corporation Analysis system
US4641528A (en) * 1985-09-16 1987-02-10 American Hospital Supply Corp. Specimen analysis instrument assembly
US5121320A (en) * 1988-10-17 1992-06-09 Hitachi Software Engineering Co., Ltd. System for reading and displaying an edit-processed DNA pattern
US5273632A (en) * 1992-11-19 1993-12-28 University Of Utah Research Foundation Methods and apparatus for analysis of chromatographic migration patterns
US5389792A (en) * 1993-01-04 1995-02-14 Grumman Aerospace Corporation Electron microprobe utilizing thermal detector arrays
US5417923A (en) * 1991-04-24 1995-05-23 Pfizer Inc. Assay tray assembly
US5541064A (en) * 1985-11-04 1996-07-30 Cell Analysis Systems, Inc. Methods and apparatus for immunoploidy analysis
US5560811A (en) * 1995-03-21 1996-10-01 Seurat Analytical Systems Incorporated Capillary electrophoresis apparatus and method
US5580728A (en) * 1994-06-17 1996-12-03 Perlin; Mark W. Method and system for genotyping
US5581631A (en) * 1994-09-20 1996-12-03 Neopath, Inc. Cytological system image collection integrity checking apparatus
US5583973A (en) * 1993-09-17 1996-12-10 Trustees Of Boston University Molecular modeling method and system
US5680514A (en) * 1994-09-23 1997-10-21 Hughes Electronics Multiple elastic feature net and method for target deghosting and tracking
US5695937A (en) * 1995-09-12 1997-12-09 The Johns Hopkins University School Of Medicine Method for serial analysis of gene expression
US5757954A (en) * 1994-09-20 1998-05-26 Neopath, Inc. Field prioritization apparatus and method
US5837475A (en) * 1997-01-30 1998-11-17 Hewlett-Packard Co. Apparatus and method for scanning a chemical array
US5851769A (en) * 1995-09-27 1998-12-22 The Regents Of The University Of California Quantitative DNA fiber mapping
US5853979A (en) * 1995-06-30 1998-12-29 Visible Genetics Inc. Method and system for DNA sequence determination and mutation detection with reference to a standard
US5865975A (en) * 1995-06-06 1999-02-02 Academy Of Applied Science Automatic protein and/or DNA analysis system and method
US5876933A (en) * 1994-09-29 1999-03-02 Perlin; Mark W. Method and system for genotyping
US5916747A (en) * 1995-06-30 1999-06-29 Visible Genetics Inc. Method and apparatus for alignment of signals for use in DNA based-calling
US5945284A (en) * 1997-05-27 1999-08-31 The Perkin-Elmer Corporation Length determination of nucleic acid repeat sequences by discontinuous primer extension
US5970164A (en) * 1994-08-11 1999-10-19 Sophisview Technologies, Ltd. System and method for diagnosis of living tissue diseases
US5981190A (en) * 1997-01-08 1999-11-09 Ontogeny, Inc. Analysis of gene expression, methods and reagents therefor
US5989835A (en) * 1997-02-27 1999-11-23 Cellomics, Inc. System for cell-based screening
US6127129A (en) * 1999-05-04 2000-10-03 Wisconsin Alumni Research Foundation Process to create biomolecule and/or cellular arrays on metal surfaces and product produced thereby
US6185561B1 (en) * 1998-09-17 2001-02-06 Affymetrix, Inc. Method and apparatus for providing and expression data mining database
US6226542B1 (en) * 1998-07-24 2001-05-01 Biosense, Inc. Three-dimensional reconstruction of intrabody organs
US6245517B1 (en) * 1998-09-29 2001-06-12 The United States Of America As Represented By The Department Of Health And Human Services Ratio-based decisions and the quantitative analysis of cDNA micro-array images
US6251601B1 (en) * 1999-02-02 2001-06-26 Vysis, Inc. Simultaneous measurement of gene expression and genomic abnormalities using nucleic acid microarrays
US6263092B1 (en) * 1996-07-10 2001-07-17 R2 Technology, Inc. Method and apparatus for fast detection of spiculated lesions in digital mammograms
US6301378B1 (en) * 1997-06-03 2001-10-09 R2 Technology, Inc. Method and apparatus for automated detection of masses in digital mammograms
US6303301B1 (en) * 1997-01-13 2001-10-16 Affymetrix, Inc. Expression monitoring for gene function identification
US6345115B1 (en) * 1997-08-07 2002-02-05 Imaging Research, Inc. Digital imaging system for assays in well plates, gels and blots
US6362832B1 (en) * 1999-09-01 2002-03-26 Packard Bioscience Company Method and system for overlaying at least three microarray images to obtain a multicolor composite image
US6362004B1 (en) * 1999-11-09 2002-03-26 Packard Biochip Technologies, Llc Apparatus and method for using fiducial marks on a microarray substrate
US6381058B2 (en) * 1996-08-16 2002-04-30 Imaging Research, Inc. Digital imaging system for assays in well plates, gels and blots
US20020052882A1 (en) * 2000-07-07 2002-05-02 Seth Taylor Method and apparatus for visualizing complex data sets
US6498863B1 (en) * 2000-09-20 2002-12-24 Media Cybernetics Inc. Method, system, and product for analyzing a digitized image of an array to create an image of a grid overlay
US6537749B2 (en) * 1998-04-03 2003-03-25 Phylos, Inc. Addressable protein arrays

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6103479A (en) * 1996-05-30 2000-08-15 Cellomics, Inc. Miniaturized cell array methods and apparatus for cell-based screening
EP1007737A4 (en) * 1997-07-25 2002-07-03 Affymetrix Inc System for providing a polymorphism database

Patent Citations (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4550084A (en) * 1982-01-15 1985-10-29 Allied Corporation Analysis system
US4641528A (en) * 1985-09-16 1987-02-10 American Hospital Supply Corp. Specimen analysis instrument assembly
US5541064A (en) * 1985-11-04 1996-07-30 Cell Analysis Systems, Inc. Methods and apparatus for immunoploidy analysis
US5121320A (en) * 1988-10-17 1992-06-09 Hitachi Software Engineering Co., Ltd. System for reading and displaying an edit-processed DNA pattern
US5417923A (en) * 1991-04-24 1995-05-23 Pfizer Inc. Assay tray assembly
US5273632A (en) * 1992-11-19 1993-12-28 University Of Utah Research Foundation Methods and apparatus for analysis of chromatographic migration patterns
US5389792A (en) * 1993-01-04 1995-02-14 Grumman Aerospace Corporation Electron microprobe utilizing thermal detector arrays
US5583973A (en) * 1993-09-17 1996-12-10 Trustees Of Boston University Molecular modeling method and system
US5580728A (en) * 1994-06-17 1996-12-03 Perlin; Mark W. Method and system for genotyping
US5970164A (en) * 1994-08-11 1999-10-19 Sophisview Technologies, Ltd. System and method for diagnosis of living tissue diseases
US5581631A (en) * 1994-09-20 1996-12-03 Neopath, Inc. Cytological system image collection integrity checking apparatus
US5757954A (en) * 1994-09-20 1998-05-26 Neopath, Inc. Field prioritization apparatus and method
US5680514A (en) * 1994-09-23 1997-10-21 Hughes Electronics Multiple elastic feature net and method for target deghosting and tracking
US5876933A (en) * 1994-09-29 1999-03-02 Perlin; Mark W. Method and system for genotyping
US5560811A (en) * 1995-03-21 1996-10-01 Seurat Analytical Systems Incorporated Capillary electrophoresis apparatus and method
US5865975A (en) * 1995-06-06 1999-02-02 Academy Of Applied Science Automatic protein and/or DNA analysis system and method
US5853979A (en) * 1995-06-30 1998-12-29 Visible Genetics Inc. Method and system for DNA sequence determination and mutation detection with reference to a standard
US5916747A (en) * 1995-06-30 1999-06-29 Visible Genetics Inc. Method and apparatus for alignment of signals for use in DNA based-calling
US5695937A (en) * 1995-09-12 1997-12-09 The Johns Hopkins University School Of Medicine Method for serial analysis of gene expression
US5851769A (en) * 1995-09-27 1998-12-22 The Regents Of The University Of California Quantitative DNA fiber mapping
US6263092B1 (en) * 1996-07-10 2001-07-17 R2 Technology, Inc. Method and apparatus for fast detection of spiculated lesions in digital mammograms
US6498690B2 (en) * 1996-08-16 2002-12-24 Imaging Research, Inc Digital imaging system for assays in well plates, gels and blots
US6441973B1 (en) * 1996-08-16 2002-08-27 Imaging Research, Inc. Digital imaging system for assays in well plates, gels and blots
US6381058B2 (en) * 1996-08-16 2002-04-30 Imaging Research, Inc. Digital imaging system for assays in well plates, gels and blots
US5981190A (en) * 1997-01-08 1999-11-09 Ontogeny, Inc. Analysis of gene expression, methods and reagents therefor
US6303301B1 (en) * 1997-01-13 2001-10-16 Affymetrix, Inc. Expression monitoring for gene function identification
US5837475A (en) * 1997-01-30 1998-11-17 Hewlett-Packard Co. Apparatus and method for scanning a chemical array
US5945679A (en) * 1997-01-30 1999-08-31 Hewlett-Packard Company Apparatus for scanning a chemical array
US5989835A (en) * 1997-02-27 1999-11-23 Cellomics, Inc. System for cell-based screening
US5945284A (en) * 1997-05-27 1999-08-31 The Perkin-Elmer Corporation Length determination of nucleic acid repeat sequences by discontinuous primer extension
US6301378B1 (en) * 1997-06-03 2001-10-09 R2 Technology, Inc. Method and apparatus for automated detection of masses in digital mammograms
US6345115B1 (en) * 1997-08-07 2002-02-05 Imaging Research, Inc. Digital imaging system for assays in well plates, gels and blots
US6537749B2 (en) * 1998-04-03 2003-03-25 Phylos, Inc. Addressable protein arrays
US6226542B1 (en) * 1998-07-24 2001-05-01 Biosense, Inc. Three-dimensional reconstruction of intrabody organs
US6185561B1 (en) * 1998-09-17 2001-02-06 Affymetrix, Inc. Method and apparatus for providing and expression data mining database
US6245517B1 (en) * 1998-09-29 2001-06-12 The United States Of America As Represented By The Department Of Health And Human Services Ratio-based decisions and the quantitative analysis of cDNA micro-array images
US6251601B1 (en) * 1999-02-02 2001-06-26 Vysis, Inc. Simultaneous measurement of gene expression and genomic abnormalities using nucleic acid microarrays
US6127129A (en) * 1999-05-04 2000-10-03 Wisconsin Alumni Research Foundation Process to create biomolecule and/or cellular arrays on metal surfaces and product produced thereby
US6362832B1 (en) * 1999-09-01 2002-03-26 Packard Bioscience Company Method and system for overlaying at least three microarray images to obtain a multicolor composite image
US6362004B1 (en) * 1999-11-09 2002-03-26 Packard Biochip Technologies, Llc Apparatus and method for using fiducial marks on a microarray substrate
US20020052882A1 (en) * 2000-07-07 2002-05-02 Seth Taylor Method and apparatus for visualizing complex data sets
US6498863B1 (en) * 2000-09-20 2002-12-24 Media Cybernetics Inc. Method, system, and product for analyzing a digitized image of an array to create an image of a grid overlay

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040001623A1 (en) * 2002-06-19 2004-01-01 Commissariat A L'energie Atomique Image analysis process for measuring the signal on biochips
US7136517B2 (en) * 2002-06-19 2006-11-14 Commissariat A L'energie Atomique Image analysis process for measuring the signal on biochips
US20100204057A1 (en) * 2009-02-10 2010-08-12 Samsung Electronics Co., Ltd. Substrate for microarray, method of manufacturing microarray using the same and method of obtaining light data from microarray

Also Published As

Publication number Publication date
US6674882B1 (en) 2004-01-06
US6349144B1 (en) 2002-02-19
US6577956B1 (en) 2003-06-10

Similar Documents

Publication Publication Date Title
US6674882B1 (en) Automated DNA array image segmentation and analysis
US6990221B2 (en) Automated DNA array image segmentation and analysis
US7317820B2 (en) System and method for automatically identifying sub-grids in a microarray
US6731781B1 (en) System and method for automatically processing microarrays
JP4420821B2 (en) Histological assessment
Seo et al. Interactively exploring hierarchical clustering results [gene identification]
US6633659B1 (en) System and method for automatically analyzing gene expression spots in a microarray
US6980677B2 (en) Method, system, and computer code for finding spots defined in biological microarrays
Balagurunathan et al. Simulation of cDNA microarrays via a parameterized random signal model
Bajcsy Gridline: automatic grid alignment DNA microarray scans
US7599090B2 (en) Method and apparatus for automatically segmenting a microarray image
US20060173628A1 (en) Method and system for determining feature-coordinate grid or subgrids of microarray images
US6754375B1 (en) Method and system for interactively developing at least one grid pattern and computer-readable storage medium having a program for executing the method
US7099502B2 (en) System and method for automatically processing microarrays
US6993173B2 (en) Methods for estimating probe cell locations in high-density synthetic DNA microarrays
US20020150935A1 (en) Dotted micro-array data extraction method
US20030198385A1 (en) Method apparatus for image analysis
US20030087289A1 (en) Image analysis of high-density synthetic DNA microarrays
Karthik et al. Automatic gridding of noisy microarray images based on coefficient of variation
US20040181342A1 (en) System and method for automatically analyzing gene expression spots in a microarray
CN114758104B (en) Automatic interference-dispersion-resistant bar arrangement marking method suitable for three-dimensional design steel bar drawing
US20080123898A1 (en) System and Method for Automatically Analyzing Gene Expression Spots in a Microarray
Barra Robust segmentation and analysis of DNA microarray spots using an adaptative split and merge algorithm
US7555154B2 (en) Fully automated segmentation of genetic micro-array images
Kuklin et al. High throughput screening of gene expression signatures

Legal Events

Date Code Title Description
AS Assignment

Owner name: BIODISCOVERY, INC., CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SHAMS, SOHEIL;REEL/FRAME:015704/0731

Effective date: 20000728

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED