Brevet FR2614450A1 — publié le 19881028

-- 1 4--5 PROCEDE ET APPAREIL D'EXAMEN DE MILIEUX PAR ECHOGRAPHIE ULTRA- SONORE La présente invention concerne un procédé d'examen de milieux par échographie ultrasonore, comprenant d'une part l'émission d'ondes ultrasonores vers le milieu à explorer à l'aide d'une structure transductrice excitée à cet effet avec des fréquence, durée et rythme d'émission déterminés, et d'autre part la réception et le traitement des signaux écho- graphiques correspondant aux échos renvcyés vers ladite struc- ture transductrice par les obstacles rencontrés dans le milieu exploré, ledit traitement des signaux échographiques associés à chaque ligne échographique du plan de coupe dans le milieu exploré comprenant une amplification desdits signaux, leur correction des effets de l'atténuation en fonction du temps ou de la profondeur d'exploration, et leur visualisation pour constituer une image de référence. L'invention concerne égale- ment un appareil pour la mise en oeuvre de ce procédé. Le but de l'invention est de proposer un procédé et un appareil d'échographie ultrasonore faisant appel à un mode de traitement original des signaux échographiques de ré- ception, en vue de l'obtention d'images plus facilement ex- ploitables des milieux explorés. A cet effet, l'invention concerne un procédé re- marquable en ce que ledit traitement comprend également les étapes suivantes, réalisées successivement-pour chaque ligne échographique ou pour des segments spatialement voisins d'un ensemble de lignes échographiques juxtaposées: (a) une conversion analogique-numérique des si- gnaux échographiques corrigés, avec mise en mémoire des si- gnaux ainsi numérisés; (b) une transformation de signaux échographiques comprenant successivement pour chaque ligne échographique le -2- découpage de celle-ci en n éléments de p points, avec réarran- gement desdits n x p points en une matrice X de n lignes et p colonnes et mémorisation de cette matrice, le calcul et la mé- morisation du produit matriciel de X par sa matrice transposée X', après détermination de la moyenne de chaque colonne de X et soustraction de chacun des termes de la matrice de la mo- yenne correspondant à sa colonne, le calcul et la mémorisation des valeurs et vecteurs propres dudit produit matriciel, avec classement des valeurs propres selon un ordre croissant ou dé- croissant, et, avant envoi vers la deuxième mémoire-tampon, la reconstitution des éléments de lignes échographiques par sé- lection de quelques facteurs principaux impliquant une recons- titution de la matrice somme des séries de produits matriciels unitaires relatifs auxdits facteurs principaux sélectionnés; (c) après mémorisation des signaux ainsi reconsti- tués, une conversion numérique-analogique; (d) une visualisation des signaux convertis, de façon similaire à celle de l'image de référence. Malgré les avantages de l'échographie ultrasonore, tels que l'absence de radiations, la robustesse des équipe- ments actuels, etc..., la complexité des interactions entre les ondes ultrasonores et les milieux explorés -qui sont des tissus biologiques, dans le cas de l'échographie à application médicale- et le nombre élevé des paramètres physiques entrant en jeu empêchent certainement de tirer pleinement profit des possibilités de cette méthode pour l'analyse et la connais- sance desdits milieux. Le procédé proposé conformément à l'invention a pour objectif de contourner cette limitation. Il permet en effet d'extraire des signaux de réception complexes des infor- mations choisies à volonté, jugées importantes par le prati- cien, ces informations conduisant à des images de même nature que les images classiques obtenues jusqu'à présent mais d'in- terprétation plus facile. Les particularités et avantages de l'invention - 3 - apparaîtront maintenant de façon plus précise dans la descrip- tion qui suit et dans les dessins annexés, donnés à titre d'exemple non limitatif et dans lesquels: - la figure 1 montre le découpage d'une ligne échographique et son réarrangement en une matrice de n x p points; - la figure 2 montre un exemple de réalisation d'un appareil selon l'invention. L'analyse factorielle constitue une méthode de statistique descriptive multidimensionnelle qui utilise des calculs d'ajustement faisant appel à l'algèbre linéaire. Une telle méthode permet de fournir des représentations synthé- tiques d'ensembles de valeurs numériques, les objets à décrire devenant en effet des points dans un espace à plusieurs dimen- sions. Le principe essentiel de l'analyse factorielle est la recherche de la reconstitution d'un ensemble de n x p valeurs xij à partir d'un plus petit nombre de ces valeurs. Subsidiairement, une telle méthode d'analyse per- met non pas la reconstitution du tableau complet des n x p va- leurs, mais la détermination de q x (n+p) d'entre elles, avec bien entendu la relation q x (n+p) est suffisamment petit. Le principe général est le suivant. Etant donné une matrice X à n lignes et p colonnes que l'on veut recons- tituer à partir d'un nombre de valeurs numériques plus petit que n x p, ladite reconstitution de X avec q x (n+p) valeurs est possible si l'on peut trouver q vecteurs-colonnes à n composantes et q vecteurs-colonnes à p composantes, notés respectivement u et v, tels que X = uv. On a dans ce cas reconstitué les np valeurs de X avec seulement n+p valeurs numériques, ce qui est la décomposition la plus simple possible. Le problème ainsi résolu-peut recevoir une inter- prétation géométrique dans un espace à p dimensions RP. -- 4 - -4- Les n lignes de X peuvent être considérées comme les coordon- nées d'un ensemble de n points de cet espace, qu'on appellera nuage. Si ce nuage de point est en fait entièrement contenu dans un sous-espace vectoriel Rq de RP, il est possible de retrouver les positions des n points à partir de leurs coor- données sur les q nouveaux axes de ce sous-espace et, bien en- tendu, à partir des composantes de ces nouveaux axes. Le processus d'ajustement par un sous-espace vec- toriel dans RP est réalisable de la manière suivante. On re- cherche d'abord la droite F1 passant par l'origine des p axes qui définissent RP et qui, simultanément, ajuste au mieux le nuage des n points représentant dans RP les lignes de X. Si M1 est le point de RP qui représente la première ligne de X, et si u est le vecteur unitaire porté par la droite F1 et tel i = p que u'u = 1 ou encore E u2 = 1, alors les n lignes du i=1 i vecteur-colonne Xu sont les n produits scalaires de ces points avec u, et représentent par conséquent les longueurs des seg- ments déterminés par projection de ces n points sur la droite Fi. Pour que la droite ajuste au mieux le nuage des n points, il faut donc minimiser la somme de ces segments E MjH2 (j variant de 1 à n), expression dans laquelle M1, M2,..., Mn sont les n points de RP représentant les lignes de X, et H1, H2,..., Hn leur projection sur la droite FI. Puisque l'expression E OM2 est fixée (0 étant l'origine des p j axes dans RP), rendre l'expression E MjH2 minimale revient à rendre l'expression E OH? maximale, c'est-à-dire à recher- cher le vecteur unitaire u qui rend maximale la quantité (Xu)'Xu = u'X'Xu (X' = transposée de X). On pourrait montrer de façon détaillée -mais on l'admettra-,que ce sous-espace vectoriel optimal à une di- mension, appelé alors ul, qui ajuste au mieux au sens des - 5 - moindres carrés le nuage des n points, existe. On montrerait de même que le meilleur sous-espace vectoriel à deux dimen- sions réalisant cet ajustement existe également et qu'il con- tient ul. On trouve ce sous-espace en cherchant le vecteur unitaire u2 orthogonal à ul (c'est-à-dire tel que u'2u1 = O et u'2u2 = 1) qui rend maximale la forme quadratique u'2X'Xu2. De façon analogue, en raisonnant par récurrence, on voit que le sous-espace vectoriel optimal à q dimensions (avec q inférieur ou égal à p) est engendré par des vecteurs ul, u2, u3,..., uq_1, Uq o uq est orthogonal à ul, u2, u3,...., Uq-1 et rend maximale la forme quadratique u'rX'Xuq, avec u'quq = 1. On démontre également queu1 est le vecteur propre de la matrice X'X qui a la plus grande valeur propre À1, et de même pour u2, u3, etc... avec des valeurs propres A2, A3,... Aq, en considérant que A2 est la valeur propre immédiatement inférieure ou égale à A1, A3 la valeur propre immédiatement inférieure ou égale à A2, et ainsi de suite. Finalement, les q vecteurs propres normés correspondant aux q plus grandes va- leurs propres de la matrice symétrique X'X constituent une base orthonormée du sous-espace vectoriel Rq de RP qui s'ajuste au mieux au nuage des n points. La démonstration précédente peut être reprise de façon similaire dans un espace à n dimensions Rn, les p co- lonnes de X pouvant cette fois être considérées comme les coordonnées d'un ensemble, ou nuage, de p points de cet es- pace. On constate de même qu'une base orthonormée du sous- espace vectoriel à q dimensions s'ajustant au mieux, au sens des moindres carrés, au nuage des p points est constituée par q vecteurs propres vl, v2, v3,..., Vq-1, vq correspondant aux q plus grandes valeurs propres P1, P2, P3,..., Pq-1, pq de la matrice symétrique XX', le vecteur cherché devant rendre maximale la forme quadratique v'XX' v avec v'v = 1. On démontre enfin que, quel que soit a inférieur ou égal au rang r de X et donc de XX', on a les relations sui- vantes: - 6 - Ua =; X'V" Va = X vU va = P uc Une reconstitution approchée X* du tableau des données numé- riques X est alors donnée par l'expression: a =q a=q X a Da vau a= 1 (Xua) UIa Cette express:on montre qu'en fait le terme Xu. représente le poids du vecteur u. dans la reconstitution approchée de X. Le calcul des différents termes Xua permet donc de con- naître l'importance relative de tel ou tel facteur dans ladite reconstitution. L'analyse factorielle ainsi décrite va être main- tenant appliquée à l'échographie ultrasonore. En effet, les diverses lois de la physique qui sont impliquées dans un pro- cessus échographique conduisent malheureusement à des équa- tions très difficilement exploitables du fait de l'extrême im- brication des différents paramètres concernés. L'exploitation des images échographiques repose alors essentiellement sur la grande expérience des observateurs, médecins, etc... L'appli- cation de l'analyse factorielle à l'échographie permet de pro- poser un processus plus objectif de traitement des données échographiques. On sait en effet que l'interaction d'ondes ultra- sonores longitudinales émises par un transducteur (ou une structure transductrice) ultrasonore avec des tissus biolo- giques rend possible la formation d'images de ces tissus et surtout des interfaces entre zones relativement homogènes de -7- ces tissus. Les échos renvoyés par ces interfaces qui marquent des ruptures d'impédance acoustique peuvent en effet être dé- tectés, puis traités de façon à localiser les zones échogènes, compte tenu notamment de la direction des ondes et du temps d'arrivée des fronts d'ondes. Il suffit alors d'explorer sé- quentiellement un plan de coupe des tissus pour obtenir la formation d'une image échographique à partir de lignes écho- graphiques successives enregistrées en déplaçant progressi- vement le transducteur ultrasonore sur le patient. Plus précisément, une image ultrasonore est tradi- tionnellement reconstituée, pour chacune des lignes qui la composent, à partir de l'enveloppe des cdonnées correspondant à chacune desdites lignes échographiques. Pour éviter la perte d'informations d'amplitude et de phase résultant de la déter- mination des enveloppes, on va ici appliquer l'analyse facto- rielle non pas à l'image globale, mais à chacune des lignes échographiques avant la prise d'enveloppe, ladite image glo- bale étant reconstituée ultérieurement à partir des lignes traitées par analyse factorielle pour en extraire les infor- mations essentielles. La méthode de traitement selon l'invention est la suivante. Chacune des lignes échographiques, qui représente (selon le mode de représentation dit mode A) l'amplitude des échos reçus en fonction du temps écoulé ou de la distance ou profondeur des tissus, est découpée en éléments de même lon- gueur, qui peuvent soit être juxtaposés, soit se chevaucher. Par exemple, une ligne est échantillonnée sur 2048 points, correspondant à une profondeur d'exploration de 6,4 centimè- tres, et découpée en 32 éléments juxtaposés de 64 points. Le principe de l'analyse factorielle appliqué dans ce cas con- siste à construire à partir d'un tel découpage le tableau de données sur lequel portera l'analyse et à effectuer cette analyse. La figure I illustre ce découpage des lignes écho- graphiques et la constitution du tableau des données. Le dé- coupage en n = 32 éléments de p = 64 points chacun donne lieu - 8 - à la constitution d'un tableau de n lignes de p points cha- cune. On prend comme point de départ du tableau le premier point de la ligne échographique. La longueur de 64 points par ligne matricielle du tableau de données a été ici retenue comme étant de l'ordre de la longueur de l'impulsion ultra- sonore incidente, mais une longueur inférieure ou, de' préfé- rence, supérieure peut bien entendu être choisie. On adopte par ailleurs, dans le cas présent, un pas de ligne égal au nombre de points par élément, ici 64. Il suffira alors de juxtaposer les lignes de la matrice traitée par analyse fac- torielle pour reconstituer une image correspondant à la tota- lité des n x p points de la ligne échographique d'origine. On a vu plus haut que la reconstitution du tableau de données initial est donnée par une expression comprenant des produits du type Xuu', la reconstitution d'un des n élé- ments à partir de la nouvelle base à q dimensions étant en effet réalisée par formation de la somme des produits Xuu' re- latifs aux valeurs propres retenues. La suppression, au cours de ladite reconstitution, de l'un ou de plusieurs de ces pro- duits revient à sélectionner de façon corrélative l'informa- tion que l'on désire plus spécialement conserver. En traitant ensuite les lignes (ainsi reconstituées) comme des lignes échographiques normales, on aboutit à la formation d'une image de même nature qu'une image classique, mais dont certaines composantes ont été volontairement éliminées. L'application d'une telle méthode permet par exemple d'éliminer des signaux échographiques d'origine tout ce qui correspond à du bruit, ce qui permet d'obtenir des images plus facilement interprétables et s'avère particulièrement utile dans le cadre d'applications médicales. La méthode de traitement de signaux échogra- phiques peut alors, de façon avantageuse, être mise en oeuvre selon les étapes suivantes: (1) réception et échantillonnage d'une ligne écho- graphique avec repérage de sa position, et mise en mémoire des -9- signaux résultant de l'échantillonnage et des paramètres de position; (2) réarrangement des données selon le découpage choisi, de façon à constituer la matrice X, et mise en mémoire de cette matrice; (3) calcul de la moyenne de chaque colonne de la matrice X et soustraction de chacun des termes de la matrice de la moyenne correspondant à sa colonne; (4) calcul et mise en mémoire du produit matriciel X'X; (5) calcul et mise en mémoire des valeurs propres et des vecteurs propres associés, dans le cas présent à l'aide d'un algorithme itératif utilisant la méthode dite de Jacobi, la mise en mémoire n'étant effectuée qu'après classement, par ordre de valeurs propres décroissantes par exemple; (6) reconstitution des éléments de lignes échogra- phiques à partir de quelques facteurs principaux choisis, par reconstitution de la matrice somme des produits Xuu' relatifs aux facteurs choisis, addition à chaque élément de la matrice ainsi obtenue de la moyenne calculée précédemment pour chaque colonne, et réarrangement des lignes de la matrice de manière à reconstituer la ligne échographique; (7) prise d'enveloppe de la ligne échographique et mise en mémoire; (8) visualisation à partir desdites enveloppes, avec notamment interpolation entre lesdites enveloppes; Pour la mise en oeuvre de cette méthode, un mode de réalisation préférentiel d'échographe ultrasonore va être maintenant décrit. L'appareil d'examen de milieux selon l'in- vention comprend, dans l'exemple décrit, d'une part des élé- ments de type classique, tels que la structure transductrice, l'étage d'émission et, dans l'étage de réception et de trai- tement, le dispositif de visualisation, et d'autre part les éléments originaux, à savoir un dispositif de traitement des signaux échographiques qui est prévu dans ledit étage de 261.... - 10 - réception, en parallèle sur le dispositif de visualisation. Plus précisément, l'appareil selon l'invention comprend ici un transducteur ultrasonore 101 dont l'excitation est commandée par un étage d'émission 102 composé essentiel- lement d'un circuit d'excitation de transducteur 104 et d'un circuit d'aiguillage 103. Une horloge interne à l'étage d'émission définit les rythme et fréquence d'émission des ondes ultrasonores par le transducteur, et définit les durées des périodes d'émission et de réception. A l'émission, le cir- cuit d'excitation 104 et le transducteur 101 sont reliés par l'intermédiaire du circuit d'aiguillage 103 qui occupe sa po- sition émission. A la réception, les signaux échographiques fournis par le transducteur en réponse aux échos qu'il reçoit sont envoyés, par l'intermédiaire du circuit d'aiguillage 103 venu en position réception, vers l'étage de réception com- prenant successivement un amplificateur 106, un circuit 107 de correction de l'atténuation des ondes ultrasonores dans les milieux traversés (en fait, un amplificateur à gain croissant en fonction du temps), et, dans deux voies en parallèle, le dispositif de visualisation de type connu et le dispositif de traitement selon l'invention. Le dispositif de visualisation, référencé 200, comprend ici un redresseur 201, un filtre passe-bas 202, et un étage 203 incluant tous les circuits nécessaires à la consti- tution d'une image à partir des données fournies pour chacune des lignes échographiques concernant le plan de coupe balayé. Un dispositif de visualisation de ce type est classique dans un échographe, et lesdits circuits ne seront donc pas décrits plus en détail. Le dispositif de traitement selon l'invention com- prend un dispositif de calcul 400 suivi d'un dispositif de vi- sualisation 300 semblable au dispositif 200 et composé donc de circuits 301, 302, 303 identiques aux circuits 201, 202, 203. Ce dispositif de calcul 400 comprend d'abord un convertisseur analogiquenumérique rapide 401, par exemple un convertisseur - 1I - 8 bits/20 mégahertz tel le convertisseur 1048 de la société TRW, La Jolla, CA 92038, USA. Ce convertisseur 401 est suivi d'une première mémoiretampon 402 qui permet de stocker avant leur traitement l'ensemble des signaux concernant les lignes échographiques du plan de coupe balayé. Une unité de calcul programmable 403 reçoit alors la sortie de la mémoiretampon 402. Autour de cette unité de calcul, qui est par exemple un microprocesseur 68020 de la Société Motorola, 1303 East Algonquin Road, Schaumburg, Illinois 60196, USA, sont prévues une m-moire de travail 404, une mémoire 405 de stockage des programmes de traitement, une interface 406 praticien/unité de calcul qui constitue un or- gane d'entrée-sortie pour le choix des paramètres intervenant dans le traitement, une mémoire 410 de stockage de ces para- mètres, une première mémoire 407 à adressage matriciel pour le stockage sous forme matricielle de chaque ligne échographique en cours de traitement, une mémoire 411 de stockage des va- leurs moyennes de chaque colonne de la matrice stockée dans la mémoire 407, une mémoire intermédiaire 408 de stockage des valeurs propres et des vecteurs propres relatifs au traitement d'une ligne échographique, une deuxième mémoire 409 à adres- sage matriciel pour le stockage des signaux résultant des re- combinaisons partielles de lignes échographiques effectuées pour reconstituer celles-ci à partir de quelques facteurs principaux choisis par le praticien selon son expérience de l'interprétation des images et selon celles des informations qu'il veut mettre en évidence au vu de l'image de référence, et une deuxième mémoire-tampon 412 qui permet de stocker les lignes échographiques reconstituées, après leur traitement. En sortie de cette mémoire-tampon 412 est enfin prévu un conver- tisseur numérique-analogique 413, en sortie duquel les signaux sont envoyés vers le dispositif de visualisation 300. Le fonctionnement de l'appareil ainsi décrit est essentiellement le suivant. Après un balayage mécanique ou électronique correspondant à un plan de coupe déterminé dans - 12 - le milieu à explorer, les signaux échographiques disponibles en sortie du circuit 107 de correction des effets de l'atté- nuation suivent deux chemins parallèles. D'une part ces si- gnaux sont envoyés vers le dispositif de visualisation 200, de type classique comme on l'a vu, en vue de la formation d'une image traditionnelle de référence. D'autre part ces mêmes si- gnaux sont reçus par le dispositif de calcul 400 qui assure d'abord la numérisation des lignes échographiques, dans le convertisseur 401, puis leur stockage, dans la première mé- moire-tampon 402. Après ce traitement préliminaire, les lignes échographiques sont traitées l'une après l'autre. Grâce au logiciel de traitement stocké dans la mé- moire 405 de stockage des programmes de traitement, chaque ligne échographique est d'abord découpée, puis replacée sous forme matricielle X dans la première mémoire 407 à adressage matriciel. Ce logiciel stocké dans la mémoire 405 réalise, comme fonctions essentielles, la soustraction de chaque élé- ment de X de la moyenne correspondant à la colonne de X o se trouve cet élément (cette moyenne est stockée dans la mémoire 411), le produit matriciel XX', le calcul des valeurs et vec- teurs propres de la matrice ainsi obtenue, le classement de ces valeurs propres (par ordre de valeurs propres décrois- santes par exemple), et le stockage de ces valeurs et vecteurs propres. En fonction des numéros des facteurs principaux choi- sis a priori par le praticien (et stockés dans la mémoire 410 de stockage de paramètres), le logiciel exécute alors égale- ment les produits Xuu' correspondants, puis l'addition de ces produits et leur commande de stockage dans la deuxième mémoire 409 à adressage matriciel. A chaque élément de la nouvelle matrice ainsi constituée est alors rajoutée la moyenne stockée en mémoire 411 et correspondant à sa colonne. Enfin chaque ligne échogra- phique est reconstituée par une opération de réalignement in- verse de celles du découpage initial et de la mise sous forme matricielle, puis mémorisée dans la deuxième mémoire-tampon - 13 - 412 avant d'être envoyée par l'intermédiaire du convertisseur numériqueanalogique 413 vers le dispositif de visualisation 300. * Bien entendu, la présente invention n'est pas li- mitée à l'exemple de réalisation décrit, à partir duquel des variantes peuvent être proposées sans pour cela sortir du cadre de l'invention. En effet, l'appareil décrit ici, et qui fonctionne avec l'assistance d'un microcalculateur pour l'exé- cution de certaines fonctions, ne constitue pas la seule réa- lisation envisageable. On peut évidemment utiliser, en lieu et place de ce logiciel, les circuits numériques qui réalisent chacune de ses étapes, mais une telle utilisation est au prix d'une complexité accrue. Par ailleurs, les premier et deuxième dispositifs de visualisation de l'étage de réception et de traitement peuvent être regroupés dans un même ensemble ou étage de visualisation. Enfin, on remarquera bien que l'inven- tion a été décrite dans le cas du traitement d'une ligne écho- graphique, découpée ensuite en segments, mais que lesdits seg- ments peuvent également ne pas provenir de la même ligne écho- graphique, et être par exemple des segments spatialement voi- sins d'un ensemble de lignes échographiques successives. 261 4450 - 14 - REVENDICATIONS 1. Procédé d'examen de milieux par échographie ultra- sonore, comprenant d'une part l'émission d'ondes ultrasonores vers le milieu à explorer à l'aide d'une structure transduc- trice excitée à cet effet avec des fréquence, durée et rythme d'émission déterminés, et d'autre part la réception et le traitement des signaux échographiques correspondant aux échos renvoyés vers ladite structure transductrice par les obstacles rencontrés dans le milieu exploré, ledit traitement des si- gnaux échographiques associés à chaque ligne échographique du plan de coupe dans le milieu exploré comprenant une amplifica- tion desdits signaux, leur correction des effets de l'atténua- tion en fonction du temps ou de la profondeur d'exploration, et leur visualisation pour constituer une image de référence, caractérisé en ce que ledit traitement comprend également les étapes suivantes, réalisées successivement pour chaque ligne échographique ou pour des segments spatialement voisins d'un ensemble de lignes échographiques juxtaposées: (a) une conversion analogique-numérique des si- gnaux échographiques corrigés, avec mise en mémoire des si- gnaux ainsi numérisés; (b) une transformation de signaux échographiques comprenant successivement pour chaque ligne échographique le découpage de celle-ci en n éléments de p points, avec réarran- gement desdits nxp points en une matrice X de n lignes et p colonnes et mémorisation de cette matrice, le calcul et la mé- morisation du produit matriciel de X par sa matrice transposée X', après détermination de la moyenne de chaque colonne de X et soustraction de chacun des termes de la matrice de la moyenne correspondant à sa colonne, le calcul et la mémorisa- tion des valeurs et vecteurs propres dudit produit matriciel, avec classement des valeurs propres selon un ordre croissant ou décroissant, et, avant envoi vers la deuxième mémoire-tam- pon (412), la reconstitution des éléments de lignes échogra- phiques par sélection de quelques facteurs principaux impli- quant une reconstitution de la matrice somme des séries de - 15 - produits matriciels unitaires relatifs auxdits facteurs prin- cipaux sélectionnés; (c) après mémorisation des signaux ainsi reconsti- tués, une conversion numérique-analogique; (d) une visualisation des signaux convertis, de façon similaire à celle de l'image de référence. 2. Appareil d'examen de milieux par échographie ul- trasonore pour la mise en oeuvre du procédé selon la revendi- cation 1, comprenant une structure transductrice (101) asso- ciée à un étage d'émission (102) pour l'envoi d'ondes ultra- sonores vers le milieu à explorer et à un étage de réception et de traitement des échos renvoyés vers ladite structure transductrice par les obstacles rencontrés par les ondes ul- trasonores dans le milieu exploré, ledit étage de réception et de traitement recevant les signaux échographiques associés à chaque ligne échographique du plan de coupe dans le milieu exploré et comprenant au moins un amplificateur (106), un cir- cuit (107) de correction des effets de l'atténuation des ondes ultrasonores dans le milieu exploré, et un premier dispositif de visualisation (200) recevant les signaux corrigés présents en sortie dudit circuit, caractérisé en ce que ledit étage de réception comprend également en parallèle sur le premier dis- positif de visualisation un dispositif de calcul (400) suivi d'un deuxième dispositif de visualisation (300), ledit dispo- sitif de calcul comprenant lui-même un convertisseur analo- gique-numérique (401), une première mémoire-tampon (402) de stockage desdits signaux échographiques corrigés, des moyens (403 à 411) de transformation ligne échographique par ligne échographique de ces signaux échographiques corrigés, une deuxième mémoire-tampon (412) de stockage des signaux écho- graphiques transformés, un convertisseur numérique-analogique (413), et un deuxième dispositif de visualisation (300), la- dite transformation de signaux échographiques comprenant suc- cessivement pour chaque ligne échographique ou pour une ligne constituée de segments spatialement voisins d'un ensemble de - 16 - lignes échographiques juxtaposées le découpage de celle-ci en n éléments de p points, avec réarrangement desdits nxp points en une matrice X de n lignes et p colonnes et mémorisation de cette matrice, le calcul et la mémorisation du produit matri- ciel de X par sa matrice transposée X', après détermination de la moyenne de chaque colonne de X et soustraction de chacun des termes de la matrice de la moyenne correspondant à sa co- lonne, le calcul et la mémorisation des valeurs et vecteurs propres dudit produit matriciel, avec classement des valeurs propres selon un ordre croissant ou décroissant, et, avant en- voi vers la deuxième mémoire-tampon (412), la reconstitution des éléments de lignes échographiques par sélection de quel- ques facteurs principaux impliquant une reconstitution de la matrice somme des séries de produits matriciels unitaires re- latifs auxdits facteurs principaux sélectionnés. 3. Appareil selon la revendication 2, caractérisé en ce que les premier et deuxième dispositifs de visualisation sont inclus dans un même étage de visualisation.