FR 2862803 A1 20050527 FR 0313754 A 20031124 Pour contrôler la qualité des cordons de soudure et des points de soudure d'étanchéité, on acquiert une image radiographique numérique des zones correspondantes d'un crayon à contrôler. Ensuite, on compare l'image acquise avec une image de référence stockée. Cette image de référence a été obtenue à partir d'un crayon de combustible réputé ne présenter aucun défaut. On soustrait donc l'image de référence à l'image acquise, puis on traite numériquement l'image en résultant afin de faire apparaître plus nettement les éventuels défauts, et on affecte une note globale pour le cordon de soudure, une note pour le point de soudure d'étanchéité et une note pour les inclusions de tungstène. Au-delà de valeurs seuils, le crayon est considéré comme non- conforme. Du fait des incertitudes de positionnement des crayons et des tolérances de fabrication, la position d'un crayon et ses dimensions sur l'image acquise sont différentes de celles du crayon de l'image de référence. Cela impose donc un traitement complexe de l'image issue de la soustraction pour pouvoir détecter les éventuels défauts de fabrication. En outre, on a constaté que la fiabilité de ce contrôle restait faible, malgré ce traitement de l'image issue de la soustraction, ce qui impose de rejeter, par sécurité, un grand nombre de crayons comme étant nonconforme. Un but de l'invention est de résoudre ce problème en fournissant un procédé plus fiable de contrôle non-destructif d'un élément pour réacteur nucléaire. A cet effet, l'invention a pour objet un procédé de contrôle nondestructif d'un élément pour réacteur nucléaire, procédé comprenant les étapes 25 de: a) acquisition d'une image radiographique et numérique d'au moins une zone de l'élément, b) création d'une image de référence par traitement numérique de l'image acquise, et c) comparaison de l'image acquise, éventuellement traitée, à l'image de référence pour détecter la présence de défaut(s). Selon des modes particuliers de réalisation, le procédé peut comprendre l'une ou plusieurs des caractéristiques suivantes prise(s) isolément ou selon toutes les combinaisons techniquement possibles: - l'élément est une partie d'un assemblage de combustible nucléaire, - l'élément est un crayon de combustible nucléaire s'étendant le long d'un axe longitudinal et comprenant une gaine fermée par des bouchons inférieur et supérieur et contenant un combustible nucléaire, - la zone comprend un cordon de soudure entre un des bouchons et la gaine, - la zone comprend un point de soudure obturant un canal traversant un bouchon, - l'étape b) comprend une sous-étape b1) d'ouverture ou de fermeture de l'image par adjonction d'un élément structurant, - l'élément structurant a une forme allongée le long de l'axe longitudinal du crayon, - l'élément structurant est un segment de p pixels, p étant un nombre entier non nul, - un défaut à détecter est une porosité ou une sous-pénétration du cordon de soudure, et la sous-étape b1) est une sous-étape d'ouverture de l'image par adjonction de l'élément structurant, - un défaut à détecter est une inclusion de tungstène, et la sous-étape b1) est une sous-étape de fermeture de l'image par adjonction de l'élément structurant, - l'étape b) comprend, préalablement à la sous-étape b1), une sous- étape b0) de lissage de l'image par un convolueur, - le convolueur est un carré de n pixels de côté, n étant un nombre entier non nul, - l'étape c) comprend une sous-étape cl) de calcul de la différence entre l'image acquise, éventuellement traitée, et l'image de référence, et de division de cette différence par l'image acquise, éventuellement traitée, ou par l'image de référence, - le procédé comprend, après la sous-étape cl), une sous-étape c2) de multiplication de l'image par un coefficient correspondant sensiblement à l'intensité lumineuse maximale du dispositif de prise de vue utilisé pour acquérir l'image radiographique, - le procédé comprend, après la sous-étape c1), une sous-étape de lissage de l'image grâce à un convolueur, - le convolueur est un carré de q pixels de côté, q étant un nombre entier, le convolueur est allongé transversalement par rapport à l'axe longitudinal, - le procédé comprend, après la sous-étape cl), une étape de binarisation de l'image, - le défaut à détecter étant un défaut d'étanchéité au droit du point de soudure, l'étape b) comprend une sous-étape b1) de projection de l'image le long de l'axe longitudinal et de reconstruction de l'image à partir de sa projection le long de cet axe, - l'étape b) comprend, préalablement à la sous-étape b1), une sous-étape b0) de lissage de l'image acquise par un convolueur, - le convolueur est un carré de t pixels de côté, t étant un nombre entier, - l'étape c) comprend une sous-étape c1) de soustraction de l'image de référence à l'image acquise, éventuellement traitée, - le procédé comprend, après la sous-étape c1), une sous-étape de binarisation de l'image, - le procédé comprend une étape d) de détection automatique et de détermination de caractéristiques d'une région de l'image issue de l'étape c) correspondant à un défaut, - une des caractéristiques est la position du défaut détecté dans l'image, - une des caractéristiques est représentative de la dimension du défaut détecté, - le procédé comprend une étape d) de détermination de l'épaisseur axiale minimale du point de soudure, - les étapes a) à d) sont effectuées pour plusieurs angles de vue, - le procédé comprend une étape e) de reconstruction des défauts détectés dans les images correspondant à différents angles de vues, - l'étape e) comprend une sous-étape el) de détermination des positions qu'un défaut détecté dans une première image correspondant à un premier angle de vue peut occuper dans une deuxième image correspondant à un deuxième angle de vue, une sous-étape e2) de comparaison des positions ainsi déterminées avec des positions de défaut(s) effectivement détecté(s) dans la deuxième image afin de déterminer si le défaut a été détecté dans la deuxième image, et, si le défaut a été détecté dans la deuxième image, une sous-étape e3) de calcul d'une dimension du défaut à partir des caractéristiques représentatives des dimensions du défaut déterminées dans la première et la deuxième images, - le procédé comprend une étape e) de sommation des caractéristiques représentatives de la dimension déterminées pour plusieurs angles de vue et de comparaison de la somme avec une valeur seuil afin de rendre un verdict sur la conformité de l'élément avec des critères prédéterminés de fabrication, - le procédé comprend une étape e) de calcul d'une moyenne des épaisseurs minimales déterminées pour plusieurs angles de vue et de comparaison avec une valeur seuil afin de rendre un verdict sur la conformité de l'élément avec des critères prédéterminés de fabrication. L'invention sera mieux comprise à la lecture de la description qui va 20 suivre, donnée uniquement à titre d'exemple, et faite en se référant aux dessins annexés, sur lesquels: - la figure 1 est une vue schématique d'une installation de contrôle non- destructif pour la mise en oeuvre d'un procédé selon l'invention, - la figure 2 est une vue partielle et en coupe longitudinale, illustrant 25 l'extrémité supérieure d'un crayon de combustible nucléaire avant remplissage de la gaine par un gaz inerte, la figure 3 est une vue partielle, en coupe longitudinale et agrandie, illustrant l'extrémité supérieure du bouchon supérieur du crayon de la figure 2, après fermeture étanche du canal de remplissage, - la figure 4 illustre une image acquise lors de la mise en oeuvre du procédé de contrôle non-destructif selon l'invention, - la figure 5 est un diagramme en trois dimensions illustrant les niveaux de gris de la zone de l'image de la figure 4 correspondant au cordon de soudure circulaire, - la figure 6 est une vue partielle agrandie de la figure 5 correspondant à une région présentant une porosité, - la figure 7 est une vue analogue à la figure 6 après lissage, - la figure 8 est une vue analogue à la figure 6 après érosion, - la figure 9 est un diagramme illustrant l'effet d'une érosion sur une courbe, - la figure 10 est une vue analogue à la figure 6 après dilatation, - la figure 11 est une vue analogue à la figure 9 illustrant l'effet d'une dilatation, - la figure 12 est une vue analogue à la figure 9 illustrant l'effet d'un top-hat , - la figure 13 est une vue analogue à la figure 6 après soustraction logarithmique, - la figure 14 est une vue analogue à la figure 6 après un nouveau lissage, - la figure 15 représente l'image de la zone sélectionnée après ces différents traitements, - la figure 16 est un organigramme illustrant des étapes mises en oeuvre pour la reconstruction des porosités détectées, - la figure 17 est une vue illustrant l'affichage sur l'écran de l'ordinateur de l'installation après reconstruction des porosités, - la figure 18 est une vue analogue à la figure 15, après la mise en oeuvre des différents traitements d'images en vue de détecter les sous- pénétrations dans le cordon de soudure circulaire, - la figure 19 est un diagramme en trois dimensions illustrant les niveaux de gris de la zone de l'image de la figure 4 correspondant au point de soudure fermant le canal de remplissage du bouchon supérieur, la figure 20 est une figure analogue à la figure 19 après projection et reconstruction, - la figure 21 est une vue analogue à la figure 19 après soustraction, - la figure 22 représente l'image de la zone correspondant au point de soudure après la mise en oeuvre de ces différents traitements, et - les figures 23 à 25 sont des vues illustrant l'affichage sur l'écran de l'ordinateur de l'installation de la figure 1 après la mise en oeuvre du procédé de contrôle non-destructif. La figure 1 illustre schématiquement une installation de contrôle nondestructif de crayons de combustible nucléaire 3 pour réacteur à eau pressurisée. Comme illustré plus particulièrement par la figure 2, un tel crayon 3 s'étend le long d'une direction longitudinale L et comprend une gaine tubulaire 5 qui renferme un combustible nucléaire, par exemple sous forme de pastilles 7 empilées. La gaine 5 est fermée à ses extrémités par un bouchon 9 supérieur et un bouchon 9 inférieur, qui n'est pas visible sur la figure 2 mais est disposé à gauche sur la figure 1. La structure générale des bouchons 9 est analogue et seule celle du bouchon 9 supérieur sera décrite par la suite. Le bouchon 9 supérieur comprend une partie inférieure 11 qui est engagée dans l'extrémité supérieure de la gaine 5, pratiquement sans jeu. Cette partie inférieure 11 est reliée par un épaulement 13 à la partie courante du bouchon 9. Un cordon de soudure circulaire 15 s'étend entre cet épaulement 13 et la tranche supérieure de la gaine 5. Ce cordon 15 assure la fixation étanche du bouchon 9 supérieur à la gaine 5. Le bouchon 9 comporte en outre une partie amincie 17 permettant l'engagement de mors de préhension. En outre, et à la différence du bouchon 9 inférieur, le bouchon 9 supérieur comporte un canal 19 de remplissage de la gaine 5 par un gaz inerte sous pression. Ce canal est parfois dénommé canal de pressurisation. Le canal 19 comporte plusieurs tronçons successifs dont les diamètres sont généralement décroissants depuis la partie 11 engagée dans la gaine 5 vers la partie d'extrémité extérieure 21 du bouchon 9. A son extrémité supérieure ou extrémité externe, le canal 19 comporte un tronçon terminal 23 ou queusot. Ce tronçon 23 débouche dans la surface 25 d'extrémité supérieure du bouchon 9 via une ouverture d'entrée circulaire dont le centre se trouve sur l'axe longitudinal L. Après avoir rempli la gaine 5 des pastilles 7, on vient la fermer par les bouchons 9 inférieur et supérieur que l'on soude à la gaine 5. On notera qu'un ressort (nonreprésenté) est alors classiquement interposé entre la pastille 7 supérieure et le bouchon 9 supérieur. Ensuite, on remplit par l'intermédiaire du canal 19 la gaine 5 d'un gaz inerte sous pression, tel que de l'hélium, puis on assure la fermeture étanche du 35 canal 19 dans son tronçon supérieur 23. Comme illustré par la figure 3, cette fermeture étanche est obtenue grâce à un point de soudure 27 obtenu par fusion du matériau du bouchon 9 au voisinage du tronçon supérieur 23. En revenant à la figure 1, l'installation 1 comprend principalement: - une source 29 de rayons X, - un dispositif 31 de conversion de rayons X en rayons lumineux, - un dispositif 33 d'acquisition d'images numériques, des moyens 35 de déplacement d'un crayon de combustible 3, et - un système informatique 37. La source 29, le dispositif de conversion 31 et le dispositif d'acquisition 33 sont disposés à l'intérieur d'une enceinte 39 de protection biologique qui reçoit l'extrémité du crayon 3 à examiner au travers d'une ouverture 41 susceptible d'être fermée, en l'absence de crayon 3, par un clapet 43. Sur la figure 1 c'est l'extrémité supérieure du crayon 3 qui est introduite dans l'enceinte 39 et qui va donc être contrôlée. La source 29 comprend un tube 45 de création de rayons X et un collimateur 47 pour orienter les rayons X radialement vers le bouchon 9 supérieur dont il faut contrôler l'étanchéité des soudures 15 et 27. La source 29 comprend en outre un obturateur 49 qui peut être commandé, tout comme le tube 45, par le système informatique 37. La source 29 comprend également un indicateur de qualité de l'image (IQI) 51. Classiquement un tel indicateur comprend une plaque percée de deux trous calibrés. La comparaison des dimensions réelles de ces trous avec celles qu'ils ont sur l'image acquise permet de s'assurer de la qualité de l'image acquise. Le bouchon 9 supérieur est engagé dans un bloc de compensation 53 qui a une forme sensiblement parallélépipédique et présente une cavité intérieure de forme sensiblement complémentaire de celle de l'extrémité supérieure du crayon 3 qu'elle reçoit. Ce bloc 53 est réalisé par exemple en alliage de zirconium. Les rayons X issus du tube 45 traversent le bloc de compensation 53 et l'extrémité supérieure du crayon 3 radialement par rapport à sa direction longitudinale L. Ensuite, les rayons X viennent frapper le dispositif 31 qui est par exemple un écran comprenant une pastille d'aluminium recouverte d'oxysulfure de gadolinium dopé au terbium. Les rayons X sont alors transformés en photons qui sont renvoyés, par l'intermédiaire d'un miroir 55, vers le dispositif 33 d'acquisition d'images. Ce dispositif 33 est par exemple une caméra CCD. La caméra 33 est reliée au système informatique 37 pour lui fournir les 5 images numériques acquises et être pilotées par ce dernier. L'image ainsi acquise LAI est représentée sur la figure 4. La lettre A signifie qu'il s'agit de l'image acquise et le chiffre 1, qu'il s'agit de la première prise de vue. La bande verticale sombre au centre correspond à un cache en plomb que l'on a placé dans le bloc de compensation 53 pour pallier la forte luminosité que devrait induire la partie amincie 17 du bouchon 9. Les moyens 35 de déplacement comprennent notamment des galets motorisés 57 pour déplacer longitudinalement le crayon 3 et une pince tournante 59 permettant de le faire tourner autour de son axe L, afin d'acquérir des images selon différents angles de vue. Les galets motorisés 57, qui sont pilotés par le système informatique 37, permettent de faire pénétrer l'extrémité d'un crayon 3 à examiner dans l'enceinte 39, de l'engager dans le bloc de compensation 53 et de l'extraire de l'enceinte 39 après mise en oeuvre du procédé de contrôle décrit ultérieurement. Les moyens 35 de déplacement comprennent en outre un capteur 61 de mesure de la rotation du crayon 3 raccordé au système informatique 37. Ce capteur 61 comprend par exemple une roue en appui sur le crayon 3 pour être entraînée en rotation par le crayon 3 lorsqu'il tourne. Le système informatique 37 comprend un ordinateur comprenant lui-même un ou plusieurs processeurs, des moyens de stockage, des moyens d'entrée/sortie et éventuellement des moyens d'affichage. Dans l'exemple représenté, il s'agit d'un micro-ordinateur qui comprend notamment un écran 63, une unité centrale 65 et un clavier 67. Un programme est stocké dans l'ordinateur 37 et permet, en association avec une librairie de traitement d'images la mise en oeuvre du 30 procédé de contrôle non-destructif. De telles librairies de traitement d'images sont couramment disponibles dans le commerce. Les étapes mises en oeuvre par le procédé de contrôle non-destructif sont les suivantes. L'ordinateur 37 va tout d'abord détecter et analyser les éventuelles 35 porosités présentes dans le cordon de soudure circulaire 15. De telles porosités doivent se traduirent par une moindre absorption des rayons X et donc par des régions plus claires. Pour cela, on sélectionne une zone IA1 C de l'image acquise IA1 qui correspond à la zone du crayon 3 contenant le cordon de soudure 15, comme illustré par la figure 4. La lettre C signifie que cette zone est relative au cordon de soudure 15. La figure 5 représente les différents niveaux de gris de la zone IA1 C. Ces niveaux de gris ont été portés selon l'axe Oz du repère Oxyz représenté sur cette figure. On a également matérialisé les axes Ox et Oy sur la figure 4. Comme cela apparaîtra ultérieurement, le cordon de soudure 15 présente une porosité 71 qui n'est pas décelable sur l'image acquise. La figure 6 est une vue agrandie de la partie de la figure 5 où la porosité 71 est située. Pour pouvoir détecter cette porosité 71, on va créer par traitement 15 numérique de la zone IA1 C, une image de référence IR1 C correspondant à un cordon 15 qui ne comporterait pas de défaut. Pour cela, on va soumettre la zone IA1 C à une opération de lissage par un convolueur g qui est un carré de n pixels de côté, n étant un nombre entier, de préférence compris entre 3 et 9. Le convolueur g est centré sur son pixel central. L'opération à laquelle la zone lA1 C est soumise correspond donc à la formule suivante: +00 +00 lxg(x,y)= E Eg(v,w)xl(x v,y w) V= co W= co Dans cette formule et celles qui suivent, I désigne l'image en cours de traitement. Cette opération va donc remplacer chaque pixel par la moyenne des pixels environnants contenus dans le convolueur centré sur le pixel considéré. Par exemple avec un convolueur g de 5 pixels de côté, chaque pixel serait remplacé par la moyenne du pixel considéré et des 24 pixels environnants. La zone IA1 C issue de l'étape de lissage est notée ILI C. La figure 7 est une vue analogue à la figure 6 à l'issue de l'étape de lissage. Ensuite, on va soumettre la zone IL1 C lissée à une ouverture par adjonction d'un élément structurant h qui est un segment, parallèle à l'axe Ox, et centré sur son pixel central. Ce segment a une longueur de p pixels, p étant un nombre entier, de préférence compris entre 5 et 21. L'ouverture par adjonction est une érosion par un élément structurant suivi d'une dilatation adjointe, c'est à dire d'une dilatation par le même élément structurant. L'ouverture par adjonction de l'élément structurant h peut donc être définie par la formule: (Ieh)O+h L'érosion correspond à la première partie de cette formule et consiste à passer l'élément structurant h sur la zone ILI C, comme un convolueur, et à affecter la valeur minimum au pixel central. Ainsi, l'érosion correspond à la formule (I e h) (x) = Min (I(v)) v E h (x) Plus de détails sur les opérations de morphologie mathématique telles que les ouvertures, les top-hat , ainsi que sur les opérations de traitement d'images telles que les lissages... peuvent être trouvés dans l'ouvrage de J. SERRA: Image Analysis and mathematical Morphology Ac. Press, Tome 1 (1982), Tome 2 (1988). La figure 8 correspond à la figure 6 après cette opération d'érosion. La zone traitée est alors dénommée IE1 C. La figure 9 illustre schématiquement et de manière générale l'effet d'une telle érosion produite sur les variations des niveaux de gris d'une image le long de l'axe Ox. La courbe 73 correspond à une image avant érosion et la courbe 75 à une image après érosion. L'érosion garde donc le fond de l'image en éliminant les pics de niveau de gris dont la largeur est inférieure à celle de l'élément structurant. La dilatation, qui est la seconde partie de l'opération d'ouverture, a pour but de ramener les niveaux de gris de l'image où ils se trouvaient avant l'érosion. La dilatation consiste à passer l'élément structurant h sur l'image, comme un convolueur, et à affecter la valeur maximum au pixel central, conformément à la formule: (I O h) (x) = Max (I(v)) v eg(x) Comme illustré par la figure 10, les niveaux de gris de l'image reviennent à ceux de l'image lissée, mais la porosité 71 a disparu. Sur la figure 11, la courbe 77 correspond à une image après dilatation. L'image obtenue à l'issue de cette opération d'ouverture par adjonction constitue donc une image de référence IR1 C puisque les défauts éventuels en ont été gommés. On notera que la position et les dimensions des différentes parties du crayon 3 sur cette image correspondent à celles de la zone IA1 C de l'image acquise puisqu'elle a été obtenue par traitement numérique de cette dernière. Ensuite, pour détecter et analyser les défauts éventuels, on va comparer l'image lissée ILI C à l'image de référence IR1 C. Pour cela, on va effectuer l'opération suivante: ILRIC IRICx K IR1C Cela signifie que l'on va soustraire au niveau de gris de chaque pixel de l'image ILI C, le niveau de gris du pixel correspondant de l'image de référence IR1 C, puis diviser par le niveau de gris du pixel de l'image de référence IR1 C, et enfin multiplier par le coefficient K et réaffecter le résultat ainsi obtenu au pixel considéré. La première partie de cette formule correspond à une soustraction dite logarithmique . En effet, l'image de référence IRIC n'est pas seulement soustraite à l'image lissée ILIC, mais cette différence est en outre divisée par un des termes de la soustraction. On pourra trouver plus de renseignements sur les soustractions 25 logarithmiques dans la publication: A model for logarithmic image processing by Michel Jourlin and Jean- Charles Pinoli Journal of Microscopy, Vol. 149, Pt 1, January 1988, pp. 21-35 (The Royal Microscopical Society). Une telle soustraction logarithmique permet de prendre en compte 30 l'absorption des rayons X par la matière. En effet, un défaut vu au travers de peu de matière aurait tendance à apparaître plus gros que ce même défaut vu au travers de plus de matière. Cette soustraction dite logarithmique permet donc de compenser cet effet. Le coefficient multiplicateur K est un paramètre représentatif de la caméra 33 et peut correspondre sensiblement à la dynamique ou résolution de l'image acquise par la caméra 33, c'est à dire à la valeur maximale de l'intensité lumineuse. Cette multiplication par le coefficient K permet de ramener les niveaux de gris à ceux de l'image IA1 C. La combinaison des opérations d'ouverture par adjonction de l'élément structurant h puis de soustraction peut être en fait apparentée à un top-hat ou chapeau haut de forme , bien que dans de tels top hat la soustraction ne soit pas logarithmique . Sur la figure 12 qui illustre de manière générale l'effet d'un tel top-hat , la courbe 78 correspond à la courbe 73 après ouverture et la courbe 79 correspond à la soustraction de la courbe 78 à la courbe originale 73. La figure 13 correspond à la figure 6 à l'issue du traitement de tophat appliqué par l'ordinateur 37. L'image est alors dénommée ISI C. L'opération de top-hat a donc fait ressortir la porosité 71. Ensuite, on va soumettre l'image ISI C à une opération de moyennage, c'est à dire à un lissage, grâce à un convolueur j qui est également un carré, mais de dimensions plus importantes que le convolueur g. Ainsi, ce carré a par exemple un carré de q pixels de côté, q étant un nombre entier, de préférence compris entre 5 et 21. Comme on le voit sur la figure 14 illustrant l'image à l'issue de cette seconde étape de lissage, la porosité 71 devient plus compacte et les parasites sont atténués. L'image obtenue est alors dénommée IMICP. La lettre P précise que cette image est dédiée à la recherche de porosités. L'image IMI CP va alors être soumise à une étape de binarisation. Une telle étape consiste à afficher en blanc les pixels de l'image dont le niveau de gris est au-dessus d'un seuil de binarisation déterminé, et afficher en noir ceux dont le niveau de gris est en dessous de ce seuil. Le seuil de binarisation, qui peut être modifié, correspond à la sensibilité souhaitée de l'analyse. L'image finale à l'issue de ces traitements apparaît donc comme illustré par la figure 15 et est désignée IFICP. La porosité 71 y apparaît sous forme d'une tâche blanche ou singularité. Ensuite, l'ordinateur 37 va détecter et analyser les défauts dans l'image finale IFICP. Pour cela, l'ordinateur 37 utilise l'outil de recherche de singularités ou blobs de la librairie de traitement d'images qui en outre détermine des caractéristiques des singularités de l'image IFICP, c'est à dire des porosités ainsi détectées. En particulier, l'ordinateur 37 calcule la position des porosités et leur diamètre. Ensuite, l'ordinateur 37 commande la rotation du crayon 3 d'un angle de 120 , acquiert une image IA2 puis procède à l'ensemble des opérations précédemment décrites sur sa zone IA2C, ce qui permet d'obtenir une deuxième vue IF2CP de la porosité 71. Ici, le chiffre 2 signifie qu'il s'agit des images relatives à la deuxième vue. Enfin, l'ordinateur 37 provoque à nouveau une rotation de 120 du crayon 3 autour de son axe longitudinal L et acquiert une image IA3 dans cette nouvelle position dont il traite la zone IA3C de la manière indiquée précédemment. Ensuite, l'ordinateur 37 procède à une reconstruction du ou des défauts détectés, en l'occurrence, la porosité 71. Les opérations alors mises en oeuvre par l'ordinateur 37 sont illustrées par l'organigramme de la figure 16. Dans une première étape 91, l'ordinateur 37 sélectionne une porosité détectée dans la vue IF1 CP. L'ordinateur 37 calcule ensuite, à partir des coordonnées de la porosité dans la vue IF1 CP, les coordonnées qu'elle pourrait occuper dans les vues IF2CP et IF3CP. C'est l'objet de l'étape 93. Ensuite, dans l'étape 95, l'ordinateur 37 vérifie si ces coordonnées calculées correspondent à celles d'une porosité effectivement détectée dans les vues IF2CP et IF3CP. Si ce n'est pas le cas, la porosité détectée seulement dans la vue IF1 CP est ignorée et l'étape 91 est effectuée avec une autre porosité de la vue IF1 CP. Si le résultat de l'étape de comparaison 95 est positif, la porosité a alors été détectée dans au moins deux vues et l'ordinateur 37 peut déterminer dans l'étape 97 ses coordonnées exactes et évaluer son diamètre réel. Si la porosité n'a été détectée que dans deux vues, le diamètre réel est calculé à partir, par exemple, d'une formule du type Am2 + Bm + c où m est la moyenne des diamètres déterminés dans les deux vues dans laquelle la porosité a été détectée. Si la porosité a été détectée dans trois vues, le diamètre réel estcalculé à partir d'une formule du type A' m + B' où m est le diamètre moyen des diamètres déterminés dans les trois vues IFI CP, IF2CP et IF3CP. Ainsi, l'ordinateur 37 connaît les positions exactes des porosités 5 détectées et a déterminé leurs diamètres réels. Il est alors en mesure de représenter en coupe transversale ces porosités. Cela est illustré par la figure 17 qui représente l'affichage sur l'écran 63 du cordon de soudure 15 en coupes transversales et où la seule porosité 71 10 détectée est affichée. Chacune des coupes transversales correspond à un angle de prise de vue. L'ordinateur 37 peut enfin comparer les diamètres des porosités détectées à une valeur maximale autorisée, ainsi que comparer la somme de ces diamètres à une valeur seuil. Si les valeurs mesurées et la somme sont inférieures à la valeur maximale et à la valeur seuil, le crayon 3 est déclaré conforme pour ce qui concerne les porosités du cordon de soudure 15. L'installation 1 décrite ci-dessus permet ensuite de détecter et d'analyser la présence de sous-pénétration(s) dans le cordon de soudure 15. On rappelle qu'une telle sous-pénétration correspond à une épaisseur radiale locale insuffisante du cordon 15 de sorte que l'étanchéité entre le bouchon 9 considéré et la gaine 5 peut être compromise. Une telle sous-pénétration, qui correspond à une absence de matière, apparaît sous forme d'une tâche plus claire sur les radiographies. Pour détecter de telles sous-pénétrations, l'ordinateur 37 traite les 25 images IA1 C, IA2C et IA3C. Les opérations de traitements sont les mêmes que celles décrites cidessus pour ce qui concerne la recherche des porosités, sauf que la seconde étape de lissage ou étape de moyennage est effectuée avec un convolueur k, qui est d'une part de dimensions plus importantes que dans le cas de la recherche des porosités et qui, d'autre part, est allongé le long de la direction Oy, c'est-à-dire transversalement à l'axe longitudinal L. Il peut ainsi s'agir d'un rectangle centré sur son pixel central. Ce rectangle peut avoir par exemple r pixels de côté selon la direction Ox, r étant un nombre entier compris entre 5 et 15, et s pixels de côté selon la direction Oy, s étant un nombre entier par exemple compris entre 10 et 20. L'utilisation d'un convolueur k allongé le long de la direction Oy permet de faire ressortir les défauts ayant une telle forme, ce qui est le cas des sous-pénétrations qui ont une certaine étendue angulaire, mais pas le cas des porosités qui s'apparentent à des bulles. De même, le seuil de binarisation, c'est-à-dire le seuil de niveau de gris à partir duquel les pixels sont affichés en blanc, est plus élevé que pour la recherche des porosités. Les images finales obtenues à l'issue de ces traitements sont référencées IF1 CS, IF2CS, IF3CS. La figure 18 illustre l'image IFI CS. Cette image est issue d'une succession d'images correspondant à chacune des étapes du traitement, à savoir les images IA1 C, IL1 C, IEI C, IDIC, ISIC, IM1 CS et IFCS. La lettre S indique que ces images sont utilisées pour la recherche des sous-pénétrations. On peut également déterminer grâce à la librairie de traitement d'images et afficher les coordonnées et les surfaces des singularités 99 mises en évidence, en l'occurrence les sous-pénétrations. Ces sous-pénétrations 99 sont au nombre de deux sur l'image IFI CS. Ces mêmes étapes sont mises en oeuvre à partir des images IF2CS et 20 IF3CS. Ensuite, l'ordinateur 37 peut calculer une note relative aux souspénétrations du cordon de soudure 15. Cette note peut être la somme de toutes les surfaces des sous-pénétrations 99 détectées, éventuellement pondérées. La comparaison de cette note de sous-pénétration avec une note seuil permet de rendre un verdict sur la conformité du crayon 3 analysé pour ce qui concerne les sous-pénétrations du cordon de soudure 15. L'installation 1 permet ensuite de vérifier l'étanchéité au niveau du point de soudure 27 du canal de remplissage 19 de la gaine 5. Pour cela on sélectionne une zone IA1Q de l'image acquise IA1, comme illustré par la figure 4. La lettre Q signifie que cette zone est relative à la vérification de l'étanchéité au droit du queusot 23. La figure 19 illustre les niveaux de gris dans cette zone IA1 Q. Ensuite on soumet cette image IA1 Q à une opération de lissage par un convolueur I qui est par exemple un carré de t pixels de côté centré sur son pixel 35 central, t étant un nombre entier par exemple compris entre 3 et 9. La taille de ce convolueur est réduite pour modifier le moins possible l'image IA1 Q. L'image à l'issue de ce traitement est dénommée ILI Q. Ensuite on soumet l'image ILI Q à une étape de projection selon l'axe Ox et de reconstruction. Plus précisément, cette étape de projection et de reconstruction consiste à affecter au pixel de chaque colonne s'étendant le long de l'axe Oy, la valeur moyenne des niveaux de gris des pixels de cette colonne. Ce traitement permet de privilégier le fond de l'image ILIQ est donc de gommer les éventuelles sous-pénétrations longitudinales locales du point de soudure. L'image ainsi obtenue est une image de référence IRIQ dont les niveaux de gris varient uniquement le long de la direction Ox, comme illustré par la figure 20. Ensuite, on soustrait l'image de référence IR1 Q à l'image lissée ILI Q. II s'agit ici d'une soustraction simple ILI Q-IR1 Q et non pas d'une soustraction logarithmique. La figure 21 illustre les niveaux de gris à l'issue de cette opération. Ensuite, on binarise l'image IS1 Q issue de la soustraction afin d'obtenir une image finale IFIQ où tous les pixels ayant un niveau de gris inférieur à une valeur seuil se voient affecter la valeur zéro, et tous les autres se voient affecter la valeur maximale en niveau de gris. La figure 22 illustre l'image IF1 Q. On peut y distinguer à gauche la région 101 du canal 19 sous le point de soudure 27 et à droite la cavité 103 formée au-dessus du point de soudure 27. On notera qu'en cas de soudure TIG, le point soudure 27 ne forme pas de cavité mais présente une surface extérieure convexe. Ensuite, l'ordinateur 37 va détecter, grâce à l'outil de recherche de singularités de la librairie de traitement d'images, les zones 101 ou 103. Si, aucune zone ou une seule zone traversant horizontalement l'image IA1 Q est détectée, le queusot 23 est déclaré non-soudé et le crayon 3 nonconforme. Ce dernier cas signifie en effet que les zones 101 ou 103 se rejoignent. Si deux zones 101 et 103 sont détectées, l'ordinateur 37 calcule, grâce à la librairie de traitement d'images, la distance entre ces deux zones 101 et 103 qui apparaissent en blanc sur l'image IFIC. Pour cela, la librairie trace par exemple le bord 105 de la cavité 103 sur l'image IF1 C puis va construire le bord 107 de la zone 101. Pour cela l'ordinateur 37 parcourt, pour chaque ligne, l'image IFIQ vers la droite en partant du bord gauche de l'image IFIQ jusqu'à trouver le bord 105. Ensuite, l'ordinateur 37 parcourt la ligne considérée en sens inverse jusqu'à détecter le premier pixel blanc. S'il détecte un tel pixel blanc, celui-ci est le point du bord 107 de la zone 101 pour la ligne considérée. S'il ne détecte pas de pixel blanc, ce qui correspond en fait au cas où seule une zone 103 est présente à droite de l'image IFIQ, l'ordinateur 37 considère que le bord 107 est le bord gauche de l'image IFIQ. Ainsi, en procédant ligne par ligne, l'ordinateur 37 va déterminer les coordonnées des points des bords 105 et 107 des deux zones claires 101 et 103. Ensuite, à partir de l'ensemble de ces coordonnées, l'ordinateur 37 va calculer la distance minimale entre ces bords 105 et 107, distance qui correspond à la pénétration axiale du point de soudure 27 dans le canal 19. L'ordinateur 37 va ensuite procéder ainsi pour les vues IA2 et IA3 et fournir, à partir des zones IA2Q et IA3Q des images finales IF2Q et IF3Q, puis déterminer les pénétrations axiales du point de soudure 27 dans le canal 19 pour chacune de ces vues. L'ordinateur 37 calcule enfin la pénétration axiale réelle du point de soudure 27 comme étant la moyenne des pénétrations axiales déterminées à partir des trois vues IF1Q, IF2Q et IF3Q. Il compare enfin cette valeur avec une valeur seuil, en-deçà de laquelle il considère que le point de soudure 27 ne pénètre pas suffisamment dans le canal 19 pour assurer l'étanchéité du crayon 3. On notera également que les images IA1 Q, IA2Q et IA3Q peuvent être soumises à un traitement pour la recherche de porosités, ce traitement étant analogue à celui décrit précédemment pour le cordon de soudure 15. Enfin, l'installation 1 permet également, dans le cas où le cordon 15 et/ou le point 17 est/sont réalisé(s) par soudure TIG de déterminer la présence d'inclusion(s) de tungstène. En effet, ces inclusions apparaissent sous forme de zones plus sombres sur les radiographies numériques acquises par la caméra 33. On peut alors, grâce à des traitements analogues à ceux décrits pour les porosités ou les sous-pénétrations, détecter et analyser de tels défauts et rendre un verdict à partir de critères de conformité prédéterminés. On utilise alors à la place de l'opération d'ouverture par adjonction, une opération de fermeture par un élément structurant h. Cette opération est définie par la formule Une telle fermeture permet d'éliminer les pics négatifs et ainsi de former des images de référence satisfaisantes. L'élément structurant h est alors de préférence analogue à celui décrit précédemment. On notera que les défauts recherchés apparaissant sous forme de pics négatifs dans les niveaux de gris, la binarisation est inversée par rapport à celle décrite précédemment. En variante, c'est l'opération de soustraction qui peut être inversée, c'est-à-dire que l'on soustrait l'image lissée à l'image érodée, le résultat étant divisé par l'image érodée et multipliée par le coefficient K. II est alors possible de visualiser sur l'écran 63 l'ensemble des images acquises et traitées pour chaque angle de vue. La figure 23 illustre ainsi l'affichage de l'écran 63 pour le premier angle de vue. La succession supérieure d'images correspond aux images relatives à la recherche de porosités à savoir IA1 C, ILI C, IEIC, IDIC, ISIC, IMICP et IFICP. Les images IS1C et IMICP sont bien présentes entre les images IDIC et IFICP mais ne sont pas visibles du fait de leurs faibles niveaux de gris. La succession intermédiaire d'images correspond à une zone de l'image LAI dans lequel les orifices de l'indicateur de qualité de l'image (IQI) 51 sont visibles. Ces images ont été soumises exactement aux mêmes traitements que l'image IA1 C, y compris les opérations de détermination des positions et des diamètres des singularités. La succession inférieure d'images inférieures correspond à la succession d'images relatives au point de soudure 27, à savoir IA1 Q, IL1 Q, IR1 Q, IS1 Q et IF1 Q. On notera que dans l'exemple de la figure 23, le queusot 23 n'est pas soudé. La partie inférieure de l'écran 63 peut fournir des informations relatives à la position et au diamètre des porosités 71 détectées, à la position et au diamètre des taches 109 détectées dans la zone correspondant à l'indicateur de qualité d'image (IQI) 51 ainsi que des informations relatives à la pénétration axiale du point de soudure 27 pour l'angle de vue considéré. Les dimensions des taches détectées pour les images correspondant à l'IQI permettent de vérifier la fiabilité des mesures effectuées. L'écran 63 peut également afficher dans sa partie supérieure la succession d'images correspondant à la recherche de sous-pénétrations. Cela est illustré par la figure 24. Enfin, et comme illustré par la figure 25, l'écran peut également afficher simultanément: - dans une partie supérieure les images IA1C, IA2C, IA3C, IFICP, IF2CP, IF3CP, IFI CS, IF2CS, IF3CS, relatives au contrôle du cordon 15, - dans une région intermédiaire, les images initiales et finales pour les trois vues concernant l'indicateur de qualité IQI, - dans une région inférieure, les images IA1 Q, IA2Q, IA3Q, IF1 Q, IF2Q et IF3Q, relatives au contrôle du point de soudure 27. Ainsi, l'installation 1 a permis de déterminer automatiquement si le bouchon supérieur 9 du crayon 3 présentait des porosités ou des souspénétrations pour le cordon 15 ou un défaut d'étanchéité au droit du point de soudure 27, et éventuellement des inclusions de tungstène, justifiant le rejet du crayon 3. Ensuite, on procède de manière analogue pour le bouchon 9 inférieur. On a pu constater que le procédé de contrôle non-destructif décrit est particulièrement fiable et d'exécution rapide, puisqu'il permet de rendre un verdict pour chaque bouchon 9 en moins de 26 secondes. Le nombre de crayons 3 indûment rejetés est également réduit. En outre, le procédé décrit permet de détecter et caractériser les défauts séparément et d'en reconstruire certains, c'est à dire de déterminer leurs positions exactes et certaines de leurs caractéristiques. Cette fiabilité et cette rapidité d'exécution est due à la réalisation des images de référence, notamment IRIC, IR2C, IR3C, IRIQ, IR2Q et IR3Q à partir des images acquises IA1 C, IA2C, IA3C, IA1 Q, IA2Q, et IA3Q. En effet, cela permet de garantir une correspondance exacte entre les contours sur ces images, ce qui simplifie les traitements nécessaires pour isoler ensuite les défauts recherchés. De manière plus générale, d'autres traitements numériques des images acquises peuvent être effectués pour former des images de référence. Pour le cordon de soudure 25, on utilisera de préférence une ouverture par adjonction qui, combinée avec l'opération suivante de soustraction, forme un top-hat . L'utilisation d'un élément structurant h sous forme d'un segment parallèle à l'axe L s'avère particulièrement avantageuse. En effet, un tel élément structurant permet de limiter les aberrations introduites dans les images du fait de la présence du canal 19. Plus généralement, l'élément structurant h peut ne pas être un segment mais avoir une forme allongée le long de l'axe longitudinal L. Il peut ainsi s'agir d'un rectangle. De manière plus générale, les principes ci-dessus, à savoir la création d'une image de référence par traitement d'une image acquise de l'élément à contrôler, peuvent être mis en oeuvre pour contrôler d'autres éléments de réacteurs nucléaires que des crayons de combustible nucléaire, en particulier leurs soudures. 22 REVENDICATIONS 1. Procédé de contrôle non-destructif d'un élément (3) pour réacteur nucléaire, procédé comprenant les étapes de: a) acquisition d'une image radiographique et numérique (IA1C, IA2C, IA3C, IA1 Q, IA2Q, IA3Q) d'au moins une zone (15, 27) de l'élément (3), b) création d'une image de référence (IR1 C, IR2C, IR3C, IRIQ, IR2Q, IR3Q) par traitement numérique de l'image acquise, et c) comparaison de l'image acquise (IA1 C, IA2C, IA3C, IA1 Q, IA2Q, IA3Q), éventuellement traitée, à l'image de référence pour détecter la présence de défaut(s). 2. Procédé selon la revendication 1, caractérisé en ce que l'élément (3) est une partie d'un assemblage de combustible nucléaire. 3. Procédé selon la revendication 2, caractérisé en ce que l'élément est un crayon de combustible nucléaire (3) s'étendant le long d'un axe longitudinal (L) et comprenant une gaine (5) fermée par des bouchons (9) inférieur et supérieur et contenant un combustible nucléaire (7). 4. Procédé selon la revendication 3, caractérisé en ce que la zone comprend un cordon de soudure (15) entre un des bouchons (9) et la gaine (5). 5. Procédé selon la revendication 3, caractérisé en ce que la zone comprend un point de soudure (27) obturant un canal (19) traversant un bouchon (9). 6. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'étape b) comprend une sous-étape b1) d'ouverture ou de fermeture de l'image par adjonction d'un élément structurant (h). 7. Procédé selon les revendications 3, 4 et 6 prises ensembles, caractérisé en ce que l'élément structurant (h) a une forme allongée le long de l'axe longitudinal (L) du crayon (3). 8. Procédé selon la revendication 7, caractérisé en ce que l'élément structurant (h) est un segment de p pixels, p étant un nombre entier non nul. 9. Procédé selon la revendication 7 ou 8, caractérisé en ce qu'un défaut à détecter est une porosité (71) ou une sous-pénétration (99) du cordon de soudure (15), et en ce que la sous-étape b1) est une sous-étape d'ouverture de l'image par adjonction de l'élément structurant (h). 10. Procédé selon la revendication 7 ou 8, caractérisé en ce qu'un défaut à détecter est une inclusion de tungstène, et en ce que la sous-étape b1) est une sous-étape de fermeture de l'image par adjonction de l'élément structurant (h). 11. Procédé selon l'une quelconque des revendications 6 à 10, caractérisé en ce que l'étape b) comprend, préalablement à la sous-étape b1), une sous-étape b0) de lissage de l'image par un convolueur (g). 12. Procédé selon la revendication 11, caractérisé en ce que le convolueur (g) est un carré de n pixels de côté, n étant un nombre entier non nul. 13. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'étape c) comprend une sous-étape cl) de calcul de la différence entre l'image acquise, éventuellement traitée, et l'image de référence, et de division de cette différence par l'image acquise, éventuellement traitée, ou par l'image de référence. 14. Procédé selon la revendication 13, caractérisé en ce qu'il comprend, après la sous-étape cl), une sous-étape c2) de multiplication de l'image par un coefficient correspondant sensiblement à l'intensité lumineuse maximale du dispositif de prise de vue (33) utilisé pour acquérir l'image radiographique. 15. Procédé selon la revendication 13 ou 14, caractérisé en ce qu'il comprend, après la sous-étape cl), une sous-étape de lissage de l'image grâce à un convolueur (j, k). 16. Procédé selon la revendication 15, caractérisé en ce que le convolueur est un carré de q pixels de côté, q étant un nombre entier. 17. Procédé selon la revendication 15, caractérisé en ce que le convolueur (k) est allongé transversalement par rapport à l'axe longitudinal (L). 18. Procédé selon l'une des revendications 13 à 17, caractérisé en ce qu'il comprend, après la sous-étape cl), une étape de binarisation de l'image. 19. Procédé selon la revendication 5, caractérisé en ce que le défaut à détecter étant un défaut d'étanchéité au droit du point de soudure (27), l'étape b) comprend une sous-étape b1) de projection de l'image le long de l'axe longitudinal (L) et de reconstruction de l'image à partir de sa projection le long de cet axe (L). 20. Procédé selon la revendication 19, caractérisé en ce que l'étape b) 5 comprend, préalablement à la sous-étape b1), une sous-étape b0) de lissage de l'image acquise par un convolueur (I). 21. Procédé selon la revendication 20, caractérisé en ce que le convolueur (I) est un carré de t pixels de côté, t étant un nombre entier. 22. Procédé selon l'une des revendications 19 à 21, caractérisé en ce que l'étape c) comprend une sous-étape cl) de soustraction de l'image de référence à l'image acquise, éventuellement traitée. 23. Procédé selon la revendication 22, caractérisé en ce qu'il comprend, après la sous-étape c1), une sous-étape de binarisation de l'image. 24. Procédé selon l'une des revendications précédentes, caractérisé en ce que le procédé comprend une étape d) de détection automatique et de détermination de caractéristiques d'une région (71; 99) de l'image issue de l'étape c) correspondant à un défaut. 25. Procédé selon la revendication 24, caractérisé en ce qu'une des caractéristiques est la position du défaut détecté dans l'image. 26. Procédé selon la revendication 24 ou 25, caractérisé en ce qu'une des caractéristiques est représentative de la dimension du défaut détecté. 27. Procédé selon l'une des revendications 19 à 23, caractérisé en ce qu'il comprend une étape d) de détermination de l'épaisseur axiale minimale du point de soudure (27). 28. Procédé selon l'une des revendications 24 à 27, caractérisé en ce que les étapes a) à d) sont effectuées pour plusieurs angles de vue. 29. Procédé selon les revendications 25, 26 et 28 prises ensemble, caractérisé en ce qu'il comprend une étape e) de reconstruction des défauts détectés dans les images correspondant à différents angles de vues. 30. Procédé selon la revendication 29, caractérisé en ce que l'étape e) comprend une sous-étape el) de détermination des positions qu'un défaut détecté dans une première image correspondant à un premier angle de vue peut occuper dans une deuxième image correspondant à un deuxième angle de vue, une sous-étape e2) de comparaison des positions ainsi déterminées avec des positions de défaut(s) effectivement détecté(s) dans la deuxième image afin de déterminer si le défaut a été détecté dans la deuxième image, et, si le défaut a été détecté dans la deuxième image, une sousétape e3) de calcul d'une dimension du défaut à partir des caractéristiques représentatives des dimensions du défaut déterminées dans la première et la deuxième images. 31. Procédé selon les revendications 26 et 28 prises ensemble, caractérisé en ce qu'il comprend une étape e) de sommation des caractéristiques représentatives de la dimension déterminées pour plusieurs angles de vue et de comparaison de la somme avec une valeur seuil afin de rendre un verdict sur la conformité de l'élément (3) avec des critères prédéterminés de fabrication. 32. Procédé selon les revendications 27 et 28 prises ensembles, caractérisé en ce qu'il comprend une étape e) de calcul d'une moyenne des épaisseurs minimales déterminées pour plusieurs angles de vue et de comparaison avec une valeur seuil afin de rendre un verdict sur la conformité de l'élément (3) avec des critères prédéterminés de fabrication.