EP1461769A1 - Method and arrangement for reducing noise in a microfluid device - Google Patents

Method and arrangement for reducing noise in a microfluid device

Info

Publication number
EP1461769A1
EP1461769A1 EP02793749A EP02793749A EP1461769A1 EP 1461769 A1 EP1461769 A1 EP 1461769A1 EP 02793749 A EP02793749 A EP 02793749A EP 02793749 A EP02793749 A EP 02793749A EP 1461769 A1 EP1461769 A1 EP 1461769A1
Authority
EP
European Patent Office
Prior art keywords
image
detection area
detection
area
substance
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.)
Withdrawn
Application number
EP02793749A
Other languages
German (de)
French (fr)
Inventor
Tobias SÖDERMAN
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.)
Gyros Patent AB
Original Assignee
Gyros AB
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
Priority claimed from SE0104461A external-priority patent/SE0104461D0/en
Priority claimed from US10/062,258 external-priority patent/US20030054563A1/en
Application filed by Gyros AB filed Critical Gyros AB
Publication of EP1461769A1 publication Critical patent/EP1461769A1/en
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N21/6452Individual samples arranged in a regular 2D-array, e.g. multiwell plates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6402Atomic fluorescence; Laser induced fluorescence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/75Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
    • G01N21/76Chemiluminescence; Bioluminescence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N35/00Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor
    • G01N35/00029Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor provided with flat sample substrates, e.g. slides
    • G01N35/00069Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor provided with flat sample substrates, e.g. slides whereby the sample substrate is of the bio-disk type, i.e. having the format of an optical disk
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01LCHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
    • B01L3/00Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
    • B01L3/50Containers for the purpose of retaining a material to be analysed, e.g. test tubes
    • B01L3/502Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
    • B01L3/5027Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1765Method using an image detector and processing of image signal

Definitions

  • the present invention generally concerns the field of microfluidic devices. More particularly, the present invention concerns a method and an arrangement for reducing noise in a substance raw data image containing noise. The present invention also relates to a method and an arrangement for determining a measure of at least one substance.
  • the substance raw data image has been obtained from a search area that comprises a detection area associated with a detection microcavity that is part of a microchannel structure of a microfluidic device detectable signals.
  • the substance raw data image reflects from a substance that is present in the detection microcavity after one or more liquid aliquots have been processed through the microchannel structure that comprises the detection microcavity. See further below.
  • the present invention belongs to the field of miniaturization of processes which comprise sample treatment, assay protocols, chemical and/or biochemical synthesis etc., within medicine, chemistry, biochemistry, molecular biology and the like.
  • one important goal within this field is to reduce the costs for these processes, for instance to reduce the amount of reagents needed per assay, reduce time per assay, etc.
  • One route has been to increase the degree of parallelity, for instance by integrating as many as possible of similar process runs in one and the same device in order to carry out all the runs in parallel.
  • large numbers of research groups and companies are involved in developing technology that will solve the numerous problems encountered.
  • One problem is related to the optimal way of configuring the detector in relation to the microdevice used for performing the processes while maintaining an acceptable sensitivity and reproducibihty. This problem may become particularly pronounced if the measuring step is performed by continuously moving the detector unit and the detection areas of a microfluidic device relative to each other during the measurement operation.
  • microfluidic device is in the form of a disc which is skewed because then it becomes difficult to maintain the optical focus in the right position relative to the detection areas. Without proper arrangement skewed discs will reduce sensitivity. This problem in particular applies to discs made of plastic material.
  • Another problem is related to maintaining an acceptable sensitivity and reproducibihty when changing sample volumes from the ⁇ l-range to the nl-range and performing the process protocols with a high degree of parallelity within the same device.
  • the inventors have found that under these circumstances the materials from which the microfluidic devices are fabricated and the various treatments during the manufacturing and conditioning of the devices easily introduce signal artifacts that are of the same kind and of comparable or larger size as the desired signals.
  • a more recent problem relates to the fact that the present assignee recently has managed to control the liquid flow in microfluidic devices containing a plurality of microchannel structures in such a way that the inter-channel variation for a device with respect to flow becomes insignificant.
  • This progress has enabled the assignee to quantify with a low inter-assay variation and a high sensitivity analytes, such as antigens, in the subfemtomole range in nl- volumes by carrying out the solid phase reaction of a heterogeneous sandwich immunoassay under flow conditions in small columns (nl-columns).
  • the present invention is directed to a method which meets the discussed problems.
  • a first object is to reduce the electronic, chemical, optical and dust noise that is present in the stored measurement data of a search area, i.e. that is present in the stored image of a detection area.
  • Electronic noise is introduced by the detector arrangement instrument.
  • the optical noise arises from the design of the disc, the lid of the disc (if present), the detector head (e.g. lenses) etc.
  • the chemical noise arise from chemicals and materials used in the manufacture of the disc and from undesired reactions including undesired adsorption related to the various reagents and solutions that have been used for placing the signal creating substance in the detection microcavity. Peak noise may derive from the manufacturing or the processing of the disc. Different kinds of dust, like particles, fibres, hair, etc may also introduce peak noise.
  • the invention typically reduces noise-from steps related to chemical noise, e.g. disturbing fluorescence in the case the substance raw data image is based on fluorescence.
  • a second object is to provide a method, software and arrangement for treating radiation data from one or more detection areas of the kind of microfluidic device described herein.
  • a third object is to provide arrangement, software and methods enabling accurate integration of radiation deriving from a desired substance as a function of subareas of individual detection areas. In particular this object aims at avoiding the problems discussed herein which the inventors have found may appear when measuring low amounts of substances with a high accuracy in microfluidic devices.
  • Figure 1 illustrates a schematic view on an arrangement of the invention and its main parts.
  • Figure 2 illustrates a detector head placed above a disc (microfluidic device) (cross- sectional view).
  • Figure 3 illustrates a set of microchannel structures that can be used in a circular disc.
  • Figure 4 illustrates the search and detection areas.
  • Figure 5 is a flow chart illustrating a first variant, ⁇ of a method step ⁇ according to the invention.
  • Figure 6 is a flowchart illustrating a second variant, ⁇ ", of the method step ⁇ according to the invention.
  • Figure 7 is a flowchart illustrating a second variant, ⁇ ", of the method step ⁇ according to the invention.
  • Figure 8 is a flowchart illustrating a method step ⁇ according to the invention.
  • Figure 9 illustrates an edge curve
  • Figure 10 is a flowchart illustrating a method step ⁇ according to the invention.
  • Figure 11 is a flowchart illustrating a method step ⁇ according to the invention.
  • Figure 12 is a flowchart illustrating a method step ⁇ according to the invention.
  • Figure 13 illustrates the search and detection areas and the surrounding areas A1-A3.
  • rotating refers to spinning. Yet further, the term “rotating” may also include, but is not limited to a step-wise rotation of a disc.
  • reagent includes, but is not limited to an analyte.
  • the term "in a circular manner” refers to around the center of the disc (circumferential direction).
  • a plurality of microchannel structures means two, three or more microchannel structures.
  • plural means > 10, such > 50 or > 100 microchannel structures.
  • removing noise and “reducing noise” means that the noise contribution in the measured data (values) is reduced to zero, or is negligible.
  • reduce and “remove” will be used alternately in the text, but there is no difference in the meaning of the two words in this specification. They should be regarded as equivalent.
  • a detector arrangement is typically adapted for measuring radiation from a plurality of detection areas (317a.b.c. etc., in figure 3) each of which is associated with a detection microcavity (306a,b,c, etc., in figure 3) in a microfluidic disc (101 in figure 1, 301 in figure 3).
  • the arrangement comprises:
  • a detector head (102) with a focal area and a disc holder (105) which are linked to a means enabling for the detector head (102), i.e., the focal area to transverse, the surface of the disc (101) in an essentially circular (means I) and/or radial manner (means II) when the disc is placed in the disc holder (105)
  • an angular aligning system (108,109) for recognizing the angular position of the part area which at a particular time is covered by the focal area
  • controller (112) e. , computer with software, which controls
  • the disc holder (105) is linked to at least one of the means for enabling the detector head and the disc to move relative to each other so that the focal area transverses the surface of the disc.
  • the disc holder (105) is linked to either means I or means II or both means I and means II.
  • means I, means II, and the detector head can be linked directly or indirectly to the controller.
  • Means I comprises three main variants with respect to how the disc and detector moves:
  • means I preferably comprises a spinner motor (103) and a shaft (104) carrying the disk holder (105). See figure 1.
  • the detector head (focal area) is circularly moving around the axis of symmetry of the disc.
  • means I preferably comprises a spinner motor and a shaft with an arm that carries the detector head. This shaft typically has the same direction as the axis of symmetry of the disc when placed on the disc holder.
  • Means II comprises three main variants with respect to how the disc and detector head moves:
  • means II preferably comprises a motor for laterally moving the disc holder.
  • means II preferably comprises a linear frame (106) that carries the detector head, and a drive unit (107) for radial movement of the focal area, over a disc which is placed in the disc holder (105). See Figure 1.
  • the detector head (102) and the motor (103) e.g., a spinner
  • the motor (103) controls the rotating speed that can be varied, e.g., within an interval between 0 - 15,000 rpm, such as within an interval between 60-5,000 rpm.
  • the rotation of the disc may be stepwise. When the disc is rotating, the focal area of the detector head will successively scan angularly adjacent part areas of the disc surface.
  • Means I in figure 1 is according to variant (1) above and comprises the motor (103) and the shaft (104).
  • the disc holder (105) is preferably a plate on which the disc can be placed.
  • the disc holder could also be a device that holds the disc at its periphery.
  • the side (114) of the plate facing the disc may comprise a system of evenly distributed uncovered shallow grooves or channel openings that are connected to a vacuum system by which the disc can be sucked to the plate. See for instance PCT/SE02/01682 (Gyros AB).
  • Detection areas may be present in one or both sides of the disc. If the disc holder is in the form of a plate as illustrated in figure 1 and the microfluidic disc has detection areas on the side opposing, it is important to secure the plate so as not to disturb the collection of radiation.
  • the plate may thus have a smaller diameter than the disc with the detector areas being located in an annular zone that is not covered by the plate.
  • the plate may be in a material that is translucent for the radiation utilized for the measurement, at least at the detector areas.
  • the angular aligning system may comprise:
  • a home position mark detector (108) which is able to detect when a home position mark (305) on a rotating disc (101) placed on the disc holder (105) is passing.
  • a home position mark (305 in figure 3) is preferably placed in an outer circumferential zone outside the detection areas or in some other position, which can be detected with high accuracy.
  • the position coordinates of each specific spot of the surface of the disc is given as the angular position relative to the home position mark and as the radial position relative to the circumference or axis of symmetry or relative to any other arbitrary fixed position on the disc.
  • a home position mark detector (108) typically has a fixed position outside the disc, for instance on the frame structure (113).
  • An accurate and preferred alternative for determining when a predetermined angular position is in front of the objective is to include an encoder that progressively gives the angular distance from the home position mark while the disc is rotating.
  • This kind of encoder ( 109) is typically associated with means I, for instance the motor (103), the shaft (104) or the disc holder (105). By associating the encoder directly with the disc (101) it is likely that the most accurate determination will be accomplished.
  • the encoder typically divides each revolution of the shaft into a large number of grades, for instance > 5 000, such as > 10 000 or > 20 000 or > 30 000.
  • a simple but less accurate alternative is to include calculating means that calculates the time needed from a preset rotation speed and the angular distance between the predetermined position and the home position mark (i.e., from the preset rotation speed and the angular position coordinate).
  • This kind of calculating means may be associated with the controller.
  • the angular aligning system should be able to give the angular position coordinate for the part of the disc which is in front of the objective with an accuracy of ⁇ 1°, such as within ⁇ 0.1° or within ⁇ 0.01° (provided there are 360° per revolution).
  • the exact accuracy needed will depend on the size of the disc, radial position of the detection area, the required sensitivity, size of detection area, etc.
  • the detector head (102) is guided on a linear frame (106) that may be the upper part of the frame structure (113) for linear displacement and positioning in a first plane PI, transversely through the central axis CL of shaft (104) and running in a radial direction thereto.
  • the linear frame (106) prohibits uncontrolled movement of the detector head in any other direction relative to this linear displacement.
  • the drive unit (107) for this displacement may be in the form of a translational responder for incrementally changing the position of the detector head (102) in the first plane PI (radial movement) and for enabling scanning of radially adjacent subareas of a microfluidic disc device placed in the disc holder (105).
  • Means II comprises the linear frame (106) and the drive unit (107) and is in figure 1 according to variant (2).
  • the drive unit (107) is associated with a unit (110) for determining the linear displacement and thus also the radial position coordinate of the part of the disc which is in front of the objective of the detector head (102).
  • This unit (110) may be in the form of an encoder that gives a translational position and movement of the focal area (objective of the detector head) in relation to a translational home position (111) on the responder. This home position in turn may be associated with a unique radial position in the disc.
  • the measuring unit (110) should be able to translate a translational position and movement of the detector head (focal area) into a radial position coordinate of the disc used with a high accuracy, typically within ⁇ 10 ⁇ m such as ⁇ 1 ⁇ m or ⁇ 0.1 ⁇ m.
  • the drive unit (107) and the vertical height of plane PI may be adjustable for focusing purposes.
  • Controller Control means for instance electronic and programmable control means
  • a computer program/computer program product i.e. a computer software application in a computer software application and computer program storage 120, and reading and accessing software applications from the same.
  • the present invention may be implemented as a computer program product, i.e. a computer software application, directly loadable into an internal memory storage of a processing unit within the controller 120.
  • Said computer program product comprises the software code means for performing the steps of the invented methods.
  • the present invention relates to a computer program product stored on a computer usable medium, comprising readable program for causing a processing unit in a computer means, such as the controller 112, to control an execution of the steps of the invented methods.
  • the computer usable medium is a record medium, a computer memory, a Read-Only
  • Different parts of the arrangement may communicate (1 15) with the controller (1 12).
  • the controller will in the preferred variants instruct the detector head to successively collect radiation from distinct and preselected parts of the surface of the disc.
  • the controller is programmed to start collecting radiation at a position, primarily an angular and/or a radial position, which is prior to an intended detection area, and to stop collecting at a position, which is after the same detection area.
  • the start and stop positions define a search area that at least should encompass the detection area or the relevant part of the detection area in which there is radiation from the substance. See further figure 4. If the radiation requires that the substance is irradiated, which is the case if-fluorescence is measured, the control means also defines the settings for the start and stop positions for irradiation. These latter settings are typically essentially the same as for collecting radiation.
  • the start and stop signals for collecting radiation is preferably directly linked to the angular positions in the disc at which collection is to start and end. respectively. This also includes that due account is taken for delays that may be inherent in the system or preset, i.e., the start and stop signals may have to be initiated before the focal area is positioned in front of the start and stop position, respectively.
  • the angular aligning system comprises an encoder, the encoder signals corresponding to a start position and a stop position are used to define the time
  • the start and stop for collecting radiation is linked to a preset rotating speed, i.e., the controller calculates from a preset rotating speed the time at which the start and stop position should be in front of the objective.
  • the controller may be programmed to change the radial position of the detector head (focal area) after a predetermined number of revolutions of the disc, for instance after 1. 2 or more revolutions with preference for 1.
  • the controller may be capable of changing the radial position of the objective (focal area) during a revolution of the disc.
  • the objective for these variants, one can envisage that radiation is collected for all relevant subareas at a common angular position before the disc is rotated (in a single step) to a subsequent angular position.
  • the objective for this variant the objective (focal area) is transversing the disc surface in a spiral-like manner, i.e , the radial position is changed successively during a revolution.
  • the collected radiation data is stored in a data storage 130 in a form that is retrievable for each individual subarea, for instance in the control unit.
  • the proper treatment of the data takes into account the overlapping effect and creates a true image of the radiation associated with different parts (subareas) of a detection area.
  • the microfluidic devices used in the present invention typically require detection of very low absolute amounts or concentrations of substances in the detection microcavity. It is therefore imperative in many variants of the invention that the detection principle shall enable detection and quantification of substance amounts that are ⁇ 10 " " mole per detection area/detection microcavity, such as ⁇ 10 " '” " ' mole or ⁇ 10 "1 mole per detection area/detection microcavity.
  • Typical detection principles may include, for example, spectrometric detection, such as radioactivity, absorption of radiation (e.g. light), fluorescence, chemiluminescence, bioluminescence, scattered light etc.
  • a change in wave-length(s), polarization, life-time, scattering, intensity, etc., between irradiation and radiation may be measured as a function of the presence of the substance of interest in the detection microcavity.
  • Typical examples include, but are not limited to luminescence and fluorescence principles, with laser induced fluorescence (LIF) being preferred.
  • the size of the focal area of the detector head typically is less than the size of a search and/or a detection area.
  • the width of a focal area in one direction is typically essentially ⁇ 1/5, such as of the ⁇ 1/10, of the corresponding width of a search and/or detection area and in a perpendicular direction (direction 2) within the same limits or larger.
  • a suitable detector head in the form of a pick-up head (200) is illustrated in figure 2 and comprises a laser source (201) whose beam is reflected on a dichroic mirror (202) and focused through an objective (205) to a part of a detection microcavity (203) positioned in front of the head.
  • the epi- fluorescent light is passed through the dichroic mirror and through a band-pass filter (206), selective for the flourochrome at hand and is finally focused onto the entrance of a photo multiplicator tube (PMT) (207) by means of an aspheric lens (208).
  • Pin-holes (209 and 210) are positioned between the entrance of the laser beam and the dichroic mirror, and in front of the PMT (207).
  • the size and position of the pin-hole (210) are adapted so that the focal area of the laser beam is inside the detection microcavity.
  • the size of pin-hole (209) is adapted so that preferentially emission light emanating from the focal area is passed into the PMT (207).
  • a mirror spot on a glass disc or the like may replace the dichroic mirror.
  • the microfluidic device used in the various aspects of the invention comprises a plurality of microchannel structures in which aliquots of liquids are transported and/or processed.
  • the devices typically are disc-shaped.
  • the structures are covered in the sense that their interior is in contact with ambient atmosphere primarily via separate inlet and/or outlet openings and/or vents.
  • Each microchannel structure comprises one or more detection microcavities and possibly also one or more reaction microcavities, and microconduits connecting these parts with each other.
  • a reaction microcavity may coincide with a detection microcavity.
  • the result of the processing in a microchannel structure is measured as radiation from a detection area which is directly or indirectly associated with a detection microcavity.
  • a disc typically has an axis of symmetry (C ⁇ ) perpendicular to the disc plane where n is an integer > 2, 3, 4, 5 5 preferably ⁇ (C ⁇ ).
  • the disc may thus be rectangular or have a circular shape.
  • inertia force may be used, for instance by spinning the disc (centrifugal force).
  • Other forces are capillary forces, electrokinetic forces. hydrodynamic forces etc.
  • Microfluidic devices that have an axis of symmetry and are intended for rotation may have a home position mark as discussed above.
  • microchannel contemplate that a channel structure comprises one or more cavities and/or channels/conduits that have a cross-sectional dimension that is ⁇ 10 ⁇ m, preferably ⁇ 10 " ⁇ m.
  • the lower limit is typically significantly larger than the size of the largest reagents and constituents of aliquots that are to pass through a microchannel.
  • the volumes of microcavities/microchambers are typically ⁇ 1000 nl, such as ⁇ 500 nl or ⁇ 100 nl or ⁇ 50 nl or ⁇ 25 nl, which in particular applies to the detection microcavities. Chambers/cavities directly connected to inlet ports for liquids may be considerably larger, e.g..
  • Microformat means that one, two, three or more liquid aliquots that are transported within the device have a volume in the ⁇ l-range, i.e., ⁇ 1000 ⁇ l such as ⁇ 100 ⁇ l or ⁇ 50 ⁇ l including but not limited to the nl -range (nanoformaf). such as ⁇ 1000 nl or ⁇ 500 nl or ⁇ 100 nl or ⁇ 50 nl.
  • the microfluidic device may be made from different materials, such as plastics, glass, silicone polymers, etc.
  • the detection area should be transparent/translucent for the detection principle utilized by the detector. From the manufacturing point of view plastic material is many times preferred because the costs for the material are normally low and mass production can easily be done, for instance by replication. Typical manufacturing processes involving plastic material are replication by embossing, moulding, etc., followed by attaching a top lid covering the open microchannel structures so obtained. See for instance WO 91 16966
  • the preferred disc material is plastic material, such as polycarbonates and plastic material based on monomers which consist of a polymerisable carbon-carbon double or triple bonds and saturated branched straight or cyclic alkyl and/or alkylene groups. Typical examples are Zeonex and Zeonor I M from Nippon Zeon, Japan, with preference for the latter.
  • silicone polymers such as poly dimethyl siloxane (PDMS) and the like are not considered to be plastic material.
  • black plastic material for instance containing graphite powder or carbon black, absorbs light and therefore has a low auto-fluorescence. Transport of light within black plastic material is prevented. Black plastic material will be very efficient for microfluidic devices when fluorescence and luminiscence measurements are relied upon. Black plastic material should in this case be avoided in the detection areas.
  • plastic lids of different or • iigin may be used, for instance MelinexTM 12 PET and MelinexTM 17 OPP (Du Pont, U.S.A.), etc
  • the plurality of dejection microcavities and the corresponding search and/or detection areas are preferably arranged in subgroups such that all members of a subgroup are positioned at the same radial distance and/or at the same angular position and/or have equal length and/or cross-sectional dimensions.
  • a detection area in the inventive arrangement typically has a size within the range of lxlO 2 - 2x10° ⁇ m 2 , such as lxlO - lO 3 ⁇ m 2 .
  • Their length and/or breadth are typically within the range of 0.5x10 - 5xl0 4 ⁇ m, such as 1x10 -10 4 ⁇ m.
  • Design microchannel structures (301a,b.c. etc.) that are adapted to a circular disc are illustrated in figure 3.
  • the structures are linked together by a common distribution channel (302) and a common waste channel (303).
  • the orientation of the microchannel structures around a common axis of symmetry is apparent.
  • the circumference (304) of the disc has a home position mark (305).
  • Each of the combined reaction/detection microcavities (306a,b,c, etc.) is communicating in the downstream direction with the common waste channel (303) and in the upstream direction via separate connections with the common distribution channel (302) and separate volume measuring units (307a,b,c, etc.).
  • a surface detection area (317a,b,c, etc.) is associated with each detection microcavity.
  • the common distribution channel (302) carries at one of its ends and at an intermediary position inlet ports (308 and 309. respectively). Another kind of inlet port (310) is located at each volume measuring unit (307a,b,c, etc.). Each microchannel structure (301a,b,c, etc.) also has an outlet to the common waste channel (303) and an outlet port (318) at the remaining end of the common distribution channel (302).
  • An inlet vent (31 1) to ambient atmosphere is connected to the common distribution channel via a common venting channel (3E2).
  • Other vents (313 and 314) to ambient atmosphere are placed in the common waste channel (303) and in the connecting microconduit between each reaction/detection microcavity (306a,b,c. etc.) and the common waste channel (303).
  • Appropriate valvings are positioned at 315 and 316 in each microchannel structure (301a,b,c, etc.).
  • the diameter of the preferred circular discs is about the same as conventional CDs but may be larger or smaller, for instance up to 300 % or larger and down to 10 % or less.
  • the processes that are carried out within the individual microchannel structures comprise assay protocols, organo-themical or biochemical synthesis protocols, etc.
  • the protocols comprise introduction of one or more liquid aliquots containing the necessary reagents/reactants into a microchannel structure.
  • one of the aliquots is a sample which is uncharacterized with respect to at least one feature, e.g.. type, form and/or amount of an analyte.
  • the processes comprises that the substance, which is associated with the radiation to be collected, is formed and/or retained in the detection microcavity under static conditions or under flow conditions.
  • the reaction system for retaining may be homogeneous or heterogeneous, i.e., with or without the desired substance being partitioned between a liquid phase and a solid phase.
  • the flow direction in the detection microcavity may be towards the circumference (outwards) or towards the center of the disc (inwards), or essentially parallel with the circumference of the disc. Also other directions may be utilized.
  • Process protocols may utilize specific reactions between reactants having mutual affinity to each other leading to a (a) formation of an affinity complex that is immobilized to a solid phase in a detection microcavity or (b) one or more other reaction products that may be soluble or insoluble in the detection microcavity.
  • the reaction conditions including selection of reactants it can often be arranged so that the product obtained and/or a reagent in excess are detectable with a signal that can be (a) measured from the above-mentioned detection areas and (b) related to one or more features of a starting liquid aliquot introduced into a microchannel structure. Typical such features are kind, form and amount including activity. etc., of a particular reactant (e.g.
  • an analyte including for instance an affinity reactant such as an enzyme etc.
  • the term "can be related to one or more features" includes also the determination of the manner in which reaction variables such as pH. ionic strength, detergents, etc., might influence the reaction used for forming the reaction product.
  • reaction variables such as pH. ionic strength, detergents, etc.
  • detection principles based on radioactivity, fluorescence, chemiluminescence, bioluminescence, enzymatic activity, chromogens, light scattering (turbidometry), etc., for instance by utilizing a reactant that carries a group providing detectability, either by being detectable as such or by being transformable to a detectable group.
  • the utilized process protocol means that a detectable reactant is incorporated into a complex or into some other reaction product.
  • Detectable products, reagents, etc , that can be formed and/or retained and measured iq a microcavity via a detection area are collectively called "substance" in other parts of this specification.
  • Typical reactants in this context include members of affinity pairs such as (a) antigen/hap ten and the corresponding antibody including its antibody active fragments, (b) lectin and the corresponding carbohydrate structure, (c) native ligands and the corresponding receptors.
  • complementary nucleic acids including synthetic variants such as synthetic oligonucleotides.
  • synthetic variants such as synthetic oligonucleotides.
  • Ig(Fc)-binding proteins and Protein A, Protein G and other Ig(Fc)-receptors ion pairs of opposite charges, enzyme and the substrate, inhibitor, cofactor. coenzyme etc., that can bind to the enzyme, etc.
  • Synthetic variants more or less mimicking a native affinity interaction are also included.
  • One aspect of the invention is a method and corresponding arrangement, for determining the amount of a substance in a detection microcavity of a microfluidic device and comprises collecting radiation associated with the substance from a detection area associated with the detection microcavity.
  • the method is characterized in comprising the steps of:
  • a microfluidic device e.g., in the form of a disc, comprising
  • each of the detection areas being (1 ) associated with one of said detection microcavities, (2) present in the surface of said device and (3) translucent/transparent for said radiation, and
  • a detector arrangement which is capable of collecting radiation from individual subareas of each of said search and detection areas;
  • step (b) scanning the search and/or detection areas associated with the detection microcavities that are part of microchannel structures in which step (b) has been carried out to obtain radiation from individual subareas of each scanned search or detection area, said scanning being performed by the use of said detector arrangement;
  • step (d) characterizing for each of the amounts obtained in step (d) a reaction variable that has been included in the process protocol used for each microchannel structure.
  • the characterization includes e.g., that the amount of the substance in absolute or relative terms in each of the detection microcavities is determined from each of the amounts obtained in step d).
  • the substance in this context may be the substance from which the radiation is directly derived or a substance, such as an analyte, the amount of which in one of the liquid aliquots used in step (b) is related to the amount of radiation determined in step (d). Since it is the result of the reaction that is measured the analyte may or may not be present in a detection microcavity after step (b).
  • Step (c) comprises collection of radiation intensities which refer to changes in the radiation due to the amount of a particular substance in the detection microcavity. This includes true intensity values but also other changes, such as changes in refraction, in scattering, in the plane of polarization and life time of fluorescent light (including half life, e.g.), etc.
  • step (c) also comprises non-circular scanning for instance scanning by lateral movement of a linear detector head comprising one or more rows of detector elements over the detection areas. Also imaging by a CCD camera is included.
  • the microfluidic disc may be circular or have some other geometric form, for instance triangular, rectangular, etc., including also irregular forms.
  • Step (d) means integration over each detection area, i.e., primarily over subareas which have radiation values that deviate from the values obtained for surface parts of the device that surrounds the detection area. Alternatively one may exclude all or selected parts of a detection area which corresponds to parts of a detector microcavity in which the presence of the substance is insignificant.
  • the integrating step includes the substeps of (a) finding a start and/or stop position at edges of the detection area, e.g., at an end corresponding to the inlet end of the detection microcavity, and (b the factual integrating.
  • substep (a) is carried out by determining the inflection point for the amount of radiation per subarea versus position along the detection area.
  • substep (a) comprises determining a threshold that segments detection area pixels from the background, e.g., if the detection microcavity is filled with a particle bed a threshold that segments the particle bed pixels from background pixels. This can be done with optimal thresholding or determining median or mean background or any other way of determining background (for "optimal thresholding” see for instance reference [1]: “Digital Image Processing", 2nd edition. Editors: Gonzales R C et al, page 354).
  • substep (b) will mean integrating a selection of pixels, i.e., those pixels which have radiation values above the threshold and belong to the main group (detection area) and excluding noise pixels above threshold that do not belong to the main group of connected pixels.
  • the integrating typically starts from pixels corresponding to one end of the detection area, e.g., the inlet end.
  • substeps (a) and (b) there are preferably additional substeps for refining the method, such as
  • the factual integrating (main step d) will in this case mean integrating the radiation values for the binary high pixels of the main group.
  • step (d) may also comprise substeps for removing peak noise.
  • One way of doing this is la Place filtering, point detection and the like. See for instance see reference [1]: “Digital Image Processing", 2nd edition, Editors: Gonzales R C et al, pages 333, 334 and 339 (1987).
  • the width and the position of each peak is calculated for instance by including edge detection or edge linking based on a local area etc.. See for instance reference [1]: "Digital Image Processing”. 2nd edition, Editors: Gonzales R
  • Step (e) is conventional and typically includes that the integrated value is compared with the value for one or more standards.
  • a standard value is typically the integrated value for a known amount of a standard substance, which in most cases is the same as the substance under investigation.
  • the scanning step (c) and integrating step (d) may in certain innovative variants of one aspect of the invention be optional.
  • these steps are preferably included when the desired substance is unevenly distributed within the detection microcavity. In a typically case, this may happen when the substance is retained within the detection microcavity under flow conditions, for instance by being captured to a solid phase introduced into the microcavity prior to the fluorescent substance.
  • the scanning step (c) and integrating step (d) may be replaced with collecting the radiation intensities for selected subareas of a detection area and letting these intensities represent the total amount of radiation from the detection area, for instance as a mean or maximum value.
  • Step (e) can then be carried out on these values in the same manner as for values obtained by scanning and integrating. This way of performing the method is also a part of the present invention.
  • This algorithm may be started for each scanned search area/detection area/particle bed of interest.
  • microfluidic device in a disc format, comparable to Compact Disc (CD), having an axis of symmetry perpendicular to the disc plane and with a radial flow direction from the axis of symmetry for each of the detection microcavities (see figure 3: reference number 306 a,b,c ... ) that comprises a bed of particles.
  • This kind of microfluidic device will henceforth be called CD.
  • the signal detected and measured may be fluorescence from an antibody fluorochrom conjugate (antibody-fluorochrom).
  • Each bed defines a detection area (see figure 3 : reference number 317 a,b,c ... ) in the surface of the disc.
  • the detection area is part of a search area comprising also a defined minor zone around the detection area.
  • the principles outlined in this example may be adapted also to other devices. The principles are also applicable to other kinds of radiation as discussed elsewhere in this text.
  • a controller means (fig. 1 : ref. No. 1 12) is programmed to control the scanning and detection.
  • Each search area is scanned by a detector head that is capable of detecting a detectable physical feature deriving from a substance present in the microcavity. like radiation caused by radioactivity, fluorescence, chemiluminescene, bioluminescence, enzymatic activity, scattering of light, absorption of light, reflection of light etc.
  • the measured data are stored in a data storage, an image data storage, connected to or connectable to an processing unit, which may either be comprised in the controller means, or as in another possible embodiment, incorporated in a separate computer system.
  • Each search area is identified by an identifier, which may be used as address. The identifier makes it possible for the system and controller means to locate each single search area and its detection area on a CD and link the corresponding stored measured data in the image storage.
  • each search area 402 has a pre-defined radial start position, Xstart, a pre-defined angular start position.
  • Ystart. a pre-defined angular stop position.
  • the start of the collecting of data is started at a position, primarily an angular and/or a radial position that is prior to an intended detection area, and to end the collecting at a position that is after the same detection area.
  • the search area will therefore cover an outer area 406 surrounding the detection area 404.
  • the search area 402 is scanned in radial direction by the detector disc lap by disc lap.
  • the radial start position Xstart is the number of disc laps/revolutions, counted from a spinning start position, Spinstart (not shown).
  • the radial stop position is the number of laps, counted from the spinning start position, when the last angular scanning of the search area is to be performed.
  • the scanning of the search area in the angular direction starts during the lap, which is defined by the radial start position, at the angular start position in the angular direction, which is defined by the spinning direction of the disc and the detector position in relation to a Home Mark on the CD.
  • the scanning in angular direction is stopped at the angular stop position.
  • Figure 4 is schematic in the sense that the search area of a rotatable device typically is a) slightly narrowing when going from an outer position inwards towards the centre of the CD, and b) slightly arc-shaped. Compare for instance that the Ystart position for each lap are always at the same angular position (which also applies to the Ystop position).
  • a stream of collected data from that search area will generate a corresponding stream of measured data that are stored into the storage. This is repeated lap by lap from the radial start value to the radial stop value.
  • the processing unit of the controller 120 links measured data to the correct identified detection area, and there will be no mixing of measured data from different search and/or detection areas. Each scanned search area has a corresponding set of collected measured data.
  • the processing unit will organize the data in such an order, position by position in each set of measured data, and it is therefore possible for the image processor to locate a certain measured value in the set.
  • the set of measured data has therefore similarities to an electronic and/or digitally stored image.
  • An image is defined in very small elements, called pixels. Each pixel corresponds to a position of a measurement area in the search area.
  • the measured data of a search area is possible to plot by use of printer/plotter software program used by the processing unit in a three dimensional diagram, wherein the horizontal plane is defined by the two dimensions x and y, in this case the radial direction, x-axis, and the angular direction, y-axis, which is perpendicular to the x-axis.
  • Positive x-direction is defined as the main flow direction (downstream) through the detection area, from the inlet 408 of the detection area 404 to the outlet of the same.
  • Positive y-direction is defined as the scanning direction which is orthogonal/perpendicular to the radius (or flow direction).
  • the third dimension, z-axis is defined by use of printer/plotter software program used by the processing unit in a three dimensional diagram, wherein the horizontal plane is defined by the two dimensions x and y, in this case the radial direction, x-axis, and the angular direction, y-axis, which is perpen
  • an image is the stored detector measurement data belonging to a search area/detection area/particle bed.
  • a pixel is here defined as a discrete measurement area within the search area/detection area and each pixel may have a measured intensity value, a pixel value.
  • the resolution of the image is depending of the number of discrete measurement points in the y-direction and the number of feps in the x-direction within the search area. The resolution increases with increasing density of measure points and scanning laps, and consequently, the area of a pixel will decrease.
  • the electromagnetic-emitting source is a laser
  • the possible minimum pixel area is limited by the minimum cross-section of the laser beam with the search area.
  • the pixel intensities are determined by the detection process and stored as digital pixel data representing raw data images of the search areas.
  • the whole or a part of the volume may represent a wanted measure for a chemical process, i.e. a reaction variable.
  • This value is possible to calculate by using a mathematical method, e. g. addition or integration of measurement data over the whole area of the detection area, or one sub-area 410 (fig. 4) or a number of sub-areas of the detection area.
  • the object is to calculate an amount for an subarea, here called a normalized area, i.e. a predefined amount of pixels, inside the detection area. This is done to minimize the variance between detected structures.
  • the variance depends among other things on background radiation, such as fluorescence, from the detection area.
  • a substance raw data image comprises measured substance data and noise data.
  • the noise cause irregularities, like noise peaks, in the surface, and evidently, a calculation of the volume in an area having one or many noise peaks will not result in a true measure for a chemical process.
  • the present invented method provides a method for reducing noise contribution in a substance raw data image derived from a detection area of a microfluidic device of the kind described above.
  • the method comprises a number of main-steps for reducing different kind of disturbances and noise within an ima • 'gtoe 1 -.
  • the method removes noise from substance raw data information in substance raw data images.
  • Each image is associated to one detection microcavity and one substance that is present in the detection microcavity of a microfluidic structure of a microfluidic device.
  • the substance raw data information is obtained by use of a detector arrangement capable of scanning a search area associated to a detection area, which is associated to the detection microcavity.
  • the detector arrangement is capable of detecting at least one of a number detectable physical features discussed above, e.g. fluorescence, luminescence, etc. emitting radiation from the search area, and generate output that is proportional to the magnitude of detectable activity.
  • the output is converted to measure data storable as substance raw data in a set in a data storage.
  • the method may comprise following steps in any order for removing different noise signals contributions step-by- step in any order in a substance raw data image:
  • step ⁇ - reducing background radiation
  • step ⁇ - reducing peak disturbance
  • step ⁇ - locating the detection area within the search area
  • step ⁇ - moving binary artefacts
  • step ⁇ removing unwanted areas of the detection area
  • step ⁇ applying default detection area in noisy images
  • Each of the steps ⁇ - ⁇ comprises a number of sub-steps.
  • Steps ⁇ - ⁇ may be regarded as an algorithm, but is not necessary that all of the steps ⁇ - ⁇ are carried out for all images.
  • the order and kind of steps is typically preselected and depends among others of the sensitivity and accuracy required in each individual step. In general, inclusion of more steps is required when high sensitivity is desired.
  • Step ⁇ Reducing background radiation and construction of background processed images.
  • the object is to remove the electronic, optical and dust noise that may be visible in the image.
  • This part of the invented method mainly comprises two variants.
  • the first variant a background radiation image is used.
  • the first variant is presented schematically in a flowchart. Following steps may be used:
  • step ⁇ ' 1 accessing a background radiation image or. if not available, generating a background radiation image (step ⁇ ' 1); - if available, accessing corresponding substance raw data image or, if not available, generating substance raw data image (step ⁇ '2):
  • the first step, step ' 1 is to scan the disc and measure the radiation, e.g. the intensity of the fluorescence, without presence of the substance giving raise to the measured signal from the detection areas, e.g. before step b) in the method for determining the amount of a substance, retaining the radiation associated substance before the detection areas/particle beds of the disc have been treated, e.g. with antibody-fluorochrom. All search areas/detection areas are scanned and each of the corresponding sets of measured data is stored in a background radiation image. When using antibody-fluorochrom, the corresponding image is called background fluorescence image. The image contains the measured data of the background radiation noise signal and other noise signals.
  • the radiation e.g. the intensity of the fluorescence
  • the scanning for the background radiation image is carried out as close as possibly to the step in which the substance causing the radiation is formed and/or retained in the detection microcavities.
  • the manufacturer may provide a background radiation image that for instance may be based on the average of many CDs of the same kind.
  • step ⁇ '2 is in this first variant the same as step c) in the method for determining the amount of a substance, above.
  • the substance is now present in the detection microcavities, for instance retained in particle beds that may have been washed, e.g. after the particle beds have been treated with e.g. antibody-fluorochrom.
  • This image is the substance raw data image, described above.
  • step ⁇ ' 1 and step ⁇ '2 instead of generating them again by the scanning and detecting processes.
  • the two images should be detected under the same conditions, e,g. same laser power and same detector sensitivity if the detection principles is based on fluorescence.
  • step ⁇ '3 is to correlate the two images in two dimensions, that define a correlation image, by use of a two-dimensional correlation.
  • the use of two-dimensional correlation is well-known in the art of signal-processing, see reference [1], page 90-92.
  • the step has resulted in a correlation image.
  • step ⁇ '4 the position of the pixel, or pixels, having the largest intensity value in the correlation image is/are located. This could be done by using a simple comparison algorithm, wherein the recently largest value is stored and compared with the next data, e.g.intensity value, that hasn ' t been compared with any of the earlier pixel data. The positions of the pixels having the largest pixel values, largest intensity values, are stored and used in the next step, step ⁇ '5.
  • step ⁇ '5 the background radiation image is moved to the correct position of the correlation image by using the identified pixel positions.
  • the parts of the image that overlap is subtracted pixel by pixel.
  • the substance raw data image comprise the wanted substance radiation data. If not any new noise have been introduced by the introduction of the substance giving rise to the measured signal into the microcavities (part of step b). the process of noise removing may be stopped.
  • step ⁇ determining a Median value for background radiation data (step ⁇ ' " 2); - subtracting pixel values of substance raw data image with the median value, and store background processed image (step ⁇ "3).
  • This second variant for removing background radiation, like fluorescence, especially for removing electronic and optical noise that exists in the search area outside the detection area, is possible to use instead of the five steps described above.
  • This alternative variant works very well under the condition that there is no unwanted object of large size outside the detection area and the background signal intensity is constant inside the image.
  • This variant has the advanta *gfc e that the disc is scanned onlv once, which is described hereinafter.
  • the first step, step ⁇ " 1 in this alternative variant, should start after the radiation associated substance has been incorporated into the detection microcavities.
  • the method for determining the amount of a substance above this for instance corresponds to after the detection areas have been treated with the antibody-fluorochrom, i.e. in step c.
  • the pixel values of the substance raw data images are detected and determined, and stored as digital pixel data. If the substance raw data images exist already and are availably stored, there is consequently only to copy/read the images, instead of generating them again by scanning and detecting processes.
  • step ⁇ "2 a median value of the pixel values for at least a portion of the search area outside the detection area is calculated and used as background radiation value, e.g. as a value for the background fluorescence.
  • the median value is an estimated value for the background radiation (intensity).
  • step ⁇ "3 the substance raw data is decreased with the median value contribution, e.g. by subtraction, pixel pixel, until all pixels are processed and a background processed image is provided.
  • both methods, ⁇ ' and ⁇ " provide a background processed image for each scanned search area, i.e. an image with reduced background radiation, such as fluorescence, and can be used separately or in combination in other steps of the innovative method, with the former being preferred at the moment.
  • Step ⁇ Reducing peak disturbances.
  • the method of this aspect of the invention may comprise one or more steps in which peak disturbances originating from dust, chemistry and other artefacts are reduced, e.g. reducing peaks from fibres.
  • the disturbance peaks have typically a relatively high amplitude and a relatively narrow distribution.
  • the image of an ideal detection area is smooth and does not contain any disturbance peaks.
  • step ⁇ for reducing the peak disturbances is presented schematically in a flowchart. Following steps may be used:
  • step ⁇ l - Generating a Laplacian processed image from the earlier stored image, typically a background processed image, by using a Laplacian filter in two dimensions (step ⁇ l);
  • step ⁇ 5 Determining a minimum rectangular area surrounding the start points (step ⁇ 5);
  • step ⁇ 6 Calculating an interpolated value within each peak area and substituting the original data of the peaks, and storing the disturbance filtered image (step ⁇ 6).
  • the process for reducing peak disturbances will now be described in more detail.
  • a radiation intensity image for instance a background processed image as discussed above, may be processed by a Laplacian filter in two dimensions.
  • the result of this step is that peaks in a background processed image are amplified.
  • the result is normalized with the maximum radiation (intensity) in the background image.
  • This step provides a Laplacian processed image.
  • a disturbance threshold limit is calculated by use of the standard statistic functions average and standard deviation of the stochastic background noise, calculated from the image used in step ⁇ 1 . This step enhances the disturbance peaks of the image.
  • step ⁇ 3 the Laplacian processed image is binary encoded.
  • pixel intensities/values equal to or higher than a threshold limit are set to one of the two corresponding binary values, for instance "1 " or "positive” and pixel intensities less than the disturbance threshold limit are set to the other of the binary values, for instance "0" or "negative”.
  • the determined disturbance threshold limit determined in step ⁇ 2 may be used in step ⁇ 3 for encoding. This step creates an image wherein the high intensity pixels of the peaks appears as groups, even called clusters.
  • the Laplacian processed image may then be processed for finding the pixel in each peak having the largest value by using the position information from the clusters in the binary Laplacian processed image.
  • Each pixel having the largest value in a cluster is defined to be a center and start point. The other pixels of the clusters are set low in the binary Laplacian processed image. This step results in a cluster image.
  • step ⁇ $ the start points defined in the cluster image is used for determining a minimum rectangular area by iterative testing using different area sizes.
  • the iteration that works on the Laplacian processed image, is stopped when a minimum rectangular area contains only pixels having a negative intensity value or zero value according to the definition in previous steps and as a result of the earlier performed Laplacian filtering. Said pixels surround the corresponding start point and defines the position, the area of the corresponding disturbance peak and the boarder pixels which is outside the peak area.
  • step ⁇ 6 the disturbance peaks of the background processed image are removed by substituting the original data within each peak area with an interpolated value between the boarder pixels of the corresponding peak.
  • the interpolated values are calculated by use of an appropriate interpolation function, both a linear and a non-linear function is a possible choice.
  • the process has now generated a "disturbance filtered image" of the radiation intensity, e.g. fluorescence intensity, of the detection area.
  • Step ⁇ Determining a global threshold value
  • the object of this part of the invention is to make a radiation intensity image, such as a disturbance filtered image, binary by determining a global threshold that separates the measured substance signals, i.e. the signals from the substance in the detection area, from the background signals.
  • a global threshold that separates the measured substance signals, i.e. the signals from the substance in the detection area, from the background signals.
  • the detection area pixels should be equal to binary high, like "1", and the rest of the image equal to binary low, like "0".
  • a global threshold is valid when only a few pixels in a part of the search area that is not connected to the detection area have an intensity higher than the pixels in the detection area.
  • step ⁇ for determining a global threshold value is presented schematically in a flowchart. Following steps may be used:
  • step ⁇ 4 Determining radial directed edges of detection area (step ⁇ 4); Generating a x-summation graph (step ⁇ 5);
  • step ⁇ 9 Determining the start and end positions of the detection area
  • step ⁇ lO Determining a start value of the global threshold value
  • step ⁇ l pixel values are summarized in the radius direction, pixel values, line by line, in a radiation intensity image, typically the disturbance filtered image, to get a one-dimensional y-summation graph in the y-direction of the radiation intensity, e.g. fluorescence intensity (angular direction on a spinning disc).
  • a radiation intensity image typically the disturbance filtered image
  • a one-dimensional y-summation graph in the y-direction of the radiation intensity, e.g. fluorescence intensity (angular direction on a spinning disc).
  • Each y-coordinate on the graph will be the sum of the pixel values in radial direction having the same y-coordinat.
  • step ⁇ 2 the y-summation graph is filtered to remove disturbing peaks, e.g. with a median filter function.
  • step ⁇ 3 the centre of the detection area in y-direction (angular direction on a spinning disc) is calculated and located by use of a correlation, and with the filtered y-summation graph and the known detection area width as input value to the correlation length. This should normally give a graph with one large peak.
  • step ⁇ 4 the centre of the largest peak is set to the centre of the detection area. It's now possible to calculate the positions and the radial extension of the particle bed edges by use of the known data for the centre and detection area/particle bed width. Hereby, the radial directed boarder/edge pixels are known.
  • step ⁇ 5 the radiation intensities . are summarized in the image in y- direction (angular direction on a spinning disc) using the limits (boarders) determined in the previous step to get a x-summation graph (radial summation graph) in x-direction (radius direction) of the radiation intensities.
  • the summation of pixel values is performed in the same way, but in the orthogonal direction in the same plane, as in step ⁇ l .
  • step ⁇ 6 the x-summation graph is filtered with a median filter function to remove disturbing peaks.
  • step ⁇ 7 the start of the detection area in the x-direction is determined by use of correlation between a step function and the filtered x-summation graph.
  • the mathematical process results in an edge curve, see figure 9.
  • step ⁇ 8 disturbance peaks, close to edge that is sitting opposite to the inlet, are restrained by multiplying the edge curve with a ramp function.
  • step ⁇ 9. The purpose of step ⁇ 9. is to find the start position of the detection area by finding the highest positive peak centre in the filtered edge curve (from step ⁇ 8) and using a minimum detection area length to calculate the end position of the detection area.
  • step ⁇ lO the global threshold value is now determined by use of optimal threshold technique for image processing (reference [1], starting at page 360). It is assumed that the background and particle bed peaks are normal distributed with the same variance. Calculate from the disturbance filtered image the estimated background intensity as max of (median background intensity on the areas A1+A3 or median intensity background on the area A2)). The positions and orientations of the areas Al . A2 and A3 is illustrated in figure 13. Calculate the estimated substance intensity as median of the detection area. Calculate the start value for the iteration step as average of the estimated background intensity and the estimated substance intensity.
  • step ⁇ l 1 the global threshold value is iterative adjusted, between a maximum intensity value, defined by the estimated substance intensity, and a min value, defined by the estimated background intensity, so that the average number of pixels in the y-direction (the angular direction) inside the approximate position of the detection area is equal to the detection area width.
  • step ⁇ l2 the calculated global threshold value is used for binary encoding a stored image, in particular the stored disturbance filtered image into a binary coded disturbance filtered image.
  • Step ⁇ Remove binary objects outside the detection area
  • the object of the following steps ⁇ l - ⁇ 2 is to remove any disturbances outside the detection area.
  • the following algorithm processes the data of the binary coded disturbance filtered image, in which all objects are set to 1, i.e. binary high. Further, it is assumed that there is no object to be removed inside the boarders of the approximately calculated detection area position, and the detection area is the largest area of the image.
  • step ⁇ for removing binary objects outside the detection area is presented schematically in a flowchart.
  • step ⁇ l labelling of all the objects in the binary coded disturbance filtered image, wherein the number of connected pixels, i.e. the area, of each object is determined. Labelling is a known method in image processing.
  • step ⁇ 2 generating a binary detection area image by removing all but the largest objects that are inside the approximate detection area position.
  • Step ⁇ Removing unwanted areas of the detection area
  • the object is to calculate an amount for a ""normalized " area. i.e. a predefined amount of pixels, inside the detection area. This is done to minimize the variance between detection areas on the same microfluidic device.
  • the variance depends among other things on background radiation, such as fluorescence, from the detection area. It is assumed that the 3D distribution of an intensity signal inside the detection area is independent of small variations in the detection area length.
  • the algorithm used is a process for generating a new binary image that only comprises pixels within the normalized area ((detection area width) multiplied by (min detection area length)) in the binary detection area image starting from the inner part of the CD.
  • step ⁇ for removing unwanted areas of the detection area is presented schematically in a flowchart.
  • step ⁇ l generating a normalized area image by setting all binary high indicated pixels outside the location of the normalized area ((detection area width) multiplied by (min detection area length)) to binary "low” in the binary detection area image.
  • step ⁇ 2 multiplying each corresponding pixel in the disturbance filtered image, generated and stored earlier in the invented process, and the normalized area image with each other and perform a calculation, e.g. integration, on the normalized area.
  • Step ⁇ Apply default detection area in noisy images
  • the described method above will not work properly, or not at all, for very noisy images hiding the substance signal.
  • the following algorithm works, if all processed images originate from the same set of detection areas, e.g. the same CD.
  • the detection areas are then well-defined in angular position and quite well in radius direction.
  • step ⁇ for applying a default detection area in noisy images is presented schematically in a flowchart.
  • the method for applying default detection area in noisy images may comprise following steps: Step ⁇ l : Saving automatically all determined detection area positions (step ⁇ ) fulfilling high confidence criterion.
  • Step ⁇ 2 Calculating an average detection area position.
  • Step ⁇ 3 Using automatically the average detection area position if one or more predetermined condition is fulfilled.
  • the average position for a detection area will automatically be applied for all calculated detection areas that fulfil all following conditions:
  • the number of binary high pixels, in the normalized area image, in angular direction for all radiuses is less than [amount] multiplied by [detection area width] multiplied by [min detection area length], wherein the [amount] is chosen 0-1.
  • Steps d) and e) in the method for determining the amount of a substance that is associated with the collected radiation as well as any of the other individual method steps or their combinations discussed are typically performed by the appropriate software for instance included in the controller or elsewhefe, for instance not in physical association with the innovative arrangements described herein.
  • the invention also comprises computer program-related aspects for treating radiation data that have been assembled by steps (a)-(c) of the innovative method.
  • One such aspect is a computer program product that (1) comprises program code elements corresponding to a sequence variant of step (d) above comprising one or more of the substeps described for this step, and
  • step (c) when installed on the appropriate hardware is capable of causing the hardware (computer) to execute the sequence of substeps on data which have been obtained by performing step (c) above on any kind of microfluidic device irrespective of being spinnable or not.
  • code elements may be included, for instance corresponding to step (e) in order to execute the sequence (d)-(e).
  • Another computer program-related aspect is the computer program product stored at a computer program readable means which, when the product is loaded, makes it possible for a computer to perform the sequence of steps corresponding to the code elements in the stored computer program product.
  • a third computer program- related aspect is a carrier having at least one of the innovative computer program products thereon.
  • the carrier may be a computer memory, a Read-Only Memory or an electrical signal carrier.
  • the present invention may be implemented as a computer program product directly loadable into an internal memory storage of a processing unit within the controller in the detector arrangement, comprising the software code means for performing the steps of any claims 1-16.
  • the present invention relates to a computer program product stored on a computer usable medium, comprising readable program for causing a processing unit in a computer means to control an execution of the steps of any of the claims 1-16.
  • the computer usable medium is a record medium, a computer memory, a Read-Only Memory or an electrical carrier signal.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Plasma & Fusion (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Image Processing (AREA)

Abstract

The present invention concerns a method and an arrangement for reducing noise in a substance raw data image containing noise. The present invention also relates to a method and an arrangement for determining a measure of at least one substance. Each image is associated to one detection microcavity and one substance that is present in the detection microcavity of a microfluidic structure of a microfluidic device. The substance raw data information is obtained by use of a detector arrangement capable of scanning a search area associated to a detection area, which is associated to the detection microcavity. The method may comprise following steps in any order for removing different noise signals contributions step-by-step in any order in a substance raw data image : - reducing background radiation (step α); - reducing peak disturbance (step β); - locating the detection area within the search area (step χ); - moving binary artefacts (step Δ); - removing unwanted areas of the detection area (step ϵ); - applying default detection area in noisy images (step φ).

Description

Method and arrangement for reducing noise in a microfluid device
BACKGROUND OF THE INVENTION
I. Field of the Invention
The present invention generally concerns the field of microfluidic devices. More particularly, the present invention concerns a method and an arrangement for reducing noise in a substance raw data image containing noise. The present invention also relates to a method and an arrangement for determining a measure of at least one substance. The substance raw data image has been obtained from a search area that comprises a detection area associated with a detection microcavity that is part of a microchannel structure of a microfluidic device detectable signals.
The substance raw data image reflects from a substance that is present in the detection microcavity after one or more liquid aliquots have been processed through the microchannel structure that comprises the detection microcavity. See further below.
II. Related Art
All patents and patent applications cited in this specification are incorporated by reference in their entirety.
A. Background Publications
Detector arrangements for measuring radiation signals from individual detection areas on a circular substrate have been described in a number of previous publications. See for instance: EP 392475 (Idemitsu Petrochemical Co, Yamaji Kazutaka et al.) US 5994150 (Imafion Corp, Challener et al), US 5892577 (The University Court of the University of Glasgow. Gordon), WO 9721090 (Gamera, Mian et al), Duffy et al., "Microfabricated Centrifugal Microfluidic systems: Characterization and multiple Enzymatic Assays" (Anal. Chem. 71 (1999) 4669-4678), WO 0040857 (Amersham Pharmacia Biotech AB. Bjδrkesten et al)
B. Background Technology and Problems
The present invention belongs to the field of miniaturization of processes which comprise sample treatment, assay protocols, chemical and/or biochemical synthesis etc., within medicine, chemistry, biochemistry, molecular biology and the like. At present one important goal within this field is to reduce the costs for these processes, for instance to reduce the amount of reagents needed per assay, reduce time per assay, etc. One route has been to increase the degree of parallelity, for instance by integrating as many as possible of similar process runs in one and the same device in order to carry out all the runs in parallel. At present, large numbers of research groups and companies are involved in developing technology that will solve the numerous problems encountered.
One problem is related to the optimal way of configuring the detector in relation to the microdevice used for performing the processes while maintaining an acceptable sensitivity and reproducibihty. This problem may become particularly pronounced if the measuring step is performed by continuously moving the detector unit and the detection areas of a microfluidic device relative to each other during the measurement operation.
Another problem arises if the microfluidic device is in the form of a disc which is skewed because then it becomes difficult to maintain the optical focus in the right position relative to the detection areas. Without proper arrangement skewed discs will reduce sensitivity. This problem in particular applies to discs made of plastic material.
Another problem is related to maintaining an acceptable sensitivity and reproducibihty when changing sample volumes from the μl-range to the nl-range and performing the process protocols with a high degree of parallelity within the same device. The inventors have found that under these circumstances the materials from which the microfluidic devices are fabricated and the various treatments during the manufacturing and conditioning of the devices easily introduce signal artifacts that are of the same kind and of comparable or larger size as the desired signals.
During recent years it has become popular to fabricate microfluidic devices in plastic material. This kind of material is typically highly fluorescent ("auto-fluorescent") with emission wavelengths covering most of the wavelengths normally utilized in fluorescent measurements. Compared to microtitre wells and other uncovered microstructures the problem becomes more severe for the kind of covered microchannel structures used in the present invention, because the exciting and emitted radiation has to pass through plastic material. For transparent plastic material there is also a problem with "cross-talks" between the detection area/detection microcavities. Similar problems may also be at hand for spectroscopic methods in which the radiation to be measured is created within the detection microcavity (for instance chemiluminescence, bioluminescence, etc.).
A more recent problem relates to the fact that the present assignee recently has managed to control the liquid flow in microfluidic devices containing a plurality of microchannel structures in such a way that the inter-channel variation for a device with respect to flow becomes insignificant. This progress has enabled the assignee to quantify with a low inter-assay variation and a high sensitivity analytes, such as antigens, in the subfemtomole range in nl- volumes by carrying out the solid phase reaction of a heterogeneous sandwich immunoassay under flow conditions in small columns (nl-columns). This has raised the question about measuring the amount of an affinity complex such as an immune complex as a function of position along the flow direction in a column. See WO 02075312 (Gyros AB) and assignee's poster presented on September 17, 2001 at Proteomic Forum September 16-19, 2001, Munich, Germany.
BRIEF SUMMARY OF THE INVENTION
The present invention is directed to a method which meets the discussed problems.
A first object is to reduce the electronic, chemical, optical and dust noise that is present in the stored measurement data of a search area, i.e. that is present in the stored image of a detection area. Electronic noise is introduced by the detector arrangement instrument. The optical noise arises from the design of the disc, the lid of the disc (if present), the detector head (e.g. lenses) etc. The chemical noise arise from chemicals and materials used in the manufacture of the disc and from undesired reactions including undesired adsorption related to the various reagents and solutions that have been used for placing the signal creating substance in the detection microcavity. Peak noise may derive from the manufacturing or the processing of the disc. Different kinds of dust, like particles, fibres, hair, etc may also introduce peak noise. The invention typically reduces noise-from steps related to chemical noise, e.g. disturbing fluorescence in the case the substance raw data image is based on fluorescence.
A second object is to provide a method, software and arrangement for treating radiation data from one or more detection areas of the kind of microfluidic device described herein. A third object is to provide arrangement, software and methods enabling accurate integration of radiation deriving from a desired substance as a function of subareas of individual detection areas. In particular this object aims at avoiding the problems discussed herein which the inventors have found may appear when measuring low amounts of substances with a high accuracy in microfluidic devices.
These objects in particular apply to radiation data that have been collected while the detection area(s) is(are) moved in front of the detector head that is used for collection of the data, or vice versa, e.g. the detection area(s) is/are present in the surface of a circular disc that is spun in front of the detector head.
The foregoing has outlined rather broadly the features and technical advantages of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of the invention will be described hereinafter which form the subject of the claims of the invention. It should be appreciated by those skilled in the art that the conception and specific embodiment disclosed may be readily utilized as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims. The novel features which are believed to be characteristic of the invention, both as to its organization and method of operation, together with further objects and advantages will be better understood from the following description when considered in connection with the accompanying figures. It is to be expressly understood, however, that each of the figures is provided for the purpose of illustration and description only and is not intended as a definition of the limits of the present invention.
BRIEF DESCRIPTION OF THE DRAWINGS
For a more complete understanding of the present invention, reference is now made to the following descriptions taken in conjunction with the accompanying drawing, in which:
Figure 1 illustrates a schematic view on an arrangement of the invention and its main parts.
Figure 2 illustrates a detector head placed above a disc (microfluidic device) (cross- sectional view). Figure 3 illustrates a set of microchannel structures that can be used in a circular disc.
Figure 4 illustrates the search and detection areas.
Figure 5 is a flow chart illustrating a first variant, α\ of a method step α according to the invention.
Figure 6 is a flowchart illustrating a second variant, α", of the method step α according to the invention.
Figure 7 is a flowchart illustrating a second variant, α", of the method step α according to the invention.
Figure 8 is a flowchart illustrating a method step χ according to the invention.
Figure 9 illustrates an edge curve.
Figure 10 is a flowchart illustrating a method step δ according to the invention.
Figure 11 is a flowchart illustrating a method step ε according to the invention.
Figure 12 is a flowchart illustrating a method step φ according to the invention.
Figure 13 illustrates the search and detection areas and the surrounding areas A1-A3.
DETAILED DESCRIPTION OF THE INVENTION
It is readily apparent to one skilled in the art that various embodiments and modifications can be made to the invention disclosed in this Application without departing from the scope and spirit of the invention.
As used herein, the use of the word "a" or "an" when used in conjunction with the term "comprising" in the sentences and/or the specification may mean "one," but it is also consistent with the meaning of "one or more," "at least one," and "one or more than one."
As used herein, the term "rotating" refers to spinning. Yet further, the term "rotating" may also include, but is not limited to a step-wise rotation of a disc. As used herein, the term "reagent" includes, but is not limited to an analyte.
As used herein, the term "in a circular manner" refers to around the center of the disc (circumferential direction).
The term "a plurality of microchannel structures" means two, three or more microchannel structures. Typically the term "plurality" means > 10, such > 50 or > 100 microchannel structures.
The terms "removing noise" and "reducing noise" means that the noise contribution in the measured data (values) is reduced to zero, or is negligible. The words "reduce" and "remove" will be used alternately in the text, but there is no difference in the meaning of the two words in this specification. They should be regarded as equivalent.
As illustrated in figure 1 , a detector arrangement is typically adapted for measuring radiation from a plurality of detection areas (317a.b.c. etc., in figure 3) each of which is associated with a detection microcavity (306a,b,c, etc., in figure 3) in a microfluidic disc (101 in figure 1, 301 in figure 3). With reference to figure 1, the arrangement comprises:
(a) a detector head (102) with a focal area, and a disc holder (105) which are linked to a means enabling for the detector head (102), i.e., the focal area to transverse, the surface of the disc (101) in an essentially circular (means I) and/or radial manner (means II) when the disc is placed in the disc holder (105)
(b) an angular aligning system (108,109) for recognizing the angular position of the part area which at a particular time is covered by the focal area, and
(c) an optional radial aligning system (1 10,11 1) for recognizing the radial position of the part area which at a particular time is covered by the focal area, and
(d) a controller (112), e. , computer with software, which controls
(i) means I and means II causing the focal area to transverse the detection areas (317a',b',c\ etc.) in an annular zone of the disc, and (ii) the detector head (102) successively collects radiation in a preselected manner from individual subareas of essentially the same size as the focal area within at least one of the detection areas in said annular zone.
In specific aspects, it is contemplated that the disc holder (105) is linked to at least one of the means for enabling the detector head and the disc to move relative to each other so that the focal area transverses the surface of the disc. Thus, the disc holder (105) is linked to either means I or means II or both means I and means II.
Yet further, it is also understood by those of skill in the art that means I, means II, and the detector head can be linked directly or indirectly to the controller.
Means I comprises three main variants with respect to how the disc and detector moves:
(1) The disc is rotating around its axis of symmetry. In this variant, means I preferably comprises a spinner motor (103) and a shaft (104) carrying the disk holder (105). See figure 1.
(2) The detector head (focal area) is circularly moving around the axis of symmetry of the disc. In this variant, means I preferably comprises a spinner motor and a shaft with an arm that carries the detector head. This shaft typically has the same direction as the axis of symmetry of the disc when placed on the disc holder.
(3) A combination of (1) and (2).
Means II comprises three main variants with respect to how the disc and detector head moves:
(1) The disc is laterally moving in front of the detector head. In this variant, means II preferably comprises a motor for laterally moving the disc holder.
(2) The detector head (focal area) is laterally moving over the disc surface. In this variant, means II preferably comprises a linear frame (106) that carries the detector head, and a drive unit (107) for radial movement of the focal area, over a disc which is placed in the disc holder (105). See Figure 1.
(3) A combination of (1) and (2). The arrangement is illustrated in figure 1.
I. Means I and the Angular Aligning System for the Detector Head
In a typical variant, the detector head (102) and the motor (103) (e.g., a spinner) with a rotatable shaft (104) carrying a disc holder (105) are supported on a frame structure (113). The motor (103) controls the rotating speed that can be varied, e.g., within an interval between 0 - 15,000 rpm, such as within an interval between 60-5,000 rpm. The rotation of the disc may be stepwise. When the disc is rotating, the focal area of the detector head will successively scan angularly adjacent part areas of the disc surface. Means I in figure 1 is according to variant (1) above and comprises the motor (103) and the shaft (104).
The disc holder (105) is preferably a plate on which the disc can be placed. The disc holder could also be a device that holds the disc at its periphery. In order to reduce wobbling of the disc (if the disc is skewed), the side (114) of the plate facing the disc may comprise a system of evenly distributed uncovered shallow grooves or channel openings that are connected to a vacuum system by which the disc can be sucked to the plate. See for instance PCT/SE02/01682 (Gyros AB).
Detection areas may be present in one or both sides of the disc. If the disc holder is in the form of a plate as illustrated in figure 1 and the microfluidic disc has detection areas on the side opposing, it is important to secure the plate so as not to disturb the collection of radiation. The plate may thus have a smaller diameter than the disc with the detector areas being located in an annular zone that is not covered by the plate. Alternatively, the plate may be in a material that is translucent for the radiation utilized for the measurement, at least at the detector areas.
The angular aligning system may comprise:
(1) a device that will enable the determination of when a predetermined angular position of a disc placed on the disc holder is in front of the objective of the detector head (102) (i.e., covered by the focal area), and/or
(2) a home position mark detector (108) which is able to detect when a home position mark (305) on a rotating disc (101) placed on the disc holder (105) is passing.
A home position mark (305 in figure 3) is preferably placed in an outer circumferential zone outside the detection areas or in some other position, which can be detected with high accuracy. The position coordinates of each specific spot of the surface of the disc is given as the angular position relative to the home position mark and as the radial position relative to the circumference or axis of symmetry or relative to any other arbitrary fixed position on the disc.
A home position mark detector (108) typically has a fixed position outside the disc, for instance on the frame structure (113).
An accurate and preferred alternative for determining when a predetermined angular position is in front of the objective is to include an encoder that progressively gives the angular distance from the home position mark while the disc is rotating. This kind of encoder ( 109) is typically associated with means I, for instance the motor (103), the shaft (104) or the disc holder (105). By associating the encoder directly with the disc (101) it is likely that the most accurate determination will be accomplished. The encoder typically divides each revolution of the shaft into a large number of grades, for instance > 5 000, such as > 10 000 or > 20 000 or > 30 000. A simple but less accurate alternative is to include calculating means that calculates the time needed from a preset rotation speed and the angular distance between the predetermined position and the home position mark (i.e., from the preset rotation speed and the angular position coordinate). This kind of calculating means may be associated with the controller.
The angular aligning system should be able to give the angular position coordinate for the part of the disc which is in front of the objective with an accuracy of ± 1°, such as within ± 0.1° or within ± 0.01° (provided there are 360° per revolution). The exact accuracy needed will depend on the size of the disc, radial position of the detection area, the required sensitivity, size of detection area, etc.
II. Means II and the Radial Aligning System for the Detector Head
The detector head (102) is guided on a linear frame (106) that may be the upper part of the frame structure (113) for linear displacement and positioning in a first plane PI, transversely through the central axis CL of shaft (104) and running in a radial direction thereto. The linear frame (106) prohibits uncontrolled movement of the detector head in any other direction relative to this linear displacement. The drive unit (107) for this displacement may be in the form of a translational responder for incrementally changing the position of the detector head (102) in the first plane PI (radial movement) and for enabling scanning of radially adjacent subareas of a microfluidic disc device placed in the disc holder (105). Means II comprises the linear frame (106) and the drive unit (107) and is in figure 1 according to variant (2).
The drive unit (107) is associated with a unit (110) for determining the linear displacement and thus also the radial position coordinate of the part of the disc which is in front of the objective of the detector head (102). This unit (110) may be in the form of an encoder that gives a translational position and movement of the focal area (objective of the detector head) in relation to a translational home position (111) on the responder. This home position in turn may be associated with a unique radial position in the disc. The measuring unit (110) should be able to translate a translational position and movement of the detector head (focal area) into a radial position coordinate of the disc used with a high accuracy, typically within ± 10 μm such as ± 1 μm or ± 0.1 μm.
The drive unit (107) and the vertical height of plane PI may be adjustable for focusing purposes.
III. Controller Control means, for instance electronic and programmable control means
(schematically illustrated by reference numeral (112)) with operator's interface and software, not further disclosed, may be assigned to the detector arrangement among others for
a) recognizing one or more pairs of start/stop-positions (angular and/or radial) for irradiating if the detection principle utilized requires irradiation and/or for collecting radiation,
b) identifying individual subareas in detection areas or elsewhere in the surface of the disc,
c) controlling the simultaneous rotating of the disc and the incremental radial displacement of the detector head (102),
d) collecting radiation data from the detection areas/detection microcavities.
e) treatment and presentation of the collected data, and/or
f) determining the time at which a particular angular position is in front of the objective of the detector head from the rotational speed. g) storing and reading collected data in a data storage 130,
h) storing a computer program/computer program product, i.e. a computer software application in a computer software application and computer program storage 120, and reading and accessing software applications from the same.
The present invention may be implemented as a computer program product, i.e. a computer software application, directly loadable into an internal memory storage of a processing unit within the controller 120. Said computer program product comprises the software code means for performing the steps of the invented methods.
Further, the present invention relates to a computer program product stored on a computer usable medium, comprising readable program for causing a processing unit in a computer means, such as the controller 112, to control an execution of the steps of the invented methods.
The computer usable medium is a record medium, a computer memory, a Read-Only
Memory or an electrical carrier signal.
Different parts of the arrangement may communicate (1 15) with the controller (1 12). The controller will in the preferred variants instruct the detector head to successively collect radiation from distinct and preselected parts of the surface of the disc. Typically the controller is programmed to start collecting radiation at a position, primarily an angular and/or a radial position, which is prior to an intended detection area, and to stop collecting at a position, which is after the same detection area. The start and stop positions define a search area that at least should encompass the detection area or the relevant part of the detection area in which there is radiation from the substance. See further figure 4. If the radiation requires that the substance is irradiated, which is the case if-fluorescence is measured, the control means also defines the settings for the start and stop positions for irradiation. These latter settings are typically essentially the same as for collecting radiation.
The start and stop signals for collecting radiation is preferably directly linked to the angular positions in the disc at which collection is to start and end. respectively. This also includes that due account is taken for delays that may be inherent in the system or preset, i.e., the start and stop signals may have to be initiated before the focal area is positioned in front of the start and stop position, respectively. If the angular aligning system comprises an encoder, the encoder signals corresponding to a start position and a stop position are used to define the time
'period during which radiation is to be collected. In an alternative, the start and stop for collecting radiation is linked to a preset rotating speed, i.e., the controller calculates from a preset rotating speed the time at which the start and stop position should be in front of the objective.
The controller may be programmed to change the radial position of the detector head (focal area) after a predetermined number of revolutions of the disc, for instance after 1. 2 or more revolutions with preference for 1.
The controller may be capable of changing the radial position of the objective (focal area) during a revolution of the disc. For these variants, one can envisage that radiation is collected for all relevant subareas at a common angular position before the disc is rotated (in a single step) to a subsequent angular position. In an alternative variant the objective (focal area) is transversing the disc surface in a spiral-like manner, i.e , the radial position is changed successively during a revolution.
In the preferred variants, the collected radiation data is stored in a data storage 130 in a form that is retrievable for each individual subarea, for instance in the control unit. This means that after collecting of radiation, it will be possible to represent the collected data as a 3-D image of each search area and/or detection area showing the amount of radiation from each individual subarea. In the case of overlapping subareas, the proper treatment of the data takes into account the overlapping effect and creates a true image of the radiation associated with different parts (subareas) of a detection area.
IV. Detection Principles and the Detector Head
The microfluidic devices used in the present invention typically require detection of very low absolute amounts or concentrations of substances in the detection microcavity. It is therefore imperative in many variants of the invention that the detection principle shall enable detection and quantification of substance amounts that are < 10" " mole per detection area/detection microcavity, such as < 10"'""' mole or < 10"1 mole per detection area/detection microcavity. Typical detection principles may include, for example, spectrometric detection, such as radioactivity, absorption of radiation (e.g. light), fluorescence, chemiluminescence, bioluminescence, scattered light etc. For variants in which radiation requires irradiation a change in wave-length(s), polarization, life-time, scattering, intensity, etc., between irradiation and radiation may be measured as a function of the presence of the substance of interest in the detection microcavity. Typical examples include, but are not limited to luminescence and fluorescence principles, with laser induced fluorescence (LIF) being preferred.
Details about the detector head are given in copending applications PCT/SE02/01678 and corresponding US application S.N. 10/062,258.
The size of the focal area of the detector head typically is less than the size of a search and/or a detection area. The width of a focal area in one direction (direction 1) is typically essentially < 1/5, such as of the < 1/10, of the corresponding width of a search and/or detection area and in a perpendicular direction (direction 2) within the same limits or larger. A suitable detector head in the form of a pick-up head (200) is illustrated in figure 2 and comprises a laser source (201) whose beam is reflected on a dichroic mirror (202) and focused through an objective (205) to a part of a detection microcavity (203) positioned in front of the head. The epi- fluorescent light is passed through the dichroic mirror and through a band-pass filter (206), selective for the flourochrome at hand and is finally focused onto the entrance of a photo multiplicator tube (PMT) (207) by means of an aspheric lens (208). Pin-holes (209 and 210) are positioned between the entrance of the laser beam and the dichroic mirror, and in front of the PMT (207). The size and position of the pin-hole (210) are adapted so that the focal area of the laser beam is inside the detection microcavity. The size of pin-hole (209) is adapted so that preferentially emission light emanating from the focal area is passed into the PMT (207). A mirror spot on a glass disc or the like may replace the dichroic mirror.
V. Microfluidic Device
The microfluidic device used in the various aspects of the invention comprises a plurality of microchannel structures in which aliquots of liquids are transported and/or processed.
The devices typically are disc-shaped. The structures are covered in the sense that their interior is in contact with ambient atmosphere primarily via separate inlet and/or outlet openings and/or vents. Each microchannel structure comprises one or more detection microcavities and possibly also one or more reaction microcavities, and microconduits connecting these parts with each other. A reaction microcavity may coincide with a detection microcavity. The result of the processing in a microchannel structure is measured as radiation from a detection area which is directly or indirectly associated with a detection microcavity.
A disc typically has an axis of symmetry (Cπ) perpendicular to the disc plane where n is an integer > 2, 3, 4, 5 5 preferably ∞ (C). The disc may thus be rectangular or have a circular shape.
Different principles may be utilized for transporting the liquid aliquots within the microfluidic device/microchannel structures. Thus inertia force may be used, for instance by spinning the disc (centrifugal force). Other forces are capillary forces, electrokinetic forces. hydrodynamic forces etc.
Microfluidic devices that have an axis of symmetry and are intended for rotation may have a home position mark as discussed above.
The terms "microchannel", "microconduit", etc., contemplate that a channel structure comprises one or more cavities and/or channels/conduits that have a cross-sectional dimension that is < 10 μm, preferably < 10" μm. The lower limit is typically significantly larger than the size of the largest reagents and constituents of aliquots that are to pass through a microchannel. The volumes of microcavities/microchambers are typically < 1000 nl, such as < 500 nl or < 100 nl or < 50 nl or < 25 nl, which in particular applies to the detection microcavities. Chambers/cavities directly connected to inlet ports for liquids may be considerably larger, e.g.. microchambers/microcavities intended for application of sample and/or washing liquids. Microformat means that one, two, three or more liquid aliquots that are transported within the device have a volume in the μl-range, i.e., ≤ 1000 μl such as < 100 μl or < 50 μl including but not limited to the nl -range (nanoformaf). such as < 1000 nl or < 500 nl or < 100 nl or < 50 nl.
The microfluidic device may be made from different materials, such as plastics, glass, silicone polymers, etc. The detection area should be transparent/translucent for the detection principle utilized by the detector. From the manufacturing point of view plastic material is many times preferred because the costs for the material are normally low and mass production can easily be done, for instance by replication. Typical manufacturing processes involving plastic material are replication by embossing, moulding, etc., followed by attaching a top lid covering the open microchannel structures so obtained. See for instance WO 91 16966
(Pharmacia Biotech AB, Ohman & Ekstrδm). However, plastic materials may interfere with several sensitive detection principles. Their high auto-fluorescence is disadvantageous for normal fluorescence techniques in case low absolute amounts of fluorescent substances are to be measured. This points to the fact that it is important to match the material in the device with the detection principle used. At the priority date of this invention the preferred disc material is plastic material, such as polycarbonates and plastic material based on monomers which consist of a polymerisable carbon-carbon double or triple bonds and saturated branched straight or cyclic alkyl and/or alkylene groups. Typical examples are Zeonex and ZeonorI M from Nippon Zeon, Japan, with preference for the latter. See for instance WO 0056808 (Gyros AB, Larsson, Ocklind and Derand) which is hereby incorporated by reference. In this context silicone polymers such as poly dimethyl siloxane (PDMS) and the like are not considered to be plastic material.
It is known that black plastic material, for instance containing graphite powder or carbon black, absorbs light and therefore has a low auto-fluorescence. Transport of light within black plastic material is prevented. Black plastic material will be very efficient for microfluidic devices when fluorescence and luminiscence measurements are relied upon. Black plastic material should in this case be avoided in the detection areas.
In case a lid is needed as for most discs obtained by replication, plastic lids of different or iigin may be used, for instance Melinex™ 12 PET and Melinex™ 17 OPP (Du Pont, U.S.A.), etc
From the auto-fluorescent point of view an optimal combination of transparent plastic material appears to be ZeonorI M for replication and Melinex I M 17 OPP as the lid. This combination of material is likely to be useful for excitation wavelengths in the interval 480-650 nm.
The plurality of dejection microcavities and the corresponding search and/or detection areas are preferably arranged in subgroups such that all members of a subgroup are positioned at the same radial distance and/or at the same angular position and/or have equal length and/or cross-sectional dimensions. Within each subgroup there may be at least two, three or more detection microcavities (detection areas), such as > 10 or > 25 or > 50 detection microcavities (detection areas). A detection area in the inventive arrangement typically has a size within the range of lxlO2 - 2x10° μm2, such as lxlO - lO3 μm2. Their length and/or breadth are typically within the range of 0.5x10 - 5xl04 μm, such as 1x10 -104 μm.
Design microchannel structures (301a,b.c. etc.) that are adapted to a circular disc are illustrated in figure 3. The structures are linked together by a common distribution channel (302) and a common waste channel (303). The orientation of the microchannel structures around a common axis of symmetry is apparent. The circumference (304) of the disc has a home position mark (305). Each of the combined reaction/detection microcavities (306a,b,c, etc.) is communicating in the downstream direction with the common waste channel (303) and in the upstream direction via separate connections with the common distribution channel (302) and separate volume measuring units (307a,b,c, etc.). A surface detection area (317a,b,c, etc.) is associated with each detection microcavity. The common distribution channel (302) carries at one of its ends and at an intermediary position inlet ports (308 and 309. respectively). Another kind of inlet port (310) is located at each volume measuring unit (307a,b,c, etc.). Each microchannel structure (301a,b,c, etc.) also has an outlet to the common waste channel (303) and an outlet port (318) at the remaining end of the common distribution channel (302). An inlet vent (31 1) to ambient atmosphere is connected to the common distribution channel via a common venting channel (3E2). Other vents (313 and 314) to ambient atmosphere are placed in the common waste channel (303) and in the connecting microconduit between each reaction/detection microcavity (306a,b,c. etc.) and the common waste channel (303). Appropriate valvings are positioned at 315 and 316 in each microchannel structure (301a,b,c, etc.).
The diameter of the preferred circular discs is about the same as conventional CDs but may be larger or smaller, for instance up to 300 % or larger and down to 10 % or less.
VI. Processes to be Preformed within the Microfluidic Device The processes that are carried out within the individual microchannel structures comprise assay protocols, organo-themical or biochemical synthesis protocols, etc. Typically the protocols comprise introduction of one or more liquid aliquots containing the necessary reagents/reactants into a microchannel structure. In the case of assay protocols, one of the aliquots is a sample which is uncharacterized with respect to at least one feature, e.g.. type, form and/or amount of an analyte. The processes comprises that the substance, which is associated with the radiation to be collected, is formed and/or retained in the detection microcavity under static conditions or under flow conditions. The reaction system for retaining may be homogeneous or heterogeneous, i.e., with or without the desired substance being partitioned between a liquid phase and a solid phase. In case of flow conditions and microfluidic devices in the form of discs, the flow direction in the detection microcavity may be towards the circumference (outwards) or towards the center of the disc (inwards), or essentially parallel with the circumference of the disc. Also other directions may be utilized.
Process protocols may utilize specific reactions between reactants having mutual affinity to each other leading to a (a) formation of an affinity complex that is immobilized to a solid phase in a detection microcavity or (b) one or more other reaction products that may be soluble or insoluble in the detection microcavity. By properly selecting the reaction conditions including selection of reactants, it can often be arranged so that the product obtained and/or a reagent in excess are detectable with a signal that can be (a) measured from the above-mentioned detection areas and (b) related to one or more features of a starting liquid aliquot introduced into a microchannel structure. Typical such features are kind, form and amount including activity. etc., of a particular reactant (e.g. an analyte) including for instance an affinity reactant such as an enzyme etc. The term "can be related to one or more features" includes also the determination of the manner in which reaction variables such as pH. ionic strength, detergents, etc., might influence the reaction used for forming the reaction product. Typically, one makes use of detection principles based on radioactivity, fluorescence, chemiluminescence, bioluminescence, enzymatic activity, chromogens, light scattering (turbidometry), etc., for instance by utilizing a reactant that carries a group providing detectability, either by being detectable as such or by being transformable to a detectable group. Typically the utilized process protocol means that a detectable reactant is incorporated into a complex or into some other reaction product. See for instance WO 02075312 (Gyros AB). Detectable products, reagents, etc , that can be formed and/or retained and measured iq a microcavity via a detection area are collectively called "substance" in other parts of this specification.
Typical reactants in this context include members of affinity pairs such as (a) antigen/hap ten and the corresponding antibody including its antibody active fragments, (b) lectin and the corresponding carbohydrate structure, (c) native ligands and the corresponding receptors.
(d) complementary nucleic acids including synthetic variants such as synthetic oligonucleotides. (e) Ig(Fc)-binding proteins and Protein A, Protein G and other Ig(Fc)-receptors, (e) ion pairs of opposite charges, enzyme and the substrate, inhibitor, cofactor. coenzyme etc., that can bind to the enzyme, etc. Synthetic variants more or less mimicking a native affinity interaction are also included.
One aspect of the invention is a method and corresponding arrangement, for determining the amount of a substance in a detection microcavity of a microfluidic device and comprises collecting radiation associated with the substance from a detection area associated with the detection microcavity. The method is characterized in comprising the steps of:
a) providing
(i) a microfluidic device, e.g., in the form of a disc, comprising
A) a plurality of microchannel structures, each of which has an inlet port, a detection microcavity and a microconduit connecting the inlet port with the detection microcavity, and
B) a plurality of detection areas and associated search areas as discussed elsewhere in this specification, each of the detection areas being (1 ) associated with one of said detection microcavities, (2) present in the surface of said device and (3) translucent/transparent for said radiation, and
ii) a detector arrangement which is capable of collecting radiation from individual subareas of each of said search and detection areas;
b) processing one or more liquid aliquots in at least one of said plurality of microchannel structures so that the substance is retained in the detection microcavity of each of said at least one of said plurality of microchannel structures, provided that said substance and/or one or more reagents that are necessary for the substance to be retained in a detection microcavity are present in at least one of said one or more aliquots;
c) scanning the search and/or detection areas associated with the detection microcavities that are part of microchannel structures in which step (b) has been carried out to obtain radiation from individual subareas of each scanned search or detection area, said scanning being performed by the use of said detector arrangement;
d) integrating radiation as a function of the subareas of each scanned search area or detection area to obtain the amount of radiation from each search area or detection area;
e) characterizing for each of the amounts obtained in step (d) a reaction variable that has been included in the process protocol used for each microchannel structure.
Further details about characterization of reaction variables are discussed above and in WO 02075312 (Gyros AB). The characterization includes e.g., that the amount of the substance in absolute or relative terms in each of the detection microcavities is determined from each of the amounts obtained in step d). The substance in this context may be the substance from which the radiation is directly derived or a substance, such as an analyte, the amount of which in one of the liquid aliquots used in step (b) is related to the amount of radiation determined in step (d). Since it is the result of the reaction that is measured the analyte may or may not be present in a detection microcavity after step (b).
Step (c) comprises collection of radiation intensities which refer to changes in the radiation due to the amount of a particular substance in the detection microcavity. This includes true intensity values but also other changes, such as changes in refraction, in scattering, in the plane of polarization and life time of fluorescent light (including half life, e.g.), etc.
The microfluidic disc and steps b) and c) are illustrated elsewhere in this specification. In addition to circular scanning, e.g., by rotating a disc, step (c) also comprises non-circular scanning for instance scanning by lateral movement of a linear detector head comprising one or more rows of detector elements over the detection areas. Also imaging by a CCD camera is included. The microfluidic disc may be circular or have some other geometric form, for instance triangular, rectangular, etc., including also irregular forms.
Step (d) means integration over each detection area, i.e., primarily over subareas which have radiation values that deviate from the values obtained for surface parts of the device that surrounds the detection area. Alternatively one may exclude all or selected parts of a detection area which corresponds to parts of a detector microcavity in which the presence of the substance is insignificant. The integrating step includes the substeps of (a) finding a start and/or stop position at edges of the detection area, e.g., at an end corresponding to the inlet end of the detection microcavity, and (b the factual integrating. In one variant substep (a) is carried out by determining the inflection point for the amount of radiation per subarea versus position along the detection area. Due the fact that the delineating part of the detection area may be curved, at the upstream end, the invention also suggests that the integration should account for curvatures in the circumference of a detection area. In a preferred variant, substep (a) comprises determining a threshold that segments detection area pixels from the background, e.g., if the detection microcavity is filled with a particle bed a threshold that segments the particle bed pixels from background pixels. This can be done with optimal thresholding or determining median or mean background or any other way of determining background (for "optimal thresholding" see for instance reference [1]: "Digital Image Processing", 2nd edition. Editors: Gonzales R C et al, page 354). In this case substep (b) will mean integrating a selection of pixels, i.e., those pixels which have radiation values above the threshold and belong to the main group (detection area) and excluding noise pixels above threshold that do not belong to the main group of connected pixels. The integrating typically starts from pixels corresponding to one end of the detection area, e.g., the inlet end.
Between the preferred variant of substeps (a) and (b), there are preferably additional substeps for refining the method, such as
i) Creating a binary image from the calculated threshold.
ii) Labelling the high binary pixels into different groups (= labeling the image). Each group will consist of pixels that are bordering to each other (close to each other). The binary high pixels from the detection area will define the main group.
iii) Optionally removing all binary high pixels that do not belong to the main group.
The factual integrating (main step d) will in this case mean integrating the radiation values for the binary high pixels of the main group.
The inventors have found that when working with very small volumes and amounts, the material from which a device is manufactured as well as the pretreatment procedures may introduce radiation artifacts in form of peak noise that is comparable to the radiation coming from subareas outside the peak. Therefore step (d) may also comprise substeps for removing peak noise. This typically means a first substep in which the deviating radiation (peaks) are made more apparent. One way of doing this is la Place filtering, point detection and the like. See for instance see reference [1]: "Digital Image Processing", 2nd edition, Editors: Gonzales R C et al, pages 333, 334 and 339 (1987). In a subsequent substep the width and the position of each peak is calculated for instance by including edge detection or edge linking based on a local area etc.. See for instance reference [1]: "Digital Image Processing". 2nd edition, Editors: Gonzales R
C et al, pages 334 and 344 (1987). In the next substep the peak noise is removed by interpolating from surrounding subareas. An alternative way to the whole process of removing peaks is morphological opening, un-linear filtering operating in local histogram domains etc. These substeps for removing peak noise is performed prior to the factual integration.
Step (e) is conventional and typically includes that the integrated value is compared with the value for one or more standards. A standard value is typically the integrated value for a known amount of a standard substance, which in most cases is the same as the substance under investigation.
The scanning step (c) and integrating step (d) may in certain innovative variants of one aspect of the invention be optional. Thus, these steps are preferably included when the desired substance is unevenly distributed within the detection microcavity. In a typically case, this may happen when the substance is retained within the detection microcavity under flow conditions, for instance by being captured to a solid phase introduced into the microcavity prior to the fluorescent substance.
In the that case the substance is homogenously distributed within the detection microcavity, the scanning step (c) and integrating step (d) may be replaced with collecting the radiation intensities for selected subareas of a detection area and letting these intensities represent the total amount of radiation from the detection area, for instance as a mean or maximum value. Step (e) can then be carried out on these values in the same manner as for values obtained by scanning and integrating. This way of performing the method is also a part of the present invention. *
Homogenous distribution of the substance in the detection microcavity typically is at hand in case the substance is present in the detection microcavity in an equilibrated solution and/or when a reaction is going in the microcavity between homogeneously distributed reactants either forming or producing the substance that is to be detected. The innovative method of the aspect discussed in the preceding paragraphs includes that steps (d) and/or (e) are performed in close connection to the preceding steps [(a)-(c) or (a)-
(d), respectively] and/or that the data from the scanning and/or the integrating have been obtained at an earlier time, for instance at a different geographical location, and/or by different individuals.
In the described method for determining an amount of a substance in a detection microcavity of a microfluidic device, it may be necessary to reduce the noise contribution in the detected radiation intensity from scanning step (c) in the integration step (d).
A preferred embodiment of the invention is now described.
Method for reducing noise in a substance raw data image containing noise.
An algorithm, or rather an automated method, for that is detectable by use of a particle bed is presented hereinafter. This algorithm may be started for each scanned search area/detection area/particle bed of interest.
The invention will now be further detailed with respect to a particular example comprising a microfluidic device (see figure 3) in a disc format, comparable to Compact Disc (CD), having an axis of symmetry perpendicular to the disc plane and with a radial flow direction from the axis of symmetry for each of the detection microcavities (see figure 3: reference number 306 a,b,c ... ) that comprises a bed of particles. This kind of microfluidic device will henceforth be called CD. The signal detected and measured may be fluorescence from an antibody fluorochrom conjugate (antibody-fluorochrom). Each bed defines a detection area (see figure 3 : reference number 317 a,b,c ... ) in the surface of the disc. See for instance WO 02075312 (Gyros AB). The detection area is part of a search area comprising also a defined minor zone around the detection area. The principles outlined in this example may be adapted also to other devices. The principles are also applicable to other kinds of radiation as discussed elsewhere in this text.
Generation of a substance raw data image
As described in section III, Controller, of this description, a controller means (fig. 1 : ref. No. 1 12) is programmed to control the scanning and detection. Each search area is scanned by a detector head that is capable of detecting a detectable physical feature deriving from a substance present in the microcavity. like radiation caused by radioactivity, fluorescence, chemiluminescene, bioluminescence, enzymatic activity, scattering of light, absorption of light, reflection of light etc. The measured data are stored in a data storage, an image data storage, connected to or connectable to an processing unit, which may either be comprised in the controller means, or as in another possible embodiment, incorporated in a separate computer system. Each search area is identified by an identifier, which may be used as address. The identifier makes it possible for the system and controller means to locate each single search area and its detection area on a CD and link the corresponding stored measured data in the image storage.
Now with reference to figure 4. each search area 402 has a pre-defined radial start position, Xstart, a pre-defined angular start position. Ystart. a pre-defined angular stop position. Ystop, and a pre-defined radial stop position, Xstop. The start of the collecting of data is started at a position, primarily an angular and/or a radial position that is prior to an intended detection area, and to end the collecting at a position that is after the same detection area. The search area will therefore cover an outer area 406 surrounding the detection area 404. The search area 402 is scanned in radial direction by the detector disc lap by disc lap. The radial start position Xstart is the number of disc laps/revolutions, counted from a spinning start position, Spinstart (not shown). When the scanning of a search area starts, and the radial stop position is the number of laps, counted from the spinning start position, when the last angular scanning of the search area is to be performed. The scanning of the search area in the angular direction starts during the lap, which is defined by the radial start position, at the angular start position in the angular direction, which is defined by the spinning direction of the disc and the detector position in relation to a Home Mark on the CD. The scanning in angular direction is stopped at the angular stop position.
Figure 4 is schematic in the sense that the search area of a rotatable device typically is a) slightly narrowing when going from an outer position inwards towards the centre of the CD, and b) slightly arc-shaped. Compare for instance that the Ystart position for each lap are always at the same angular position (which also applies to the Ystop position).
During a lap, from the angular start position to the angular stop position, a stream of collected data from that search area will generate a corresponding stream of measured data that are stored into the storage. This is repeated lap by lap from the radial start value to the radial stop value. The processing unit of the controller 120 links measured data to the correct identified detection area, and there will be no mixing of measured data from different search and/or detection areas. Each scanned search area has a corresponding set of collected measured data.
The processing unit will organize the data in such an order, position by position in each set of measured data, and it is therefore possible for the image processor to locate a certain measured value in the set. The set of measured data has therefore similarities to an electronic and/or digitally stored image. An image is defined in very small elements, called pixels. Each pixel corresponds to a position of a measurement area in the search area.
When the scanning, detecting and storing are finished, the measured data of a search area is possible to plot by use of printer/plotter software program used by the processing unit in a three dimensional diagram, wherein the horizontal plane is defined by the two dimensions x and y, in this case the radial direction, x-axis, and the angular direction, y-axis, which is perpendicular to the x-axis. Positive x-direction is defined as the main flow direction (downstream) through the detection area, from the inlet 408 of the detection area 404 to the outlet of the same. Positive y-direction is defined as the scanning direction which is orthogonal/perpendicular to the radius (or flow direction). The third dimension, z-axis. perpendicular to the horizontal plane defined by the x- and y-axis, is the magnitude of the measured value, e.g. measured value of the fluorescence intensity. The measured data defines a surface of measured data, having peaks and valleys. As mentioned above, the three dimensional figure has essential similarities with an image having an intensity value defined for each pixel, and the surface of measured data is therefore called an image of a search area/detection area/particle bed comprising the substance from which the measured signal derives. In the description of this application, an image is the stored detector measurement data belonging to a search area/detection area/particle bed. A pixel is here defined as a discrete measurement area within the search area/detection area and each pixel may have a measured intensity value, a pixel value. The resolution of the image is depending of the number of discrete measurement points in the y-direction and the number of feps in the x-direction within the search area. The resolution increases with increasing density of measure points and scanning laps, and consequently, the area of a pixel will decrease. If the electromagnetic-emitting source is a laser, the possible minimum pixel area is limited by the minimum cross-section of the laser beam with the search area. During the scanning of the search areas the pixel intensities are determined by the detection process and stored as digital pixel data representing raw data images of the search areas.
The space between the surface of measured data for the search area/detection area and the base plane, defined by the two axes in the base plane, constitutes a volume. The whole or a part of the volume may represent a wanted measure for a chemical process, i.e. a reaction variable. This value is possible to calculate by using a mathematical method, e. g. addition or integration of measurement data over the whole area of the detection area, or one sub-area 410 (fig. 4) or a number of sub-areas of the detection area. Sometimes the object is to calculate an amount for an subarea, here called a normalized area, i.e. a predefined amount of pixels, inside the detection area. This is done to minimize the variance between detected structures. The variance depends among other things on background radiation, such as fluorescence, from the detection area.
In an ideal image of a detection surface, the area is smooth, because of the absence of noisy background contribution in the measured raw data signals from the detector. This image includes only measured data from a substance signals, i.e. the substance contribution to the detected and measured raw data signal. In reality, a substance raw data image comprises measured substance data and noise data. The noise cause irregularities, like noise peaks, in the surface, and evidently, a calculation of the volume in an area having one or many noise peaks will not result in a true measure for a chemical process.
It is therefore necessary to reduce the noise contribution in a substance raw data image as much as possible before calculation of the measure.
The present invented method provides a method for reducing noise contribution in a substance raw data image derived from a detection area of a microfluidic device of the kind described above. The method comprises a number of main-steps for reducing different kind of disturbances and noise within an ima 'gtoe1-.
The method removes noise from substance raw data information in substance raw data images. Each image is associated to one detection microcavity and one substance that is present in the detection microcavity of a microfluidic structure of a microfluidic device. The substance raw data information is obtained by use of a detector arrangement capable of scanning a search area associated to a detection area, which is associated to the detection microcavity. The detector arrangement is capable of detecting at least one of a number detectable physical features discussed above, e.g. fluorescence, luminescence, etc. emitting radiation from the search area, and generate output that is proportional to the magnitude of detectable activity. The output is converted to measure data storable as substance raw data in a set in a data storage. The method may comprise following steps in any order for removing different noise signals contributions step-by- step in any order in a substance raw data image:
- reducing background radiation (step α);
- reducing peak disturbance (step β); - locating the detection area within the search area (step χ);
- moving binary artefacts (step δ);
- removing unwanted areas of the detection area (step ε);
- applying default detection area in noisy images (step φ).
Each of the steps α - φ comprises a number of sub-steps. Steps α - φ may be regarded as an algorithm, but is not necessary that all of the steps α - φ are carried out for all images. The order and kind of steps is typically preselected and depends among others of the sensitivity and accuracy required in each individual step. In general, inclusion of more steps is required when high sensitivity is desired.
The different steps will now be described in more detail.
Step α: - Reducing background radiation and construction of background processed images.
The object is to remove the electronic, optical and dust noise that may be visible in the image. This part of the invented method mainly comprises two variants.
In the first variant, a background radiation image is used. In figure 5, the first variant is presented schematically in a flowchart. Following steps may be used:
- if available, accessing a background radiation image or. if not available, generating a background radiation image (step α' 1); - if available, accessing corresponding substance raw data image or, if not available, generating substance raw data image (step α'2):
correlating the two images and generating a correlation image (step α'3);
- locate the positions of the largest pixel values (step α'4);
- positioning the images correctly, subtracting corresponding pixel values and storing the generated background processed image(step α'5).
This first variant is now discussed in more detail.
In the first variant, the first step, step ' 1, is to scan the disc and measure the radiation, e.g. the intensity of the fluorescence, without presence of the substance giving raise to the measured signal from the detection areas, e.g. before step b) in the method for determining the amount of a substance, retaining the radiation associated substance before the detection areas/particle beds of the disc have been treated, e.g. with antibody-fluorochrom. All search areas/detection areas are scanned and each of the corresponding sets of measured data is stored in a background radiation image. When using antibody-fluorochrom, the corresponding image is called background fluorescence image. The image contains the measured data of the background radiation noise signal and other noise signals. In preferred variants the scanning for the background radiation image is carried out as close as possibly to the step in which the substance causing the radiation is formed and/or retained in the detection microcavities. Although not preferred, the manufacturer may provide a background radiation image that for instance may be based on the average of many CDs of the same kind.
The second step, step α'2, is in this first variant the same as step c) in the method for determining the amount of a substance, above. The substance is now present in the detection microcavities, for instance retained in particle beds that may have been washed, e.g. after the particle beds have been treated with e.g. antibody-fluorochrom. This image is the substance raw data image, described above.
If the background radiation images and corresponding substance raw data images exist already and are availably stored, there is consequently only to access (e.g. by copying or reading) the images in step α' 1 and step α'2 instead of generating them again by the scanning and detecting processes.
The two images should be detected under the same conditions, e,g. same laser power and same detector sensitivity if the detection principles is based on fluorescence.
The third step of the first variant, step α'3, is to correlate the two images in two dimensions, that define a correlation image, by use of a two-dimensional correlation. The use of two-dimensional correlation is well-known in the art of signal-processing, see reference [1], page 90-92. Thus, the step has resulted in a correlation image.
In the fourth step, step α'4, of the variant, the position of the pixel, or pixels, having the largest intensity value in the correlation image is/are located. This could be done by using a simple comparison algorithm, wherein the recently largest value is stored and compared with the next data, e.g.intensity value, that hasn't been compared with any of the earlier pixel data. The positions of the pixels having the largest pixel values, largest intensity values, are stored and used in the next step, step α'5.
In step α'5, the background radiation image is moved to the correct position of the correlation image by using the identified pixel positions. The parts of the image that overlap is subtracted pixel by pixel.
These five steps have resulted in a new image, the background processed image.
The substance raw data image comprise the wanted substance radiation data. If not any new noise have been introduced by the introduction of the substance giving rise to the measured signal into the microcavities (part of step b). the process of noise removing may be stopped.
In the second variant-is a Median value of the background radiation data used. In figure 6, the first variant is presented schematically in a flowchart. Following steps may be used:
- accessing or generating a substance raw data image (step α"' 1);
determining a Median value for background radiation data (step α' "2); - subtracting pixel values of substance raw data image with the median value, and store background processed image (step α"3).
The second variant of reducing background radiation and construction of background processed images will now be described in more detail.
This second variant for removing background radiation, like fluorescence, especially for removing electronic and optical noise that exists in the search area outside the detection area, is possible to use instead of the five steps described above. This alternative variant works very well under the condition that there is no unwanted object of large size outside the detection area and the background signal intensity is constant inside the image. This variant has the advanta *gfc e that the disc is scanned onlv once, which is described hereinafter.
The first step, step α" 1 in this alternative variant, should start after the radiation associated substance has been incorporated into the detection microcavities. In the method for determining the amount of a substance above this for instance corresponds to after the detection areas have been treated with the antibody-fluorochrom, i.e. in step c. As described above, during the scanning of the search areas for the radiation, the pixel values of the substance raw data images are detected and determined, and stored as digital pixel data. If the substance raw data images exist already and are availably stored, there is consequently only to copy/read the images, instead of generating them again by scanning and detecting processes.
In the second step, step α"2. a median value of the pixel values for at least a portion of the search area outside the detection area is calculated and used as background radiation value, e.g. as a value for the background fluorescence. The median value is an estimated value for the background radiation (intensity).
In step α"3, the substance raw data is decreased with the median value contribution, e.g. by subtraction, pixel pixel, until all pixels are processed and a background processed image is provided.
Potentially the step for determining Median value ma)' be replaced by determining some other form of an average value and use this value instead of the Median in the subsequent steps. Both methods, α' and α", provide a background processed image for each scanned search area, i.e. an image with reduced background radiation, such as fluorescence, and can be used separately or in combination in other steps of the innovative method, with the former being preferred at the moment.
Step β: Reducing peak disturbances.
The method of this aspect of the invention may comprise one or more steps in which peak disturbances originating from dust, chemistry and other artefacts are reduced, e.g. reducing peaks from fibres. The disturbance peaks have typically a relatively high amplitude and a relatively narrow distribution. However, the image of an ideal detection area is smooth and does not contain any disturbance peaks.
In figure 7, step β for reducing the peak disturbances is presented schematically in a flowchart. Following steps may be used:
- Generating a Laplacian processed image from the earlier stored image, typically a background processed image, by using a Laplacian filter in two dimensions (step βl);
Calculating a disturbance threshold value from the image created in the previous step, typically from a background processed image (step β2);
- Generating binary encoding Laplacian processed image by applying the disturbance threshold limit to the image used for calculating the treshold value (step β3);
- Finding the start points in clusters (step β4);
Determining a minimum rectangular area surrounding the start points (step β5);
- Calculating an interpolated value within each peak area and substituting the original data of the peaks, and storing the disturbance filtered image (step β6). The process for reducing peak disturbances will now be described in more detail.
Thus in a first step, step β 1 , a radiation intensity image, for instance a background processed image as discussed above, may be processed by a Laplacian filter in two dimensions. This is a standard method earlier known from reference [1]. The result of this step is that peaks in a background processed image are amplified. The result is normalized with the maximum radiation (intensity) in the background image. This step provides a Laplacian processed image.
In a second step, step β2, a disturbance threshold limit is calculated by use of the standard statistic functions average and standard deviation of the stochastic background noise, calculated from the image used in step β 1 . This step enhances the disturbance peaks of the image.
In the third step, step β3, the Laplacian processed image is binary encoded. In the encoding process pixel intensities/values equal to or higher than a threshold limit are set to one of the two corresponding binary values, for instance "1 " or "positive" and pixel intensities less than the disturbance threshold limit are set to the other of the binary values, for instance "0" or "negative". The determined disturbance threshold limit determined in step β2 may be used in step β3 for encoding. This step creates an image wherein the high intensity pixels of the peaks appears as groups, even called clusters.
In a fourth step, step β4, the Laplacian processed image may then be processed for finding the pixel in each peak having the largest value by using the position information from the clusters in the binary Laplacian processed image. Each pixel having the largest value in a cluster is defined to be a center and start point. The other pixels of the clusters are set low in the binary Laplacian processed image. This step results in a cluster image.
In a fifth step, step β$, the start points defined in the cluster image is used for determining a minimum rectangular area by iterative testing using different area sizes. The iteration, that works on the Laplacian processed image, is stopped when a minimum rectangular area contains only pixels having a negative intensity value or zero value according to the definition in previous steps and as a result of the earlier performed Laplacian filtering. Said pixels surround the corresponding start point and defines the position, the area of the corresponding disturbance peak and the boarder pixels which is outside the peak area.
In a sixth step, step β6, the disturbance peaks of the background processed image are removed by substituting the original data within each peak area with an interpolated value between the boarder pixels of the corresponding peak. The interpolated values are calculated by use of an appropriate interpolation function, both a linear and a non-linear function is a possible choice.
The process has now generated a "disturbance filtered image" of the radiation intensity, e.g. fluorescence intensity, of the detection area.
Step χ: Determining a global threshold value
The object of this part of the invention is to make a radiation intensity image, such as a disturbance filtered image, binary by determining a global threshold that separates the measured substance signals, i.e. the signals from the substance in the detection area, from the background signals. For instance, the detection area pixels should be equal to binary high, like "1", and the rest of the image equal to binary low, like "0".
The purpose of the following steps is to determine a valid global threshold value. This is done by an iteration that tests different threshold values. A global threshold is valid when only a few pixels in a part of the search area that is not connected to the detection area have an intensity higher than the pixels in the detection area.
In figure 8, step χ for determining a global threshold value is presented schematically in a flowchart. Following steps may be used:
Generating a y-sujnmation graph (step χl);
Median filtering of the y-summation graph (step χ2);
- Determining centre position of detection area (step χ3);
Determining radial directed edges of detection area (step χ4); Generating a x-summation graph (step χ5);
- Median filtering of the x-summation graph (step χ6);
Determining an edge curve (correlation) (step χ7);
Restraining disturbance peaks in the edge curve (step χ8);
- Determining the start and end positions of the detection area (step χ9);
- Determining a start value of the global threshold value (step χlO);
Adjusting the global threshold value to the value giving correct detection area width (step χl 1);
Generating a binary encoded disturbance filtered image by use of the global threshold value (step χl 2).
This process will now be described in more detail.
In a first step, step χl , pixel values are summarized in the radius direction, pixel values, line by line, in a radiation intensity image, typically the disturbance filtered image, to get a one-dimensional y-summation graph in the y-direction of the radiation intensity, e.g. fluorescence intensity (angular direction on a spinning disc). Each y-coordinate on the graph will be the sum of the pixel values in radial direction having the same y-coordinat.
In a second step, step χ2,the y-summation graph is filtered to remove disturbing peaks, e.g. with a median filter function.
In a third step, step χ3. the centre of the detection area in y-direction (angular direction on a spinning disc) is calculated and located by use of a correlation, and with the filtered y-summation graph and the known detection area width as input value to the correlation length. This should normally give a graph with one large peak.
In a fourth step, step χ4, the centre of the largest peak is set to the centre of the detection area. It's now possible to calculate the positions and the radial extension of the particle bed edges by use of the known data for the centre and detection area/particle bed width. Hereby, the radial directed boarder/edge pixels are known.
In a fifth step, step χ5, the radiation intensities. are summarized in the image in y- direction (angular direction on a spinning disc) using the limits (boarders) determined in the previous step to get a x-summation graph (radial summation graph) in x-direction (radius direction) of the radiation intensities. The summation of pixel values is performed in the same way, but in the orthogonal direction in the same plane, as in step χl .
In the next step, step χ6, the x-summation graph is filtered with a median filter function to remove disturbing peaks.
In the seventh step, step χ7, the start of the detection area in the x-direction is determined by use of correlation between a step function and the filtered x-summation graph. The mathematical process results in an edge curve, see figure 9.
In step χ8, disturbance peaks, close to edge that is sitting opposite to the inlet, are restrained by multiplying the edge curve with a ramp function.
The purpose of step χ9. is to find the start position of the detection area by finding the highest positive peak centre in the filtered edge curve (from step χ8) and using a minimum detection area length to calculate the end position of the detection area.
An approximate position of the detection area edges, boarder pixel is now known.
In the tenth step, step χlO, the global threshold value is now determined by use of optimal threshold technique for image processing (reference [1], starting at page 360). It is assumed that the background and particle bed peaks are normal distributed with the same variance. Calculate from the disturbance filtered image the estimated background intensity as max of (median background intensity on the areas A1+A3 or median intensity background on the area A2)). The positions and orientations of the areas Al . A2 and A3 is illustrated in figure 13. Calculate the estimated substance intensity as median of the detection area. Calculate the start value for the iteration step as average of the estimated background intensity and the estimated substance intensity. In step χl 1, the global threshold value is iterative adjusted, between a maximum intensity value, defined by the estimated substance intensity, and a min value, defined by the estimated background intensity, so that the average number of pixels in the y-direction (the angular direction) inside the approximate position of the detection area is equal to the detection area width.
In step χl2, the calculated global threshold value is used for binary encoding a stored image, in particular the stored disturbance filtered image into a binary coded disturbance filtered image.
Step δ: Remove binary objects outside the detection area
The object of the following steps δl - δ2, is to remove any disturbances outside the detection area. The following algorithm processes the data of the binary coded disturbance filtered image, in which all objects are set to 1, i.e. binary high. Further, it is assumed that there is no object to be removed inside the boarders of the approximately calculated detection area position, and the detection area is the largest area of the image.
In figure 10, step δ for removing binary objects outside the detection area is presented schematically in a flowchart.
The first step, step δl : labelling of all the objects in the binary coded disturbance filtered image, wherein the number of connected pixels, i.e. the area, of each object is determined. Labelling is a known method in image processing.
The second step, step δ2: generating a binary detection area image by removing all but the largest objects that are inside the approximate detection area position.
Step ε: Removing unwanted areas of the detection area
Sometimes the object is to calculate an amount for a ""normalized" area. i.e. a predefined amount of pixels, inside the detection area. This is done to minimize the variance between detection areas on the same microfluidic device. The variance depends among other things on background radiation, such as fluorescence, from the detection area. It is assumed that the 3D distribution of an intensity signal inside the detection area is independent of small variations in the detection area length. The substance signal deriving from a portion of the detection area that may be closests to the inlet of the corresponding detection microcavity. The algorithm used is a process for generating a new binary image that only comprises pixels within the normalized area ((detection area width) multiplied by (min detection area length)) in the binary detection area image starting from the inner part of the CD.
In figure 11, step ε for removing unwanted areas of the detection area is presented schematically in a flowchart.
The first step, step εl : generating a normalized area image by setting all binary high indicated pixels outside the location of the normalized area ((detection area width) multiplied by (min detection area length)) to binary "low" in the binary detection area image.
The second step, step ε2: multiplying each corresponding pixel in the disturbance filtered image, generated and stored earlier in the invented process, and the normalized area image with each other and perform a calculation, e.g. integration, on the normalized area.
Step φ: Apply default detection area in noisy images
The described method above will not work properly, or not at all, for very noisy images hiding the substance signal. The following algorithm works, if all processed images originate from the same set of detection areas, e.g. the same CD. The detection areas are then well-defined in angular position and quite well in radius direction.
During the algorithm to locate the position of the detection area within the search area, all positions for detectio areas having high confidence are stored.
In figure 12, step φ for applying a default detection area in noisy images is presented schematically in a flowchart.
The method for applying default detection area in noisy images may comprise following steps: Step φl : Saving automatically all determined detection area positions (step χ) fulfilling high confidence criterion.
Step φ2: Calculating an average detection area position.
Step φ3: Using automatically the average detection area position if one or more predetermined condition is fulfilled.
The average position for a detection area will automatically be applied for all calculated detection areas that fulfil all following conditions:
• It's angular start position differ from the average start position more than [amount] multiplied by [detection area width], wherein the [amount] is chosen 0-1.
• The average number of binary high pixels, in the normalized area image. in angular direction for all radiuses is less than [amount] multiplied by [detection area width], wherein the [amount] is chosen 0-1.
• The difference between angular stop and start positions is greater than [amount] multiplied by [detection area width], wherein the [amount] is chosen 1-10.
• The number of binary high pixels, in the normalized area image, in angular direction for all radiuses is less than [amount] multiplied by [detection area width] multiplied by [min detection area length], wherein the [amount] is chosen 0-1.
Steps d) and e) in the method for determining the amount of a substance that is associated with the collected radiation as well as any of the other individual method steps or their combinations discussed are typically performed by the appropriate software for instance included in the controller or elsewhefe, for instance not in physical association with the innovative arrangements described herein.
The invention also comprises computer program-related aspects for treating radiation data that have been assembled by steps (a)-(c) of the innovative method. One such aspect is a computer program product that (1) comprises program code elements corresponding to a sequence variant of step (d) above comprising one or more of the substeps described for this step, and
(2) when installed on the appropriate hardware is capable of causing the hardware (computer) to execute the sequence of substeps on data which have been obtained by performing step (c) above on any kind of microfluidic device irrespective of being spinnable or not.
Other code elements may be included, for instance corresponding to step (e) in order to execute the sequence (d)-(e). Another computer program-related aspect is the computer program product stored at a computer program readable means which, when the product is loaded, makes it possible for a computer to perform the sequence of steps corresponding to the code elements in the stored computer program product. A third computer program- related aspect is a carrier having at least one of the innovative computer program products thereon. The carrier may be a computer memory, a Read-Only Memory or an electrical signal carrier.
The present invention may be implemented as a computer program product directly loadable into an internal memory storage of a processing unit within the controller in the detector arrangement, comprising the software code means for performing the steps of any claims 1-16.
Further, the present invention relates to a computer program product stored on a computer usable medium, comprising readable program for causing a processing unit in a computer means to control an execution of the steps of any of the claims 1-16.
The computer usable medium is a record medium, a computer memory, a Read-Only Memory or an electrical carrier signal.
Certain innovative aspects of the invention is defined in more detail in the appending claims. Although the present invention and its advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the invention as defined by the appended claims. Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the disclosure of the present invention, processes, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present invention. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps.

Claims

1. Method for reducing noise in an image, typically a substance raw data image, comprising substance and noise data which is associated to one detection microcavity and one substance is present in the detection microcavity of a microchannel e structure of a microfluidic device, the substance raw data information obtained by use of a detector arrangement capable of scanning and collecting radiation data in a search area associated to a detection area, the collected radiation data are stored as measure data for the amount of the substance as substance raw data in a set in a data storage, said method is characterized in that the method comprises following step: - reducing background radiation data either by use of a image comprising background radiation data and no substance data information, a background radiation image, (step α') or, by use of a calculated median value of the background radiation data (step α").
2. Method according to claim 1, characterized in that the method comprises the step of: - if available, accessing a background radiation image or, if not available, generating a background radiation image (step α' l);
- if available, accessing corresponding substance raw data image or, if not available, generating a substance raw data image (step α'2);
- correlating the two images and generating a correlation image (step α'3); - locate the positions of the largest pixel values (step α'4);
- positioning the images correctly, subtracting corresponding pixel values and storing the generated background processed image(step α'5).
3. Method according to claim 1, characterized in that the method comprises the step of:
- if available, accessing corresponding substance raw data image or, if not available, generating a background radiation image substance raw data image (step α" 1);
- determining a Median value for background radiation (step α"2);
- subtracting pixel values of substance raw data image with the median value contribution, and store background processed image (step α! ,3).
4. Method according to claim 1, characterized in that the method comprises the step of: Reducing peak disturbances (Step β).
5. Method according to claim 4, characterized in that the step of reducing peak disturbances (Step β) comprises following steps:
- Generating a Laplacian processed image from the earlier stored image, typically a background processed image, by using a Laplacian filter in two dimensions (step βl); - Calculating a disturbance threshold value from the earlier stored image, typically a background processed image (step β2);
- Generating binary encoding Laplacian processed image by use of the disturbance threshold limit (step β3);
- Finding the start points in clusters (step β4); - Determining a minimum rectangular area surrounding the start points (step β5);
- Calculating an interpolated value within each peak area and substituting the original data of the peaks, and storing the disturbance filtered image (step β6).
6. Method according to claim 1. characterized in that the method comprises the step of: Determining a global threshold value (step χ).
7. Method according to claim 6, characterized in that the step of determining a global threshold value (step χ) comprises following steps:
- Generating a y-summation graph (step χl);
- Median filtering of the y-summation graph (step χ2);
- Determining centre position of detection area (step χ3); - Determining radial directed edges of detection area (step χ4);
- Generating a x-summation graph (step χ5);
- Median filtering of the x-summation graph (step χ6);
- Determining an edge curve (correlation) (step χ7);
- Restraining disturbance peaks in the edge curve (step χ8); - Determining the start and end positions of the detection area (step χ9);
- Determining a start value of the global threshold value (step χlO);
- Adjusting the global threshold value to the value giving correct detection area width (step χl 1);
- Generating a binary encoded disturbance filtered image by use of the global threshold value(step χ 12).
8. Method according to claim 1, characterized in that the method comprises the step of: removing binary objects outside the detection area (step δ).
9. Method according to claim 8, characterized in that the step of removing binary objects outside the detection area (step δ) comprises following steps: - labelling of all the objects in the binary coded disturbance filtered image, wherein the number of connected pixels of each object is determined (step δl);
- generating a binary detection area image by removing all but the largest objects that are inside the approximate detection area position (step δ2).
10. Method according to claim 1, characterized in that the method comprises the step of: removing unwanted areas of the detection area (step ε).
11. Method according to claim 10, characterized in that the step of removing unwanted areas of the detection area (step ε) comprises following steps:
- generating a normalized area image by setting all binary high indicated pixels outside the location of the normalized area ((detection area width) multiplied by (min detection area length)) to binary '"low" in the binary detection area image (step εl);
- multiplying each corresponding pixel in the disturbance filtered image, generated and stored earlier in the invented process, and the normalized area image with each " other and perform the calculation, e.g. integration, on the normalized area (step ε2).
12. Method according to claim 1. characterized in that the method comprises the step of: Applying default detection area in noisy images (step φ).
13. Method according to claim 12, characterized in that the step of applying default detection area in noisy images comprises following steps:
- Saving automatically all detection area positions fulfilling high confidence criterion (step φ 1 ); ; - Calculating an average detection area position (step φ2);
- Using automatically the average detection area position if one or more predetermined condition is fulfilled (step φ3).
14. A method for determining a measure of at least one substance comprising the steps of: (a) providing a microfluidic device and a detector arrangement, wherein said microfluidic device comprises a plurality of microchannel structures, each of which has an inlet port and a detection microcavity, and a plurality of detection areas, each of which is associated with one of said detection microcavities; and wherein said detector arrangement is capable of collecting radiation from individual subareas of each of said detection areas;
(b) processing one or more liquid aliquots in at least one of said plurality of microchannel structures so that said substance is retained in the detection microcavity of each of said at least one of said plurality of microchannel structures;
(c) scanning the detection areas associated with the detection microcavities that are part of microchannel structures in which step (b) has been carried out to obtain radiation from individual subareas (pixels) of each scanned detection area, said scanning being performed by the use of said detector arrangement;
(d) integrating radiation as a function of subarea for each scanned detection area to obtain the amount of radiation from each detection area; and
(e) characterizing for each of the amounts obtained in step (d) a reaction variable that has been included in the process protocol used for each microchannel structure.
15. The method of claim 14, wherein at least one of the steps in claims 1-13 is carried out in step (d).
16. The method of claim 15, wherein integration only is on a portion of each detection area, a normalized area.
17. A computer program product directly loadable into an internal memory storage of a processing unit within the computer means, such as a controller of a detector arrangement, comprising the software code means for performing the steps of any claims 1-13 and/or 14-16.
18. The computer program product stored on a computer usable medium, comprising readable program for causing a processing unit in a computer means, such as a controller of a detector arrangement, to control an execution of the steps of any of claims 1-13 and/or 14-16.
19. The computer program product according to claim 18, characterised in that the computer usable medium is a record medium.
20. The computer program product according to claim 18, characterised in that the computer usable medium is a computer memory.
21. The computer program product according to claim 18, characterised in that the computer usable medium is a Read-Only Memory.
22. The computer program product according to claim 18, characterised in the computer usable medium is an electrical carrier signal.
23. Arrangement for reducing noise in an image, preferably a substance raw data image, comprising substance and noise data which is associated to one detection microcavity and one substance retained in the detection microcavity of a microfluidic structure of a microfluidic device, the substance raw data information obtained by use of a detector arrangement capable of scanning and collecting radiation data in a search area associated to a detection area, which is associated to the detection microcavity. the collected radiation data are stored as measure data for the amount of the substance as substance raw data in a set in a data storage, the arrangement also comprises a processing unit for controlling the arrangement, characterized in that the method comprises following step:
- means for reducing background radiation data either by use of an image comprising background radiation data and no substance data information, typically a background radiation image, (step α') or, by use of a calculated median value of the background radiation data (step α").
24. Arrangement according to claim 23, characterized in that the arrangement comprises: - means for accessing a background radiation image and means for generating a background radiation image (step α' l ); - means for accessing corresponding substance raw data image and means for generating a substance raw data image (step α'2);
- means for correlating the two images and generating a correlation image (step α'3);
- means for locate the positions of the largest pixel values (step α'4); - means for positioning the images correctly, means for subtracting corresponding pixel values and means for storing the generated background processed image(step α'5).
25. Arrangement according to claim 23, characterized in that the arrangement comprises:
- means for accessing corresponding substance raw data image and means for generating a background radiation image substance raw data image (step α' " 1 );
- means for determining a Median value for background radiation data (step α'"2);
- means for subtracting pixel values of substance raw data image with the median value contribution, and means for storing the background processed image (step α"3).
26. Arrangement according to claim 23, characterized in that the arrangement comprises means for reducing peak disturbances (Step β).
27. Arrangement according to claim 26, characterized in that the means for reducing peak disturbances (Step β) comprises:
- means for generating a Laplacian processed image from the earlier stored background processed image by using a Laplacian filter in two dimensions (step β 1 );
- means for calculating a disturbance threshold value from the background processed image (step β2);
- means for generating binary encoding Laplacian processed image by use of the disturbance threshold limit (step β3); - means for findingithe start points in clusters (step β4);
- means for determining a minimum rectangular area surrounding the start points (step β5);
- means for calculating an interpolated value within each peak area and substituting the original data of the peaks, and storing the disturbance filtered image (step β6).
28. Arrangement according to claim 23, characterized in that the arrangement comprises means for determining a global threshold value (step χ).
29. Arrangement according to claim 28, characterized in that means for determining a global threshold value (step χ) comprises: - means for generating a y-summation graph (step χl);
- means for Median filtering of the y-summation graph (step χ2);
- means for determining centre position of detection area (step χ3);
- means for determining radial directed edges of detection area (step χ4);
- means for generating a x-summation graph (step χ5); - means for median filtering of the x-summation graph (step χ6);
- means for determining an edge curve (correlation) (step χ7);
- means for restraining disturbance peaks in the edge curve (step χ8);
- means for determining the start and end positions of the detection area (step χ9);
- means for determining a start value of the global threshold value (step χl O); - means for adjusting the global threshold value to the value giving correct detection area width (step χl 1);
- means for generating a binary encoded disturbance filtered image by use of the global threshold value (step χl2).
30. Arrangement according to claim 23, characterized in that the arrangement comprises means for removing binary objects outside the detection area (step δ).
31. Arrangement according to claim 30, characterized in that the means for removing binary objects outside the detection area (step δ) comprises:
-means for labelling of all the objects in the binary coded disturbance filtered image, wherein the number of connected pixels of each object is determined (step δl ); - means for generating a binary detection area image by removing all but the largest objects that are inside the approximate detection area position (step δ2).
32. Arrangement according to claim 23, characterized in that the arrangement comprises means for removing unwanted areas of the detection area (step ε).
33. Arrangement according to claim 32, characterized in that means for removing unwanted areas of the detection area (step ε) comprises:
- means for generating a normalized area image by setting all binary high indicated pixels outside the location of the normalized area ((detection area width) multiplied by (min detection area length)) to binary "low" in the binary detection area image
(step εl);
- means for multiplying each corresponding pixel in the disturbance filtered image, generated and stored earlier in the invented process, and the normalized area image with each other and perform the calculation, e.g. integration, on the normalized area (step ε2).
34. Arrangement according to claim 23. characterized in that the arrangement comprises means for applying default detection area in noisy images (step φ).
35. Arrangement according to claim 34, characterized in that the means for applying default detection area in noisy images comprises: - means for having automatically all detection area positions fulfilling high confidence criterion (step φl);
- means for calculating an average detection area position (step φ2);
- means for using automatically the average detection area position if one or more predetermined condition is fulfilled (step φ3).
36. An arrangement for determining a measure of at least one substance, the arrangement comprises:
(a) a microfluidic device and a detector arrangement, wherein said microfluidic device comprises a plurality of microchannel structures, each of which has an inlet port and a detection microcavity, and a plurality of detection areas, each ofiwhich is associated with one of said detection microcavities: and wherein said detector arrangement is capable of collecting radiation from individual subareas of each of said detection areas;
(b) means for integrating radiation as a function of subarea for each scanned detection area to obtain the amount of radiation from each detection area; and (c) means for characterizing for each of the amounts obtained in step (d) a reaction variable that has been included in the process protocol used for each microchannel structure.
37. The arrangement according to claim 36, at least one of the method steps α - φ in claims 1-13 is carried out in step (d).
38. The arrangement according to claim 37, characterized in that the subarea in step d) is a normalized area.
EP02793749A 2001-12-31 2002-12-30 Method and arrangement for reducing noise in a microfluid device Withdrawn EP1461769A1 (en)

Applications Claiming Priority (7)

Application Number Priority Date Filing Date Title
US62258 1998-04-17
SE0104461A SE0104461D0 (en) 2001-09-17 2001-12-31 Detector arrangement for microfluidic devices
SE0104461 2001-12-31
US10/062,258 US20030054563A1 (en) 2001-09-17 2002-01-31 Detector arrangement for microfluidic devices
PCT/SE2002/001678 WO2003025548A1 (en) 2001-09-17 2002-09-17 Detector arrangement for microfluidic devices
WOPCT/SE02/01678 2002-10-17
PCT/SE2002/002455 WO2003056517A1 (en) 2001-12-31 2002-12-30 Method and arrangement for reducing noise in a microfluid device

Publications (1)

Publication Number Publication Date
EP1461769A1 true EP1461769A1 (en) 2004-09-29

Family

ID=27354789

Family Applications (1)

Application Number Title Priority Date Filing Date
EP02793749A Withdrawn EP1461769A1 (en) 2001-12-31 2002-12-30 Method and arrangement for reducing noise in a microfluid device

Country Status (4)

Country Link
EP (1) EP1461769A1 (en)
JP (1) JP4423039B2 (en)
AU (1) AU2002359227A1 (en)
WO (1) WO2003056517A1 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE0300822D0 (en) 2003-03-23 2003-03-23 Gyros Ab A collection of Micro Scale Devices
SE0301945D0 (en) 2003-06-30 2003-06-30 Gyros Ab Confidence determination
US7718133B2 (en) 2003-10-09 2010-05-18 3M Innovative Properties Company Multilayer processing devices and methods
CN100343655C (en) * 2004-06-18 2007-10-17 中国科学院上海有机化学研究所 Online self-calibration laser induced fluorescence detection method based on electric charge coupling apparatus
EP1779111B1 (en) 2004-07-16 2010-07-07 Gyros Patent Ab Grading of immune responses
JP5021486B2 (en) * 2004-11-10 2012-09-05 ユィロス・パテント・アクチボラグ Liquid detection and reliability determination
WO2007069940A1 (en) 2005-12-12 2007-06-21 Gyros Patent Ab Microfluidic assays and microfluidic devices

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5892539A (en) * 1995-11-08 1999-04-06 Alpha Innotech Corporation Portable emission microscope workstation for failure analysis
WO2000007026A1 (en) * 1998-07-27 2000-02-10 Caliper Technologies, Corp. A distributed database for analytical instruments

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9418981D0 (en) * 1994-09-21 1994-11-09 Univ Glasgow Apparatus and method for carrying out analysis of samples
GB9828785D0 (en) * 1998-12-30 1999-02-17 Amersham Pharm Biotech Ab Sequencing systems
WO2000067635A1 (en) * 1999-05-07 2000-11-16 Applied Spectral Imaging Ltd. Spectral bio-imaging of the eye
US6545758B1 (en) * 2000-08-17 2003-04-08 Perry Sandstrom Microarray detector and synthesizer

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5892539A (en) * 1995-11-08 1999-04-06 Alpha Innotech Corporation Portable emission microscope workstation for failure analysis
WO2000007026A1 (en) * 1998-07-27 2000-02-10 Caliper Technologies, Corp. A distributed database for analytical instruments

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of WO03056517A1 *

Also Published As

Publication number Publication date
AU2002359227A1 (en) 2003-07-15
JP2005513681A (en) 2005-05-12
WO2003056517A1 (en) 2003-07-10
JP4423039B2 (en) 2010-03-03

Similar Documents

Publication Publication Date Title
US7221783B2 (en) Method and arrangement for reducing noise
US20030054563A1 (en) Detector arrangement for microfluidic devices
US6992278B2 (en) Homing process
US8877129B2 (en) Method and device for optical detection of substances in a liquid or gaseous medium
EP0782705B8 (en) Apparatus and method for carrying out analysis of samples
US7428200B2 (en) Method for triggering through disc grooves and related optical analysis discs and system
US20050037484A1 (en) Optical bio-discs including spiral fluidic circuits for performing assays
US20080138247A1 (en) Liquid Detection and Confidence Determination
US20050002827A1 (en) Optical discs including equi-radial and/or spiral analysis zones and related disc drive systems and methods
AU2003209372A1 (en) Method for triggering through disc grooves and related optical analysis discs and system
King et al. Optical detection strategies for centrifugal microfluidic platforms
WO2003056517A1 (en) Method and arrangement for reducing noise in a microfluid device
US8600684B2 (en) Spinner home sequence
JP6883248B2 (en) Particle measuring machine, analyzer, and analysis method
US7616791B2 (en) Confidence determination
CN1556924A (en) Optical analysis disc and related drive assembly for performing interactive centrifugation

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20040507

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SI SK TR

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: GYROS PATENT AB

17Q First examination report despatched

Effective date: 20061115

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20160330