METHOD FOR PHOTOCOPYING A SEQUENCE OF CUT SURFACES INSIDE A LIGHT-SCATTERING OBJECT WITH IMPROVED SCANNING

Abstract

The invention relates to a free-beam interferometric method for illuminating a sequence of sectional areas in the interior of the light-scattering object. The method makes it possible for the user to select a larger image field and/or a higher image resolution than previously possible with the occurrence of self-interference of the specimen light from a scattering specimen.

Claims

1. A method for illuminating a sequence of sectional areas in the interior of a light-scattering object, the method comprising: provision of a light source, which emits light with a predetermined central wavelength .sub.0 and a coherence length less than 25 microns; splitting the light of the light source into specimen light and reference light; illumination of the object over the surface with the specimen light; imaging of the specimen light scattered by the object with numerical aperture NA and magnification factor M onto an electronic camera with pixels with a pitch P along at least one axis in the camera plane; causing interference of reference light and specimen light on the camera by establishing a path length profile and a phase gradient of the reference light along the predetermined axis in the camera plane; displacement of the path length profile of the reference beam at a time-dependent rate; and detection of further camera images in each case at least indexed with a dimension for the time-dependent displacement of the path length profile, wherein the condition $\frac{\sqrt{2}}{3+\sqrt{2}}.Math.{}_0<\frac{2.Math.\mathrm{NA}*P}{M}\frac{\sqrt{2}}{1+\sqrt{2}}.Math.{}_0$ is met.

2. The method according to claim 1, wherein the phase gradient along the predetermined axis in the camera plane is established on half the Nyquist frequency.

3. The method according to any one of the preceding claims, wherein a surrounding interval of the indexing is predetermined for each indexed camera image of at least a partial number of the detected camera images and a weighted stack mean value of all the camera images of the surrounding interval is calculated and the calculated stack mean value of the indexed camera image is deducted to calculate the structurally improved camera image.

4. The method according to claim 3, further comprising calculation of the images of the sectional areas in the interior of the object by two-dimensional Fourier filtering of the structurally improved camera images.

5. The method according to claim 1, wherein the light-scattering object is the retina of a living eye.

6. The method according to claim 5, further comprising biometric verification of the identity of a user.

Description

[0029] The single FIG. 1 also serves for the following explanation of the invention. In the FIGURE:

[0030] FIG. 1 shows a sketch of the Fourier space related to a quadratic camera with the definition ranges of the Fourier components a) of the cross-correlated and autocorrelated interference signals according to the prior art and b) solely of the cross-correlated useful signal according to the present invention.

[0031] WO 2017/029160 A1 states that the imaging of the object is to be set up such that average speckle diameter D is greater than two pixel pitches P, and that furthermore the phase gradient in the camera is to be set from the interval between 2/D and /P. This is equivalent to the requirement on angle of incidence a of the reference beam onto the camera plane:

[00003]$\frac{\mathrm{NA}}{1.22*M}<\sin .Math.<\frac{{}_0}{2*P}$

[0032] Here, M denotes the magnification factor of the image, i.e. the ratio of image size to object size, so that D=M*D.sub.0 applies with

[00004]$D_0=\frac{1.2.Math.2*{}_0}{\mathrm{NA}}$

[0033] as the speckle size related to the object plane. On account of the random distribution of the light intensity in a speckle field, the diameter of the speckle can only be given as a statistical magnitude, which can be calculated from the autocorrelation function. For a circular aperture, an Airy function ensues (see Dainty J., The statistics of speckle patterns, Progress in Optics, E. Wolf, ed. (Elsevier), pp. 1-46, 1977). The first zero point is usually defined as a measure for the average speckle size, which then also corresponds to the diameter of the Airy disc of the point spread function with a diffraction-limited resolution.

[0034] The right-hand side of the above inequality requires the scanability of the stripe pattern with the given camera (Nyquist condition), and the left-hand side of the inequality requires that the stripe pattern in the individual speckles must be able to be detected on the camera.

[0035] With regard to the Fourier coefficients of the cross-correlated, a stripe pattern with an amount of the wave number vector {right arrow over (k)}.sub.F in the vicinity of the Nyquist frequency is however not expedient, because then fine structures in the sought sectional areas lead to frequencies which would be located beyond the Nyquist frequency, i.e. beyond the limit of the scanability.

[0036] According to the invention, it is on the contrary advantageous rather to approach half the Nyquist frequency with |{right arrow over (k)}.sub.F| and to favourably establish the amount for the half bandwidth k of the useful signal by selecting average speckle diameter D.

[0037] Speckle diameter D is limited downwards by two conditions: [0038] (i) D must be at least as large as the fringe spacing of the stripe pattern generated by the phase gradient on the camera. [0039] (ii) D must be so large that the bandwidth of the cross-correlated useful signal modulated with the phase gradient does not contain any frequency components beyond the Nyquist frequency of the camera.

[0040] If condition (i) is not met, it is not then possible to determine the phase positions of the speckles, and if condition (ii) is not met, useful signal information is lost.

[0041] For half bandwidth k of the useful signal, condition (ii) signifies an upper limit. This can be seen from FIG. 1 a), in which the Fourier space of a quadratic camera is sketched. The represented circles describe the definition ranges of the non-zero Fourier coefficients of the cross-correlated signal modulated with k.sub.F=|{right arrow over (k)}.sub.F| top left and bottom right) and of the autocorrelated signal (by k=0), wherein the circle radii are the half bandwidths of respectively k and 2 k. The circle centre points are displaced in relation to one another precisely by |{right arrow over (k)}.sub.F|=3 k. In order that the definition ranges on the one hand do not overlap and on the other hand lie completely in the Fourier space of the camera,

[00005]$.Math.k\frac{\sqrt{2}}{3+\sqrt{2}}.Math.\frac{}{P}$

must be selected.

[0042] If the autocorrelated interference signal is not however regarded as an interference signal at all or is simply ignored, use can be made, as sketched in FIG. 1 b), of a broadband up to

[00006]$.Math..Math.k=\frac{\sqrt{2}}{1+\sqrt{2}}.Math.\frac{}{P},$

at which the Fourier components do not overlap with the complex conjugates in the definition range (except at k=0).

[0043] According to the invention, the aim is to make accessible the broadband interval

[00007]$\frac{\sqrt{2}}{3+\sqrt{2}}.Math.\frac{}{P}<.Math..Math.k=2.Math..Math.\frac{\mathrm{NA}}{{}_0}\frac{\sqrt{2}}{1+\sqrt{2}}.Math.\frac{}{P}$

hitherto problematic for the image evaluation. This can be directly reformulated into a requirement on the measurement setup, namely into

[00008]$\frac{\sqrt{2}}{3+\sqrt{2}}.Math.{}_0<2.Math.\mathrm{NA}*P\frac{\sqrt{2}}{1+\sqrt{2}}.Math.{}_0$

as a setup requirement for the camera parameters in relation to the wavelength used. The phase gradient can then be also selected at

[00009]$.Math.{\stackrel{.fwdarw.}{k}}_F.Math.<3.Math..Math..Math.k=6.Math..Math.\frac{\mathrm{NA}}{{}_0}.$

In order to actually maximise the half bandwidth of the useful signal, the phase gradient should preferably be put precisely at the half Nyquist frequency. In many practical cases, however, a maximisation is not sought, but only an improvement. Other phase gradients, which are still smaller than |{right arrow over (k)}.sub.F|=3 k, can then also be set by another selection of the angle of incidence of the reference light onto the camera.

[0044] It should be noted that in FIG. 1 the Nyquist frequency amounts to

[00010]$\frac{\sqrt{2}.Math.}{P}$

along the diagonal, because a quadratic camera with quadratic pixels arranged in a checker-board manner and pixel pitch P is adopted there. Along the diagonal, the effective pixel pitch is then shortened by the factor {square root over (2)}, i.e. there is a higher scanning density, which is known per se.

[0045] The speckle diameter on the camera must always simultaneously meet the two conditions (i) and (ii). Nothing is changed to condition (i) by the present invention. By the use of smaller phase gradients than in the prior art, condition (ii) permits greater bandwidths according to the invention and therefore smaller speckle diameters than previously. The speckle size can be reduced by a factor up to:

[00011]$\frac{\sqrt{2}+3}{\sqrt{2}+1}=2.Math.\sqrt{2}-11.8.Math.2$

[0046] Even if speckle size D is necessarily selected intentionally somewhat larger than according to conditions (i) and (ii), in order for example to have a safety distance from the edge of the Fourier space of the camera, a markedly improved scanning of the cross-correlated is already achieved with the given camera.

[0047] Approximately three times the Fourier space can be utilised with the invention.

[0048] Three facts now have to be stated and emphasised: [0049] I) The method of WO 2017/029160 A1 is precisely specified in a development by the present invention. In particular, greater freedom can thus be taken in the selection of the angle of incidence of the reference light onto the camera. [0050] II) The method according to the invention gains more structural information of internal sectional areas with higher image resolution and/or enlarged image field, because it is possible to transfer to smaller speckle diameters and larger bandwidths of the useful signal than previously possible in the prior art [0051] III) The actual presence of the autocorrelated in the detected individual images is simply ignored during the actual image detection. No isolated background images (without reference) are detected, and the autocorrelated also do not have to appear the same over all the specimen depths.

[0052] It could at first be surprising to the person skilled in the art that ignoring the autocorrelated would be possible without harm.

[0053] The invention relates from the outset to a sequence of camera images which is obtained by changing the path length of the reference light in the manner of a Time-Domain-OCT. Not a single one of these images is accessible to direct filtering to separate background and useful signal. But in the combined view of the images indexed with the path length, it appears that a numerical averaging procedure according to the invention virtually erases the useful signal in a mean value image. A still strongly structured image of the mean values then remains, which is independent of the reference light and as such can be perceived as a background image. This background image essentially comprises the autocorrelated for the specimen depth, from the surrounding area of which the path length interval is selected for averaging the indexed individual images, and furthermore also all the static structures from the instrumental set-up, e.g. dust particles on lenses and suchlike.

[0054] It is clear that a path length interval can be defined for each individual image detected and indexed with path length l.sub.0 of the reference light. A distinction can be made between two cases:

[0055] If the path lengths of the interval differ at most by coherence length l.sub.c of the light, e.g.

[00012]$\left[l_0-\frac{l_c}{2},l_0+\frac{l_c}{2},\right]$

the individual images comprise speckles in respect of indexes from this interval, which change from constructive into destructive interference with a variation of the path length of the reference light by .sub.0/2. If the reference mirror is displaced with a technical drive by several wavelengths between two individual images, the brightness of these speckles is virtually a random magnitude distributed around 0.5. If a plurality of such images are added, all the speckles assume on average the same value of 0.5, and the structure of the useful signal becomes invisible. The Fourier coefficients of the cross-correlated becomewith the exception of the integral, i.e. of the coefficient k=0zero.

[0056] For the measurement of larger specimen volumes, however, it may be more favourable to carry out only two or fewer measurements per coherence length in order to reduce the measurement time. Since the measurements in this case no longer originate from the same coherence volume, the speckle patterns in the cross-correlated are usually no longer correlated in the case of measurements in scattering fabric. Here, therefore, the phases and also the spatial positions of the speckles are randomly distributed. In this case, too, the cross-correlated can be erased by averaging over a none too large number of imagestypically up to 20 images.

[0057] A path length interval of the form [l.sub.0l, l.sub.0+l,] should be designated as a surrounding interval of path length of l.sub.0. A partial number of detected camera images is assigned to the surrounding interval. The partial number can also be denoted as a stack of individual images. If these are for example images at the path lengths l.sub.j[l.sub.0l, l.sub.0+l], j=1, . . . , N, a weighted stack mean value reads

[00013]$B_{a.Math.v}\left(x,y,l_0\right)=\frac{1}{N}.Math.\underset{j=1}{\overset{N}{.Math.}}.Math.\left[B\left(x,y,l_j\right)-B_0\left(l_j\right)\right]G\left(l_j\right)$

wherein x, y denote pixel coordinates and B.sub.av, B, B.sub.0[0,1] denote brightness values (or grey scale values) on the pixels. Parameter B.sub.0 can be used to standardise the average brightness of all the individual images if this appears necessary. Parameter G Can be used as a weighting factor once again if required to standardise the contrast breadth of the images of the stack. They are in addition suitable for leaving out of account individuale.g. obviously erroneousimages and/or for limiting the effective range of the summation, e.g. by diminishing weights with the distance to l.sub.0.

[0058] Ultimately, the user will, in view of a detected sequence of camera images, determine in post-processing the precise selection of the parameters of the stack averaging. The determination of the parameters can also be automated, in particular can be carried out with software on a conventional computer, if the properties of the image detection device are known to the user and he specifies specific requirements on the background images. An illustrative criterion to the effect that the calculated stack mean value represents a good approximation for a background image independently of the reference arm length is its slow variability of all the image components in the transition from l.sub.0 to a nearby depth plane, i.e. the following should hold B.sub.av(x, y, l.sub.0)B.sub.av(x,y,l.sub.0+dl).

[0059] After the determination of the background images for all the detected camera images, the differences

B(x,y,l.sub.0)=B(x,y,l.sub.0)B.sub.av(x,y,l.sub.0)

can be calculated in the simplest and preferred case and the structurally improved images B can be further processed. It is expedient to then subject the structurally improved images to a Fourier transformation and to eliminate any remaining coefficients close to k=0 and the complex conjugates (Fourier filtering). The Fourier reverse transformation then leads to an image of an inner sectional area of the illuminated object indexed with the path length displacement l.sub.0 of the reference light, i.e. the sectional area can be assigned to a depth position in the object.

[0060] It is thus possible, with the method according to the invention for the illumination of inner sectional areas of scattering objects, which relates specifically to a sequence of camera images at different path lengths of the reference light and observes the condition according to the invention for the illumination parameters, to arrive at a dataset, which in post-processing generates sectional images with a higher resolution and larger image field than in the prior art. This post-processing is optionally chronologically downstream of the data acquisition within the scope of the disclosure of WO 2017/029160 A1, it can even take place not until very much later by a third party and to this extent does not have to be an integral component of the illumination method according to the invention or of a measurement device. The generation of the raw dataset is itself also a valuable achievement in its own right, the usefulness of which is not in doubt according to the above explanations.

[0061] The previously described invention is preferably used for the illumination of biological scattering specimens, particularly preferably living tissue, very particularly preferably the illumination of the retina of a living eye. Apart from the medical applications in ophthalmology, the invention can also preferably be used by a user for the biometric verification of the identity of the user.