FR2982364A1 - Procede et dispositif d'estimation d'un taux de porosite d'un echantillon de materiau a partir d'au moins une image codee en niveaux de gris - Google Patents

Abstract

Le procédé selon l'invention comporte : - une étape (E40) d'évaluation d'un taux intermédiaire (T(S)) de l'échantillon pour chaque valeur (S) d'une pluralité de valeurs de seuil comprises entre deux valeurs limites déterminées (S1,S2), ce taux intermédiaire étant égal au ratio entre le nombre de pixels de ladite au moins une image ayant un niveau de gris borné par cette valeur de seuil et le nombre total de pixels de ladite au moins une image ; et - une étape d'estimation (E50-E60) du taux de porosité de l'échantillon en analysant les variations du taux intermédiaire en fonction des valeurs de seuil.

Description

Arrière-plan de l'invention L'invention se rapporte au domaine général des matériaux. Elle concerne plus particulièrement l'estimation du taux de porosité d'un matériau (taux surfacique ou volumique), tel que par exemple un matériau composite. De façon connue, la porosité d'un matériau caractérise la faculté de ce matériau à comporter des pores, c'est-à-dire des vides interstitiels connectés entre eux ou non. L'invention a ainsi une application privilégiée mais non limitative dans le domaine de l'aéronautique.

En effet, il est courant aujourd'hui d'utiliser des structures primaires en matériau composite dans la production d'aéronefs (ex. aubes de turboréacteur, etc.). Ces structures sont soumises à des contrôles qualité stricts au cours desquels, le taux de porosité volumique des matériaux composites est étroitement surveillé. La présence de porosités dans le matériau pouvant handicaper sa bonne tenue mécanique, il s'agit par ces contrôles de s'assurer que le taux de porosité volumique du matériau ne dépasse pas une valeur limite prédéfinie. Pour déterminer le taux de porosité d'un matériau composite, il est connu de recourir à une technique de dégradation ou de dissolution (ex. chimique par attaque acide ou par calcination) de la matrice du matériau. Conformément à cette technique, on réalise un relevé de masses d'un échantillon du matériau composite avant et après dissolution de la matrice. Ce relevé de masses permet, à partir de la connaissance des densités des fibres et de la matrice du matériau, de calculer aisément le taux de porosité volumique du matériau composite. Toutefois, celte technique présente un certain nombre d'inconvénients, Tout d'abord, elle dépend rottement de la précision du relevé ffinsique de la connaissance des fibres COrl 2,5 ,-)L1 LIC composants étrangers dans le matériau, tels que par exemple des inclusions métalliques ou de la fibre de verre, a une forte incidence sur le calcul du taux de porosité volumique. Enfin, cette technique est difficilement applicable de manière industrielle sur des matériaux métalliques et sur des matériaux composites à matrice céramique ou métallique. Le document de Y. Ledru et al. intitulé «Quantification 2-D et 3D de la porosité par analyse d'images dans les matériaux composites stratifiés aéronautiques », JNC 16, Toulouse, 2009, propose une technique d'estimation du taux de porosité d'un matériau composite non destructive, basée sur l'analyse d'images codées en niveaux de gris. Elle propose plus particulièrement d'isoler dans ces images les pixels correspondant aux porosités des pixels correspondant à la matière. Le taux de porosité volumique du matériau composite est alors déduit du nombre de pixels correspondant aux porosités ainsi isolés. Toutefois cette technique repose sur la fixation par un opérateur d'un seuil de niveau de gris pour différencier les pixels. L'intervention de l'opérateur rend ainsi l'analyse conduite dans ce document subjective et donc difficilement contrôlable et reproductible.

Objet et résumé de l'invention La présente invention permet de remédier notamment à cet inconvénient en proposant un procédé d'estimation d'un taux de porosité d'un échantillon de matériau à partir d'au moins une image codée en niveaux de gris représentant cet échantillon, ce procédé comportant : une étape d'évaluation d'un taux intermédiaire pour chaque valeur d'une pluralité de valeurs de seuil de niveaux de gris comprises entre deux valeur-7, limites déterminées., ce taux intermédiaire étant égal au ratio entre le nombre de piels de ladite au moins une image ayant un niveau de gris borné par cette voleur de seuil et le nombre total de CC-. bidte au mciin inc in ag ; et ilion en ars d'au moins une image codée en niveaux de gris représentant cet échantillon, ce dispositif comportant - des moyens pour évaluer un taux intermédiaire pour chaque valeur d'une pluralité de valeurs de seuil de niveaux de gris comprises entre deux valeurs limites déterminées, ce taux intermédiaire étant égal au ratio entre le nombre de pixels de ladite au moins une image ayant un niveau de gris borné par cette valeur de seuil et le nombre total de pixels de ladite au moins une image ; et des moyens pour estimer le taux de porosité de l'échantillon en analysant les variations du taux intermédiaire en fonction des valeurs de seuil. Par « pixels ayant un niveau de gris borné par une valeur de seuil », on entend au sens de l'invention les pixels qui ont des niveaux de gris respectifs situés d'un même côté de la valeur de seuil, c'est-à-dire soit bornés à droite par la valeur de seuil, soit bornés à gauche. Le choix de l'un ou l'autre des côtés dépend de la façon dont les porosités et/ou la matière sont représentées sur les images codées en niveaux de gris, i.e. selon si l'on attribue sur les images les niveaux de gris les plus faibles aux porosités ou à la matière.

Autrement dit, le taux intermédiaire évalué conformément à l'invention est égal - à la proportion de pixels des images ayant un niveau de gris inférieur à la valeur de seuil si sur les images les niveaux de gris représentant les porosités de l'échantillon sont inférieurs aux niveaux de gris représentant la matière ; ou à la proportion de pixels des images ayant un niveau de gris supérieur à la valeur de seuil si sur les images, les niveaux de gris représentant es porosités de l'échantillon sont supérieurs aux niveaux de gris I eprésentant la matière.

L'invention propose ainsi une méthode non destructive, simple taLiI. de pour rOterminer CCvc.,![..irnique ppliqu H ent à tout est très complexe d'utiliser des méthodes classiques telles que des méthodes de dissolution par voie acide. Le procédé selon l'invention peut être entièrement automatisé et s'affranchir de l'intervention d'un opérateur. Ainsi, non seulement les résultats obtenus sont reproductibles et précis, mais l'utilisation du procédé selon l'invention résulte également en un gain de temps non négligeable par rapport aux techniques de dissolution (gain supérieur à un facteur 3 pour un échantillon de dimensions 10mmx10mmx10mm, en incluant l'acquisition des images codées en niveaux de gris).

On notera par ailleurs que cette technique permet de traiter des échantillons de tailles plus importantes que les techniques de dissolution précédemment décrites. Préférentiellement, les images codées en niveaux de gris considérées sont des images représentant l'échantillon en trois dimensions. Ces images sont obtenues par exemple par tomographie. Les pixels sont alors des pixels en trois dimensions (chaque pixel associe par exemple un niveau de gris aux coordonnées spatiales d'un point de l'image), aussi connus sous le nom de voxels (contraction de « volumetric pixels »). On obtient ainsi directement, grâce au procédé d'estimation selon l'invention, un taux de porosité volumique. Il n'est de ce fait pas nécessaire pour obtenir le taux de porosité volumique de l'échantillon de procéder à une conversion d'un taux surfacique : l'estimation obtenue est par conséquent plus précise. La courbe représentant les variations du taux intermédiaire en fonction des valeurs de seuil est très caractéristique : elle présente en effet de part et d'autre d'un point de « rupture », deux parties sensiblement linéaires de pentes très distinctes. Or, les inventeurs ont avantageusement constaté (et vérifie à l'aide dune technique de dissolution) que le point de rupture reflétant le changement de pente de la courbe donne une estimation précise du taux de porosité volumique du itériai; , au cours de - on estime le taux de porosité de l'échantillon à partir de l'ordonnée de ce point. Par changement de pente significatif, on entend ici qui marque une rupture franche dans la pente de la courbe, c'est-à-dire que la pente de la courbe de part et d'autre de ce point est très différente (ou de façon équivalente, le rapport des pentes de part et d'autre de ce point est supérieur à un seuil prédéfini). Ainsi par exemple, pour détecter ce point, on s'intéresse à la dérivée seconde de la courbe : le point reflétant un changement de pente significatif de la courbe correspond en effet à un maximum de la dérivée seconde. En variante, on s'assure d'un changement significatif de pente en vérifiant que le rapport (ratio) des pentes de part et d'autre de ce point est supérieur à un seuil prédéfini.

Dans un mode de réalisation, les valeurs limites sont déterminées à partir des informations contenues dans au moins une dite image codée en niveaux de gris. Par exemple, ces valeurs limites peuvent être extraites automatiquement en observant les valeurs prises par les niveaux de gris dans l'image.

Dans un autre mode de réalisation, le procédé comprend en outre avant l'étape d'évaluation, une étape de classement des niveaux de gris de ladite au moins une image codée en niveaux de gris par ordre croissant ou décroissant. Cette étape facilite le traitement des images et le comptage des pixels de niveaux de gris identiques. Par ailleurs, il n'est pas utile grâce à cette étape d'analyser la totalité des pixels des images : on peut en effet se limiter à l'analyse des pixels ayant un niveau de gris compris entre les Jet< valeurs limites des valeurs de seuil. Dans un mode particulier de réalisation, les différentes étapes du procédé d'estimation du taux de porosité d'un échantillon de matériau sont d tuiminc ictions de [-)rncirarrunes d'urcliné3teurs. LLi 1-111\., ut( c nro(iiian'me instructions adaptées à la mise en oeuvre des étapes d'un procédé d'estimation tel que décrit ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur. D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres 20 moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en 25 question. On peut également envisager, dans d'autres modes de réalisation, que le procédé et/ou le dispositif d'estimation d'un taux de porosit.:: selon l'invention présentent en combinaison tout ou partie des caractéristiques orécitées 30 1P,A211t1(.)P1 la figure 1 représente, dans son environnement, un dispositif d'estimation d'un taux de porosité d'un matériau conforme à l'invention, dans un mode particulier de réalisation ; - la figure 2 représente, sous forme d'ordinogramme, les principales étapes d'un procédé d'estimation conforme à l'invention, dans un mode particulier de réalisation dans lequel il est mis en oeuvre par le dispositif de la figure 1; - la figure 3 représente schématiquement les variations de niveau de gris sur une image ; et la figure 4 donne un exemple de courbe représentant les variations du taux intermédiaire en fonction des niveaux de seuil. Description détaillée de l'invention La figure 1 représente, dans son environnement un dispositif 1 d'estimation d'un taux de porosité TE d'un échantillon E de matériau, conforme à l'invention, dans un mode de réalisation particulier de l'invention. Dans l'exemple envisagé ici, l'échantillon E est un échantillon de matériau composite à matrice organique, céramique ou métallique.

Toutefois, l'invention ne se limite pas aux matériaux composites et peut être appliquée à des échantillons de matériau quelconque, tels que par exemple à un échantillon de matériau homogène (ex. métallique, plastique, etc.). Conformément à l'invention, le dispositif d'estimation 1 est apte 25 à estimer le taux de porosité de l'échantillon E à partir d'une ou de plusieurs images II., I2,..., IN (N?.1, ex. N=1000) représentant l'échantillon E et codées en niveaux de gris. Dans l'exemple envisagé ici, le vide est représenté sur les imades II, IN par des niveaux de gris faibles » par rapport aux niveaux de gris représentant la matière : 30 autrement dit, sur les images II, IN, les niveaux de gris représentant- des porosités de 'échantilion E sont infE.,',rinut-s au: niveau>. de aria ru unta ut la matière. n end mitative, l'inventi, Le nombre N d'images considérées pour estimer le taux de porosité peut dépendre de divers paramètres, comme notamment de la présence importante ou non de porosités sur les images (ex. si une image représente un nombre important de porosités de l'échantillon de matériau, cette seule image peut s'avérer suffisante pour estimer de façon fiable le taux de porosité de l'échantillon conformément à l'invention), d'un compromis entre précision de l'estimation et ressources en termes de puissance calcul et/ou de mémoire requises pour l'estimation, etc. Dans le mode de réalisation décrit ici, les images représentent l'échantillon E en trois dimensions (i.e. sous la forme d'une pluralité de voxels, c'est-à-dire de pixels en trois dimensions), de sorte que le taux de porosité estimé par l'invention est un taux de porosité volumique. Les images Il, I2,...,IN sont obtenues ici par tomographie, par 15 exemple à rayons X. De façon connue en soi, les techniques de tomographie à rayons X permettent une reconstruction volumique de l'objet traité en utilisant des rayons X. La qualité des images tomographiques sera choisie suffisante pour permettre de différentier la matière du vide, ce qui ne pose en pratique aucune difficulté aujourd'hui 20 compte tenu des techniques utilisées en tomographie. Toutefois ces hypothèses ne sont pas limitatives. En effet, l'invention permet également d'estimer un taux de porosité à partir d'images représentant l'échantillon en deux dimensions (sous la forme d'une pluralité de pixels en deux dimensions). Le taux de 25 porosité obtenu sera alors un taux de porosité surfacique. Ce taux pourra être si nécessaire converti en taux de porosité volumique, de façon connue de l'homme du métier. A titre d'exemple, une technique de conversion est proposée dans le document de Y. Ledru et al. cité pré :écfeni ment. 30 Ains sens de l'invention, le Ur-me pr..e.1 désigne aussi bien en CIEL flFi t2T) BOLS Clill]CHsiOnS voxel fonction CIL tvp Cl!rages cod.,2,' considérées pour d'imagerie par résonance magnétique (IRM), dès lors que les images fournies sont de qualité suffisante pour distinguer le vide de la matière. Le dispositif d'estimation 1 dispose ici de l'architecture matérielle d'un ordinateur. Il comporte notamment un processeur 2, une mémoire vive 3 et une mémoire morte 4, connus en soi. La mémoire morte 4 du dispositif d'estimation 1 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 2 et sur lequel est enregistré un programme d'ordinateur conforme à l'invention, comportant des instructions pour l'exécution des étapes d'un procédé d'estimation d'un taux de porosité selon l'invention décrites maintenant en référence à la figure 2. La figure 2 représente, sous forme d'ordinogramme, les principales étapes d'un procédé d'estimation d'un taux de porosité conforme à l'invention, dans un mode particulier de réalisation.

Dans l'exemple envisagé ici, les N images H, I2,..., IN de l'échantillon E codées en niveaux de gris et représentant l'échantillon E en trois dimensions sont transmises au dispositif d'estimation 1 (étape E10). Ces images sont stockées dans la mémoire morte 4 du dispositif d'estimation 1 sous forme de fichiers F1,...,FN, chaque fichier Fn, n=1,...,N, associant à chaque pixel 3D (voxel) de l'image In (identifié par ses coordonnées) un niveau de gris codé sur k bits (par exemple k=16 bits, ce qui permet de coder 65536 niveaux de gris). Par souci de simplification dans la suite de la description, on fera référence aux pixels de l'image In, qui sont en fait dans l'exemple envisagé ici, des voxels.

Dans le mode de réalisation décrit ici, le dispositif d'estimation 1 crée un fichier F à partir des informations contenues dans les fichiers dans lequel il répertorie les niveaux de gris de tous les pixels representés sur les images I1,...,IN. En outre, dans ce fichier, les niveaux de gris sont classés par ordre croissant à l'aide d'un algorithme de tri connu en ni afin de faciliter leur tri.-iiten-k-,nt. entendu, cette hynothèseImitative, et un autre type (Jc. (7135 correspondant au fichier sélectionné (une ligne correspond à des pixels d'ordonnées constantes dans le fichier sélectionné). Préférentiellement, le fichier sélectionné par le dispositif d'estimation 1 comporte au moins une ligne traversant une porosité. Pour identifier un tel fichier, on peut par exemple considérer, en plus de l'échantillon dont on cherche à déterminer le taux de porosité, un échantillon étalon comportant une porosité connue, et sélectionner un fichier comportant une ligne traversant cette porosité connue de l'échantillon étalon.

En variante, le dispositif 1 d'estimation sélectionne un fichier au hasard ou un fichier prédéterminé, par exemple Fi, ou encore un fichier identifié par un opérateur parmi les fichiers F1,...,FN comme comportant au moins une ligne traversant une porosité. La figure 3 illustre un exemple de courbe représentant les variations de niveaux de gris sur une ligne de l'image ainsi sélectionnée (x représentant l'abscisse des pixels). A partir de l'analyse de ces variations, le dispositif 1 extrait deux valeurs limites de niveaux de gris, Si et S2 (étape E30). Ces deux valeurs limites définissent un intervalle de niveaux de gris dans lequel se situe le seuil de niveau de gris qui correspond à la frontière entre le vide et la matière, autrement dit ici, en deçà duquel un pixel est représentatif d'une porosité de l'échantillon E (i.e. vide) et au-dessus duquel un pixel correspond à de la matière. En variante, si sur les images Il, 12,..., IN, les niveaux de gris 25 représentant les porosités sont supérieurs aux niveaux de gris représentant la matière, le seuil de niveau de gris correspondant à la frontière entre le vide et la matière correspond au seuil au-dessus duquel un est représentatif d'une porosité éc l'échantillon E et en deçà correspond à de là mat.ière. 30 Plus précisément il chc.i)st_ici pour Sitaleur minimale prise ne -(360 et pour 'Ç)0O sur cette Par ailleurs, en variante, on peut extraire 51 et 52 des variations de niveaux de gris correspondant à plusieurs images, afin notamment de confirmer les valeurs sélectionnées sur une image. Dans une autre variante encore de réalisation de l'invention, le dispositif d'estimation 1 détermine comme valeurs limites de niveaux de gris les deux valeurs extrêmes S1=0 et 52=2k de niveaux de gris pouvant être prises par les pixels de l'image (soit ici 52=65536, k étant égal à 16 bits). Les valeurs Si et S2 ainsi fixées constituent des valeurs de seuil limites au sens de l'invention. Le dispositif 1 initialise alors une valeur de seuil S à Si. et évalue un taux intermédiaire T(S) de l'échantillon E de matériau (étape E40) à l'aide de l'équation suivante T(S)=AnfiAtot où Anf désigne le nombre de pixels représentant l'échantillon E répertoriés dans le fichier F (autrement dit de l'ensemble des images I1,..,IN) ayant un niveau de gris inférieur à la valeur de seuil S et Atot désigne le nombre total de pixels répertoriés dans le fichier F représentant l'échantillon E (i.e. somme des nombres de pixels de chaque image représentant l'échantillon E). En variante, si sur les images Il, I2,..., IN, les niveaux de gris représentant les porosités sont supérieurs aux niveaux de gris représentant la matière, le nombre Ah-if désigne le nombre de pixels répertoriés dans le fichier F ayant un niveau de gris supérieur à la valeur de seuil S. On notera que les niveaux de gris des pixels étant classés par ordre croissant, l'identification des pixels ayant un niveau de gris inférieur au seuil S est facilitée. Par ailleurs, de façon avantageuse, il n'est pas nécessaire d'analyser les niveaux de gris de tous les pixels stockés dans le fichier F pour déterminer les nombres A.. et A. l_dy,_,Ieur de tau j nte éd ai e r(S) ainsi obtenue est ;ri C l iernDI ; dispcsI 1 I ii;sociation avec la intermédiaire pour cette nouvelle valeur de seuil, etc., jusqu'à la valeur S=S2. Dans le mode de réalisation décrit ici, l'incrément incr est constant sur la plage de valeurs [S1,S2]. En variante, cet incrément peut dépendre des valeurs de seuil S considérées. Ainsi, par exemple, en référence à la figure 3, pour des valeurs de seuil S comprises entre 7000 et 9000 on pourra choisir un incrément relativement grossier (ex. de l'ordre de 500), peu de points étant compris entre ces deux valeurs sur l'image considérée. En revanche, pour des valeurs de seuil S supérieures à 9000, on choisira un incrément plus petit, afin d'identifier plus précisément le niveau de gris représentatif de la frontière réelle entre vide et matière (ex. incrément de l'ordre de 100 ou 50). Le taux de porosité TE de l'échantillon E est ensuite obtenu en analysant les variations du taux intermédiaire T(S) en fonction des valeurs de seuil S. Comme mentionné précédemment, la courbe C représentant les variations du taux intermédiaire en fonction des valeurs de seuil est très caractéristique. Un exemple d'une telle courbe est donné sur la figure 4. Comme cela est apparent sur cette figure, la courbe C présente en effet de part et d'autre d'un point de « rupture » P, deux portions de courbe Cl et C2 sensiblement linéaires et de pentes très distinctes. Or, les inventeurs ont avantageusement constaté que le point de rupture P reflétant ce changement de pente significatif de la courbe C donne une estimation précise du taux de porosité volumique TE du matériau.

Pour estimer le taux de porosité volumique TE, le dispositif 1 identifie donc dans un premier temps le point de rupture P (étape E50). Par définition, le point P correspond à un maximum de la dérivée seconde de la courbe. Par conséquent, dans le mode de réalisation ciecrit ici, le Lfispositif 1 evalue pour cheque valeur de seuil S, la dérivée seconde du ai.rs': intermédiaire T!H, é Je moyens connus en soi, tels que par Les valeurs de la dérivée seconde T"(S) sont ensuite stockées dans la mémoire vive 3 du dispositif 1 en association avec les valeurs de seuil S correspondantes. Puis, le dispositif 1 recherche parmi les valeurs stockées la valeur maximale de la dérivée seconde, et extrait la valeur de seuil correspondante notée Smax. A partir de cette valeur de seuil Smax, le dispositif 1 obtient la valeur du taux intermédiaire T(Smax), stockée au cours de l'étape E40 dans la mémoire vive 3. Autrement dit, le point de rupture P a pour coordonnées Smax et 10 T(Smax). Le dispositif 1 estime alors le taux de porosité volumique TE de l'échantillon E à partir de l'ordonnée du point (étape E60), soit ici TE= T(Smax). En variante, si les images sont des représentations en deux 15 dimensions de l'échantillon E, le taux intermédiaire T(Smax) est un taux surfacique. Le taux de porosité volumique de l'échantillon E peut alors être obtenu par conversion de ce taux surfacique en un taux volumique, de façon connue en soi. Dans le mode de réalisation décrit ici, le point de rupture est 20 identifié en recherchant le maximum de la dérivée seconde de la courbe C, autrement dit du taux intermédiaire. Dans un autre mode de réalisation, les portions de courbe Cl et C2 étant quasiment linéaires, le point de rupture est obtenu en approximant chaque portion de courbe par une portion de droite à l'aide d'une régression linéaire, puis en recherchant 25 l'intersection des deux portions de droite ainsi obtenues. Les techniques d'application d'une régression linéaire à un ensemble de points et de recherche d'une intersection de deux droites étant connues en soi, elles ne seront pas ck±tailiées davantage ici. Le noint d'intersection ainsi obtenu est ensuite appro.inié par un point de la courbe C (i.e. on recherche le point 30 de la courbe C le plus proche du point d'intersection obtenu), l'ordonnée de ce oint constituant une eqtimation ( taux de poro;ité volumique. On notera que 1;.prochc v osce da 11'; 'H VU-1111:D7', l'invention (à titre indicatif, une précision d'estimation de l'ordre de 0.1% a été obtenue lors des tests réalisés).

Claims (4)

1. REVENDICATIONS1. Procédé d'estimation d'un taux de porosité d'un échantillon de matériau à partir d'au moins une image (I1,...,IN) codée en niveaux de gris représentant cet échantillon, ledit procédé comportant : une étape (E40) d'évaluation d'un taux intermédiaire (T(S)) pour chaque valeur (S) d'une pluralité de valeurs de seuil de niveaux de gris comprises entre deux valeurs limites déterminées (S1,S2), ce taux intermédiaire étant égal au ratio entre le nombre de pixels de ladite au moins une image ayant un niveau de gris borné par cette valeur de seuil et le nombre total de pixels de ladite au moins une image ; - une étape d'estimation (E50-E60) du taux de porosité (TE) de l'échantillon en analysant les variations du taux intermédiaire en fonction des valeurs de seuil.
2. 2. Procédé d'estimation selon la revendication 1 dans lequel au cours de l'étape d'estimation du taux de porosité: - on identifie (E50), sur une courbe représentant les variations du taux intermédiaire en fonction des valeurs de seuil, un point (P) reflétant un changement de pente significatif de la courbe ; et - on estime (E60) le taux de porosité de l'échantillon à partir de l'ordonnée de ce point.
3. 3. Procédé d'estimation selon la revendication 2, dans lequel le point (P) reflétant un changement de pente significatif de la courbe correspond à un maximum de la dérivée seconde de la courbe.
4. 4. Procédé d'estimation selon l'une quelconque des revendications I 0 5, dans lequel ladite au moins une image (11,...,1 codée en niveau>: de gus est Lune image représentent l'échantillon en trois dimensions, et' taux de oornsfté .es u n ta u e usité \y/Glu-1 u6. Procédé d'estimation selon l'une quelconque des revendications 1 à 5 dans lequel les valeurs limites (S1,S2) sont déterminées à partir des informations contenues dans au moins une dite image codée en niveaux de gris. 7. Procédé d'estimation selon l'une quelconque des revendications 1 à 6 dans lequel le matériau est un matériau composite. 10 8. Procédé d'estimation selon l'une quelconque des revendications 1 à 7 comprenant en outre, avant l'étape d'évaluation, d'une étape (E20) de classement des niveaux de gris de ladite au moins une image codée en niveaux de gris par ordre croissant ou décroissant. 15 9. Programme d'ordinateur comportant des instructions pour l'exécution des étapes du procédé d'estimation selon l'une quelconque des revendications 1 à 8 lorsque ledit programme est exécuté par un ordinateur. 20 10. Dispositif (1) d'estimation d'un taux de porosité d'un échantillon de matériau à partir d'au moins une image codée en niveaux de gris représentant cet échantillon, ledit dispositif comportant : des moyens pour évaluer un taux intermédiaire (T(S)) pour chaque valeur (S) d'une pluralité de valeurs de seuil de niveaux de gris 25 comprises entre deux valeurs limites (S1,S2) déterminées, ce taux intermédiaire étant égal au ratio entre le nombre de pixels de ladite au moins une image ayant un niveau de gris borné par cette valeur de seuil et le nombre total de pixels de ladite au moins une image S moyens pour estimer le taux de porosité (TE) de l'échantillon en 30 analysant les variations du taux intermédiaire en fonction des valeurs de uil