WO2016158719A1 - Procédé de traitement d'images, programme de commande, et dispositif de traitement d'images - Google Patents

Procédé de traitement d'images, programme de commande, et dispositif de traitement d'images Download PDF

Info

Publication number
WO2016158719A1
WO2016158719A1 PCT/JP2016/059549 JP2016059549W WO2016158719A1 WO 2016158719 A1 WO2016158719 A1 WO 2016158719A1 JP 2016059549 W JP2016059549 W JP 2016059549W WO 2016158719 A1 WO2016158719 A1 WO 2016158719A1
Authority
WO
WIPO (PCT)
Prior art keywords
spheroid
circular
separation
image
degree
Prior art date
Application number
PCT/JP2016/059549
Other languages
English (en)
Japanese (ja)
Inventor
治郎 津村
藤本　博己
Original Assignee
株式会社Ｓｃｒｅｅｎホールディングス
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社Ｓｃｒｅｅｎホールディングス filed Critical 株式会社Ｓｃｒｅｅｎホールディングス
Priority to JP2017509890A priority Critical patent/JP6470399B2/ja
Publication of WO2016158719A1 publication Critical patent/WO2016158719A1/fr

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
    • C12M1/00—Apparatus for enzymology or microbiology
    • C12M1/34—Measuring or testing with condition measuring or sensing means, e.g. colony counters
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00—General purpose image data processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes

Definitions

  • the present invention relates to a technique for evaluating a spheroid included in an image by image processing based on an original image obtained by imaging the spheroid.
  • cell clusters in which a large number of cells gather in a spherical shape are cultured and observed.
  • a spheroid with higher cell activity has a shape closer to a sphere. From this, it is performed to evaluate how the spheroid has a shape similar to a circle in an image obtained by imaging the spheroid.
  • Examples of techniques for evaluating the roundness of an object such as a cell included in an image include the techniques described in Patent Documents 1 to 3. In these techniques, the roundness of an object included in an image is evaluated by the roundness obtained from the ratio of the area of the object to the circumference.
  • the above-described conventional evaluation method can achieve a certain result.
  • the object is a spheroid
  • good results are not always obtained.
  • the reason is as follows. Spheroids are composed of a large number of cells, and irregularities due to the contours of individual cells appear on the outer periphery of the spheroids. Therefore, the circumference tends to be estimated longer than a circle of the same size. As a result, the calculated roundness does not necessarily indicate the roundness of the spheroid.
  • This invention is made in view of the said subject, and it aims at providing the technique which can obtain
  • the image processing method specifies a target region including an object corresponding to the spheroid from an original image in which a spheroid is captured, and a circular separability filter operation to determine the target.
  • a step of detecting a circular region having a maximum degree of separation in the region, and a value proportional to the degree of separation of the circular region calculated by the circular separation degree filter calculation, the spheroid specified by the circular region And outputting as an index value indicating roundness.
  • the image processing apparatus includes an image acquisition unit that acquires an original image in which a spheroid is captured, and a target area that includes an object corresponding to the spheroid in the original image.
  • a detection means for executing a circular separability filter calculation to detect a circular area having a maximum value of the separability, and a value proportional to the separability of the circular area calculated by the circular separability filter calculation
  • Output means for outputting as an index value indicating the roundness of the spheroid specified by the region.
  • Another aspect of the present invention is a control program that causes a computer to execute each step of the above-described image processing method.
  • the circle having the maximum degree of separation among the circles drawn in the target region of the original image is specified by the circular separation degree filter calculation.
  • This circle is considered to have a high probability of corresponding to the outline of the spheroid existing in the target region. From this, it can be said that the circular region specified by this circle approximately represents the spheroid region. If the spheroid is close to a circle in the image, the value of the degree of separation is large, and if the deviation from the circle is large, the value of the degree of separation is small. Therefore, the value of the separation degree itself is information that quantitatively indicates how close the spheroid is in the shape of a circle in the image.
  • a value proportional to the value of the degree of separation obtained for the circular region estimated to be a spheroid is output as an index value indicating the roundness of the spheroid.
  • the index value obtained in this way indicates the roundness of the spheroid quantitatively and objectively.
  • an index value that quantitatively represents the roundness of a spheroid can be obtained by performing a circular separability filter operation using an image obtained by imaging the spheroid.
  • FIG. 1 is a diagram illustrating a schematic configuration of an imaging apparatus as an embodiment of an image processing apparatus. It is a 1st figure which shows typically the external appearance of a spheroid. It is a 2nd figure which shows the external appearance of a spheroid typically. It is a 3rd figure which shows the external appearance of a spheroid typically. It is a flowchart which shows operation
  • the example of the original image obtained by imaging a well is shown.
  • region cut out from the original image is shown.
  • region cut out from the original image is shown.
  • region cut out from the original image is shown.
  • FIG. 1 is a diagram showing a schematic configuration of an imaging apparatus as an embodiment of an image processing apparatus according to the present invention.
  • the imaging apparatus 1 is an apparatus that images spheroids (cell conglomerates) that are cultured in a liquid injected into a recess called a well W formed on the upper surface of a well plate WP.
  • spheroids cell conglomerates
  • a well W formed on the upper surface of a well plate WP.
  • XYZ orthogonal coordinate axes are set as shown in FIG.
  • the XY plane represents a horizontal plane
  • the Z axis represents a vertical axis. More specifically, the (+ Z) direction represents a vertically upward direction.
  • the well plate WP is generally used in the fields of drug discovery and biological science, and a plurality of wells W are provided on the upper surface of a flat plate.
  • the cross section of the well W is formed in a substantially circular cylindrical shape, and the bottom surface is transparent and flat.
  • the number of wells W in one well plate WP is arbitrary, for example, 96 (12 ⁇ 8 matrix arrangement) can be used.
  • the diameter and depth of each well W is typically about several mm.
  • the size of the well plate and the number of wells targeted by the imaging apparatus 1 are not limited to these and are arbitrary. For example, those having 12 to 384 holes are generally used.
  • the imaging apparatus 1 can be used not only for well plates having a plurality of wells but also for imaging spheroids cultured in, for example, a flat container called a dish.
  • a predetermined amount of liquid as the medium M is injected into each well W of the well plate WP, and spheroids cultured in the liquid under predetermined culture conditions are imaging targets of the imaging apparatus 1.
  • the medium M may be added with an appropriate reagent, or may be in a liquid state and gelled after being put into the well W. Commonly used liquid volume is about 50 to 200 microliters.
  • the imaging apparatus 1 can target a spheroid cultured in the medium M, for example, on the inner bottom surface of the well W, as an imaging target.
  • the imaging device 1 includes a holder 11, an illumination unit 12, an imaging unit 13, and a control unit 14.
  • the holder 11 abuts on the peripheral edge of the lower surface of the well plate WP that holds the sample together with the liquid in each well W and holds the well plate WP in a substantially horizontal posture.
  • the illumination unit 12 is disposed above the holder 11.
  • the imaging unit 13 is disposed below the holder 11.
  • the control unit 14 includes a CPU 141 that controls operations of these units.
  • the illumination unit 12 emits appropriate diffused light (for example, white light) toward the well plate WP held by the holder 11. More specifically, for example, a combination of a white LED (Light-Emitting-Diode) light source as a light source and a diffusion plate can be used as the illumination unit 12. The spheroid in the well W provided in the well plate WP is illuminated from above by the illumination unit 12.
  • white light for example, white light
  • a white LED (Light-Emitting-Diode) light source as a light source
  • diffusion plate can be used as the illumination unit 12.
  • the spheroid in the well W provided in the well plate WP is illuminated from above by the illumination unit 12.
  • An imaging unit 13 is provided below the well plate WP held by the holder 11.
  • an imaging optical system 131 is disposed immediately below the well plate WP, and the optical axis AX of the imaging optical system 131 is oriented in the vertical direction (Z direction).
  • the imaging unit 13 images the spheroid in the well W. Specifically, light emitted from the illumination unit 12 and incident on the liquid from above the well W illuminates the imaging target. Light transmitted downward from the bottom surface of the well W is incident on the light receiving surface of the image sensor 132 via the imaging optical system 131. An image of the imaging target imaged on the light receiving surface of the imaging element 132 by the imaging optical system 131 is captured by the imaging element 132.
  • the image sensor 132 for example, a CCD sensor or a CMOS sensor can be used, and either a two-dimensional image sensor or a one-dimensional image sensor may be used.
  • the imaging unit 13 can be moved in the XYZ directions by a mechanical control unit 146 provided in the control unit 14. Specifically, the mechanical control unit 146 moves the imaging unit 13 in the X direction and the Y direction based on a control command from the CPU 141. As a result, the imaging unit 13 moves in the horizontal direction with respect to the well W. Further, focus adjustment is performed by moving the imaging unit 13 in the Z direction. When the imaging object in one well W is imaged, the mechanical control unit 146 positions the imaging unit 13 in the horizontal direction so that the optical axis AX coincides with the center of the well W.
  • the imaging device of the imaging unit 13 is a one-dimensional image sensor
  • a two-dimensional image can be captured by scanning the imaging unit 13 in a direction orthogonal to the longitudinal direction of the image sensor.
  • the mechanical control unit 146 moves the illumination unit 12 in the XY direction integrally with the imaging unit 13 as indicated by a dotted arrow in the drawing. That is, the illumination unit 12 is arranged so that the optical center thereof substantially coincides with the optical axis AX of the imaging unit 13, and moves in the XY direction in conjunction with the imaging unit 13 moving in the XY direction.
  • the center of the well W and the light center of the illumination unit 12 are always located on the optical axis of the imaging unit 13. Therefore, it is possible to maintain the imaging condition favorably while keeping the illumination condition for each well W constant.
  • the image signal output from the image sensor of the imaging unit 13 is sent to the control unit 14. That is, the image signal is input to an AD converter (A / D) 143 provided in the control unit 14 and converted into digital image data.
  • the CPU 141 executes various image processing based on the received image data.
  • the control unit 14 further includes an image memory 144 for storing and storing image data, and a memory 145 for storing and storing programs to be executed by the CPU 141 and data generated by the CPU 141. These may be an integral memory.
  • the CPU 141 executes various control processes described later by executing a control program stored in the memory 145.
  • control unit 14 is provided with an interface (IF) unit 142.
  • the interface unit 142 receives an operation input from the user and presents information such as a processing result to the user.
  • the interface unit 142 exchanges data with an external device connected via a communication line.
  • the control unit 14 may be a dedicated device including the hardware described above.
  • a control unit 14 that incorporates a control program for realizing a processing function described later in a general-purpose processing device such as a personal computer or a workstation may be used. That is, a general-purpose computer device can be used as the control unit 14 of the imaging device 1. In the case of using a general-purpose processing device, it is sufficient that the imaging device 1 has a minimum control function for operating each unit such as the imaging unit 13.
  • the imaging device 1 configured as described above has a function of imaging a spheroid cultured in the well W and a function of calculating the roundness of the spheroid based on the captured image.
  • the “roundness” here is an index value that quantitatively indicates how close the spheroid has a circular shape in the captured image.
  • FIG. 2A to 2C are diagrams schematically showing the appearance of the spheroid.
  • the spheroid Sp1 shown in FIG. 2A is an example composed of cells with relatively high viability (viability, cell activity).
  • viability viability, cell activity
  • a shape close to a sphere is formed by a large number of cells C. Therefore, in the two-dimensional image obtained by imaging the spheroid Sp1, the outline of the spheroid Sp1 is substantially circular.
  • roundness or “circularity” is generally used as an index value indicating how close the contour of the imaged object is to a circle, that is, the “roundness” of the spheroid.
  • the object has a shape that is close to a circle, depending on the ratio of the radius of the circle that circumscribes the contour of the object or the ratio of the area of the object and the circumference of the contour. Is expressed.
  • spheroids are a collection of many cells. Therefore, as shown in FIG. 2A, even if the spheroid Sp1 appears to be approximately circular, microscopically, irregularities due to the surfaces of the individual cells C appear in the outline. Therefore, there are cases where the “roundness” of the spheroid cannot be appropriately represented by an index value using roundness or circularity.
  • an actual spheroid has a three-dimensional shape, as shown in FIG. 2C, a shadow caused by illumination light may cause the spheroid Sp1 image to be shaded. In this case, the contour of the spheroid cannot be appropriately extracted, and as a result, the roundness or circularity obtained for the spheroid may have an inappropriate value.
  • the “roundness” of the spheroid is introduced as a new index value.
  • the operation of the imaging apparatus 1 from imaging of the spheroid to calculation of the roundness will be described.
  • FIG. 3 is a flowchart showing the operation of this imaging apparatus. This operation is realized by causing the CPU 141 of the control unit 14 to execute a control program prepared in advance in accordance with a user instruction and causing each unit of the imaging apparatus 1 to perform a predetermined operation defined by the program.
  • the imaging unit 13 images the well W (step S102). By performing imaging a plurality of times as necessary, each well W is imaged and an original image is acquired.
  • the object region having characteristics different from the background region is extracted from the captured original image by image processing (step S103).
  • Various methods for extracting an object region from an image by image processing are known, and an appropriate method can be applied in this embodiment as well.
  • an area in which pixel values of pixels constituting the original image are in a predetermined range can be set as an object area.
  • an area surrounded by edges extracted by edge extraction can be set as an object area.
  • step S103 may be omitted.
  • a setting input for a target region to be subjected to the following processing is received from the captured original image (step S104).
  • the setting input is given to the control unit 14 when the user operates a keyboard, a mouse, a tablet, or the like as an input device provided in the interface unit 142, for example.
  • FIG. 4A to 4D are diagrams schematically showing examples of images. More specifically, FIG. 4A shows an example of an original image obtained by imaging the well W. 4B to 4D show examples of images of target regions cut out from the original image. As shown in FIG. 4A, spheroids are distributed at various locations in the well W in the original image. Here, seven spheroids Sp3 to Sp9 are included, but the number of spheroids varies, and of course, one may be sufficient.
  • the original image can include spheroids having various characteristics. For example, spheroids having a shape close to a circle (Sp4 to Sp8), spheroids having a more irregular shape (Sp3, Sp9), spheroids (Sp4 to Sp6) that are close to each other, and the like may be included.
  • the user designates a partial region including a spheroid for which the degree of roundness is to be obtained from such an original image as a target region. Then, a partial image in the target region is cut out from the original image, and a circular separation degree filter calculation is executed for the partial image (step S105).
  • the target area is preferably selected so that a partial image obtained by cutting out the target area from the original image includes one spheroid for which the roundness is to be obtained.
  • spheroids Sp8 and Sp9 have a relatively large distance from other spheroids and are isolated.
  • the target region is set so that the partial image P includes the entire spheroid and the surrounding background region.
  • the spheroid Sp5 there are other spheroids in the vicinity.
  • the target region is selected so that the partial image P includes the entire spheroid and other spheroids are excluded as much as possible.
  • the target area is designated by the user here, the target area may be automatically set by the imaging apparatus 1 instead of or in addition to this. That is, when the object area is extracted from the original image in step S103, the position and size of the object area in the original image can be known. Therefore, it is possible to automatically set the target area so as to include the extracted object area.
  • FIG. 5 is a flowchart showing the processing contents of the circular separation degree filter calculation in this embodiment.
  • the average value of the pixel values included in the partial image corresponding to the designated target area, that is, the average pixel value is calculated (step S201).
  • a difference image is created in which the absolute value of the difference between the pixel value of each pixel of the original partial image and the average pixel value is the pixel value of the pixel (step S202).
  • the difference image is a virtual intermediate image used for subsequent processing. It is only necessary to obtain the pixel value of each pixel in the difference image, and it is not necessary to actually create or output an image.
  • the contrast between the object region and the other background region may be smaller than that of the original image.
  • the center coordinate position and diameter of a virtual circle are provisionally set to appropriate values as parameters for calculating the degree of separation (step S203). If the object area extracted in step S103 is present in the target area, the position and size are known in advance, and parameter setting can be performed based on the information. That is, it is only necessary to set a circle whose center position and size are equal to those of the object area.
  • the parameter setting value may be given by the user.
  • the degree of separation is calculated for a circle in the target area specified by the set center coordinate position and diameter (step S204).
  • the degree of separation is calculated based on the pixel value of a pixel included in an annular outer region adjacent to the outside of the circle and the pixel value of a pixel included in an annular inner region adjacent to the inside of the circle. Is required.
  • the width of these annular regions can be, for example, several pixels.
  • Steps S203 and S204 are a loop process, and the calculation of the degree of separation is repeated while variously changing the combination of the center position and the diameter of the circle.
  • the calculation may be performed for all circles that can be drawn in the target area.
  • the center position of the circle and the change range of the diameter can be limited based on the information. By doing so, the processing time can be significantly shortened compared with the case where the calculation is performed for all combinations of the center coordinate position and the diameter.
  • step S205 When the degree of separation is obtained for all the circles to be set, the loop process ends (step S205).
  • the maximum value of the degree of separation thus obtained (hereinafter referred to as “maximum degree of separation”) is stored in the memory 145 together with the center position and diameter of the circle having the maximum value (step S206). As a result, a circle having the maximum degree of separation in the target region is specified.
  • the circle thus identified can be presumed to schematically indicate the outline of the spheroid existing in the target region. That is, there is a high probability that a circular region having the circle as an outer periphery is a spheroid region in the target region. In particular, when the spheroid has a shape close to an ideal circle, the circular region should substantially match the spheroid region.
  • FIG. 6A to 6C are diagrams showing the relationship between the shape of the spheroid and the corresponding circular region.
  • the spheroid Sp has a substantially circular shape like the spheroid Sp8 shown in FIG. 4A
  • the obtained circular region CR and the actual spheroid Sp region substantially coincide as shown in FIG. 6A.
  • the obtained value of the maximum degree of separation is relatively large.
  • the target region is set so as to mainly include a single spheroid.
  • the value of the degree of separation in the circular region increases as the spheroid is closer to a circle. Accordingly, for example, a value proportional to the value of the maximum separation degree, which takes a larger value as the maximum separation degree increases, can be used as an index value representing the roundness of the spheroid. In the index value of the roundness based on the degree of separation, if the envelope shape of the spheroid is substantially circular, the influence of fine irregularities on the outer periphery is less likely to appear.
  • step S105 when a circle having the maximum degree of separation in the target area is specified, the roundness of the spheroid in the target area is calculated based on the value of the maximum degree of separation (step S106).
  • a value proportional to the value of the maximum separation can be used. For example, a value obtained by normalizing the value of the maximum degree of separation as it is or multiplying by a predetermined proportionality coefficient can be set as the value of the roundness.
  • the degree of roundness may be quantified by an appropriate calculation formula that monotonously increases with respect to the value of the maximum degree of separation.
  • the value of the maximum degree of separation obtained by the above calculation tends to increase as the concentration of the spheroid region increases if the concentration of the background region is the same. Therefore, when comparing the roundness of spheroids having greatly different concentrations, it is not appropriate to use the value of maximum separation as an index value as it is. In order to avoid this problem, for example, a value obtained by dividing the value of the maximum separation by the average pixel value in the spheroid region and reducing the influence of the density difference may be used as the roundness.
  • the roundness of the spheroid thus obtained is output in a predetermined format and presented to the user.
  • the user can be notified by displaying the result on a display (output device) provided in the interface unit 142 (step S107).
  • FIG. 7 is a diagram showing an example of a display image.
  • the display mode is not limited to that shown here, but is arbitrary. Further, the roundness calculation result may be output in a manner other than being displayed on the display device in this way, for example, transmitted and output to an external recording medium or an information processing device.
  • the original image is displayed in the box area B1 of the display image IM, and the range of the set target area is indicated by a dotted line.
  • the box area B2 an enlarged image of the target area is displayed.
  • a circular region corresponding to the maximum degree of separation specified by the circular degree of separation filter calculation is indicated by a dotted line.
  • the box area B3 the calculation results of the center coordinate position, the diameter, and the roundness of the spheroid of the circular area are displayed as text. Thus, the roundness is calculated for one target area, and the result is displayed.
  • step S108 it is determined whether another target area is set and calculation of the roundness in that area is to be continued. If it is necessary to continue the process, the process returns to step S104 and the above process is repeated.
  • the user can be inquired whether to continue the process.
  • the target area is automatically set, if there is an extracted object area whose roundness has not yet been calculated, the process for the object area is continued. Then, if the roundness calculation processing has been completed for all object regions, the processing may be terminated.
  • the circular separation degree filter calculation is performed on the target region designated to include the spheroids in the original image obtained by imaging.
  • a value proportional to the value of the separation degree when the separation degree is maximized is output as an index value indicating the roundness of the spheroid in the target region.
  • the roundness can be calculated without specifying the contour in advance. In this sense, the roundness can be calculated even if step S103 in FIG. 3 is omitted. In particular, when it is clear that one spheroid is included in the target region, it is possible to omit specifying the contour.
  • the circular separability filter calculation is not used for the purpose of searching a circular area from the original image. That is, on the premise that the position of the object corresponding to the spheroid is specified to some extent, the circular separability filter calculation is used to obtain the maximum value of the separability within the target area including this object.
  • the roundness of the spheroid obtained in this way becomes a large value when the envelope shape of the spheroid is close to a circle, and becomes a smaller value as it deviates from the circle. Since the deviation from the spheroid-shaped circle is mainly caused by the weakening of the spheroid-forming cells, the roundness value can be used for the evaluation of spheroid activity.
  • the imaging device 1 functions as the “image processing device” of the present invention.
  • the imaging unit 13 functions as the “image acquisition unit” of the present invention
  • the CPU 141 functions as the “detection unit” of the present invention
  • the interface unit 142 functions as the “output unit” and the “accepting unit” of the present invention.
  • the image processing apparatus of the above embodiment includes the imaging unit 13 as the “image acquisition unit” of the present invention.
  • an image processing apparatus that does not include an imaging function may be used. That is, the present invention is implemented in such a manner that the interface unit 142 in the control unit 14 of the above embodiment acquires an original image by receiving image data from the outside, and the roundness is calculated using the original image. May be.
  • the interface unit 142 functions as an “image acquisition unit”, and an imaging function is not required.
  • the present invention is implemented by the CPU 141 executing the control program stored in the memory 145.
  • a general-purpose computer device can be used as the control unit 14 in this embodiment. Therefore, the present invention may be provided to the user as a control program for causing the computer apparatus to execute the above-described processing on the premise that the program is read by such a computer apparatus, and in a mode in which the program is recorded on an appropriate recording medium. Is possible. Accordingly, for example, it is possible to add a new function to an imaging apparatus that is already in operation. Further, a general-purpose computer device can function as the “image processing device” of the present invention.
  • the specific embodiment has been described by way of example, and in the circular separation degree filter calculation according to the present invention, at least one of the center position and the diameter in the target region is changed and set in a plurality of stages.
  • the degree of separation may be determined. By doing so, it is possible to accurately select a circle having the maximum degree of separation from various circles that can be set in the target region.
  • a difference between the pixel value of the pixel and the average pixel value of the target area is obtained, and a circular separation filter operation is performed on the difference image using the difference as the pixel value of the pixel. Also good. By doing so, it is possible to reduce the influence of density unevenness in the object region and the background region in the target region on the calculation result of the separation degree.
  • a value obtained by normalizing the value of the degree of separation of the circular area with the average pixel value of the pixels in the circular area may be used as an index value indicating roundness.
  • concentration is different for each spheroid
  • the value of the degree of separation is different even if the spheroids have the same roundness. This problem can be solved by normalizing the degree of separation with the average pixel value of the circular area corresponding to the spheroid.
  • it may be configured such that a setting input for specifying the target area from the original image is received from the user, and the target area is specified according to the setting input. By doing so, it is possible to preferentially process an object that the user desires to evaluate and present the calculation result to the user in a short time.
  • a step of detecting an object corresponding to the spheroid from the original image by image processing may be provided, and the target region may be specified based on the detection result. By doing so, it is possible to automatically perform calculation for each object in the original image regardless of user input.
  • the image acquisition means may have a configuration having an imaging unit that captures a sample and acquires an original image.
  • the original image may be given from the outside, but by providing the imaging unit, it is possible to consistently perform from the imaging of the original image to the calculation of the roundness of the spheroid.
  • the present invention is suitable for use in evaluating spheroids for the purpose of observation, analysis, and the like, and can be used for various experiments in the medical and biological science fields, for example.
  • Imaging device image processing device
  • Imaging unit image acquisition means
  • control unit CPU (detection means)
  • Interface unit output means, accepting means

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Biotechnology (AREA)
  • Organic Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Immunology (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Microbiology (AREA)
  • Sustainable Development (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Medicinal Chemistry (AREA)
  • Genetics & Genomics (AREA)
  • Mathematical Physics (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

La présente invention consiste : en une étape de spécification d'une zone sujet incluant un objet correspondant à un sphéroïde, à partir d'une image brute capturée du sphéroïde (étape S103) ; en une étape de détection d'une zone circulaire présentant le degré de séparation maximal dans la zone sujet, au moyen d'une opération de filtre de degré de séparation circulaire (étape S105) ; et en une étape de sortie d'une valeur proportionnelle au degré de séparation circulaire de la zone circulaire calculée au moyen de l'opération de filtre de degré de séparation circulaire, la valeur étant sortie sous forme de valeur d'indice indiquant la rotondité du sphéroïde spécifié par l'aire circulaire (étape S106).
PCT/JP2016/059549 2015-03-31 2016-03-25 Procédé de traitement d'images, programme de commande, et dispositif de traitement d'images WO2016158719A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017509890A JP6470399B2 (ja) 2015-03-31 2016-03-25 画像処理方法、制御プログラムおよび画像処理装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015-071655 2015-03-31
JP2015071655 2015-03-31

Publications (1)

Publication Number Publication Date
WO2016158719A1 true WO2016158719A1 (fr) 2016-10-06

Family

ID=57005836

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/059549 WO2016158719A1 (fr) 2015-03-31 2016-03-25 Procédé de traitement d'images, programme de commande, et dispositif de traitement d'images

Country Status (2)

Country Link
JP (1) JP6470399B2 (fr)
WO (1) WO2016158719A1 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018110058A1 (fr) * 2016-12-13 2018-06-21 株式会社サイフューズ Dispositif de production d'une structure cellulaire
WO2019106945A1 (fr) * 2017-11-28 2019-06-06 パナソニックＩｐマネジメント株式会社 Dispositif de détermination d'état de culture, procédé et programme de détermination d'état de culture

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001116534A (ja) * 1999-10-15 2001-04-27 Kubota Corp 真円度判定方法、真円度演算装置及び記録媒体
JP2005196385A (ja) * 2004-01-06 2005-07-21 Toshiba Corp 画像処理装置、画像処理方法およびデジタルカメラ
JP2008217614A (ja) * 2007-03-06 2008-09-18 Fujitsu Ten Ltd 顔方位判定装置
JP2009077635A (ja) * 2007-09-25 2009-04-16 Yokohama City Univ 軸索内移動粒子の自動追跡システム
WO2014038025A1 (fr) * 2012-09-05 2014-03-13 Ｓｃｉｖａｘ株式会社 Procédé de criblage de matière agissant sur une maintenance épithéliale de cellules
JP2014109488A (ja) * 2012-11-30 2014-06-12 Dainippon Screen Mfg Co Ltd 画像処理装置、画像処理方法、および画像処理プログラム

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013137635A (ja) * 2011-12-28 2013-07-11 Dainippon Screen Mfg Co Ltd 画像表示装置および画像表示方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001116534A (ja) * 1999-10-15 2001-04-27 Kubota Corp 真円度判定方法、真円度演算装置及び記録媒体
JP2005196385A (ja) * 2004-01-06 2005-07-21 Toshiba Corp 画像処理装置、画像処理方法およびデジタルカメラ
JP2008217614A (ja) * 2007-03-06 2008-09-18 Fujitsu Ten Ltd 顔方位判定装置
JP2009077635A (ja) * 2007-09-25 2009-04-16 Yokohama City Univ 軸索内移動粒子の自動追跡システム
WO2014038025A1 (fr) * 2012-09-05 2014-03-13 Ｓｃｉｖａｘ株式会社 Procédé de criblage de matière agissant sur une maintenance épithéliale de cellules
JP2014109488A (ja) * 2012-11-30 2014-06-12 Dainippon Screen Mfg Co Ltd 画像処理装置、画像処理方法、および画像処理プログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MASAMU TAKAKI ET AL.: "Enkei Bunrido Filter ni yoru Ryushi Kenshutsu to Kinbo Ryushi no Ido Kyori Saishoka ni Motozuku Tsuiseki -Ryoiki Base no Ryushi Kenshutsu to Ryushi no Kinbo Joho o Koryo shita Tsuiseki", VIEW VISION ENGINEERING WORKSHOP KOEN RONBUNSHU, December 2010 (2010-12-01), pages 483 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018110058A1 (fr) * 2016-12-13 2018-06-21 株式会社サイフューズ Dispositif de production d'une structure cellulaire
WO2018110053A1 (fr) * 2016-12-13 2018-06-21 株式会社サイフューズ Dispositif de production de structure cellulaire et plateau à cellules
JPWO2018110053A1 (ja) * 2016-12-13 2019-03-28 株式会社サイフューズ 細胞構造体製造装置および細胞トレイ
JPWO2018110058A1 (ja) * 2016-12-13 2019-10-24 株式会社サイフューズ 細胞構造体製造装置
WO2019106945A1 (fr) * 2017-11-28 2019-06-06 パナソニックＩｐマネジメント株式会社 Dispositif de détermination d'état de culture, procédé et programme de détermination d'état de culture
US11674889B2 (en) 2017-11-28 2023-06-13 Panasonic Intellectual Property Management Co., Ltd. Culture state determination based on direction-dependent image information

Also Published As

Publication number Publication date
JPWO2016158719A1 (ja) 2017-10-26
JP6470399B2 (ja) 2019-02-13

Similar Documents

Publication Publication Date Title
US9865054B2 (en) Evaluation method of spheroid and spheroid evaluation apparatus
JP5920994B2 (ja) マイクロプレートのウェル壁境界を識別するための方法及びシステム
WO2013099045A1 (fr) Dispositif d'affichage d'image et procédé d'affichage d'image
JP2013526717A5 (fr)
EP3375859B1 (fr) Procédé de construction de classificateur et procédé de détermination de la vie ou de la mort de cellules l'employant
Wang et al. Recognition and matching of clustered mature litchi fruits using binocular charge-coupled device (CCD) color cameras
JP2010181402A (ja) 受精卵品質評価支援システム、受精卵品質評価支援装置および受精卵品質評価支援方法
JP6342810B2 (ja) 画像処理
KR102580984B1 (ko) 화상 처리 방법, 프로그램 및 기록 매체
JP6882329B2 (ja) Ｉｈｃ画像解析のための空間インデックス作成
JP6470399B2 (ja) 画像処理方法、制御プログラムおよび画像処理装置
JP6229365B2 (ja) コロニーカウント装置、コロニーカウント方法、およびコロニーカウントプログラム
Faria et al. AI-powered mobile image acquisition of vineyard insect traps with automatic quality and adequacy assessment
JP4985480B2 (ja) がん細胞を分類する方法、がん細胞を分類するための装置及びがん細胞を分類するためのプログラム
JP7006833B2 (ja) 細胞解析装置
Yin et al. A novel method of situ measurement algorithm for oudemansiella raphanipies caps based on YOLO v4 and distance filtering
US20180189606A1 (en) Method and device for supporting creation of teaching data, program and program recording medium
JP5762315B2 (ja) 画像処理方法
Niederlein et al. Image analysis in high content screening
JP6326445B2 (ja) 画像処理方法、画像処理装置、および画像処理プログラム
WO2017069035A1 (fr) Procédé de traitement d'image et procédé de création de données de référence d'ombrage
EP3779463B1 (fr) Procédé de traitement d'image, dispositif de traitement d'image, programme et support d'enregistrement
US11181463B2 (en) Image processing device, cell recognition apparatus, cell recognition method, and cell recognition program
WO2017056600A1 (fr) Procédé de traitement d'image et programme de commande
Oduntan High Throughput Phenotyping of Cross-Sectional Morphology in Plant Stalks: Sample Preparation, Image Processing, and Feature Extraction

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16772600

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017509890

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16772600

Country of ref document: EP

Kind code of ref document: A1