WO2013083563A1 - Methode d'integration de phase holographique differentielle - Google Patents
Methode d'integration de phase holographique differentielle Download PDFInfo
- Publication number
- WO2013083563A1 WO2013083563A1 PCT/EP2012/074362 EP2012074362W WO2013083563A1 WO 2013083563 A1 WO2013083563 A1 WO 2013083563A1 EP 2012074362 W EP2012074362 W EP 2012074362W WO 2013083563 A1 WO2013083563 A1 WO 2013083563A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- areas
- phase
- differential
- phase image
- Prior art date
Links
- 230000010354 integration Effects 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 claims description 7
- 238000009647 digital holographic microscopy Methods 0.000 claims description 2
- 238000000386 microscopy Methods 0.000 abstract description 6
- 230000003287 optical effect Effects 0.000 description 15
- 230000008569 process Effects 0.000 description 15
- 238000012937 correction Methods 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 7
- 238000001093 holography Methods 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 5
- 230000007547 defect Effects 0.000 description 3
- 238000005305 interferometry Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000000399 optical microscopy Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 210000004102 animal cell Anatomy 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
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/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
Definitions
- the present invention relates to a method of integrating a differential holographic representation to determine the optical phase of this representation.
- Digital holography microscopy makes it possible to obtain the quantitative optical phase of an object which constitutes an accurate measurement of its optical thickness in transmission mode or of its relief in reflective mode [1-5]. This feature contributes greatly to the interest that is brought to the microscopy in digital holography because it allows the implementation of an unprecedented quantitative phase contrast technique in conventional optical microscopy with a precision of the order of a few nanometers. 6].
- the instrumentation used to implement this mode of operation is optical interferometry.
- a Mach-Zehnder or Michelson type 2 optical channel interferometer is often used in which one of the channels, called the object channel, comprises the sample and a microscope objective which produces the image of the sample on the sensor. a camera (light beam object).
- the second channel, called reference gives rise to the reference light beam which interferes optically with the object beam.
- the optical phase is obtained by a numerical calculation performed on the recorded interference pattern.
- Placing the sample in one of the arms of the interferometer can change the overall direction of the object beam and can therefore induce misalignment between the object and reference beams. It is therefore often necessary to realign the holographic microscope after the introduction of the sample, and at each sample change.
- Two-channel interferometry is very sensitive to the relative directions of object beams and reference. If the object strongly deflects the beam that illuminates it, the fringe density of the interference pattern increases and one quickly finds oneself in a situation in which the fringes can no longer be resolved by the camera sensor.
- 1 'interferometer is insensitive to misalignment when placing the sample because it is placed at
- the interferometer must not be realigned when the object is changed because it is placed outside the interferometer;
- the interference pattern keeps fringes sufficiently spaced apart on the sensor thus avoiding the problem of resolution at the sensor level.
- the present invention aims to provide a method for calculating the optical phase of a holographic representation from its differential phase.
- the present invention relates to a method for integrating a differential phase image comprising the steps of:
- said predetermined value assigned to the zones undisturbed by the object is a constant, preferably equal to zero.
- the continuous function is determined on an integration line given by linear interpolation of the raw phase image between the two points of intersection of the integration line with the limits of the object. being made when the two points of intersection belong to the image.
- the continuous function is equal, on a given integration line, to the value of the raw phase at the intersection of said integration line with the limit of the object, this equality being applied when a single point of intersection between said integration line and the limit of the object is inside the image.
- the determination of undisturbed areas by an object is performed from the light intensity image.
- predetermined threshold and to 1 the areas where the square module is greater than the predetermined threshold.
- FIG. 1 represents the overlap of the diffracted and undiffracted beams on the sensor for an off-focus object.
- Figure 2 shows the holographic reconstruction of the image of the object and its conjugate.
- Figure 3a shows an intensity image of slightly fuzzy objects.
- Figure 3b shows a detection image of the objects.
- FIG. 3c represents the point-by-point multiplication of the images of FIGS. 3-a and fig. 3-b
- FIG. 4 Configuration for phase correction at a point (s, t).
- FIG. 5 Configuration for the phase correction at a point (s, t) belonging to an area that intercepts an edge of the image
- FIG. 6a represents a differential phase image (animal cells in suspension)
- Figure 6b shows the integrated differential phase of the image of Figure 6a using the proposed method.
- Figure 7a shows the measured phase in a direct manner.
- Figure 7b shows the integrated differential phase without correction process, corresponding to Figure 7a.
- the optical phase of an object is expressed by ⁇
- the index c indicates that it is a signal
- the corresponding differential phase can be expressed by: o Where represents the differential displacement between the two beams.
- interferometric signals are detected on a sensor consisting of a matrix of points photosensitive with a sampling step usually
- the sampled signal can be expressed according to the equation:
- N the number of points per side of the sensor.
- equation (6) The integration method expressed by equation (6) is well suited to the case of integration of the differential phase because it often happens in real situations that the differential step is not constant over the entire field of the sensor. In this case, equation (6) can be used provided that in each point the according to a prior calibration.
- the integration process can be done iteratively to eliminate these errors: after a first iteration, the error is evaluated by differentiating the integrated phase and subtracting the result at the differential phase of departure. The result of this subtraction is reintegrated using equation (6) to calculate a correction of the integrated phase initially obtained. This correction process can be performed several times until an integrated phase is reached below a previously defined error threshold.
- the integration process can not however apply alone because it requires knowledge of the integrated phase on the upper and left edges of the sensor. This is classically the problem of determining the constants of integration. Generally, in the prior art, these integration constants are set arbitrarily equal to zero, which produces artifacts near these edges, or lines that do not correspond to reality. These integration constants can not be known a priori and must be established through additional information. For this, the present invention uses the light intensity information as described in the following section. Determination of integration constants thanks to light intensity
- a differential holographic microscope provides both the differential phase information and the light intensity information of the sample. This intensity information makes it possible to determine the value of the optical phase in certain zones of the image and makes it possible to fix in most cases the local integration constants.
- the configuration shown in Figure 1 is similar to that used in the Gabor holography, namely the first type of optical holography that was invented.
- the Gabor holography allows a holographic reconstruction of the object as well as its so-called conjugated image (Fig.2). We do not get the image phase as obtained with an interferometer configuration.
- the interference pattern recorded between the beam diffracted by the image of the object and the non-diffracted beam contains information on the image of the object and its conjugated image and more particularly on the spatial frequency components of the net image of the object.
- the spatial frequency component (u, v) is defined by the spatial Fourier transform of light intensity:
- Figure 4 shows an object presence area for which the integrated phase without correction in one (s, t) was calculated.
- di is the Euclidean distance between the point (s, t) and the point (si, ti)
- d ⁇ is the Euclidean distance between the point (s, t) and the point (so, to).
- the corrected ⁇ can be calculated according to the following equation:
- the corrected phase can be calculated according to the following equation:
- the method described is effective if there are still areas outside the object that define the integration constants.
- the following figures show an example, corresponding to the image of FIG. 3-a, differential phase integration using the method described.
- FIG. 6-a is the differential image of departure with a horizontal differentiation axis. From this image, and thanks to the image of FIG. 3-b, the integrated phase is calculated to obtain the image of FIG. 6-b. This image is compared to the image of FIG. 7-a which represents the directly recorded non-differential phase. We observe that these two images are very similar. Between these two images, there is a slight overall offset in gray level. This offset has no physical meaning and therefore does not matter.
- FIG. 7-b shows a digital integration process of the direct differential phase without using the integration zones and the correction process described in the present application.
- the initial phases on the upper and left edges were taken equal to zero. This figure shows numerous integration errors that validate the utility of the proposed method.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Holo Graphy (AREA)
Abstract
La présente invention se rapporte à une méthode d'intégration d'une image de phase différentielle comprenant les étapes de : • - obtenir une image de phase différentielle utilisant la microscopie holographique digitale différentielle, pour laquelle un objet à analyser est placé avant l'interféromètre; • - obtenir une image d'intensité lumineuse de l'objet; • - à partir de ladite image d'intensité lumineuse, déterminer des zones de l'image de phase différentielle n'étant pas perturbées par un objet; • - attribuer une valeur de phase prédéterminée auxdites zones non perturbées par un objet; • - intégrer numériquement l'image de phase différentielle à partir des bords de l'image en considérant une constante d'intégration arbitraire pour ces bords pour obtenir une image de phase brute; • - corriger l'image de phase brute par soustraction d'une fonction continue telle que la valeur de l'image de phase corrigée dans les zones sans objet prenne la valeur de phase prédéterminée.
Description
METHODE D'INTEGRATION DE PHASE HOLOGRAPHIQUE DIFFERENTIELLE
Objet de l'invention
[0001] La présente invention se rapporte à une méthode d'intégration d'une représentation holographique différentielle afin de déterminer la phase optique de cette représentation . Etat de la technique
[0002] La microscopie en holographie digitale permet d'obtenir la phase optique quantitative d'un objet qui constitue une mesure précise de son épaisseur optique en mode de transmission ou de son relief en mode réflectif [1- 5] . Cette fonctionnalité contribue fortement à l'intérêt qui est porté à la microscopie en holographie digitale car elle permet 1 ' implémentation d'une technique de contraste de phase quantitative inédite en microscopie optique classique avec une précision de l'ordre de quelques nanomètres [5-6].
[0003] L' instrumentation utilisée pour implémenter ce mode de fonctionnement est 1 ' interférométrie optique. En transmission, on utilise souvent un interféromètre à 2 canaux optiques de type Mach-Zehnder ou Michelson dans lequel un des canaux, dit canal objet, comporte l'échantillon et un objectif de microscope qui réalise l'image de l'échantillon sur le senseur d'une caméra (faisceau lumineux objet) . Le deuxième canal, dit de référence, donne lieu au faisceau lumineux de référence qui
interfère optiquement avec le faisceau objet. La phase optique est obtenue par un calcul numérique effectué sur la figure d'interférence enregistrée.
[0004] En holographie digitale, on enregistre également l'intensité lumineuse de l'objet. Cette information couplée à la phase optique permet de reconstituer l'amplitude complexe provenant de l'échantillon. Grâce à l'amplitude complexe, il devient possible, par le procédé de reconstruction holographique, de pouvoir remettre nette l'image d'un objet qui a été enregistrée floue [2-6] .
[0005] L' implémentation de ce type de microscopie avec un interféromètre à deux canaux présente toutefois des inconvénients :
L'utilisation de deux faisceaux séparés nécessite 1 ' implémentation de sources lumineuses optiques de cohérence optique élevée. Même si certaines configurations permettent l'utilisation des sources de cohérence partielle [4,6], elles sont toujours plus cohérentes que celles utilisées en microscopie optique classique. Or, plus la cohérence de la source est importante, plus le bruit dit de cohérence augmente et on obtient donc généralement des images plus bruitées (speckle) qu'en microscopie classique.
Le fait de placer l'échantillon dans un des bras de 1 ' interféromètre peut changer la direction globale du faisceau objet et peut donc induire un désalignement entre les faisceaux objet et de référence. Il est donc souvent nécessaire de réaligner le microscope holographique après l'introduction de l'échantillon, et lors de chaque changement d'échantillon.
L' interférométrie à deux canaux est très sensible aux directions relatives des faisceaux objet et de
référence. Si l'objet dévie fortement le faisceau qui l'illumine, la densité de franges de la figure d' interférence augmente et on se retrouve rapidement dans une situation dans laquelle les franges ne peuvent plus être résolues par le senseur de la caméra .
[0006] Pour éviter ces inconvénients, on peut utiliser la microscopie holographique digitale dite différentielle [7,8]. Avec ce type d'instrument, l'objet est placé sur le trajet optique avant 1 ' interféromètre . En transmission, le faisceau lumineux provenant de l'objet passe au travers d'un interféromètre à deux bras qui introduit un petit décalage spatial entre les deux faisceaux lorsqu'ils sont incidents sur le senseur. La figure d'interférence des deux faisceaux permet alors de calculer la phase différentielle suivant la direction du décalage spatial entre les deux faisceaux lumineux. Le document EP 1631788 décrit ce type de dispositifs.
[0007] Cette configuration évite tous les inconvénients de la microscopie holographique tels que décrits ci-dessus :
des sources de cohérence très réduite peuvent être utilisées ;
1 ' interféromètre est peu sensible au désalignement quand on place l'échantillon car celui-ci est placé à
1 ' extérieur ;
il ne faut pas réaligner 1 ' interféromètre quand on change l'objet car il est placé à l'extérieur de 1 ' interféromètre ;
- en mesurant la phase différentielle, même quand l'objet introduit une déviation importante de la direction du faisceau lumineux, la figure d' interférence garde des franges suffisamment espacées
sur le senseur évitant ainsi le problème de résolution au niveau du senseur.
La configuration en holographie digitale différentielle est donc très intéressante pour de nombreuses applications.
[0008] Toutefois, pour des évaluations quantitatives des échantillons, il est souvent plus avantageux de mesurer la phase plutôt que la phase différentielle. Il y a donc un processus de calcul à effectuer qui permet d'obtenir la phase directe à partir de la phase différentielle.
[0009] Néanmoins, ce processus d'intégration nécessite la connaissance des conditions aux limites (constantes d'intégrations) à partir desquelles on intègre de proche en proche. Généralement, dans l'art antérieur, on considère arbitrairement une valeur constante sur deux bords de l'image, à partir desquels on intègre numériquement la phase différentielle. Cette fixation arbitraire induit alors des erreurs se propageant sur l'ensemble de l'image, en particulier lorsque des objets se trouvent sur ces bords.
Buts de 1 ' invention
[0010] La présente invention vise à fournir une méthode de calcul de la phase optique d'une représentation holographique à partir de sa phase différentielle.
Résumé de 1 ' invention
[0011] La présente invention se rapporte à une méthode d'intégration d'une image de phase différentielle comprenant les étapes de :
- déterminer des zones de l'image de phase différentielle n'étant perturbées par un objet ;
attribuer une valeur prédéterminée auxdites zones non perturbées par un objet;
intégrer l'image de phase différentielle à partir des bords de l'image en considérant une constante d' intégration arbitraire pour ces bords pour obtenir une image de phase brute ;
- corriger l'image de phase brute par soustraction d'une fonction continue telle que la valeur de l'image de phase corrigée dans les zones sans objet prenne la valeur prédéterminée.
[0012] De préférence, ladite valeur prédéterminée attribuée aux zones non perturbées par l'objet est une constante, de préférence égale à zéro.
[0013] Avantageusement, la fonction continue est déterminée sur une droite d' intégration donnée par interpolation linéaire de l'image de phase brute entre les deux points d' intersection de la droite d' intégration avec les limites de l'objet, cette détermination étant faite lorsque les deux points d' intersection appartiennent à 1 ' image .
[0014] De préférence, la fonction continue est égale, sur une droite d'intégration donnée, à la valeur de la phase brute à l'intersection de ladite droite d'intégration avec la limite de l'objet, cette égalité étant appliquée lorsqu'un seul point d'intersection entre ladite droite d'intégration et la limite de l'objet se trouve à l'intérieur de l'image.
[0015] Avantageusement, la détermination des zones non perturbées par un objet est effectuée à partir de l'image d'intensité lumineuse.
[0016] Avantageusement, la détermination des zones non perturbées par l'objet comprend les étapes suivantes : calcul de la transformée de Fourier de l'image d'intensité lumineuse I(u,v);
filtrage passe haut de I(u,v) par une fonction 1- g(u,v), pour obtenir la transformée filtrée I' (u,v), g(u,v) est une fonction gaussienne centrée sur (u,v) = (0,0) ;
inverse ;
filtrage passe bas du module carré ;
seuil prédéterminé, et à 1 les zones ou le module carré est supérieur au seuil prédéterminé.
Description des figures
[0017] La figure 1 représente le recouvrement des faisceaux diffractés et non diffractés sur le senseur pour un objet hors focus .
[0018] La figure 2 représente la reconstruction holographique de l'image de l'objet et de son conjugué.
[0019] La figure 3a représente une image d' intensité d'objets légèrement flous.
[0020] La figure 3b représente une image de détection des objets.
[0021] La figure 3c représente la multiplication point par point des images des fig. 3-a et fig. 3-b
[0022] Fig. 4 : Configuration pour la correction de phase en un point (s,t) .
[0023] Fig. 5 : Configuration pour la correction de phase en un point (s,t) appartenant à une zone qui intercepte un bord de l'image
[0024] La figure 6a représente une image de phase différentielle (cellules animales en suspension)
[0025] La figure 6b représente la phase différentielle intégrée de l'image de la figure 6a utilisant la méthode proposée.
[0026] La figure 7a représente la phase mesurée de façon directe.
[0027] La figure 7b représente la phase différentielle intégrée sans processus de correction, correspondant à la figure 7a. Description détaillée de 1 ' invention
[0028] Pour obtenir la phase optique à partir de la phase différentielle, il faut effectuer une opération numérique d'intégration. Une méthode adaptée à cette situation est décrite ci-après.
[0029] On suppose que la phase différentielle initiale ait été obtenue par un processus de calcul adéquat de manière à éliminer les sauts de phase. Un exemple d'une telle procédure adéquate est décrit par Vyacheslav V. Volkov and Yimei Zhu, dans le document "Deterministic phase unwrapping in the présence of noise", OPTICS LETTERS, 28, 2156-2158 (2003) que nous incorporons ici par référence.
continu suivant les axes x et y du senseur 2D. La phase différentielle correspondante peut s'exprimer par :
o Ùù
représente le déplacement différentiel entre les deux faisceaux.
[0031] Comme les signaux interférométriques sont détectés sur un senseur constitué d'une matrice de points
photosensibles avec un pas d' échantillonnage généralement
Avec s, t=-N/2,..., (N/2)-l, où N est le nombre de points par côté du senseur.
[0032] En supposant que les variations des fonctions sont lentes à l'échelle de la dimension caractéristique de 1 ' échantillon Δ , il vient :
[0033] L'équation (6) permet donc de calculer la phase
connaissant
et la phase différentielle
. Connaissant la phase sur le bord supérieur (s =-N/2) et sur le bord gauche (t=-N/2) du senseur, nous pouvons donc par un balayage ligne par ligne calculer la phase intégrée en tout point du senseur.
[0034] La méthode d'intégration exprimée par l'équation (6) convient bien au cas de l'intégration de la phase différentielle car il arrive souvent dans des situations réelles que le pas différentiel
ne soit pas constant sur tout le champ du senseur. Dans ce cas, l'équation (6) peut être utilisée à condition que l'on adapte en chaque point les
en fonction d'une calibration préalable.
[0035] Dans le raisonnement qui précède l'hypothèse de la variation lente de la phase différentielle entre des points d'échantillonnage voisins a été faite. On peut donc s'attendre à obtenir des erreurs d'intégration qui se cumulent .
[0036] Toutefois, le processus d'intégration peut être fait de manière itérative pour éliminer ces erreurs : après une première itération, on évalue l'erreur en différenciant la phase intégrée et en soustrayant le résultat à la phase différentielle de départ. Le résultat de cette soustraction est réintégré en utilisant l'équation (6) pour calculer une correction de la phase intégrée initialement obtenue. Ce processus de correction peut être effectué plusieurs fois jusqu'à obtenir une phase intégrée en deçà d'un seuil d'erreur préalablement défini.
[0037] Le processus d'intégration ne peut toutefois pas s'appliquer seul car il nécessite la connaissance de la phase intégrée sur les bords supérieur et gauche du senseur. Il s'agit classiquement du problème de la détermination des constantes d'intégration. Généralement, dans l'art antérieur, ces constantes d'intégration sont fixées arbitrairement égales à zéro, ce qui produit des artéfacts à proximité de ces bords, ou des lignes ne correspondant pas à la réalité. Ces constantes d' intégration ne peuvent donc pas être connues a priori et il faut les établir grâce à une information supplémentaire. Pour cela, la présente invention utilise l'information d'intensité lumineuse comme décrit à la section suivante. Détermination des constantes d'intégration grâce à l'intensité lumineuse
[0038] Un microscope holographique différentiel fournit à la fois l'information de phase différentielle et l'information d'intensité lumineuse de l'échantillon. Cette information d' intensité permet de déterminer la valeur de la phase optique en certaines zones de l'image et permet de fixer dans la plupart des cas les constantes d' intégration locales .
[0039] Dans le cas de l'enregistrement d'un hologramme différentiel d'un objet légèrement flou en transmission, l'intensité, et la figure d'interférence entre l'objet qui diffracte la lumière et le faisceau d'illumination de l'échantillon (Fig .1) sont enregistrés.
[0040] La configuration représentée par la figure 1 est similaire à celle utilisée dans l'holographie de Gabor, à savoir le premier type d'holographie optique qui a été inventé. Or, l'holographie de Gabor permet une reconstruction holographique de l'objet ainsi que son image dite conjuguée (Fig.2). On n'obtient toutefois pas l'image
de phase comme on l'obtient avec une configuration en interféromètre .
[0041] Cela signifie que la figure d'interférence enregistrée entre le faisceau diffracté par l'image de l'objet et le faisceau non diffracté contient l'information sur l'image de l'objet et son image conjuguée et plus particulièrement sur les composantes en fréquences spatiales de l'image nette de l'objet.
[0042] On peut également rencontrer des situations similaires lorsque l'objet est illuminé en réflexion. Par exemple lorsque l'échantillon est une plaque plane réfléchissante sur laquelle il y a un défaut localisé. Dans cette situation, le faisceau réfléchi par la partie plane de l'objet interférera avec le faisceau réfléchi par la zone du défaut. Dans ce qui suit, on considérera comme objet ce qui produit une perturbation non prédictible du faisceau (i.e. un défaut localisé, une particule en suspension dans un milieu transparent) tandis que les zones produisant des modifications prédictibles du faisceau seront considérées comme des zones « sans objet » (i.e. un support plan, un milieu transparent homogène dans un porte échantillon) .
[0043] Dans la configuration en transmission (fig.
1), l'intensité lumineuse enregistrée sur le senseur est
. La composante en fréquences spatiales (u,v) est définie par la transformée de Fourier spatiale de l'intensité lumineuse :
[0044] On montre que s'il n' y a pas d'objet, est significativement différent de zéro que lorsque (u,v)
est proche de (0,0) . Par contre, s'il y a la présence de l'image d'un objet proche du senseur, il y aura des composantes significativement différentes de zéro pour une gamme étendue des valeurs de u et v.
[0045] Par conséquent, en opérant un filtrage passe haut sur la transformée de Fourier, il est possible de déterminer les zones d' intensité lumineuse influencées par la présence de l'image de l'objet. [0046] L'algorithme qui permet de déterminer ces zones est le suivant :
2/ Filtrage passe haut de
par une fonction (l-g(u,v)), ou g{u,v) est une fonction gaussienne réelle centrée sur (u, v) = (0,0) .
3/ Calcul de la transformée de Fourier inverse
4/ Calcul du module carré
5/ Filtrage passe bas par une convolution avec une fonction gaussienne .
6/ Application d'une fonction seuil qui met à 1 les valeurs supérieures au seuil et à 0 les valeurs inférieures au seuil .
Ce processus s'avère très efficace pour sélectionner les zones de présence d'objets comme illustré par la figure 3.
[0047] Les zones de présence d'images d'objets étant déterminées (zones blanches sur Fig. 3-b) , on dispose d'un moyen pour fixer les constantes d' intégration dans le processus décrit précédemment.
[0048] Partout où il n'y a pas d'objet (zone noire), on considérera que la phase intégrée doit être constante. Sans nuire à la généralité, on supposera que cette constante est égale à zéro. Le processus d'intégration
précédant peut alors être appliqué pour tout objet qui est entouré par une zone sans objet (région noire) .
[0049] On applique également une correction supplémentaire directionnelle en chaque point d'une zone signalant la présence d'un objet en considérant que la
Phase intégrée suivant la direction
doit être nulle quand on entre et on sort de la zone, la figure 4 montre une zone de présence d'objet pour laquelle la phase intégrée sans correction
en un (s,t) a été calculée. On trace par le point (s,t) où on souhaite corriger la phase la droite D de direction
. On recherche les points frontières de la zone (si,ti) et (so,to) qui appartiennent à la droite D de direction et qui
passe par le point (s,t) .
[0050] Si les phases intégrées aux points (si,ti) et
Où di est la distance euclidienne entre le point (s,t) et le point (si,ti), et où d∑ est la distance euclidienne entre le point (s,t) et le point (so,to).
On constate alors que cette correction annule bien
[0051] Lorsque une zone de présence d'un objet intercepte un bord de l'image et qu'il n'est pas possible de trouver les deux points (si,ti) ou (so,to), la correction est effectuée sur le seul point (si,ti) ou (so,to) qui se trouve dans l'image (fig. 5) .
[0052] Si (si,ti) appartient à l'image, la ρηε corrigée peut être calculée selon l'équation suivante:
La méthode décrite s'avère efficace s'il subsiste des zones hors objet qui permettent de définir les constantes d' intégration. [0053] Les figures suivantes montrent un exemple, correspondant à l'image de la Fig. 3-a, d'intégration de phase différentielle en utilisant la méthode décrite.
[0054] La Fig. 6-a est l'image différentielle de départ avec un axe de différenciation horizontal. A partir de cette image, et grâce à l'image de la Fig. 3-b, on calcule la phase intégrée pour obtenir l'image de la Fig. 6-b. Cette image est à comparer à l'image de la Fig. 7-a qui représente la phase non différentielle directement enregistrée. On observe que ces deux images sont très similaires. Entre ces deux images, il y a un léger offset global en niveau de gris. Cet offset n'a aucune signification physique et n'a donc aucune importance.
[0055] Les zones abruptes noires qui apparaissent dans certaines zones claires sur les deux images des Fig. 6-b et Fig. 7-a sont dues à de sauts de phase qui proviennent de la périodicité inhérente à
1 ' interférométrie . Ces sauts sont normaux et peuvent être
éliminés en utilisant un processus de déroulage de phase usuel .
[0056] La Fig. 7-b montre un processus d'intégration numérique de la phase différentielle directe sans utiliser les zones d' intégration et le processus de correction décrits dans la présente demande. Pour l'obtention de la phase sans processus de correction (Fig. 6-b) , les phases initiales sur les bords supérieur et gauche ont été prises égales à zéro. On observe sur cette figure de nombreuses erreurs d'intégration qui valident l'utilité de la méthode proposée .
Bibliographie
1. U. Schnars and W. Jtiptner, "Direct recording of holograms by a CCD target and numerical reconstruction",
Appl. Opt. 33, 179-181 (1994).
2. I.Yamaguchi and T. Zhang, " Ph a s e-shifting digital holography", Opt. Lett. 22, 1268-1270 (1997).
3. T. Zhang and I. Yamaguchi, "Three-dimensional microscopy with phase-shi fting digital holography", Opt. Lett. 23,
1221-1223 (1998) .
4. F. Dubois, L. Joannes, J.-C. Legros, "Improved three- dimensional imaging with digital holography microscope using a partial spatial cohérent source " Appl. Opt. 38, 7085-7094 (1999)
5. E.Cuche, F.Bevilacqua and C. Depeursinge, "Digital holography for quantitative phase contrast imaging", Opt. Lett. 24, 291-293 (1999)
6. Frank Dubois, Catherine Yourassowsky, Olivier Monnom, Jean-Claude Legros, Olivier Debeir, Philippe Van Ham,
Robert Kiss, Christine Decaestecker, "Digital holographie microscopy for the three-dimensional dynamic analysis of in vitro cancer cell migration", Journal of Biomédical Optics 11-5, 054032 September/October 2006
7. F. Dubois, C. Yourassowsky, « Digital holographie microscope for 3D Imaging and procès using it », Patent EP1631788B1, 2007.
8. F. Dubois, C. Yourassowsky, "Off Axis interferometer", Patent WO 2011/042442 Al, 2011.
Claims
1. Méthode d'intégration d'une image de phase différentielle comprenant les étapes de :
obtenir une image de phase différentielle utilisant la microscopie holographique digitale différentielle, pour laquelle un objet à analyser est placé avant 1 ' interféromètre ;
obtenir une image d'intensité lumineuse de l'objet; déterminer des zones de l' image de p hase différentielle n'étant perturbées par un objet ;
attribuer une valeur prédéterminée auxdites zones non perturbées par un objet;
intégrer l'image de phase différentielle à partir des bords de l'image en considérant une constante d' intégration arbitraire pour ces bords pour obtenir une image de phase brute ;
corriger l'image de phase brute par soustraction d'une fonction continue telle que la valeur de l'image de phase corrigée dans les zones sans objet prenne la valeur prédéterminée.
2. Méthode de la revendication 1 dans laquelle la ladite valeur prédéterminée attribuée aux zones non perturbées par l'objet est une constante, de préférence égale à zéro.
3. Méthode des revendications 1 ou 2 dans laquelle la fonction continue est déterminée sur une droite d'intégration donnée par interpolation linéaire de l'image de phase brute entre les deux points d' intersection de la droite d'intégration avec les limites de l'objet, cette détermination étant faite lorsque les deux points d'intersection appartiennent à l'image.
4. Méthode de l'une quelconque des revendications précédentes dans laquelle la fonction continue est égale, sur une droite d'intégration donnée, à la valeur de la phase brute à l'intersection de ladite droite d'intégration avec la limite de l'objet, cette égalité étant appliquée lorsqu'un seul point d'intersection entre ladite droite d'intégration et la limite de l'objet se trouve à l'intérieur de l'image.
5. Méthode selon l'une quelconque des revendications précédentes, dans laquelle la détermination des zones non perturbées par un objet est effectuée à partir de l'image d'intensité lumineuse.
6. Méthode de la revendication 5 dans laquelle la détermination des zones non perturbées par l'objet comprend les étapes suivantes :
calcul de la transformée de Fourier de l' image d'intensité lumineuse I(u,v);
filtrage passe haut de I(u,v) par une fonction 1- g(u,v), pour obtenir la transformée filtrée I' (u,v), g(u,v) est une fonction gaussienne centrée sur (u,v) = (0,0) ;
inverse ;
filtrage passe bas du module carré ;
seuil prédéterminé, et à 1 les zones ou le module carré est supérieur au seuil prédéterminé.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP11192638.2 | 2011-12-08 | ||
EP11192638 | 2011-12-08 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013083563A1 true WO2013083563A1 (fr) | 2013-06-13 |
Family
ID=47471707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/EP2012/074362 WO2013083563A1 (fr) | 2011-12-08 | 2012-12-04 | Methode d'integration de phase holographique differentielle |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2013083563A1 (fr) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1524491A1 (fr) * | 2003-10-16 | 2005-04-20 | Universite Libre De Bruxelles | Appareil associant un interféromètre et un microscope |
EP1631788A1 (fr) | 2003-05-16 | 2006-03-08 | Universite Libre De Bruxelles | Microscope holographique numerique pour imagerie 3d et procede l'utilisant |
US20080309944A1 (en) * | 2006-04-21 | 2008-12-18 | Consiglio Nazionale Delle Ricerche | Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus |
WO2011042442A1 (fr) | 2009-10-08 | 2011-04-14 | Universite Libre De Bruxelles | Interféromètre indirect |
-
2012
- 2012-12-04 WO PCT/EP2012/074362 patent/WO2013083563A1/fr active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1631788A1 (fr) | 2003-05-16 | 2006-03-08 | Universite Libre De Bruxelles | Microscope holographique numerique pour imagerie 3d et procede l'utilisant |
EP1524491A1 (fr) * | 2003-10-16 | 2005-04-20 | Universite Libre De Bruxelles | Appareil associant un interféromètre et un microscope |
US20080309944A1 (en) * | 2006-04-21 | 2008-12-18 | Consiglio Nazionale Delle Ricerche | Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus |
WO2011042442A1 (fr) | 2009-10-08 | 2011-04-14 | Universite Libre De Bruxelles | Interféromètre indirect |
Non-Patent Citations (8)
Title |
---|
E.CUCHE; F.BEVILACQUA; C. DEPEURSINGE: "Digital holography for quantitative phase contrast imaging", OPT. LETT., vol. 24, 1999, pages 291 - 293, XP000823520 |
F. DUBOIS; L. JOANNES; J.-C. LEGROS: "Improved three-dimensional imaging with digital holography microscope using a partial spatial coherent source", APPL. OPT., vol. 38, 1999, pages 7085 - 7094 |
FRANK DUBOIS; CATHERINE YOURASSOWSKY; OLIVIER MONNOM; JEAN-CLAUDE LEGROS; OLIVIER DEBEIR; PHILIPPE VAN HAM; ROBERT KISS; CHRISTINE: "Digital holographic microscopy for the three-dimensional dynamic analysis of in vitro cancer cell migration", JOURNAL OF BIOMÉDICAL OPTICS, September 2006 (2006-09-01), pages 11 - 5 |
I.YAMAGUCHI; T. ZHANG: "Phase-shifting digital holography", OPT. LETT., vol. 22, 1997, pages 1268 - 1270 |
PRETZIER G ET AL: "High-accuracy differential interferometry for the investigation of phase objects", MEASUREMENT SCIENCE AND TECHNOLOGY, IOP, BRISTOL, GB, vol. 4, no. 6, 1 June 1993 (1993-06-01), pages 649 - 658, XP020065466, ISSN: 0957-0233, DOI: 10.1088/0957-0233/4/6/003 * |
T. ZHANG; I. YAMAGUCHI: "Three-dimensional microscopy with phase-shifting digital holography", OPT. LETT., vol. 23, 1998, pages 1221 - 1223, XP000783068 |
U. SCHNARS; W. JÜPTNER: "Direct recording of holograms by a CCD target and numerical reconstruction", APPL. OPT., vol. 33, 1994, pages 179 - 181, XP000353611, DOI: doi:10.1364/AO.33.000179 |
VYACHESLAV V. VOLKOV; YIMEI ZHU: "Deterministic phase unwrapping in the presence of noise", OPTICS LETTERS, vol. 28, 2003, pages 2156 - 2158 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zuo et al. | Transport of intensity equation: a tutorial | |
JP5364203B2 (ja) | 観察装置 | |
JP4772961B2 (ja) | ディジタル・ホログラムを数値的に再構成することにより、振幅コントラスト画像と定量的位相コントラスト画像を同時に形成する方法 | |
US8687253B2 (en) | Speckle noise reduction based on longitudinal shift of sample | |
US20130070251A1 (en) | Systems and Methods of Dual-Plane Digital Holographic Microscopy | |
EP2915009B1 (fr) | Système d'imagerie holographique auto-référencé | |
EP2438413A2 (fr) | Capteur d'imagerie de front d'onde | |
US10006860B2 (en) | Digital holography recording device, digital holography playback device, digital holography recording method, and digital holography playback method | |
US20140307261A1 (en) | White Light Diffraction Tomography of Unlabeled Live Cells | |
Ayoub et al. | A method for assessing the fidelity of optical diffraction tomography reconstruction methods using structured illumination | |
US20080309944A1 (en) | Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus | |
Jayakumar et al. | Sampling Moiré method: a tool for sensing quadratic phase distortion and its correction for accurate quantitative phase microscopy | |
Malik et al. | A practical criterion for focusing of unstained cell samples using a digital holographic microscope | |
Lee et al. | Two step on-axis digital holography using dual-channel Mach-Zehnder interferometer and matched filter algorithm | |
Joglekar et al. | Compact, low cost, large field-of-view self-referencing digital holographic interference microscope | |
Romero et al. | Digital holographic microscopy for detection of Trypanosoma cruzi parasites in fresh blood mounts | |
FR3050038A1 (fr) | Procede et dispositif de microscopie interferentielle plein champ en lumiere incoherente | |
Almoro et al. | Object wave reconstruction by speckle illumination and phase retrieval | |
KR100871270B1 (ko) | 광학 간섭계를 이용한 샘플 시료의 이미지 정보 생성 방법및 장치 | |
KR20170037800A (ko) | 투명 매질의 미세 결함 검출방법 및 시스템 | |
Grebenyuk et al. | Theory of imaging and coherence effects in full-field optical coherence microscopy | |
WO2013083563A1 (fr) | Methode d'integration de phase holographique differentielle | |
Liebling et al. | Local amplitude and phase retrieval method for digital holography applied to microscopy | |
Wei et al. | Rapid quantitative interferometric microscopy using fast Fourier transform and differential–integral based phase retrieval algorithm (FFT-DI-PRA) | |
JP6196841B2 (ja) | 透過波面計測装置及び透過波面計測方法 |
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: 12809139 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 12809139 Country of ref document: EP Kind code of ref document: A1 |