EP1461660B1 - Ensemble imageur a instrument optique a synthese d'ouvertures - Google Patents

Ensemble imageur a instrument optique a synthese d'ouvertures Download PDF

Info

Publication number
EP1461660B1
EP1461660B1 EP02803826A EP02803826A EP1461660B1 EP 1461660 B1 EP1461660 B1 EP 1461660B1 EP 02803826 A EP02803826 A EP 02803826A EP 02803826 A EP02803826 A EP 02803826A EP 1461660 B1 EP1461660 B1 EP 1461660B1
Authority
EP
European Patent Office
Prior art keywords
zones
transfer function
central section
modulation transfer
smaller
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
EP02803826A
Other languages
German (de)
English (en)
Other versions
EP1461660A1 (fr
Inventor
Eric Thomas
Philippe Blanc
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Thales SA
Original Assignee
Thales SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Thales SA filed Critical Thales SA
Publication of EP1461660A1 publication Critical patent/EP1461660A1/fr
Application granted granted Critical
Publication of EP1461660B1 publication Critical patent/EP1461660B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/58Optics for apodization or superresolution; Optical synthetic aperture systems

Definitions

  • the present invention relates to an imaging assembly with an optical aperture synthesis instrument.
  • optical imaging systems with aperture synthesis. These systems consist of a combination of subsystems each having a reduced size so as to obtain practically the same result as a large optical system.
  • a system composed of small subsystems is easier to achieve than a large optical system and has fewer launching problems with a satellite.
  • an imaging system of this type has a telescope function.
  • the combination of the signals obtained by each of the instruments is generally done by interferometry Fizeau or Michelson.
  • the figure 1 represents, by way of illustration, a pupil configuration of an interferometer with three telescopes.
  • the main optical characteristic of an aperture synthesis optical apparatus is determined by the diameter of the pupils 10, 12, 14 and their respective positions.
  • the pupils all have the same diameter D 'and their centers are arranged according to an equilateral triangle of base B'.
  • Such aperture synthesis devices may have a modulation transfer function (MTF) whose support is of discontinuous type when the distance separating the pupils is sufficiently large relative to their diameter, that is to say say that this function has cancellation ranges.
  • MTF modulation transfer function
  • the figure 2 represents the modulation transfer function of the instrument represented on the figure 1 .
  • the modulation transfer function of an optical device is the response of this instrument to the various input spatial frequencies.
  • the columnar spatial frequencies are represented on the abscissa u axis and the line frequencies are represented on the ordinate v axis.
  • LH means that the zone concerns low frequencies (L) in line and high frequencies (H) in column.
  • the designations HH, H'H ', H'H, HH' mean high frequencies both online and in column.
  • the images are sampled with a spatial sampling frequency that can be different in line and in column.
  • the central rectangle 20 0.0 of the useful part of the spectrum is replicated to form a paving with identical blocks to the central pad 20 constituting 0.0 replicas 20 1.0, 20 0.1, 20 1.1, 20 0, -1 , 20 -1,1 , etc.
  • the sampling frequency can not fall below the values corresponding to those of the figure 3 because, for lower sampling frequencies, the FTM zones of the replicas overlap the MTF zones of the middle section 0,0 , which would result in a spectrum fold, i.e. a deterioration of the spectrum. 'information.
  • the number of pixels (picture elements) needed to sample the image is, of course, a direct function of the sampling frequency. It As a result, the higher the sampling frequency, the higher the number of pixels required.
  • the invention allows a substantially lower sampling frequency, without spectrum folding.
  • the invention makes it possible, at a constant field, to reduce the number of pixels necessary to restore images and to increase the quality of the image (signal-to-noise ratio).
  • the invention takes advantage of the fact that the central zone 0,0 has important parts in which the modulation transfer function is zero.
  • in-line and column sampling rates are chosen whose values are such that replicas overlap this central section, but avoid intersections between non-zero modulation transfer function (circular) zones (ie that is, non-zero FTM zones of replicas fill zones of zero value of the central section of the MTF), and in the central section thus filled, those which are closest to the origin of the frequencies and which make it possible to reconstruct the modulation transfer frequency are selected and reconfigured.
  • this new central section smaller in size than the original central section, so that it corresponds, for image reconstruction, to the original central section.
  • the reconstruction of the image is preferably done on the ground.
  • the in-line sampling step can be increased by a factor of about 1.62 and the column sampling step can be increased by a factor of about 1.67.
  • the total number of pixels necessary for the sampling of the image is, with a constant field, decreased by a factor of 2.7.
  • the invention generally relates to an imaging assembly comprising an aperture synthesis instrument having the features summarized in claim 1.
  • the values of the sampling frequencies are such that at least in one direction, in line or in column, the areas of the reduced central section of the modulation transfer function are tangent.
  • the assembly is for example of the interferometer type, including telescope.
  • a method for determining the optical sampling rate of an imaging assembly comprising a multi-pupil aperture synthesis type instrument, such as a telescope.
  • interferometer which is characterized in that this sampling frequency is chosen at a value lower than that imposed by the Shannon criterion and such that the replicas close to the central section of the modulation transfer function overlap this central section of in such a way that the non-zero transfer function zones of these replicas are intercalated in zones with zero transfer function of the central section, without intersection with a non-zero transfer function zone, and in such a way that a reduced section inside the central section contains all the useful zones to reconstruct the modulation transfer function.
  • the line and column sampling frequency is selected so that at least in one direction, in line or in column, the areas of the reduced section that are adjacent are tangent.
  • the invention also relates to an image reconstruction method comprising the features summarized in claim 6.
  • the section (replica) 20 '1.0 for which the center of the central region LL 1 to 0 abscissa and ordinate, overlaps the central section 20 0.0.
  • the areas H 'H', H and H'H of this replica 20 ' 1.0 are within the central section 0.0 .
  • the area H 1.0 is between the areas LL 0.0 and LH 0.0 .
  • the indices for the zones correspond to the coordinates of the center of the central zone LL of the corresponding replica.
  • the LH- 1.0 zone lies between the zones H 0.0 and LL 0.0 . It can also be seen that the zone H'H ' 0.1 is between the zones LH -1.0 and H'H 0.0 and is tangent to these two zones. Likewise, the zone HH ' 0.1 is between the zones HH 0.0 and L'H 1.0 and is tangent to these two zones. In a symmetrical manner, the zone H'H 0, -1 is located between the zones LH -1.0 and H'H ' 0.0 and is tangent to these two zones; finally, the zone HH 0.1 is between the zones H 1.0 and HH ' 0.0 and is tangent to these two zones.
  • a central section 22 0 of smaller dimension than the central section 20 0.0 include the seven areas of optical transfer function: LL, L'H, LH, HH ', H'H, HH and H'H ', the LL zone from the central section 0.0 while the zones H'H' and HH 'come from the replica 0.1, the zones H'H and HH come from the replica 0, - 1, the H area is from the 1.0 replica and the LH area is from the -1.0 replica.
  • D 'and B' constitutes the maximum value of D 'which makes it possible to use the method according to the invention. Indeed, for larger values of D ', it is no longer possible to insert replicas in the central zone.
  • the gain on the signal-to-noise ratio is 4.52
  • the gain on the signal-to-noise ratio is 4 , 52 that is 2.13.
  • Such an interferometer is represented on the figure 6 . It comprises four circular pupils 30, 32, 34 and 36 of diameter D 'whose centers are arranged along a diagonal square B'.
  • B ' 4D', which is the minimum value of D 'making it possible to use the method according to the invention.
  • this FTM has two zones whose centers are on the abscissa axis, namely the areas LH and L'H.
  • the center of area LH has abscissa B and that of area H has abscissa -B.
  • two areas HL and H'L present their centers on the y-axis.
  • the ordinate of the center of the zone HL is + B, while the ordinate of the zone H'L is -B.
  • LH means a low frequency in line and a high frequency in column
  • HL a high frequency in line and a low frequency in column
  • the same column and in-line sampling frequency can be chosen and this sampling frequency can be chosen with a value such that, in the reduced central section, all circular areas are tangent.
  • this sampling frequency in column and in line is 6D.
  • a reduced central section 40 ' 0,0 comprising a central zone LL 0,0 and eight other peripheral circular zones HH, H'H, H'L, LH, H, H'H', HL and HH 'from the central section and replicas.
  • the indices corresponding to the coordinates of the center of the area LL of the corresponding replica have been assigned.
  • the reduced central section 40 ' 0,0 comprises, in addition to the area LL 0,0 , the zones H'H 0,0 , HH 0,0 , H'H' 0,0 and HH ' 0,0 .
  • the area HL 0, -1 is inserted between the zones H'H '0,0 and HH' 0,0 and is tangential to these areas. This zone is also tangent to the central zone LL 0,0 .
  • the area H 1.0 is inserted between the zones HH 0,0 and HH ' 0,0 and their is tangent as well as the zone LL 0,0 .
  • the zone H'L 0.1 is interspersed between the zones H'H 0,0 and HH 0,0 and their is tangent as well as the central zone LL 0,0 .
  • the LH -1.0 zone is interspersed between the zones H'H 0,0 and H'H ' 0,0 and is tangent to them as well as to the central zone LL 0,0 .
  • the value (6D) of the sampling frequency is equal to five-thirds of the value of the sampling frequency corresponding to the Shannon criterion.
  • the gain on the signal-to-noise ratio is between 4.66 and 2.16.
  • the invention is, of course, not limited to a number of pupils equal to three or four. It applies regardless of the number of pupils of the interferometer or, more generally, of the instrument for synthesizing apertures. However, regardless of the embodiment, the ratio between the diameter of each pupil and the distance between these pupils must have a maximum value that can be easily determined so that the zones with zero transfer function are sufficiently wide to allow inserting non-zero replica transfer function media into the central section of this MTF.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
  • Microscoopes, Condenser (AREA)
  • Mounting And Adjusting Of Optical Elements (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Fire-Detection Mechanisms (AREA)
  • Optical Elements Other Than Lenses (AREA)
  • Non-Portable Lighting Devices Or Systems Thereof (AREA)

Description

  • La présente invention est relative à un ensemble imageur à instrument optique à synthèse d'ouvertures.
  • Il est envisagé, notamment dans l'industrie spatiale, d'utiliser des systèmes optiques imageurs à synthèse d'ouvertures. Ces systèmes consistent en une combinaison de sous-systèmes ayant chacun une taille réduite de façon à obtenir pratiquement le même résultat qu'un système optique de grandes dimensions. Un système composé de sous-systèmes de petites dimensions est plus aisé à réaliser qu'un système optique de grande dimension et pose moins de problèmes de lancement à l'aide d'un satellite. A titre d'exemple, un système imageur de ce type a une fonction de télescope.
  • La combinaison des signaux obtenus par chacun des instruments s'effectue en général par interférométrie de Fizeau ou de Michelson.
  • Le document Applied Optics, Vol. 12, No. 3, pages 487-492 décrit un système à synthèse d'ouvertures destiné à être utilisé dans le domaine infrarouge.
  • La figure 1 représente, à titre d'illustration, une configuration de pupilles d'un interféromètre à trois télescopes.
  • La caractéristique optique principale d'un appareil optique à synthèse d'ouvertures est déterminée par le diamètre des pupilles 10, 12, 14 et leurs positions respectives. Dans cet exemple, les pupilles ont toutes le même diamètre D' et leurs centres sont disposés selon un triangle équilatéral de base B'.
  • On sait que de tels appareils à synthèse d'ouvertures peuvent présenter une fonction de transfert de modulation (FTM) dont le support est de type discontinu lorsque la distance séparant les pupilles est suffisament grande par rapport à leur diamètre, c'est-à-dire que cette fonction présente des plages d'annulation. La figure 2 représente la fonction de transfert de modulation de l'instrument représenté sur la figure 1.
  • La fonction de transfert de modulation d'un appareil optique est la réponse de cet instrument pour les diverses fréquences spatiales d'entrée. Dans le diagramme de la figure 2, les fréquences spatiales en colonne sont représentées sur l'axe u des abscisses et les fréquences en ligne sont représentées sur l'axe v des ordonnées. Le support de la fonction de transfert de modulation de l'instrument représenté sur la figure 1 présente ainsi sept zones circulaires, chacune ayant un diamètre 2D (avec D = λʹ ,
    Figure imgb0001
    λ étant la longueur d'onde) : une zone LL "centrale" pour les basses fréquences (L) en colonne et en ligne, et six zones circulaires périphériques : LH, HH, H'H, L'H, H'H' et HH'. La désignation LH signifie que la zone concerne des basses fréquences (L) en ligne et des hautes fréquences (H) en colonne. De même, les désignations HH, H'H', H'H, HH' signifient des hautes fréquences tant en ligne qu'en colonne.
  • Les centres des zones circulaires périphériques se trouvent sur un cercle de diamètre 2B ( B = λ
    Figure imgb0002
    ) centré à l'origine.
  • La ligne joignant les centres des cercles HH et H'H' forme un angle de 60° avec l'axe des abscisses et, de même, la ligne joignant les centres des cercles H'H et HH' forme un angle de 120° avec cet axe des abscisses.
  • En pratique, les images sont échantillonnées avec une fréquence d'échantillonnage spatial qui peut être différente en ligne et en colonne. Pour éviter le repliement du spectre, c'est-à-dire une perte d'information, il faut respecter le théorème de Shannon, c'est-à-dire que la fréquence d'échantillonnage en colonne doit être supérieure ou égale à 2 fois la fréquence maximale du spectre à restituer soit 2 fois la distance 1 (figure 2) telle que 1 = B + D et, en ligne, la fréquence d'échantillonnage doit être égale à deux fois la distance µ (figure 2) qui vaut : μ = 3 B / 2 + D .
    Figure imgb0003
  • Ainsi, quand la fréquence d'échantillonnage en colonne est 2B + 2D et la fréquence d'échantillonnage en ligne est 3 B / 2 + D ,
    Figure imgb0004
    on obtient une image dont le support du spectre est du type de celui représenté sur la figure 3 avec une section centrale de spectre comprenant les sept zones circulaires LL, L'H, LH, HH, H'H', H'H et HH' et délimitée sur la figure 3, par un rectangle 200,0. Le spectre comporte aussi un ensemble de répliques identiques à la section centrale qui sont décalées sur l'axe des abscisses d'un nombre entier de fréquences d'échantillonnage en colonne, et sur l'axe des ordonnées, d'un nombre entier de fréquences d'échantillonnage en ligne.
  • Ainsi, comme on le voit sur la figure 3, le rectangle 200,0 central de la partie utile du spectre est répliqué pour former un pavage avec des pavés identiques au pavé central 200,0 constituant des répliques 201,0, 200,1, 201,1, 200,-1, 20-1,1, etc.
  • Jusqu'à présent, on a considéré que la fréquence d'échantillonnage ne pouvait pas descendre au-dessous des valeurs correspondant à celles de la figure 3 car, pour des fréquences d'échantillonnage inférieures, les zones de FTM des répliques chevauchent les zones FTM de la section centrale 200,0, ce qui aboutirait à un repliement du spectre, c'est-à-dire à une détérioration de l'information.
  • Le nombre de pixels (éléments d'image) nécessaire à l'échantillonnage de l'image est, bien entendu, une fonction directe de la fréquence d'échantillonnage. Il en résulte que plus la fréquence d'échantillonnage est élevée et plus le nombre de pixels nécessaires est élevé.
  • L'invention permet une fréquence d'échantillonnage sensiblement plus faible, sans repliement de spectre. Autrement dit, l'invention permet, à champ constant, de diminuer le nombre de pixels nécessaires pour restituer des images et d'augmenter la qualité de l'image (rapport signal à bruit).
  • L'invention tire parti du fait que la zone centrale 200,0 présente des parties importantes dans lesquelles la fonction de transfert de modulation est nulle.
  • Le procédé consiste :
    • à conférer aux fréquences d'échantillonnage en colonne et en ligne des valeurs inférieures aux fréquences d'échantillonnage qui sont déterminées par le critère de Shannon, ces fréquences d'échantillonnage de valeurs inférieures fournissant des répliques occupant des parties à fonction de transfert nulle de la section centrale du spectre, sans former d'intersection avec les zones à fonction de transfert non nulle de cette section centrale, l'ensemble des zones utiles différentes les plus rapprochés de l'origine formant alors une nouvelle section centrale de dimension plus faible que la section centrale d'origine, et
    • pour reconstruire l'image, à reconfigurer selon la section centrale d'origine la section centrale de faible dimension.
  • Ainsi, on choisit des fréquences d'échantillonnage en ligne et en colonne dont les valeurs sont telles que des répliques chevauchent cette section centrale, mais évitent des intersections entre les zones (circulaires) à fonction de transfert de modulation non nulle (c'est-à-dire que les zones à FTM non nulle des répliques remplissent des zones de valeur nulle de la section centrale de la FTM), et on sélectionne dans la section centrale ainsi remplie celles des zones qui sont les plus rapprochées de l'origine des fréquences et qui permettent de reconstruire la fréquence de transfert de modulation, et on reconfigure cette nouvelle section centrale, de dimension plus faible que la section centrale d'origine, de façon qu'elle corresponde, pour la reconstruction d'image, à la section centrale d'origine.
  • Dans le cas d'un instrument embarqué, la reconstruction de l'image s'effectue de préférence au sol.
  • On montre que dans le cas du système interférométrique à trois pupilles représenté sur la figure 1, le pas d'échantillonnage en ligne peut être augmenté d'un facteur d'environ 1,62 et le pas d'échantillonnage en colonne peut être augmenté d'un facteur d'environ 1,67. Ainsi, le nombre total de pixels nécessaires à l'échantillonnage de l'image est, à champ constant, diminué d'un facteur 2,7.
  • Ainsi, l'invention concerne de façon générale un ensemble imageur comprenant un instrument à synthèse d'ouvertures comportant les caractéristiques résumées dans la revendication 1.
  • Dans une réalisation, les valeurs des fréquences d'échantillonnage sont telles qu'au moins selon une direction, en ligne ou en colonne, les zones de la section centrale réduite de la fonction de transfert de modulation sont tangentes.
  • L'instrument peut comprendre trois pupilles de diamètre D' dont les centres sont disposés selon un triangle équilatéral de côté B', tels que : 2 = 3 / 4 ,
    Figure imgb0005
    la fréquence d'échantillonnage en colonne étant égale à 3B/2 et la fréquence d'échantillonnage en ligne étant égale à 6D, avec B = λ
    Figure imgb0006
    et D = λʹ ,
    Figure imgb0007
    λ étant une longueur d'onde.
  • L'instrument comprend, dans une variante, quatre pupilles dont chacune présente un diamètre D' et dont les centres sont disposés selon un carré de diagonale B' tels que : 2D' = B'/2, les fréquences d'échantillonnage en colonne et en ligne étant de 6D, avec D = λʹ ,
    Figure imgb0008
    λ étant une longueur d'onde.
  • L'ensemble est par exemple du type interféromètre, notamment télescope.
  • Il est aussi décrit un procédé pour déterminer la fréquence d'échantillonnage optique d'un ensemble imageur comprenant un instrument du type à synthèse d'ouvertures de plusieurs pupilles, tel qu'un télescope interférométrique, qui est caractérisé en ce qu'on choisit cette fréquence d'échantillonnage à une valeur inférieure à celle imposée par le critère de Shannon et telle que les répliques voisines de la section centrale de la fonction de transfert de modulation chevauchent cette section centrale de façon que les zones à fonction de transfert non nulle de ces répliques s'intercalent dans des zones à fonction de transfert nulle de la section centrale, sans intersection avec une zone à fonction de transfert non nulle, et de façon telle qu'une section réduite à l'intérieur de la section centrale contienne l'ensemble des zones utiles pour reconstruire la fonction de transfert de modulation.
  • Dans une réalisation, on choisit la fréquence d'échantillonnage en ligne et en colonne de façon telle qu'au moins selon une direction, en ligne ou en colonne, les zones de la section réduite qui sont voisines sont tangentes.
  • L'invention concerne également un procédé de reconstruction d'images comprenant les caractéristiques résumées dans la revendication 6.
  • D'autres caractéristiques et avantages de l'invention apparaîtront avec la description de certains de ses modes de réalisation, celle-ci étant effectuée en se référant aux dessins ci-annexés sur lesquels :
    • la figure 1, déjà décrite, représente un interféromètre à trois pupilles,
    • la figure 2, déjà décrite, représente le support de la fonction de transfert de modulation avant échantillonnage de l'interféromètre représenté sur la figure 1,
    • la figure 3, également déjà décrite, représente, pour l'interféromètre représenté sur la figure 1, le support du spectre après échantillonnage correspondant au critère de Shannon,
    • la figure 4 représente, pour l'interféromètre représenté sur la figure 1, le support du spectre résultant d'une fréquence d'échantillonnage déterminée avec le procédé conforme à l'invention,
    • les figures 5a, 5b et 5c sont des schémas illustrant un traitement pour retrouver le spectre d'origine en vue de la reconstruction de l'image après utilisation du procédé représenté sur la figure 4, et
    • les figures 6, 7, 8, 9, 10a, 10b et 10c sont des figures analogues aux figures, respectivement 1, 2, 3, 4, 5a, 5b et 5c, pour un interféromètre à quatre pupilles représenté sur la figure 6.
  • On se réfère tout d'abord à la figure 4.
  • Dans cet exemple correspondant à l'interféromètre à trois pupilles représenté sur la figure 1, on choisit une fréquence d'échantillonnage en colonne ayant pour valeur 3B/2 et une fréquence d'échantillonnage en ligne de valeur 6D.
  • Dans ce cas, comme on peut le voir sur la figure 4, la section (réplique) 20'1,0, pour laquelle le centre de la zone centrale LL a 1 pour abscisse et 0 pour ordonnée, chevauche la section centrale 200,0. Ainsi, les zones H'H', L'H et H'H de cette réplique 20'1,0 se trouvent à l'intérieur de la section centrale 200,0. En particulier, on observe que la zone L'H1,0 se trouve entre les zones LL0,0 et LH0,0. Les indices pour les zones correspondent aux coordonnées du centre de la zone centrale LL de la réplique correspondante.
  • De même, la zone LH-1,0 se trouve entre les zones L'H0,0 et LL0,0. On voit aussi que la zone H'H'0,1 se trouve entre les zones LH-1,0 et H'H0,0 et est tangente à ces deux zones. De même, la zone HH'0,1 se trouve entre les zones HH0,0 et L'H1,0 et est tangente à ces deux zones. De façon symétrique, la zone H'H0,-1 se trouve entre les zones LH-1,0 et H'H'0,0 et est tangente à ces deux zones ; enfin, la zone HH0,1 se trouve entre les zones L'H1,0 et HH'0,0 et est tangente à ces deux zones.
  • On voit ainsi que dans une section centrale 220 de dimension plus réduite que la section centrale 200,0 on retrouve les sept zones de fonction de transfert de modulation : LL, L'H, LH, HH', H'H, HH et H'H', la zone LL provenant de la section centrale 0,0 alors que les zones H'H' et HH' proviennent de la réplique 0,1, les zones H'H et HH proviennent de la réplique 0,-1, la zone L'H provient de la réplique 1,0 et la zone LH provient de la réplique -1,0.
  • Ainsi, bien qu'on ait choisi des fréquences d'échantillonnage en ligne et en colonne inférieures à celles correspondant au critère de Shannon, il ne se produit pas de repliement de spectre car on obtient dans une section centrale l'ensemble des sept zones de la fonction de transfert de modulation et ces zones ne se chevauchent pas.
  • Ensuite, pour reconstruire l'image, on procède comme représenté sur les figures 5a, 5b et 5c, c'est-à-dire qu'après avoir obtenu la section centrale réduite 220 (figures 4 et 5a), on procède à un réarrangement des zones LL, LH, L'H, HH, H'H, HH' et H'H' de façon que ces zones soient disposées comme représenté dans la section 200 de la figure 3, c'est-à-dire comme sur la figure 5c. A cet effet, on permute HH et HH', on permute H'H' et H'H, et on permute LH et L'H'. On obtient ainsi la position présentée sur la figure 5b. Ensuite, il suffit d'éloigner de la zone centrale LL les six zones périphériques H'H, HH, L'H, LH, HH' et H'H' pour obtenir la configuration représentée sur la figure 5c, c'est-à-dire la configuration de départ de la section centrale du spectre.
  • On a observé qu'avec ce procédé, bien que les parties à fonction de transfert nulle contiennent de l'énergie provenant du bruit présent à toutes les fréquences lors de l'échantillonnage, le signal obtenu est de même qualité, c'est-à-dire sans perte d'information, et sensiblement de même amplitude que dans le cas où l'on effectue l'échantillonnage correspondant au critère de Shannon.
  • Lorsque 2 = 3 / 4 ,
    Figure imgb0009
    le rapport entre la fréquence d'échantillonnage en colonne quand on respecte le critère de Shannon (figure 3) et la fréquence d'échantillonnage déterminée avec le procédé de l'invention (figure 4) a pour valeur 8 + 3 / 6 1 , 62.
    Figure imgb0010
    Ainsi, on peut augmenter le pas entre les pixels dans ce rapport.
  • La relation indiquée ci-dessus entre D' et B' constitue la valeur maximale de D' qui permet d'utiliser le procédé conforme à l'invention. En effet, pour des valeurs de D' plus importantes, il n'est alors plus possible d'intercaler des répliques dans la zone centrale.
  • Ce même rapport est, en ligne, de 5/3 (toujours, dans l'hypothèse où : 2 = 3 / 4 ,
    Figure imgb0011
    ). Ainsi le pas entre les pixels en ligne peut être augmenté d'un facteur d'environ 1,67.
  • Le nombre total de pixels peut donc être réduit d'un facteur 1,62 x 1,67 = 2,7 pour un champ constant.
  • Étant donné qu'avec cet échantillonnage, pour un champ donné, on diminue le nombre de pixels, on peut donc augmenter la taille des pixels. Dans le cas d'un système d'observation de type "push-broom", pour lequel on utilise une seule ligne de détecteurs, le temps d'intégration peut être augmenté. Dans ces conditions dans le cas d'un système de type push-broom, le flux collecté par pixel est multiplié par : 1,62 x 1,672 (1,672 provient du fait que l'intégration est effectuée le long des colonnes). Ainsi, le flux collecté par pixel est multiplié par un facteur d'environ 4,52 par rapport à un échantillonnage respectant le critère de Shannon. Dans le cas le plus favorable (qui intervient quand les bruits sont indépendants de la surface du détecteur), le gain sur le rapport signal à bruit est de 4,52, et dans le cas défavorable (qui intervient quand le bruit dépend de la surface du détecteur), le gain sur le rapport signal à bruit est de 4 , 52
    Figure imgb0012
    soit 2,13.
  • On va maintenant décrire en relation avec les figures 6, 7, 8, 9 et 10a, 10b, 10c un exemple d'application de l'invention dans le cas où l'interféromètre comporte quatre pupilles.
  • Un tel interféromètre est représenté sur la figure 6. Il comprend quatre pupilles 30, 32, 34 et 36 circulaires de diamètre D' dont les centres sont disposés selon un carré de diagonale B'. Dans cet exemple, B'=4D', qui est la valeur minimale de D' permettant de faire appel au procédé conforme à l'invention.
  • Le support de la fonction de transfert de modulation d'un tel interféromètre est représentée sur la figure 7. Cette fonction de transfert de modulation comporte neuf zones circulaires ayant toute un diamètre 2D ( D = λ
    Figure imgb0013
    ), soit B/2 ( B = λ
    Figure imgb0014
    ). En plus de la zone centrale LL, dont le centre est à l'origine des fréquences en colonne et en ligne, cette FTM comporte deux zones dont les centres sont sur l'axe des abscisses, à savoir les zones LH et L'H. Le centre de la zone LH a l'abscisse B et celui de la zone L'H a l'abscisse -B. En outre, deux zones HL et H'L présentent leurs centres sur l'axe des ordonnées. L'ordonnée du centre de la zone HL est +B, tandis que l'ordonnée de la zone H'L est -B.
  • Les centres des zones HL, LH, H'L et L'H forment ainsi un carré et au centre de chaque côté de ce carré se trouvent les centres des quatre autres zones circulaires, à savoir : HH, HH', H'H' et H'H. Comme dans le diagramme de la figure 2, LH signifie une basse fréquence en ligne et une haute fréquence en colonne, et HL une haute fréquence en ligne et une basse fréquence en colonne.
  • Lorsque la fréquence d'échantillonnage correspond au critère de Shannon avec une valeur minimale 2B + 2D tant en colonne qu'en ligne, on obtient pour la FTM (figure 8) une section centrale carrée 400,0 et des répliques de cette section carrée, à savoir 401,0, 401,-1, 401,1, 400,1, etc., les indices ayant la même signification que dans le cas du diagramme de la figure 3.
  • Dans ce cas, quand on applique le procédé conforme à l'invention, on peut choisir la même fréquence d'échantillonnage en colonne et en ligne et on peut choisir cette fréquence d'échantillonnage avec une valeur telle que, dans la section centrale réduite, toutes les zones circulaires sont tangentes.
  • Dans l'exemple représenté sur la figure 9, cette fréquence d'échantillonnage en colonne et en ligne est de 6D. Ainsi, on obtient une section centrale réduite 40'0,0 comportant une zone centrale LL0,0 et huit autres zones circulaires périphériques HH, H'H, H'L, LH, L'H, H'H', HL et HH' provenant de la section centrale et de répliques. Comme dans le diagramme de la figure 4, sur le diagramme de la figure 9, on a affecté à chaque zone circulaire les indices correspondant aux coordonnées du centre de la zone LL de la réplique correspondante.
  • La section centrale réduite 40'0,0 comporte, outre la zone LL0,0, les zones H'H0,0, HH0,0, H'H'0,0 et HH'0,0. En outre, la zone HL0,-1 s'intercale entre les zones H'H'0,0 et HH'0,0 et est tangente à ces zones. Cette zone est également tangente à la zone centrale LL0,0. La zone L'H1,0 s'intercale entre les zones HH0,0 et HH'0,0 et leur est tangente ainsi qu'à la zone LL0,0. De même, la zone H'L0,1 s'intercale entre les zones H'H0,0 et HH0,0 et leur est tangente ainsi qu'à la zone centrale LL0,0. Enfin, la zone LH-1,0 est intercalée entre les zones H'H0,0 et H'H'0,0 et leur est tangente ainsi qu'à la zone centrale LL0,0.
  • Pour la reconstruction de l'image, il est nécessaire de reconfigurer la section centrale réduite 40'0,0 de façon telle que les neuf zones soient réparties comme les zones de la section centrale 400,0 représentée sur la figure 8. A cet effet, partant de la section 40'0,0, on effectue des permutations pour que la zone HL se trouve du côté des ordonnées positives et la zone H'L du côté des ordonnées négatives, et de même on permute les zones LH et L'H. On obtient ainsi la configuration de la figure 10b. Il suffit ensuite de positionner les zones HL, LH, L'H et H'L pour obtenir la configuration représentée sur la figure 10c correspondant à la section centrale 400,0 de la figure 8.
  • La valeur (6D) de la fréquence d'échantillonnage est égale aux cinq tiers de la valeur de la fréquence d'échantillonnage correspondant au critère de Shannon. Ainsi, le pas entre pixels peut être augmenté d'un facteur d'environ 1,67 et le nombre total de pixels nécessaires à l'échantillonnage peut, à champ constant, être réduit de 1,672 = 2,79. Etant donné que la taille des pixels peut être augmentée et que, de même, le temps d'intégration peut aussi être augmenté dans le cas d'un système push-broom, le flux collecté par pixel est multiplié par un facteur : 1,672 x 1,67 = 4,66. Le gain sur le rapport signal sur bruit est compris entre 4,66 et 2,16.
  • L'invention n'est, bien entendu, pas limitée à un nombre de pupilles égal à trois ou quatre. Elle s'applique quel que soit le nombre de pupilles de l'interféromètre ou de façon plus générale, de l'instrument à synthèse d'ouvertures. Cependant, quel que soit le mode de réalisation, le rapport entre le diamètre de chaque pupille et la distance entre ces pupilles doit présenter une valeur maximale qui peut être déterminée aisément afin que les zones à fonction de transfert nulle soient suffisamment étendues pour permettre d'intercaler des supports de fonction de transfert de réplique non nulle dans la section centrale de cette FTM.

Claims (7)

  1. Ensemble imageur comprenant un instrument à synthèse d'ouvertures comportant une pluralité de pupilles séparées (10, 12, 14 ; 30, 32, 34, 36) avec des moyens de détection par échantillonnage, dans lequel les valeurs des fréquences d'échantillonnage en colonne et en ligne sont inférieures aux valeurs correspondant au critère de Shannon, et sont telles que des zones de fonction de transfert de modulation de répliques résultant de l'échantillonnage s'insèrent dans des zones à fonction de transfert de modulation nulle de la section centrale de la fonction de transfert de modulation, de façon telle que ces zones ne forment pas d'intersection et que, dans une section centrale de dimension plus réduite que la section centrale de la fonction de transfert de modulation, soit présent l'ensemble des zones représentant la section centrale de la fonction de transfert de modulation, l'ensemble imageur comportant en outre des moyens pour reconfigurer le spectre de l'image acquise afin de la reconstruire, caractérisé en ce que :
    lesdits moyens permettent de positionner les zones de la section centrale de dimension réduite (40'0,0 ; 220) de façon qu'elles correspondent aux positions des zones de la section centrale de la fonction de transfert de modulation d'origine afin de rétablir la fonction de transfert de modulation à partir des zones se trouvant dans la section réduite.
  2. Ensemble selon la revendication 1 caractérisé en ce que les valeurs des fréquences d'échantillonnage sont telles qu'au moins selon une direction, en ligne ou en colonne, les zones de la section centrale réduite de la fonction de transfert de modulation sont tangentes.
  3. Ensemble selon la revendication 1 ou 2 comprenant trois pupilles (10, 12, 14) de diamètre D' dont les centres sont disposés selon un triangle équilatéral de côté B', tels que : 2 = 3 / 4 ,
    Figure imgb0015
    la fréquence d'échantillonnage en colonne étant égale à 3B/2 et la fréquence d'échantillonnage en ligne étant égale à 6D, avec B = λ
    Figure imgb0016
    et D = λʹ ,
    Figure imgb0017
    λ étant une longueur d'onde.
  4. Ensemble selon la revendication 2 comprenant quatre pupilles (30, 32, 34, 36) dont chacune présente un diamètre D' et dont les centres sont disposés selon un carré de diagonale B' tels que : 2D' = B'/2, les fréquences d'échantillonnage en colonne et en ligne étant de 6D, avec D = λʹ ,
    Figure imgb0018
    λ étant une longueur d'onde
  5. Ensemble selon l'une quelconque des revendications précédentes caractérisé en ce qu'il est du type interféromètre, notamment télescope.
  6. Procédé de reconstruction d'images dans un ensemble imageur comprenant un instrument du type à synthèse d'ouvertures de plusieurs pupilles séparées (10,12,14 ;30,32,34,36) avec des moyens de détection par échantillonnage, dans lequel on choisit la fréquence d'échantillonnage à une valeur inférieure à celle imposée par le critère de Shannon et telle que les répliques voisines de la section centrale de la fonction de transfert de modulation chevauchent cette section centrale de.façon que les zones à fonction de transfert non nulle de ces répliques s'intercalent dans des zones à fonction de transfert nulle de la section centrale de la fonction de transfert de modulation, sans intersection avec une zone à fonction de transfert non nulle; et de façon telle que dans une section centrale de dimension plus réduite que la section centrale de la fonction de transfert de modulation, et située à l'intérieur de la section centrale de la fonction de transfert de modulation, soit présent l'ensemble des zones utiles pour reconstruire la fonction de transfert de modulation, caractérisé en ce que l'on positionne les zones de la section centrale de dimension réduite (40'0,0 ; 220) de façon qu'elles correspondent aux positions des zones de la section centrale de la fonction de transfert de modulation d'origine afin de rétablir la fonction de transfert de modulation à partir des zones se trouvant dans la section centrale de dimension réduite.
  7. Procédé selon la revendication 6, caractérisé en ce qu'on choisit la fréquence d'échantillonnage en ligne et en colonne de façon telle qu' au moins selon une direction, en ligne ou en colonne, les zones de la section réduite qui sont voisines sont tangentes.
EP02803826A 2001-11-27 2002-11-21 Ensemble imageur a instrument optique a synthese d'ouvertures Expired - Lifetime EP1461660B1 (fr)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR0115303A FR2832814B1 (fr) 2001-11-27 2001-11-27 Ensemble imageur a instrument optique a synthese d'ouvertures
FR0115303 2001-11-27
PCT/FR2002/003989 WO2003046642A1 (fr) 2001-11-27 2002-11-21 Ensemble imageur a instrument optique a synthese d'ouvertures

Publications (2)

Publication Number Publication Date
EP1461660A1 EP1461660A1 (fr) 2004-09-29
EP1461660B1 true EP1461660B1 (fr) 2010-05-12

Family

ID=8869820

Family Applications (1)

Application Number Title Priority Date Filing Date
EP02803826A Expired - Lifetime EP1461660B1 (fr) 2001-11-27 2002-11-21 Ensemble imageur a instrument optique a synthese d'ouvertures

Country Status (7)

Country Link
US (1) US7180597B2 (fr)
EP (1) EP1461660B1 (fr)
AT (1) ATE467857T1 (fr)
DE (1) DE60236371D1 (fr)
ES (1) ES2343241T3 (fr)
FR (1) FR2832814B1 (fr)
WO (1) WO2003046642A1 (fr)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2858425B1 (fr) * 2003-07-31 2005-09-30 Cit Alcatel Instrument d'observation a synthese d'ouverture optique a pupilles en mouvement et procede associe
FR2879761B1 (fr) * 2004-12-20 2007-07-20 Cit Alcatel Procede de correction de la fonction de transfert d'un interferometre imageur a pupilles multiples
FR2919933B1 (fr) * 2007-08-06 2010-06-25 Onera (Off Nat Aerospatiale) Dispositif optique multi-pupilles aminci

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3556630A (en) * 1968-08-14 1971-01-19 Ibm Method and apparatus for obtaining,by a series of samples,the intensity distribution across sources of incoherent electromagnetic waves to produce a single composite picture
US4136954A (en) * 1976-12-29 1979-01-30 Jamieson John A Imaging apparatus including spatial-spectral interferometer

Also Published As

Publication number Publication date
US20040257263A1 (en) 2004-12-23
FR2832814A1 (fr) 2003-05-30
ES2343241T3 (es) 2010-07-27
EP1461660A1 (fr) 2004-09-29
WO2003046642A1 (fr) 2003-06-05
US7180597B2 (en) 2007-02-20
DE60236371D1 (de) 2010-06-24
ATE467857T1 (de) 2010-05-15
FR2832814B1 (fr) 2004-02-27

Similar Documents

Publication Publication Date Title
EP2567537B1 (fr) Procede d'imagerie polychrome
CA3059260A1 (fr) Dispositif et procede d'imagerie multispectrale dans l'infrarouge
EP1472563B1 (fr) Instrument d'optique d'observation depuis un satellite, a deux telescopes
EP3824269A2 (fr) Procédé et systèmes de caracterisation optique non invasive d'un milieu hétérogène
FR3007534A1 (fr) Realisation d'un filtre spectral nanostructure
EP1461660B1 (fr) Ensemble imageur a instrument optique a synthese d'ouvertures
WO2017042494A1 (fr) Caméra plénoptique
EP4154046B1 (fr) Instrument d'imagerie
King et al. The current ability of HST to reveal morphological structure in medium-redshift galaxies
EP2767860A1 (fr) Dispositif de dépolarisation amelioré
EP2176633B1 (fr) Procédé d'estimation d'au moins une déformation du front d'onde d'un système optique ou d'un objet observé par le système optique et dispositif associé
EP2567538A1 (fr) Procede d'imagerie
WO2021156581A1 (fr) Appareil pour l'observation d'un astre
WO2016139284A1 (fr) Procédé et dispositif d'étalonnage au sol d'un instrument optique
EP3839814B1 (fr) Instrument d'observation optique embarque a resolution spatiale et spectrale variables
EP3798710B1 (fr) Telescope de type cassegrain a plan de focal segmente
Breckinridge et al. Innovative pupil topographies for sparse aperture telescopes and SNR
EP1544667B1 (fr) Complément fréquentiel d'images acquises par un dispositif d'acquisition à synthèse d'ouverture optique dilué, par combinaison d'images fréquentiellement complémentaires
WO2020169905A1 (fr) Spectrometre a imagerie pupillaire
FR3110003A1 (fr) Ecran pour l’affichage d’une image projetée comportant un revêtement structuré.
FR3133086A1 (fr) Telescope a cinq miroirs
FR2879761A1 (fr) Procede de correction de la fonction de transfert d'un interferometre imageur a pupilles multiples
WO2014086539A1 (fr) Systeme optique comportant un filtre optique a masque de phase et/ou d'amplitude a motif periodique
FR2683928A1 (fr) Procede et dispositif de formation d'images par traitement des donnees fournies par des capteurs optroniques.
EP2012526A1 (fr) Procédé et dispositif d'observation

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20040628

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SK TR

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: ALCATEL LUCENT

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: THALES

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SK TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REF Corresponds to:

Ref document number: 60236371

Country of ref document: DE

Date of ref document: 20100624

Kind code of ref document: P

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2343241

Country of ref document: ES

Kind code of ref document: T3

REG Reference to a national code

Ref country code: NL

Ref legal event code: T3

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

REG Reference to a national code

Ref country code: IE

Ref legal event code: FD4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100813

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100913

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

Ref country code: IE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

26N No opposition filed

Effective date: 20110215

BERE Be: lapsed

Owner name: THALES

Effective date: 20101130

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 60236371

Country of ref document: DE

Effective date: 20110214

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20101130

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20101130

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20101130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20101130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20101121

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100512

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100812

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 14

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 15

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 16

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 20171115

Year of fee payment: 16

Ref country code: DE

Payment date: 20171114

Year of fee payment: 16

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 17

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: ES

Payment date: 20181203

Year of fee payment: 17

Ref country code: FR

Payment date: 20181026

Year of fee payment: 17

Ref country code: GB

Payment date: 20181119

Year of fee payment: 17

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 60236371

Country of ref document: DE

REG Reference to a national code

Ref country code: NL

Ref legal event code: MM

Effective date: 20181201

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20181201

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20190601

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20191121

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191130

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191121

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20210528

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191122