WO2016084420A1 - ディジタルホログラフィ装置及びディジタルホログラム生成方法 - Google Patents

ディジタルホログラフィ装置及びディジタルホログラム生成方法 Download PDF

Info

Publication number
WO2016084420A1
WO2016084420A1 PCT/JP2015/070019 JP2015070019W WO2016084420A1 WO 2016084420 A1 WO2016084420 A1 WO 2016084420A1 JP 2015070019 W JP2015070019 W JP 2015070019W WO 2016084420 A1 WO2016084420 A1 WO 2016084420A1
Authority
WO
WIPO (PCT)
Prior art keywords
wavelength
light
hologram
max
wavelengths
Prior art date
Application number
PCT/JP2015/070019
Other languages
English (en)
French (fr)
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 CN201580064772.4A priority Critical patent/CN107003638B/zh
Priority to US15/529,538 priority patent/US10281877B2/en
Priority to EP15863320.6A priority patent/EP3223083A4/en
Priority to JP2016561420A priority patent/JP6388036B2/ja
Publication of WO2016084420A1 publication Critical patent/WO2016084420A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0465—Particular recording light; Beam shape or geometry
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G03H2001/005—Adaptation of holography to specific applications in microscopy, e.g. digital holographic microscope [DHM]
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • G03H2001/0447—In-line recording arrangement
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/26—Processes or apparatus specially adapted to produce multiple sub- holograms or to obtain images from them, e.g. multicolour technique
    • G03H1/2645—Multiplexing processes, e.g. aperture, shift, or wavefront multiplexing
    • G03H2001/266—Wavelength multiplexing
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H2222/00—Light sources or light beam properties
    • G03H2222/10—Spectral composition
    • G03H2222/13—Multi-wavelengths wave with discontinuous wavelength ranges
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H2222/00—Light sources or light beam properties
    • G03H2222/34—Multiple light sources
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H2226/00—Electro-optic or electronic components relating to digital holography
    • G03H2226/02—Computing or processing means, e.g. digital signal processor [DSP]

Definitions

  • the present invention relates to a digital holography device and a digital hologram generation method for calculating a reproduced image of an object from a hologram in which interference fringes between an object wave and a reference wave are recorded, and in particular, to calculate phase information from a plurality of recorded holograms.
  • the present invention relates to a digital holography apparatus and a digital hologram generation method.
  • a predetermined light beam is generated based on an interference pattern (hologram) formed on a detection surface of an image sensor or the like by object light reflected or transmitted by the light source from the light source and reference light directly reaching from the same light source.
  • a reproduced image of the object is obtained by performing arithmetic processing.
  • a hologram is an image representing the intensity distribution of interference light, and this intensity distribution is determined based on the amplitude difference and phase difference between the object wave, which is a complex wavefront, and the reference wave. However, the phase that is the imaginary component of the wave cannot be recorded on the hologram.
  • a false image (conjugate image) having an amplitude complex with the real image is calculated at a position symmetrical to the real image (the desired reconstructed image) with respect to the detection surface. Is done. When these two images are superimposed, the visibility of the real image is reduced, and a clear image to be observed may not be obtained.
  • Patent Document 1 describes a method of imaging a plurality of holograms with the phase of a reference wave shifted and calculating phase information by linear calculation.
  • the method described in this document requires a piezoelectric element for performing phase shift with high accuracy and an optical system such as a mirror and a lens, which complicates the device configuration and increases the size of the device itself.
  • Non-Patent Document 1 describes a phase recovery method using an in-line type digital holography apparatus that does not require such a complicated optical system.
  • a highly light-transmitting object for example, a cell
  • a plurality of holograms having different distances from the object to the detection surface are generated.
  • phase information is calculated by repeatedly calculating propagation / reverse propagation of the light wave between the object surface and the detection surface based on the plurality of captured holograms.
  • a driving mechanism such as a piezoelectric element for moving the detector with high accuracy is still necessary.
  • Non-Patent Document 2 a phase recovery method in which the wavelength of light emitted from a light source is changed in a plurality of stages and light wave propagation is repeatedly calculated based on holograms imaged at the respective wavelengths. According to this method, phase recovery is possible with a simple apparatus configuration that does not require a complicated optical system or drive system.
  • the biomicroscope which is one of the applications of digital holography technology, has so far focused on making it possible to observe intracellular structures or objects that are much smaller than that, that is, improving resolution.
  • phase contrast microscope Although such confirmation is made by microscopic observation, although the phase contrast microscope is used because the cells themselves are difficult to visually recognize, it is not easy to judge the occurrence of a defect, which is a heavy burden on the observer. It has become. Therefore, it is considered to automate this by using digital holography technology.
  • the above-mentioned defects that occur during the culture of induced pluripotent stem cells are judged based on the state of the macro structure (specifically, the shape and texture of the colony including thickness) compared to the intracellular structure.
  • the state of the macro structure specifically, the shape and texture of the colony including thickness
  • the above-described conjugate image is superimposed on the target region and the visibility of the real image is reduced. Confirmed by calculation. That is, it is difficult to say that the configuration described in Non-Patent Document 2 is suitable for determining good or bad in the culture of induced pluripotent stem cells.
  • Non-Patent Document 2 Even if the method described in Non-Patent Document 2 is used, it is difficult to completely recover the phase, and it is difficult to avoid conjugation of conjugate images in inline digital holography. This is because a numerical value serving as a clear index indicating the degree of phase recovery is not defined.
  • a correlation coefficient of estimated amplitude values between holograms with a plurality of wavelengths is used as an evaluation index of approximation accuracy.
  • the correlation coefficient as the evaluation index can be used as a measure of convergence of the iterative calculation, but is not accurate as an index representing the accuracy of phase recovery.
  • the present invention has been made in view of the above circumstances, and an object of the present invention is to provide a digital holography device and a digital hologram generation method capable of reproducing a clear image for an observation object having a different size according to the purpose. Is to provide.
  • the phase information calculation process is an algorithm that estimates the phase of an object wave by light wave propagation calculation based on a plurality of captured holograms. An example is shown in FIG. Details of each step will be described later. It should be noted that the term “light wave propagation” in this specification includes the concepts of both propagation and back propagation described above unless otherwise specified.
  • E (x, y, 0) is the complex amplitude distribution on the object plane
  • E (x, y, z) is the complex amplitude distribution on the detection plane
  • z corresponds to the propagation distance
  • K represents the wave number.
  • the above equation (1) corresponds to obtaining an object image as a sum of diffraction gratings having various pitches by diffraction calculation.
  • the diffraction angle ⁇ of the light incident on the diffraction grating is expressed by the following formula (2).
  • ⁇ represents the wavelength
  • the fine pitch of the diffraction grating corresponds to the fine structure in the image.
  • a fine structure for example, intracellular structure
  • a large structure for example, the shape of a colony
  • FIGS. 7A and 7B are conceptual diagrams illustrating a phenomenon in which a real image and a conjugate image reproduced by light wave propagation calculation are superimposed.
  • a conjugate image having a fine structure corresponding to a high-frequency grating is greatly spread (blurred) during back-propagation to the object plane, and thus has an effect on the visibility of a real image superimposed on the object plane.
  • Small FIG. 7A
  • the conjugate image having a large structure corresponding to the low-frequency grating has a small spread as described above, and thus reaches the object plane while maintaining almost the same size as the real image (FIG. 7B). In this case, a similar line is superimposed on the object plane in the vicinity of each boundary line of the real image, and the visibility of the real image is greatly reduced.
  • the wavelength condition of illumination light for recovering the phase of the low frequency component will be described next.
  • the diffraction angle of the low-frequency component is small compared to the high-frequency component, when the optical path length difference between the detection surfaces H1 and H2 is small (d1 in the figure), the amount of interference fringe deviation due to diffraction is very small.
  • the amount of shift at d1 of the low frequency component is so small that it cannot be shown with respect to the amount of shift g1 related to the high frequency component.
  • the difference between the two holograms imaged on the detection surfaces H1 and H2 is difficult to appear for the low frequency component, and the phase recovery is difficult. Therefore, in order to easily recover the phase of the low frequency component, it is only necessary to increase the amount of interference fringe by increasing the optical path length difference (d2 in the figure) like the detection surfaces H1 and H3 in the figure. (G2 in d2).
  • these multiple detection surfaces can be regarded as conceptually showing detection surfaces at the same position for detecting illumination light of different wavelengths, and detection surfaces close to the object surface have short wavelengths and detections far away.
  • the surface corresponds to long wavelengths.
  • An optical path length difference ⁇ z AB generated when two wavelengths ⁇ A and ⁇ B ( ⁇ A ⁇ B ) are used is expressed by the following equation (3).
  • z represents the propagation distance as in the above formula (1). According to the above equation (3), the greater the difference between the two wavelengths ⁇ A and ⁇ B , the greater the optical path length difference ⁇ z AB, and thus the phase recovery of the low frequency component becomes easier.
  • the structure to be observed corresponds to a lattice having a pitch width p lim .
  • the influence on the visibility of the superimposed real image is sufficiently small if the wavefront of the conjugate image spreads more than X times the pitch width p lim on the object plane.
  • the following expression (4) Holds. 2z is twice the propagation distance z and corresponds to the calculated distance from the conjugate image to the object plane.
  • the following formula (5) is derived from the above formula (4).
  • phase of the component corresponding to the grating having the pitch width p lim that satisfies the above equation (5) may be recovered.
  • the conjugate image of the fine structure spreads larger than the coarse structure, if phase recovery of the component corresponding to the lattice of p lim that is the observation target is performed, It is considered that the influence of the high frequency component having a pitch width p smaller than p lim is sufficiently small.
  • the difference between the holograms is most noticeable, that is, phase recovery of the component corresponding to the pitch width p lim is facilitated by two interferences.
  • the digital holography device made to solve the above-mentioned problems is a) a light source that emits coherent light; b) a detector for imaging a hologram that is an interference fringe between an object wave and a reference wave caused by light emitted from the light source; c) wavelength setting means for setting a plurality of wavelengths of illumination light for generating a hologram to be imaged by the detector, and
  • the wavelength setting means enlarges the conjugate image set by the user as not to hinder the visual recognition of the real image and the conjugate image reproduced by the predetermined computing means for the structure to be observed when superimposed on the corresponding real image.
  • the shortest wavelength ⁇ min and the longest wavelength ⁇ max of the plurality of wavelengths set by the wavelength setting means are ⁇ max / ⁇ min ⁇ (1 / X + 1) (9) It is characterized by setting to satisfy.
  • “does not disturb the visual recognition of the superimposed real image” means, in other words, that the boundary line (edge) particularly noticed in the real image is clear.
  • this conjugate image does not hinder the real image from being visually recognized. Be taken care of.
  • the hologram which is an interference fringe of the object wave and reference wave resulting from the coherent light radiate
  • the wavelength setting means for determining the plurality of wavelengths sets the plurality of wavelengths so that the shortest wavelength ⁇ min and the longest wavelength ⁇ max satisfy the above formula.
  • the phase recovery accuracy for the structure to be observed is improved. Further, even if a conjugate image remaining without phase recovery is superimposed on the real image, the conjugate image is enlarged by X times or more of the real image, so that the visual recognition of the real image is not hindered. Therefore, a clear image can be obtained even for a relatively large observation target such as a cell colony.
  • the calculation of the phase information may be performed by arithmetic means included in the digital holography apparatus, or may be performed on an external computer for control / analysis connected to the digital holography apparatus.
  • the number of wavelengths to be used may be appropriately determined from the viewpoint of hardware, calculation cost, required image quality, and the like. Increasing the number of wavelengths used increases the calculation cost, but enables phase recovery in a wider frequency range.
  • the wavelength setting means of the present invention may be one of the functions provided in the light source, or may be realized on the detector side.
  • a typical example of the former is a light source capable of switching the wavelength of emitted light, and a configuration in which a detector is realized by a hyperspectral camera and images of a plurality of wavelengths are individually extracted from a captured image is considered as the latter example.
  • a switching element connected to a plurality of light sources by an optical fiber or the like may be provided, and a configuration in which an object is irradiated with light emitted from which light source can be switched.
  • a digital holography device made to solve the above-mentioned problems is as follows. a) a light source that emits coherent light; b) a detector for imaging a hologram that is an interference fringe between an object wave and a reference wave caused by light emitted from the light source; c) wavelength setting means for setting a plurality of wavelengths of illumination light for generating a hologram to be imaged by the detector, and
  • the wavelength setting means is configured such that the shortest wavelength ⁇ min and the longest wavelength ⁇ max of the plurality of wavelengths are ⁇ max / ⁇ min ⁇ 1.3 (10) It is characterized by setting to satisfy.
  • the wavelength setting means includes: 1.3 ⁇ ⁇ max / ⁇ min ⁇ 2.0 (11)
  • the shortest wavelength ⁇ min and the longest wavelength ⁇ max are set so as to satisfy the above condition.
  • a digital hologram generating method made to solve the above-described problem is a) a light emitting step of emitting coherent light from a light source; b) an imaging step in which a detector images a hologram that is an interference fringe between an object wave and a reference wave caused by light emitted from the light source; c) a wavelength setting step for setting a plurality of wavelengths of illumination light for generating a hologram to be imaged by the detector, and Based on the magnification X of the conjugate image set by the user as not to hinder the visual recognition of the real image and the conjugate image reproduced by the predetermined calculation means for the structure to be observed when superimposed on the corresponding real image, A plurality of wavelengths set in the wavelength setting step, the shortest wavelength ⁇ min and the longest wavelength ⁇ max are ⁇ max / ⁇ min ⁇ (1 / X + 1) (12) It is characterized by setting to satisfy.
  • the wavelengths of the plurality of light sources are set to 1.3 ⁇ ⁇ max / ⁇ min ⁇ 2.0 as described above. This is also useful in the phase unwrapping method.
  • multiple wavelengths are set to 1.3 ⁇ ⁇ max / ⁇ min ⁇ 2.0 It is a preferable condition even when the phase unwrapping method is applied.
  • the ratio between the maximum value and the minimum value of the wavelength used is 1.3 or more.
  • a relatively clear image can be obtained for the target structure.
  • the method according to the present invention can be used when applying the phase unwrapping method.
  • FIG. 1 is a block diagram showing a schematic configuration of a measurement system including a digital holography device according to an embodiment of the present invention.
  • the flowchart which shows the flow of the process which the measurement system shown in FIG. 1 performs.
  • the flowchart which shows an example of the flow of a hologram imaging process.
  • the flowchart which shows an example of the flow of a phase information calculation process.
  • Explanatory drawing which shows the difference in the diffraction angle of incident light by the difference in the pitch width of a diffraction grating.
  • the conceptual diagram for demonstrating the phenomenon in which the real image reproduced
  • FIG. 1 is a block diagram showing a schematic configuration of a measurement system including a digital holography device according to an embodiment of the present invention.
  • the measurement system includes a digital holography device 100 and a workstation 1 that is communicably connected to the digital holography device 100.
  • the digital holography apparatus 100 is a microscope, and includes N laser diodes (LD) 101 (1) to 101 (N), a switching element 102 (corresponding to the wavelength setting means of the present invention), an irradiation unit 103, a detector 104, and an interface. (I / F) 105 is provided.
  • Each of the LDs 101 (1) to 101 (N) is a light source that oscillates and emits coherent light, and the oscillation wavelengths ⁇ 1 to ⁇ N are set to become longer in that order.
  • These LDs 101 (1) to 101 (N) are connected to the switching element 102 through optical fibers.
  • the switching element 102 switches the LDs 101 (1) to 101 (N) used as the light source of the illumination light 120 in accordance with an instruction from the workstation 1.
  • the irradiation unit 103 emits the illumination light 120 determined as described above toward the object 110. In an actual measurement, the illumination light 120 passes through a plate, a medium, and the like in addition to the object 110, so these should also be made of a light transmissive material.
  • the detector 104 images the interference fringes generated by the illumination light 120 emitted from the irradiation unit 103 as a hologram.
  • a light wave diffracted by the object 110 is converted into an object wave (arc-shaped line on the right side of the object 110 in the figure), and a light wave (including transmitted light) that is not diffracted is a reference wave (a line on the right side of the object 110).
  • the interference fringes produced by these are recorded.
  • the detector 104 is realized by a CCD image sensor, for example.
  • the actual state of the workstation 1 is a computer, a CPU (Central Processing Unit) 10 which is a central processing unit, a memory 12, a monitor 14 including an LCD (Liquid Crystal Display), an input unit 16 including a keyboard and a mouse, and a memory.
  • the parts 20 are connected to each other.
  • the memory 12 is a volatile storage device such as a RAM (Random Access Memory), and the storage unit 20 is a ROM (Read only Memory), a flash memory, an EPROM (Erasable Programmable ROM), an EEPROM (registered trademark) (Electrically EPROM). ) ⁇ HDD (Hard Disc Drive) ⁇ SSD (Solid State Drive) and other nonvolatile storage devices.
  • the storage unit 20 is provided with an imaging control / data analysis program 21. Each element described later included in the imaging control / data analysis program 21 is a functional unit that is realized by the CPU 10 reading this program into the memory 12 and executing it.
  • the storage unit 20 also stores an OS (Operating System) 29.
  • the workstation 1 includes an interface (I / F) 18 for direct connection with an external device and connection with a network such as a LAN (Local Area Network) with the external apparatus.
  • the I / F 18 is connected to the digital holography device 100 via a network cable NW (or wireless LAN).
  • NW or wireless LAN
  • a plurality of digital holography apparatuses 100 connected to the workstation 1 may be provided.
  • the workstation 1 and the digital holography device 100 may be directly connected by a USB cable or the like.
  • the imaging control / data analysis program 21 controls hologram imaging by the digital holography device 100, reproduces an image of the object 110 by a predetermined calculation process based on the captured hologram, and displays the reproduced image on the monitor 14.
  • an imaging parameter setting unit 31 As in the imaging control / data analysis program 21, an imaging parameter setting unit 31, an imaging instruction unit 32, a hologram acquisition unit 33, a phase information calculation unit 34, an image generation unit 35, a display control unit 36, and a hologram A storage unit 37 is shown.
  • the imaging control / data analysis program 21 is not necessarily a single program.
  • the digital holography apparatus 100 may include some of the elements described above as functions.
  • FIGS. 2 to 4 are flowcharts.
  • the program 21 When the program 21 is activated, the program 21 first prompts the user to input a value of the center wavelength ⁇ mid for hologram imaging.
  • the input of the value of the center wavelength ⁇ mid may allow input of an arbitrary value, but any one of N LDs 101 (1) to 101 (N) provided in advance in the digital holography apparatus 100 may be allowed. It is also possible to make a selection.
  • the imaging control / data analysis program 21 next prompts the input of the value of the conjugate image area magnification X 2 (or magnification X). .
  • the imaging parameter setting unit 31 of the imaging control / data analysis program 21 receives the input center wavelength ⁇ mid.
  • a plurality of LD light sources satisfying the equation (9) are selected from the LDs 101 (1) to 101 (N) (step S13).
  • the imaging target object 110 is a colony of induced pluripotent stem cells such as ES cells and iPS cells.
  • the value of the area ratio X 2 of conjugate image input by the user is described as being 10 ( ⁇ 10 in magnification X). This is because, when the object 110 to be observed is a cell colony, it is sufficient that the area of the conjugate image superimposed on the real image is about 10 times that of the real image in order to obtain a clear reproduction image of the shape. This is determined from the results of preliminary experiments. On the other hand, if the wavelength difference between the longest wavelength ⁇ max and the shortest wavelength ⁇ min is too large, as described above, part of the diffracted light protrudes from the detection surface, which causes a decrease in phase recovery accuracy.
  • the wavelength input by the user may be the longest wavelength ⁇ max or the shortest wavelength ⁇ min instead of the center wavelength ⁇ mid .
  • the center wavelength ⁇ mid is a visible region (about 600 nm) that avoids the low wavelength region in which cell toxicity is a concern and the far infrared region in which heat generation due to light absorption is a concern.
  • this value is not a problem in the following explanation.
  • the switching element 102 switches the light source of the illumination light 120 to the j-th LD 101 (j) (step S105).
  • the irradiation unit 103 emits the illumination light 120 toward the object 110 (step S106).
  • the interference fringes between the object wave diffracted by the object 110 and the reference wave that has not been diffracted are imaged by the detector 104 as a hologram (step S107).
  • the imaged hologram data is transmitted to the I / F 18 of the workstation 1 via the I / F 105 (step S108).
  • the hologram acquisition unit 33 acquires the received hologram data and stores it in the hologram storage unit 37 (step S109).
  • the hologram storage unit 37 stores a plurality of (J2 ⁇ J1 + 1 in this embodiment) hologram data used for phase information calculation processing and object image reproduction, which will be described later, for each measurement object designated by the user. is there.
  • the imaging parameter setting unit 31 When the set hologram data is stored, the imaging parameter setting unit 31 then increments j (step S110). If j does not exceed the maximum value J2 (No in step S111), the process proceeds to step S102. Returning to the front, steps S102 to S111 are executed for the next wavelength ⁇ j .
  • step S111 if j exceeds J2 as a result of the increment in step S110 (Yes in step S111), the determination result is that hologram data corresponding to all the light wavelengths from ⁇ J1 to ⁇ J2 set in step S100 is obtained. It means that they are ready.
  • step S ⁇ b> 15 the phase information calculation unit 34 next executes the phase information calculation process for each hologram data. Details of the processing in this step will be described later with reference to FIG. If phase recovery is performed for each hologram data in step S15, the image generator 35 reproduces the object image based on the hologram data after the phase recovery (step S16). The reproduced object image (hereinafter referred to as “reproduced image”) is displayed on the screen of the monitor 14 by the display control unit 36 (step S17). This is the end of the rough processing by the measurement system.
  • the phase information calculation unit 34 converts each hologram into an amplitude image (step S201). Since the hologram has a distribution of intensity values, it cannot be applied to Fourier transform used in propagation calculation described later. Therefore, each intensity value is converted into an amplitude value in this step. Conversion to an amplitude image is performed by calculating the square root of each pixel value.
  • the initial value of the phase image can be arbitrarily determined. For example, all pixel values may be set to 0, or each pixel value may be set at random.
  • j is an identifier of the LD 101 that is the light source of the illumination light 120 (J1 ⁇ j ⁇ J2)
  • a is a direction value that takes a value of 1 or ⁇ 1
  • n (1 ⁇ n) is calculated as described above. Is the number of repetitions.
  • phase information calculation unit 34 calculates back propagation to the object plane (step S204).
  • the phase information calculation unit 34 determines whether or not the value of j + a falls within the range of J1 to J2 (step S205).
  • the determination result in Step S205 is Yes, and j is incremented (Step S207).
  • the phase information calculation unit 34 inverts the sign of a (step S206), and decrements j (step S207).
  • step S207 If j decreases to 1 by repeating decrement in step S207, j ⁇ 1 becomes 0. Therefore, the determination result in step S205 is No, and the sign of a is reversed again in step S206. Thus, in this flowchart, the positive / negative inversion of a and the increment / decrement of j by it are repeated.
  • the phase information calculation unit 34 updates the phase of the object wave based on ⁇ j (step S208). Specifically, in the complex wavefront on the object plane calculated in step S204, the phase is converted to that of the next wavelength by calculation (the amplitude is not updated). In this way, propagation to the detection surface is calculated in a state where only the phase is converted to the next wavelength (step S209), and the difference (that is, error) between the calculation result and the square root of each intensity value of the measured hologram. ) Is larger than the threshold ⁇ (No in step S210), the phase information calculation unit 34 increments n (step S211) and repeats the above processing.
  • step S210 the phase information calculation unit 34 considers that sufficient phase recovery has been performed, and ends the phase information calculation process.
  • FIG. 5 shows a simulation result of the reproduced object image based on the above method.
  • the center wavelength ⁇ mid (or the longest wavelength ⁇ max or the shortest wavelength ⁇ min ) input by the user from among a number of laser diodes (LD) 101 (1) to 101 (N) provided in advance.
  • the appropriate light sources LD101 (J1) to 101 (J2) are selected according to the magnification ratio X, and they are used by switching them.
  • it is not necessary to set such degrees of freedom, and the object to be photographed (observed) is If it is determined in advance, only a plurality of light sources set in advance according to the equation (9) may be provided in the digital holography apparatus 100, and only these may be used. In this case, the calculation time can be shortened by setting the number of light sources (the number of wavelengths) to 2, and in the case of the inspection of the defect of the cell colony as described above, a sufficiently clear image for practical use. Recovery is also possible.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Holo Graphy (AREA)

Abstract

　光源からコヒーレント光を出射する光出射工程と、前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを検出器が撮像する撮像工程と、前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定工程とを含むディジタルホログラム生成方法において、観察対象となる構造について所定の演算手段により再生される実像及び共役像について、対応する実像に重畳したときにその視認を妨げないとして使用者により設定された共役像の拡大率Xに基き、前記波長設定工程にて設定する複数の波長を、その最短の波長λminと最長の波長λmaxとが、λmax/λmin≧(1/X + 1)を満たすように設定する。

Description

ディジタルホログラフィ装置及びディジタルホログラム生成方法
　本発明は、物体波と参照波との干渉縞を記録したホログラムから物体の再生像を算出するためのディジタルホログラフィ装置及びディジタルホログラム生成方法に関し、特に、複数記録されたホログラムから位相情報を算出する方式のディジタルホログラフィ装置及びディジタルホログラム生成方法に関する。
　ディジタルホログラフィでは、光源からの光が物体表面で反射又は透過してくる物体光と同一光源から直接到達する参照光とがイメージセンサ等の検出面で形成する干渉パターン（ホログラム）に基き、所定の演算処理を実施することで物体の再生像を得る。ホログラムは干渉光の強度分布を表す画像であり、この強度分布は、ともに複素波面である物体波と参照波との振幅差及び位相差に基いて決定される。
　但し、ホログラムには波の虚数成分である位相を記録することができない。そのため、ホログラムに基いて物体の再生像を得る際、検出面に対して実像（所期の再生像）と対称な位置に、該実像と複素共役な振幅を持つ偽像（共役像）が算出される。これら2つの像が重畳すると、実像の視認性が低下し、観察対象の鮮明な像を得られない場合がある。
　この問題を解決するために、これまで様々な手法による物体波の位相回復（位相情報の推定）が試みられてきた。例えば特許文献１には、参照波の位相をシフトさせた複数のホログラムを撮像し、線形演算により位相情報を算出する方法が記載されている。しかし同文献に記載の方法では、位相シフトを高精度に行うための圧電素子や、ミラーやレンズ等の光学系が必要であり、装置構成が複雑になり、装置自体も大型になる。
　一方、こうした複雑な光学系を必要としないインライン型のディジタルホログラフィ装置を用いた位相回復方法が非特許文献１に記載されている。この方法では、光透過性の高い物体（例えば細胞）によって回折した光を物体波、回折することなく該物体を透過した光を参照波とし、物体から検出面までの距離が異なる複数のホログラムを撮像する。その後、撮像した複数のホログラムに基いて物体面と検出面との間で光波の伝搬・逆伝搬を反復的に計算することにより位相情報を算出する。しかしながら、この方法でも検出器を高精度に移動させるための駆動機構（圧電素子等）は依然必要である。
　これに対し近年、光源からの出射光波長を複数段階に変化させ、それぞれの波長で撮像したホログラムに基く光波伝搬の反復計算による位相回復方法が提案されている（非特許文献２）。この方法によれば、複雑な光学系や駆動系を必要としない簡易な装置構成で位相回復が可能である。
特開平10-268740号公報
Yan Zhang, et al., "Whole optical wave field reconstruction from double or multi in-line holograms by phase retrieval algorithm", Optic Express, The Optical Society of America，2003，Vol. 11，No. 24，pp. 3234-3241 Peng Bao, et al., "Lensless phase microscopy using phase retrieval with multiple illumination wavelengths", Applied Optics, The Optical Society of America, 2012，Vol. 51, No. 22, pp. 5486-5494
　ディジタルホログラフィ技術の応用先の一つとされる生体顕微鏡では、これまで細胞内構造或いはそれ以上にごく微小な物体を観察可能とすること、すなわち解像度の向上に主眼が置かれてきた。
　ところで、昨今の再生医療分野では、人工多能性幹細胞（ＥＳ細胞やｉＰＳ細胞等）を用いた研究が盛んである。この分野の研究では多能性を維持しつつ良質な細胞を培養する必要があるため、培養状態を慎重に確認する必要がある。例えば、コロニー内の或る細胞に不良（細胞の扁平化及び培地への貼り付き、並びに特定の細胞種への分化等）が生じると、この不良がコロニー内の他の細胞に伝播して、最終的にはコロニー内の細胞全てが不良となってしまう。従って、観察者は上述したような不良が発生していないかを日々確認し、不良が発生していればこれを取り除く必要がある。このような確認は顕微鏡観察によってなされるが、細胞自体が光学的に視認しにくいことから位相差顕微鏡が用いられるものの、不良の発生を判断することは容易ではなく、観察者にとっては大きな負担となっている。そこで、ディジタルホログラフィ技術を用いてこれを自動化することが考えられている。
　人工多能性幹細胞の培養中に生じる上述のような不良現象は、細胞内構造と比較するとずっとマクロな構造の状態（具体的には厚みを含むコロニーの形状及び質感等）に基いて判断される。ここで本発明者は、非特許文献２に記載の構成でこうしたマクロな構造を観察しようとした場合に、対象領域に上述した共役像が重畳して実像の視認性が低下してしまうことを計算により確認している。すなわち、人工多能性幹細胞の培養における良不良の判定においては、非特許文献２に記載の構成は好適とは言い難い。
　また、非特許文献２に記載の方法を用いても位相の完全な回復は難しく、インライン型のディジタルホログラフィにおいて共役像の重畳は避け難いものである。これは、位相回復の程度を表す明確な指標となる数値が定義されていないことによる。同文献では、複数波長によるホログラム間における推定振幅値の相関係数が近似精度の評価指標として用いられている。しかしながら、例えば複数のホログラムが撮像されたときの照明光波長がほぼ同一であれば、位相回復が全く行われていない状態であっても相関係数の値は大きいため、良好な近似が行われていると判定されることになる。従って、評価指標としての相関係数は反復計算の収束の目安とすることはできるが、位相回復の精度を表す指標として正確ではない。
　本発明は上記の事情に鑑みてなされたものであり、その目的とするところは、目的に応じて大きさの異なる観察対象について鮮明な像の再生を可能とするディジタルホログラフィ装置及びディジタルホログラム生成方法を提供することにある。
　ここで、物体波の位相回復により観察対象の画質を向上させる手法について概要を先に説明する。
　位相情報算出処理は、撮像された複数のホログラムに基き、光波伝搬計算により物体波の位相を推定するアルゴリズムである。図４に一例を示す。各ステップの詳細については後述する。なお、本明細書での「光波伝搬」なる文言には、特に断りがない限り、上述した伝搬と逆伝搬との両方の概念が含まれるものとする。
　光波伝搬の計算手法としてはフレネルの回折計算法や角スペクトル法等が知られており、このうち角スペクトル法による計算式は例えば以下の式（１）で表される。
Figure JPOXMLDOC01-appb-M000001
 　E(x,y,0)は物体面における複素振幅分布、E(x,y,z)は検出面における複素振幅分布であり、zは伝搬距離に相当する。また、kは波数を表す。
　上記の式（１）は、様々なピッチの回折格子の総和としての物体像を回折計算により求めることに相当する。回折格子に入射した光の回折角度φは以下の式（２）で表される。λは波長、pは格子のピッチ幅を表す。
Figure JPOXMLDOC01-appb-M000002
 　上記式（２）より、pが小さい、すなわち回折格子のピッチが細かい（p＝p1）と回折角度Φが大きくなり、波面（物体波）は検出面に対し垂直方向に大きく広がっていく（図６の左側）。一方、pが大きい、すなわち回折格子のピッチが粗い（p＝p2＞p1）と回折角度Φが小さくなり、検出面に対する波面の垂直方向の広がりは僅かである（同図の右側）。この性質は、共役像が検出面及び物体面へと逆伝搬する場合においても同様である。
　なお、本明細書で参照する各図面において、各部材の寸法比は説明の簡略化のため適宜変更されており、実際とは異なる場合がある。
　ここで、回折格子のピッチの細かさは画像中の構造の細かさに相当する。換言すれば、画像中の微細な構造（例えば細胞内構造）はピッチの細かい高周波の格子に、大きな構造（例えばコロニーの形状）はピッチの粗い低周波の格子と看做すことができる。
　図７Ａおよび図７Ｂは、光波伝搬計算により再生された実像及び共役像が重畳する現象を説明する概念図である。同図に示すように、高周波の格子に相当する微細な構造の共役像は物体面への逆伝搬中に大きく広がる（ぼける）ため、該物体面上で重畳する実像の視認性に与える影響が小さい（図７Ａ）。これに対し、低周波の格子に相当する大きな構造の共役像は上述のとおり広がりが小さいため、実像とほぼ同じ大きさを保ったまま物体面に到達する（図７Ｂ）。この場合、物体面において実像の各境界線のごく近傍に類似した線が重畳することになり、実像の視認性が大きく低下する。
　以上説明した原理に基づけば、観察対象となる構造が比較的大きい構造である場合には、低周波の格子によって回折した複素波面（以下「低周波成分」と称す。同様に、高周波の格子によって回折した複素波面を「高周波成分」と称す。）の位相を回復することで、実像の視認性低下を抑制することができる。
　そこで、次に低周波成分の位相を回復するための照明光の波長条件を述べる。図８に示すように、高周波成分と比較すると低周波成分の回折角度は小さいため、検出面H1及びH2間の光路長差が小さい場合（図中d1）、回折による干渉縞のずれ量はごく僅かである（同図の場合、高周波成分に係るずれ量g1に対し、低周波成分のd1でのずれ量は図示が不可能なほど微量である。）。この条件では、低周波成分については検出面H1及びH2で撮像された2つのホログラム間での差異が現れにくく、位相回復が難しい。従って、低周波成分の位相を回復しやすくするためには、同図の検出面H1及びH3のように光路長差を大きくして（図中d2）干渉縞のずれ量を増大させればよい（同d2におけるg2）。
　ここで、これら複数の検出面は、異なる波長の照明光をそれぞれ検出する同一位置の検出面を概念的に示したものと捉えることができ、物体面に近い検出面は短波長に、遠い検出面は長波長に対応している。2つの波長λA及びλB（λA＜λB）を使用した場合に生じる光路長差ΔzABは、以下の式（３）で表される。zは上記式（１）と同様に伝搬距離を示す。
Figure JPOXMLDOC01-appb-M000003
 　上記式（３）によれば、2つの波長λAとλBとの差が大きいほど光路長差ΔzABは大きくなるため、低周波成分の位相回復が容易となる。
　図７Ａに示すように、物体面上に再生された共役像の波面の広がりが大きいほど、対応する実像の視認性への影響が小さくなる。そこで、観察対象となる構造がピッチ幅plimの格子に相当する場合を考える。このとき、物体面において共役像の波面がピッチ幅plimのX倍以上に広がれば、重畳する実像の視認性への影響が十分小さいと仮定したとき、式（２）から下記の式（４）が成り立つ。なお、2zは伝搬距離zの2倍であり、算出された共役像の位置から物体面への距離に相当する。
Figure JPOXMLDOC01-appb-M000004
 　さらに、上記式（４）から以下の式（５）が導かれる。
Figure JPOXMLDOC01-appb-M000005
　すなわち、上記式（５）を満たすピッチ幅plimの格子に相当する成分の位相を回復すればよいことになる。なお、図７Ａおよび図７Ｂを参照して説明したように、細かい構造の共役像は粗い構造よりも大きく広がるため、観察対象であるplimの格子に相当する成分の位相回復がなされれば、plimよりもピッチ幅pが小さい高周波成分の影響は十分小さいと考えられる。
　さらに、2つの波長λAとλBとでホログラムを撮像したとき、ホログラム間の差異が最も顕著になる、すなわちピッチ幅plimに相当する成分の位相回復が容易になるのは、2つの干渉縞がピッチ幅plimのちょうど半分ずれた場合である。従って、式（２）及び式（３）より、以下の式（６）が成り立つ。
Figure JPOXMLDOC01-appb-M000006
 　さらに、上記式（６）から以下の式（７）が導かれ、この式（７）を満たすλA及びλBが、plimに相当する成分の位相回復に適した照明光の波長となる。
Figure JPOXMLDOC01-appb-M000007
 　式（５）及び式（７）より、以下の式（８）が導かれる。
Figure JPOXMLDOC01-appb-M000008
　以上説明した原理に基き、上記課題を解決するために成された第１の発明に係るディジタルホログラフィ装置は、
　a)コヒーレント光を出射する光源と、
　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像するための検出器と、
　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定手段と
　を備え、
　前記波長設定手段は、観察対象となる構造について所定の演算手段により再生される実像及び共役像について、対応する実像に重畳したときにその視認を妨げないとして使用者により設定された共役像の拡大率Xに基き、前記波長設定手段が設定する複数の波長を、その最短の波長λminと最長の波長λmaxとが
　λmax/λmin ≧ (1/X ＋ 1)　…（９）
　を満たすように設定することを特徴とする。
　ここで、「重畳した実像の視認を妨げない」とは、換言すれば、実像中で特に注目される境界線（エッジ）が明瞭であることを意味している。例えば、実像のエッジの近傍領域内で重畳する共役像の形状が該実像と大きく異なっており、両者の別が明確に判断可能である場合には、この共役像は実像の視認を妨げないと看做される。
　上記の構成によれば、光源から出射されたコヒーレント光に起因する物体波と参照波との干渉縞であるホログラムが検出器によって撮像される。ここで、上記複数の波長を定める波長設定手段は、最短の波長λminと最長の波長λmaxとが上記式を満たすよう、当該複数の波長を設定する。
　この関係式は上記の式（８）から導かれるものである。上述のとおり式（８）は、観察対象がピッチ幅plimの回折格子に相当する構造である場合に、実像の視認性を低下させない共役像の拡大率をXとした場合（この共役像の拡大率Xの値は、使用者が観察対象の特性を勘案して上記判断より定め、本発明の装置に与える。）の照明光波長λA及びλB（λA＜λB）の比率を定義している。従って、波長設定手段が定める最短波長λminと最長波長λmaxとの比が上記式（９）を満たせば、観察対象となる構造についての位相回復精度が向上する。また、位相回復されずに残存した共役像が実像に重畳しても、該共役像は実像のX倍以上に拡大されているため、該実像の視認は妨げられない。よって、細胞コロニーのように比較的大きな観察対象についても鮮明な像を得ることができる。なお、位相情報の算出はディジタルホログラフィ装置が備える演算手段が行ってもよいし、ディジタルホログラフィ装置と接続された制御・解析用の外部コンピュータ上で行われてもよい。
　上記「複数の波長」について、何種類の波長を使用するかについてはハードウェア、計算コスト及び求められる画質等の観点から適宜に決定すればよい。使用する波長の数を増やすことにより、計算コストは増大するものの、より広い周波数領域の位相回復が可能となる。
　また、本発明の波長設定手段は光源が備える機能の一つであってもよいし、検出器側で実現されてもよい。前者の典型例は出射光の波長を切り替え可能な光源であり、後者の例としては検出器をハイパースペクトルカメラで実現し、撮像画像から複数の波長の像を個別に抽出する構成が考えられる。さらに別の例として、光ファイバ等で複数の光源と接続されたスイッチング素子を設け、いずれの光源からの出射光を物体に照射するかを切り替える構成とすることもできる。
　上記課題を解決するために成された第２の発明に係るディジタルホログラフィ装置は、
　a)コヒーレント光を出射する光源と、
　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像するための検出器と、
　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定手段と
　を備え、
　前記波長設定手段は、前記複数の波長を、その最短の波長λminと最長の波長λmaxとが
　λmax/λmin ≧ 1.3　…（１０）
　を満たすように設定することを特徴とする。
　上記数値は、共役像の面積が実像の10倍である場合のものである。共役像と実像の面積比がこの程度であれば、一般的には実像の視認性は十分であるとの本発明者によるシミュレーション結果が得られている。このとき、√Xがほぼ0.3となり、式（９）の左辺がほぼ1.3となるため、上記式（１０）が導かれる。
　好ましい構成として、前記波長設定手段は、
　1.3 ≦ λmax/λmin ≦ 2.0　…（１１）
　を満たすよう、前記最短の波長λminと前記最長の波長λmaxとを設定する。
　式（３）に示したように、理論上はλminとλmaxとの波長差が大きいほど低周波成分の位相回復は容易となるが、この効果は一定以上のレベルにおいては高周波成分の位相回復精度とトレードオフとなる。具体的には、図８に示すように、波長差が大きいと、例えば検出面H3に対する図中の高周波成分のように、回折光の一部が検出面からはみ出してしまい、位相回復精度の低下要因となる。そこで、これら2波長の比の妥当な上限値を定めるために以下の基準を採用した。すなわち、λminとλmaxとの間の光路長差Δzが伝搬距離zに等しければ、低周波成分の位相回復は十分になされると考えられる。従って式（３）よりλmax/λmin＝2.0が上限値として妥当である。
　上記課題を解決するために成された第３の発明に係るディジタルホログラム生成方法は、
　a)光源からコヒーレント光を出射する光出射工程と、
　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを検出器が撮像する撮像工程と、
　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定工程と
　を含み、
　観察対象となる構造について所定の演算手段により再生される実像及び共役像について、対応する実像に重畳したときにその視認を妨げないとして使用者により設定された共役像の拡大率Xに基き、前記波長設定工程にて設定する複数の波長を、その最短の波長λminと最長の波長λmaxとが
　λmax/λmin ≧ (1/X + 1)　…（１２）
　を満たすように設定することを特徴とする。
　この場合においても、実際的には、
　1.3 ≦ λmax/λmin ≦ 2.0　…（１３）
としておけば、十分明瞭な実像回復を行うことができる。
　なお、ディジタルホログラフィ装置又はディジタルホログラム生成方法においてこのように、複数の光源の波長を
　1.3 ≦ λmax/λmin ≦ 2.0
とすることは、位相アンラッピング法においても有用なものとなる。
　すなわち、2つの波長λA、λBの位相画像から擬似的に長波長の位相画像を構成する位相アンラッピング法において、擬似波長λABは、
　λAB＝（λA×λB）/（λB－λA）＝λB/（λB/λA－1）
で計算される。この擬似波長λABが大きくなると位相の折り返しを生じることなく位相画像を表示することができるというメリットがあるが、位相の分解能が低下するというデメリットがある。そのため、撮影した物体の位相レンジに合わせて適切な波長λA、λBを選択することが望ましい。図９に合成前の両波長λA、λB（λB/λA）と合成後の波長λABの関係を示すが、λB/λA＝1（λA＝λB）の場合、λABは無限に発散し、λB/λA＝2の場合（λB＝2×λA）、λAB＝λBとなり、合成前の波長と同じになる。
　従って、上記のように複数の波長を
　1.3 ≦ λmax/λmin ≦ 2.0
となるように設定することは、位相アンラッピング法を適用する場合においても好適な条件である。
　本発明によれば、大きさの異なる構造を含む観察対象について、目的に応じた構造について、鮮明な像の再生が可能となる。
　例えば細胞の集合体である細胞コロニー等の観察において、コロニー自体の不良等の比較的大きな構造を撮影（観察）する場合、その使用波長の最大値と最小値の比を1.3以上とすることにより、その目的構造について比較的鮮明な画像を得ることができる。
　また、本発明に係る方法は、位相アンラッピング法を適用する場合においても使用することができる。
本発明の一実施形態に係るディジタルホログラフィ装置を含む測定システムの概略構成を示すブロック図。 図１に示す測定システムが実行する処理の流れを示すフローチャート。 ホログラム撮像処理の流れの一例を示すフローチャート。 位相情報算出処理の流れの一例を示すフローチャート。 図１に示す測定システムで再生される各種サイズの観察対象の物体画像のシミュレーション結果。 回折格子のピッチ幅の違いによる入射光の回折角度の違いを示す説明図。 光波伝搬計算により再生された実像及び共役像が重畳する現象を説明するための概念図。 光波伝搬計算により再生された実像及び共役像が重畳する現象を説明するための概念図。 位相回復に必要な光路長が高周波成分と低周波成分とで異なることを示す説明図。 位相アンラッピングにおける元波長λA、λBと合成波長λABの関係のグラフ。
　図１は本発明の一実施形態に係るディジタルホログラフィ装置を含む測定システムの概略構成を示すブロック図である。この測定システムは、ディジタルホログラフィ装置１００と、該ディジタルホログラフィ装置１００と通信可能に接続されたワークステーション１とを含む。
　〔ディジタルホログラフィ装置１００の構成〕
　ディジタルホログラフィ装置１００は顕微鏡であり、N個のレーザーダイオード（ＬＤ）１０１(1)～１０１(N)、スイッチング素子１０２（本発明の波長設定手段に相当）、照射部１０３、検出器１０４及びインターフェース（Ｉ／Ｆ）１０５を備える。
　ＬＤ１０１(1)～１０１(N)はいずれもコヒーレント光を発振・出射する光源であり、その発振波長λ1～λNは、その順に長くなるように設定されている。これらＬＤ１０１(1)～１０１(N)は光ファイバを介してスイッチング素子１０２に接続されている。
　スイッチング素子１０２は、ワークステーション１からの指示に従い、照明光１２０の光源として使用するＬＤ１０１(1)～１０１(N)を切り替える。
　照射部１０３は、上記のようにして定められた照明光１２０を物体１１０に向けて出射するものである。なお、実際の測定では照明光１２０は物体１１０以外にプレートや培地等も透過するので、これらも光透過性の材料にて構成すべきである。
　検出器１０４は、照射部１０３から出射された照明光１２０によって生じた干渉縞をホログラムとして撮像するものである。このホログラムは、物体１１０によって回折した光波を物体波（同図中において物体１１０の右側の円弧型の線）、回折しなかった光波（透過光を含む）を参照波（物体１１０の右側の線分）とし、これらによって生じた干渉縞を記録したものである。検出器１０４は例えばＣＣＤイメージセンサによって実現される。
　〔ワークステーション１の構成〕
　ワークステーション１の実態はコンピュータであり、中央演算処理装置であるＣＰＵ（Central Processing Unit）１０にメモリ１２、ＬＣＤ（Liquid Crystal Display）等から成るモニタ１４、キーボードやマウス等から成る入力部１６及び記憶部２０が互いに接続されている。このうち、上記のメモリ１２はＲＡＭ（Random Access Memory）等の揮発性記憶装置、記憶部２０はＲＯＭ（Read only Memory）・フラッシュメモリ・ＥＰＲＯＭ（Erasable Programmable ROM）・ＥＥＰＲＯＭ（登録商標）（Electrically EPROM）・ＨＤＤ（Hard Disc Drive）・ＳＳＤ（Solid State Drive）等の不揮発性記憶装置によって構成される。記憶部２０には撮像制御・データ解析プログラム２１が設けられている。撮像制御・データ解析プログラム２１が備える後述の各要素は、ＣＰＵ１０がこのプログラムをメモリ１２に読み出して実行することで実現される機能手段である。記憶部２０にはまた、ＯＳ（Operating System）２９も記憶されている。
　ワークステーション１は、外部装置との直接的な接続や、外部装置等とのＬＡＮ（Local Area Network）などのネットワークを介した接続を司るためのインターフェース（Ｉ／Ｆ）１８を備えており、該Ｉ／Ｆ１８よりネットワークケーブルＮＷ（又は無線ＬＡＮ）を介してディジタルホログラフィ装置１００に接続されている。なお、ワークステーション１と接続されるディジタルホログラフィ装置１００は複数台であっても構わない。また、ワークステーション１とディジタルホログラフィ装置１００とはＵＳＢケーブル等によって直接接続されても構わない。
　撮像制御・データ解析プログラム２１は、ディジタルホログラフィ装置１００によるホログラム撮像を制御するとともに、該撮像されたホログラムに基き、所定の演算処理によって物体１１０の像を再生し、該再生した像をモニタ１４の画面上に画像として表示させるためのアプリケーションソフトウェアである。
　図１においては、撮像制御・データ解析プログラム２１に係るように、撮像パラメータ設定部３１、撮像指示部３２、ホログラム取得部３３、位相情報算出部３４、画像生成部３５、表示制御部３６及びホログラム記憶部３７が示されている。なお、撮像制御・データ解析プログラム２１は必ずしも単体のプログラムである必要はなく、例えば上述した要素の一部をディジタルホログラフィ装置１００が機能として備えていてもよい。
　〔測定システムによる処理の流れ〕
　以下、フローチャートである図２～図４を参照しつつ、本実施形態に係るディジタルホログラフィ装置１００を含む測定システムが実行する処理の流れについて説明を行う。
　まず、図２により、ワークステーション１において撮像制御・データ解析プログラム２１が行う基本的な処理の流れについて説明する。該プログラム２１は、起動すると、最初に使用者に、ホログラム撮像の中心波長λmidの値の入力を促す。この中心波長λmidの値の入力は、任意の値の入力を許してもよいが、本ディジタルホログラフィ装置１００に予め備えられているN個のＬＤ１０１(1)～１０１(N)のうちのいずれかを選択させるという方法でもよい。入力部１６より使用者が中心波長λmidを入力する（ステップＳ１１）と、次に撮像制御・データ解析プログラム２１は共役像の面積倍率X2（又は、拡大率X）の値の入力を促す。これも、任意の値の入力を許してもよいが、予め定めた値をいくつかモニタ１４の画面上に提示し、その中から選択させるようにしてもよい。使用者が入力部１６より面積倍率X2（又は拡大率X）の値を入力する（ステップＳ１２）と、撮像制御・データ解析プログラム２１の撮像パラメータ設定部３１が、入力された中心波長λmidを中心とし、前記式（９）を満たすような複数のＬＤ光源をＬＤ１０１(1)～１０１(N)の中から選択する（ステップＳ１３）。すなわち、ここで選択された光源をＬＤ１０１(J1)～１０１(J2)とする（1 ≦ J1 < J2 ≦ N）と、ＬＤ１０１(J1)の波長λJ1～ＬＤ１０１(J2)の波長λJ2は
　λJ2/λJ1 ≧ 1/X + 1、
　average（λJ1, …, λJ2）≒ λmid　
を満たす。
　本実施形態では、撮像対象の物体１１０は、ＥＳ細胞やｉＰＳ細胞等の人工多能性幹細胞のコロニーとする。そして、使用者の入力した共役像の面積倍率X2の値が10（拡大率Xで√10）であるとして説明する。これは、観察対象となる物体１１０が細胞コロニーであるときにその形状の鮮明な再生像を得るためには、実像と重畳する共役像の面積が実像の10倍程度あれば十分であるとの予備実験の結果から定められたものである。一方、最長波長λmaxと最短波長λminの波長差が大きすぎると、前述の通り回折光の一部が検出面からはみ出してしまい、位相回復精度の低下要因となる。撮像パラメータ設定部３１はこれらを勘案してステップＳ１３において1.3（＝ 1/√10 + 1）≦ λmax /λmin ＝ λJ2/λJ1 ≦ 2.0となるように光源ＬＤ１０１(J1)～１０１(J2)を選択する。なお、ここで使用者の入力する波長は中心波長λmidではなく、最長波長λmax 又は最短波長λminとしてもよい。また、細胞観察の例の場合、中心波長λmidとしては細胞への毒性が懸念される低波長領域や、光の吸収による発熱が懸念される遠赤外領域を避けた可視領域（600nm程度)が適切であるが、以下の説明ではこの値は特に問題とはならない。
　こうして使用する光源が選択された後、ホログラム撮像処理が行われる（ステップＳ１４）。これを図３により詳しく説明する。
　まず、j = J1（J1 ≦ j ≦ J2）とし（ステップＳ１０１）、撮像パラメータ設定部３１は照射部１０３から出射される照明光１２０（図１参照）の波長をλjに設定する（ステップＳ１０２）。続いて撮像指示部３２がディジタルホログラフィ装置１００に対し、λjのホログラム撮像を指示する（ステップＳ１０３）。撮像指示部３２による上記指示は、撮像指示信号としてＩ／Ｆ１８からディジタルホログラフィ装置１００のＩ／Ｆ１０５に対し送信される（ステップＳ１０４）。
　ディジタルホログラフィ装置１００が備えるＩ／Ｆ１０５が上記の撮像指示信号を受信すると、スイッチング素子１０２は照明光１２０の光源をj番目のＬＤ１０１(j)に切り替える（ステップＳ１０５）。次に、照射部１０３が物体１１０に向けて照明光１２０を出射する（ステップＳ１０６）。そして、物体１１０により回折した物体波と回折しなかった参照波との干渉縞はホログラムとして検出器１０４に撮像される（ステップＳ１０７）。撮像されたホログラムデータはＩ／Ｆ１０５を介してワークステーション１のＩ／Ｆ１８に送信される（ステップＳ１０８）。
　ワークステーション１が備えるＩ／Ｆ１８が上記のホログラムデータを受信すると、ホログラム取得部３３が該受信したホログラムデータを取得し、ホログラム記憶部３７に保存する（ステップＳ１０９）。ホログラム記憶部３７は、後述の位相情報算出処理や物体画像の再生に用いる複数（本実施形態ではJ2－J1＋1個）のホログラムデータを、ユーザが指定した測定対象物ごとにまとめて保存するものである。
　設定されたホログラムデータが保存されると、次に撮像パラメータ設定部３１はjをインクリメントし（ステップＳ１１０）、jが最大値J2を超えていなければ（ステップＳ１１１でNo）、処理はステップＳ１０２の手前に戻り、次の波長λjについてステップＳ１０２～Ｓ１１１が実行される。
　一方、ステップＳ１１０におけるインクリメントの結果、jがJ2を上回った場合（ステップＳ１１１でYes）、この判定結果はステップＳ１００で設定されたλJ1～λJ2までの全ての光波長に対応するホログラムデータが揃ったことを意味する。
　図２に戻り、こうしてホログラム撮像処理（ステップＳ１４）が終了すると、次に位相情報算出部３４が各ホログラムデータに対し位相情報算出処理を実行する（ステップＳ１５）。本ステップにおける処理の詳細は図４を参照して後述する。
　ステップＳ１５にて各ホログラムデータについて位相回復がなされれば、該位相回復後のホログラムデータに基いて、画像生成部３５が物体画像を再生する（ステップＳ１６）。再生された物体画像（以下「再生像」と称す）は表示制御部３６によってモニタ１４の画面上に表示される（ステップＳ１７）。以上で測定システムによる大まかな処理は終了となる。
　〔位相情報算出処理の流れ〕
　上述のステップＳ１５における位相情報算出処理の流れの一例を図４に示す。この位相情報算出処理は、ステップＳ１０９にてホログラム記憶部３７に保存されたJ2－J1+1個のホログラムデータに基いて行われる。
　まず、位相情報算出部３４は、各ホログラムを振幅像に変換する（ステップＳ２０１）。ホログラムは強度値の分布であるため、後述の伝搬計算で用いるフーリエ変換に適用することができない。そのため、本ステップで各強度値を振幅値に変換する。振幅像への変換は各画素値の平方根を算出することでなされる。
　次に、位相情報算出部３４は、j = 1、a = 1、n = 1とし、検出面における位相像の初期値を設定する（ステップＳ２０２）。位相像の初期値は任意に定めることができ、例えば全ての画素値を0としてもよいし、各画素値をランダムに設定してもよい。なお、jは前記同様、照明光１２０の光源となるＬＤ１０１の識別子であり（J1 ≦ j ≦ J2）、aは1又は-1の値をとる方向値であり、n（1 ≦ n）は演算の繰り返し回数である。
　次に、位相情報算出部３４はλjの振幅像を更新する（ステップＳ２０３）。具体的にはステップＳ２０１にて強度値からの変換で求めた振幅像を代入する。同図に示す処理ではjの更新と対応付けて「更新」との表現を用いているが、j＝J1の場合は「代入」との表現がより相応しいであろう。
　次に、位相情報算出部３４は物体面への逆伝搬を計算する（ステップＳ２０４）。
　そして位相情報算出部３４は、j＋aの値がJ1～J2の範囲に収まるか否かを判定する（ステップＳ２０５）。ここで、J2の値を仮に5とした場合、最初の試行ではj＋aの値は2であるため、ステップＳ２０５における判定結果はYesとなり、jがインクリメントされる（ステップＳ２０７）。このようにインクリメントが繰り返された結果jの値が5となった場合、換言すればjがJ2に等しくなった場合、j＋1は6となりJ2を上回ることとなるため、ステップＳ２０５における判定結果はNoとなる。このとき、位相情報算出部３４はaの正負を反転させ（ステップＳ２０６）、jをデクリメントする（ステップＳ２０７）。ステップＳ２０７におけるデクリメントの繰り返しによりjの値が1まで減少すればj－1は0となるため、ステップＳ２０５における判定結果はNoとなり、ステップＳ２０６にて再びaの正負が反転する。このように、本フローチャートではaの正負反転とそれによるjのインクリメント及びデクリメントが繰り返される。
　上述したように、j＝1のときステップＳ２０７ではjはインクリメントされて2となる。続いて、位相情報算出部３４はλjによる物体波の位相を更新する（ステップＳ２０８）。具体的には、ステップＳ２０４で計算した物体面における複素波面において、計算により位相を次の波長のものに変換する（振幅は更新しない）。このようにして位相だけを次の波長に変換した状態で、検出面への伝搬を計算し（ステップＳ２０９）、該計算結果と実測値であるホログラムの各強度値の平方根との差分（すなわち誤差）の総和が閾値εよりも大きい場合（ステップＳ２１０でNo）、位相情報算出部３４はnをインクリメントし（ステップＳ２１１）、上記処理を繰り返す。
　一方、誤差の総和が閾値εよりも小さくなっていれば（ステップＳ２１０でYes）、位相情報算出部３４は十分な位相回復がなされたと看做し、位相情報算出処理を終了する。
　上記方法に基づき、再生される物体画像のシミュレーション結果を図５に示す。これらはいずれもUSAFチャートホログラムモデルの、位相回復・逆伝播のシミュレーション結果であり、シミュレーション結果５１は位相回復なしで逆伝播(400nm)した場合、シミュレーション結果５２は400nmと520nm（λmax/λmin = 1.3）で位相回復して逆伝播した場合、シミュレーション結果５３は400nmと800nm（λmax/λmin = 2.0）で位相回復して逆伝播した場合、シミュレーション結果５４は400nmと520nmと800nm（λmax/λmin = 2.0、但し3波長）で位相回復して逆伝播した場合である。λmax/λmin ≧ 1.3とした場合、顕著な位相回復効果が見られる。また、その範囲内で波長の数を増やすことによる効果も見られる。
　〔変更例〕
　上記実施形態では、予め備えられた多数のレーザーダイオード（ＬＤ）１０１(1)～１０１(N)の中から、使用者が入力する中心波長λmid（或いは最長波長λmax又は最短波長λmin）と拡大率Xに応じて適切な光源ＬＤ１０１(J1)～１０１(J2)を選択し、それらを切り替えて使用するとしたが、そのような自由度の設定が必要でなく、撮影（観察）対象が予め定まっている場合には、前記式（９）に従って予め設定した複数の光源のみをディジタルホログラフィ装置１００に設けておき、それらのみを使うようにしてもよい。この場合、光源の数（波長の数）を2とすることにより計算時間を短縮することができ、前記のような細胞コロニーの不良の検査等の目的の場合には、実用上十分鮮明な画像回復も可能となる。
１０…ＣＰＵ
１２…メモリ
１４…モニタ
１６…入力部
１８…Ｉ／Ｆ
２０…記憶部
２１…撮像制御・データ解析プログラム
３１…撮像パラメータ設定部
３２…撮像指示部
３３…ホログラム取得部
３４…位相情報算出部
３５…画像生成部
３６…表示制御部
３７…ホログラム記憶部
１００…ディジタルホログラフィ装置
１０１(1)～１０１(N)…レーザーダイオード（ＬＤ）
１０２…スイッチング素子
１０３…照射部
１０４…検出器
１０５…Ｉ／Ｆ
１１０…物体
１２０…照明光

Claims (9)

  1. 　a)コヒーレント光を出射する光源と、
    　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像するための検出器と、
    　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定手段と
    　を備え、
    　前記波長設定手段は、観察対象となる構造について所定の演算手段により再生される実像及び共役像について、対応する実像に重畳したときにその視認を妨げないとして使用者により設定された共役像の拡大率Xに基き、前記波長設定手段が設定する複数の波長を、その最短の波長λminと最長の波長λmaxとが
    　λmax/λmin ≧ (1/X + 1)
    　を満たすように設定することを特徴とするディジタルホログラフィ装置。
  2. 　a)コヒーレント光を出射する光源と、
    　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像するための検出器と、
    　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定手段と
    　を備え、
    　前記波長設定手段は、前記複数の波長を、その最短の波長λminと最長の波長λmaxとが
    　λmax/λmin ≧ 1.3
    　を満たすように設定することを特徴とするディジタルホログラフィ装置。
  3. 　前記波長設定手段は、前記複数の波長を、その最短の波長λminと最長の波長λmaxとが
    　1.3 ≦ λmax/λmin ≦ 2.0
    　を満たすように設定することを特徴とする請求項２に記載のディジタルホログラフィ装置。
  4. 　前記波長設定手段は、光ファイバを介して複数の前記光源と接続されたスイッチング素子であり、該接続された複数の光源を切り替えて使用することを特徴とする請求項１～３のいずれかに記載のディジタルホログラフィ装置。
  5. 　a)共にコヒーレント光を出射する2つの光源であって、各発光波長λ1及び波長λ2がλ2/λ1≧1.3を満たす光源と、
    　b)前記両光源のそれぞれについて、各出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像するための検出器と、
    　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定手段と
    　を備えることを特徴とするディジタルホログラフィ装置。
  6. 　a)光源からコヒーレント光を出射する光出射工程と、
    　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを検出器が撮像する撮像工程と、
    　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定工程と
    　を含み、
    　観察対象となる構造について所定の演算手段により再生される実像及び共役像について、対応する実像に重畳したときにその視認を妨げないとして使用者により設定された共役像の拡大率Xに基き、前記波長設定工程にて設定する複数の波長を、その最短の波長λminと最長の波長λmaxとが
    　λmax/λmin ≧ (1/X + 1)
    　を満たすように設定することを特徴とするディジタルホログラム生成方法。
  7. 　a)コヒーレント光を出射する光出射工程と、
    　b)前記光源からの出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像するための撮像工程と、
    　c)前記検出器に撮像させるホログラムを生成する照明光の波長を複数設定する波長設定工程と
    　を含み、
    　前記波長設定工程において、前記複数の波長を、その最短の波長λminと最長の波長λmaxとが
    　λmax/λmin ≧ 1.3
    　を満たすように設定することを特徴とするディジタルホログラム生成方法。
  8. 　前記波長設定工程において、前記複数の波長を、その最短の波長λminと最長の波長λmaxとが
    　1.3 ≦ λmax/λmin ≦ 2.0
    　を満たすように設定することを特徴とする請求項７に記載のディジタルホログラム生成方法。
  9. 　共にコヒーレント光を出射する2つの光源であって、各発光波長λ1及び波長λ2がλ2/λ1≧1.3を満たす光源を用意し、
    　前記両光源のそれぞれについて、各出射光に起因する物体波と参照波との干渉縞であるホログラムを撮像する
    　ことを特徴とするディジタルホログラム生成方法。
PCT/JP2015/070019 2014-11-27 2015-07-13 ディジタルホログラフィ装置及びディジタルホログラム生成方法 WO2016084420A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201580064772.4A CN107003638B (zh) 2014-11-27 2015-07-13 数字全息摄像装置以及数字全息图生成方法
US15/529,538 US10281877B2 (en) 2014-11-27 2015-07-13 Digital holography device and digital hologram generation method
EP15863320.6A EP3223083A4 (en) 2014-11-27 2015-07-13 Digital holography device and digital hologram generation method
JP2016561420A JP6388036B2 (ja) 2014-11-27 2015-07-13 ディジタルホログラフィ装置及びディジタルホログラム生成方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014-239506 2014-11-27
JP2014239506 2014-11-27

Publications (1)

Publication Number Publication Date
WO2016084420A1 true WO2016084420A1 (ja) 2016-06-02

Family

ID=56074004

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/070019 WO2016084420A1 (ja) 2014-11-27 2015-07-13 ディジタルホログラフィ装置及びディジタルホログラム生成方法

Country Status (5)

Country Link
US (1) US10281877B2 (ja)
EP (1) EP3223083A4 (ja)
JP (1) JP6388036B2 (ja)
CN (1) CN107003638B (ja)
WO (1) WO2016084420A1 (ja)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108513031A (zh) * 2017-02-28 2018-09-07 株式会社岛津制作所 细胞观察系统
GB2561964A (en) * 2017-02-28 2018-10-31 Shimadzu Corp Cell observation apparatus
WO2018235476A1 (ja) * 2017-06-22 2018-12-27 ソニー株式会社 情報処理装置、情報処理方法、及びプログラム
WO2019171453A1 (ja) * 2018-03-06 2019-09-12 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
WO2019171546A1 (ja) * 2018-03-08 2019-09-12 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
WO2020105534A1 (ja) * 2018-11-22 2020-05-28 株式会社島津製作所 逐次近似計算方法、逐次近似計算装置およびプログラム
WO2020262551A1 (ja) * 2019-06-26 2020-12-30 株式会社島津製作所 細胞機能の評価方法及び細胞解析装置
US11106178B2 (en) 2017-03-03 2021-08-31 Shimadzu Corporation Cell observation device
US11609537B2 (en) 2017-03-02 2023-03-21 Shimadzu Corporation Cell analysis method and cell analysis system using a holographic microscope
US11686721B2 (en) 2018-03-20 2023-06-27 Shimadzu Corporation Cell image analysis apparatus, cell image analysis system, method of generating training data, method of generating trained model, training data generation program, and method of producing training data

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017006370A1 (ja) * 2015-07-07 2017-01-12 オリンパス株式会社 デジタルホログラフィック撮影装置
JP6995975B2 (ja) 2018-03-12 2022-02-04 富士フイルム株式会社 判定方法
CN110895389B (zh) * 2018-09-10 2021-08-24 英属开曼群岛商音飞光电科技股份有限公司 全像影像片体、全像影像记录方法与重建方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010092739A1 (ja) * 2009-02-13 2010-08-19 国立大学法人京都工芸繊維大学 干渉計測装置および干渉計測方法
JP2012531584A (ja) * 2009-06-25 2012-12-10 フェイズ ホログラフィック イメージング ペーホーイー アーベー デジタルホログラフィック撮像による卵又は胚の分析
US20140268105A1 (en) * 2013-03-15 2014-09-18 Zygo Corporation Optical defect inspection system

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3471556B2 (ja) 1997-03-27 2003-12-02 理化学研究所 位相シフトディジタルホログラフィ装置
US7312875B2 (en) * 2003-04-23 2007-12-25 Ut-Battelle Llc Two-wavelength spatial-heterodyne holography
US7978336B2 (en) * 2009-03-16 2011-07-12 Ut-Battelle, Llc Three wavelength quantitative imaging systems
DE102011016660B4 (de) * 2011-04-04 2012-10-25 Universität Stuttgart Verfahren und Anordnung zur Kurz-Kohärenz-Holografie
US20130250383A1 (en) * 2012-03-21 2013-09-26 Coherix, Inc. Apparatus for multi-wavelength holographic imaging
FR2997518B1 (fr) * 2012-10-30 2014-12-05 Centre Nat Rech Scient Systeme d'imagerie holographique auto-reference
US10126709B2 (en) * 2013-12-02 2018-11-13 Imec Vzw Apparatus and method for performing in-line lens-free digital holography of an object
EP2879004A1 (en) * 2013-12-02 2015-06-03 IMEC vzw Apparatus and method for performing in-line lens-free digital holography of an object

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010092739A1 (ja) * 2009-02-13 2010-08-19 国立大学法人京都工芸繊維大学 干渉計測装置および干渉計測方法
JP2012531584A (ja) * 2009-06-25 2012-12-10 フェイズ ホログラフィック イメージング ペーホーイー アーベー デジタルホログラフィック撮像による卵又は胚の分析
US20140268105A1 (en) * 2013-03-15 2014-09-18 Zygo Corporation Optical defect inspection system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MEMMOLO P.: "Multi-wavelength digital holography: reconstruction, synthesis and display of holograms using adaptive transformation", OPTICS LETTERS, vol. 37, no. 9, pages 1445 - 1447, XP001575318, DOI: doi:10.1364/OL.37.001445 *
See also references of EP3223083A4 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2561964A (en) * 2017-02-28 2018-10-31 Shimadzu Corp Cell observation apparatus
CN108513031A (zh) * 2017-02-28 2018-09-07 株式会社岛津制作所 细胞观察系统
CN108513031B (zh) * 2017-02-28 2020-09-25 株式会社岛津制作所 细胞观察系统
US10798333B2 (en) 2017-02-28 2020-10-06 Shimadzu Corporation Cell observation system
US11609537B2 (en) 2017-03-02 2023-03-21 Shimadzu Corporation Cell analysis method and cell analysis system using a holographic microscope
US11106178B2 (en) 2017-03-03 2021-08-31 Shimadzu Corporation Cell observation device
WO2018235476A1 (ja) * 2017-06-22 2018-12-27 ソニー株式会社 情報処理装置、情報処理方法、及びプログラム
JPWO2018235476A1 (ja) * 2017-06-22 2020-04-23 ソニー株式会社 情報処理装置、情報処理方法、及びプログラム
JPWO2019171453A1 (ja) * 2018-03-06 2021-02-04 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
WO2019171453A1 (ja) * 2018-03-06 2019-09-12 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
JP7070656B2 (ja) 2018-03-06 2022-05-18 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
WO2019171546A1 (ja) * 2018-03-08 2019-09-12 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
JPWO2019171546A1 (ja) * 2018-03-08 2021-02-04 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
JP7163955B2 (ja) 2018-03-08 2022-11-01 株式会社島津製作所 細胞画像解析方法、細胞画像解析装置、及び学習モデル作成方法
US11978211B2 (en) 2018-03-08 2024-05-07 Shimadzu Corporation Cellular image analysis method, cellular image analysis device, and learning model creation method
US11686721B2 (en) 2018-03-20 2023-06-27 Shimadzu Corporation Cell image analysis apparatus, cell image analysis system, method of generating training data, method of generating trained model, training data generation program, and method of producing training data
JPWO2020105534A1 (ja) * 2018-11-22 2021-09-30 株式会社島津製作所 逐次近似計算方法、逐次近似計算装置およびプログラム
WO2020105534A1 (ja) * 2018-11-22 2020-05-28 株式会社島津製作所 逐次近似計算方法、逐次近似計算装置およびプログラム
JPWO2020262551A1 (ja) * 2019-06-26 2020-12-30
WO2020261455A1 (ja) * 2019-06-26 2020-12-30 株式会社島津製作所 細胞機能の評価方法及び細胞解析装置
WO2020262551A1 (ja) * 2019-06-26 2020-12-30 株式会社島津製作所 細胞機能の評価方法及び細胞解析装置

Also Published As

Publication number Publication date
US10281877B2 (en) 2019-05-07
CN107003638B (zh) 2019-08-20
JP6388036B2 (ja) 2018-09-12
US20170329281A1 (en) 2017-11-16
JPWO2016084420A1 (ja) 2017-07-06
EP3223083A4 (en) 2017-12-27
EP3223083A1 (en) 2017-09-27
CN107003638A (zh) 2017-08-01

Similar Documents

Publication Publication Date Title
JP6388036B2 (ja) ディジタルホログラフィ装置及びディジタルホログラム生成方法
JP4772961B2 (ja) ディジタル・ホログラムを数値的に再構成することにより、振幅コントラスト画像と定量的位相コントラスト画像を同時に形成する方法
KR101810637B1 (ko) 티코그래피에서 프로브의 보정
EP2592459B1 (en) Holographic microscope and method for processing microscopic subject hologram image
JP5802110B2 (ja) 光変調制御方法、制御プログラム、制御装置、及びレーザ光照射装置
JP4594114B2 (ja) 画像処理装置および屈折率分布測定装置
KR101605178B1 (ko) 다중 기준 영상들을 이용한 ３차원 물체 형상 복원 장치 및 방법
JP2013504780A (ja) 波動場の位相を回復する方法及び装置
KR20160019308A (ko) 홀로그래픽 영상에 대한 품질 측정 장치 및 방법
JP4499331B2 (ja) 波面の位相情報を回復するためのシステムおよび方法
Kocsis et al. Single exposure lensless subpixel phase imaging: optical system design, modelling, and experimental study
Ma et al. Fast algorithm for reliability-guided phase unwrapping in digital holographic microscopy
US11430144B2 (en) Device and process for the contemporary capture of standard images and plenoptic images via correlation plenoptic imaging
Xu et al. Enhanced multiple-plane phase retrieval using a transmission grating
Petrov et al. Image reconstruction using measurements in volume speckle fields formed by different wavelengths
Kuś et al. Advances in design and testing of limited angle optical diffraction tomographysystem for biological applications
Besaga et al. Monitoring of photochemically induced changes in phase-modulating samples with digital holographic microscopy
Huang et al. Infrared Digital Holography
Besaga et al. Performance evaluation of digital holographic microscopy for rapid inspection
WO2022224722A1 (ja) 評価値の取得方法
Li et al. Wavefront sensing and optical surface measurement method based on reference light modulation
Loparev et al. A high-speed modulation interference microscope for biomedical studies
Bae et al. Beam propagation analysis on thickness measurements in quantitative phase microscopy
Stuerwald et al. Imaging micro-optical components with short coherent digital holographic microscopy
Haouat et al. Visualizing the fine structure and dynamics of living cells with temporal polychromatic digital holographic microscopy

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: 15863320

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016561420

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15529538

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

REEP Request for entry into the european phase

Ref document number: 2015863320

Country of ref document: EP