FR 3161954 A1 20251107 FR 2404751 A 20240506 Procédé de mesure par tomographie par impédance électrique L’invention porte sur un procédé de reconstruction de la distribution des propriétés électriques d’au moins un matériau d’un corps comportant une partie cylindrique contenant un fluide, utilisant un réseau de neurones et des données de valeurs électriques du corps mesurées au préalable par tomographie par impédance électrique à l’aide d’électrodes agencées autour d’une périphérie de la partie cylindrique du corps, chaque électrode ayant été excitée par un potentiel de forme prédéfinie, le réseau de neurones comportant : un auto-encodeur comprenant plusieurs niveaux (L 1 ,…,L J ) comportant chacun au moins une couche de convolution et au moins une couche d’abandon, et un module d’attention comprenant des portes d’attention (AG 1 ,…, AG J-1 ), ayant chacune un signal d’attention (g j ) associé. Figure pour l’abrégé : Fig. 2 Procédé de reconstruction de la distribution des propriétés électriques des matériaux utilisant la tomographie par impédance électrique La présente invention concerne le domaine de la tomographie par impédance électrique. L’invention porte plus particulièrement sur un procédé de reconstruction de la distribution des propriétés électriques des matériaux d’un corps comportant une partie cylindrique contenant un fluide, utilisant un réseau de neurones et des données de propriétés électriques du corps mesurées au préalable par tomographie par impédance électrique. L’invention concerne également un produit programme d’ordinateur configuré pour mettre en œuvre ce procédé. L’application principale visée par la présente invention concerne la surveillance d’écoulements de fluides susceptibles de varier de façon abrupte, comme cela peut être le cas pour des fluides circulant sous haute pression et à haute température. Une application d’intérêt particulier concerne la surveillance de conduits d’installations nucléaires, mais d’autres applications peuvent être envisagées dans le cadre de l’invention. La Tomographie par Impédance Electrique (TIE) est une technique non invasive, non destructive, qui permet de visualiser, en temps réel et en continu, l’intérieur d’un objet en mesurant les propriétés électriques (potentiel et courant électrique) de l’objet à sa surface. Cette approche, robuste, est particulièrement adaptée à la réalisation de mesures non intrusives dans des environnements sous hautes pressions et/ou hautes températures. Par rapport à l'IRM et au scanner, la TIE offre une fréquence d'images rapide de plusieurs kHz, avec une résolution spatiale plus faible en raison de la non-linéarité du problème inverse. La TIE consiste plus précisément à injecter des courants ou potentiels électriques au moyen d’un ensemble d’électrodes non-intrusives disposées à la surface de l’objet surveillé et ensuite à mesurer les potentiels ou courants électriques à la surface de l’objet. Les électrodes peuvent être uniquement en contact avec la surface extérieure de l’objet. Toutefois, dans le cas où la surface de l’objet est en métal, les électrodes doivent traverser la paroi et être en contact avec le fluide. Le modèle d'excitation appliqué à travers les paires d'électrodes, avec les mesures de tension, résulte en une cartographie entre le courant et la tension. La carte de l’impédance à l’intérieur de l’objet est reconstruite en résolvant le problème inverse associé, afin de retrouver la distribution interne de la conductivité et de la permittivité, conformément à la loi d'Ohm. Initialement, les algorithmes utilisés pour résoudre le problème inverse non linéaire de la TIE étaient principalement des théories mathématiques itératives, telles que la régularisation de Tikhonov (TR), la rétroprojection linéaire (LBP), la méthode de Gauss-Newton (GNM), la méthode d'erreur en une étape de Newton (NOSER). Ces méthodes utilisent l'approximation linéaire de premier ordre pour conjuguer la cartographie non linéaire, ce qui conduit à une faible résolution et à une grande sensibilité au bruit pendant les mesures, sans compter que le temps de calcul exponentiel augmente avec les données. Au cours des trente dernières années, les méthodes utilisant les réseaux neuronaux artificiels ont permis des avancées majeures dans le domaine de la TIE, en raison de leur capacité à cartographier des relations non linéaires complexes avec une bonne convergence et de faibles erreurs. Cela se traduit par une résolution spatiale plus élevée et un temps d'exécution plus court, non seulement pour le problème inverse, mais aussi pour le débruitage, la super-résolution et la segmentation des images. La demande FR 3 121 234 décrit une méthode de reconstruction d'image utilisant l'algorithme NOSER et mettant en œuvre un multiplexage en fréquence dans lequel les signaux d’excitation sont imposés simultanément à toutes les électrodes. Pour permettre la discrimination des signaux, chaque électrode est excitée par un signal de forme trigonométrique. Le fait d’exciter toutes les électrodes simultanément évite la redondance que l’on retrouve dans les données obtenues en excitant séquentiellement les électrodes appairées. Cependant, la résolution spatiale obtenue n’est parfois pas suffisante pour détecter certains objets. L’article [1] propose le réseau ADALINE, basé sur un élément linéaire adaptatif. Les travaux décrits dans l’article [2] ont appliqué les réseaux de rétropropagation, et ceux de l’article [3] ont étudié les perceptrons multicouches bayésiens. Ces études ont été lancées sur des opérateurs de reconstruction de formation linéaire, en mesurant la différence entre les mesures de tension et en évitant la non-linéarité entre la tension et la conductivité. Pour surmonter la non-linéarité de la TIE, la méthode de l’article [4] utilise le réseau neuronal Radial-Basis-Functions, et celle de l’article [5] a exploré un réseau neuronal dense avec des données simulées d'EIDORS. Les études récentes sont davantage axées sur les architectures complexes, comme décrit dans l’article [6], où est appliquée une méthode d'auto-encodeur finement ajustée à la surveillance pulmonaire EIT, ou dans l’article [7] qui propose un auto-encodeur multicouche. Les chercheurs de l’article [8] ont étudié l'utilisation d’antécédents basés sur l'énergie avec un réseau neuronal informé par la physique pour améliorer l'apprentissage du modèle, tandis que ceux de l’article [9] ont utilisé l'apprentissage hybride-fusion pour la reconstruction à haute résolution. La demande WO 2022/77866 décrit une méthode d'imagerie par impédance électrique utilisant un ensemble original de données de tension mesurées sur une zone à tester. Une séquence de distribution de conductivité initiale de la zone forme un ensemble de données d'entraînement correspondant, auquel du bruit est ajouté. Un auto-encodeur variationnel est entraîné avec ces données d'entraînement. En fonction de l'ensemble de données de tension d’entrée, un encodeur de l'auto-encodeur variationnel est entraîné pour obtenir une caractéristique de données de tension d'entrée, et une relation de mise en correspondance est établie entre l'ensemble de données de tension d’entrée et une séquence de distribution de conductivité correspondante en utilisant un décodeur de l'auto-encodeur variationnel entraîné. Avec les progrès continus de l'apprentissage profond dans la TIE, de meilleurs résultats sont explorés par rapport aux algorithmes traditionnels. De façon connue, le problème inverse est résolu sans tenir compte de l'intérêt de donner des indications à l'intelligence artificielle (IA) sur les points à privilégier pour améliorer la résolution spatiale du résultat. L'amélioration de la reconstruction de l'image est pourtant primordiale. L'identification de la meilleure correspondance entre le courant ou le potentiel injecté et la distribution de la conductivité peut conduire à une faible résolution spatiale en raison du caractère mal posé du problème. Il existe donc un besoin pour proposer un procédé de traitement de mesures par TIE qui remédie aux inconvénients de l’art antérieur, notamment pour améliorer la résolution spatiale des images. Le but de l’invention est de répondre au moins partiellement à ce besoin. Pour ce faire, l’invention a pour objet, selon un de ses aspects, un procédé de reconstruction de la distribution des propriétés électriques d’un corps comportant une partie cylindrique au sein de laquelle s’écoule un fluide, utilisant un réseau de neurones et des données de valeurs électriques du corps mesurées au préalable par tomographie par impédance électrique à l’aide d’électrodes agencées autour d’une périphérie de la partie cylindrique du corps, chaque électrode ayant été excitée par un potentiel de forme prédéfinie, le réseau de neurones comportant : un auto-encodeur comprenant plusieurs niveaux L 1 ,…,L J comportant chacun au moins une couche de convolution et au moins une couche d’abandon, et un module d’attention comprenant des portes d’attention AG 1 ,…, AG J-1 , ayant chacune un signal d’attention g j associé, le procédé comportant les étapes suivantes : i) les données de valeurs électriques sont traitées successivement par chaque niveau L 1 ,…,L J de couches de convolution et de couches d’abandon de l’auto-encodeur du réseau de neurones de sorte à coder les données pour obtenir des données encodées x 1 ,…, x J pour chaque niveau de couches L 1 ,…,L J , ii) afin de calculer le coefficient αj d’une porte d’attention AG j , cette dernière est configurée pour additionner les données encodées x 1 ,…, x J au signal d’attention g j de la porte d’attention AG j correspondante, fonction au moins de la concaténation des données de sortie j+1 de la porte d’attention de rang supérieur AG j+1 et des données décodées provenant du niveau de rang supérieur L j+2 à celui de la porte d’attention de rang supérieur AG j+1 , le résultat de cette addition étant transformé par au moins une fonction mathématique de sorte à obtenir le coefficient α j , iii) les données de sortie j d’un niveau de couches L j de l’auto-encodeur sont multipliées par le coefficient α j de la porte d’attention correspondante AG j pour obtenir des données de sortie j de la première porte d’attention AG j , iv) l’étape précédente est renouvelée jusqu’à une première porte d’attention AG 1 afin d’obtenir les données de sortie 1 de la première porte d’attention AG 1 , v) lesdites données de sortie 1 de la première porte d’attention AG 1 sont concaténées avec les données décodées provenant du deuxième niveau L 2 afin d’obtenir une matrice de données représentative d’une image de la distribution des propriétés électriques d’au moins un matériau du corps. Grâce à l’invention, il est possible de reconstruire la distribution des propriétés électriques des matériaux, notamment la conductivité et la permittivité, avec une meilleure résolution spatiale. Le procédé selon l’invention utilise un réseau neuronal de bout en bout, pouvant être nommé Autoencoder Improved Attention-Net (AIA-Net). Le premier segment est un auto-encodeur utilisé comme architecture d'extraction de caractéristiques de remodelage pour l'entrée de tension, et le second segment est un module d’attention amélioré permettant de résoudre le problème inverse et fournissant l'image de distribution des propriétés électriques des matériaux sous test. Le potentiel d’excitation des électrodes étant de forme prédéfinie connue, il est possible de reconstruire la conductivité électrique et/ou la permittivité électrique. Grâce aux portes d’attention de l’invention, ayant la connaissance de la correspondance entre les valeurs électriques mesurées et les propriétés électriques des matériaux, la résolution spatiale de l’image en sortie est très satisfaisante. Auto-encodeur Les couches d’abandon, ou de « dropout », servent, de façon connue, à activer ou non les neurones. Cette technique permet de réduire l’overfitting lors de l’entraînement du modèle. Certains neurones peuvent être temporairement désactivés dans le réseau, ainsi que toutes les connexions entrantes et sortantes correspondantes. Le choix des neurones à désactiver est avantageusement aléatoire. Dans un mode de réalisation préféré, un échantillonnage par valeur maximale est mis en œuvre d’un niveau de couches L 1 ,…,L J à l’autre afin de compresser les données de valeurs électriques. Ce processus est appelé, de façon connue, « max pooling ». Le procédé de max-pooling s'effectue en passant une fenêtre, ou filtre, de dimension prédéterminée sur l'ensemble de l'entrée. La fenêtre se déplace sur l'image et sélectionne le maximum de chaque portion de l'image, appelée fenêtre ou région de pooling. La fenêtre glisse sur l'entrée et, pour chaque position de la fenêtre, la plus grande valeur est sélectionnée et placée dans la matrice de sortie. Des connexions de saut sont avantageusement réalisées entre les niveaux de couches L 1 ,…,L J pour obtenir les données encodées x 1 ,…, x J à l’entrée des portes d’attention. De façon connue, cette technique permet aux réseaux neuronaux convolutifs de contourner certaines couches et de se connecter directement à des couches plus profondes ou moins profondes. De façon préférée, à l’étape ii), les données encodées (x 1 ,…, x J ) sont, selon leur dimension spatiale, convoluées ou transposées pour s’adapter à la dimension spatiale du signal d’attention considéré. Module d’attention Les portes d'attention du module d’attention dans la partie décodage sont avantageusement utilisées comme filtres pour les informations croisées entre la partie codage, en transportant les caractéristiques à travers les connexions de saut, et la partie décodage. Ces informations croisées peuvent être améliorées en ajoutant aux couches supérieures les caractéristiques extraites de chaque couche dans la partie décodage. Des méthodes connues de tomographie par impédance électrique utilisent l’intelligence artificielle en se concentrant sur la correspondance non linéaire entre la tension et la conductivité, sans donner d'indications au réseau neuronal sur l'endroit où chercher. Le mécanisme d'attention a été introduit dans l’article [10], où il est formulé comme suit pour l'attention dite douce : De façon préférée, le coefficient αj de la porte d’attention AG j correspondante est calculé grâce aux équations suivantes : Où . est le quotient d'attention utilisé pour déterminer les coefficients d'attention α j , et sont respectivement des fonctions d'activation d’unité linéaire rectifiée, notée ReLU, et sigmoïde, , et sont des transformations linéaires destinées à réduire la consommation de paramètres et de ressources informatiques, combinées aux termes de biais , , and pour former un ensemble de paramètres . Le module d'attention permet de mettre en évidence les régions importantes pour l'extraction des caractéristiques lors de l'exécution des couches de convolution, et de supprimer les régions non pertinentes de l'image. De préférence, l'opération d'addition des données encodées est effectuée sur toutes les données d'entrée avec chaque signal d’attention , principalement afin d’augmenter les poids alignés et diminuer les poids non alignés, ce qui réduit l'information spatiale et met en évidence les caractéristiques importantes. Les caractéristiques extraites de l'entrée sont avantageusement sélectionnées grâce aux informations contenues dans les signaux d’attention . Les caractéristiques d'entrée x 1 ,…, x J sont avantageusement mises à l'échelle par rapport aux coefficients d'attention α j , afin d'identifier les zones significatives de l'image, après être passées par une interpolation trilinéaire , et . Dans un mode de réalisation préféré, la transformation ψ produit une matrice de poids qui est donnée à la fonction sigmoïde d’activation afin de limiter les valeurs des coefficients d'attention ( ) ∈ [0, 1]. Un rééchantillonneur peut être utilisé pour étendre les dimensions de la matrice de poids. À la sortie de chaque porte d’attention, une opération de concaténation peut être effectuée entre le signal d’attention g j et les données de sortie , étant donné qu'ils sont de dimensions différentes. La multiplication par éléments entre les caractéristiques d'entrée et les coefficients d'attention conduit avantageusement aux données de sortie ( ). Mesures préalables Dans un mode de réalisation préféré, la mesure tomographique par impédance électrique du corps, comportant une partie cylindrique contenant un fluide, afin d’obtenir les données de valeurs électriques du corps mesurées au préalable, comporte les étapes suivantes : i. agencement d’un nombre n e d’électrodes autour d’une périphérie de la partie cylindrique du corps, ii. excitation simultanée de chacune des n e électrodes, chaque électrode étant excitée par un potentiel V n exc ayant la forme : où A est une amplitude de signal, θ n est la position angulaire de l’électrode E n , f m = m * f 0 est une fréquence d’oscillation, f 0 est une fréquence fondamentale choisie telle que f m soit inférieure à la fréquence de Nyquist du système pour tout m, iii. mesure des propriétés électriques du corps à l’aide des n e électrodes, iv. traitement des données issues de l’étape iii de mesure, comportant les sous-étapes suivantes : a) pour chaque électrode E n , calcul des points de données M n définis par : où R est la valeur de la résistance utilisée pour la mesure de V n meas avec V n meas = R I n aux bornes de la résistance, P est le nombre de points d’une séquence discrète de mesure du courant I n , p est le temps discret, k est un coefficient de Fourier compris entre 1 et (n e – 1) et β p = (2πp/P). Utilisation du procédé L’invention concerne également l’utilisation du procédé qui vient d’être décrit pour la reconstruction de la distribution des propriétés électriques d’un conduit d’installation nucléaire, un écoulement diphasique traversant notamment le conduit. Produit programme d’ordinateur L’invention concerne également, sous un autre de ses aspects, un produit programme d’ordinateur comportant un support et, enregistrées sur ce support, des instructions lisibles par un processeur pour que, lorsqu’exécutées, elles permettent de mettre en œuvre le procédé de reconstruction selon l’invention. Dispositif L’invention concerne enfin un dispositif pour la mise en œuvre du procédé selon l’invention, comportant : un système d’acquisition comportant au moins un réseau logique programmable, un module de génération de signaux analogiques et un module de mesure de signaux analogiques, une pluralité d’électrodes étant notamment connectées au système d’acquisition, un ordinateur configuré pour contrôler le système d’acquisition, un réseau de neurones intégré ou relié à l’ordinateur, le réseau de neurones comportant : un auto-encodeur comprenant plusieurs niveaux L 1 ,…,L J comportant chacun au moins une couche de convolution et au moins une couche d’abandon, et un module d’attention comprenant des portes d’attention AG 1 ,…, AG J-1 , ayant chacune un signal d’attention g j associé. Les caractéristiques énoncées ci-dessus pour le procédé s’appliquent à l’utilisation, au produit programme d’ordinateur et au dispositif, et vice versa. FIG. 1 La FIG. 1 représente un diagramme illustrant un exemple d’étapes de mise en œuvre du procédé selon l’invention, FIG. 2 La FIG. 2 illustre un exemple d’architecture pour la mise en œuvre du procédé selon l’invention. FIG. 3 [ FIG. 3 ) Les figures 3a et 3b représentent en détails un exemple de portes d’attention utilisées dans l’invention, FIG. 4 La FIG. 4 illustre des résultats comparatifs entre le procédé selon l’invention et l’art antérieur, et FIG. 5 La FIG. 5 représente une courbe illustrant l’erreur quadratique moyenne obtenue lors de l’utilisation du procédé selon l’invention. Description détaillée d’un exemple On a illustré en FIG. 1 un diagramme représentant les étapes d’un exemple de mise en œuvre du procédé selon l’invention. Lors d’une première étape, des mesures expérimentales par tomographie par impédance électrique sont réalisées sur un corps comportant une partie cylindrique à l’intérieur de laquelle s’écoule un fluide, à l’aide d’électrodes agencées autour d’une périphérie de la partie cylindrique du corps, chaque électrode ayant été excitée par un potentiel de forme prédéfinie. Les électrodes sont disposées de façon non-intrusive sur une périphérie circulaire du corps, comme visible à la FIG. 2 . Les électrodes sont dans cet exemple angulairement réparties de manière régulière autour de la périphérie du corps. Les électrodes sont reliées à un circuit imprimé, lui-même connecté à un système d’acquisition de données. Un écran permet de visualiser les données et les images produites à partir de ces données. Le système d’acquisition de données contient le système d’exploitation linux (HÔTE) qui contrôle un réseau logique programmable FPGA, également contenu dans le système d’acquisition de données. Le système d’acquisition permet de générer les signaux analogiques d’excitation et de mesurer les signaux analogiques de mesure en provenance des électrodes. Des données de valeurs électriques M n (k) sont ainsi obtenues. Dans une deuxième étape, ces données sont traitées par l’auto-encodeur du réseau de neurones selon l’invention afin de les remodeler pour obtenir des données encodées x 1 ,…, x J pour chaque niveau de couches L 1 ,…,L J de l’auto-encodeur, comme visible à la FIG. 2 . Chaque niveau comporte au moins une couche de convolution et au moins une couche d’abandon, dite de « dropout ». Les données sont traitées par couches et concaténées entre elles. Comme visible à la FIG. 2 , un échantillonnage par valeur maximale, dit « max-pooling », est mis en œuvre d’un niveau de couches L 1 ,…,L J à l’autre afin de compresser les données de valeurs électriques. Dans une troisième étape, le module d’attention selon l’invention est utilisé pour la résolution du problème inverse. Ce module d’attention comprend des portes d’attention AG 1 ,…, AG J-1 , ayant chacune un signal d’attention g j associé. Comme représenté aux figures 3a et 3b, afin de calculer le coefficient αj d’une porte d’attention AG j , cette dernière est configurée pour additionner les données encodées x 1 ,…, x J au signal d’attention g j de la porte d’attention AG j correspondante, fonction au moins de la concaténation des données de sortie j+1 de la porte d’attention de rang supérieur AG j+1 et des données décodées provenant du niveau de rang supérieur L j+2 à celui de la porte d’attention de rang supérieur AG j+1 . Le résultat de cette addition est transformé par au moins une fonction mathématique de sorte à obtenir le coefficient α j . De préférence et dans l’exemple considéré, le coefficient α j de la porte d’attention AG j correspondante est calculé grâce aux équations suivantes : où est le quotient d'attention utilisé pour déterminer les coefficients d'attention αj, et sont respectivement des fonctions d'activation d’unité linéaire rectifiée, notée ReLU, et sigmoïde, , et sont des transformations linéaires destinées à réduire la consommation de paramètres et de ressources informatiques, combinées aux termes de biais , , and pour former un ensemble de paramètres . La transformation ψ produit une matrice de poids qui est donnée à la fonction sigmoïde d’activation afin de limiter les valeurs des coefficients d'attention, un rééchantillonneur est utilisé pour étendre les dimensions de la matrice de poids. Dans l’exemple illustré, comme visible à la FIG. 2 , des connexions de saut sont réalisées entre les niveaux de couches L 1 ,…,L J pour obtenir les données encodées x 1 ,…, x J à l’entrée des portes d’attention. Lors de l’étape du calcul des coefficients d’attention, les données encodées x 1 ,…, x J sont, selon leur dimension spatiale, convoluées ou transposées pour s’adapter à la dimension spatiale du signal d’attention considéré, comme représenté aux figures 3a et 3b. Par exemple, dans la porte d’attention AG 1 , la dimension spatiale de x 1 est supérieure, du double, à celle du signal d’attention g 1 . Pour pouvoir les additionner, il convient donc de diviser la dimension de x 1 par 2 (Conv2D avec stride = 2). En ce qui concerne les autres signaux x j en entrée, leur dimension spatiale est inférieure à celle du signal d’attention g 1 , sauf pour x 2 qui est au même niveau de dimension. Il convient ainsi d’effectuer une convolution inverse avec des pas différents (stride = 1 pour x 2 , 2 pour x 3 et 4 pour x 4 ) afin d’additionner les signaux avec g 1 . Le même principe se poursuit pour AG 2 , AG 3 , AG 4 . Sur cette figure, F, H, W, D, F int , sont respectivement les caractéristiques, la hauteur, la largeur, le canal d'entrée ou la profondeur, dans le cas de données 3D, et les caractéristiques intermédiaires. Les données de sortie x J d’un niveau de couches L j de l’auto-encodeur sont multipliées par le coefficient α j de la porte d’attention correspondante AG j , ainsi calculé, afin d’obtenir des données de sortie j de la première porte d’attention AG j , tel qu’illustré aux figures 3a et 3b. Cette étape est renouvelée jusqu’à une première porte d’attention AG 1 afin d’obtenir les données de sortie 1 de la première porte d’attention AG 1 . Comme visible à la figure 2, les données de sortie 1 de la première porte d’attention (AG 1 ) sont concaténées avec les données décodées provenant du deuxième niveau (L 2 ) afin d’obtenir, dans la dernière étape de la FIG. 1 , une matrice de données représentative d’une image de la distribution des propriétés électriques d’au moins un matériau du corps. Dans cet exemple, l’image est une matrice 32 x 32. La FIG. 4 présente une partie des résultats obtenus grâce à la mise en œuvre du procédé selon l’invention et comparés à d'autres architectures connues : 2D-CNN, Variational-Auto-Encoder (VAE), Generative DNN with U-Net, Dense ResNet, Conv ResNet, U-Net3++, Attention-Net. Les images sur la FIG. 4 correspondent à des tests réalisés sur un corps ayant une partie cylindrique qui peut comporter, de gauche à droite: - deux bâtons de section circulaire (de 1 cm et 2 cm de diamètre), - deux bâtons de section circulaire (de 1 cm de diamètre), - un bâton de section circulaire (de 2 cm de diamètre), - un bâton de section rectangulaire (de 2.5 × 1 cm de surface), et - deux bâtons de section circulaire (de 1 cm et 2 cm de diamètre). Il est à noter que, dans un souci de concision, seuls les résultats de certains tests ont été illustrés mais que les conclusions sont les mêmes pour tous les tests réalisés. Nous pouvons constater qu'avant U-Net3++, les architectures connues n'étaient pas suffisamment optimisées pour reconstruire la conductivité, principalement en raison du problème du gradient de fuite, puisqu'il n'y a pas de connexions sautées. L'exception serait l'architecture DNN+U-Net dans laquelle U-Net inclut ces connexions sautées mais l'extraction des caractéristiques est initiée avec un réseau de neurones profond, qui n'est pas optimisé pour les données d'image par rapport à un réseau de neurones convolutif. Pour U-Net3++, en regardant la deuxième image de gauche à droite, le réseau manquait d'identification des caractéristiques spatiales pour une reconstruction précise, ce qui a été amélioré avec le réseau classique Attention-Net, mais ce réseau de neurones a eu des problèmes avec la conductivité de l'eau en arrière-plan et peut être vu dans l'image N°5 de la FIG. 4 . Le réseau utilisé dans l’invention montre de très bons résultats de reconstruction avec une meilleure identification des bords spatiaux. Le tableau 1 montre que les mesures donnent de meilleurs résultats avec l’invention. Même si le R² a une valeur élevée pour l'architecture proposée, il peut conduire à de mauvaises interprétations. Les valeurs de RMSE et SER sont assez faibles, ce qui est le signe d'une bonne régression, et se situent dans la plage de la fonction de perte MSE, comme illustré à la FIG. 5 . Métriques R² RMSE SER Architectures (%) (μS/cm) (μS/cm) 2D-CNN 68.3 0.3001 0.3240 VAE 70.9 0.2769 0.2831 DNN + U-Net Dense ResNet Conv ResNet U-Net3++ Attention-Net AIA-Net 76.3 80.4 85.9 91.7 92.3 98.5 0.2547 0.2010 0.1832 0.1521 0.1308 0.1012 0.2342 0.1912 0.1447 0.1272 0.1187 0.1009 La FIG. 5 montre que la fonction de perte de l’erreur quadratique moyenne est très faible, ce qui indique une bonne correspondance entre les données d'entrée de tension et les données de sortie de conductivité. La perte de validation présente de faibles variations, ce qui indique que le modèle du réseau de neurones mis en œuvre est stable. En outre, l'absence de surajustement ou de sous-ajustement est expliquée par les valeurs proches entre les deux pertes. L’invention n’est pas limitée aux exemples qui viennent d’être décrits ; on peut notamment combiner entre elles des caractéristiques des exemples illustrés au sein de variantes non illustrées. D’autres variantes et améliorations peuvent être envisagées sans pour autant sortir du cadre de l’invention. En particulier, le procédé selon l’invention peut être mis en œuvre au moyen d’un réseau de neurones différent de celui décrit, comportant notamment des couches supplémentaires. L’invention peut être utilisée dans le domaine nucléaire, notamment comme méthode de prévention lors de présence de bulles dans le circuit primaire, et de cavitation dans les pompes des différents circuits. L’invention peut être appliquée au domaine médical, en particulier pour la supervision des poumons et la détection d’anomalies, en mettant notamment en œuvre une excitation simultanée. Le procédé selon l’invention peut trouver une application dans le domaine agroalimentaire, notamment pour le comptage de fruits et légumes, dans le domaine pétrolier, en particulier pour la détection d’encrassement dans les puits de production, ou encore dans le domaine pharmaceutique, par exemple pour la détection de corps étrangers lors de l’usinage de médicaments. Liste des références cités R. Guardo & al., "A Neural Network Approach To Image Reconstruction In Electrical Impedance Tomography", IEEE, 1991 A. Nejatali & al., "An iterative algorithm for electrical impedance imaging using neural networks", IEEE, 1998 J. Lampinen & al., "Application of Bayesian neural network in electrical impedance tomography", IEEE, 1999 Chao Wang & al., "RBF neural network image reconstruction for electrical impedance tomography", IEEE, 2004 Xiuyan Li & al., "An image reconstruction framework based on deep neural network for electrical impedance tomography", IEEE, 2017 Jin K. S. & al., "A Learning-Based Method for Solving Ill-Posed Nonlinear Inverse Problems: ASimulation Study of Lung EIT", SIAM J. Imaging sciences, 2019 Xiaoyan C. & al., "Deep Autoencoder Imaging Method for Electrical Impedance Tomography", IEEE, 2021 Akarsh P. & al., "Improved Training of Physics-Informed Neural Networks Using Energy-Based Priors: a Study on Electrical Impedance Tomography", ICLR, 2023 Hao Y. & al., "High-resolution conductivity reconstruction by electrical impedance tomography using structure-aware hybrid-fusion learning", Sciencedirect, 2024 Ozan Oktay & al., “AttentionU-Net : Learning Where to look for the Pancreas,” 1 st Conference on MIDL, 2018 Procédé de reconstruction de la distribution des propriétés électriques d’un corps comportant une partie cylindrique au sein de laquelle s’écoule un fluide, utilisant un réseau de neurones et des données de valeurs électriques du corps mesurées au préalable par tomographie par impédance électrique à l’aide d’électrodes agencées autour d’une périphérie de la partie cylindrique du corps, chaque électrode ayant été excitée par un potentiel de forme prédéfinie, le réseau de neurones comportant : un auto-encodeur comprenant plusieurs niveaux (L 1 ,…,L J ) comportant chacun au moins une couche de convolution et au moins une couche d’abandon, et un module d’attention comprenant des portes d’attention (AG 1 ,…, AG J-1 ), ayant chacune un signal d’attention (g j ) associé, le procédé comportant les étapes suivantes : i) les données de valeurs électriques sont traitées successivement par chaque niveau (L 1 ,…,L J ) de couches de convolution et de couches d’abandon de l’auto-encodeur du réseau de neurones de sorte à coder les données pour obtenir des données encodées (x 1 ,…, x J ) pour chaque niveau de couches (L 1 ,…,L J ), ii) afin de calculer le coefficient (α j ) d’une porte d’attention (AG j ), cette dernière est configurée pour additionner les données encodées (x 1 ,…, x J ) au signal d’attention (g j ) de la porte d’attention (AG j ) correspondante, fonction au moins de la concaténation des données de sortie j+1 ) de la porte d’attention de rang supérieur (AG j+1 ) et des données décodées provenant du niveau de rang supérieur (L j+2 ) à celui de la porte d’attention de rang supérieur (AG j+1 ), le résultat de cette addition étant transformé par au moins une fonction mathématique de sorte à obtenir le coefficient (α j ), iii) les données de sortie (x J ) d’un niveau de couches (L j ) de l’auto-encodeur sont multipliées par le coefficient (α j ) de la porte d’attention correspondante (AG j ) pour obtenir des données de sortie j ) de la première porte d’attention (AG j ), iv) l’étape précédente est renouvelée jusqu’à une première porte d’attention (AG 1 ) afin d’obtenir les données de sortie ( 1 ) de la première porte d’attention (AG 1 ), v) lesdites données de sortie ( 1 ) de la première porte d’attention (AG 1 ) sont concaténées avec les données décodées provenant du deuxième niveau (L 2 ) afin d’obtenir une matrice de données représentative d’une image de la distribution des propriétés électriques d’au moins un matériau du corps. Procédé selon la revendication 1, dans lequel un échantillonnage par valeur maximale est mis en œuvre d’un niveau de couches (L 1 ,…,L J ) à l’autre afin de compresser les données de valeurs électriques. Procédé selon l’une quelconque des revendications précédentes, dans lequel des connexions de saut sont réalisées entre les niveaux de couches (L 1 ,…,L J ) pour obtenir les données encodées (x 1 ,…, x J ) à l’entrée des portes d’attention. Procédé selon l’une quelconque des revendications précédentes, dans lequel, à l’étape ii), les données encodées (x 1 ,…, x J ) sont, selon leur dimension spatiale, convoluées ou transposées pour s’adapter à la dimension spatiale du signal d’attention considéré. Procédé selon l’une quelconque des revendications précédentes, dans lequel le coefficient (α j ) de la porte d’attention (AG j ) correspondante est calculé grâce aux équations suivantes : où est le quotient d'attention utilisé pour déterminer les coefficients d'attention αj, et sont respectivement des fonctions d'activation d’unité linéaire rectifiée, notée ReLU, et sigmoïde, , et sont des transformations linéaires destinées à réduire la consommation de paramètres et de ressources informatiques, combinées aux termes de biais , , and pour former un ensemble de paramètres . Procédé selon la revendication précédente, dans lequel, la transformation (ψ) produisant une matrice de poids qui est donnée à la fonction sigmoïde d’activation afin de limiter les valeurs des coefficients d'attention, un rééchantillonneur est utilisé pour étendre les dimensions de la matrice de poids. Utilisation du procédé selon l’une des revendications précédentes pour la reconstruction de la distribution des propriétés électriques d’un conduit d’installation nucléaire, un écoulement diphasique traversant notamment le conduit. Produit programme d’ordinateur comportant un support et, enregistrées sur ce support, des instructions lisibles par un processeur pour que, lorsqu’exécutées, elles permettent de mettre en œuvre le procédé de reconstruction selon l’une quelconque des revendications 1 à 6. Dispositif pour la mise en œuvre du procédé selon l’une quelconque des revendications 1 à 6, comportant : un système d’acquisition comportant au moins un réseau logique programmable, un module de génération de signaux analogiques et un module de mesure de signaux analogiques, une pluralité d’électrodes étant notamment connectées au système d’acquisition, un ordinateur configuré pour contrôler le système d’acquisition, un réseau de neurones intégré ou relié à l’ordinateur, le réseau de neurones comportant : un auto-encodeur comprenant plusieurs niveaux (L 1 ,…,L J ) comportant chacun au moins une couche de convolution et au moins une couche d’abandon, et un module d’attention comprenant des portes d’attention (AG 1 ,…, AG J-1 ), ayant chacune un signal d’attention (g j ) associé.