FR 2851676 A1 20040827 FR 0450325 A 20040220 PROCÉDÉ D'OBTENTION D'UNE IMAGE A RÉSOLUTION AUGMENTÉE PAR L'UTILISATION D'UNE PLURALITÉ D'IMAGES A FAIBLE RÉSOLUTION La présente invention porte sur un procédé de traitement d'images et, plus particulièrement, sur un procédé d'obtention d'une image à résolution augmentée par l'utilisation d'une pluralité d'images à faible résolution. En ce qui concerne l'arrière-plan technologique du procédé de l'invention, nous nous référons aux documents suivants : [1] M. Irani et S. Peleg "Améliorer la Résolution par Enregistrement d'Images (Improving Resolution by Image Registration)" CVGIP . Graphical Models and Image Proc., 1991, Vol. 53, pages 231-239 ; [2] R.Y. Tsai et T.S. Huang "Restauration et Enregistrement d'Images Multicadres (Multiframe Image Restoration and Registration)" dans Advances in Computer Vision and Image Processing, Vol. 1 (T.S. Huang, ed.), Greenwich, CT: Jai Press, 1984, pages 317-339 ; [3] P. Cheeseman, B. Kanefsky, R. Kruft, J. Stutz et R. Hanson "Reconstruction de Surface Super-Résolue à partir d'Images Multiples (Super-Resolved Surface Reconstruction from Multiple Images)", NASA Technical Report FIA-94-12,1994 ; [4] A.M. Tekalp, M. K. Ozkan et M. I. Sezan "Reconstruction d'Images à Haute Résolution Pour des Séquences d'Images à Résolution Inférieure et Restauration d'Images Variant dans l'Espace (High-Resolution Image Reconstruction for Lower-Resolution Image Sequences and Space-Varying Image Restoration) ", IEEE International Conference on Acoustics, Speech, and Signal Processing, San Francisco, CA, 1992, Vol. III, pages 169-172 ; [5] R. Gonzalez et R. E. Woods, Digital Image Processing, Addison-Wesley, Reading, MA, 1992 ; [6] W. K. Pratt, Digital Image Processing, 2ème Ed., Wiley, New York, 2001 ; [7] Brevet américain US NO. 6 330 344 ; [8] Brevet américain US NO. 5 969 848. En raison des contraintes de l'environnement et de la résolution des détecteurs d'images, parfois nous pouvons seulement obtenir des images à faible résolution. Pour améliorer la qualité et la résolution d'images vues par les yeux humains, il faut plus d'une image d'entrée. Avec des séquences d'images, une scène floue, une figure indistincte ou un objet peu clair de qualité médiocre peut être reconstruit en une image de sortie à résolution augmentée et peut alors être facilement observé et reconnu. Les recherches antérieures concernant la reconstruction d'une image à résolution augmentée par l'utilisation d'images à faible résolution sont principalement divisées en méthodes itératives [1], méthodes à domaine de fréquence [2] et méthodes statistiques de Bayesian [3]. Dans les méthodes précitées, jusqu'ici, l'algorithme itératif qui a été développé par Irani (1] en 1991 et reconstruit une image à résolution augmentée principalement par enregistrement d'images, est encore le plus fiable dans les techniques antérieures concernant la reconstruction d'une image à résolution augmentée. La méthode itérative consiste principalement en trois phases : estimation initiale, procédé d'imagerie et procédé de reconstruction. Les procédures des trois phases de la méthode itérative seront décrites en détail dans la description détaillée de l'invention. Cependant, il est souligné que la méthode itérative consommera davantage de temps de calcul alors que le facteur de grossissement prédéterminé dans la méthode itérative deviendra plus grand, c'est-à-dire que la dimension de l'image reconstruite deviendra plus grande. Typiquement, le temps d'exécution de la reconstruction d'image par la méthode itérative est de l'ordre d'heures et dépend de la performance du système d'ordinateur. Par conséquent, un objectif de l'invention est de proposer un procédé pour reconstruire une image à résolution augmentée avec des algorithmes perfectionnés de résolution augmentée, qui est basé sur la méthode itérative d'Irani et emploie une interpolation initiale, une sélection d'image automatique et un enregistrement d'image robuste que nous proposons. De plus, l'image à résolution augmentée reconstruite par le procédé de l'invention peut être améliorée par un post-procédé d'amélioration d'image pour augmenter la qualité d'image de celle-ci. Alors que les systèmes classiques de la reconstruction d'images à résolution augmentée fonctionnent bien seulement alors que les séquences d'images à faible résolution sont capturées par le déplacement d'une caméra stationnaire dans un déplacement constant en relation avec la scène entière, c'est-à-dire que les cibles, qui ont besoin d'être reconstruites, sont associées à des scènes stationnaires, un autre objectif de l'invention est de proposer un procédé pour la reconstruction d'une image à résolution augmentée, qui fonctionne bien non seulement pour les conditions alors que les cibles, qui ont besoin d'être reconstruites, sont associées à des scènes stationnaires, mais encore pour les conditions alors que les cibles, qui ont besoin d'être reconstruites, sont associées à des objets en mouvement. Un objectif de l'invention est de proposer un procédé de reconstruction d'une image à résolution augmentée avec des algorithmes perfectionnés de résolution augmentée. De plus, le procédé de l'invention, basé sur la méthode itérative, utilise une interpolation initiale, une sélection d'image automatique et un enregistrement d'image robuste, que nous proposons. De plus, l'image à résolution augmentée reconstruite par le procédé de l'invention peut être améliorée par un post-procédé d'amélioration d'image pour augmenter la qualité d'image de celle-ci. Un autre objectif de l'invention est de proposer un procédé pour reconstruire une image à résolution augmentée, qui fonctionne bien non seulement pour les conditions alors que les cibles, qui ont besoin d'être reconstruites, sont associées à des scènes stationnaires, mais encore pour les conditions alors que les cibles, qui ont besoin d'être reconstruites, sont associées à des objets en mouvement. La présente invention a donc pour objet un procédé d'obtention d'une image à résolution augmentée, caractérisé par le fait qu'il comprend l'utilisation de M premières images à faible résolution associées à une scène, chacune des M premières images étant représentée par un ensemble de pixels uniformément échantillonnés et quantifiés, une relation d'un facteur de grossissement existant entre l'image à résolution augmentée et les M premières images, ledit procédé comprenant les étapes consistant à : (a) stocker les M premières images ; (b) sélectionner l'une parmi les M premières images en tant qu'image prototype, et désigner les (M-1) premières images non-sélectionnées comme étant les (M- 1) secondes images ; (c) sur la base du facteur de grossissement, interpoler des pixels supplémentaires dans l'ensemble de pixels de l'image prototype, puis déduire les valeurs des pixels supplémentaires conformément aux valeurs de pixels voisins de ceux-ci pour produire une image prototype interpolée, et calculer une translation respective entre chacune des (M-l) secondes images et l'image prototype interpolée ; (d) diviser les translations des (M-1) secondes images par le facteur de grossissement, respectivement, afin d'obtenir le module en ce qui concerne les translations des (M-l) secondes images, sélectionner, sur la base d'un critère, l'une parmi les secondes images dont les modules apparentés sont les mêmes, et désigner les secondes images sélectionnées conjointement avec le reste des secondes images dont les modules apparentés ne sont pas les mêmes comme étant les N troisièmes images, où N est égal à ou inférieur à (M-l) ; (e) sous-échantillonner l'image prototype interpolée N fois conformément à la translation respective entre chacune des N troisièmes images et l'image prototype interpolée pour produire N quatrièmes images qui correspondent chacune à l'une des N troisièmes images ; (f) calculer la différence entre chacune des N troisièmes images et la quatrième image correspondante de celle- ci ; (g) ajuster les valeurs des pixels de l'image prototype interpolée conformément à une moyenne des différences calculées à l'étape (f) ; et (h) répéter l'étape (e) à l'étape (g) jusqu' à ce que les valeurs des pixels de l'image prototype interpolée convergent vers un résultat satisfaisant, et désigner l'image prototype interpolée dont les valeurs de pixels convergent vers le résultat satisfaisant comme étant l'image à résolution augmentée. L'étape (e) peut être effectuée conformément à la formule suivante : O kn) \Tk l.l (n) l hl L S , dans laquelle : gk(n) est la kème quatrième image au nème procédé de sous- échantillonnage ; (ni est l'image prototype interpolée étant ajustée n fois ; - h est l'opérateur de création de flou défini par une fonction de diffusion de points ; - Tk est l'opérateur de transformation en ce qui concerne la translation entre la kème troisième image et l'image prototype interpolée ; - #s est l'opérateur de sous-échantillonnage. Les étapes (f) et (g) peuvent être effectuées conformément à la formule suivante : /(n+l) = /(11) + f,Tk-l (((gk - g/n ) Î s). p), dans laquelle : f(n+l) est l'image prototype interpolée pour être ajustée à (n+1) fois ; - K est le nombre des troisièmes images étant calculé ; gk est la kème troisième image ; p est l'opérateur de retrait du flou ; Tk-1 est l'opérateur de transformation inverse en ce qui concerne la translation entre la kème troisième image et l'image prototype interpolée ; Ts est l'opérateur de sur-échantillonnage. A l'étape (c), la déduction des valeurs des pixels supplémentaires peut être basée sur l'une choisie dans le groupe consistant en une interpolation de premier ordre, une interpolation de troisième ordre et une interpolation de cinquième ordre. A l'étape (c) , la translation (LTi) entre la ième seconde image et l'image prototype interpolée peut être calculée conformément aux étapes suivantes : (c-1) assigner au moins un point intéressant P1 (x,y) sur la ième seconde image, chaque point intéressant correspondant à une fenêtre locale (w) et à un point (u,v) sur l'image prototype interpolée ; (c-2) calculer une différence absolue (LADi(x,y;u,v)) entre chaque point intéressant Pi(x,y) et le point correspondant(u,v) de celui-ci conformément à la formule : LADt(x,y;u,v)= LII(x+m,y+n)-Io(u+m,v+n)1 (m,n)ew (c-3) calculer un minimum (LRi(x,y)) de la différence absolue (LADi(x,y ;u,v) calculée à l'étape (c-2) conformément à la formule suivante : LR(x,y) = a,rg min LAD, (x, y; u, v) (u,v) (c-4) calculer la translation (LTi) en ce qui concerne chaque point intéressant Pi(x,y) conformément à la formule suivante : LTi (x,y) = LRi(x,Y)-(X,Y)8 facteur de grossissement. A l'étape (c), la translation (LTi) entre la ième seconde image et l'image prototype interpolée peut être calculée conformément aux étapes suivantes : (c-1) calculer une différence absolue (GAD(u,v)) en ce qui concerne la ième seconde image entière conformément à la formule suivante : GAD(u,v)= LII(x,y) -Io(u + x, v+ y)1 (x,y)E dans laquelle le point (u,v) est un point, sur l'image prototype interpolée, correspondant à un point de départ de la ième seconde image ; (c-2) calculer un minimum (GR(i)) de la différence absolue (GAD(u,v)) calculée à l'étape (c-1) conformément à la formule suivante : GR(i) = arg min GAD, (u,v) (u,v) et considérer GR(i) comme LT1. Conformément à un mode de réalisation préféré de l'invention, M premières images à faible résolution sont utilisées pour reconstruire une image à résolution augmentée. Les M premières images à faible résolution sont associées à une scène. Chacune des M premières images est représentée par un ensemble de pixels uniformément échantillonnés et quantifiés. Une relation d'un facteur de grossissement existe entre l'image à résolution augmentée et chacune des M premières images. Tout d'abord, les M premières images sont stockées. Ensuite, l'une est sélectionnée parmi les M premières images en tant qu'image prototype, et les (M-l) premières images non sélectionnées sont désignées comme étant les (M-l) secondes images. Sur la base du facteur de grossissement, des pixels supplémentaires sont interpolés dans l'ensemble de pixels de l'image prototype, puis chaque valeur de pixels supplémentaires est calculée conformément aux valeurs de pixels voisins de ceux-ci pour produire une image prototype interpolée. Les translations respectives, qui existent chacune entre chacune des (M-l) secondes images et l'image prototype interpolée, sont également calculées. Les translations des (M-l) secondes images sont divisées par le facteur de grossissement, respectivement, pour obtenir le module en ce qui concerne les translations des (M-l) secondes images. Sur la base d'un critère, l'une est sélectionnée à partir des secondes images dont les modules apparentés sont les mêmes. Les secondes images sélectionnées conjointement avec le reste des secondes images, dont les modules apparentés ne sont pas les mêmes, sont désignées comme étant les N troisièmes images où N est égal à ou inférieur à (M-1). L'image prototype interpolée est sous-échantillonnée N fois conformément aux translations respectives entre chacune des N troisièmes images et l'image prototype interpolée pour produire N quatrièmes images qui correspondent chacune à l'une des N troisièmes images. La différence entre chacune des N troisièmes images et la quatrième image correspondante de celle-ci est calculée. Les valeurs des pixels de l'image prototype interpolée sont ajustées conformément à une moyenne des différences entre chacune des N troisièmes images et la quatrième image correspondante de celle-ci. Les étapes de sous-échantillonnage et d'ajustement des valeurs des pixels de l'image prototype interpolée sont répétées jusqu'à ce que les valeurs des pixels de l'image prototype interpolée convergent vers un résultat satisfaisant. Finalement, l'image prototype interpolée dont les valeurs de pixels convergent vers le résultat satisfaisant est désignée comme étant l'image à résolution augmentée. Les avantages et l'esprit de l'invention peuvent être compris par les éléments suivants conjointement avec les dessins annexés. BRÈVE DESCRIPTION DES DESSINS ANNEXÉS - la FIGURE 1 est un organigramme montrant les phases de travail du procédé selon la présente invention ; - la FIGURE 2A et la FIGURE 2B sont des diagrammes schématiques montrant le calcul de la valeur P par l'utilisation d'un mode d'interpolation du troisième ordre ; - la FIGURE 3 évalue la performance de modes d'interpolation de différents ordres par le Rapport Signal de Pic à Bruit (Peak Signal-to-Noise Ratio - PSNR) entre l'image initiale et les images reconstruites ; - la FIGURE 4A montre schématiquement l'enregistrement d'image à l'aide d'une technique d'appariement locale ; - la FIGURE 4B montre schématiquement l'enregistrement d'image à l'aide d'une technique d'appariement globale ; les FIGURES 5A à 5D montrent les résultats sur le premier cas de reconstruction d'image à résolution augmentée conformément au procédé de l'invention ; les FIGURES 6A à 6D montrent les résultats sur le second cas de reconstruction d'image à résolution augmentée conformément au procédé de l'invention ; la FIGURE 7 montre le PSNR concernant le second cas alors que le nombre d'itérations croît ; les FIGURES 8A à 8D montrent les résultats sur le troisième cas de reconstruction d'image à résolution augmentée conformément au procédé de l'invention. L'invention consiste à proposer un procédé de reconstruction d'une image à résolution augmentée avec des algorithmes perfectionnés de résolution augmentée, qui est basé sur la méthode itérative d'Irani et emploie une interpolation initiale, une sélection d'image automatique et un enregistrement d'image robuste que nous proposons. L'image à résolution augmentée, reconstruite par le procédé de l'invention, peut en outre être améliorée par un postprocédé d'amélioration d'image pour augmenter la qualité d'image de celle-ci. Dans ce qui suit, les étapes du procédé selon l'invention et des modes de réalisation préférés de l'invention seront décrits en détail pour illustrer l'esprit, les caractéristiques et les avantages de l'invention. Seront également illustrées dans ce qui suit la manière de reconstruire une image à résolution augmentée associée à une scène stationnaire et la manière de reconstruire une image à résolution augmentée associée à un objet en mouvement. Si l'on se réfère à la FIGURE 1, on peut voir que celle-ci est un organigramme montrant les phases de travail du procédé selon la présente invention. On soulignera que M premières images à faible résolution sont utilisées pour produire une image à résolution augmentée dans le procédé selon l'invention. De plus, les M premières images à faible résolution sont associées à la même scène, c'est-àdire que les M premières images à faible résolution sont toutes capturées à partir de la même scène. Chacune des M premières images est représentée par un ensemble de pixels uniformément échantillonnés et quantifiés. Une relation d'un facteur de grossissement désiré par l'utilisateur existe entre l'image à résolution augmentée et chacune des M premières images. Au départ, le procédé selon l'invention effectue l'étape S10 pour stocker les M premières images. Le procédé de l'invention effectue alors l'étape S12 pour sélectionner l'une des M premières images comme image prototype, et pour désigner les (M-1) premières images non sélectionnées comme étant les (M-l) secondes images. Par la suite, le procédé de l'invention effectue l'étape S14. La phase d'évaluation initiale du procédé de l'invention est réalisée pendant l'étape S14. L'étape S14 consiste à interpoler des pixels supplémentaires dans l'image prototype sur la base du facteur de grossissement. A l'étape S14, les valeurs des pixels supplémentaires sont calculées conformément aux valeurs de pixels voisins de ceux-ci pour produire une image prototype interpolée, et les translations respectives qui existent chacune entre chacune des (M-1) secondes images et l'image prototype interpolée sont également calculées. Après l'étape S14, l'étape S16 est effectuée. La phase de sélection d'image automatique du procédé de l'invention est réalisée pendant l'étape S16. L'étape S16 consiste à diviser les translations des (M-l) secondes images par le facteur de grossissement, respectivement, pour obtenir le module en ce qui concerne les translations des (M-1) secondes images. L'étape 16 consiste également à sélectionner, sur la base d'un critère, l'une des secondes images dont les modules apparentés sont les mêmes. Les secondes images sélectionnées conjointement avec le reste des secondes images, dont les modules apparentés ne sont pas les mêmes, sont désignées comme étant les N troisièmes images où N est égal ou inférieur à (M-l). Après cela, l'étape S18 est effectuée pour sous- échantillonner, conformément à la translation respective entre chacune des N troisièmes images et l'image prototype interpolée, l'image prototype interpolée N fois pour produire N quatrièmes images. Chacune des N quatrièmes images correspond à l'une des N troisièmes images. Après l'étape S18, l'étape S20 est effectuée pour calculer la différence entre chacune des N troisièmes images et la quatrième image correspondante de celle-ci. Ensuite, l'étape S22 est effectuée pour ajuster les valeurs des pixels de l'image prototype interpolée conformément à une moyenne des différences calculées à l'Etape 20. Finalement, l'étape S24 est effectuée pour déterminer si les valeurs des pixels de l'image prototype interpolée convergent vers un résultat satisfaisant. Si NON à l'étape S24, les étapes S18 à S 22 sont répétées. Dans le résultat, l'image prototype interpolée, dont les valeurs de pixels convergent vers le résultat satisfaisant, est désignée comme étant l'image à résolution augmentée. Dans un mode de réalisation préférée, l'étape S18 est mise en #uvre conformément à la formule suivante : gkn> Tk J () h J S , dans laquelle : gk(n) est la kème quatrième image au nème procédé de sous- échantillonnage ; f(n) est l'image prototype interpolée étant ajustée n fois ; h est l'opérateur de création de flou ( blurring ) défini par une fonction de diffusion de points ; Tk est l'opérateur de transformation en ce qui concerne la translation entre la kème troisième image et l'image prototype interpolée ; #s est l'opérateur de sous-échantillonnage. Dans un mode de réalisation préféré, l'Etape S18 et l'Etape S22 sont mises en #uvre conformément à la formule suivante : (n+1) - en) +, 1 'k (((gk - ) ) ' K k=1 dans laquelle : f(n+1) est l'image prototype interpolée pour être ajustée à (n+1) fois ; - K est le nombre des troisièmes images étant calculé ; - gk es la kème troisième image ; - p est l'opérateur de retrait de flou ( deblurring ) ; Tk-1 est l'opérateur de transformation inverse en ce qui concerne la translation entre la kème troisième image et l'image prototype interpolée ; - Ts est l'opérateur de sur-échantillonnage. On soulignera que la phase d'évaluation initiale affectera de façon significative la performance du système d'ordinateur étant donné que la méthode itérative est utilisée pour reconstruire une image à résolution augmentée. Autrement dit, si une meilleure évaluation initiale est appliquée, une grande quantité de temps de calcul sera économisée. En raison du fait que l'évaluation initiale est effectuée simplement une fois au début du procédé, la complexité de la méthode itérative globale ne dépend pas de la complexité de l'évaluation initiale, qui est basée sur les techniques d'interpolation. La différence entre ces modes d'interpolation de différents ordres est juste que les nombres de pixels voisins utilisés pour le calcul du procédé de calcul sur ordinateur sont différents. Par la suite, un mode d'interpolation du troisième ordre sera pris comme exemple pour illustrer les procédures des modes d'interpolation. Le mode d'interpolation du troisième ordre, ou appelé le mode d'interpolation cubique, considère 4 variables inconnues. Si l'on suppose que la fonction d'interpolation est : y = /3(x) = ax3 +bx2 +cx + d et que les pixels voisins connus comprennent (-l,A), (0, B) , (1, C) ,and (2,D), alors les formules suivantes peuvent être conduites : A -1 1 -1 1 a B 0 0 0 1 b C 1 1 1 le D D] J 8 4 2 1 d a -1 1 -1 IV1 N b 0 0 0 1 B c 1 1 1 1 C d 8 4 2 1 D - 0.1667 0.5 - 0.5 0.1667 A 0.5 -1 0.5 0 B - 0.3333 - 0.5 1 -0.1667 C 0 1 0 0 D Si l'on se réfère aux FIGURES 2A et 2B, les procédures du mode d'interpolation du troisième ordre sont encore illustrées. Comme représenté sur la Figure 2A, pour interpoler la valeur du pixel P dans une image à deuxdimensions, les valeurs de pixels A', B', C' et D' sont calculées au préalable. Les valeurs connues de pixels Al, B1, Cl et Dl sont utilisées pour déterminer les coefficients de la fonction du troisième ordre f3(x). Ensuite, la valeur de pixel A' est interpolée par l'utilisation d'une interpolation à une dimension. De façon analogue, les valeurs de pixels B', C' et D' sont interpolées conformément aux valeurs de pixels Ai, Bi, Ci et Di, alors que i=2,3 et 4. Finalement, la valeur de pixel P est calculée avec l'interpolation à une dimension dans la direction verticale comme démontré sur la FIGURE 2A. Comme représenté sur la FIGURE 2B, le mode d'interpolation du troisième ordre avec une interpolation à une dimension est utilisé pour déterminer la valeur de pixel P conformément aux valeurs connues de pixels A', B', C' et D' . De la même manière, en algorithme d'interpolation à deux dimensions, les coefficients de fn(x) concernant ces modes d'interpolation d'autres ordres peuvent également être résolus. Il a été observé que les modes d'interpolation de différents ordres conduisent à des images d'évaluation initiale différentes et à différentes vitesses de convergence de la qualité d'image alors que le nombre d'itérations croît. En sélectionnant l'ordre d'interpolation le plus approprié, les meilleurs résultats de la méthode itérative seront obtenus, étant donné que la phase d'évaluation initiale a une grande influence sur la performance de l'enregistrement d'image et sur le nombre nécessaire d'itérations pour parvenir au résultat de pic d'image. Dans une expérience, les performances au sujet de ces modes d'interpolation de différents ordres sont évaluées par le Rapport Signal de Pic à Bruit (PSNR) entre l'image initiale et les images reconstruites. Les résultats expérimentaux sont représentés sur la FIGURE 3. Comme représenté sur la FIGURE 3, dans la plupart des situations, le mode d'interpolation du troisième ordre est au rang du meilleur choix d'évaluation initiale si à la fois la complexité et la qualité d'image reconstruite sont concernées. Le procédé selon la présente invention préfère utiliser seulement l'interpolation du premier ordre, l'interpolation du troisième ordre ou l'interpolation du cinquième ordre. Dans un mode de réalisation préféré, alors que les (M-1) secondes images sont les séquences d'images à faible résolution associées à un objet en mouvement, la translation (LTi) à l'étape S14, entre la ième seconde image et l'image prototype interpolée est calculée conformément aux étapes suivantes consistant à : (a) assigner au moins un point intéressant Pi(x,y) sur la ième seconde image, chaque point intéressant correspondant à une fenêtre locale (w) et un point (u,v) sur l'image prototype interpolée ; (b) calculer une différence absolue (LADi(x,y;u,v)) entre chaque point intéressant P1(x,y) et le point correspondant (u,v) de celui-ci, conformément à la formule suivante : LAD,(x,y;u,v)= \l,(x + m,y + ri)- I0(u + m,v+ ri)\ (m,n)Ew (c) calculer un minimum (LRi(x,y) de la différence absolue LADi(x,y;u,v) calculée à l'étape (b) conformément à la formule suivante : LR= (x, y) = arg min LADt (x, y; u, v) ; et (u,v) (d) calculer la translation (LT1), en ce qui concerne chaque point intéressant Pi(x,y) conformément à la formule suivante : LT1(x,y)= LRi(x,y)-(x,y)#facteur de grossissement. Dans ce mode de réalisation préféré, pour parvenir à un enregistrement d'image plus précis, puis reconstruire l'image à résolution élevée d'un objet en mouvement, les points intéressants sont choisis selon les contraintes suivantes : a. le gradient au point intéressant doit être plus grand qu'un seuil ; b. pour chaque point intéressant sur une image à faible résolution, la complexité correspondant au point correspondant sur l'image à résolution augmentée simulée doit être supérieure à la complexité des points périphériques ; c. la translation entre chaque point intéressant et le point correspondant ne doit pas être zéro. La technologie mentionnée ci-dessus est également appelée la technologie d'appariement locale comme représenté sur la FIGURE 4A. Dans un autre mode de réalisation préféré, alors que les (M-l) secondes images sont les séquences d'images à faible résolution associées à un objet stationnaire, la translation (LTi), à l'étape S14, entre la ième seconde image et l'image prototype interpolée est calculée conformément aux étapes suivantes consistant à : (b) calculer une différence absolue (GAD(u,v)) en ce qui concerne la ième seconde image conformément à la formule : GAD (u, v) ~ 7,()-(M+v+ (,)6; dans laquelle le point (u,v) est un point de position correspondante de l'image prototype interpolée qui correspond à un point de départ de la ième seconde image ; (b) calculer un minimum (GR(i)) de la différence absolue (GAD (u,v)) calculée à l'étape (a) conformément à la formule suivante : GR(i) = arg min GAD, (u, v) (u,v) et considérer GR (i) comme LTi. La technologie mentionnée ci-dessus est également appelée technologie d'appariement globale, comme représenté sur la FIGURE 4B. A l'étape S16, l'objectif de la sélection d'image automatique est de sélectionner un nombre minimal d'images à faible résolution utiles pour les procédés de traitement d'image tel que l'enregistrement d'image. Ce mode peut réduire un calcul non nécessaire pour accélérer la reconstruction de l'image à résolution augmentée. Les principes de sélection d'images à faible résolution utiles à partir des images à faible résolution initiales, associées à une scène stationnaire ou à un objet en mouvement, sont différents. Les différents principes seront illustrés comme suit. Dans un mode de réalisation, alors que les (M-l) secondes images sont les séquences d'images à faible résolution associées à un objet en mouvement, le principe de sélectionner des images à faible résolution utiles à l'étape S16 est que, pour deux secondes images i, j ayant le même module et (ui,vi)=Ti, (uj, vj)=Tj, la seconde image i est sélectionnée si GADi(u1,v1) Dans un mode de réalisation, alors que les (M-l) secondes images sont les séquences d'images à faible résolution associées à une scène stationnaire, le principe de sélection d'image à faible résolution à l'étape S16 est que, pour deux secondes images j et j ayant le même module, et (ui,vi)=Ti et (uj,vj)=Tj, la seconde image i est sélectionnée si 62 Si l'on se réfère au Tableau 1, on peut voir que le Tableau 1 indique les performances d'un système appliquant le procédé selon l'invention. Ce système est équipé d'une unité centrale CPU Pentium Intel III et de mémoires vives RAM de 128MB. Dans ce cas, cinq images à faible résolution 62 x 62 sont utilisées pour reconstruire une image à résolution augmentée sur la base d'un facteur de grossissement 3. Dans le Tableau 1, les résultats de la sélection du système exécutant une sélection d'image automatique ou non sont également énumérés. Tableau 1 Temps de Calcul PSNR (db) (secondes) ~~~~~~~ Technique Avec Sélection 155,4 26,78 d'Appariement d'Image Locale Sans Sélection 582,4 26,66 Locale d'Image 582,4 26,66 d'Image , Avec Sélection ~c -, -# Technique Avec Sélection 75,8 26,78 d'Appariement d'Image Globale Sans Sélection 496,2 26,66 d'Image ' Conformément aux résultats du Tableau 1, il apparaît à l'évidence que le système exécutant la sélection d'images automatique proposé par l'invention peut réduire le temps de calcul de façon significative. Certains cas de reconstruction d'image utilisant le procédé de l'invention seront décrit dans ce qui suit pour améliorer encore la performance du système appliquant le procédé de l'invention. Si l'on se réfère aux FIGURES 5A à 5D, un premier cas est illustré. La FIGURE 5A montre l'une des images à faible résolution. La FIGURE 5B montre le résultat après la phase d'évaluation initiale. La FIGURE 5C montre l'image reconstruite après 100 itérations. La FIGURE 5D montre l'image à résolution augmentée reconstruite avec un procédé de filtrage passe-haut pour rendre nettes les bordures de celle-ci. Si l'on se réfère aux FIGURES 6A à 6D, un second cas est illustré. Dans le second cas, une image initiale est prise comme scène stationnaire et est sous- échantillonnée par une caméra simulée pour produire une pluralité de séquences d'images à faible résolution. A l'aide de la caméra simulée, ces points sur l'image initiale, auxquels les points de commencement des images à faible résolution correspondent, sont différents, c'est-àdire que la caméra simulée se déplace lors de la capture des images à faible résolution. Ensuite, les images à faible résolution sont utilisées pour reconstruire une image à résolution augmentée avec un facteur de grossissement de 4, conformément au procédé de l'invention. La FIGURE 6A montre l'une des images à faible résolution dans le second cas. La FIGURE 6B montre le résultat après la phase d'évaluation initiale. La FIGURE 6C montre l'image reconstruite après 100 itérations. La FIGURE 6D montre l'image améliorée reconstruite avec un procédé postimage. Comme représenté sur la FIGURE 7, le PSNR entre l'image initiale et l'image reconstruite dans le second cas converge rapidement alors que le nombre d'itérations croît. Par le résultat de la FIGURE 7, il est manifeste que le système appliquant le procédé de l'invention fournit une bonne performance. Si l'on se réfère aux FIGURES 8A à 8D, un troisième cas est illustré. Dans le troisième cas, plusieurs images à faible résolution, associées à un objet en mouvement par rapport à un arrière-plan stationnaire, sont utilisées pour reconstruire une image à résolution augmentée sur la base d'un facteur de grossissement de 2. La FIGURE 8A montre l'une des séquences d'images à faible résolution. La FIGURE 8B montre le résultat après la phase d'évaluation initiale. La FIGURE 8C montre l'image reconstruite après 100 itérations. La FIGURE 8D montre l'image améliorée reconstruite avec un procédé de postimage. De la FIGURE 8A à la FIGURE 8D, comme cela apparaît, étant donné que le nombre d'itérations augmente, l'object devient plus distinct alors que l'arrière-plan devient flou et les bordures de mots sont davantage discernables. Pour résumer les modes de réalisation ci-dessus et les cas pratiques, l'invention propose un procédé qui reconstruit des images à résolution augmentée avec des algorithmes perfectionnés de résolution augmentée, et emploie l'interpolation initiale, la sélection d'image automatique et l'enregistrement d'image robuste que nous proposons. Conformément au procédé de l'invention, la reconstruction d'une image à résolution augmentée peut être rapidement obtenue non seulement par des séquences d'images à faible résolution associées à une scène stationnaire mais encore par des séquences d'images à faible résolution associées à un objet en mouvement. Les images à résolution augmentée reconstruites peuvent encore effectuer le postprocédé pour augmenter la qualité d'image. Il va de soi que le procédé de l'invention satisfait aux conditions de nouveauté et d'activité inventive et satisfait les directives sur les demandes de brevets. Avec l'exemple et les explications ci-dessus, les caractéristiques et l'esprit de l'invention seront bien décrits. L'homme du métier observera facilement que de nombreuses modifications et variantes du dispositif qui peuvent être effectuées tout en conservant l'enseignement de l'invention. En conséquence, la description ci-dessus doit être considérée comme limitée seulement par la portée des revendications annexées. REVENDICATIONS 1 - Procédé d'obtention d'une image à résolution augmentée, caractérisé par le fait qu'il comprend l'utilisation de M premières images à faible résolution associées à une scène, chacune des M premières images étant représentée par un ensemble de pixels uniformément échantillonnés et quantifiés, une relation d'un facteur de grossissement existant entre l'image à résolution augmentée et les M premières images, ledit procédé comprenant les étapes consistant à : (a) stocker les M premières images (S10) ; (b) sélectionner l'une parmi les M premières images en tant qu'image prototype (S12), et désigner les (M-1) premières images non-sélectionnées comme étant les (M- 1) secondes images ; (c) sur la base du facteur de grossissement, interpoler des pixels supplémentaires dans l'ensemble de pixels de l'image prototype, puis déduire les valeurs des pixels supplémentaires conformément aux valeurs de pixels voisins de ceux-ci pour produire une image prototype interpolée, et calculer une translation respective entre chacune des (M-1) secondes images et l'image prototype interpolée (S14) ; (d) diviser les translations des (M-1) secondes images par le facteur de grossissement, respectivement, afin d'obtenir le module en ce qui concerne les translations des (M-1) secondes images, sélectionner, sur la base d'un critère, l'une parmi les secondes images dont les modules apparentés sont les mêmes (S16), et désigner les secondes images sélectionnées conjointement avec le reste des secondes images dont les modules apparentés ne sont pas les mêmes comme étant les N troisièmes images, où N est égal à ou inférieur à (M-l) ; (e) sous-échantillonner l'image prototype interpolée N fois conformément à la translation respective entre chacune des N troisièmes images et l'image prototype interpolée pour produire N quatrièmes images (S18) qui correspondent chacune à l'une des N troisièmes images ; (f) calculer la différence entre chacune des N troisièmes images et la quatrième image correspondante de celle- ci (S20) ; (g) ajuster les valeurs des pixels de l'image prototype interpolée conformément à une moyenne des différences calculées à l'étape (f) (S22) ; et (h) répéter l'étape (e) à l'étape (g) jusqu'à ce que les valeurs des pixels de l'image prototype interpolée convergent vers un résultat satisfaisant (S24), et désigner l'image prototype interpolée dont les valeurs de pixels convergent vers le résultat satisfaisant comme étant l'image à résolution augmentée. 2 - Procédé selon la revendication 1, caractérisé par le fait que l'étape (e) est effectuée conformément à la formule suivante : gkn) Tk \J (n) h J S , dans laquelle : gk(n) est la kème quatrième image au nème procédé de sous- échantillonnage ; en) est l'image prototype interpolée étant ajustée n fois ; h est l'opérateur de création de flou défini par une fonction de diffusion de points ; Tk est l'opérateur de transformation en ce qui concerne la translation entre la kème troisième image et l'image prototype interpolée ; s est l'opérateur de sous-échantillonnage. 3 - Procédé selon la revendication 2, caractérisé par le fait que les étapes (f) et (g) sont effectuées conformément à la formule suivante : (n+1) - en) + fTk-1 (((gk - gk(n) S p), K k=1 dans laquelle : f(n+1) est l'image prototype interpolée pour être ajustée à (n+l)fois ; - K est le nombre des troisièmes images étant calculées ; - gk est la kçme troisième image ; - p est l'opérateur de retrait du flou ; Tk-1 est l'opérateur de transformation inverse en ce qui concerne la translation entre la Kème troisième image et l'image prototype interpolée ; - #s est l'opérateur de sur-échantillonnage. 4 - Procédé selon la revendication 3, caractérisé par le fait qu'à l'étape (c), la déduction des valeurs des pixels supplémentaires est basée sur l'une choisie dans le groupe consistant en une interpolation de premier ordre, une interpolation de troisième ordre et une interpolation de cinquième ordre. 5 - Procédé selon la revendication 4, caractérisé par le fait qu'à l'étape (c), la translation (LTi) entre la ième seconde image et l'image prototype interpolée est calculée conformément aux étapes suivantes : (c-1) assigner au moins un point intéressant Pi(x,y) sur la ième seconde image, chaque point intéressant correspondant à une fenêtre locale (w) et à un point (u,v) sur l'image prototype interpolée ; (c-2) calculer une différence absolue (LADi(x,y;u,v)) entre chaque point intéressant Pi(x,y) et le point correspondant(u,v) de celui-ci conformément à la formule : LADi(x,y;u,v)= LII(x+ m,y+n)-Io(u+ m,v+ n)1 (m,n)#w (c-3) calculer un minimum (LRi(x,y)) de la différence absolue (LADi(x,y;u,v) calculée à l'étape (c-2) conformément à la formule suivante : LR(x,y) = arg min LAD (x, y; u, v) (u,v) (c-4) calculer la translation (LT1) en ce qui concerne chaque point intéressant Pi(x,y) conformément à la formule suivante : LTi(x,y)= LRi(x,y)-(x,y)#facteur de grossissement. 6 - Procédé selon la revendication 4, caractérisé par le fait qu'à l'étape (c), la translation (LTi) entre la ième seconde image et l'image prototype interpolée est calculée conformément aux étapes suivantes : (c-1) calculer une différence absolue (GAD(u,v)) en ce qui concerne la ième seconde image entière conformément à la formule suivante : G.()= /,()-(M++ (x,y)E dans laquelle le point (u,v) est un point, sur l'image prototype interpolée, correspondant à un point de départ de la ième seconde image ; (c-2) calculer un minimum (GR(i)) de la différence absolue (GAD(u,v)) calculée à l'étape (c-1) conformément à la formule suivante : GR(i) = arg min GAD1 (u, v) (u,v) et considérer GR(i) comme LT1.