WO1997005451A1 - Procede d'acquisition d'une image par balayage pousse-balai - Google Patents

Procede d'acquisition d'une image par balayage pousse-balai Download PDF

Info

Publication number
WO1997005451A1
WO1997005451A1 PCT/FR1996/001188 FR9601188W WO9705451A1 WO 1997005451 A1 WO1997005451 A1 WO 1997005451A1 FR 9601188 W FR9601188 W FR 9601188W WO 9705451 A1 WO9705451 A1 WO 9705451A1
Authority
WO
WIPO (PCT)
Prior art keywords
detectors
sampling
pixels
processing
acquiring
Prior art date
Application number
PCT/FR1996/001188
Other languages
English (en)
Inventor
Bernard Rouge
Christophe Latry
Jean-Claude Favard
Gilbert Pauc
Original Assignee
Centre National D'etudes Spatiales
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
Priority claimed from FR9509242A external-priority patent/FR2737375B1/fr
Priority claimed from FR9509241A external-priority patent/FR2737326B1/fr
Priority claimed from FR9509263A external-priority patent/FR2737376B1/fr
Application filed by Centre National D'etudes Spatiales filed Critical Centre National D'etudes Spatiales
Priority to DE69600785T priority Critical patent/DE69600785T2/de
Priority to EP96927095A priority patent/EP0782695B1/fr
Priority to JP9507289A priority patent/JP3012865B2/ja
Priority to US08/849,902 priority patent/US6215522B1/en
Publication of WO1997005451A1 publication Critical patent/WO1997005451A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/02—Picture taking arrangements specially adapted for photogrammetry or photographic surveying, e.g. controlling overlapping of pictures
    • G01C11/025—Picture taking arrangements specially adapted for photogrammetry or photographic surveying, e.g. controlling overlapping of pictures by scanning the object
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformation in the plane of the image
    • G06T3/40—Scaling the whole image or part thereof
    • G06T3/4084—Transform-based scaling, e.g. FFT domain scaling
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/48—Increasing resolution by shifting the sensor relative to the scene
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N3/00—Scanning details of television systems; Combination thereof with generation of supply voltages
    • H04N3/02—Scanning details of television systems; Combination thereof with generation of supply voltages by optical-mechanical means only
    • H04N3/08—Scanning details of television systems; Combination thereof with generation of supply voltages by optical-mechanical means only having a moving reflector

Definitions

  • the present invention relates to a method for acquiring a satellite observation image of the earth by "push-broom” scanning (Push-broom according to English terminology).
  • -saxonne from a bar or matrix of charge transfer type (CCD) detectors scrolling opposite the observed area.
  • CCD charge transfer type
  • FIG. 1 The principle of a sweep-type sweep has been illustrated in FIG. 1 in the case of a strip 1 of detectors.
  • the bar 1 performs the successive observation, during the scrolling of the satellite which carries it, of lines LI, L2, ..., Ln perpendicular to the movement (arrow D).
  • Instrumentation optics 2 at all times form the image of a landscape line on the line of detectors, the bar 1 being placed in the focal plane of optics 2, perpendicular to the speed vector of the satellite.
  • the landscape scrolls past each detector which integrates the light flux during the exposure time and transforms it into a proportional electric charge.
  • FIG. 2 illustrates a conventional processing chain for the images thus identified.
  • This processing chain schematically comprises a unit 3 for processing and amplification at the output of the detectors of the array 1, an analog-digital encoder 4 receiving the signal at the output of the unit 3, means 5 for transmission to the from the satellite of the digital images thus recorded, a unit 6 for the reconstruction on the ground of the images.
  • the unit 3 comprises in particular a shift register in which the information integrated and stored in each detector of the strip 1 in the form of a load is transferred at the end of an exposure time. This register then ensures an electronic transfer of the charges, which are converted into a succession of voltages proportional to the light fluxes received and integrated.
  • Unit 6 on the ground reconstructs the images, in particular by implementing deconvolution treatments intended to compensate for the instrumental defects, as well as, if necessary, interpolation treatments to reconstruct certain pixels of the image.
  • the acquisition time between two successive lines, or sampling time is such that the point on the ground located vertically above the satellite has moved a distance equal to the dimensions of an elementary detector projected on the ground .
  • the Modulation Transfer Function of an optical system is a function of frequency space which translates the ability of the system to transmit the various frequencies. It is characteristic of the restitution of the contrasts of the landscape by the system.
  • the Modulation Transfer Function depends mainly on the system optics, the yarn effect (sagging of contrast due to displacement), as well as detectors. Observation systems are characterized by the cutoff frequency fc, above which the FTM is negligible.
  • a cutoff frequency can be associated with each of the effects contributing to the Global Modulation Transfer Function (optical, spun and integration in the photosensitive zone), the overall cutoff frequency being the smallest of the previous three values.
  • the Modulation Transfer Function associated with integration on the photosensitive zone is canceled for a frequency equal to the inverse of the size of an elementary detector which constitutes as a first approximation the corresponding cutoff frequency.
  • An object of the invention is to propose a treatment method which allows:
  • oversampling is carried out by:
  • the superposition of the two bars can be obtained by means of an optical line divider or by using the scrolling of the satellite by shifting in the field of the second bar of an integer number of sampling steps in the direction of advancement and half a step in the direction of the bar.
  • the offset of the two sampling grids can be different from the half pixel, provided that it is taken into account in post-processing, the optimal value of the offset is however half the size of an elementary detector in row and in column.
  • FIG. 4 illustrates a matrix of pixels obtained with such processing. Pixels XB1 are those acquired with strip B1, pixels XB2 those acquired with strip B2.
  • Each of the two bars B1 and B2 therefore generates a rectangular sampling grid, the pitch of which along the lines is equal to the size of an elementary detector and the pitch along the columns of half this value, the two grids being offset by a distance equal to half of an elementary detector.
  • a square sampling grid is reconstituted, the line pitch of which coincides with the column pitch and is half the size of an elementary detector.
  • FIG. 5 Illustrated in FIG. 5 is the representation in the Fourier plane of a Transfer Function of Realistic modulation (FTM) for an observation instrument of the "PUSH-BROOM" type, used according to the acquisition method described in FR 2 678 460.
  • FTM Transfer Function of Realistic modulation
  • This FTM is expressed as the product of an ftml frequency transfer function normalized lines fx and a transfer function ftmc of the column frequencies fy, these two functions being given by the following formulas in the case of the SPOT 5 instrument:
  • the FTM is normalized at the origin and the frequencies fx and fy are normalized compared to twice the classical sampling frequency and the different curves represent level curves in steps of 0.01 up to 'at the value 0.1.
  • the Modulation Transfer Functions are weak over a large part of the frequency origin spectrum.
  • An object of the invention is therefore to propose an acquisition method which overcomes this drawback: it makes it possible to maintain a negligible spectral aliasing while optimizing the acquisition throughput.
  • the invention provides a method of acquiring and processing a satellite observation image of the earth by means of at least one strip or matrix of detectors of the charge transfer type scrolling opposite an observed area. , in which an acquisition of a plurality of pixels is carried out, which is sampled in staggered rows, characterized in that the row and column steps of this staggered sampling are such that the spectrum of the image is mainly located in an area of the Fourier plane where the Modulation Transfer Function is significant and where aliasing is limited.
  • a method for acquiring and processing a satellite observation image of the earth by means of at least two arrays or arrays of detectors of the same pitch, extending juxtaposed perpendicularly to their directions of movement and an integer plus a fraction of an elementary step parallel to their direction of movement, said bars or dies being offset one with respect to the other of an elementary step fraction, perpendicular to their direction of movement, characterized in that the sampling times of these two bars are equal to the inverse of the cutoff frequency of the elementary detectors.
  • a square sampling grid is then generated in the geometry defined by the speed of the sub-satellite point and its perpendicular, the sampling step in these two directions being equal to p / Vl + n2.
  • the sampling is then densified by a factor vl + n 2 at the cost of reducing the swath by the same factor.
  • a method for acquiring and processing a satellite observation image of the earth by means of a single array or matrix of detectors, with a time of sampling corresponding to the inverse of the cutoff frequency of said detectors, characterized in that the pixels thus obtained are agglomerated by group of at least four neighboring pixels so as to carry out staggered sampling, the step between two lines of the sampling thus carried out corresponding to the pitch between two successive lines of the initial pixel matrix.
  • the sampling thus obtained from these different implementation variants are transmitted to the ground where a staggered interpolation is performed on the pixels received. Thanks to this agglomeration and staggered subsampling, it is possible to reconstruct an image on the ground without folding spectra.
  • FIG. 1 illustrates the principle of acquisition by sweeping
  • - Figure 2 illustrates a processing chain of images identified by push-broom scanning
  • FIG. 3 illustrates the acquisition according to the method described in FR 2 678 460;
  • Figure 4 represents a matrix of pixels obtained by implementing the acquisition of Figure
  • - Figure 5 is a representation in the Fourier plane of the contour lines of the Modulation Transfer Function with rectangular symmetry; - Figure 6 illustrates a possible implementation for the invention;
  • Figure 7 illustrates a pixel matrix acquired by implementing the method of Figure 6;
  • Figure 9 is a representation in the Fourier plane of the spectrum of an image acquired by implementing the method of Figure 6;
  • FIG. 6 illustrates the acquisition processing carried out by two linear bars Bil and B12 in accordance with a possible embodiment for the invention.
  • Bars Bil and B12 have an arrangement similar to that of bars Bl and B2: they are parallel and juxtaposed and are offset in the direction in which they extend, by a fraction of a pixel, for example by a half-pixel (that is to say of p / 2, where p is the pitch between two elementary detectors of the same strip).
  • the sampling time of each of the bars Bil, B12 is chosen equal to 1 / fc, where fc is the cut-off frequency corresponding to the FTM of the apparatus. In practice, the sampling time 1 / fc and the step p between detectors coincide.
  • the pixel matrix then obtained is that illustrated in FIG. 7, in which the pixels acquired on the strip Bil are referenced XB11, those acquired by the strip B12 being referenced by XB12.
  • pixels 0 are restored to the ground by an interpolation from the pixels XB11 and XB12 in their vicinity.
  • interpolation processing making it possible to reconstruct these missing pixels is for example a Fourier interpolation, conventionally known to those skilled in the art and which will not be described here in more detail.
  • Other interpolation treatments are of course possible.
  • the step p generally corresponding to 1 / fc, the acquisition carried out is therefore equivalent to a staggered sampling, also called centered hexagonal, with a step p.V2 of the same image rotated by 45 °.
  • the sampling proposed by the invention can also be obtained, as illustrated in FIG. 8, by means of a single bar 11 parallel to the zone which it observes and oriented in its plane of displacement. an angle ⁇ of 45 ° relative to its direction of movement D.
  • sampling time te is then chosen so as to correspond to p / V2, where p is the sampling step along the bar 1, the sampling thus generated corresponding to a square sampling of steps p / v ⁇ .
  • the spectral aliasing is outside the area occupied by the spectrum (FIG. 10b). It is the same in the case of the staggered processing proposed by the invention in which the aliasing occurs outside the zone Zi
  • the optimality property of the proposed sampling comes from the collapse of the FTM (therefore of the information) along the diagonal precisely in the direction where the staggered sampling has the coarsest geometric resolution. This allows deconvolution to be carried out under good conditions.
  • row and column resolutions are identical to those obtained with conventional sampling, that is to say by taking all the points which cancel each other out at the Shannon frequency.
  • FIG. 12 Another possible embodiment for the invention has been illustrated in FIG. 12.
  • a first sampling is carried out conventionally by means of a single strip of detectors and with a sampling time corresponding to 1 / fc.
  • the pixels X of three lines t lf t 2 , t 3 shown with this sampling have been shown in thin lines by means of a strip of CCD detectors.
  • the pixels X are agglomerated, in groups of at least four, while retaining the same sampling frequency.
  • each Ag pixel being the sum of four neighboring pixels X distributed in square.
  • the pixels Ag of the same line of the agglomerated frame are juxtaposed in a contiguous manner, the pixels
  • Ag of two successive lines of the agglomerated frame being offset by the width of a half pixel Ag and overlapping the height of a line of pixels X.

Abstract

Procédé d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins une barrette ou une matrice de détecteurs de type à transfert de charge défilant en regard d'une zone observée, dans lequel on réalise une acquisition d'une pluralité de pixels, que l'on échantillonne en quinconce, caractérisé en ce que les pas en ligne et en colonne de cet échantillonnage quinconce sont tels que le spectre de l'image est principalement situé dans une zone du plan de Fourier où la Fonction de Transfert de Modulation est significative et où le repliement du spectre est limité.

Description

PROCEDE D'ACQUISITION D'UNE IMAGE PAR BALAYAGE POUSSE-BALAI La présente invention est relative à un procédé d'acquisition d'une image d'observation satellitaire de la terre par balayage "Pousse-balai" (Push-broom selon la terminologie anglo-saxonne) à partir d'une barrette ou d'une matrice de détecteurs de type à transfert de charge (CCD) défilant en regard de la zone observée.
Le principe d'un balayage de type "pousse-balai" a été illustré sur la figure 1 dans le cas d'une barrette 1 de détecteurs.
La barrette 1 réalise l'observation successive, au cours du défilement du satellite qui la porte, de lignes LI, L2,..., Ln perpendiculaires au déplacement (flèche D). Une optique 2 d'instrumentation forme à chaque instant l'image d'une ligne de paysage sur la ligne des détecteurs, la barrette 1 étant placée dans le plan focal de l'optique 2, perpendiculairement au vecteur vitesse du satellite. Le paysage défile devant chaque détecteur qui intègre le flux lumineux pendant le temps de pose et le transforme en une charge électrique proportionnelle.
On a illustré sur la figure 2 une chaîne de traitement classique des images ainsi relevées.
Cette chaîne de traitement comporte de façon schématique une unité 3 de traitement et d'amplification en sortie des détecteurs de la barrette 1, un codeur 4 analogique-numérique recevant le signal en sortie de l'unité 3, des moyens 5 pour la transmission au sol à partir du satellite des images numériques ainsi relevées, une unité 6 pour la reconstitution au sol des images . L'unité 3 comporte en particulier un registre à décalage- dans lequel les informations intégrées et stockées dans chaque détecteur de la barrette 1 sous forme de charge sont transférées à l'issue d'un temps de pose. Ce registre assure ensuite un transfert électronique des charges, qui sont converties en une succession de tensions proportionnelles aux flux lumineux reçus et intégrés.
L'unité 6 au sol reconstruit les images en mettant notamment en oeuvre des traitements de déconvolution destinés à compenser les défauts instrumentaux, ainsi que le cas échéant des traitements d'interpolation pour reconstituer certains pixels de l'image.
Usuellement, le temps d'acquisition entre deux lignes successives, ou temps d'échantillonnage, est tel que le point au sol situé à la verticale du satellite s'est déplacé d'une distance égale aux dimensions d'un détecteur élémentaire projeté au sol.
En se ramenant dans le plan focal de l'instrument, on dispose alors d'un échantillonnage de même fréquence en ligne et en colonne, égale à l'inverse de la taille d'un détecteur élémentaire dans le plan focal.
Toutefois, cette approche ne tient pas compte de la Fonction de Transfert de Modulation (FTM) instrumentale .
On rappelle que la Fonction de Transfert de Modulation d'un système optique est une fonction de l'espace fréquentiel qui traduit l'aptitude du système à transmettre les diverses fréquences. Elle est caractéristique de la restitution des contrastes du paysage par le système. Dans le cas d'une acquisition de type pousse- balai, la Fonction de Transfert de Modulation dépend principalement de l'optique du système, de l'effet de filé (affaissement de contraste dû au déplacement), ainsi que des détecteurs. Les systèmes d'observation sont caractérisés par la fréquence de coupure fc, au delà de laquelle la FTM est négligeable. On peut associer une fréquence de coupure à chacun des effets contribuant à la Fonction de Transfert de Modulation globale (optique, filé et intégration sur la zone photosensible), la fréquence de coupure globale étant la plus petite des trois valeurs précédentes.
La Fonction de Transfert de modulation associée à l'intégration sur la zone photosensible s'annule pour une fréquence égale à l'inverse de la taille d'un détecteur élémentaire qui constitue en première approximation la fréquence de coupure correspondante.
En général, c'est l'effet d'intégration sur la zone photosensible qui impose la fréquence de coupure globale. Dans le cas d'une acquisition classique où la fréquence d'échantillonnage fe est égale à la fréquence de coupure fc, la condition de Shannon (fe > 2.fc) n'est donc pas respectée et il en résulte un fort repliement de spectre introduisant des artefacts et rendant difficile toute tentative de déconvolution et d'interpolation.
Un but de l'invention est de proposer un procédé de traitement qui permet :
- de limiter le repliement du spectre, afin d'autoriser au sol une restauration et/ou interpolation convenable de l'image ;
- de minimiser la redondance de l'information transmise au sol de façon à pouvoir utiliser au mieux les capacités de transmission de données des instruments.
Il a déjà été proposé dans la demande de brevet FR 2 *678 460 un procédé de traitement permettant d'enregistrer et de retransmettre des images avec un pas d'échantillonnage en ligne et en colonne deux fois plus petit que pour un instrument de balayage classique, et ce avec un télescope et des détecteurs élémentaires identiques.
Selon ce procédé, qui a été illustré sur la figure 3, on réalise un suréchantillonnage en :
- divisant le temps d'intégration et le temps d'échantillonnage par deux, la fréquence d'échantillonnage en colonne devenant alors égale au double de la fréquence de coupure ;
- superposant, pour les lignes, les relevés des deux barrettes Bl et B2 décalées l'une par rapport à l'autre d'une distance égale à la moitié de la taille d'un détecteur élémentaire.
La superposition des deux barrettes peut être obtenue par l'intermédiaire d'un diviseur optique de ligne ou en utilisant le défilement du satellite par décalage dans le champ de la deuxième barrette d'un nombre entier de pas d'échantillonnage dans le sens d'avancement et d'un demi pas dans le sens de la barrette. Le décalage des.deux grilles d'échantillonnage peut être différent du demi pixel, moyennant sa prise en compte dans le post traitement, la valeur optimale du décalage est toutefois la moitié de la taille d'un détecteur élémentaire en ligne et en colonne.
On a illustré sur la figure 4 une matrice de pixels obtenus avec un tel traitement. Les pixels XB1 sont ceux acquis avec la barrette Bl, les pixels XB2 ceux acquis avec la barrette B2.
Chacune des deux barrettes Bl et B2 génère donc une grille d'échantillonnage rectangulaire dont le pas le long des lignes vaut la taille d'un détecteur élémentaire et le pas le long des colonnes la moitié de cette valeur, les deux grilles étant décalées d'une distance égale à la moitié d'un détecteur élémentaire. Par entrelacement tel qu'indiqué sur la figure 4, on reconstitue bien une grille d'échantillonnage carrée dont le pas en ligne coïncide avec le pas en colonne et vaut la moitié de la taille d'un détecteur élémentaire.
On dispose ainsi d'un procédé respectant la condition d'échantillonnage de Shannon au prix d'un quadruplement du débit d'information.
On a illustré sur la figure 5 la représentation dans le plan de Fourier d'une Fonction de Transfert de Modulation (FTM) réaliste pour un instrument d'observation de type "PUSH-BROOM", utilisé selon le procédé d'acquisition décrit dans FR 2 678 460. Cette FTM s'exprime comme le produit d'une fonction de transfert ftml des fréquences lignes normalisées fx et d'une fonction de transfert ftmc des fréquences colonnes fy, ces deux fonctions étant données par les formules suivantes dans le cas de l'instrument SPOT 5 :
e(-3-43i fx) sin (2ττfx) ftml =
2πfx e(-3.03fy) sin (2rrfy) sin(7τfy) ftmc =
2ττ2fy2 Sur cette représentation, la FTM est normalisée à l à l'origine et les fréquences fx et fy sont normalisées par rapport au double de la fréquence d'échantillonnage classique et les différentes courbes représentent des courbes de niveaux par pas de 0,01 jusqu'à la valeur 0,1. On voit sur cette représentation que les Fonctions de Transfert de Modulation sont faibles sur une grande partie du spectre d'origine fréquentiel.
Les fonctions ftml et ftmc s'annulent pour fx = fy = 0.5, et la condition de Shannon est donc satisfaite. En fait, toute l'information significative est essentiellement contenue dans le triangle supérieur Is, dont les extrémités correspondent à l'origine du plan de Fourier et aux points (0 ; 0,5) et (0,5 ; 0) .
Cette zone est en effet tangentée par la courbe FTM de niveau 0,03, valeur en dessous de laquelle on considère que le signal est noyé dans le bruit. Au delà de cette courbe, c'est-à-dire en dehors de la zone Is précitée, le signal est inexploitable. Par conséquent, l'échantillonnage proposé dans la technique décrite dans la demande FR 2 678 460 n'est pas optimal. Il transmet inutilement une grande partie du plan fréquentiel . Un but de l'invention est donc de proposer un procédé d'acquisition qui pallie cet inconvénient : il permet de conserver un repliement spectral négligeable tout en optimisant le débit d'acquisition.
On se place dans tout le présent texte dans le cas où la Fonction de Transfert de Modulation est à symétrie rectangulaire .
L'invention propose un procédé d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins une barrette ou une matrice de détecteurs de type à transfert de charge défilant en regard d'une zone observée, dans lequel on réalise une acquisition d'une pluralité de pixels, que l'on échantillonne en quinconce, caractérisé en ce que les pas en ligne et en colonne de cet échantillonnage quinconce sont tels que le spectre de l'image est principalement situé dans une zone du plan de Fourier où la Fonction de Transfert de Modulation est significative et où le repliement du spectre est limité.
Par rapport au procédé de FR 2 678 460 décrit précédemment, on divise donc par deux le débit de télémesure correspondant à une même zone imagée au sol.
Plus particulièrement, dans un mode de mise en oeuvre avantageux, il est proposé un procédé d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins deux barrettes ou matrices de détecteurs de même pas, s'étendant de façon juxtaposée perpendiculairement à leurs directions de déplacement et d'un nombre entier plus une fraction de pas élémentaire parallèlement à leur direction de déplacement, lesdites barrettes ou matrices étant décalées l'une par rapport à l'autre d'une fraction de pas élémentaire, perpendiculairement à leur direction de déplacement, caractérisé en ce que les temps d'échantillonnage de ces deux barrettes sont égaux à l'inverse de la fréquence de coupure des détecteurs élémentaires.
Dans un autre mode de mise en oeuvre avantageux, il est proposé un procédé d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins une barrette ou matrices de détecteurs, dont le pas p correspond à l'inverse de la fréquence de coupure desdits détecteurs, caractérisé en ce que ladite barrette est orientée dans son plan de déplacement d'un angle α = 1/n par rapport à sa direction de déplacement, n étant un nombre entier, et en ce que le temps d'échantillonnage des détecteurs élémentaires correspond à un déplacement de p/Vl+n2. on génère alors une grille d'échantillonnage carrée dans la géométrie définie par la vitesse du point sub-satellite et sa perpendiculaire, le pas d'échantillonnage dans ces deux directions valant p/Vl+n2. On densifiê ainsi l'échantillonnage d'un facteur vl+n2 au prix d'une réduction de la fauchée du même facteur.
Dans un autre mode de mise en oeuvre avantageux encore, il est proposé un procédé d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'une seule barrette ou matrice de détecteurs, avec un temps d'échantillonnage correspondant à l'inverse de la fréquence de coupure desdits détecteurs, caractérisé en ce qu'on agglomère les pixels ainsi obtenus par groupe d'au moins quatre pixels voisins de façon à réaliser un échantillonnage quinconce, le pas entre deux lignes de l'échantillonnage ainsi réalisé correspondant au pas entre deux lignes successives de la matrice de pixels initiale. L'échantillonnage ainsi obtenu à partir de ces différentes variantes de mise en oeuvre est transmis au sol où l 'on réalise au sol une interpolation quinconce sur les pixels reçus . Grâce à cette agglomération et au sous- échantillonnage en quinconce, on peut reconstituer au sol une image sans repliement de spectres .
D ' autres caractéristiques et avantages de l ' invention ressortiront encore de la description qui suit . Cette description est purement illustrative et non limitative . Elle doit être lue en regard des dessins annexés sur lesquels seul le quart du spectre est représenté , le complément s ' en déduisant par des considérations de symétrie : la figure 1 illustre le principe d ' une acquisition par balayage pousse-balai ; - la figure 2 illustre une chaîne de traitement d' images relevées par balayage pousse-balai ;
- la figure 3 illustre l ' acquisition selon le procédé décrit dans FR 2 678 460 ;
- la figure 4 représente une matrice de pixels obtenue en mettant en oeuvre l ' acquisition de la figure
3 ;
- la figure 5 est une représentation dans le plan de Fourier des courbes de niveau de la Fonction de Transfert de Modulation à symétrie rectangulaire ; - la figure 6 illustre un mode de mise en oeuvre possible pour l'invention ;
- la figure 7 illustre une matrice de pixels acquise en mettant en oeuvre le procédé de la figure 6 ;
- la figure 8 illustre un autre mode de mise en oeuvre possible pour l'invention ;
-" la figure 9 est une représentation dans le plan de Fourier du spectre d'une image acquise en mettant en oeuvre le procédé de la figure 6 ;
- les figures 10a, 10b, 10c illustrent les repliements de spectre pour chacune de ces acquisitions ; - la figure 11 illustre un autre mode de mise en oeuvre de l'invention semblable à celui de la figure 8 ;
- la figure 12 illustre un autre mode de traitement possible pour l'invention. On a illustré sur la figure 6 le traitement d'acquisition que réalisent deux barrettes linéaires Bil et B12 conformément à un mode de mise en oeuvre possible pour l'invention.
Les barrettes Bil et B12 présentent une disposition semblable à celle des barrettes Bl et B2 : elles sont parallèles et juxtaposées et sont décalées dans la direction selon laquelle elles s'étendent, d'une fraction de pixel, par exemple d'un demi-pixel (c'est-à- dire de p/2, où p est le pas entre deux détecteurs élémentaires d'une même barrette).
Le temps d'échantillonnage de chacune des barrettes Bil, B12 est choisi égal à 1/fc, où fc est la fréquence de coupure correspondant à la FTM de 1'appareillage. Dans la pratique, le temps d'échantillonnage 1/fc et le pas p entre détecteurs coïncident.
La matrice de pixels alors obtenue est celle illustrée sur la figure 7, sur laquelle les pixels acquis sur la barrette Bil sont référencés XB11, ceux acquis par la barrette B12 étant référencés par XB12.
Ainsi qu'on peut le voir sur la figure 7, on dispose ainsi de deux images correspondant aux pas d'acquisition classiques de 1/fc et p, mais décalées l'une par rapport à l'autre de l/2fc en ligne et de p/2 en colonne.
Par conséquent, par rapport à l'acquisition décrite en référence aux figures 3 et 4, on a enregistré un pixel sur deux en ligne et en colonne, c'est-à-dire que l'on a réalisé un échantillonnage quinconce. Sur la matrice de la figure 7, les pixels de la matrice de la figure 4 qui n'ont pas été acquis ont été remplacés par des 0.
Ces pixels 0 sont restitués au sol par une interpolation à partir des pixels XB11 et XB12 de leur voisinage.
Le traitement d'interpolation permettant de reconstituer ces pixels manquant est par exemple une interpolation de Fourier, classiquement connue de l'Homme du Métier et qui ne sera pas ici décrite plus en détails. D'autres traitements d'interpolation sont bien entendu possibles. A titre d'exemple, on pourra également utiliser un traitement d'interpolation du type de celui décrit dans la publication : "Analyse multirésolution pour les images avec un facteur de résolution V2" - J.C. FEAUVEAU - Traitement du signal - Volume 7 - n°2 - 1990.
Le pas p correspondant généralement à 1/fc, l'acquisition réalisée est donc équivalente à un échantillonnage quinconce, également appelée hexagonal centré, avec un pas p.V2 de la même image tournée de 45°.
Par conséquent, l'échantillonnage proposé par l'invention peut également être obtenu, ainsi qu'illustré sur la figure 8, au moyen d'une seule barrette 11 parallèle à la zone qu'elle observe et orientée dans son plan de déplacement d'un angle α de 45° par rapport à sa direction de déplacement D.
Le temps te d'échantillonnage selon la direction de déplacement est alors choisi de façon à correspondre à p/V2, où p est le pas d'échantillonnage le long de la barrette 1, l'échantillonnage ainsi généré correspondant à un échantillonnage carré de pas p/v^.
Le spectre d'un tel échantillonnage a été illustré sur la figure 9 et correspond à celui de la figure 5 tourné de 45°. On voit sur cette figure que la zone du plan de Fourier occupée par le spectre de l'image est le triangle tangent à la courbe FTM0,o3«
On se réfère maintenant aux figures 10a à 10c. Les zones grisées représentées sur ces figures sont les zones du plan de Fourier où le repliement de spectre est supérieur à 3 %, pour chacun des trois modes d'acquisition.
On voit sur la figure 10a que dans le cas d'une acquisition classique, seule une petite partie de l'information est peu polluée (moins de 15 %).
Dans le cas d'une acquisition selon le traitement décrit dans FR 2 678 460, les repliements spectraux sont en dehors de la zone occupée par le spectre (figure 10b). II en est de même dans le cas du traitement quinconce proposé par l'invention dans lequel le repliement de spectre intervient en dehors de la zone Zi
(figure 10c).
Ces résultats sont par ailleurs indiqués dans le tableau I ci-dessous sur lequel on a porté :
- la zone du plan de Fourier polluée par le repliement de spectres, exprimée en fraction de la surface de la zone a priori accessible compte tenu de 1'échantillonnage, - la zone du plan de Fourier pour lequel le contenu fréquentiel n'est pas trop atténué par la fonction de transfert, caractérisée par une valeur de FTM supérieure au seuil de 3% ;
- la quantité de données produites pour imager une même zone, normalisée par le débit "classique". TABLEAU I
Figure imgf000014_0001
On voit sur ce tableau que le traitement décrit dans FR 2 678 460 ne présente pas de repliement spectral mais génère des spectres en grande partie vides. A l'inverse, le traitement classique utilise une zone du plan de Fourier qui est trop restreinte pour contenir toute l'information transmise par la FTM instrumentale, ce qui induit un fort niveau de repliement des spectres et interdit tout déconvolution. L'échantillonnage proposé par l'invention permet quant à lui d'obtenir un spectre qui occupe une zone où la Fonction de Transfert de Modulation est suffisante pour que le signal soit exploitable et où le repliement du spectre est négligeable. Le débit de données est deux fois moins important qu'avec l'échantillonnage de FR 2 678 460.
La propriété d'optimalité de l'échantillonnage proposé provient de l'affaissement de la FTM (donc de l'information) suivant la diagonale précisément dans la direction où l'échantillonnage quinconce possède la résolution géométrique la plus grossière. Ceci permet d'effectuer une déconvolution dans de bonnes conditions.
A l'inverse, les résolutions ligne et colonne sont identiques à celles que l'on obtient avec un échantillonnage classique, c'est-à-dire en prenant tous les points qui s'annulent à la fréquence Shannon.
De façon plus générale, ainsi qu'illustré sur la figure 11, une FTM significative et un repliement de spectre limité sont avantageusement obtenus au moyen d'une barrette de détecteurs parallèle à la zone qu'elle observe et orientée dans son plan de déplacement d'un angle α par rapport à sa direction de déplacement, tel que : Tg α = 1/n, avec n entier ; le temps d'échantillonnage correspond à un déplacement de
Figure imgf000015_0001
Le suréchantillonnage illustré sur la figure 10 correspond au cas n = 2. Le cas particulier où α = 45° et le temps d'échantillonnage est égale à p/V2 correspond à n = 1.
Bien entendu, on se place dans ces différents cas dans l'hypothèse où p correspond à 1/fc.
On a illustré sur la figure 12 un autre mode de mise en oeuvre possible pour l'invention.
Dans ce mode de mise en oeuvre, un premier échantillonnage est réalisé de façon classique au moyen d'une seule barrette de détecteurs et avec un temps d'échantillonnage correspondant à 1/fc. Sur la figure 12 précitée, on a représenté en traits fins les pixels X de trois lignes tlf t2, t3 relevées avec cet échantillonnage au moyen d'une barrette de détecteurs CCD.
Selon le traitement conforme à cette mise en oeuvre de l'invention, on agglomère les pixels X, par groupe d'au moins quatre, tout en conservant la même fréquence d'échantillonnage.
Les pixels Ag ainsi obtenus sont représentés en traits gras sur la figure 12, chaque pixel Ag étant la somme de quatre pixels X voisins répartis en carré.
Les pixels Ag d'une même ligne de la trame agglomérée sont juxtaposés de façon jointive, les pixels
Ag de deux lignes successives de la trame agglomérée étant décalés de la largeur d'un demi-pixel Ag et se recouvrant sur la hauteur d'une ligne de pixels X.
On dispose ainsi d'un suréchantillonnage quinconce analogue à celui de la figure 7 :
- deux lignes successives sont séparées de 1/fc, c'est-à- dire de deux fois 1/Fc où Fc correspond à la fréquence de coupure de quatre détecteurs agglomérés de la barrette utilisée pour l'acquisition initiale,
- deux pixels Ag d'une même ligne sont séparés d'une distance correspondant à la largeur élémentaire desdits pixels Ag.
A la réception , on réalise s ur ce suréchantillonnage un traitement d 'interpolation quinconce et une déconvolution.
Avec une tel traitement , il est possible de restituer quasiment la résolution initiale de l ' image , alors que le débit est divisé par deux.

Claims

REVENDICATIONS
1. Procédé d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins une barrette ou une matrice de détecteurs de type à transfert de charge défilant en regard d'une zone observée, dans lequel on réalise une acquisition d'une pluralité de pixels, que l'on échantillonne en quinconce, caractérisé en ce que les pas en ligne et en colonne de cet échantillonnage quinconce sont tels que le spectre de l'image est principalement situé dans une zone du plan de Fourier où la Fonction de Transfert de Modulation est significative et où le repliement du spectre est limité.
2. Procédé selon la revendication 1 d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins deux barrettes ou matrices de détecteurs de même pas, s'étendant de façon juxtaposée perpendiculairement à leurs directions de déplacement et décalées d'un nombre entier plus une fraction de pas élémentaire parallèlement à leur direction de déplacement, lesdites barrettes ou matrices étant décalées l'une par rapport à l'autre d'une fraction de pas élémentaire, perpendiculairement à leur direction de déplacement, caractérisé en ce que les temps d'échantillonnage de ces deux barrettes sont égaux à l'inverse de la fréquence de coupure des détecteurs élémentaires.
3. Procédé selon la revendication 1 d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'une seule barrette ou matrice de détecteurs, dont le pas p correspond à l'inverse de la fréquence de coupure desdits détecteurs, caractérisé en ce que ladite barrette est orientée dans son plan de déplacement d'un angle α = 1/n par rapport à sa direction de déplacement, n étant un nombre entier, et en ce que le temps d'échantillonnage des détecteurs élémentaires est de p/Vl+n2.
4. Procédé selon la revendication 3 , caractérisé en ce que l ' angle entre la barrette ou matrice de détecteurs et la direction de déplacement est de 45°, le temps d'échantillonnage des détecteurs élémentaires étant de ρ/V2.
5. Procédé selon la revendication 1 d'acquisition et de traitement d'une image d'observation satellitaire de la terre au moyen d'au moins une barrette ou matrice de détecteurs, avec un temps d'échantillonnage correspondant à l'inverse de la fréquence de coupure desdits détecteurs, caractérisé en ce qu'on agglomère les pixels ainsi obtenus par groupe d'au moins quatre pixels voisins de façon à réaliser un échantillonnage quinconce, le pas entre deux lignes de l'échantillonnage ainsi réalisé correspondant au pas entre deux lignes successives de la matrice de pixels initiale.
6. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'échantillonnage ainsi obtenu est transmis au sol et en ce qu'on réalise au sol une interpolation quinconce sur les pixels reçus.
PCT/FR1996/001188 1995-07-28 1996-07-26 Procede d'acquisition d'une image par balayage pousse-balai WO1997005451A1 (fr)

Priority Applications (4)

Application Number Priority Date Filing Date Title
DE69600785T DE69600785T2 (de) 1995-07-28 1996-07-26 Bildaufnahmeverfahren mit "push-broom" abtastung
EP96927095A EP0782695B1 (fr) 1995-07-28 1996-07-26 Procede d'acquisition d'une image par balayage pousse-balai
JP9507289A JP3012865B2 (ja) 1995-07-28 1996-07-26 プッシュブルーム走査を用いた画像取り込み方法
US08/849,902 US6215522B1 (en) 1995-07-28 1996-07-26 Method of acquiring an image by push-broom scanning

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
FR95/09241 1995-07-28
FR9509242A FR2737375B1 (fr) 1995-07-28 1995-07-28 Procede d'acquisition d'une image par balayage pousse-balai
FR95/09263 1995-07-28
FR95/09242 1995-07-28
FR9509241A FR2737326B1 (fr) 1995-07-28 1995-07-28 Procede pour le traitement ou la transmission d'une image acquise sous la forme d'une matrice de pixels
FR9509263A FR2737376B1 (fr) 1995-07-28 1995-07-28 Procede et dispositif pour l'acquisition d'une image par echantillonnage par une barrette ou matrice de detecteurs optiques elementaires a transfert de charge

Publications (1)

Publication Number Publication Date
WO1997005451A1 true WO1997005451A1 (fr) 1997-02-13

Family

ID=27253082

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR1996/001188 WO1997005451A1 (fr) 1995-07-28 1996-07-26 Procede d'acquisition d'une image par balayage pousse-balai

Country Status (7)

Country Link
US (1) US6215522B1 (fr)
EP (1) EP0782695B1 (fr)
JP (1) JP3012865B2 (fr)
AT (1) ATE172295T1 (fr)
DE (1) DE69600785T2 (fr)
ES (1) ES2124636T3 (fr)
WO (1) WO1997005451A1 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6493021B1 (en) 1998-03-20 2002-12-10 Centre National D'etudes Spatiales Image acquisition by push-broom scanning
EP2012526A1 (fr) 2007-07-06 2009-01-07 Thales Procédé et dispositif d'observation

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2763199B1 (fr) * 1997-05-12 1999-07-30 Centre Nat Etd Spatiales Procede pour l'acquisition d'une image par balayage pousse-balai surechantillonne
JP4097966B2 (ja) * 2002-03-22 2008-06-11 オリンパス株式会社 画像取得装置
US7381211B2 (en) * 2002-08-21 2008-06-03 Kci Licensing, Inc. Medical closure screen device and method
FR2899344B1 (fr) * 2006-04-03 2008-08-15 Eads Astrium Sas Soc Par Actio Procede de restitution de mouvements de la ligne de visee d'un instrument optique
CN101257034B (zh) * 2008-04-10 2011-05-18 中国科学院长春光学精密机械与物理研究所 可提高分辨率的ccd像元
FR3005226B1 (fr) * 2013-04-25 2015-04-10 Astrium Sas Saisie d'image avec addition de signaux d'accumulation pour des photodetecteurs adjacents
FR3015740B1 (fr) * 2013-12-20 2016-02-05 Centre Nat Etd Spatiales Procede, dispositif et systeme de traitement de correction d'images

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0368354A2 (fr) * 1988-11-10 1990-05-16 Sony Corporation Circuit de traitement de signal par un dispositif de prise d'image à l'état solide
WO1991003125A1 (fr) * 1989-08-25 1991-03-07 Thomson Consumer Electronics S.A. Camera et dispositif de visualisation
US5067016A (en) * 1990-03-07 1991-11-19 Industrial Technology Research Institute Progressive scan system with field-offset sampling
FR2678460A1 (fr) * 1991-06-27 1992-12-31 Centre Nat Etd Spatiales Procede d'enregistrement et de traitement d'une image et instrument de prise de vues pour la mise en óoeuvre de ce procede.
DE4307833A1 (de) * 1993-03-12 1994-09-15 Thomson Video Europ Gmbh Videokamera mit einem Halbleitersensor

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4683496A (en) 1985-08-23 1987-07-28 The Analytic Sciences Corporation System for and method of enhancing images using multiband information
CA2021052C (fr) 1990-07-12 1995-08-29 Clifford D. Anger Imageur spectrographique a barrette de capteurs
US5307427A (en) 1991-09-04 1994-04-26 Westinghouse Electric Corp. Non-integral pixel compression for second generation FLIR sensors
US5251037A (en) 1992-02-18 1993-10-05 Hughes Training, Inc. Method and apparatus for generating high resolution CCD camera images
NL9300312A (nl) 1993-02-19 1994-09-16 Hollandse Signaalapparaten Bv Optische aftastinrichting.
US5692062A (en) * 1994-10-03 1997-11-25 Recon/Optical, Inc. Electro-optical imaging array with profiled foward motion compensation
US5668593A (en) * 1995-06-07 1997-09-16 Recon/Optical, Inc. Method and camera system for step frame reconnaissance with motion compensation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0368354A2 (fr) * 1988-11-10 1990-05-16 Sony Corporation Circuit de traitement de signal par un dispositif de prise d'image à l'état solide
WO1991003125A1 (fr) * 1989-08-25 1991-03-07 Thomson Consumer Electronics S.A. Camera et dispositif de visualisation
US5067016A (en) * 1990-03-07 1991-11-19 Industrial Technology Research Institute Progressive scan system with field-offset sampling
FR2678460A1 (fr) * 1991-06-27 1992-12-31 Centre Nat Etd Spatiales Procede d'enregistrement et de traitement d'une image et instrument de prise de vues pour la mise en óoeuvre de ce procede.
DE4307833A1 (de) * 1993-03-12 1994-09-15 Thomson Video Europ Gmbh Videokamera mit einem Halbleitersensor

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6493021B1 (en) 1998-03-20 2002-12-10 Centre National D'etudes Spatiales Image acquisition by push-broom scanning
EP2012526A1 (fr) 2007-07-06 2009-01-07 Thales Procédé et dispositif d'observation

Also Published As

Publication number Publication date
DE69600785T2 (de) 1999-06-10
EP0782695B1 (fr) 1998-10-14
EP0782695A1 (fr) 1997-07-09
ATE172295T1 (de) 1998-10-15
JP3012865B2 (ja) 2000-02-28
US6215522B1 (en) 2001-04-10
DE69600785D1 (de) 1998-11-19
ES2124636T3 (es) 1999-02-01
JPH10508172A (ja) 1998-08-04

Similar Documents

Publication Publication Date Title
EP0782695B1 (fr) Procede d'acquisition d'une image par balayage pousse-balai
WO2011138541A1 (fr) Procede d'imagerie polychrome
FR2779852A1 (fr) Appareil de multiplication de lignes video procurant une finesse de details amelioree, et procede associe
EP0983684B9 (fr) Perfectionnement a l'acquisition d'image par balayage pousse-balai
FR2513014A1 (fr) Dispositif de detection de radiations avec dispositifs de lecture et de transfert ainsi que dispositif et procede de production d'indications electriques liees au balayage d'image
WO2009030698A1 (fr) Procédé d'augmentation de résolution d'images multi-spectrales
EP2567538A1 (fr) Procede d'imagerie
EP0604302B1 (fr) Procédé et dispositif d'analyse, par radiograhpie, d'objects en défilement
WO2015091518A1 (fr) Procédé, dispositif et système de traitement de correction d'images
FR2585527A1 (fr) Dispositif pour la restitution et/ou l'analyse d'images en couleurs utilisant un tube cathodique ligne a ecran sur fibres optiques
FR2737375A1 (fr) Procede d'acquisition d'une image par balayage pousse-balai
Reichenbach et al. Restoration and reconstruction from overlapping images for multi-image fusion
EP3777129B1 (fr) Capteur d'images aeroporte realisant des prises de vue matricielle par decalage temporel et sommation multi-spectrales
FR3024782A1 (fr) Systeme et procede d'imagerie infrarouge
FR3032543A1 (fr) Procede de determination d'un axe de rotation d'un objet en tomographie et procede de caracterisation par tomographie
Zhang et al. Spatio-spectral modulation using a binary photomask for compressive chromotomography
EP2002226B1 (fr) Dispositif de formation de spectre sur un capteur optique a rejet spatial
FR2737376A1 (fr) Procede et dispositif pour l'acquisition d'une image par echantillonnage par une barrette ou matrice de detecteurs optiques elementaires a transfert de charge
WO2023046551A1 (fr) Procede pour l'acquisition d'images multispectrales et d'imagettes panchromatiques
FR2737326A1 (fr) Procede pour le traitement ou la transmission d'une image acquise sous la forme d'une matrice de pixels
Omrane Super-resolution for unregistered satellite images
FR3104744A1 (fr) Instrument d'observation optique embarque a resolution spatiale et spectrale variables
FR3066268A1 (fr) Capteur d'images
FR2683928A1 (fr) Procede et dispositif de formation d'images par traitement des donnees fournies par des capteurs optroniques.

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): JP US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 1996927095

Country of ref document: EP

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 08849902

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 1996927095

Country of ref document: EP

WWG Wipo information: grant in national office

Ref document number: 1996927095

Country of ref document: EP