FR 2754362 A1 19980410 FR 9712418 A 19971006 PROCEDE DE VISUALISATION D'UNE IMAGE D'UNE SOURCE D'ONDES La présente invention concerne un procédé de visualisation d'une image d'une source d'ondes qui est basé sur des données de me- sure d'hologramme obtenues par observation d'un hologramme d'ondes radioélectriques. Des procédés de visualisation d'ondes radioélectriques qui utili- sent l'interférence d'ondes comprennent à la fois des techniques de vi- sualisation actives comme en radar à synthèse d'ouverture, en radar à neige et en radar souterrain, et des techniques de visualisation passives, comme la radioastronomie et la radiométrie [(1) Masaru Matsuo, Ku- niyoshi Yamane: Radar Holography, publié par Electronic Communication Society, 1980. (2) Yoshinao Aoki: Wave Signal Processing, Morikita Shuppan, 1986.]. De tels procédés ont déjà été proposés pour des appli- cations dans la détection de sources d'ondes produisant des émissions d'ondes radioélectriques indésirables [(3) Jun'ichi Kikuchi, Motoyuki Sato, Yoji Nagasawa, Risaburo Sato; A proposai for Searching the Electroma- gnetic Wave Sources by Using a Synthetic Aperture Technique, IECE transactions on B, Vol. J68-B, n 10, pages 1194-1201, octobre 1985.]. Des algorithmes à haute résolution spectrale faisant appel aux procédés dits MLM ("Maximum Likelihood Method", c'est-à-dire procédé à vraisem- blance maximale), MEM ("Maximum Entropy Method", c'est-à-dire procédé à entropie maximale) et MUSIC ("Multiple Signal Classification", c'est-à- dire classification de signaux multiples) 1(4) S. Kesler: Modem Spectrum Analysis II, IEEE Press, 1986)] ou ESPRIT amélioré ("Estimation of Si- gnal Parameters via Rotational Invariance Techniques", c'est-à-dire esti- mation de paramètres de signal par des techniques d'invariance en rota- tion) [(5) R. Roy et T. Kailath: ESPRIT-Estimation of Signal Parameters via Rotational Invariance Techniques, IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-37, n 7, pages 984-995, juillet 1989.] ont également été proposés pour obtenir des résolutions plus élevées dans la détection d'une source d'ondes éloignée. De plus, des exemples de l'application d'algorithmes MEM à la détection d'une source d'ondes proche ont également été signalés [(2) Yoshinao Aoki: Wave Signal Processing, Morikita Shuppan, 1986. (6) Jun'ichi Kikuchi, Yoji Nagasawa: Search for Electromagnetic Wave Sources by Using the Maximum Entropy Method, IECE Transactions on B, Vol. J69-B, n 9, pa- ges 949-957, septembre 1986.]. Cependant, les algorithmes MLM, MUSIC et ESPRIT sont limités à l'estimation de sources ponctuelles, et en MEM on peut sélectionner soit un mode de reconstruction de source d'ondes ponctuelle, soit un mode de reconstruction de source d'ondes ayant un étalement corres- pondant à une distribution normale, conformément au nombre de termes du filtre d'erreur de prédiction. Cependant, lorsqu'on a affaire à une émission indésirable pro- venant d'un dispositif réel, la source d'ondes présente rarement une dis- tribution correspondant à un nombre limité de sources quasi-ponctuelles, et elle est plus probablement une source d'ondes plane étendue. En d'autres termes, on considère que les algorithmes mentionnés ci-dessus sont inappropriés pour visualiser la distribution de sources d'ondes d'une configuration qui n'est pas nette. La visualisation qui a pour objet une émission indésirable doit donc considérer la distribution de sources d'on- des comme une distribution de sources d'ondes continue, plutôt que comme une combinaison d'un nombre limité de sources d'ondes ponc- tuelles. Les inventeurs de la présente invention ont également proposé un algorithme SPIM ("Spectrum Phase Interpolation Method", c'est-à-dire procédé d'interpolation de phase de spectre) et son application pour la reconstruction d'une image de sources d'ondes correcte d'une configura- tion dans laquelle des sources d'ondes ponctuelles coexistent avec une source d'ondes qui s'étale sur une certaine étendue [(7) Hitoshi Kitayoshi: High Resolution Technique for Short Time Frequency Spec- trum Analysis, IEICE Transactions on A, Vol. J76-A, n 1, pages 78-81, janvier 1993; (8) Hitoshi Kitayoshi: High Resolution Techniques for 2-D Complex Spectrum Analysis, IEICE Transactions on A, Vol. J76-A, n 4, pages 687-689, avril 1993]. Néanmoins, ce procédé exige également la fixation d'une valeur de seuil pour estimer une source d'ondes ponctuelle, et il est donc inap- proprié pour la visualisation d'une émission indésirable. On a signalé un exemple dans lequel, à titre de technique de visualisation pour travailler avec une distribution de sources d'ondes continue, la relation entre le champ observé et une source d'ondes d'une forme arbitraire est représentée sous la forme d'une intégrale de diffrac- tion de Fresnel-Kirchhoff [(9) Yoshinao Aoki, Shigeki Ishizuka: Numerical Two-Dimensional Fresnel Transform Methods, IECE Transactions on B, Vol. J57-B, n 8, pages 511-518, août 1974]. Cependant, ce procédé im- pose une limite à la région qui peut être visualisée, du fait qu'il prend pour son objet une région de Fresnel. Ce procédé souffre en outre du défaut consistant en une mauvaise résolution aux bords de la plage vi- sualisée, du fait que l'on utilise une quasi-transformation de Fourier et que les résultats de reconstruction deviennent un espace de nombres d'ondes. Enfin, des recherches concernant des mesures de champ qui prennent en considération la directivité de l'antenne comprennent une mesure de champ au voisinage de l'antenne [(10) A.D. Yaghjian: An Overview of Near-Field Antenna Measurements, IEEE Transactions on Antennas and Propagation, Vol. AP-34, n 1, pages 30-45, janvier 1986; et (11) Tasaku Teshirogi, Near-Field Antenna Measurement, IECE Maga- zine, Vol. 62, n 10, pages 1145-1153, octobre 1979.]. Cependant, dans ces cas, la distribution de champ est mesurée en utilisant une surface de mesure (D x D) suffisamment grande par rapport à l'antenne (L x L), qui est l'objet d'observation, dans des conditions dans lesquelles L dans la mesure est un dispositif qui est destiné à l'étude du champ loin- tain, et par conséquent le procédé ne peut pas être utilisé sans modifica- tion dans des conditions dans lesquelles L > D. On expliquera ensuite un procédé de reconstruction numérique d'hologramme radioélectrique de l'art antérieur, qui n'utilise pas l'ap- proximation de Fresnel ou l'approximation de Fraunhofer. Une source de courant de surface J(R') existe à un point de source d'ondes x', y', z' dans les coordonnées orthogonales qui sont re- présentées dans le modèle d'observation d'hologramme de la figure 1, et le champ au point d'observation x, y; z = 0 qui est ainsi créé est E(R). De plus, la plage de distribution de la source d'ondes J(R') est un plan bidi- mensionnel fini à z = z', avecÀ -L/2 L/2 L.L/2 E(R) = -Jỏ)t0 " j G (R - R'). J (R') dx'dy' (1) J-L12 J-L/2 Cependant, si la longueur de vecteur effective d'une antenne qui est utilisée pour observer le vecteur de champ E est fixée à Il et si la longueur d'onde est X, la tension V pour la réception dans la région r= R - R' " > est' V = g Il. E (2) Ici, g est une constante. En d'autres termes, on peut considérer Il comme une fonction indiquant la directivité et la sensibilité qui ne dé- pend pas de la distance r. D'après les équations (1) et (2), la tension de sortie d'une antenne qui est placée au point d'observation R peut être donnée par: L/2 1/2 V(R) = - jolotg I11(R-R'). G(R-R').J(R')dx'dy' (3) J-L/2 -L/2 L'équation (3) donne la tension de réception d'une antenne quelconque en relation avec n'importe quelle distribution de courant de surface. Cependant, si on considère le système de balayage d'antenne dans une mesure d'hologramme réelle, on peut aisément réaliser une antenne qui reçoit le champ dans les directions x et y, mais une réception précise du champ dans la direction z est extrêmement difficile. Malgré la présence d'un vecteur de courant de la composante z sur la surface du dispositif qui est l'objet de la mesure, presque rien n'est reçu dans le plan de mesure d'hologramme. L'explication qui est présentée ci-dessous est limitée à un cas dans lequel le courant du plan a seulement une composante horizontale et une composante verticale, et dans lequel la réception est effectuée par des antennes de polarisation horizontale et de polarisation verticale. Des antennes de polarisation horizontale et de polarisation ver- ticale sont placées dans le plan de mesure d'hologramme, et les vecteurs de tension qui sont reçus par ces antennes sont LVh (R)1 [h1 V(R) = V (R) = v. E(R) (4) Ici, Vh et Vv sont les tensions reçues des antennes de polarisa- tion horizontale et de polarisation verticale, et Ih et Ilv sont les Ion- gueurs de vecteur effectives de ces antennes. Le dispositif utilise deux antennes de réception ayant la même forme, la sensibilité de polarisation principale est A0, et la sensibilité de polarisation croisée est A+. De plus, I'angle zénithal et l'angle d'azimut du point R' sont Oh et h, lorsqu'on ob- serve ce point à partir de l'antenne de réception de polarisation horizon- tale qui est placée au point R, et l'angle zénithal et l'angle d'azimut du point R' sont Ov et 'v lorsqu'on observe ce point à partir de l'antenne de réception de polarisation verticale. Les angles d'azimut *h et *v repré- sentent ici des angles qui sont respectivement mesurés par rapport à la direction de l'axe x et de l'axe y. Si la partie incluse dans le membre de droite de l'équation (3) qui représente la sensibilité de réception due à la source de courant ponctuelle, est représentée par: Ih(R -R) B(R-R')=-j. ogliG(R - R') [A@(h,h)sin6h A,0(0h 'hi)Sin1v exp(-j2Xr/l) (5) [A#(eO)si$nh AO(ev,,, )sinJ r exp(-j2xr/a) =A(R-R') r On peut trouver Vh et Vv par l'équation suivante, basée sur les équations (3) et (4): L12 1.L/2 V(R) = B(R - R'). J(R') dx'dy' (6) -L/2 2-L2 D'autre part, les équations (5) et (6) donnent J(R') = A1-(R-R').V(R). exp(j27rr/X).rdxdy (7) Conformément à l'équation (4), les vecteurs de tension reçue sont: V(R) = [Vh (R) = Vh( y) (8)) LVv (R)j LVv(,, y]) Et d'après l'équation (7), seulement deux composantes de vec- teurs de courant de source d'ondes sont visualisées pour ces vecteurs de tension: ix (R)] rJx(x', Y" Z)l J(R') = (R') LJy(x, y', z)(9) L'art antérieur décrit ci-dessus présent les problèmes suivants D'après l'équation (6), V(R) comprend N x N éléments de don- nées et le point visualisé de J(R') comprend M x M éléments, ce qui exige un calcul de l'ordre de N2 x M2 et exige en outre un temps considérable pour obtenir les résultats de visualisation. On ne peut pas mesurer de façon précise le vecteur de courant de source d'ondes du fait que, d'après l'équation (9), on ne prend pas en considération la directivité correspond à la composante Jz d'une antenne- sonde qui est utilisée dans la mesure de V(R). De plus, on ne peut visua- liser que deux composantes des vecteurs de courant de source d'ondes. Le but de la présente invention est de procurer un procédé de visualisation d'une image d'une source d'ondes qui permette de réaliser une holographie par le calcul plus rapide, en utilisant une transformation de Fourier rapide (ou FFT pour "Fast Fourier Transform") avec division partielle, et qui soit capable de réduire à un niveau pratique le temps né- cessaire pour obtenir des résultats de visualisation. Un autre but de la présente invention est de procurer un procé- dé de visualisation d'une image d'une source d'ondes qui permette la vi- sualisation de trois composantes de vecteurs de courant de source d'on- des, à partir de données de mesure d'hologramme. Le procédé de visualisation d'une image d'une source d'ondes de la présente invention est mis en oeuvre par holographie par le calcul, qui consiste en un procédé de reconstruction numérique de données de mesure d'hologramme qui sont obtenues par l'observation d'un holo- gramme radioélectrique; et il comprend: une étape de segmentation de la région de reconstruction d'une image reconstruite d'hologramme, par la taille de la surface de mesure d'hologramme; une étape d'exécution d'une première transformation de Fourier rapide (FFT) après multiplication des données de mesure d'hologramme par une fonction de pondération; une étape d'exécution d'une seconde FFT après avoir donné la distance focale de reconstruction et une paire de nombres désignant un segment d'une image reconstruite, et de multiplication par une fonction de pondération d'une valeur qui est obtenue en multipliant une matrice inverse de sensibilité d'antenne pour chaque segment, par une fonction de propagation inverse; une étape d'exécution d'une transformée de Fourier rapide in- verse (ou IFFT pour "Inverse Fast Fourier Transform"), après la multipli- cation de l'information de sortie de la première FFT par la seconde FFT; une étape de compensation, au niveau de chaque segment, de la fonction de pondération qui est utilisée dans la multiplication des don- nées de mesure d'hologramme; une étape de compensation, au niveau de la région de recons- truction, de la fonction de pondération qui est utilisée dans la multiplica- tion de la valeur obtenue en multipliant la matrice inverse de sensibilité d'antenne, pour chaque segment, par la fonction de propagation inverse; une étape consistant à visualiser l'information de sortie de cette étape à une position qui désigne le segment, et à faire de l'information de sortie de cette étape une image reconstruite au niveau du segment; et l'accomplissement de toutes ces opérations pour tous les seg- ments, et la formation d'une image reconstruite au niveau de l'holo- gramme. L'utilisation de la FFT avec division partielle réduit ainsi le temps nécessaire pour obtenir les résultats de visualisation, à partir des données de mesure d'hologramme. De plus, la segmentation peut être effectuée en chevauche- ment. En d'autres termes, par la segmentation de l'étendue de re- construction d'une image reconstruite au niveau de l'hologramme, dans laquelle des segments se chevauchent, on peut éviter l'utilisation de données situées à des bords, ou une erreur de coupure est grande, et une image reconstruite peut être générée sous la forme de segments dans laquelle la précision s'améliore avec la proximité vis-à-vis du centre de chaque segment. Le procédé de visualisation d'une image d'une source d'ondes de la présente invention permet la visualisation de composantes de vec- teurs de courant de sources d'ondes à partir de données de mesure d'hologramme qui sont obtenues par l'observation d'un hologramme ra- dioélectrique, et il comprend: une étape d'accomplissement d'une observation d'hologramme radioélectrique sur deux surfaces de balayage, à savoir une première surface de balayage et une seconde surface de balayage, de mise en place d'au moins deux antennes-sondes, à savoir une antenne à polari- sation horizontale et une antenne à polarisation verticale, sur chaque surface de balayage, et d'enregistrement de la distribution de vecteurs de tension qui est reçue par les antennes; une étape de calcul des caractéristiques directionnelles totales des antennes-sondes, en utilisant un procédé de moment, et de détermi- nation de la matrice de sensibilité d'antenne de réception pour l'antenne de réception à polarisation horizontale, et de la matrice de sensibilité d'antenne de réception pour l'antenne de réception à polarisation verti- cale, en relation avec des vecteurs de courant de source d'ondes; et une étape de comparaison des déterminants des matrices de sensibilité d'antenne de réception, de détermination de la grandeur de l'erreur de la matrice de sensibilité d'antenne de réception pour l'antenne de réception à polarisation verticale, et de l'erreur de la matrice de sen- sibilité d'antenne de réception pour l'antenne de réception à polarisation horizontale, et d'utilisation sélective des vecteurs de tension reçus et des matrices inverses de sensibilité d'antenne de réception, pour trouver et visualiser trois composantes de vecteurs de courant de source d'ondes. Les trois composantes de vecteurs de courant de source d'on- des peuvent donc être visualisées à partir de données de mesure d'holo- gramme, en accomplissant chacune des étapes décrites ci-dessus. L'invention sera mieux comprise à la lecture de la description qui va suivre de modes de réalisation, donnés à titre d'exemples non li- mitatifs. La suite de la description se réfère aux dessins annexés, dans lesquels: La figure I montre un modèle d'observation d'hologramme; La figure 2 est un organigramme montrant le premier mode de réalisation de la présente invention; La figure 3 montre la structure et la distribution de courant d'un élément d'antenne de réception à balayage, pour l'utilisation dans la me- sure d'hologramme, avec une analyse utilisant un procédé de moment, et elle consiste en une photographie d'une image tramée qui est représen- tée sur un dispositif de visualisation; La figure 4 est une représentation graphique montrant une comparaison de la résolution d'une image reconstruite en fonction de la position d'une source d'ondes sur l'axe x; La figure 5 est une représentation graphique montrant la dé- pendance de la valeur d'amplitude d'une image reconstruite vis-à-vis de la position de la source d'ondes; La figure 6 est une représentation graphique montrant la dé- pendance de la résolution de reconstruction d'hologramme vis-à-vis de la distance z' depuis la surface de mesure jusqu'au point de source d'on- des; La figure 7 est une représentation de structure schématique montrant un dispositif de mesure préparé sur une base expérimentale pour observer un hologramme radioélectrique; La figure 8 est une représentation graphique montrant les ré- sultats expérimentaux de valeurs d'amplitude obtenues dans la recons- truction d'image d'hologramme d'une antenne dipôle placée sur un axe horizontal; La figure 9 est un tableau montrant la résolution de visualisa- tion de sources d'ondes ponctuelles; La figure 10 est une photographie d'une image tramée qui est représentée sur un dispositif de visualisation d'images d'ondes électro- magnétiques d'un équipement sous test, observées à 10 GHz; La figure 11 est une photographie d'un équipement sous test qui est utilisé dans des expériences de visualisation d'images d'ondes électromagnétiques; et La figure 12 est une représentation schématique de structure d'un dispositif de mesure préparé sur une base expérimentale pour ob- server des hologrammes radioélectriques, correspondant au second mode de réalisation de la présente invention. Le premier mode de réalisation de la présente invention utilise un procédé d'élaboration du procédé de reconstruction numérique d'holo- gramme de la présente invention, qui n'utilise pas une approximation de Fresnel ou une approximation de Fraunhofer, à partir d'une transforma- tion de Fourier faisant appel à une division partielle, et il comprend une compensation de la directivité d'antennes de réception. Premièrement, dans l'équation (6), la plage d'observation E(R) est -D/2 L = MD (10) la modification qui est indiquée dans l'équation suivante est possible M-1 M-1 D/2 D/2 V(R) = 2 JD/2 (RR'mn).J(R'mn)dsdt (11) m=O n=O On a ici R'mn = (S + xm, t + Yn, z'), et xm et yn sont donnés par: Xm = (m - (M - 1)/2)D (12) Yn = (n - (M - 1)/2)D (13) Dans ce cas, Bm, (s, t)= A(R'mn) exp(- j2r'mn/) (14) r'm,n avec r'm,n = I R'm,n I. Comme on peut le comprendre d'après l'équation (5) et l'équa- tion (14), I'équation (11) est une intégrale de convolution de J(R'm,n) et Bmn(s, t), et on modifie donc l'équation (11) en représentant la transfor- mée de Fourier bidimensionnelle (s, t) par F{.}, et la transformée de Fou- rier inverse par F-1{.}. M-1 M-1 V(R) = Y 1 F'1{F{Bmn(s, t)}.F{J(R'm,n)}} (15) m=O n=O De plus, en utilisant la concordance de la taille de la surface d'observation de V(R) et de la taille de la surface de (s, t), la transforma- tion de Fourier et la multiplication par F{B-1a,b (s, t)} des deux membres de l'équation (15) donnent: F{B-1a,b(S, t)}.F{V(R)} = M-1 M-1 F{B-1a,b(s, t)}. ( 1 Y F{Bmn(S, t)}.F{J(R'mn)}) (16) m=O n-:O Ici, Ba,b(S, t) = '(R'a,b)r'ab exp(J2tr'ab/X) (17) et A-1(R'a,b) est la matrice inverse de A(R'a,b). Sur la base des expressions (14) et (17) qui figurent dans le membre de droite de l'équation (16), on note que lorsque le changement d'amplitude de A(R'mn) par rapport à s et t est suffisamment progressif, et en outre, lorsque le changement de phase est suffisamment faible en comparaison avec 2ir'm,n/I, c'est-à-dire lorsque la directivité d'amplitude et de phase d'une antenne de réception ne sont pas complexes, F{B-1a,b(s, t)}.F{Bmn(S, t)} devient une matrice unité seulement lorsque on a à la fois m = a et n = b, et dans le cas contraire elle est pratique- ment égale à zéro. Par conséquent, l'équation (16) devient F{B-la,b(S, t)}.F{V(R)} = F{J(R'a b)} (18) Une transformation de Fourier inverse du membre de droite donne: J(R'a b) = F-1{F{B-1a,b(S, t)}. F{V(R)}} (19) On trouve ensuite la distribution de sources d'ondes J(R') sur la base des données de mesure d'hologramme V(R). On a ici J(R') = Jx (R') (20) [Jy (R')] De cette manière, parmi les surfaces d'objet qui sont divisées par M x M, le numéro a, b de la surface qui doit être visualisée est subs- titué dans l'équation (19), grâce à on peut estimer la source d'ondes de cette surface. De plus, on peut réaliser un calcul rapide en utilisant une FFT bidimensionnelle pour la transformation de Fourier et la transforma- tion de Fourier inverse qui sont incluses dans l'équation (19). Si le nombre d'éléments de données pour V(R) est N x N et si le nombre d'éléments de données de J(R') à visualiser est M x M, I'utilisa- tion de l'équation (19) réduit le volume de calcul à une valeur de l'ordre de 20(M/N)2 x N 10og2 (2N), et le temps nécessaire pour obtenir des ré- sultats de visualisation est également réduit. Par exemple, dans un cas dans lequel N = 64 et M = 640, le volume de calcul peut être réduit ap- proximativement au millième de celui du procédé de l'art antérieur. La figure 2 est un organigramme montrant le premier mode de réalisation de la présente invention pour accomplir le calcul de l'équation (19). On expliquera ensuite une modification du premier mode de réalisa- tion de la présente invention, en se référant à la figure 2. Premièrement, on segmente la région de reconstruction d'une image d'hologramme reconstruite, par la taille (D x D) d'une surface de mesure d'hologramme. Ensuite, on multiplie les données de mesure d'hologramme V(R) par une fonction de pondération (par exemple cos(ulRI/21/2D))112.(S1), après quoi on accomplit une transformation de Fourier rapide (FFT) (S2). Ensuite, on affecte des valeurs à a et b, qui désignent un seg- ment de l'image reconstruite, et à une distance focale de reconstruction z (S3), et on multiplie la matrice inverse de sensibilité d'antenne A- 1(R'a b) par la fonction de propagation inverse r'a,b exp(j2lrr'ab/X), pour chaque segment, pour trouver B-lab(s,t) (S4). On multiplie à son tour ce produit par une fonction de pondération (par exemple coslRI/2c)) (S5), et on ef- fectue une transformation de Fourier rapide (FFT) (S6). Après avoir mul- tiplié les informations de sortie de S2 et S6 (S7), on effectue une trans- formation de Fourier rapide inverse (S8), après quoi on effectue au niveau de segment une compensation de la fonction de pondération qui est utilisée pour multiplier les données de mesure d'hologramme V(R) (S9), et on effectue au niveau de région de reconstruction une compensation de la fonction de pondération qui est utilisée pour multiplier B'lab(s,t) (S10). L'information de sortie de S10 est enregistrée à la position du segment a, b et elle est visualisée et prise pour l'image reconstruite du segment Sl 1. Les opérations décrites ci-dessus sont accomplies pour tous les segments de façon à produire l'image reconstruite d'hologramme (S12). Les segments divisés de la région de reconstruction de l'image reconstruite d'hologramme peuvent se chevaucher. Ceci élimine l'utilisa- tion des données aux bords coupés o l'erreur tend à être grande, ce qui permet la génération d'une image reconstruite de chaque segment dont la précision s'améliore avec la proximité vis-à-vis du centre du segment. On décrira ensuite les résultats d'une étude concernant le do- maine d'application ainsi que la résolution et la précision de ce mode de réalisation de la présente invention, qui a été effectuée par simulation par ordinateur. L'antenne de réception qui est utilisée dans la mesure d'holo- gramme utilise une antenne dipôle sur laquelle est fixée une plaque ré- fléchissante, pour réduire des changements de la directivité de l'antenne sous l'effet du système de balayage de l'antenne. L'antenne qui est réellement utilisée dans l'expérimentation comporte un câble coaxial se- mi-rigide qui fait saillie sur exactement X/14 à partir du centre d'une pla- que de base mesurant 2,3X x 2,33, qui est utilisée à titre de réflecteur. Le conducteur extérieur du câble coaxial en saillie est fendu en deux par deux fentes d'une longueur de./4, une extrémité du conducteur exté- rieure étant connectée au conducteur coaxial central, et l'autre extrémité étant simplement connectée à une antenne dipôle, remplissant ainsi la fonction d'un transformateur symétrique/dissymétrique. La figure 3 montre la structure et la distribution de courant d'un élément d'antenne de réception à balayage pour la mesure d'hologramme (f = 10 GHz), qui comprend un transformateur symétrique/dissymétrique et qui est analysée par un procédé de moment. La directivité totale de I'antenne a été obtenue à partir de cette distribution de courant, et la si- mulation par ordinateur d'une observation d'hologramme et d'une recons- truction d'image a été effectuée. Dans la simulation par ordinateur qui est décrite ci-dessous, la taille de la surface de mesure d'hologramme est D = 19,21, I'intervalle d'échantillonnage est 0,3X, le nombre d'échantillons est 64 x 64 = 4096, et l'étendue de la source d'ondes est L = 5D. La figure 4 montre les résultats d'une simulation par ordinateur qui a été accomplie pour comparer la résolution de l'image reconstruite, en ce qui concerne la position de la source d'ondes sur l'axe x, pour le procédé qui utilise une quasi-transformation de Fourier au moyen de l'approximation de Fresnel, employant la quasi-transformation de Fourier de la référence (9), et pour le procédé conforme à la transformation de Fourier avec division partielle du mode de réalisation de la présente in- vention. Dans le tableau, "Fresnel" désigne le procédé employant la quasi-transformation de Fourier avec l'approximation de Fresnel, "Complet" désigne le procédé par transformation de Fourier avec division partielle conforme au mode de réalisation de la présente invention, et le rapport d'amplitude de l'étalement d'image reconstruite est indiqué entre crochets. Ici, la source d'ondes est une source d'ondes ponctuelle dans la direction horizontale, placée à xs, Ys- Dans cette simulation, I'intervalle z' entre la surface de mesure d'hologramme et la surface de la source d'ondes est de 66,7X. Une largeur à -3 dB et une largeur à -20 dB de l'image de la source d'ondes sont indiquées pour la résolution. Cepen- * dant, dans un procédé utilisant une quasi-transformation de Fourier avec I'approximation de Fresnel, des données de mesure d'hologramme sont multipliées par exp(jU(x2 + y2)/Xz'), et la transformation de Fourier est effectuée, I'image reconstruite devenant ainsi une image dans l'espace des nombres d'ondes (spectre de Fourier), et l'erreur de coupure (apparition de données parasites à cause d'une discontinuité) des don- nées de mesure augmentant notablement et affectant l'image recons- truite. Les données de mesure d'hologramme doivent donc être multi- pliées par une fonction de pondération suffisamment lisse, par exemple (1+cos(2 1/2r IRI/D)). Pour la même raison, dans un procédé utilisant une quasi-transformation de Fourier, I'aire effective devient approximative- ment 64% de l'aire de surface de mesure d'hologramme D x D, et une résolution de -3 dB au voisinage de l'axe z devient AI = 1,25z'/D. On peut également comprendre d'après la figure 4 que, dans un cas dans lequel une source d'ondes ponctuelle est déplacée par rapport à l'axe z de fa- çon à avoir xS > z'/2, I'étalement d'une image à une résolution de -20 dB devient particulièrement grand. Dans ce procédé, la sensibilité d'antenne de la mesure d'hologramme qui est donnée par l'équation (5) n'est évi- demment pas prise en considération, et il en résulte que l'amplitude ab- solue de l'image reconstruite est inexacte, et à xS = 40X (120 cm), elle est approximativement 50% plus faible qu'à xs = 0. Dans le procédé du premier mode de réalisation de la présente invention, la pondération de l'ordre de (cos(rIRI/2112 D))1/2 pour des don- nées sur les bords, o l'erreur de coupure est élevée, est suffisante. Ce- pendant, du fait que cette pondération fait apparaître une région dans laquelle la reconstruction de l'image est impossible, on exécute un trai- tement de chevauchement pour la région de transformation, en étendant a et b de l'équation (19) à des nombres non entiers. Par ce procédé, l'aire effective devient 92% de l'aire de mesure d'hologramme D x D, et la résolution à -3 dB d'une source d'ondes ponctuelle au voisinage de l'axe z devient AI = 1,04k)z'/D. D'autre part, comme on peut le voir d'après la figure 4, dans des régions dans lesquelles une source d'ondes ponctuelle s'écarte de l'axe z, de façon à avoir xs > z'12, la résolution à -20 dB AI est améliorée d'environ 1/2 par rapport à un cas dans lequel on utilise l'approximation de Fresnel. La figure 5 montre les résultats de la simulation par ordinateur qui est accomplie pour étudier la probabilité d'amplitude d'une image re- construite par le procédé du premier mode de réalisation de la présente invention. Dans cette simulation également, on peut voir que même lors- que z' = 66,7X et xS = 40. (120 cm), I'erreur de niveau d'amplitude de l'image reconstruite ne dépasse pas 7%. Dans les résultats de simulation qui sont représentés sur la figure 5, on peut considérer que l'erreur de niveau est relativement grande lorsque la position de la source d'ondes est au voisinage de x5 = Ys = 40> (120 cm), du fait que dans la conver- sion faisant passer de l'équation (16) à l'équation (18), le changement de la directivité de l'antenne de réception dans la plage de + D/2 dans la surface de la source d'ondes, lorsqu'on observe à partir du centre de la surface de mesure d'hologramme, ne peut pas être considéré comme né- gligeable par rapport au changement dans exp(j2itr'mn/X). La figure 6 montre les résultats de la simulation par ordinateur qui est accomplie pour étudier la relation entre la résolution de la re- construction d'hologramme et la distance z' depuis la surface de mesure jusqu'au point de source d'ondes, dans le procédé du premier mode de réalisation de la présente invention. D'après ces résultats de simulation, on peut voir que AI est proportionnel à z' lorsque z' > D. Par exemple, à une résolution de -3 dB, AI est pratiquement égal à Xz'/D. Ces résultats indiquent que la résolution angulaire est pratiquement fixe. De plus, le changement de résolution par rapport à z' diminue lorsque z' que cet état conduit à l'entrée dans une région dans laquelle la compen- sation de l'influence de la directivité de l'antenne n'est plus possible. On décrira ensuite les résultats expérimentaux de l'observation d'hologramme. En utilisant un dispositif d'observation d'hologramme qui permet un certain niveau de fluctuation dans la fréquence d'observation au cours du temps, dans le but d'observer une fuite d'ondes électromagnétiques à partir d'un dispositif réel, on a observé l'émission par une antenne dipôle, on a effectué une confirmation expérimentale de la plage de visualisation ainsi que de la résolution et de la précision, et en outre, on décrit les ré- sultats de l'accomplissement d'expériences de visualisation de la fuite d'une source d'ondes électromagnétiques à partir d'un dispositif réel (armoire métallique). La figure 7 montre la structure d'un dispositif de mesure qui a été préparé en vue de l'expérimentation, pour observer des hologrammes radioélectriques réels. Sur la figure 7, la mesure d'hologramme est ac- complie à une surface d'antenne à balayage. Des signaux qui sont reçus par une antenne à balayage 3 font l'objet d'une sélection de fréquence dans un analyseur de spectre (Advantest R3271), ils sont convertis en une Fi (fréquence intermédiaire) de 21,4 MHz, et ils sont filtrés à une bande de 1 MHz. Pour annuler une fluctuation de phase au cours du temps pendant une mesure d'hologramme, une onde de référence est re- çue dans un analyseur de spectre supplémentaire, complètement asservi en phase. De plus, I'onde de référence est obtenue à partir d'une an- tenne fixe 4 ou d'une source de génération de bruit 6 directe. Le signal de mesure d'hologramme à 21,4 MHz et le signal de référence qui sont reçus aux deux analyseurs de spectre sont convertis en une fréquence intermédiaire et ils sont en outre convertis respective- ment à 15,4 MHz et 14,4 MHz, et ils sont multipliés dans un multiplieur. Dans le signal de sortie de ce multiplieur, les résultats de la multiplica- tion de la composante de fréquence du signal de mesure d'hologramme et de la composante de fréquence identique du signal de référence sont tous convertis en une seule fréquence de 1 MHz et, de plus, la phase de ce signal à 1 MHz contient une information qui est proportionnelle à la différence de phase entre le signal de mesure d'hologramme et le signal de référence. Ce signal à 1 MHz est échantillonné avec une horloge d'échan- tillonnage qui est synchronisée dans tous les oscillateurs locaux, et la distribution d'amplitude et la distribution de phase relative, par rapport à l'onde de référence à la surface de mesure d'hologramme, peuvent être obtenues par détection synchrone en utilisant une intégrale de Fourier. En prenant en considération comme il convient la résolution de recons- truction de l'image d'hologramme, la gamme de fréquence expérimentale de ce dispositif était de 2 GHz à 26,5 GHz. De plus, la plage de balayage de mesure d'hologramme était de 57 cm par 57 cm, et les pas de mouve- ment de l'antenne étaient inférieurs à X1/2, sur la base de la théorie de l'échantillonnage. La figure 8 montre les résultats de l'observation par antenne dipôle qui a été effectuée pour confirmer expérimentalement la précision de visualisation d'une source d'ondes électromagnétiques par le premier mode de réalisation de la présente invention. Les observations expéri- mentales ont été accomplies en utilisant des fréquences de 9 GHz et 10 GHz et des ondes modulées en fréquence d'une bande de 1 MHz. On a utilisé pour la source d'ondes une antenne dipôle, inclinée de 45 degrés par rapport à l'axe x, à une distance de z' = 200 cm de la surface d'ob- servation, et décalée à des intervalles de 30 cm, de -90 cm à 120 cm sur l'axe x. Sur la figure 8, la mesure d'hologramme et la reconstruction de source d'ondes sont effectuées, et les valeurs d'amplitude de l'image de source d'ondes sont portées pour chacune des diverses positions. Dans ces observations expérimentales, le signal de référence qui est employé a été produit en divisant un signal appliqué à une antenne dipôle. Ces résultats expérimentaux présentent une dépendance consi- dérable vis-à-vis de la fréquence et de la polarisation d'observation, et la précision des expériences n'est pas nécessairement élevée. L'erreur d'amplitude maximale est de 26%, et il apparaît une erreur considérable- ment plus élevée que dans les résultats de simulation qui sont représen- tés sur la figure 5. Cependant, en considérant les données pour 9 GHz, les résultats de simulation concordent dans une large mesure. On peut considérer que les sources du niveau d'erreur élevé comprennent une aberration dans la directivité de l'antenne dipôle qui est utilisée à titre d'objet d'observation, ainsi qu'une aberration dans la directivité de l'an- tenne qui est utilisée pour la mesure d'hologramme. Comme représenté sur la figure 9, les résultats expérimentaux concordent bien avec des ré- sultats de simulation en ce qui concerne la résolution de visualisation d'une source d'ondes ponctuelle. Par exemple, à une résolution de -20 dB, on a Al = 2,3 à 3,4 (Xz'/D). La figure 10 montre les résultats de la visualisation expérimen- tale dans le cas o l'on considère une fuite d'ondes électromagnétiques à partir d'un dispositif réel relativement grand (EST, c'est-à-dire Equipe- ment Sous Test). On a effectué la mesure d'hologramme à la fois pour un état dans lequel toutes les portes de l'EST étaient fermées, et pour un état dans lequel une seule porte était ouverte. La figure 11 montre une photographie de l'EST qui est utilisé à titre d'objet d'observation. Les mesures qui sont représentées sur la figure 10 ont été obtenues en pla- çant à l'intérieur de l'EST une source d'ondes à 10 GHz modulée en fré- quence, avec une bande de 1 MHz, et en effectuant ensuite des observa- tions d'hologramme à une position à l'intérieur de la pièce et à une dis- tance de 2 m de l'EST. De plus, la plage de mesure d'hologramme était de 57 cm sur 57 cm, les pas de balayage étaient de 0,9 cm, et la plage de visualisation était de 285 cm sur 285 cm. On a utilisé pour l'onde de référence une antenne fixe placée près de l'EST. Comme on peut le voir d'après les résultats de visualisation qui sont représentés sur la figure 10, dans l'état dans lequel une porte est ouverte (côté droit), la composante Jx ainsi que la composante Jy sont largement réparties dans l'ouverture de la porte, tandis que dans l'état dans lequel la porte est fermée, seule la composante Jx est concentrée dans l'interstice de la porte. On peut voir en outre qu'une partie ouverte est présente dans la partie inférieure de l'EST pour faire sortir des câ- bles, et que des émissions à partir de cette région sont particulièrement fortes lorsque la porte est fermée. On peut en outre voir que des émis- sions provenant de la partie ouverte dans la partie inférieure de l'EST sont réfléchies par la surface du plancher. Comme on l'a expliqué ci-dessus, la présente invention a l'effet de réduire le temps nécessaire pour obtenir des résultats de visualisation à partir de données de mesure d'hologramme V(R), par l'utilisation de la FFT avec division partielle. Par exemple, si le nombre d'éléments de données V(R) est N = 64 et si le nombre de points de visualisation est M = 640, le volume de calcul est approximativement le millième de celui qu'exige le procédé de l'art antérieur décrit ci-dessus, et le temps nécessaire pour obtenir les résultats de visualisation peut également être réduit à un niveau prati- que. De plus, en établissant un chevauchement dans la segmenta- tion de la région de reconstruction de l'image d'hologramme reconstruite, la présente invention a l'effet d'éviter l'utilisation de données aux bords, ayant un niveau élevé d'erreur de coupure, ce qui permet la génération d'une image reconstruite segment par segment, dans laquelle la précision s'améliore avec la proximité vis-à-vis du centre d'un segment. Pour confirmer l'efficacité de la présente invention, on a effec- tué une étude par simulation par ordinateur, concernant le domaine d'ap- plication ainsi que la résolution et la précision de la présente invention, dans laquelle un procédé de reconstruction d'hologramme numérique qui n'utilise pas l'approximation de Fresnel ou l'approximation de Fraunhofer a été élaboré à partir d'une transformation de Fourier par division par- tielle qui comprend une compensation de directivité de l'antenne de ré- ception, afin d'obtenir une résolution relativement élevée et une plage de visualisation étendue, même lorsque la surface de mesure d'hologramme est petite par rapport au dispositif qui constitue l'objet d'observation. Ensuite, en utilisant un dispositif d'observation d'hologramme qui permet un certain niveau de fluctuation de la fréquence d'observation au cours du temps, dans le but d'observer une fuite d'ondes électroma- gnétiques à partir d'un dispositif réel, on a observé l'émission d'une an- tenne dipôle et on a effectué une confirmation expérimentale concernant la plage de visualisation ainsi que la résolution et la précision, pour montrer ainsi finalement les résultats expérimentaux de visualisation de la fuite de source d'ondes électromagnétiques à partir d'un dispositif réel (armoire métallique), et pour prouver l'efficacité de la présente invention. On décrira ensuite le second mode de réalisation de la présente invention en se référant aux dessins annexés. Comme représenté sur la figure 12, I'observation d'hologramme est accomplie sur deux surfaces de balayage, à savoir une surface de balayage (1) et une surface de balayage (2). Deux antennes-sondes, comprenant une antenne-sonde à polarisation horizontale et une an- tenne-sonde à polarisation verticale, sont placées sur chaque surface de balayage, et les vecteurs de tension qui sont reçus par ces antennes sont: V (R) = Vh (R)1 =Vhl(X, y) (21) V1(R) = LVvl (R)/ LVvl(x, Y)] pour la surface de balayage (1); et V2(R) = Vh2 (R) = LVh2(x y)1 (22) LVv2 (R)J Vv2(X Y) pour la surface de balayage (2). Sur la figure 12, Ax et Ay indiquent la valeur de décalage de la surface de balayage (1) et de la surface de balayage (2). Les angles de coordonnées polaires entre chaque antenne pla- cée au point R = (x, y; z = 0) sur la surface de balayage (1) et le vecteur de courant de source d'ondes du point R' = (x', y', z'), sont eY = -- + t (24) eO= =-tan1[]2 (25) )= tan-1 [Y (26) eh = - QX (27) év = Oc- e (28) oh= ±tan (29) Z2 +(y _,) + ta-' Y-Y (24) Vy 2 + a-1FZ2 + (XX2 Ir z oz + tan(-' y) (25) fvz =+tan1[ Y-Y] (2630) Les angles de coordonnées polaires entre chaque antenne pla- cée au point Rh = 0) sur la surface de balayage (2) et le vecteur eV = 0 y (28) óh = -+ tan-'1yY,(9 2 L 4)v = - + tan-'1 (30) Les angles de coordonnées polaires entre chaque antenne pla- cée au point R = (x, y; z = 0) sur la surface de balayage (2) et le vecteur de courant de source d'ondes du point R' = (x', y', z') sont obtenus en remplaçant x par x + Ax et en remplaçant y par y + A dans les équations allant de l'équation (12) à l'équation (19), ce qui donne eOx', Oy', ez', z', oh' ev', Xh' et kv'- Les caractéristiques directionnelles globales des antennes- sondes qui sont utilisées dans l'observation d'hologramme sont calculées par exemple par un procédé de moment, et si la sensibilité de polarisa- tion principale de ces antennes-sondes est AO et la sensibilité de polari- sation croisée est A+, la matrice de sensibilité d'antenne de réception de I'antenne de réception à polarisation horizontale, conformément au vec- teur de courant de source d'onde, est AO(Oh, *h)àWIz AeOh. h)"Oy (A%(OhCO5 z±1(6hVsw)n4)Si Aoe-y- IA#^( ^bSh)^s A(O tJ 4.6r,)sinOz)iS (31O,)x3A( T)I)fO (31) et la matrice de sensibilité d'antenne de réception de l'antenne de récep- tion à polarisation verticale qui est due aux vecteurs de courant de source d'ondes, est [Ae(6Sh')sinOZ A,(81,'41L)sinO,,(Ae(8"'óh)c r + A, (O"'ó)sIn')]O Aoe-R') = JA(e,,^,)sine A.i(8r*,)siDe (A,(6v55)oe*+(vgiy+ t0)sin0g)SilB I Afl".4C)sin9S1 A(O',.S)sin9, (A 4) + 4 r) (32) En étendant à trois dimensions l'équation (7) de l'art antérieur, on ob- tient: J(R') = S -1(R - R'). V (R).exp(j2rr/X).rdxdy (33) à partir de quoi on peut trouver les trois composantes de vecteur de cou- rant de source d'ondes et réaliser la visualisation. Jx(R') Jx(x', y', z')1 J(R) = Jy (R') = Jy(x', y', z') (34) LJz(R')JLJz(x', y', z')j Dans ce cas, si IAH(R-R') I > I v(R-R') | (35) alors l'erreur de la matrice de sensibilité d'antenne de réception de l'an- tenne de réception à polarisation horizontale sera inférieure à l'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation verticale, ce qui fait que l'on a -I(R-R')-,=,H-I(R-R') (36) et: Vhl(R) [Vhl (X Y) V(R) = Vvl(R) =VvW(XY) (37) Vh2(R) Vh2(x+ Ax y + Ay).exp(j2i(r'-r) / X). r'/r D'autre part, si IAH(R-R') I alors l'erreur de la matrice de sensibilité d'antenne de réception de l'an- tenne de réception à polarisation verticale est inférieure à l'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation horizontale, ce qui fait que l'on a A-1(R-R') _- v-1(R-R') (39) et: Vhl (R) Vh1 (X, Y) V(R) = V (R)= Vv (x Y) (40) Vv2(R) Vv2(x+ AX, y + Ay).exp(j2it(r'-r) /).rlr Cependant, dans les équations (37) et (38) r= J(x-x')2+(y_y,)2 +z2 (41) et: r'= j(x+Ax-x')2+(y+Ay-y')2+z2 (42) Comme décrit ci-dessus, en comparant les déterminants des matrices de sensibilité d'antenne de réception, en trouvant la grandeur d'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation verticale et d'erreur de la matrice de sensibi- lité d'antenne de réception de l'antenne de réception à polarisation hori- zontale, en utilisant sélectivement les vecteurs de tension reçus V(R) et la matrice inverse de sensibilité d'antenne de réception A-1(R-R'), on peut trouver les trois composantes du vecteur de courant de source d'on- des J(R'), et réaliser la visualisation. Du fait que la composante z, soit Jz, du vecteur de courant de source d'ondes J(R') ne peut pas être observée par l'avant, cette opération sélective est d'importance primordiale pour éviter le zéro de sensibilité correspondant à Jz. Le calcul du déterminant de la matrice de sensibilité d'antenne de réception dépend seulement du vecteur de position R' du point de source d'ondes, et il peut donc être déterminé sans effectuer le calcul pour tous les intervalles d'intégration. En d'autres termes, si IAH(0- R') I > IAv(O-R') I (43) alors l'erreur de la matrice de sensibilité d'antenne de réception de l'an- tenne de réception à polarisation horizontale est inférieure à l'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation verticale, ce qui fait que l'on a: A (R-R')- H1 (R-R') (44) ce qui conduit à Vhl (R) Vh 1 (yX Y) V(R) = Vvl alors l'erreur de la matrice de sensibilité d'antenne de réception de l'an- tenne de réception à polarisation verticale est inférieure à l'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation horizontale, et il en résulte que: 1 (R-R') _ Av- (R-R') (47) ce qui conduit à E Vhl(R).Vh1 (X, Y) V(R) = V(R)=Vv((XY) (48) Vv2(R) Vv2(x+ Ax y + ày).exp(j2r(r'-r) /).r'/r Comme décrit dans l'explication précédente, la présente inven- tion comprend: une étape pour accomplir une observation d'hologramme radioélectrique à deux surfaces de balayage, à savoir une première sur- face de balayage et une seconde surface de balayage, pour placer une antenne-sonde à polarisation horizontale et une antenne-sonde à polari- sation verticale sur chaque surface de balayage, et pour enregistrer la distribution de vecteurs de tension que reçoivent les antennes; une étape pour calculer toutes les caractéristiques directionnelles des antennes- sondes en utilisant un procédé de moment, pour trouver la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisa- tion horizontale, et la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation verticale, qui sont dues aux vec- teurs de courant de source d'ondes; et une étape pour comparer les dé- terminants des matrices de sensibilité d'antenne de réception, pour trou- ver la grandeur de l'erreur de la matrice de sensibilité d'antenne de ré- ception de l'antenne de réception à polarisation verticale et de l'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation horizontale, et pour utiliser sélectivement les vecteurs de tension reçus et la matrice inverse de sensibilité d'antenne de réception, et trouver les trois composantes de vecteurs de courant de source d'on- des et effectuer une visualisation, grâce à quoi la présente invention a l'effet de permettre la visualisation des trois composantes de vecteurs de courant de source d'ondes. En étendant de deux composantes à trois composantes la vi- sualisation de vecteur de courant de source d'ondes, la présente inven- tion a l'effet d'améliorer de + 6 dB à + 1 dB la précision de la visualisa- tion d'ondes. Il va de soi que de nombreuses modifications peuvent être ap- portées au procédé décrit et représenté, sans sortir du cadre de l'inven- tion. REVENDICATIONS 1. Procédé de visualisation d'une image d'une source d'ondes, mis en oeuvre par transformation de Fourier rapide (ou FFT) avec divi- sion partielle, pour accélérer l'holographie par le calcul, qui est un pro- cédé de reconstruction numérique de données de mesure d'hologramme, au moyen de l'observation d'un hologramme radioélectrique; caractérisé en ce qu'il comprend les étapes suivantes: on segmente la région de re- construction d'une image d'hologramme reconstruite, par la taille d'une surface de mesure d'hologramme; on effectue une première FFT (S2) après avoir multiplié des données de mesure d'hologramme par une fonction de pondération (S1); on exécute seconde FFT (S6) après avoir donné une distance focale de reconstruction et une paire de nombres dé- signant un segment de l'image reconstruite, et avoir multiplié par une fonction de pondération (S5) une valeur obtenue en multipliant une ma- trice inverse de sensibilité d'antenne pour chacun des segments par une fonction de propagation inverse (S4); on effectue une transformation de Fourier rapide inverse (S8) après avoir multiplié (S7) I'information de sortie de la première FFT par celle de la seconde FFT; on compense (S9), au niveau de chaque segment, la fonction de pondération qui est utilisée dans la multiplication des données de mesure d'hologramme; on compense (S10), au niveau de la région de reconstruction, la fonction de pondération qui est utilisée dans la multiplication de la valeur obtenue en multipliant la matrice inverse de sensibilité d'antenne pour chaque seg- ment par la fonction de propagation inverse; et on visualise (S11) I'infor- mation de sortie de cette étape à une position qui désigne le segment considéré, et on fait de l'information de sortie de cette étape une image de segment reconstruite; et on effectue pour tous les segments toutes les opérations d'écrites ci-dessus, et on produit une image d'hologramme re- construite (S12). 2. Procédé de visualisation d'une image d'une source d'ondes selon la revendication 1, caractérisé en ce qu'il est mis en oeuvre en uti- lisant une FFT par division partielle pour accélérer l'holographie par le calcul, la segmentation produisant des segments en chevauchement. 3. Procédé de visualisation d'une image d'une source d'ondes qui visualise trois composantes de vecteurs de courant de source d'ondes à partir de données de mesure d'hologramme qui sont obtenues au moyen d'une observation d'hologramme radioélectrique, caractérisé en ce qu'il comprend les étapes suivantes: on effectue une observation d'holo- gramme radioélectrique sur deux surfaces de balayage, à savoir une première surface de balayage et une seconde surface de balayage; on place au moins deux antennes-sondes (3), à savoir une antenne de ré- ception à polarisation horizontale et une antenne de réception à polari- sation verticale, sur chacune des surfaces de balayage, et on enregistre la distribution de vecteurs de tension qui est reçue par les antennes (3); on calcule des caractéristiques directionnelles totales des antennes- sondes (3), en utilisant un procédé de moment, et on trouve une matrice de sensibilité d'antenne de réception de l'antenne de réception à polari- sation horizontale, et une matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation verticale, qui sont produites par des vecteurs de courant de source d'ondes; et on compare des détermi- nants des matrices de sensibilité d'antenne de réception, on trouve une grandeur d'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polarisation verticale, et d'erreur de la matrice de sensibilité d'antenne de réception de l'antenne de réception à polari- sation horizontale, et on utilise sélectivement les vecteurs de tension re- çus et des matrices inverses de sensibilité d'antenne de réception, pour trouver et visualiser trois composantes de vecteurs de courant de source d'ondes.