Procédé (100) de contrôle non destructif d'une pièce en matériau composite par squelettisation de torons (Tx, Ty) de la pièce dans une image tridimensionnelle de ladite pièce, les torons étant associés à une même direction prédéterminée (X,Y), l'image tridimensionnelle comprenant des tranches (TRx, TRy) successives orthogonales à la direction prédéterminée, ledit procédé comprenant: Déterminer (102), pour chaque tranche (TRx, TRy) de l'image tridimensionnelle, des positions de centres (Cx, Cy) de sections (SECx, SECy) de torons (Tx, Ty), Pour au moins une tranche (TRx, TRy) parmi les tranches de l'image tridimensionnelle, corriger (103) au moins une position de centre (Cx, Cy) de section (SECx, SECy) de toron parmi les positions de centres (Cx, Cy) de sections de torons déterminées pour ladite au moins une tranche (TRx, TRy) ladite correction comprenant une segmentation des sections de torons dans une région de ladite au moins une tranche ayant un taux de compaction de torons inférieur à un seuil prédéfini, Déterminer (104) une pluralité de squelettes (Sx, Sy) de torons (Tx, Ty) en liant, pour plusieurs tranches (Tx, Ty) de l’image tridimensionnelle, un centre (Cx, Cy) de section de toron d’une tranche courante de l’image tridimensionnelle à un centre de section (Cx, Cy) de toron d’une tranche de l’image tridimensionnelle parmi des tranches voisines de la tranche courante, Emettre (111) une alerte lorsqu’une position d'au moins un squelette déterminé s'écarte d'une position d'un squelette prédéterminée. Procédé (100) selon la revendication précédente selon lequel, la correction de la au moins une position de centre (Cx, Cy) de section comprend : Localiser une région de ladite au moins une tranche (TRx, TRy) dont le taux de compaction de torons est inférieur au seuil prédéfini, Déterminer une pluralité de contours de sections de torons (Tx, Ty) dans la région localisée, Calculer une position d’un barycentre de chaque contour de section de toron déterminé, Pour chaque section de toron dont le contour a été déterminé, remplacer la position du centre de la section de toron par la position calculée du barycentre du contour de ladite section. Procédé (100) selon l'une quelconque des revendications précédentes selon lequel la liaison du centre (Cx, Cy) de section de toron de la tranche courante de l’image tridimensionnelle au centre de section (Cx, Cy) de toron de la tranche de l’image tridimensionnelle parmi les tranches voisines de la tranche courante est réalisée en minimisant une distance entre la position du centre de section de toron de la tranche courante de l’image tridimensionnelle et la position des centres de section de torons d’au moins une tranche de l'image tridimensionnelle parmi les tranches voisines de la tranche courante . Procédé (100) selon l'une quelconque des revendications précédentes selon lequel chaque squelette (Sx, Sy) de toron s'étend entre une première extrémité du squelette (Sx, Sy) et une deuxième extrémité du squelette (Sx, Sy), le procédé (100) comprenant en outre, pour chaque squelette (Sx, Sy) de toron déterminé et avant l’émission de l'alerte : Déterminer (105) à partir des centres (Cx, Cy) de sections de torons (Tx, Ty) dudit chaque squelette, un chemin géodésique entre la première extrémité du squelette et la deuxième extrémité dudit chaque squelette; et Corriger (106) ledit chaque squelette à partir du chemin géodésique déterminé. Procédé (100) selon la revendication précédente dans lequel la détermination (105) du chemin géodésique est réalisée par une transformation de distance géodésique. Procédé (100) selon l'une quelconque des revendications précédentes dans lequel l'écart entre la position d'au moins un des squelettes obtenus et la position prédéterminée d'un squelette est estimé via : un indice global, appelé indice global de confiance sur la propagation du squelette, calculé à partir d’un écart entre une direction de propagation du squelette obtenu et la direction prédéterminée selon laquelle s'étend le toron associé audit squelette et/ou; un indice global, appelé indice global de confiance sur la périodicité du squelette, calculé à partir d’un écart entre une périodicité dudit squelette et une périodicité prédéterminée et/ou; un indice local de confiance sur les niveaux de gris du squelette, calculé à partir d'un écart entre des niveaux de gris de voxels formant le squelette obtenu et des niveaux de gris prédéterminés et/ou, un indice local de confiance sur la proximité du squelette avec ses voisin, calculé à partir d'une distance du squelette obtenu par rapport à des squelettes voisins. Procédé (100) selon l'une quelconque des revendications précédentes selon lequel la direction prédéterminée correspond à une direction de chaine ou de trame. Procédé (100) selon l'une quelconque des revendications précédentes selon lequel l'image tridimensionnelle est une image par tomographie à rayons X comprenant une pluralité de voxels, chaque voxel étant associé à un niveau de gris N respectif. Procédé (100) selon l'une quelconque des revendications précédentes selon lequel les torons (Tx, Ty) sont noyés dans une résine et chaque voxel est associé à une classe respective indiquant si ledit voxel correspond à de la résine ou non et le taux de compaction d’une région d’une tranche de l’image tridimensionnelle est déterminé à partir d'un taux de voxels correspondant à de la résine parmi les voxels de ladite région. Procédé (100) selon la revendication précédente et selon la revendication 2 selon lequel la localisation de la région de ladite tranche dont le taux de compaction de torons est inférieur au seuil prédéfini comprend: déterminer, pour au moins un groupe de voxels voisins de ladite tranche, un taux de compaction respectif à partir des classes associées aux voxels dudit au moins un groupe, et comparer le taux de voxels déterminé au seuil prédéfini. Procédé (100) selon l'une quelconque des revendications précédentes caractérisé en ce que la détermination (102) de la pluralité de centres de section de torons comprend : Appliquer (1024) un filtre de différences de Gaussiennes tridimensionnel sur l'image tridimensionnelle, pour obtenir une première image filtrée tridimensionnelle, Appliquer (1025) un filtre Laplacien de Gaussiennes bidimensionnel puis appliquer (1026) un seuillage sur chaque tranche (TRx, TRy) de la première image filtrée, pour localiser sur chaque tranche (TRX, TRy), une pluralité de zones correspondant respectivement à une pluralité de sections de torons (SECx, SECy), pour obtenir une deuxième image filtrée tridimensionnelle, Pour chaque zone correspondant à une section de toron, déterminer (1027), sur la deuxième image filtrée, une position de barycentre correspondant à la position du centre de ladite section de toron. Procédé (100) selon les revendications 8, 9 et 11 selon lequel un niveau de gris des voxels de l’image tridimensionnelle ne correspondant pas à de la résine est inférieur à un niveau de gris de référence, le procédé comprenant en outre: inverser le niveau de gris N associé à chaque voxel de l'image tridimensionnelle, et remplacer le niveau de gris de chaque voxel de l'image tridimensionnelle correspondant à de la résine par un niveau de gris minimal des voxels de l'image tridimensionnelle pour obtenir une image tridimensionnelle modifiée. Procédé (100) selon la revendication 2 ou l'une quelconque des revendications 3 à 12 en combinaison avec la revendication 2, selon lequel la détermination de la pluralité de contours de sections de torons dans chaque région localisée est réalisée à partir d'une méthode de segmentation par contours actifs ou par une méthode de segmentation par croissance de régions, ou une méthode de segmentation par ligne de partage des eaux. Procédé (100) selon l'une quelconque des revendications précédentes selon lequel la pièce est une pièce aéronautique. Dispositif (200) de contrôle non destructif d'une pièce en matériau composite par squelettisation de torons (Tx, Ty) de la pièce dans une image tridimensionnelle de ladite pièce, les torons étant associés à une même direction prédéterminée (X,Y), l'image tridimensionnelle comprenant des tranches (TRx, TRy) successives orthogonales à la direction prédéterminée, le dispositif étant configuré pour: Déterminer, pour chaque tranche (TRx, TRy) de l'image tridimensionnelle, des positions de centres (Cx, Cy) de sections (SECx, SECy) de torons (Tx, Ty), Pour au moins une tranche (TRx, TRy) parmi les tranches de l'image tridimensionnelle, corriger au moins une position de centre (Cx, Cy) de section (SECx, SECy) de toron parmi les positions de centres (Cx, Cy) de sections de torons déterminées pour ladite au moins une tranche (TRx, TRy) ladite correction comprenant une segmentation des sections de torons localement peu compactés dans l’image tridimensionnelle, Déterminer une pluralité de squelettes (Sx, Sy) de torons (Tx, Ty) en reliant, pour plusieurs tranches (Tx, Ty) de l’image tridimensionnelle, un centre (Cx, Cy) de section de toron d’une tranche courante de l’image tridimensionnelle à un centre de section (Cx, Cy) de toron d’une tranche de l’image tridimensionnelle parmi les tranches voisines de la tranche courante, Emettre une alerte lorsqu’une position d'au moins un squelette déterminé s'écarte par rapport à une position d'un squelette prédéterminée. Produit programme d'ordinateur comportant des instructions pour mettre en œuvre le procédé selon l’une des revendications 1 à 14 lorsque ce programme est exécuté par un processeur.