Method and device for wide-field astronomical imaging

Abstract

According to one aspect, the present description relates to a wide-field astronomical imaging device of a scene comprising: an optical system (101) including an optical axis (270) movable in rotation about at least two axes of rotation; a two-dimensional optical detector (102) configured for the acquisition of raw images of portions of said scene with said optical system (101); a control unit (104) configured for the control of said optical axis (270) and the acquisition, over time, of a plurality of raw images defining a plurality of different portions of said scene, two successive portions of said plurality of portions having at least one common part, said plurality of raw images corresponding respectively to a plurality of positions of said optical axis;and a processing unit (105) configured to calculate, for each acquisition of a raw image, an extended image (240) larger than an extended image (210) calculated previously. FIG. 2;

Description

Title of the invention: Method and device for wide-field astronomical imaging Technical field

[0001] The present description relates to a method and a device for wide-field astronomical imaging. State of the art

[0002] An astronomical imaging instrument generally comprises a two-dimensional detector and an optical system configured to generate an image of a scene, for example a constellation in the sky.

[0003] The object field of such an astronomical imaging instrument, called the “intrinsic” field of the instrument, corresponds to the set of points in the object space of this instrument of which an image can be formed, and depends on several parameters, including the size of the detection surface of the detector and the characteristics of the optical system (e.g. focal length, aperture).

[0004] In the field of astronomy, the objects studied (the stars) often cover a very large portion of the object space (the sky) which exceeds the dimensions of the intrinsic field of the instrument, which prevents the observation of the entirety of these objects.

[0005] In this context, it is useful to increase the object field of the instrument beyond its intrinsic field to obtain images corresponding to a larger object field than what the instrument could natively capture. In the present description, these images are called "expanded field" images.

[0006] Rather than increasing the intrinsic field of the instrument by changing the characteristics of the optical system, it is possible to obtain an expanded field by using the so-called "mosaicing" technique, described for example in [REF. 1].

[0007] In the mosaicking technique, illustrated for example in [Fig.l], the object studied (in this example a constellation 10 in the sky 15) is imaged in several stages. The instrument acquires images 11, 12, 13 of portions of the object (“raw images”) whose dimensions are limited by the intrinsic field of the instrument. The raw images 11, 12, 13 are then assembled by an image registration technique to reconstruct an enlarged field image 18 which contains the object studied 10 and which corresponds to a field larger than the intrinsic field of the instrument with which the raw images 11, 12, 13 are obtained.

[0008] To carry out the registration of the images, it is known to acquire successive raw images 11, 12, 13 having common parts 16, 17 (overlap zone) which allow these images to be compared.

[0009] Typically, in the mosaicking technique, the overlap rate, defined as the size of the common part 16 between two successive raw images 11, 12 relative to the size of one of the raw images 11, 12, is between approximately 5% and approximately 40%.

[0010] Furthermore, it is known to increase the signal-to-noise ratio (SNR) of an astronomical image by using the “stacking” technique described for example in [REF 2].

[0011] For this, a plurality of images of the same scene are acquired and summed together, the resulting image (called stacked image in the present description) having an SNR greater than the SNR of each initial image; this method can be seen as an averaging of the noise of the initial images to reduce the noise level.

[0012] It would be possible to combine mosaicking and stacking by acquiring a plurality of raw images for each portion of the scene and then summing them. The stacked images of each portion thus obtained then have an improved SNR compared to the SNR of the raw images of each portion. The stacked images are then realigned to reconstruct the wide-field image, according to the known mosaicking method. A mosaic of stacks is thus produced. In such a method, the SNR of the wide-field image is in practice similar to the SNR of each stacked image.

[0013] The present description provides a wide-field imaging method that allows an increase in SNR compared to the SNR obtained with a stack mosaicking method as described above. Summary of the invention

[0014] In this description, the term "comprise" means the same as "include" or "contain", and is inclusive or open and does not exclude other elements not described or shown.

[0015] Furthermore, in the present description, the term "approximately" or "substantially" is synonymous with (means the same as) a lower and / or upper margin of 10%, for example 5%, of the respective value.

[0016] According to a first aspect, the present description relates to a device for wide-field astronomical imaging of a scene comprising: - an optical system comprising an optical axis movable in rotation along at least two axes of rotation; - a two-dimensional optical detector comprising a detection surface substantially perpendicular to said optical axis, said detection surface defining, with the optical system, a given intrinsic field, said two-dimensional detector being configured for the acquisition, by means of said optical system, of raw images of portions of said scene of dimensions equal to the dimensions of the intrinsic field; - a control unit configured to control the rotation of said optical axis, so as to allow the acquisition, over time, by the optical detector, of a plurality of raw images of a plurality of different portions of said scene, the superposition of said portions forming said widened field, two successive portions of said plurality of portions having at least one common zone, said plurality of raw images corresponding respectively to a plurality of positions of said optical axis; and a processing unit configured to calculate, for each acquisition of a raw image at a given time, an extended image of a size larger than a previously calculated extended image, said extended image being calculated from said previously calculated extended image and said raw image.

[0017] In the present description, the imaging device acquires raw images of "portions" of the scene, the geometric size of the portions being limited by the intrinsic field of the instrument.

[0018] All the portions of the scene acquired over time form a “field path”.

[0019] Two successive portions are called two portions acquired successively by the optical detector, the acquisitions of the images of two successive portions being separated by a time interval, for example inversely proportional to an acquisition rate of the detector.

[0020] In the present description, the two-dimensional detector is generally formed from a two-dimensional matrix of elementary detectors and a matrix of digital values ​​is called a "raw image", each digital value depending on a number of photons received by one or more elementary detectors of the two-dimensional detector when an image of a portion of the scene of dimensions equal to the dimensions of the intrinsic field is formed by the imaging system on the detection surface. The elements of the matrix are called "image points", and the digital value associated with a point of the image is called the "value of the image point".

[0021] In an astronomical imaging device according to the first aspect, by virtue of the fact that for each acquisition of a raw image, a new extended image is calculated from said raw image and a previously calculated extended image, a “sliding” stack is created which is likely to considerably increase the SNR due to the fact that a point of the extended image obtained at the end of the acquisition corresponds to a point of the object field acquired a very large number of times.

[0022] According to one or more exemplary embodiments, said common part between two successive portions of said plurality of portions has a surface area greater than or equal to approximately 85% of the intrinsic field. A coverage rate is thus defined, equal to the ratio between the common part between two successive portions of said plurality of portions and the intrinsic field. An astronomical imaging device according to the first aspect makes it possible to work, due to the sliding stack, with a high overlap rate, and to obtain an increase in the SNR compared to known devices.

[0023] In practice, it will be sought that two successive portions of said plurality of portions have an overlap strictly less than 100% of said intrinsic field of the optical system, in order to have successive raw images that are always distinct and to increase the size of the extended image after each acquisition of a raw image. Of course, it may happen during acquisition and due to experimental hazards that two successive portions are substantially confused, although this is not desired.

[0024] According to one or more examples, said extended image is calculated from the previously calculated extended image: - by summing the point values ​​of the previously calculated extended image to the point values ​​of the raw image which correspond to identical points of the object field; and - by completing the previously calculated extended image with points from the raw image which correspond to points in the object field not included in the previously calculated extended image.

[0025] By "summation" of point values ​​of several images together, it is understood within the meaning of the present description, any linear or non-linear combination operation between these point values.

[0026] In the device thus described, and unlike a device which would implement a simple stack mosaicking technique, point values ​​of the previously calculated extended image and of the new raw image are summed, when these points correspond to the same point of the object field.

[0027] This results in a greater increase in the SNR, in particular for a given observation time and / or a given expanded field size.

[0028] According to one or more exemplary embodiments, the summation comprises a simple sum, the point values ​​of the previously calculated extended image being added to the point values ​​of the raw image which correspond to identical points of the object field.

[0029] According to one or more exemplary embodiments, the summation comprises a weighted sum, the values ​​of the points of each image being able to be multiplied by one or more factors before being added. Said factors can for example be used to correct an exposure defect, or the effect of an optical aberration for only a few points of the image.

[0030] According to one or more exemplary embodiments, the summation of image point values ​​comprises a conditional sum, criteria being able to be defined to determine which points of the images must be summed.

[0031] According to one or more examples, the criteria comprise a maximum deviation of the value of a point of the image with respect to a median or an average calculated on values ​​of points of other previously acquired images, said points corresponding to the same point of the object field.

[0032] According to one or more examples, the control of said optical axis comprises compensation, at least partial, of the effects of the rotation of the Earth on the position of said portions in the scene.

[0033] This makes it possible to partially or totally compensate for the effect of the movement of the celestial vault linked to the rotation of the Earth and to determine a field path by controlling the position of each portion of the scene by rotation of the optical axis. Alternatively, the positions of the portions of the scene result from the movement of the celestial vault linked to the rotation of the Earth.

[0034] According to one or more examples, said plurality of positions of said optical axis is calculated from a predetermined observation time.

[0035] The observation time within the meaning of the present description (or “measurement time”) is the acquisition time of said plurality of raw images of the plurality of portions of the scene by the optical detector to obtain an enlarged field image of said scene.

[0036] This advantageously allows a user to choose the time it will take to acquire an image in wide field.

[0037] In this configuration, the user specifies the expanded field of which he wishes to obtain an image and the observation time. The processing unit then determines a field path from these two parameters and the successive positions of the optical axis to acquire the images of the portions of the field path.

[0038] In other exemplary embodiments, the observation time is not predetermined. For example, the user can decide to interrupt the measurement at the desired time.

[0039] According to one or more examples, the control unit is configured to control the rotation of said optical axis, so as to allow successive acquisitions of pluralities of raw images of portions of said scene, in order to scan the enlarged field several times. Thus, the portions of the scene acquired over time and forming the field path can be acquired cyclically, the portions of the scene for each cycle being able to be merged or offset.

[0040] This allows the user to obtain a wide-field image without time constraints. The device scans the wide-field several times to obtain a wide-field image. wide field whose SNR increases as the scans progress. The acquisition of the wide field image stops, for example, when the user interrupts the measurement.

[0041] According to a second aspect, the present description relates to a method of wide-field astronomical imaging of a scene comprising: - the acquisition, over time, by means of an astronomical imaging device, of a plurality of raw images of a plurality of different portions of said scene, each portion having dimensions equal to the dimensions of an intrinsic field of the imaging device, the superposition of said portions forming said enlarged field, two successive portions of said plurality of portions having at least one common part; and - the calculation, for each acquisition of a raw image at a given time, of an extended image of a size larger than a previously calculated extended image, said extended image being calculated from said previously calculated extended image and said raw image.

[0042] According to one or more examples, said extended image is calculated from the previously calculated extended image: - by summing the values ​​of the points of the previously calculated extended image to the values ​​of the points of the raw image which correspond to identical points of the object field; and - by completing the previously calculated extended image with points from the raw image which correspond to points in the object field not included in the previously calculated extended image.

[0043] According to one or more examples, the method further comprises, after the acquisition of each raw image, a step of filtering said raw image before the calculation of the extended image, to improve the quality of the extended image.

[0044] The filtering of said raw image comprises for example the suppression of possible artifacts from the raw images, eg the suppression of image points with abnormal values, the suppression of artifacts due to the detector and the optical system, color correction, suppression of the sky background and the suppression of certain stars.

[0045] According to one or more examples, the method comprises successive acquisitions of pluralities of raw images of portions of said scene, so as to scan the enlarged field several times.

[0046] In other words, the portions of the scene acquired over time and forming the field path can be acquired cyclically, the portions of the scene for each cycle being able to be merged or offset.

[0047] This allows the user to obtain a wide-field image without time constraints. The device scans the wide-field several times to obtain a wide-field image. wide field whose SNR increases as the scans progress. The acquisition of the wide field image stops, for example, when the user interrupts the measurement.

[0048] According to one or more examples, said common portion between two successive portions of said plurality of portions has a surface area greater than or equal to approximately 85% of the intrinsic field. In other words, this means that the overlap rate between two successive portions of said plurality of portions is greater than or equal to approximately 85%.

[0049] This allows for an increase in the SNR over a large part of each extended image, by maximizing the extent of the common parts between each new raw image acquired at a given time and each extended image calculated at a previous time.

[0050] This also makes it possible to improve the calculation of a new extended image at a given time by facilitating the comparison between a raw image acquired at said given time and an extended image calculated at a previous time because a greater number of image points of the two images correspond to the same points of the object field.

[0051] In practice, it will be sought that two successive portions of said plurality of portions have an overlap strictly less than 100% of said intrinsic field of the optical system, in order to have successive raw images that are always distinct and to increase the size of the extended image after each raw image acquisition.

[0052] The above-mentioned characteristics and advantages, as well as others, will appear on reading the detailed description which follows, of exemplary embodiments of the device and of the proposed method. This detailed description refers to the attached drawings. Brief description of the figures

[0053] Other advantages and characteristics of the invention will appear on reading the description, illustrated by the following figures:

[0054] [Fig.l], already described, represents a simplified diagram illustrating a mo-saiking method according to the prior art.

[0055] [Fig.2], represents a simplified diagram illustrating a wide-field astronomical imaging device according to the present description;

[0056] [Fig.3A], represents a simplified diagram of an object to be imaged according to the present description;

[0057] [Fig.3B], represents a simplified diagram illustrating an example of path portions field defined by keypoints according to this description;

[0058] [Fig.3C], represents a simplified diagram illustrating an example of a field path according to the present description;

[0059] [Fig.4], represents a simplified diagram illustrating a method of wide-field astronomical imaging according to the present description. Detailed description of the invention

[0060] In the figures, the elements are not shown to scale for better visibility.

[0061] [Fig. 2] illustrates a wide-field astronomical imaging device according to the present description. The device 100 comprises an optical system 101 with a rotating optical axis 270, a two-dimensional detector 102, for example arranged in a focal plane of the optical system 101, a control unit 104 for the rotation of the optical axis, and a processing unit 105. The device 100 may comprise a base 103 for holding said optical system 102 and supporting its movement.

[0062] The device 100 makes it possible to acquire several raw images of a scene, generally the sky, and to combine these images by means of the processing unit 105 configured to implement steps of the method according to the present description, in order to form an enlarged field image larger than one of the raw images.

[0063] For a given position of the optical axis, the optical system 101 and the detector 102 allow the acquisition of an image in said scene according to a given field, also called “intrinsic field” of the device in the present description. The intrinsic field depends mainly on the characteristics of the components used (for example the types of lenses, the vergences of the lenses, the number of detection elements of the detector, and the size of the detection surface) and the geometric arrangements between elements (for example the distances between lenses).

[0064] The detector 102 is formed for example from a two-dimensional set of detection elements and the image provided by the detector 102 is a matrix of digital values, each element of the matrix (or “point”) having a digital value which depends on the number of photons received by one or more detection elements.

[0065] The optical axis 270 of the optical system 101 is rotatable along at least two axes. It is possible to pivot the optical axis 270 along said axes so that the optical system 101 successively points at several portions of said scene and allows the successive acquisition of images of the successive portions on the detector 102.

[0066] Although in the diagram of [Fig. 2], it is shown that the optical axis 270 can be pivoted to a position 280 to successively acquire the images of the portions 205, 215, it will be understood that the optical axis is generally pivoted to allow the acquisition of a large number of raw images. Although this is not visible in [Fig. 2], in practice the overlap rate between the portions 205 and 215 is advantageously greater than 85%.

[0067] The successive pointing of several portions of the scene to acquire several images is called scanning the scene, or scanning the wide field, in the present description.

[0068] The succession of successive portions traversed during scanning is called a field path in the present description. Each portion has a size in the object field of the device which corresponds to the intrinsic field of the device.

[0069] The control unit 104 is configured for controlling the optical axis 270, that is to say in particular the rotation of the optical axis 270 and the parts integral with the optical axis (the optical system 101 and the detector 102) to perform a scan of a scene along a field path.

[0070] The control unit 104 is further configured to control the acquisition of the raw images by means of the detector 102 during the scanning of the scene by rotation of the optical axis 270. The control unit 104 makes it possible to define, for example, an exposure time for each acquisition of a raw image, an overlap rate between two successive portions of the field path or an observation time, that is to say the time taken to measure the image in an expanded field, in order to determine a particular field path.

[0071] According to one or more exemplary embodiments, the recommended observation time depends on the object observed and is on average greater than approximately 30 minutes. There is generally no maximum recommended time because the SNR of the image increases with the observation time.

[0072] The processing unit 105 is configured to process the raw images acquired by the detector 102.

[0073] In particular, the processing unit 105 is configured to calculate, after each acquisition of a raw image, a new extended image from an extended image calculated during the previous acquisition and from said raw image.

[0074] Generally, the processing unit 105 and the control unit 104 referred to in the present description may be part of the same system comprising

[0075] one or more physical entities, for example one or more computers. When in the present description, reference is made to calculation or processing steps for the implementation in particular of method steps, it is understood that each calculation or processing step can be implemented by software, hardware, firmware, microcode or any appropriate combination of these technologies. When software is used, each calculation or processing step can be implemented by computer program instructions or software code. These instructions can be stored or transmitted to a storage medium readable by the processing unit and the control unit and / or be executed by the control unit and the processing unit in order to implement these calculation or processing steps.

[0076] To scan a scene with the device according to the present description, several field paths are possible. Choosing a particular field path means choosing a way to scan the entire expanded field by acquiring images of portions of that field.

[0077] Generally, the field path can be defined by a list of coordinates defining each portion of the field path, for example a list of the coordinates of the successive centers of the portions of the field path.

[0078] According to one or more examples, to determine a field path from an expanded field, the control unit 104 may use a portion defined as the initial portion, a portion defined as the final portion, and then interpolate intermediate positions between the positions of the initial and final portions to define intermediate portions. The device then acquires a raw image of each intermediate portion.

[0079] According to one or more examples, it is possible to further define particular points (or "key points" in the present description) in the expanded field. The field path is then calculated to pass in an orderly manner through these key points, a succession of portions of the field path between successive key points being determined by interpolation according to the procedure described previously.

[0080] [Fig.3A], [Fig.3B] and [Fig.3C] schematically illustrate an example of constructing a field path in which the scan of an expanded field includes an object to be imaged 301, for example a galaxy.

[0081] [Fig.3A] illustrates the object to be imaged 301 of which an image is to be made with the device according to the present description. For this purpose, the device will scan the object to be imaged 301 along a field path.

[0082] In the example of [Fig.3B], 6 key points (k0, k1, k2, k3, k4 and k5) are defined. The field path starts from the reference key point k0, then passes successively through the key points k1, k2, k3, k4 and ends at the key point k5. The portions of the field path are here identified by the coordinates of their center. In the non-limiting example of [Fig.3B], the movement of the field path describes a curved shape, for example a spiral shape and the points k0 and k4 are merged.

[0083] [Fig.3C] illustrates, by way of example, how a straight line connecting the points k0 and kl is interpolated to define a plurality of intermediate positions i0 to i5, which correspond to intermediate field portions of which images will be acquired by the device according to the present description, between the image of the portion centered at k0 and the image of the portion centered at kl.

[0084] In the same way, it is possible to define intermediate portions between the other key points (from kl to k5) in order to fully define the field path making it possible to obtain the wide-field image comprising the object 301.

[0085] Although a particular field path has been described in detail using the [Fig.3A], [Fig.3B] and [Fig.3C], there are other embodiments of the device according to the present description in which the field path scans the expanded field in another manner. For example, it is possible to choose a field path such that the passage from one field portion to a next portion comprises either a movement in a first direction or a movement in a second direction orthogonal to the first direction.

[0086] The device according to the present description can operate in several modes of use.

[0087] In a first embodiment, called “time mode” in the present description, the acquisition of the wide-field image is carried out in a limited time, called observation time r0, which is chosen by the user.

[0088] The observation time is a time elapsing during the acquisition of the wide-field image, that is to say between the instant at which the user starts the acquisition and the instant at which the acquisition ends and where the user can observe the wide-field image that he wishes to obtain with the device.

[0089] The observation time can for example be selected by the user after estimating the size of the object studied (i.e. the enlarged field to be imaged) relative to the intrinsic field of the device.

[0090] For example, if the object studied covers an area of ​​the sky whose dimensions are four times greater than the dimensions of an area of ​​the sky that it is possible to image with the intrinsic field of the device, then the user can choose an observation time greater than four times the optimal time for acquiring a raw image, advantageously an observation time greater than or equal to 10 times this value.

[0091] In the time mode, the user of the device selects the observation time Lo and the expanded field to be imaged, then the field path is calculated by the control unit 104 as a function of the observation time and the expanded field.

[0092] More precisely, the coordinates of the portions of the field path can be calculated before the acquisition of the raw images, from the enlarged field, the observation time and auxiliary parameters preselected during the construction of the device.

[0093] In time mode, auxiliary parameters include for example: - an exposure time of a raw image Tp, that is to say the time interval during which the detector 102 will accumulate the photons coming from a portion of the scene in order to acquire an image of this portion; and - a size of a raw image which is equal to the intrinsic field of the instrument.

[0094] The user can, if desired, modify the auxiliary parameters to precisely adjust the characteristics of the acquisition of the raw images. The path of field is then determined based on these parameters and any key points as described previously.

[0095] According to a second embodiment, called “continuous mode” in the present description, the device performs a continuous scanning of the widened field to obtain a widened field image whose SNR increases with time.

[0096] In the continuous mode, the observation time is not limited, the user only chooses the wide field. The scanning of the wide field stops, for example, when the user chooses to interrupt the acquisition of the wide field image.

[0097] In the continuous mode, the field path can be calculated by the control unit 104 to scan the expanded field several times in different ways. This means that at each scan, the device acquires images of portions different from the portions already scanned at a previous scan. This makes it possible to obtain an expanded image whose SNR increases over time because, at each new scan, a greater number of points of the expanded image come from an overlapping area and can therefore be summed.

[0098] In the continuous mode, the coordinates of the portions of the field path can be calculated, before the acquisition of the raw images, from the expanded field and auxiliary parameters.

[0099] In continuous mode, the auxiliary parameters include for example: - the exposure time of a raw TP image; - the size of a raw image ( b^ by); - a recovery rate Rm; and - a shift between scans §y).

[0100] The overlap rate Rm is defined by the ratio between the common part between two successive portions during a sweep of the expanded field and the intrinsic field. Each sweep is generally carried out with the same overlap rate Rm.

[0101] The offset between scans is defined, for example, by a pair of coordinates equal to the differences between the coordinates of the position of the first portion of a scan and the coordinates of the position of the first portion of a successive scan. Choosing an offset between two successive scans for which the coordinates of the pair of coordinates are both non-zero makes it possible to obtain images of different portions at each scan.

[0102] The user can, if desired, modify the auxiliary parameters. In particular, the user can modify the overlap rate Rm if he wishes, for example, for the scanning of the expanded field to be more or less rapid.

[0103] [Fig.4] illustrates steps of a method of wide-field astronomical imaging according to the present description, which can for example be implemented with the device of [Fig.2].

[0104] In particular, [Fig.4] illustrates a step carried out iteratively during the acquisition of a plurality of raw images corresponding to portions of a field path in order to obtain an expanded field image.

[0105] In [Fig.4], at a given time T+1, a raw image 215 is acquired by the device 100 to increase the field of an extended image 210 obtained at an earlier time T.

[0106] The method according to the present description being iterative, the extended image 210 obtained at an earlier time T is itself calculated from a raw image acquired at a time T and an extended image obtained at a time T1. The first step (T = 0, not illustrated in [Fig.4]) corresponds to the acquisition of a first raw image, the extended image then being merged with this first raw image.

[0107] It will generally be sought that the two raw images acquired respectively at T and T+1 are two images acquired successively, that is to say separated by a minimum time interval, equal for example to the inverse of the acquisition rate of the imaging device. Nevertheless, it may happen that the time interval is sometimes greater under certain conditions, for example in the presence of an anomaly during the acquisition of the images or an error making one of the raw images unusable.

[0108] In a registration step 220, the raw image 215 acquired at time T+1 is compared to the extended image 210 to identify a common part. The registration of two images is a matching of the two images by searching for a geometric transformation of one of the images making it possible to superimpose, in the two images, the points corresponding to identical points of the scene.

[0109] During the comparison, it is a matter of establishing a correspondence between certain points of the raw image 215 and certain points of the extended image 210 to determine the points of the raw image 215 which correspond to identical points of the object field in the extended image, called common points 226 in the present description.

[0110] The points of the raw image 215 which are not common points 226 (and which are called new points 227 in the present description) are therefore points which correspond to points of the object field which have not yet been imaged.

[0111] Different types of geometric transformations of the raw image 215 may be used to match points of the raw image 215 and points of the extended image 210, for example rotations, translations, homotheties and homographies. These transformations are known to those skilled in the art.

[0112] In an exemplary embodiment for a simple case illustrated in [Fig.4], the imaged scene is considered to be planar and the movement of the scene relative to the focal plane of the optical system is the composition of a translation and a rotation, a geometric transformation model using a translation and a rotation can therefore be used to match the raw image 215 and the extended image 210.

[0113] In more complex embodiments, where the imaged scene is not planar and / or where the movement of the scene relative to the focal plane of the optical system comprises several rotations and translations, a geometric transformation model using, for example, combinations of translations, rotations and homotheties can be used to match raw images and extended images.

[0114] Once the common points 226 and the new points 227 of the raw image 215 have been identified, a completion step 230 is carried out to obtain at time T+1, a new extended image 240 from the raw image 215 and the extended image 210 at time T.

[0115] On the one hand, the numerical values ​​of the common points 226 of the raw image 215 are summed with the values ​​of the points of the extended image 210 to which they correspond, that is to say with the points of the extended image 210 which correspond to the same points of the object field. The summation may comprise for example a simple summation, a weighted summation, a conditional summation or a combination of these operations.

[0116] On the other hand, the new points 227 of the raw image 215 are added to those of the extended image 210 to form a new extended image 240 of a size greater than the extended image 210. It is said that the extended image 210 is “completed” with the new points 227 of the raw image 215, which are points which correspond to points of the object field not included in the extended image 210.

[0117] According to certain embodiments, after calculation of the extended image, the extended image can further be completed by a plurality of points to which a given value (for example zero) is assigned in order to obtain a rectangular extended image 241. This makes it possible to facilitate the numerical calculation which is more suited to rectangular images.

[0118] By virtue of the method according to the present description, after having acquired several raw images it is possible to obtain an extended image of size greater than each of the raw images and having an SNR greater than the SNR of each of the raw images on the overlapping zones.

[0119] For example, in an extended image calculated after acquiring two raw images of given SNR SNRx, the SNR of the points corresponding to an overlapping area between the two raw images will have a theoretical SNR equal to SNRb = SNRx-

[0120] Thus, advantageously when a high overlap rate is used (for example greater than 85%), the method according to the present description makes it possible to obtain an image in wide field whose SNR is increased compared to the SNR of the images raw SNRx for a large number of points. The SNR increase is close to the theoretical limit achievable, i.e. the limit obtained for a 100% coverage rate.

[0121] In contrast, in prior art methods, overlapping areas are generally not used to increase the SNR of the wide-field image. Thus, once two stacked images of different portions of the scene are registered, information from only one of the two images is used to obtain the points of the wide-field image corresponding to the overlapping area between the two portions.

[0122] Although described through a certain number of exemplary embodiments, the methods and devices for wide-field astronomical imaging according to the present description include various variations, modifications and improvements which will be obvious to those skilled in the art, it being understood that these various variations, modifications and improvements are part of the scope of the invention as defined by the following claims. References

[0123] [REF 1] SZELISKI, Richard. Image alignment and stitching. In : Handbook of ma-thematical models in computer vision. Springer, Boston, MA, p. 273-292, 2006

[0124] [REF 2] POPOWICZ, Adam, KUREK, Aleksander, POLLO, Agnieszka, et al. “Beyond the current noise limit in imaging through turbulent medium”. Optics letters, 2015, vol. 40, no 10, p. 2181-2184.

Claims

Claims

1. Device (100) for wide-field astronomical imaging of a scene comprising: an optical system (101) comprising an optical axis (270) movable in rotation along at least two axes of rotation; a two-dimensional optical detector (102) comprising a detection surface substantially perpendicular to said optical axis, said detection surface defining, with the optical system, a given intrinsic field, said two-dimensional detector (102) being configured for the acquisition, by means of said optical system (101), of raw images of portions of said scene of dimensions equal to the dimensions of the intrinsic field;a control unit (104) configured to control the rotation of said optical axis, so as to allow the acquisition, over time, by the optical detector (102), of a plurality of raw images of a plurality of different portions of said scene, the superposition of said portions forming said enlarged field, two successive portions of said plurality of portions having at least one common part, said plurality of raw images corresponding respectively to a plurality of positions of said optical axis; and a processing unit (105) configured to calculate, for each acquisition of a raw image (215) at a given time, an extended image (240) of a size larger than a previously calculated extended image (210), said extended image (240) being calculated from said previously calculated extended image (210) and said raw image (215).;

2. Device according to claim 1, wherein said extended image is calculated from the previously calculated extended image: by summing the point values ​​of the previously calculated extended image with the point values ​​of the raw image which correspond to identical points of the object field; and by supplementing the previously calculated extended image with points of the raw image which correspond to points of the object field not included in the previously calculated extended image.

3. Device according to any one of the preceding claims, in which said common part between two successive portions of said plurality of portions has a surface area greater than or equal to approximately 85% of the intrinsic field.

4. Device according to any one of the preceding claims, in which the control of said optical axis comprises compensation, at least partially, for the effects of the rotation of the Earth on the position of the portions in the scene.

5. A method for wide-field astronomical imaging of a scene comprising: acquiring, over time, by means of an astronomical imaging device, a plurality of raw images of a plurality of different portions of said scene, each portion having dimensions equal to the dimensions of an intrinsic field of the imaging device, superimposing said portions forming said wide field, two successive portions of said plurality of portions having at least one common part; and calculating, for each acquisition of a raw image (215) at a given time, an extended image (240) of a size larger than a previously calculated extended image (210), said extended image (240) being calculated from said previously calculated extended image (210) and said raw image (215).

6. A method according to claim 5, wherein said extended image is calculated from the previously calculated extended image: by summing the values ​​of the points of the previously calculated extended image with the values ​​of the points of the raw image which correspond to identical points of the object field; and by completing the previously calculated extended image with points of the raw image which correspond to points of the object field not included in the previously calculated extended image.

7. Method according to any one of claims 5 to 6, further comprising, after the acquisition of each raw image, a step of filtering said raw image before the calculation of the extended image.

8. A method according to any one of claims 5 to 7, wherein said plurality of portions of the scene are determined from a predetermined observation time.

9. A method according to any one of claims 5 to 8, wherein said common portion between two successive portions of said plurality of portions has a surface area greater than or equal to approximately 85% of the intrinsic field.

Images