Brevet FR2855271A1 — publié le 20041126

La présente invention concerne un procédé pour l'exploration et l'analyse d'une structure volumique par un traitement approprié de signaux représentatifs d'ondes, notamment d'ondes ultrasonores réfléchies ou transmises par cette structure volumique. Elle s'applique notamment, mais non exclusivement, à la réalisation d'appareils tels que des échographes, des appareils de contrôle non destructif d'objets, des sonars ou même des radars. Les appareils classiques de ce genre font habituellement intervenir des moyens 20 d'émission qui émettent une onde incidente dans le milieu à examiner et des moyens de réception utilisant éventuellement tout ou partie des moyens d'émission qui reçoivent les ondes réfléchies par les obstacles rencontrés par l'onde incidente. Des moyens sont en outre prévus pour traiter les signaux reçus par les moyens de réception et les présenter sous une forme exploitable 25 par l'utilisateur, par exemple sous la forme d'une image permettant de localiser la position des obstacles engendrant des réflexions de l'onde incidente. La méthode la plus courante consiste à utiliser des ondes impulsionnaires selon un processus consistant à transmettre une impulsion d'ondes ultrasonores 30 dans une direction donnée, à détecter le retour des échos, à mesurer le temps écoulé entre l'émission et la réception et à en déduire la distance, compte tenu de la vitesse de propagation de l'onde ultrasonore, et donc la position de l'obstacle qui a engendré chaque écho. Ce processus est ensuite répété dans différentes directions, selon une loi de balayage prédéterminé. Il devient ensuite possible de réaliser des images mettant en évidence les obstacles détectés par les échos, dont on connaît la position. Dans le but d'éviter les inconvénients des appareils utilisant un mode de traitement séquentiel (exploration ligne par ligne), on a proposé d'émettre dans la structure volumique à explorer une onde sensiblement plane, de section 10 relativement importante engendrée par une sonde constituée d'un réseau comprenant une pluralité d'éléments émission/réception de faibles dimensions, de préférence inférieure à la longueur d'onde ultrasonore, de façon à posséder un très large diagramme de rayonnement, ces éléments étant attaqués simultanément en parallèle. A la réception, chaque élément émission/réception 15 fonctionne de façon indépendante et reçoit donc séparément les ondes réfléchies par les obstacles interceptant le faisceau d'ondes ultrasonores qui se trouve dans sa zone de réception. Après numérisation, les informations délivrées par ces éléments émission/réception (champ d'ondes réfléchies) sont mémorisées dans des mémoires dont la lecture s'effectue en sens inverse de 20 l'écriture. Le brevet US-A-4 817 434 décrit un dispositif comprenant un générateur d'adresses par élément récepteur de la sonde, qui fournit l'adresse à lire dans une mémoire de champ correspondant au point image à reconstituer. Toutefois 25 ce dispositif ne permet qu'une cadence de reconstitution d'image relativement faible. Dans le but de supprimer cet inconvénient, le Demandeur s'est basé sur la constatation que dans un processus tel que celui précédemment décrit, chaque 30 point de l'objet à explorer donne naissance à une onde qui se trouve mémorisée à des adresses d'une mémoire de champ réparties sous la forme d'un arc de courbe, pseudo-hyperbolique, dont les caractéristiques dépendent de la position du point par rapport à la sonde et du diagramme de rayonnement de chaque élément (cette pseudo-hyperbole se réduisant théoriquement aux deux asymptotes pour les points situés contre la sonde). Il a mis au point un procédé comprenant les étapes suivantes: - l'émission dans ladite structure d'une onde incidente, - la réception des ondes réfléchies ou transmises par les obstacles rencontrés par l'onde incidente à l'intérieur de ladite structure, par une pluralité 10 d'éléments de détection indépendants les uns des autres, - la mémorisation, après numérisation, des signaux délivrés par les éléments de détection dans une mémoire de champ comprenant une ligne respective par élément de détection, et - la reconstitution et/ou l'analyse de la structure volumique à partir des 15 informations lues dans la mémoire de champ, dans laquelle sont calculées pour chaque point de la structure, les positions de la mémoire de champ contenant les signaux détectés par les éléments de détection, correspondant aux ondes réfléchies ou transmises par ce point, ces positions étant calculées à l'aide d'une loi d'adressage dont les paramètres dépendent de la 20 position de ce point par rapport aux éléments de détection, et dans laquelle pour chaque point, les lignes de la mémoire de champ sont lues aux positions respectives calculées au préalable pour ce point et stockées dans des mémoires d'adressage respectivement associées auxdites lignes de la mémoire de champ. Un calcul est alors appliqué aux informations lues pour ce point afin d'obtenir un résultat représentatif de l'importance de l'onde réfléchie ou transmise par ce point, lors de ce calcul, toutes les lignes de la mémoire de champ étant lues en parallèle pour chaque point aux positions indiquées pour ce point 30 respectivement par les mémoires d'adressage associées, le calcul du résultat étant ensuite appliqué à toutes les valeurs lues dans la mémoire de champ, ce résultat étant ensuite traité ou mémorisé dans une mémoire spécifique. Ce procédé décrit dans les brevets EP 0 825 453 B 1 et EP 0 872 742 B 1 5 déposés au nom du Demandeur, permet d'obtenir des performances remarquables en termes: - de vitesse d'examen et par conséquent autorise l'analyse de pièces de grandes dimensions, - de résolution spatiale compte tenu de la grande ouverture de rayonnement, 10 - de reproductibilité puisque le champ ultrasonore est émis par une onde plane. Néanmoins, l'expérience montre que dans certaines applications, des contraintes limitent les performances du procédé, en particulier, la vitesse de 15 traitement. La très haute résolution spatiale obtenue avec ledit procédé conduit à des images très fines et si l'on veut conserver la finesse des détails, et notamment mesurer avec précision l'amplitude des signaux réfléchis, la trame d'analyse 20 doit être très serrée, et par voie de conséquence, le nombre de points à calculer est élevé, ce qui ralentit la vitesse de traitement d'une manière excessive. En d'autres termes, l'amplitude maximale de l'écho peut ne pas être détecté dans la trame d'analyse s'il se situe entre deux pas d'analyse; cette erreur peut être rédhibitoire s'il s'avère nécessaire de mesurer l'importance de l'obstacle 25 avec une grande précision. On peut être conduit, en fonction de la précision demandée, à utiliser un pas d'analyse plus petit, donc un nombre de points calculés plus important, ce qui augmente le temps de calcul, et diminue la vitesse de contrôle. Il s'avère que dans le procédé décrit dans les brevets EP 0 825 453 B 1 et EP 0 872 742 B 1 déposés au nom du Demandeur, le pas d'analyse est indépendant du nombre de voies électroniques; il est souhaitable que le pas d'analyse soit un multiple ou un sous-multiple entier du pas inter éléments de détection; ainsi, dans le cas d'une sonde comportant 32 éléments de détection, espacés de 0,8 mm, associés à 32 voies électroniques, le pas d'analyse pourra être égal à 1,6 mm, 3,2 mm, 4,8 mm, 6,4 mm, etc... A titre d'exemple, on peut considérer le cas de sondage d'une plaque métallique d'épaisseur 120 mm par une sonde de 64 éléments au pas de 1 mm et une fréquence d'échantillonnage du signal détecté de 60 MHz. 10 La vitesse du son dans le matériau à sonder est de 6000 m/s; le temps de parcours, aller et retour, est de 40 [is. Sachant que le pas d'analyse du signal détecté est fixé à la période 15 d'échantillonnage et compte tenu du nombre de lignes horizontales équivalent au nombre d'éléments de la sonde, le nombre de points analysés est de x 60 x 64, soit 153 600. Le temps total de traitement est de 153 600/60 soit 2560 bts, soit une cadence 20 de 390 Hz. Il s'avère que la résolution latérale à 3 dB du procédé, à faible profondeur, est de l'ordre de 0,25 mm et que le pas d'analyse adopté dans ledit exemple de 1 mm, peut entraîner des erreurs supérieures à 10 dB sur l'amplitude mesurée; 25 cet écart est inacceptable dans la plupart des cas. Il est donc nécessaire de réduire le pas d'analyse d'un facteur 4, la cadence de traitement se situe en deçà de 100 Hz; cette cadence peut devenir incompatible avec la vitesse de contrôle exigée de pièces de grandes 30 dimensions. L'invention a donc plus particulièrement pour but de supprimer ces inconvénients grâce à un procédé permettant d'augmenter la vitesse de traitement sans sacrifier les autres performances, en modifiant les paramètres de traitement en fonction de la situation de la zone traitée de la mémoire de champ et d'optimiser la vitesse et/ou la précision de ce traitement. L'optimisation de la résolution et du temps de traitement repose sur le décodage des informations lues dans la mémoire de champ; ainsi, la courbe dite "de décodage" est la courbe de lecture des valeurs lues dans la mémoire 10 de champ; les positions de la mémoire de champ sont calculées à l'aide d'une loi d'adressage dont les paramètres dépendent de la position du point à analyser par rapport aux éléments de détection. La courbe, dite "inscrite", correspond aux signaux détectés par les éléments de 15 détection qui correspondent aux ondes réfléchies ou transmises par le point ou l'obstacle à analyser. La présente invention se base sur la constatation que, dans un processus tel que celui décrit précédemment, lorsque la courbe de décodage ne coïncide pas 20 exactement avec la courbe inscrite dans la mémoire de champ, la courbe de décodage peut détecter successivement des signaux positifs et négatifs. En conséquence, elle propose un procédé pour l'exploration et l'analyse d'une structure volumique, à partir des informations lues dans la mémoire de champ 25 par une courbe de décodage, ce procédé comprenant la détection du signe des signaux détectés, la détection d'une zone utile de ladite courbe de décodage dans laquelle les signaux sont de même signe et le calcul de la position du point analysé à partir d'une intégration des amplitudes des signaux détectés dans ladite zone utile de ladite courbe de décodage. 30 D'une façon plus précise, ce procédé pourra comprendre: - la détection des signaux positifs et négatifs, - la détection du centre de ladite courbe de décodage, - la mémorisation du signe des signaux au centre de la courbe de décodage, - la détection d'une zone utile de la courbe de décodage dans laquelle les signaux sont de même signe, - la sommation des signaux de même signe dans la zone utile de la courbe de décodage, - la division de ladite somme par le nombre de signaux analysés de même signe ou par une valeur constante prédéterminée correspondant à 10 l'ouverture totale de la courbe de décodage, - la mémorisation du résultat obtenu pour le point analysé de la structure volumique. Optionnellement, le procédé selon l'invention pourra mettre en oeuvre, à partir 15 des informations lues dans une trame d'analyse et au cours de la phase de mémorisation dans la mémoire de champ: - l'utilisation d'un pas horizontal variable en fonction de la profondeur d'analyse, et/ou - la modification du pas horizontal au voisinage de l'obstacle, et/ou - le traitement des signaux détectés préalablement à la phase de mémorisation. Un mode de mise en oeuvre du procédé selon l'invention sera décrit ciaprès, à titre d'exemple non limitatif, en référence aux dessins associés dans lesquels 25 La figure 1 représente les signaux détectés sur la courbe de décodage et le contenu de la mémoire de champ de part et d'autre de la courbe de décodage, dans le cas de la coïncidence de la courbe de décodage avec la courbe inscrite; La figure 2 représente les signaux détectés sur la courbe de décodage et le contenu de la mémoire de champ de part et d'autre de la courbe de décodage, dans le cas de la non coïncidence de la courbe de décodage avec la courbe inscrite; La figure 3 représente le résultat du décodage avec la prise en compte des signaux positifs et négatifs détectés sur la courbe de décodage; La figure 4 représente le résultat du décodage avec la prise en compte 10 uniquement des signaux positifs détectés sur la courbe de décodage; La figure 5 représente une courbe théorique de l'image d'un obstacle dans une trame d'analyse sans aucun dispositif de correction; La figure 6 représente l'enveloppe d'un signal numérisé en fonction du temps, lue suivant une ligne de champ mémorisé, avant traitement; La figure 7 représente l'enveloppe d'un signal numérisé en fonction du temps, lue suivant une ligne de champ mémorisé, après traitement; 20 La figure 8 représente un dispositif de traitement du signal numérisé, lu suivant une ligne de champ mémorisé, dont les résultats sont illustrés par les figures 6 et 7. Dans un premier exemple illustré par les figures 1 à 4, l'étape d'optimisation de la résolution spatiale, à partir des informations lues dans la mémoire de champ, comprend l'évaluation de l'écart horizontal entre la courbe de décodage et la courbe inscrite, et à corriger l'amplitude des signaux mesurés en fonction de cet écart; ainsi, lorsque la courbe de décodage ne coïncide pas 30 exactement avec la courbe inscrite, la courbe de décodage peut détecter successivement des signaux positifs et négatifs, la somme de ces dits signaux positifs et négatifs tendant vers 0. Réciproquement, lorsque la courbe de décodage coïncide exactement avec la courbe inscrite, la courbe de décodage détecte des signaux de même polarité. Ainsi, dans l'exemple représenté sur la figure 1, le contenu de la mémoire de champ est représenté par les amplitudes des signaux détectés et mémorisés Mll, M12, M13, M14; la courbe de décodage CD1 coïncide presque 10 exactement avec une courbe inscrite représentée par les signaux détectés et mémorisés M13 dans la mémoire de champ; ainsi, les amplitudes des signaux lus dans la mémoire de champ, représentés par M10 sont tous de même signe. Réciproquement, dans l'exemple représenté sur la figure 2, le contenu de la 15 mémoire de champ est représenté par les amplitudes des signaux détectés et mémorisés M21, M22, M23, M24; la courbe de décodage CD2 ne coïncide avec aucune courbe inscrite représentée par les signaux détectés et mémorisés dans la mémoire de champ; la courbe de décodage est décentrée et traverse successivement plusieurs courbes inscrites; ainsi, les amplitudes des signaux 20 lus dans la mémoire de champ, représentés par M20 sont tantôt positifs, tantôt négatifs. Soit S+ le nombres de valeurs positives, S- le nombre des valeurs négatives et SM le nombre total des valeurs positives et négatives sur une courbe de 25 décodage; par conséquent SM = (S+) + (S-). En supposant (cas théorique) que toutes les valeurs positives et négatives sont de même amplitude A, l'amplitude détectée, après sommation, est: A0 = SMxA. Dans le cas o la courbe de décodage ne coïncide pas avec la courbe inscrite, l'amplitude détectée, après sommation, devient A0 = A x (S+) - A x (S-); ce qui peut s'écrire: A0 = A x K, avec K =SM -2 x (S-) Le facteur K dépend de la différence entre les points de signes opposés; il sera égal à SM lorsque toutes les valeurs sont de même signe, et égal à zéro lorsqu'il y a autant de valeurs positives que de valeurs négatives. Il faut noter que le facteur K n'est pas constant du fait que les valeurs inscrites varient en fonction de leur position horizontale sur la courbe inscrite; en effet la directivité des obstacles réfléchissants provoque un maximum du signal détecté dans l'axe de visée de l'obstacle, et par conséquent un signal détecté plus faible sur la périphérie autour de l'axe de visée. 15 Ainsi, le facteur K doit être remplacé par une fonction f(k); cette fonction peut être déterminée théoriquement ou expérimentalement en déplaçant la sonde devant un obstacle connu et en mesurant pour chaque position de ladite sonde, l'amplitude du signal réfléchi et la valeur du facteur K correspondante. 20 Dans l'exemple représenté sur la figure 3, la coïncidence entre la courbe de décodage CD3 et les courbes inscrites M31, M32, M33, M34, n'est pas réalisée; en effet les signaux M30 changent de signe en s'écartant du centre de la courbe de décodage, et le nombre de signaux de même signe à partir du 25 centre est d'autant plus faible que le décalage est grand, de sorte que la somme des amplitudes décroît très rapidement. La valeur exacte de l'amplitude A peut alors être obtenue simplement en divisant A0 par K. Ce procédé a l'avantage de pouvoir être réalisé de façon 30 très simple, puisqu'il suffit de faire des additions sur le bit de signe du signal mémorisé. - 11 Un procédé très voisin, basé sur le même principe, mais dont la réalisation est légèrement plus complexe, peut être utilisé pour effectuer cette correction et réduire également les variations d'amplitude en fonction de la position des courbes de décodage par rapport aux courbes inscrites. En effet, si on additionne uniquement les signaux de même signe autour de la partie centrale, et si on divise cette somme par le nombre de signaux de même signe autour de cette partie centrale, on obtient une valeur très voisine S3, 10 sinon égale à celle de l'amplitude détectée en cas de coïncidence parfaite entre la courbe de décodage et la courbe inscrite. Ce processus peut être effectué d'une manière séquentielle, en mémorisant le signe du signal au centre de la courbe de décodage, puis en effectuant une 15 analyse de part et d'autre du centre, laquelle analyse est arrêtée en cas de détection de changement de signe du signal détecté ; ainsi, l'analyse consiste à effectuer la somme des signaux détectés divisée par le nombre de signaux analysés de même signe; cette valeur est mémorisée pour le point considéré; l'opération est ensuite effectuée de la même manière pour les autres points à 20 analyser, ce traitement peut être également effectué en parallèle par des circuits logiques. Par ailleurs, ce procédé amplifie de façon excessive des signaux correspondant à du bruit ou à des signaux réfléchis par des obstacles situés loin du point à 25 analyser, et l'image résultante risque d'être confuse. Pour éviter cet inconvénient, on considère qu'en dessous d'un certain nombre de signaux de même signe détectés à partir du centre, il ne s'agit plus d'un défaut, mais d'un bruit de fond, et la somme des ces signaux n'est plus alors 30 divisée par leur nombre, mais par une valeur constante, en général celle correspondant à l'ouverture totale de la courbe de décodage. - 12 Dans l'exemple représenté sur la figure 4, la coïncidence entre la courbe de décodage CD4 et les courbes inscrites M41, M42, M43, M44, n'est pas réalisée; dans le cas présent, seules les valeurs centrales de même signe ont 5 été prises en compte; la courbe résultante S4 présente une amplitude voisine de celle de la courbe résultante S3 et reste constante le long d'un déplacement horizontal plus important de part et d'autre du centre de la courbe de décodage. Ainsi, la précision sur l'amplitude du signal détecté est plus importante tout en adoptant un pas d'analyse plus grand. 10 Les différents processus susmentionnés contribuent à l'optimisation de la résolution spatiale, à partir des informations lues dans la mémoire de champ; l'association éventuelle de processus permettant l'augmentation de la vitesse de traitement des signaux dans la trame d'analyse contribue à l'amélioration 15 des performances du procédé décrit dans les brevets EP 0 825 453 B 1 et EP 0 872 742 B 1 déposés au nom du Demandeur, notamment dans l'analyse de pièces de grandes dimensions à cadence élevée. Dans cet exemple, selon l'invention, l'étape d'optimisation du temps de 20 traitement, à partir des informations lues dans la trame d'analyse et au cours de la phase de mémorisation dans la mémoire de champ, comprend en premier lieu l'utilisation d'un pas horizontal variable en fonction de la profondeur d'analyse. En effet, la résolution horizontale est fonction de la profondeur de l'obstacle; à titre d'exemple, une résolution de 0,25 mm à 3 mm de profondeur deviendra de 1 mm à 30 mm de profondeur. Il est donc envisageable d'adopter un pas horizontal variable en fonction de la 30 profondeur, ledit pas étant défini par la résolution acoustique à cette profondeur; ainsi, connaissant la profondeur des points à analyser, le pas - 13 d'analyse sera préalablement mémorisé et permettra un gain en temps d'analyse de la structure voisin de 2 à 3. L'étape d'optimisation du temps de traitement, à partir des informations lues 5 dans la trame d'analyse et au cours de la phase de mémorisation dans la mémoire de champ, comprend en second lieu la modification du pas horizontal au voisinage de l'obstacle. En effet, la détection de tous les signaux dont l'amplitude est supérieure au 10 bruit, avec un pas d'analyse important, permet la localisation approximative de ces dits signaux; une détection plus fine est ensuite effectuée dans la zone contenant chacun de ces dits signaux. Dans l'exemple représenté sur la figure 5, un signal S représente une courbe 15 théorique de l'image d'un obstacle dans une trame d'analyse. Ledit signal S est caractérisé par une amplitude maximale VM et un niveau de bruit d'amplitude maximale VB, ladite courbe théorique étant obtenue avec un pas d'analyse voisin de zéro. Soit VA le niveau d'analyse, inférieur à l'amplitude VM, considéré comme étant le seuil de détection des signaux à analyser; une erreur de détection inférieure à A = VM - VA, impose un pas d'analyse inférieur à P1, définit par l'intersection de la courbe S et du niveau d'analyse VA. De la même manière, pour détecter des signaux d'amplitude supérieure au niveau du bruit, le pas d'analyse sera nettement plus important, équivalent à P2, définit par l'intersection de la courbe S et du niveau de bruit d'amplitude maximale VB. Ainsi, la détection des signaux sera effectuée avec un pas d'analyse P2, définit comme étant le pas correspondant au seuil de détection situé audessus du - 14 niveau du bruit, permettant d'effectuer la recherche de signaux dont l'amplitude est supérieure au niveau d'analyse VA. Lorsqu'un obstacle est détecté à une certaine abscisse xp, l'analyse de cet 5 obstacle consiste à reprendre la détection dudit obstacle, à partir de l'abscisse Xp - P2 avec un pas d'analyse Pl. Ladite analyse est terminée au voisinage de l'obstacle, quand le signal détecté est inférieur au niveau de bruit VB et l'abscisse supérieure à xp. Il est donc envisageable d'adopter un pas horizontal variable en fonction des obstacles rencontrés, ledit pas étant relativement grand, permettant des cadences élevées; ainsi, ayant détecter un obstacle à analyser, le pas d'analyse sera réduit dans la zone de détection dudit obstacle; le temps de traitement n'est pratiquement pas affecté par l'analyse à pas réduit en regard du gain de 15 temps apporté par un pas initial élevé. L'étape d'optimisation du temps de traitement, à partir des informations lues dans la trame d'analyse et au cours de la phase de mémorisation dans la mémoire de champ, comprend en troisième lieu le traitement des signaux 20 détectés. En effet, dans chaque ligne de champ de la mémoire de champ sont mémorisés des signaux réfléchis ou transmis par un obstacle, sous forme échantillonnée; ainsi, l'amplitude des échantillons représente l'enveloppe desdits signaux 25 détectés; la fréquence d'échantillonnage est supérieure à la fréquence des signaux détectés de manière à détecter les extremums desdits signaux. Dans l'exemple représenté sur la figure 6, la période d'échantillonnage du signal détecté est environ dix fois plus faible que la période dudit signal 30 détecté. - 15 Le traitement du signal échantillonné consiste à détecter les extremums dudit signal échantillonné, à mémoriser les échantillons d'amplitude correspondante durant une demi période du signal détecté, et à mémoriser dans la ligne de champ de la mémoire de champ, non pas la totalité des échantillons dudit 5 signal échantillonné, mais uniquement les échantillons d'amplitude correspondant aux extremums. Dans l'exemple représenté sur la figure 7, le signal échantillonné est constitué des extremums du signal d'origine, représenté sur la figure 6; ainsi, le pas 10 d'analyse du signal traité peut atteindre la demi période du signal d'origine, tout en conservant la précision requise sur l'amplitude des signaux détectés. Le susmentionné traitement peut être effectué sous forme logicielle ou sous forme matérielle; dans l'exemple représenté sur la figure 8, le traitement est 15 effectué par des moyens matériels. Un bloc horloge H délivre un signal d'horloge SH qui est appliqué d'une part à un bloc registre à décalage RD à deux étages, et d'autre part à un bloc mémoire principale M2; le signal d'origine S, échantillonné à la fréquence 20 dudit signal horloge SH, est appliqué d'une part à l'entrée du bloc registre à décalage RD, et d'autre part à l'entrée d'une mémoire tampon Mi. Le bloc registre à décalage RD délivre des signaux SN et SN+1 correspondant à deux échantillons successifs dudit signal S; les susdits signaux SN et SN+, sont 25 appliqués aux deux entrées d'un comparateur C dont les deux sorties basculent de l'état 1 à l'état 0 suivant que l'échantillon N est plus grand ou plus petit que l'échantillon N+1; le basculement 1 vers 0 correspond à un extremum du signal S. - 16 Ledit basculement commande un circuit de détection D qui commande la mémorisation du signal S dans ladite mémoire tampon M1. La sortie de la mémoire M1 est ensuite appliquée à l'entrée de la mémoire principale M2. Ainsi, la mémoire principale M2 contient les valeurs des extremums à chaque demi période du signal d'origine S. L'ensemble des processus susmentionnés d'optimisation de la résolution spatiale, à partir des informations lues dans la mémoire de champ, et 10 d'optimisation du temps de traitement des signaux détectés peuvent être appliqués séparément ou d'une manière combinée; ils concourent ainsi à permettre d'accroître la vitesse de traitement des signaux détectés tout en conservant une très haute résolution spatiale, et de permettre une analyse tridimensionnelle des structures volumiques à cadence élevée. - 17 Revendications 1. Procédé pour l'exploration et l'analyse d'une structure volumique par un traitement approprié de signaux représentatifs d'ondes, notamment d'ondes 5 ultrasonores réfléchies ou transmises par cette structure volumique, ledit traitement consistant à reconstituer ou à analyser la structure volumique à partir des informations lues dans une mémoire de champ, dans laquelle sont calculés pour chaque point de structure, les positions de la mémoire de champ contenant les signaux détectés par les éléments de détection, correspondant 10 aux ondes réfléchies ou transmises par ce point, lesquelles informations dans la mémoire de champ sont lues par une courbe de décodage, caractérisé en ce qu'il comprend la détection du signe des signaux détectés, la détection d'une zone utile de ladite courbe de décodage dans laquelle les signaux sont de même signe et le calcul de la position du point analysé à partir 15 d'une intégration des amplitudes des signaux détectés dans ladite zone utile de ladite courbe de décodage. 2. Procédé selon la revendication 1, caractérisé en ce qu'il met en oeuvre: - la détection des signaux positifs et négatifs, - la détection du centre de ladite courbe de décodage, - la mémorisation du signe desdits signaux au centre de la courbe de décodage, - la détection d'une zone utile de la courbe de décodage dans laquelle les 25 signaux sont de même signe, - la sommation desdits signaux de même signe dans la zone utile de la courbe de décodage, - la division de ladite somme par le nombre de signaux analysés de même signe ou par une valeur constante prédéterminée correspondant à 30 l'ouverture totale de la courbe de décodage, - 18 - la mémorisation du résultat obtenu pour le point analysé de la structure volumique. 3. Procédé selon la revendication 1, caractérisé en ce qu'il comprend, à partir des informations lues dans une trame d'analyse et au cours de la phase de mémorisation dans la mémoire de champ - l'utilisation d'un pas horizontal variable en fonction de la profondeur, et/ou - la modification du pas horizontal au voisinage de l'obstacle détecté, et/ou - le traitement des signaux détectés préalablement à la phase de 10 mémorisation. 4. Procédé selon la revendication 3, caractérisé en ce que l'utilisation d'un pas horizontal variable en fonction de la profondeur comprend la mesure de la résolution spatiale en fonction de la 15 profondeur et la mémorisation de ladite résolution. 5. Procédé selon la revendication 3, caractérisé en ce que la modification du pas horizontal au voisinage de l'obstacle comprend une détection de tous les signaux dont l'amplitude est 20 supérieure au bruit, une localisation desdits signaux détectés, et une détection plus fine dans la zone contenant chacun desdits signaux détectés. 6. Procédé selon la revendication 3, caractérisé en ce que le traitement des signaux détectés préalablement à la 25 phase de mémorisation comprend la détection des extremums desdits signaux échantillonnés, la mémorisation des échantillons d'amplitude correspondante durant une demi période desdits signaux détectés, et la mémorisation dans la ligne de champ de la mémoire de champ des échantillons d'amplitude correspondant aux extremums.