FR 3108761 A1 20211001 FR 2003093 A 20200328 Un procédé de squelettisation des torons d’une pièce en matériau composite dans une image 3D (V) comprenant une phase détection de centres de section orientés (C X , C Y ) des torons et une liaison (E L ) de chaque centre (C X , C Y ) au centre (C X , C Y ) de même orientation le plus proche, ladite détection étant mise en œuvre au moyen d’au moins un volume de référence comportant un motif de toron centré orienté, ladite phase détection comprenant pour chaque volume de référence : une étape de détermination de portions de l’image (V), une étape de calcul d’un score de corrélation entre le volume de référence et chaque portion, associé au voxel central de chaque portion, de manière à obtenir un score de corrélation pour chaque voxel de l’image (V) et une étape de détermination des centres (C X , C Y ) des torons correspondant aux voxels dont le score de corrélation correspond à un maximum local. Figure de l’abrégé : Figure 9 Procédé et système de squelettisation des torons d’une pièce en matériau composite dans une image 3D La présente invention concerne le domaine de la détermination de la structure interne d’une pièce en matériau composite, et vise plus particulièrement un procédé et un système de squelettisation des torons d’une pièce en matériau composite dans une image 3D. On connaît dans l’art antérieur une aube en matériau composite 1, dont une portion de la paroi est illustrée sur la figure 1, comprenant des fibres de renfort F, par exemple en carbone, noyées dans une matrice G, telle qu’une matrice organique à base de résine polymère, servant de liant et conférant à l’aube 1 sa forme extérieure. De manière connue, les fibres de renfort F sont organisées en faisceaux de fibres, connus sous le terme de « torons T X , T Y », comprenant de l’ordre d’une centaine à un millier de fibres de renfort F tressées ou alignées les unes avec les autres. Les torons T X , T Y sont quant à eux agencés entre eux par tissage et/ou empilement, tel que dans l’exemple de la figure 1, un empilement de torons longitudinaux T X tissés avec des torons transversaux T Y selon un repère orthogonal (X, Y, Z) défini par la direction de l’empilement Z, la direction longitudinale X et la direction transversale Y. Pour vérifier les propriétés mécaniques d’une aube en matériau composite 1 et notamment sa résistance à un endommagement, une analyse basée sur la forme extérieure de l’aube 1 n’est pas suffisante. Il est en effet nécessaire d’analyser la position et la forme des torons T X , T Y au sein de l’aube 1 afin de détecter tout défaut interne. En particulier, il est souhaité de détecter des torons T X , T Y qui sont rompus ou endommagés ou encore des irrégularités dans leur tissage. Pour cela, comme illustré sur la figure 2, il est connu d’acquérir par tomographie une image 3D en niveaux de gris N de l’aube 1 à partir d’un grand nombre de radiographies 2D obtenues en bombardant l’aube 1 de rayons X suivant différents angles. Une telle image 3D, dite « volume tomographique V » par la suite, est divisée en voxels vox s’étendant selon un repère orthogonal i, j, k. Chaque voxel vox comprend un niveau de gris N donné, correspondant au coefficient d’absorption de rayons X moyen de la zone d’aube 1 illustrée sur le voxel vox. Les voxels vox de niveaux de gris N sombres correspondent ainsi à une zone d’aube 1 comprenant majoritairement un matériau de faible coefficient d’absorption des rayons X N2, à savoir de la matrice G. Les voxels vox de niveaux de gris N clairs correspondent quant à eux à une zone d’aube 1 comprenant majoritairement un matériau de grand coefficient d’absorption N1, à savoir des torons T X , T Y . En pratique, une aube en matériau composite 1 utilisée dans l’industrie aéronautique possède une densité de torons T X , T Y comprise entre 80% et 90%, c’est-à-dire que les torons T X , T Y sont agencés de manière très compactée au sein de la matrice G, rendant l’identification des torons T X , T Y difficile dans le volume tomographique V. En outre, la plage de niveaux de gris diffère d’une radiographie 2D à l’autre en fonction des conditions d’acquisition. Le volume tomographique V est en outre imparfait, comportant des artefacts, rendant l’identification des torons T X , T Y d’autant plus complexe. En référence à la figure 3, pour analyser la forme et la position des torons T X , T Y dans le volume tomographique V, il est connu un procédé de squelettisation des torons T X , T Y par détection E D -AA puis liaison E L -AA de centres de section orientés C X , C Y de chaque toron T X , T Y dans le volume tomographique V. Dans l’exemple de la figure 4, des centres longitudinaux C X et transversaux C Y de torons respectivement longitudinaux T X et transversaux T Y sont ainsi déterminés par la détection E D -AA détaillée ci-après puis les centres C X , C Y associés à un même toron T X , T Y sont reliés entre eux. Les courbes obtenues sont désignées « squelettes SqX, SqY des torons T X , T Y » et permettent d’accéder à la position des torons T X , T Y dans le volume tomographique V. Plus précisément, en référence à la figure 5, la détection E D -AA des centres C X , C Y des torons T X , T Y est réalisée de manière connue sur des tranches W X , W Y du volume tomographique V normales à l’axe des torons T X , T Y . Dans cet exemple, les centres longitudinaux C X et transversaux C Y sont ainsi déterminés sur des tranches W X , W Y respectivement normales à l’axe longitudinal X et l’axe transversal Y, dites « tranches longitudinales W X » et « tranches transversales W Y ». De manière connue, un opérateur annote manuellement les centres C X , C Y sur environ 1% des tranches W X , W Y par interprétation visuelle des niveaux de gris N. Pour chaque toron, les centres C X , C Y des autres tranches W X , W Y sont ensuite déterminés par interpolation à partir des centres C X , C Y déterminés manuellement. Un tel procédé de squelettisation des torons T X , T Y d’une aube 1 est laborieux et coûteux en temps et en ressources. En particulier, l’annotation manuelle nécessite environ cinq jours de travail pour déterminer le squelette Sq X , Sq Y d’environ 800 torons T X , T Y . L’annotation manuelle est en outre soumise à une éventuelle erreur et une interprétation propre à l’opérateur, ce qui affecte également la précision des centres C X , C Y interpolés. En outre, un tel procédé de squelettisation ne permet pas de détecter des défauts parmi les centres interpolés. Pour supprimer ces inconvénients, il est connu par l’article intitulé « A 3D image analysis method for fibrous microstructures : Discretization and fibers tracking » publié en 2012, de binariser le volume tomographique en affectant un premier label aux zones de niveaux de gris claires et un deuxième label aux zones de niveau de gris sombres. La détection des centres des torons est ensuite réalisée en déterminant les centres des zones de premier label, attribuées aux torons. Un tel procédé ne permet toutefois pas une squelettisation exacte et précise des torons T X , T Y d’une aube 1 de densité de torons T X , T Y importante, la binarisation conduisant alors à des zones de premier label étendues de plusieurs torons contigus indistinguables les uns des autres. L’article « Towards the 3D in-situ characterization of deformation micro-mechanisms within a compressed bundle of fibres », publié en 2011, enseigne également de binariser le volume tomographique puis de détecter, sur plusieurs coupes, les centres des zones de label attribué aux torons. Il est également connu par l’article « FiberScout : an interactive tool for exploring and analyzing fiber reinforced polymers » publié en 2014 de déterminer, directement sur le volume tomographique, le squelette des torons à partir des positions de leurs extrémités et de leur orientation locale en plusieurs points. L’orientation locale d’un toron en un point donné est obtenue en déterminant, pour un voxel donné, le voxel voisin de niveau de gris le plus proche du sien. Un tel procédé ne permet pas non plus une squelettisation exacte et précise pour une aube 1 de densité de torons T X , T Y importante, la détermination de l’orientation locale des torons pouvant être faussée par la proximité des torons voisins, en particulier contigus, et par les artefacts du volume tomographique. L’article intitulé « Composites : Part A, Fiber orientation measurements in composite materials » publié en 2006 présente également un procédé de squelettisation à partir de l’orientation des torons, en déterminant une tranche du volume tomographique de normale ni parallèle ni orthogonale à l’axe des torons, de sorte que les sections des torons sur la tranche soient sensiblement elliptiques. L’orientation des torons est alors déterminée à partir de la forme elliptique des sections. Ce procédé présente toutefois les mêmes inconvénients que précédemment. De manière incidente, on connaît par l’article « Analysis of woven fabric structure using X-ray CT images » publié en 2005 un procédé de squelettisation des fils d’un tissu dans un volume tomographique, afin d’analyser son tissage. Dans un premier temps, des tranches de normale orthogonale à l’axe des fils du tissu sont déterminées dans le volume tomographique et binarisées tel que décrit précédemment. Les centres des fils de tissu sont ensuite annotés manuellement sur une première tranche. Sur la tranche voisine suivante, les centres des fils sont déterminés par corrélation avec les centres déterminés sur la première tranche, en recherchant chaque centre d’un fil dans une région de calcul restreinte autour du centre du même fil de la première tranche. La corrélation est réalisée successivement sur chacune des tranches suivantes, de proche en proche. Un tel procédé est peu robuste, car si un centre est mal déterminé sur une tranche, l’erreur risque d’être propagée aux tranches suivantes. En outre, un tel procédé nécessite une binarisation non adaptée pour une aube 1 de densité de torons T X , T Y élevée, qui plus est dans laquelle les fibres F sont enrobées dans une matrice G contrairement aux fils du tissu. L’invention vise ainsi un procédé et un système de squelettisation des torons d’une aube en matériau composite dans une image 3D, et plus généralement d’une pièce en matériau composite, qui soit précis, fiable et rapide à mettre en œuvre. PRESENTATION DE L’INVENTION L’invention concerne un procédé de squelettisation des torons d’une pièce en matériau composite dans une image 3D de ladite pièce, ladite image comprenant une pluralité de voxels s’étendant dans un repère orthogonal et comprenant chacun un niveau de gris, ledit procédé comprenant : une phase de détection d’une pluralité de centres de section orientés des torons de la pièce dans l’image et une phase de liaison de chaque centre au centre de même orientation le plus proche, afin de relier les centres associés à un même toron, de manière à obtenir le squelette des torons de ladite pièce. L’invention est remarquable en ce que la phase de détection est mise en œuvre au moyen d’au moins un volume de référence comprenant une pluralité de voxels, lesdits voxels s’étendant dans le repère orthogonal et comprenant chacun un niveau de gris de manière à ce que le volume de référence comporte un motif de toron centré s’étendant suivant une direction donnée, ladite phase de détection comprenant, pour chaque volume de référence : une étape de détermination de portions de l’image comportant un voxel central et de dimension égale à celle du volume de référence, une étape de calcul d’un score de corrélation entre le volume de référence et chaque portion de l’image par comparaison du niveau de gris de leurs voxels respectifs, ledit score de corrélation étant associé au voxel central de chaque portion, de manière à obtenir un score de corrélation pour chaque voxel de l’image formant une matrice de scores de corrélation, et une étape de détermination des centres de section orientés des torons de la pièce, lesdits centres correspondant aux voxels de l’image dont le score de corrélation correspond à un maximum local de la matrice de scores de corrélation, l’orientation desdits centres étant définie par la direction du motif de toron du volume de référence. On précise ici que par motif de toron, on désigne un motif d’une tranche transversale de toron inscrite dans le volume de référence. Autrement dit, le motif de toron représente un toron tronqué dans le sens de la longueur et dont la section est entièrement représentée. Grâce à l’invention, la phase de détection des centres de section des torons est réalisée par corrélation d’images avec un ou plusieurs volume(s) de référence ce qui, par rapport à une annotation manuelle, permet un gain de temps et de ressources en opérateurs et limite de risque d’erreur et d’interprétation propre à un opérateur. Une telle détection est également plus fiable et plus précise car elle est basée sur une corrélation d’images avec un ou plusieurs volume(s) de référence de plusieurs voxels comprenant à minima un motif de toron, ce qui permet de bien discerner les torons et est moins sensible aux artefacts et différences locales de niveaux de gris au sein de l’image 3D de la pièce. Autrement dit, une telle détection permet une corrélation globale dans des portions de l’image 3D de la pièce, ce qui est plus pertinent que la binarisation ou la détermination de l’orientation locale des torons de l’art antérieur. Par ailleurs, une telle détection est plus précise car elle ne nécessite pas d’interpolation comme dans l’art antérieur ce qui permet de détecter des anomalies locales et/ou des discontinuités de torons qui étaient non visibles par interpolation. Selon un aspect de l’invention, le procédé de squelettisation est un procédé de squelettisation des torons d’une pièce en matériau composite comprenant une densité de torons supérieure à 50%, de préférence comprise entre 80% et 90%. De manière avantageuse, un tel procédé de squelettisation permet de bien discerner des torons compactés au sein d’une matrice, autrement dit contigus, ce que ne permettait pas la binarisation ou la détermination de l’orientation locale des torons de l’art antérieur. De préférence, l’étape de détermination de portions est mise en œuvre pour chaque voxel de l’image de la pièce, afin d’obtenir le squelette global des torons de la pièce dans l’image. Selon un aspect de l’invention, la phase de détection est mise en œuvre au moyen d’une pluralité de volumes de référence, les motifs de toron desdits volumes de référence s’étendant suivant différentes directions données. De manière avantageuse, les volumes de référence permettent de détecter les torons de toute direction dans l’image de la pièce, au moyen de différents motifs de toron d’orientation correspondante. Le procédé de squelettisation permet ainsi une squelettisation complète des torons de la pièce, en particulier de toutes directions. Selon un autre aspect de l’invention, la phase de détection est mise en œuvre au moyen d’une pluralité de volumes de référence, les motifs de toron desdits volumes de référence comportant différentes sections. On précise que la section d’un motif de toron se différencie de celle d’un autre motif de toron par sa forme et/ou sa taille. De manière avantageuse, les volumes de référence permettent de détecter les torons de toute section de l’image 3D de la pièce, en particulier de toute forme et toute taille, au moyen de différents motifs de toron de section correspondante. Le procédé de squelettisation permet ainsi une squelettisation complète des torons de la pièce, en particulier de toutes sections. Selon un aspect, au moins un volume de référence comprend autour du motif de toron un voisinage comprenant au moins un motif de toron voisin. Avantageusement, un volume de référence comprend un motif de toron central et un ou plusieurs motifs de torons voisins, entiers ou partiels, ce qui permet de réaliser une corrélation d’images qui n’est pas limitée au toron mais analyse également le voisinage. Une telle corrélation est ainsi plus fiable et plus précise. Selon un aspect de l’invention, la phase de détection est mise en œuvre au moyen d’une pluralité de volumes de référence, les motifs de toron desdits volumes de référence comportant des voisinages différents. On précise que le voisinage d’un volume de référence se différencie de celui d’un autre volume de référence par la présence ou non d’un ou plusieurs motif(s) de toron, par leur position et/ou par leur section. De manière avantageuse, les volumes de référence permettent de détecter des torons de différents voisinages de l’image 3D de la pièce, en particulier ayant des torons voisins ou non et plus ou moins compactés. Le procédé de squelettisation permet ainsi une squelettisation complète des torons de la pièce, en particulier de tous voisinages. De préférence, la phase de détection est mise en œuvre au moyen d’une pluralité de volumes de référence différents en ce que le motif de toron de chaque volume de référence s’étend suivant une direction donnée différente et/ou comporte une section différente et/ou comporte un voisinage différent. Avantageusement, une pluralité de volumes de référence comportant un panel de motifs de torons différents permet de squelettiser tous les torons de l’image 3D de la pièce, et ce de manière précise et fiable. Selon un aspect de l’invention, lorsque la phase de détection est mise en œuvre au moyen d’une pluralité de volumes de référence, le procédé de squelettisation comprend une étape de sélection du score de corrélation maximal associé à chaque voxel de l’image 3D de la pièce. Ceci permet avantageusement un gain en précision et en fiabilité en cherchant la meilleure corrélation parmi différents motifs de toron. Selon un autre aspect de l’invention, le procédé de squelettisation comprend une étape préliminaire de détermination d’au moins un volume de référence par sélection d’une portion ou d’une moyenne de portions de l’image 3D de la pièce comprenant un voxel central et comportant un motif de toron centré s’étendant suivant une direction donnée. Avantageusement, un volume de référence est adapté à l’image 3D de la pièce à analyser. En particulier, un volume de référence comporte la même plage de niveaux de gris et comporte un motif de toron similaire aux torons présents dans la pièce. Selon un autre aspect, au moins un volume de référence est obtenu par simulation informatique de tomographies. Avantageusement, un volume de référence peut être déterminé aisément sans analyse préliminaire de l’image 3D de la pièce, en particulier par un opérateur. En outre, un tel volume de référence est universel et ne correspond pas à un toron en particulier de la pièce, qui pourrait correspondre à une exception. Selon un aspect, le procédé de squelettisation comprend une étape préliminaire de filtrage de l’image 3D de la pièce, de préférence par un filtre de Gabor. De manière avantageuse, un tel filtre permet une meilleure perception des torons pouvant améliorer la détection des centres. Selon un aspect de l’invention, la phase de liaison comprend : une étape de calcul, pour chaque centre de section orienté, de la distance le séparant du centre de même orientation le plus proche et lorsque la distance est inférieure à une distance maximale, une étape de tracé de la courbe reliant chaque centre au centre de même orientation le plus proche. De manière avantageuse, une telle phase de liaison est simple et faible à réaliser. Une telle phase de liaison est en outre fiable et permet de vérifier que les centres ne sont pas trop éloignés les uns des autres, ce qui pourrait indiquer un éventuel défaut de toron, tel qu’une discontinuité. De préférence, lorsque la distance est supérieure à la distance maximale, un contrôle manuel visuel de la phase de détection des centres est réalisé par un opérateur. De manière avantageuse, l’expertise d’un opérateur pour l’analyse de zones locales où les centres sont anormalement espacés augmente la précision et la fiabilité du procédé de squelettisation. L’invention concerne également un système de squelettisation pour la mise en œuvre du procédé de squelettisation des torons d’une pièce en matériau composite dans une image 3D de ladite pièce tel que décrit précédemment, ledit système de squelettisation comprenant au moins une base de données comprenant au moins un volume de référence comportant un motif de toron centré s’étendant suivant une direction donnée, ledit système de squelettisation comprenant au moins un organe de calcul qui, pour chaque volume de référence, est configuré pour : déterminer des portions de l’image de dimension égale à celle du volume de référence, calculer un score de corrélation entre le volume de référence et chaque portion de l’image par comparaison du niveau de gris de leurs voxels respectifs, ledit score de corrélation étant associé au voxel central de chaque portion, de manière à obtenir un score de corrélation pour chaque voxel de l’image formant une matrice de scores de corrélation, et déterminer les centres de section orientés des torons de la pièce, lesdits centres correspondant aux voxels de l’image dont le score de corrélation correspond à un maximum local de la matrice de scores de corrélation, l’orientation desdits centres étant définie par la direction du motif de toron du volume de référence, ledit organe de calcul étant de plus configuré pour relier chaque centre au centre de même orientation le plus proche, afin de relier les centres associés à un même toron, de manière à obtenir le squelette des torons de la pièce. Avantageusement, un tel système de squelettisation permet d’obtenir automatiquement le squelette des torons d’une pièce, sans aide d’un opérateur, grâce à une base de données préalablement remplie. Une telle base de données peut en outre être améliorée au fur et à mesure de son utilisation, par retour d’expérience. De préférence, l’organe de calcul est configuré pour déterminer au moins un volume de référence par simulation informatique de tomographies. De préférence, l’organe de calcul est configuré pour filtrer l’image 3D de la pièce, préférentiellement par un filtre de Gabor. De préférence, dans le cas d’une pluralité de volumes de référence, l’organe de calcul est configuré pour sélectionner le score de corrélation maximal associé à chaque voxel de l’image 3D de la pièce. PRESENTATION DES FIGURES L’invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple, et se référant aux dessins annexés donnés à titre d’exemples non limitatifs, dans lesquels des références identiques sont données à des objets semblables et sur lesquels : est une représentation schématique en niveaux de gris d’une portion de la paroi d’une aube en matériau composite comprenant des torons selon l’art antérieur ; est une représentation schématique en perspective d’un volume tomographique de la portion de paroi d’aube de la figure 1 ; est une représentation schématique d’un procédé de squelettisation des torons de la portion de paroi d’aube de la figure 1 dans le volume tomographique de la figure 2 selon l’art antérieur ; est une représentation schématique d’une phase de liaison de centres de section orientés des torons selon le procédé de squelettisation de la figure 3 ; est une représentation schématique d’une phase de détection des centres de section orientés des torons selon le procédé de squelettisation de la figure 3 ; est une représentation schématique en perspective d’une portion de paroi d’aube en matériau composite comprenant des torons selon l’invention ; est une représentation schématique en niveaux de gris d’un volume tomographique de la portion de paroi d’aube de la figure 6 ; est une représentation schématique d’un procédé et d’un système de squelettisation des torons de la portion de paroi d’aube de la figure 6 dans le volume tomographique de la figure 7 selon un mode et une forme de réalisation de l’invention ; est une représentation schématique d’une phase de liaison de centres de section orientés des torons selon le procédé de squelettisation de la figure 8 ; et sont des représentations schématiques en perspective d’un volume de référence comprenant respectivement un motif de toron longitudinal et un motif de toron transversal pour la mise en œuvre d’une phase de détection des centres de section orientés des torons selon le procédé de squelettisation de la figure 8 ; , et sont des représentations schématiques d’une étape de détermination de portions du volume tomographique et d’une étape de calcul d’un score de corrélation pour trois voxels du volume tomographique selon le procédé de squelettisation de la figure 8 ; est une représentation schématique d’un procédé de squelettisation selon un autre mode de réalisation de l’invention comprenant une étape de détermination d’au moins un volume de référence, est une représentation schématique d’un procédé de squelettisation selon un autre mode de réalisation de l’invention au moyen de volumes de référence comportant des motifs de torons de différentes sections, est une représentation schématique d’une étape de sélection du score de corrélation maximal associé à un voxel du procédé de squelettisation de la figure 14, est une représentation schématique d’un procédé de squelettisation selon un autre mode de réalisation de l’invention au moyen de volumes de référence comportant des motifs de torons comportant différents voisinages et est une représentation schématique d’une étape de filtrage du volume de référence d’un procédé de squelettisation selon un mode de réalisation alternatif de l’invention. Il faut noter que les figures exposent l’invention de manière détaillée pour mettre en œuvre l’invention, lesdites figures pouvant bien entendu servir à mieux définir l’invention le cas échéant. DESCRIPTION DETAILLEE DE L’INVENTION On considère une aube en matériau composite 1, dont une portion de la paroi est représentée sur la figure 6, dont on cherche à déterminer la structure interne afin de vérifier les propriétés mécaniques de l’aube 1 et notamment sa résistance à un endommagement. De manière connue et en référence à la figure 6, une telle aube 1 comprend des fibres de renfort F, par exemple en carbone, noyées dans une matrice G, telle qu’une matrice organique à base de résine polymère, servant de liant et conférant à l’aube 1 sa forme extérieure. Les fibres de renfort F sont organisées en faisceaux de fibres, connus sous le terme de « torons T X , T Y », comprenant de l’ordre d’une centaine à un millier de fibres de renfort F tressées ou alignées les unes avec les autres. Les torons T X , T Y sont quant à eux agencés entre eux par tissage et/ou empilement, tel que dans l’exemple de la figure 6, un empilement de torons longitudinaux T X tissés avec des torons transversaux T Y selon un repère orthogonal (X, Y, Z) défini par la direction de l’empilement Z, la direction longitudinale X et la direction transversale Y. Dans l’exemple de la figure 6, les torons T X , T Y sont espacés les uns des autres à des fins de clarté, mais une aube 1 classiquement utilisée dans l’industrie aéronautique possède une densité de torons T X , T Y comprise entre 80% et 90%, c’est-à-dire que les torons T X , T Y sont agencés de manière très compactée au sein de la matrice G. Les torons T X , T Y d’une telle aube 1 sont ainsi en pratique majoritairement contigus les uns aux autres. De manière connue et en référence à la figure 7, la structure interne d’une telle aube 1, et plus précisément la forme et la position des torons T X , T Y de l’aube 1, est déterminée dans une image 3D en niveaux de gris N de l’aube 1 obtenue par acquisition tomographique à partir d’un grand nombre de radiographies 2D obtenues en bombardant l’aube 1 de rayons X suivant différents angles. Une telle image 3D, dite « volume tomographique V » par la suite, est divisée en voxels vox s’étendant selon un repère orthogonal i, j, k. Chaque voxel vox comprend un niveau de gris N donné, correspondant au coefficient d’absorption de rayons X moyen de la zone d’aube 1 illustrée sur le voxel vox. Les voxels vox de niveaux de gris N sombres correspondent ainsi à une zone d’aube 1 comprenant majoritairement un matériau de faible coefficient d’absorption des rayons X N2, à savoir de la matrice G. Les voxels vox de niveaux de gris N clairs correspondent quant à eux à une zone d’aube 1 comprenant majoritairement un matériau de grand coefficient d’absorption N1, à savoir des torons T X , T Y . A titre d’exemple, le voxel vox dans le coin supérieur gauche est un voxel de niveau de gris N sombre, car la zone d’aube 1 correspondante comprend environ 80% de matrice G (matériau de faible coefficient d’absorption N2) et 20% de torons T X , T Y (matériau de grand coefficient d’absorption N1). De manière générale, le niveau de gris N de chaque voxel vox de la figure 7 est à interpréter à partir du pourcentage de matériau de faible coefficient d’absorption N2 et de grand coefficient d’absorption N1. A noter que les voxels vox représentés sur la figure 7 sont de grande taille uniquement à des fins de clarté. Comme illustré sur les figures 8 et 9, l’invention concerne un procédé de squelettisation des torons T X , T Y d’une telle aube 1, qui est mis en œuvre dans le volume tomographique V en déterminant une courbe directrice reflétant la position de chaque toron T X , T Y , connue sous le terme de « squelette Sq X , Sq Y de toron T X , T Y ». Le procédé de squelettisation selon l’invention s’applique de manière générale à toute pièce en matériau composite et comprend : une phase de détection E D de centres de section orientés C X , C Y des torons T X , T Y dans le volume tomographique V puis une phase de liaison E L de chaque centre C X , C Y au centre C X , C Y de même orientation le plus proche, afin de relier les centres C X , C Y associés à un même toron T X , T Y , de manière à obtenir le squelette Sq X , Sq Y des torons T X , T Y . Selon l’invention et comme illustré sur les figures 10A et 10B, la phase de détection E D est mise en œuvre au moyen d’un ou plusieurs volume(s) de référence R X , R Y comprenant des voxels vox, et notamment un voxel central, qui s’étendent dans le repère orthogonal i, j, k du volume tomographique V et comprennent chacun un niveau de gris N, de manière à ce que chaque volume de référence R X , R Y comporte un motif de toron centré M X , M Y s’étendant suivant une direction donnée X, Y. Autrement dit, un volume de référence R X , R Y est une image 3D en niveaux de gris N avec un voxel central et un motif de toron M X , M Y central. On précise ici que par motif de toron, on désigne un motif d’une tranche transversale de toron inscrite dans le volume de référence. Autrement dit, le motif de toron représente un toron tronqué dans le sens de la longueur et dont la section est entièrement représentée. En reprenant le cas de l’aube 1, dans l’exemple des figures 10A et 10B, la phase de détection E D est mise en œuvre au moyen de deux volumes de référence R X , R Y , à savoir un volume de référence longitudinal R X comprenant un motif de toron longitudinal M X et un volume de référence transversal R Y comprenant un motif de toron transversal M Y. Comme cela sera présenté par la suite, les volumes de référence R X , R Y peuvent être plus nombreux et de natures différentes. Selon l’invention et comme illustré sur la figure 8, de tels volumes de référence R X , R Y permettent de mettre en œuvre la phase de détection E D en réalisant pour chaque volume de référence R X , R Y : une étape de détermination E1 de portions P X , P Y du volume tomographique V comportant un voxel central et de dimension égale à celle du volume de référence R X , R Y , une étape de calcul E2 d’un score de corrélation S X , S Y entre le volume de référence R X , R Y et chaque portion P X , P Y déterminée par comparaison du niveau de gris N de leurs voxels vox respectifs, ledit score de corrélation S X , S Y étant associé au voxel central de chaque portion P X , P Y , de manière à obtenir un score de corrélation S X , S Y pour chaque voxel vox du volume tomographique V formant une matrice de scores de corrélation, et une étape de détermination E3 des centres de section orientés C X , C Y des torons T X , T Y qui correspondent aux voxels centraux des portions P X , P Y dont le score de corrélation S X , S Y correspond à un maximum local de la matrice de scores de corrélation, leur orientation étant définie par la direction X, Y du motif de toron M X , M Y du volume de référence R X , R Y . Dans le cas de l’aube 1 de la figure 6 et comme illustré sur la figure 9, la phase de détection E D permet ainsi de déterminer des centres de section longitudinaux C X et des centres de section transversaux C Y à partir de scores de corrélation S X , S Y calculés respectivement entre le volume de référence longitudinal R X de la figure 10A et des portions longitudinales P X et entre le volume de référence transversal R Y de la figure 10B et des portions transversales P Y . En référence aux figures 8 et 9, la phase de liaison E L est quant à elle plus précisément mise en œuvre par : une étape de calcul E4, pour chaque centre de section orienté C X , C Y , de la distance d X , d Y le séparant du centre C X , C Y de même orientation le plus proche et lorsque la distance d X , d Y est inférieure à une distance maximale d Xmax , d Ymax , une étape de tracé E5 de la courbe reliant chaque centre C X , C Y au centre le plus proche C X , C Y de même orientation. Autrement dit, la phase de liaison E L permet de relier les centres de section longitudinaux C X entre eux de proche en proche, de même pour les centres de section transversaux C Y , et ainsi déterminer les squelettes de torons longitudinaux Sq X et les squelettes de torons transversaux Sq Y , dont certains sont représentés sur la figure 9. Comme illustré sur la figure 8, lorsque la distance d X , d Y calculée est supérieure à la distance maximale d Xmax , d Ymax , un opérateur H vérifie par contrôle visuel la détection des centres de section C X , C Y concernés dans le volume tomographique V et si besoin, annote manuellement des centres de section C X , C Y non détectés ou déplace des centres de section C X , C Y mal détectés. De manière avantageuse, le procédé de squelettisation selon l’invention permet un gain en temps et en ressources et limite le risque d’erreur et d’interprétation propre à l’opérateur, en particulier au niveau de la phase de détection E D des centres de section orientés C X , C Y . En effet, la phase de détection E D n’est plus opérée par annotation manuelle des centres C X , C Y sur la base d’une analyse visuelle suivie d’une interpolation entre les centres annotés manuellement, mais par corrélation d’images sur l’ensemble du volume tomographique V. Une telle corrélation d’images est en outre précise et fiable car elle compare un ou plusieurs volume(s) de référence orienté(s) R X , R Y à des portions du volume tomographique V de plusieurs voxels vox. La corrélation d’images attribue ainsi un score de corrélation S X , S Y aux voxels vox du volume tomographique V, et plus précisément aux voxels centraux des portions P X , P Y , qui prend en compte le voisinage des voxels vox, à savoir les voxels vox situés autour du voxel central de chaque portion P X , P Y. Le score de corrélation est ainsi plus pertinent car la corrélation d’images est peu sensible aux artefacts présents dans le volume tomographique V induisant des différences globales ou locales de plages de niveaux de gris N dans le volume tomographique V dues aux conditions d’acquisition. La corrélation d’images permet en outre de mettre en évidence des défauts locaux qui pouvaient être masqués par l’interpolation entre les centres annotés manuellement de l’art antérieur. En outre, prendre en compte le voisinage d’un voxel vox présente l’avantage de bien discerner des torons T X , T Y contigus et/ou compactés, ce que ne permettait pas les procédés de squelettisation de l’art antérieur utilisant la binarisation ou déterminant l’orientation locale des torons. Le procédé de squelettisation selon l’invention est donc adapté à des pièces en matériau composite de densité de torons T X , T Y importante comprise entre 80% et 90%. La fiabilité et la précision du procédé de squelettisation sont en outre renforcées par l’appel à un opérateur H en cas d’espacement anormal entre des centres C X , C Y afin qu’il vérifie localement la bonne détection des centres C X , C Y. On décrit par la suite plus précisément chacune des étapes du procédé de squelettisation selon l’invention appliqué à la détermination des torons T X , T Y de l’aube 1 de la figure 6 et mis en œuvre avec le volume de référence longitudinal R X de la figure 10A et le volume de référence transversal R Y de la figure 10B. A noter que dans l’exemple des figures 10A et 10B, les volumes de référence R X , R Y sont de très faible dimension : 3 x 3 x 1 voxels dans le repère i, j, k pour le volume de référence longitudinal R X et 3 x 1 x 3 voxels pour le volume de référence transversal R Y . En pratique, la dimension d’un volume de référence R X , R Y est plus grande dans la réalité, et ce pour une détection plus faible et plus précise des torons T X , T Y. A noter également que la description qui suit pourrait être adaptée à toute pièce en matériau composite en remplaçant l’aube 1 par ladite pièce. Les volumes de référence R X , R Y pourraient également être remplacés par d’autres volumes de référence en un nombre quelconque, de préférence choisis de sorte que leurs motifs de toron soient représentatifs des torons de la pièce. A titre d’exemple, les motifs de toron doivent être représentatifs de l’orientation des torons dans la pièce, à savoir la direction longitudinale X et la direction transversale Y des torons longitudinaux T X et transversaux T Y pour l’aube 1 de la figure 6. En référence aux figures 11A, 11B et 11C et comme décrit précédemment, le procédé de squelettisation débute par une étape de détermination E1 de portions longitudinales P X de dimension égale au volume de référence longitudinal R X et de portions transversales P Y de dimension égale au volume de référence transversal R Y . Chaque portion longitudinale P X et transversale P Y comporte un voxel central auquel sera associé le score de corrélation S X , S Y . A titre d’exemples, une première portion longitudinale P X1 de premier voxel central vox1 est représentée sur la figure 11A, une deuxième portion longitudinale P X2 de deuxième voxel central vox2 est représentée sur la figure 11B et une troisième portion longitudinale P X3 est représentée sur la figure 11C. Ces trois portions longitudinales P X1 , P X2 , P X3 sont de dimension 3 x 3 x 1 voxels dans le repère i, j, k de même que le volume de référence longitudinal R X de la figure 10A. De la même façon, les portions transversales P Y (non représentées) sont de dimension 3 x 1 x 3 voxels dans le repère i, j, k de même que le volume de référence transversal R Y de la figure 10B. Dans l’exemple des figures 10A et 10B, le volume de référence longitudinal R X et le volume de référence transversal R Y sont de dimensions différentes mais il va de soi que leurs dimensions pourraient être identiques auquel cas les portions longitudinales P X et les portions transversales P Y seraient identiques. De préférence, chaque voxel vox du volume tomographique V est le voxel central d’une portion longitudinale P X et d’une portion transversale P Y de manière à calculer un score de corrélation longitudinal S X et un score de corrélation transversal S Y en chaque voxel vox du volume tomographique V. On remarquera donc que les portions P X , P Y ne sont pas distinctes les unes des autres mais comportent chacune des voxels vox communs à d’autres portions P X , P Y . Par exemple, la première portion longitudinale P X1 (fig. 11A) et la deuxième portion transversale P X2 (fig. 11B) comprennent six voxels communs, dont le premier voxel central vox1 et le deuxième voxel central vox2. On remarquera également que pour certains voxels centraux situés aux extrémités du volume tomographique V, les portions P X , P Y du volume tomographique V associées « débordent » du volume tomographique V, c’est-à-dire qu’elles comprennent des voxels vox du volume tomographique V dits « voxels réels » et des voxels fictifs. Comme décrit précédemment et toujours en référence aux figures 11A, 11B et 11C, une fois les portions longitudinales P X et transversales P Y déterminées, le procédé de squelettisation comprend une étape de calcul E2 d’un score de corrélation longitudinal S X entre le volume de référence longitudinal R X et chaque portion longitudinale P X et d’un score de corrélation transversal S Y entre le volume de référence transversal R Y et chaque portion transversale P Y. Pour une portion P X , P Y du volume tomographique V donnée, un tel score de corrélation S X , S Y est calculé en comparant le niveau de gris N de chaque voxel vox du volume de référence R X , R Y avec celui du voxel vox correspondant de ladite portion P X , P Y suivant une fonction de corrélation normalisée f. Pour un voxel vox de coordonnées a, b, c dans le repère i, j, k et pour un volume de référence R et une portion P de dimension a1xb1xc1voxels, la fonction de corrélation normalisée f s’exprime : A noter que pour les voxels centraux de portions P X , P Y situés aux extrémités du volume tomographique V, le score de corrélation S X , S Y est calculé en comparant le niveau de gris N de chaque voxel vox réel desdites portions P X , P Y uniquement avec celui des voxels vox correspondant du volume de référence R X , R Y . Autrement dit, la corrélation est effectuée uniquement sur les voxels non fictifs dans le cas des voxels centraux situés aux extrémités du volume tomographique V. Le score de corrélation longitudinal S X et le score de corrélation transversal S Y obtenus pour chaque voxel central de portion P X , P Y sont compris entre 0 et 1, la corrélation étant d’autant plus faible que le score S X , S Y est proche de 0 et d’autant plus forte que le score S X , S Y est proche de 1. A titre d’exemples, le calcul d’un premier score de corrélation longitudinal S X1 pour le premier voxel central vox1 est illustré sur la figure 11A, celui d’un deuxième score de corrélation longitudinal S X2 pour le deuxième voxel central vox2 est illustré sur la figure 11B et celui d’un troisième score de corrélation longitudinal S X3 pour le troisième voxel central vox3 est illustré sur la figure 11C. Dans cet exemple, la troisième portion longitudinale P X3 comportant un motif de toron centré contrairement aux première et deuxième portions longitudinales P X1 , P X2 , le troisième voxel central vox3 comporte donc le plus grand score de corrélation longitudinal S X . Comme décrit précédemment, une fois les scores de corrélation longitudinal S X et transversal S Y calculés pour chaque voxel vox du volume tomographique V, le procédé de squelettisation comprend une étape de détermination E3 des centres de section longitudinaux C X et transversaux C Y des torons T X , T Y de l’aube 1. Pour cela, les maxima locaux de la matrice des scores de corrélation S X , S X1 , S X2 , S X3 , S Y du volume tomographique V sont identifiés, notamment par transformation H-maxima Hmax, comme référencée sur la figure 8. Une telle transformation H-maxima est connue en soi de l’homme du métier et une description peut être trouvée dans le chapitre 6 dans la seconde édition parue en 2003 du livre « Morphological Image Analysis : Principles and Applications » de P. Soille. En reprenant l’exemple précédent, seul le troisième score de corrélation longitudinal S X3 correspond à un maximum local parmi les trois scores de corrélation longitudinaux S X1 , S X2 , S X3 . Comme décrit précédemment, les voxels vox, vox1, vox2, vox3 du volume tomographique V dont le score de corrélation S X , S X1 , S X2 , S X3 , S Y correspond à un maximum local forment chacun un centre C X , C Y . Ces différents centres C X , C Y sont orientés, l’orientation d’un centre C X , C Y étant définie par la direction longitudinale X ou transversale Y du motif de toron M X , M Y du volume de référence R X , R Y d’où est issu le score de corrélation S X , S Y . Toujours en reprenant l’exemple précédent, seul le troisième voxel central vox3 est désigné centre de section longitudinal C X . De préférence, la transformation H-maxima Hmax est mise en œuvre au moyen d’une hauteur minimale de contraste suffisamment grande pour ne pas déterminer de centre C X , C Y correspondant à des artefacts du volume tomographique V et suffisamment petite pour déterminer l’ensemble des centres de torons T X , T Y de l’aube 1. Dans les faits, la valeur de la hauteur minimale de contraste est directement liée au choix des volumes de référence R X , R Y. Plus les volumes de référence R X , R Y comportent un motif de toron M X , M Y proche des torons T X , T Y de l’aube 1 et plus la valeur de la hauteur minimale de contraste pourra être grande, et donc la corrélation discriminante. En pratique, la valeur de la hauteur minimale de contraste est choisie par retour d’expérience. A noter que la hauteur minimale de contraste peut être différente pour chaque volume de référence R X , R Y . Dans cet exemple, il pourrait donc y avoir une hauteur minimale de contraste longitudinale et un transversale. Les trois étapes précédemment décrites permettent la détection E D des centres de section longitudinaux C X et transversaux C Y respectivement des torons longitudinaux T X et transversaux T Y de l’aube 1. Comme décrit précédemment et en référence à la figure 9, la détermination des squelettes Sq X , Sq Y des torons T X , T Y de l’aube 1 est ensuite réalisée par liaison E L des centres C X , C Y qui ont été précédemment détectés. En référence aux figures 8 et 9, la phase de liaison E L débute par une étape de calcul E 4 , pour chaque centre de section longitudinal C X , de la distance longitudinale d X le séparant du centre de section longitudinal C X le plus proche, et pour chaque centre de section transversal C Y , de la distance transversale d Y le séparant du centre de section transversal C Y le plus proche. Une telle étape de calcul E 4 est mise en œuvre par transformation de distance et permet de déterminer les centres C X , C Y appartenant à un même toron T X , T Y de manière simple et fiable. Toujours en référence aux figures 8 et 9, la phase de liaison E L se poursuit par une étape de tracé E 5 des squelettes longitudinaux Sq X des torons longitudinaux T X de l’aube 1 en reliant les centres longitudinaux C X les plus proches les uns des autres et une étape de tracé E 5 des squelettes transversaux Sq Y des torons transversaux T Y de l’aube 1 en reliant les centres transversaux C Y les plus proches les uns des autres. Toutefois, lorsque la distance longitudinale d X , séparant deux centres longitudinaux C X les plus proches est supérieure à une distance maximale longitudinale d Xmax , c’est-à-dire que les centres longitudinaux C X sont anormalement espacés, un opérateur H vérifie par contrôle visuel la détection des centres longitudinaux C X concernés dans le volume tomographique V et si besoin, annote manuellement des centres longitudinaux C X non détectés ou déplace des centres longitudinaux C X mal détectés. Une telle vérification peut également être réalisée pour des centres longitudinaux C X anormalement rapprochés. Il est réalisé de même pour les centres transversaux C Y . En pratique, la distance maximale longitudinale d Xmax et la distance maximale transversale d Ymax sont choisis en fonction de la courbure physique des torons, de la résolution du volume tomographique V et de l’orientation des torons T X , T Y dans le volume tomographique V. Il va de soi que la distance maximale longitudinale d Xmax et la distance maximale transversale d Ymax peuvent être identiques. Pour résumer, les squelettes Sq X , Sq Y des torons T X , T Y de l’aube 1 sont déterminés par phases de détection E D puis de liaison E L de centres de section longitudinaux C X et transversaux C Y respectivement des torons longitudinaux T X et transversaux T Y . Les centres longitudinaux C X sont déterminés au moyen d’un volume de référence longitudinal R X pour lequel une corrélation est effectuée avec des portions longitudinales P X du volume tomographique V. De même pour les centres transversaux C Y . Les voxels centraux vox1, vox2, vox3 des portions longitudinales P X et transversales P Y avec les meilleurs scores de corrélation S X1 , S X2 , S X3 définissent respectivement les centres longitudinaux C X et transversaux C Y , qui sont ensuite reliés de proche en proche. De manière avantageuse, la corrélation est effectuée sur la base d’un ou plusieurs volumes de référence R X , R Y comprenant un motif de toron M X , M Y central et non sur la base d’un ou plusieurs voxels de référence, ce qui rend la corrélation fiable et précise par rapport aux artefacts et dans le cas d’une densité de torons T X , T Y importante. Par ailleurs, le choix du ou des volumes de référence R X , R Y est déterminant pour la bonne détection des torons T X , T Y de l’aube 1. En effet, les volumes de référence R X , R Y doivent être représentatifs de l’ensemble des torons T X , T Y de l’aube 1 pour permettre une détection systématique et précise. Ainsi, dans l’exemple précédemment décrit, un volume de référence longitudinal R X est utilisé pour détecter les centres longitudinaux C X de torons longitudinaux T X et un volume de référence transversal R Y est utilisé pour détecter les centres transversaux C Y de torons transversaux T Y . Pour que les volumes de référence R X , R Y soient les plus représentatifs possibles de l’ensemble des torons T X , T Y de l’aube 1, selon un mode de réalisation préféré de l’invention illustré sur la figure 13, le procédé de squelettisation comprend une étape préliminaire de détermination E0 du volume de référence longitudinal R X et du volume de référence transversal R Y directement dans le volume tomographique V. Pour cela, une portion du volume tomographique V comportant un voxel central et un motif de toron M X , M Y central représentatif d’un toron longitudinal T X et d’un toron transversal T Y de l’aube 1 est sélectionnée par un opérateur H. La détermination du volume de référence longitudinal R X et du volume de référence transversal R Y est ainsi simple et facile à mettre en œuvre et permet d’obtenir des volumes de référence R X , R Y représentatifs des torons T X , T Y de l’aube 1. De manière alternative, les volumes de référence R X , R Y peuvent également être déterminés par simulation informatique de tomographies, ce qui présente un gain de temps en évitant de devoir analyser le volume tomographique V de manière préliminaire. En pr Par ailleurs, dans la description précédente, la détection E D des centres C X , C Y des torons T X , T Y est effectuée uniquement au moyen de deux volumes de référence R X , R Y , à savoir un volume de référence longitudinal R X représentatif de l’ensemble des torons longitudinaux T X et un volume de référence transversal R Y représentatif de l’ensemble des torons transversaux T Y . En pratique, les torons longitudinaux T X et les torons transversaux T Y possèdent des sections de différentes formes, notamment elliptiques, et tailles et sont plus ou moins rapprochés de torons voisins suivant leur position dans l’aube et/ou suivant la compaction locale. Ainsi, selon un mode de réalisation préféré de l’invention illustré sur la figure 14, pour que les volumes de référence R X , R Y soient les plus représentatifs possibles de l’ensemble des torons T X , T Y de l’aube 1, plusieurs volumes de référence longitudinaux R X , R X’ et plusieurs volumes de référence transversaux R Y comportant un motif de toron M X , M X’ , M Y de sections S différentes sont utilisés pour mettre en œuvre la détection E D. Les sections S peuvent différer les unes des autres par leur taille mais également par leur forme, de manière à être représentatives de chacun des torons T X , T Y de l’aube 1. Dans l’exemple de la figure 14, deux volumes de référence longitudinaux R X , R X’ sont utilisés pour détecter les centres longitudinaux C X des torons longitudinaux T X de l’aube 1. Un premier volume de référence longitudinal R X comporte un motif de toron longitudinal M X de petite section elliptique S et un deuxième volume de référence longitudinal R X’ comporte un motif de toron longitudinal M X’ de grande section elliptique S’. L’étape de détermination E1 de portions longitudinales P X , l’étape de calcul E2 de scores de corrélation longitudinaux S X et l’étape de détermination E3 de centres longitudinaux C X sont ainsi mises en œuvre pour le premier volume de référence longitudinal R X et pour le deuxième volume de référence longitudinal R X’ . Dans cet exemple, les deux volumes de référence longitudinaux R X , R X’ possèdent la même dimension donc la corrélation peut être mise en œuvre sur les mêmes portions longitudinales P X mais il va de soi que les deux volumes de référence longitudinaux R X , R X’ pourraient comporter une dimension différente. Dans tous les cas, à l’issue de l’étape de calcul E2, les voxels centraux, et notamment le troisième voxel central vox3 de la troisième portion longitudinale P X3 illustrée sur la figure 14, comportent deux scores de corrélation longitudinaux S X3 , S X3’ . De préférence, tel qu’illustré sur les figures 14 et 15, avant la mise en œuvre de l’étape de détermination E 3 , la phase de détection E D comprend une étape de sélection E2’ du score de corrélation longitudinal S X3 , S X3 maximal associé à chaque voxel vox3. Pour le cas du troisième voxel central vox3, il s’agit ainsi du score de corrélation S X3’ issu du deuxième volume de référence longitudinal R X’ comportant un motif de toron longitudinal M X’ de grande section elliptique S’. L’étape de détermination E3 des centres longitudinaux C X est ainsi mise en œuvre à partir du score de corrélation maximal S X3’ associé à chaque voxel vox3. On décrit par la suite un autre mode de réalisation préféré de l’invention pour que les volumes de référence R X , R Y soient les plus représentatifs possibles de l’ensemble des torons T X , T Y de l’aube 1. Selon ce mode de réalisation illustré sur la figure 16, plusieurs volumes de référence longitudinaux R X , R X’ et plusieurs volumes de référence transversaux R Y comportent un motif de toron M X , M X’ ’ , M Y avec différents voisinages. Le voisinage d’un motif de toron M X , M X’ ’ , M Y peut différer par la présence ou non de motifs de torons voisins M V ainsi que par la forme et la position de ces derniers, de manière à ce que les volumes de référence R X , R X’’ , R Y soient représentatifs de chacun des torons T X , T Y de l’aube 1, et en particulier de leurs différents degrés de compaction dans l’aube 1. Un volume de référence R X , R X’ , R Y comprend ainsi un ou plusieurs motif(s) de toron, à savoir à minima un motif de toron centré M X , M X’ , M Y et éventuellement un ou plusieurs motif(s) de torons voisins M V entiers ou partiels. Dans l’exemple de la figure 16, deux volumes de référence longitudinaux R X , R X’ sont utilisés pour détecter les centres longitudinaux C X des torons longitudinaux T X de l’aube 1. Un premier volume de référence longitudinal R X comporte un motif de toron longitudinal M X sans voisinage et un deuxième volume de référence longitudinal R X’ comporte un motif de toron longitudinal M X’ avec un voisinage comprenant deux motifs de torons voisins M V partiels. De même que pour le mode de réalisation précédent illustré sur la figure 14, l’étape de détermination E1 de portions longitudinales P X , l’étape de calcul E2 de scores de corrélation longitudinaux S X et l’étape de détermination E3 de centres longitudinaux C X sont ainsi mises en œuvre pour le premier volume de référence longitudinal R X et pour le deuxième volume de référence longitudinal R X’ ’ . De même que pour le mode de réalisation précédent illustré sur la figure 14, les deux volumes de référence longitudinaux R X , R X’ ’ possèdent la même dimension donc la corrélation peut être mise en œuvre sur les mêmes portions longitudinales P X mais il va de soi que les deux volumes de référence longitudinaux R X , R X’ pourraient comporter une dimension différente. De préférence, comme décrit précédemment et illustré sur les figures 15 et 16, une étape de sélection E2’ du score de corrélation S X3’’ maximal est mise en œuvre avant l’étape de détermination E3 des centres longitudinaux C X . De manière générale, les modes de réalisation présentés sur les figures 10A, 10B, 14 et 16 sont combinables entre eux. Autrement dit, en pratique, la phase de détection E D des centres de section orientés C X , C Y des torons T X , T Y de l’aube 1, et plus généralement de toute pièce en matériau composite, est mise en œuvre avec plusieurs volumes de référence R X , R X’ , R X’’ , R Y comprenant un motif de toron M X , M X’ , M X’’ , M Y selon différentes directions X, Y, de différentes sections S et de différents voisinages. Selon un mode de réalisation alternatif de l’invention illustré sur la figure 17, le procédé de squelettisation comprend de manière préliminaire une étape de filtrage E0’ du volume tomographique V, par exemple au moyen d’un filtre de Gabor G, afin de mettre en exergue les motifs de torons dans le volume tomographique V. Une telle étape de filtrage E0’ est exempte de binarisation comme dans l’art antérieur, qui ne permet pas de distinguer des torons contigus les uns des autres. L’invention concerne également un système de squelettisation permettant la mise en œuvre du procédé de squelettisation précédemment décrit. En référence à la figure 8, dans l’exemple de l’aube 1, un tel système de squelettisation comprend un organe de calcul 2, tel qu’un ordinateur, ainsi qu’une base de données 3. La base de données 3 est configurée pour stocker les volumes de référence R X , R X’ , R X’’ , R Y acquis par simulation ou directement sur le volume tomographique V, ainsi que la hauteur de contraste minimale pour la mise en œuvre de la transformation H-maxima Hmax et la distance maximale d Xmax , d Ymax . Pour chaque volume de référence R X , R X’ , R X’’ , R Y , l’organe de calcul 2 est quant à lui configuré pour : déterminer des portions P X , P Y , P X1 , P X2 , P X3 du volume tomographique V de dimension égale à celle du volume de référence R X , R X’ , R X’’ , R Y , calculer un score de corrélation S X , S Y , S X 1 , S X 2 , S X 3 entre le volume de référence R X , R X’ , R X’’ , R Y et chaque portion P X , P Y , P X1 , P X2 , P X3 du volume tomographique V par comparaison du niveau de gris N de leurs voxels vox respectifs, ledit score de corrélation S X , S Y , S X 1 , S X 2 , S X 3 étant associé au voxel central vox1, vox2, vox3 de chaque portion P X , P Y , P X1 , P X2 , P X3 , de manière à obtenir un score de corrélation S X , S Y , S X 1 , S X 2 , S X 3 pour chaque voxel vox du volume tomographique V formant une matrice de scores de corrélation, et déterminer les centres de section orientés C X , C Y des torons T X , T Y de l’aube 1, lesdits centres C X , C Y correspondant aux voxels centraux vox1, vox2, vox3 des portions P X , P Y , P X1 , P X2 , P X3 dont le score de corrélation S X , S Y , S X 1 , S X 2 , S X 3 correspond à un maximum local de la matrice de scores de corrélation, l’orientation desdits centres C X , C Y étant définie par la direction X, Y du motif de toron M X , M Y , M X’ , M X’’ du volume de référence R X , R X’ , R X’’ , R Y . L’organe de calcul 2 est de plus configuré pour relier chaque centre C X , C Y au centre C X , C Y de même orientation le plus proche, afin de relier les centres C X , C Y associés à un même toron T X , T Y , de manière à obtenir le squelette Sq X , Sq Y des torons T X , T Y de l’aube 1. Selon un aspect, l’organe de calcul 2 est configuré pour déterminer les volumes de référence R X , R X’ , R X’’ , R Y par simulation tomographique. Selon un autre aspect, l’organe de calcul 2 est configuré pour filtrer le volume de tomographique V, par exemple par un filtre de Gabor G. Selon un autre aspect, dans le cas d’une pluralité de volumes de référence R X , R X’ , R X’’ , R Y , l’organe de calcul 2 est configuré pour sélectionner le score de corrélation maximal associé à chaque voxel du volume tomographique V. Il va de soi que le système de squelettisation pourrait comprendre une pluralité d’organes de calcul 2 configurés chacun pour effectuer une partie des tâches énumérées ci-dessus, ainsi qu’une pluralité de bases de données 3, intégrées ou non à un organe de calcul 2. Procédé de squelettisation des torons (T X , T Y ) d’une pièce en matériau composite dans une image 3D (V) de ladite pièce, ladite image (V) comprenant une pluralité de voxels (vox, vox1, vox2, vox3) s’étendant dans un repère orthogonal (i, j, k) et comprenant chacun un niveau de gris (N), ledit procédé comprenant : une phase de détection (E D ) d’une pluralité de centres de section orientés (C X , C Y ) des torons (T X , T Y ) de la pièce dans l’image (V) et une phase de liaison (E L ) de chaque centre (C X , C Y ) au centre (C X , C Y ) de même orientation le plus proche, afin de relier les centres (C X , C Y ) associés à un même toron (T X , T Y ), de manière à obtenir le squelette (Sq X , Sq Y ) des torons (T X , T Y ) de ladite pièce, ledit procédé étant caractérisé par le fait que la phase de détection (E D ) est mise en œuvre au moyen d’au moins un volume de référence (R X , R X’ , R X’’ , R Y ) comprenant une pluralité de voxels (vox), lesdits voxels (vox) s’étendant dans le repère orthogonal (i, j, k) et comprenant chacun un niveau de gris (N) de manière à ce que le volume de référence (R X , R X’ , R X’’ ,R Y ) comporte un motif de toron centré (M X , M X’ , M X’’ , M Y ) s’étendant suivant une direction donnée (X, Y), ladite phase de détection (E D ) comprenant, pour chaque volume de référence (R X , R X’ , R X’’ ,R Y ) : une étape de détermination (E1) de portions (P X , P Y , P X1 , P X2 , P X3 ) de l’image (V) comportant un voxel central (vox1, vox2, vox3) et de dimension égale à celle du volume de référence (R X , R X ’ , R X ’’ , R Y ), une étape de calcul (E2) d’un score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) entre le volume de référence (R X , R X ’ , R X ’’ , R Y ) et chaque portion (P X , P Y , P X1 , P X2 , P X3 ) de l’image (V) par comparaison du niveau de gris (N) de leurs voxels (vox) respectifs, ledit score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) étant associé au voxel central (vox1, vox2, vox3) de chaque portion (P X , P Y , P X1 , P X2 , P X3 ), de manière à obtenir un score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) pour chaque voxel (vox) de l’image (V) formant une matrice de scores de corrélation, et une étape de détermination (E3) des centres de section orientés (C X , C Y ) des torons (T X , T Y ) de la pièce, lesdits centres (C X , C Y ) correspondant aux voxels (vox) de l’image (V) dont le score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) correspond à un maximum local de la matrice de scores de corrélation, l’orientation desdits centres (C X , C Y ) étant définie par la direction (X, Y) du motif de toron (M X , M Y ) du volume de référence (R X , R X ’ , R X ’’ , R Y ). Procédé de squelettisation selon la revendication 1, des torons (T X , T Y ) d’une pièce en matériau composite comprenant une densité de torons (T X , T Y ) supérieure à 50%, de préférence comprise entre 80% et 90%. Procédé de squelettisation selon l’une des revendications 1 et 2, dans lequel la phase de détection (E D ) est mise en œuvre au moyen d’une pluralité de volumes de référence (R X , R X ’ , R X ’’ , R Y ), les motifs de toron (M X , M X ’ , M X ’’ , M Y ) desdits volumes de référence (R X , R X ’ , R X ’’ , R Y ) s’étendant suivant différentes directions données (X, Y). Procédé de squelettisation selon l’une des revendications 1 à 3, dans lequel la phase de détection (E D ) est mise en œuvre au moyen d’une pluralité de volumes de référence (R X , R X ’ , R X ’’ , R Y ), les motifs de toron (M X , M X ’ , M X ’’ , M Y ) desdits volumes de référence (R X , R X ’ , R X ’’ , R Y ) comportant différentes sections (S, S’). Procédé de squelettisation selon l’une des revendications 1 à 4, dans lequel au moins un volume de référence (R X , R X ’ , R X ’’ , R Y ) comprend autour du motif de toron (M X , M X ’ , M X ’’ , M Y ) un voisinage comprenant au moins un motif de toron voisin (M V ). Procédé de squelettisation selon la revendication 5 prise seule ou en combinaison avec l’une des revendications 1 à 4, dans lequel la phase de détection (E D ) est mise en œuvre au moyen d’une pluralité de volumes de référence (R X , R X ’ , R X ’’ , R Y ), les motifs de toron (M X , M X ’ , M X ’’ , M Y ) desdits volumes de référence (R X , R X ’ , R X ’’ , R Y ) comportant des voisinages différents. Procédé de squelettisation selon l’une des revendications 1 à 6, dans lequel, lorsque la phase de détection (E D ) est mise en œuvre au moyen d’une pluralité de volumes de référence (R X , R X ’ , R X ’’ , R Y ), le procédé comprend une étape de sélection (E2’) du score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) maximal associé à chaque voxel (vox, vox1, vox2, vox3) de l’image 3D (V) de la pièce. Procédé de squelettisation selon l’une des revendications 1 à 7, comprenant une étape préliminaire de détermination (E0) d’au moins un volume de référence (R X , R X ’ , R X ’’ , R Y ) par sélection d’une portion ou d’une moyenne de portions de l’image (V) comprenant un voxel central (vox1, vox2, vox3) et comportant un motif de toron centré (M X , M X ’ , M X ’’ , M Y ) s’étendant suivant une direction donnée (X, Y). Procédé de squelettisation selon l’une des revendications 1 à 8, dans lequel la phase de liaison (E L ) comprend : une étape de calcul (E4), pour chaque centre de section orienté (C X , C Y ), de la distance (d X , d Y ) le séparant du centre (C X , C Y ) de même orientation le plus proche et lorsque la distance (d X , d Y ) est inférieure à une distance maximale (d X max , d Ymax ), une étape de tracé (E5) de la courbe reliant chaque centre (C X , C Y ) au centre (C X , C Y ) de même orientation le plus proche. Système de squelettisation pour la mise en œuvre du procédé de squelettisation des torons (T X , T Y ) d’une pièce en matériau composite dans une image 3D (V) de ladite pièce selon l’une des revendications 1 à 9, ledit système de squelettisation comprenant au moins une base de données (3) comprenant au moins un volume de référence (R X , R X ’ , R X ’’ , R Y ) comportant un motif de toron centré (M X , M X ’ , M X ’’ , M Y ) s’étendant suivant une direction donnée (X, Y), ledit système de squelettisation comprenant au moins un organe de calcul (2) qui, pour chaque volume de référence (R X , R X ’ , R X ’’ , R Y ), est configuré pour : déterminer des portions (P X , P Y , P X1 , P X2 , P X3 ) de l’image (V) de dimension égale à celle du volume de référence (R X , R X ’ , R X ’’ , R Y ), calculer un score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) entre le volume de référence (R X , R X ’ , R X ’’ , R Y ) et chaque portion (P X , P Y , P X1 , P X2 , P X3 ) de l’image (V) par comparaison du niveau de gris (N) de leurs voxels (vox) respectifs, ledit score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) étant associé au voxel central (vox1, vox2, vox3) de chaque portion (P X , P Y , P X1 , P X2 , P X3 ), de manière à obtenir un score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) pour chaque voxel (vox) de l’image (V) formant une matrice de scores de corrélation, et déterminer les centres de section orientés (C X , C Y ) des torons (T X , T Y ) de la pièce, lesdits centres (C X , C Y ) correspondant aux voxels centraux (vox1, vox2, vox3) des portions (P X , P Y , P X1 , P X2 , P X3 ) dont le score de corrélation (S X , S Y , S X1 , S X2 , S X3 , S X3’ , S X3’’ ) correspond à un maximum local de la matrice de scores de corrélation, l’orientation desdits centres (C X , C Y ) étant définie par la direction (X, Y) du motif de toron (M X , M X ’ , M X ’’ , M Y ) du volume de référence (R X , R X ’ , R X ’’ , R Y ), ledit organe de calcul (2) étant de plus configuré pour relier chaque centre (C X , C Y ) au centre (C X , C Y ) de même orientation le plus proche, afin de relier les centres (C X , C Y ) associés à un même toron (T X , T Y ), de manière à obtenir le squelette (Sq X , Sq Y ) des torons (T X , T Y ) de la pièce.