Brevet FR3033068A1 — publié le 20160826

La présente invention concerne le domaine de la reconstruction tridimensionnelle d'un objet à partir de mesures externes à l'objet, comme par exemple dans le cadre de la tomographie électronique. Le terme de tomographie désigne l'ensemble des méthodes de reconstruction exactes, ou le plus souvent approchées, de l'intérieur d'un objet à partir de mesures effectuées à l'extérieur de cet objet. Ces mesures, appelées projections par les spécialistes dans le domaine, sont liées à l'absorption par l'objet considéré d'une énergie émise par une source. Souvent, ces mesures sont effectuées par l'intermédiaire d'un capteur, et réalisées sous différents angles de mesure autour de l'objet considéré. Les sources utilisées peuvent être de plusieurs types : rayons X, rayons Gamma, ou encore électrons pour la tomographie électronique. Dans le domaine du raffinage pétrochimique, la reconstruction tridimensionnelle est utilisée pour caractériser les supports de catalyseurs. En effet, les mesures morphologiques de particules permettent de remonter aux spécificités physico-chimiques des matériaux étudiés. A l'aide par exemple d'un MET (Microscopie Electronique en Transmission), on réalise une série de projections, à intervalles angulaires réguliers. Puis une méthode de reconstruction tomographique est mise en oeuvre, par exemple par une approche de type problème inverse afin d'obtenir une image du support considéré. La reconstruction tridimensionnelle peut aussi être utilisée dans le domaine de l'exploration et de l'exploitation pétrolière, pour la caractérisation des matériaux poreux. En effet, cette technique permet d'avoir accès au réseau poreux interne d'un échantillon de roche, et par conséquent de pouvoir caractériser très précisément celui-ci. Ainsi, cette technique peut permettre d'accéder à des informations très importantes pour la caractérisation d'une roche réservoir, telles que sa porosité (densité de pores dans l'échantillon), ou encore sa perméabilité (connectivité des pores entre eux, potentiellement préférentielle dans une direction). Mais la reconstruction tridimensionnelle d'un objet peut concerner des domaines très divers : la caractérisation des matériaux en général, le contrôle non destructif, ou encore l'imagerie médicale. L'objet reconstruit, ou image de l'objet, est alors présenté sous la forme d'une maquette, généralement représentée sur ordinateur, via des niveaux d'intensité, chaque niveau d'intensité reflétant un constituant de l'objet. Par exemple dans le cas d'un support de catalyseur, des niveaux d'intensité pourront représenter le vide et d'autres la matière du support. Pour les roches, chaque constituant chimique peut être représenté par des niveaux d'intensité. Les niveaux d'intensité sont classiquement représentés en utilisant une palette de couleur en noir et blanc, conduisant à en une représentation en dégradés de gris. Le spécialiste peut alors chercher à simplifier cette information, afin de permettre des caractérisations de type morphologique ou géométrique (volume, aire, surface spécifique, longueur). On entend par simplification, une réduction, sans 3033068 2 intervention d'un opérateur, du nombre de niveaux d'intensité représentant l'objet reconstruit. Typiquement, un objet reconstruit est représenté par plusieurs milliers de niveaux d'intensité, que l'on souhaite simplifier en deux ou trois niveaux ou classes. Par exemple, dans le cas de l'estimation du volume poreux d'une roche, on pourra avantageusement simplifier l'image 5 reconstruite de l'échantillon de roche considéré en deux classes, l'une correspondant au vide et l'autre correspondant à tout constituant de la roche. Une étape de simplification peut être aussi mise en oeuvre dans le cas où l'image de l'objet comporte des intensités non entières, que l'on va regrouper en un nombre raisonnable de classes d'intensités entières. Dans le domaine technique considéré, on parle de segmentation pour désigner cette étape de simplification. La segmentation 10 consiste donc à déterminer des seuils d'intensité, délimitant des classes d'intensité représentatives d'une propriété que l'on souhaite mettre en évidence. A noter qu'une étape de segmentation peut être aussi mise en oeuvre sur une image numérique d'un objet : la projection est alors calculée numériquement, à partir de l'image, et non obtenue par des mesures expérimentales. 15 État de la technique Les documents suivants seront cités au cours de la description : 20 [Batenburg et al. 2007] K. J. Batenburg et J. Sijbers. "Optimal threshold selection for tomogram segmentation by reprojection of the reconstructed image", in: Computer Analysis of Images and Patterns, in: Lecture Notes in Computer Science, vol. 4673, Springer, Berlin/Heidelberg, 563-570, 2007. 25 [Batenburg et al. 2009] K.J. Batenburg et J. Sijbers. "Adaptive thresholding of tomograms by projection distance minimization". Pattern Recognition, 42(10), 2297-2305, 2009. [Doan D. H. et al., 2013]: Hong Doan Dinh, Pierre Delage, Jean-François Nauroy, Anh-Minh Tang, Youssef Souhail, Microstructural characterization of a Canadian oil sand, Canadian Geotechnical Journal. 01/2013; 49(10). 30 [Gilbert 1972] Gilbert P., Iterative methods for the three-dimensional reconstruction of an object from projections. Journal of Theorical Biology, 36(1) :105-117. [Huang et al. 1995] L.K. Huang et M.J. Wang, Image thresholding by minimizing the measure of fuzziness, Pattern Recognition, 41-51, 1995. [Kapur et al. 1985] J.N. Kapur, P.K. Sahoo, A.K.C. Wong. "A new method for gray-level picture 35 thresholding using the entropy of the histogram", Computer Vision, Graphics, and Image Processing, 29(3), 273-285, 1985. 3033068 3 [Lei et al. 2008] Lei He, Zhigang Peng, Bryan Everding, Xun Wang, Chia Y. Han, Kenneth L. Weiss, William G. Wee, A comparative study of deformable contour methods on medical image segmentation, Image and Vision Computing, Volume 26, Issue 2, 1 February 2008, Pages 141163. 5 [Mesanovic et al. 2013] Nihad Mesanovis, Haris Huseinagic, Samir Kamenjakovic, Automatic Region Based Segmentation and Analysis of Lung Volumes from CT Images, International Journal of Computer Science and Technology, vol. 4, Issue special, 2013. [Moreaud et al. 2008] M. Moreaud, B. Celse, F. Tihay, Analysis of the accessibility of macroporous alumino-silicate using 3D-TEM images. Material Science and Technology 2008 conference, 10 Pittsburgh USA. [Otsu 1979] N. Otsu. "A threshold selection method from gray level histograms". IEEE Trans. Systems, Man and Cybernetics 9, 62-66, 1979. [Stawiaski, 2008] Jean Stawiaski, Morphologie mathématique et graphes : application à la segmentation interactive d'images médicales, thèse, Ecole des Mines de Paris, 2008. 15 [Tran et al. 2013] VD Tran, M Moreaud, É Thiébaut, L Denis, JM Becker, Inverse problem approach for the alignment of electron tomographic series, Oil & Gas Science and TechnologyRevue d'IFP Energies nouvelles 69 (2), 279-291. [Tran 2013] V.D. Tran, Reconsutrction et segmentation d'image 3D de tomographie électronique par approche « problème inverse », Thèse, Université de Saint Etienne, 2013. 20 [Xu et al. 1996] L. Xu et M.I. Jordan, "On Convergence Properties of the EM Algorithm for Gaussian Mixture", Neural Computation, 8, 129-151, 1996. La simplification, ou segmentation, des données volumiques sans aucune intervention manuelle s'effectue classiquement avec des approches globales utilisant un histogramme (non 25 cumulé) de l'image reconstruite, comme par exemple : la maximisation de la variance interclasse [Otsu 1979], la maximisation de l'entropie [Kapur et al. 1985], l'ajustement des distributions gaussiennes multiples à l'histogramme [Xu et al. 1996], les méthodes des nuées dynamiques (ou k-moyennes) [Huang et Wang 1995]. Il est également possible d'utiliser des approches locales de type contours actifs [Lei et al. 2008], croissance de région [Mesanovic et al. 2013], ou encore issue 30 de la morphologique mathématique utilisant un opérateur de ligne de partage des eaux [Stawiaski, 2007]. Ces approches peuvent être totalement automatiques, semi-automatiques, ou guidées par l'utilisateur. Leur inconvénient majeur réside dans le fait que ces méthodes utilisent en entrée l'objet reconstruit, qui est une approximation du vrai objet, comportant des erreurs. Ces méthodes ne sont donc pas directement reliées aux informations initiales sur l'objet étudié que sont les 35 projections. 3033068 4 Pour pallier ce problème, une méthode spécifique pour la tomographie, basée sur une optimisation globale, a été récemment proposée par [Batenburg et al. 2007], puis par [Tran 2013]. Le principe consiste à déterminer les valeurs des seuils d'intensité minimisant, via une méthode d'optimisation itérative, la différence entre les projections réelles (c'est-à-dire les projections 5 réellement enregistrées par les outils de mesure de projections) et celles calculées à partir de l'objet reconstruit, puis segmenté. Cette méthode est connue par le spécialiste comme la méthode de minimisation de la distance aux projections. La Figure 1 présente le principe de cette méthode : à partir de projections réelles P, on reconstruit une image I de l'objet étudié comportant NI niveaux intensités, par exemple par reconstruction tomographique ; puis à partir de valeurs de seuils 10 d'intensité s prédéfinies, on calcule, par une méthode de segmentation S, une image segmentée IS de cette image reconstruite I, comportant NS classes d'intensité. On calcule ensuite des projections P' à partir de cette image segmentée IS. Puis, on réalise une comparaison PCP' entre les projections réelles P et des projections calculées P' sur l'image segmentée IS ; on met alors à jour les valeurs de seuils d'intensité s, jusqu'à ce que les projections réelles P et les projections 15 calculées P' soient le plus proches possibles, c'est-à-dire jusqu'à ce que leur différence soit inférieure à un seuil prédéfini. Cette méthode, si elle permet de travailler sur les informations initiales sur l'objet étudié, est très consommatrice en temps de calcul. En effet, pour chaque itération du processus d'optimisation, c'est-à-dire à chaque mise à jour s+As des valeurs de seuils d'intensité, il faut 20 recalculer les projections P' de l'objet reconstruit segmenté. La présente invention vise à réduire les temps de calculs pour obtenir un objet simplifié par la méthode de minimisation de la distance aux projections. Pour ce faire, la présente invention repose sur un pré-calcul des projections calculées P', réalisé une seule fois, permettant d'éviter le 25 calcul des projections calculées P' à chaque itération d'un processus d'optimisation. Le procédé selon l'invention 30 Ainsi, la présente invention concerne un procédé pour la segmentation d'une image d'un objet représentée par un nombre fini de niveaux d'intensité, ladite image ayant été reconstruite à partir de mesures externes à l'objet et acquises selon différents angles pris par rapport à une référence. A partir d'un nombre prédéfini de seuils d'intensité de ladite image, le procédé comporte les étapes suivantes pour au moins un desdits angles : 3033068 5 a) on détermine, pour ledit angle, une distribution représentative de la répartition desdits niveaux d'intensités de ladite image selon une direction prédéterminée dans ladite image ; b) on détermine les valeurs desdits seuils d'intensité en minimisant une fonction 5 objectif estimant la différence entre ladite mesure pour ledit angle et une mesure simulée pour ledit même angle, ladite mesure simulée étant fonction desdites valeurs desdits seuils et de ladite distribution ; c) on segmente ladite image par regroupement desdits niveaux d'intensité en classes d'intensité délimitées par lesdites valeurs desdits seuils. 10 Avantageusement, ladite direction prédéterminée peut être une droite parallèle à l'un des axes horizontaux de ladite image. Préférentiellement, ladite distribution peut être un histogramme cumulé desdits niveaux d'intensités rencontrés dans ladite image. Selon un mode de réalisation de l'invention, lesdites mesures peuvent avoir été acquises 15 dans le cadre d'une géométrie d'acquisition parallèle et ledit histogramme cumulé peut être déterminé sur une rotation de ladite image selon ledit angle. Selon un mode de mise en oeuvre de la présente invention, on peut calculer ladite mesure simulée P,4 pour un desdits angles A selon une formule du type : PA - f (TA(S1), ..., TA(Sn), ..., TA(SNS)) 20 où f est une combinaison linéaire, NS est ledit nombre prédéfini desdits seuils d'intensité, et TA(sn) est l'ensemble des éléments dudit histogramme cumulé représentatifs du nombre desdits niveaux d'intensité de ladite image inférieurs à une valeur de seuil sn, ledit histogramme étant calculé pour ledit angle A. Selon un mode de mise en oeuvre de la présente invention, on peut déterminer une 25 intensité représentative pour chacune desdites classes délimitées par lesdites valeurs desdits seuils. Selon un mode de mise en oeuvre de la présente invention, ladite intensité représentative de la dite classe peut être la valeur d'un desdits seuils délimitant ladite classe, une moyenne desdites valeurs desdits seuils délimitant ladite classe, ou toute combinaison 30 linéaire desdites valeurs desdits seuils délimitant ladite classe. Selon un mode de réalisation de l'invention, on peut appliquer les étapes a) et b) pour chacun desdits angles de mesure et on peut déterminer des valeurs de seuils représentatives de l'ensemble desdits angles en choisissant soit la médiane, soit une combinaison desdites valeurs desdits seuils déterminés pour chacun desdits angles de mesure. 3033068 6 Selon un autre mode de réalisation de l'invention, on peut déterminer des valeurs de seuils représentatives de l'ensemble desdits angles de mesure en minimisant une fonction objectif estimant les différences entre lesdites mesures et lesdites mesures simulées pour l'ensemble desdits angles de mesure. 5 De plus, l'invention concerne une utilisation du procédé selon l'invention pour la caractérisation des supports de catalyseurs dans le domaine du raffinage pétrochimique. En outre, l'invention concerne une utilisation du procédé selon l'invention pour la caractérisation des milieux poreux dans le domaine de l'exploration et de l'exploitation pétrolière. 10 Par ailleurs, l'invention concerne une utilisation du procédé selon l'invention pour la caractérisation de matériaux, pour le contrôle non destructif d'un objet, ou pour l'imagerie médicale. En outre, l'invention concerne un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou 15 exécutable par un processeur, comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon la description ci-dessus, lorsque ledit programme est exécuté sur un ordinateur. Présentation succincte des Figures 20 D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. La Figure 1 présente un schéma du principe de la méthode segmentation par minimisation de la distance aux projections. 25 Les Figures 2A à 2F présentent le principe, selon l'invention, du calcul d'une distribution représentative de la répartition des niveaux d'intensité dans une image selon une ligne horizontale, ladite distribution étant un histogramme cumulé. Les Figures 3A à 3C présentent des résultats obtenus par le procédé selon l'invention. 30 Description détaillée du procédé Les définitions suivantes sont utilisées au cours de la description de l'invention : 3033068 7 - projections : dans le domaine de la reconstruction tridimensionnelle, il s'agit de mesures déportées à l'extérieur de l'objet et acquises selon différents angles de mesures pris par rapport à une référence. Ces mesures résultent de l'émission d'une énergie émise par une 5 source, partiellement absorbée par l'objet étudié, l'énergie restante étant enregistrée par un capteur. Les sources utilisées peuvent être de plusieurs types : rayons X, rayons Gamma, ou encore électrons pour la tomographie électronique. Ces mesures sont généralement répétées pour différents angles de mesures pris par rapport à une référence. En tomographie par rayons X, il est courant d'utiliser 720 projections, à 360° autour de l'objet, par incrément de 0.5°. En 10 tomographie électronique, les projections sont réalisées de -70° et +70°, par incrément de 1°. - reconstruction tridimensionnelle: il s'agit d'une technique permettant de construire l'image volumique d'un objet à partir de projections de cet objet réalisées sous différents angles de mesure. La reconstruction conduit à une image représentée par un nombre fini de niveaux d'intensité, chaque niveau correspondant à un constituant de l'objet. En général, un constituant 15 est représenté par plusieurs niveaux d'intensité. Plus d'informations concernant la reconstruction tridimensionnelle peuvent être trouvées dans [Doan D. H. et al., 2013]. - géométrie d'acquisition parallèle/non parallèle: il s'agit de la géométrie de rayonnement de la source d'un dispositif d'acquisition de projections. On parle de géométrie non parallèle ou encore conique lorsque la source émet à partir d'un point et que les rayons traversant l'objet ne 20 sont pas parallèles. On parle de géométrie parallèle lorsque la source émet des rayons tous parallèles à partir d'une surface finie et que ces rayons traversent l'objet selon la même direction ; c'est le cas par exemple de la tomographie électronique. - classe d'intensité : il s'agit d'un regroupement de plusieurs niveaux d'intensité délimités par au moins un seuil et caractérisés par une intensité représentative des valeurs d'intensité de la 25 classe. Par exemple, une classe peut regrouper les valeurs d'intensité inférieures à un seuil donné, les valeurs d'intensité supérieures à un autre seuil, ou les valeurs comprises entre deux seuils. L'intensité représentative d'une classe donnée peut correspondre à la moyenne des valeurs d'intensité de la classe considérée, à la médiane, à l'un des seuils caractérisant la classe, ou encore être définie par tout autre calcul. 30 - segmentation : la segmentation d'une image représentée par un nombre fini de niveaux d'intensité consiste à regrouper les niveaux d'intensité en question en classes d'intensité. L'objet de la présente invention est un procédé de segmentation de l'image d'un objet représentée par un nombre fini de niveaux d'intensité. Selon l'invention, l'image considérée a été 35 reconstruite à partir de mesures externes à l'objet et acquises selon différents angles pris par rapport à une référence, qui peut être un axe vertical. Par la suite et conformément à l'usage dans le domaine, ces mesures sont appelées projections. L'objet de l'invention vise à regrouper les 3033068 8 niveaux d'intensité en un nombre prédéfini de classes d'intensité, et ce, en déterminant, de manière automatique, rapide, et en tenant compte des projections, les valeurs des seuils d'intensité délimitant les classes en question. 5 Préalablement à sa mise en oeuvre, le procédé selon l'invention peut requérir la réalisation des deux étapes suivantes : 1. acquisition des projections 10 Il s'agit ici d'acquérir les données d'entrée en vue de la reconstruction tridimensionnelle d'un objet. Le spécialiste a une parfaite connaissance des techniques d'acquisition de projections associées à un objet. Dans le cadre de la reconstruction tridimensionnelle par tomographie électronique, on peut obtenir ce type de projections en utilisant un Microscope Electronique en Transmission (MET en 15 français, TEM en anglais pour Transmission Electron Microscopy). Le MET comporte un canon à électrons qui envoie un faisceau parallèle d'électrons de haute énergie sur un objet disposé sur un porte-objet. L'objet disperse alors les électrons qui le traverse ; le faisceau sortant est ensuite refocalisé, puis enregistré par un capteur. Le signal enregistré est appelé projection ; il sera noté P dans la suite. Cette opération est répétée en inclinant le porte-objet du microscope selon différents 20 angles par rapport au faisceau d'électrons (typiquement de -70° à +70° avec un incrément de 1 °). On obtient ainsi des projections PA pour NA d'angles de mesures A. Ainsi, les projections contiennent l'information volumique sur un objet, projetée sur un plan. Le MET permet d'obtenir des informations sur un objet à l'échelle de quelques nanomètres à l'angstrôm. Classiquement, les projections peuvent être également acquises par rayons X, rayons 25 gamma ou analogues. 2. reconstruction tridimensionnelle Il s'agit ici, à partir des projections PA acquises à l'étape précédente pour différents angles de 30 mesure A, de reconstruire l'image tridimensionnelle, c'est-à-dire volumique, d'un objet. Le spécialiste a une parfaite connaissance de techniques de reconstruction tridimensionnelle d'un objet. La reconstruction peut être réalisée par plans de coupes successifs, puis le volume reconstruit final est obtenu par superposition des coupes reconstruites. Une méthode particulière de reconstruction tomographique est par exemple décrite dans [Tran et al. 2013]. On peut 3033068 9 également citer la méthode SIRT [Gilbert, 1972], couramment utilisée dans le domaine. La résolution de l'image de l'objet obtenue par reconstruction tridimensionnelle dépend notamment de la qualité de l'échantillonnage en angles de mesure (gamme d'angles de mesure considérée, et incrément entre les mesures). Ainsi, plus la couverture angulaire est complète et plus l'incrément 5 angulaire est petit, plus l'image de l'objet reconstruit est précise, c'est-à-dire plus l'image de l'objet est proche de l'objet réel. Dans le cadre de projections obtenues par MET, il est nécessaire au préalable d'appliquer une correction aux projections, en raison d'un déplacement de l'objet dans le porte-objet lors des mesures angulaires successives. 10 A partir d'un nombre de seuils d'intensité prédéfini par le spécialiste, la présente invention comporte au moins les étapes suivantes pour un angle de mesure donné : 1. Calcul d'une distribution représentative de la répartition des niveaux d'intensité : au cours de cette étape, on détermine, pour ledit angle, une distribution représentative de la 15 répartition desdits niveaux d'intensités de ladite image selon une direction prédéterminée dans ladite image ; 2. Estimation des valeurs des seuils d'intensité : au cours de cette étape, on détermine la valeur des seuils d'intensité en minimisant une fonction objectif estimant la différence entre ladite mesure pour ledit angle et une mesure simulée pour ledit même angle, ladite mesure 20 simulée étant fonction de la valeur desdits seuils et de ladite distribution ; 3. Segmentation de l'image : au cours de cette étape, on segmente ladite image par regroupement desdits niveaux d'intensité en classes délimitées par lesdites valeurs de seuil. 25 A des fins de simplification, les principales étapes du procédé selon l'invention sont présentées dans le cas de l'image d'un objet bidimensionnel, mais le même principe peut être aisément mis en oeuvre en trois dimensions par un spécialiste. Ainsi, par la suite, on considère une image I d'un objet en deux dimensions, de taille W x H, constituée de NI niveaux d'intensité. Un élément I(i,j) (ou maille) de l'image I est un élément situé au croisement de la ième colonne de I et 30 de la jème ligne de I, avec i variant entre 1 et W et j variant entre 1 et H. Un élément I(i,j) de l'image I contient une valeur d'intensité donnée. L'image I a été obtenue par reconstruction tridimensionnelle à partir de projections PA de l'objet réalisées sous un nombre NA d'angles de mesure A différents. Par ailleurs, les principales étapes du procédé selon l'invention sont présentées ci-après 35 dans le cas de la segmentation d'une image I en NS+1 classes d'intensité, c'est-à-dire dans le cas 3033068 10 général où l'on cherche à déterminer un nombre NS de valeurs de seuils d'intensité, le nombre NS de seuils étant un entier supérieur ou égal à 1, prédéfini par le spécialiste. Les principales étapes de la présente invention sont détaillées ci-après. 5 1) Calcul d'une distribution représentative de la répartition des niveaux d'intensité Au cours de cette étape, il s'agit de calculer une distribution représentative de la répartition des niveaux d'intensité rencontrés dans l'image I selon une direction prédéterminée. 10 Selon un mode préféré de réalisation de la présente invention, ladite distribution est l'histogramme cumulé T des niveaux d'intensités rencontrés dans l'image I selon une direction prédéterminée. Selon un mode de réalisation de la présente invention, la direction prédéterminée est une droite horizontale dans l'image I. Selon un mode de réalisation de la présente invention, on calcule un histogramme cumulé TA 15 pour chacun des NA angles de mesure. Selon un mode de réalisation pour lequel les projections ont été acquises via une géométrie d'acquisition parallèle, le calcul d'un histogramme cumulé TA pour un angle A donné peut se dérouler selon les deux sous-étapes suivantes : 20 1.1 Calcul d'une rotation l' de l'image I selon l'angle A : Il s'agit, dans un premier temps, de calculer une rotation l' de l'image I pour l'angle de mesure A considéré. L'image l' obtenue après rotation de l'image I a les mêmes dimensions que l'image I. Selon un mode de réalisation de la présente invention, la rotation l' d'une image I selon un 25 angle A peut être réalisée de la façon décrite ci-après. Chaque point /(i', j') de l'image après rotation l' est obtenu à partir d'un point /(i, j) de l'image I par la relation suivante : /'(i',j') = /(i,j) avec i = (i' - ic). cos(A) - (j' - jc). sin(A) + ic , j = (i' - ic). sin(A) + (j' - jc). cos(A) + jc 30 et (ic,jc) sont les coordonnées du centre de l'image I. Selon un mode de réalisation de la présente invention, lorsque les coordonnées (i,j) ne sont pas entières, on a recours à une interpolation pour estimer la valeur /(i, j) de l'image I au point (i,j). Selon un mode de mise en oeuvre de l'invention, l'interpolation utilisée peut être 3033068 11 bilinéaire, ou faire intervenir des fonctions B-spline, préférentiellement cubiques. Si l'on note Interpolation() la fonction d'interpolation choisie, ceil( ) une fonction qui associe, à un nombre réel, l'entier immédiatement supérieur, et floor( ) une fonction qui associe, à un nombre réel, l'entier immédiatement inférieur, l'expression de la valeur /(i, j) de l'image I au point(i, j) peut être 5 donnée de la façon suivante : I (i, j) = Interpolation(I (ceil (j), ceil(i)), I (ceil (j), f loor(i)), I (floor(j), ceil (i), I (f loor(j), floor(j) 1.2 Calcul d'un histogramme cumulé T des niveaux d'intensité sur l'image l' Puis, dans un deuxième temps, on procède au calcul de l'histogramme cumulé TA des niveaux d'intensité rencontrés, c, dans l'image l' obtenue après rotation de l'image I. 10 Selon un mode préférentiel de mise en oeuvre de l'invention, la direction privilégiée est une droite parallèle à l'une des lignes j de l'image l', et on calcule alors le cumul des intensités le long d'une ligne j de cette image. On définit ainsi un histogramme cumulé de dimension NI x H, et un élément TA(i, j) de l'histogramme cumulé contient le nombre de mailles de la ligne j de l'image après rotation l' dont l'intensité est inférieure ou égale à i. Un exemple d'histogramme des 15 intensités cumulées est donné en Figure 2. Les Figures 2A et 2B montrent le mécanisme de remplissage de la première ligne de l'histogramme cumulé, les Figures 2C et 2D le mécanisme de remplissage de la seconde ligne et les Figures 2E et 2F le mécanisme de remplissage de la troisième ligne. Les Figures 2A, 2C, et 2E présentent les valeurs d'intensité de l'image l', composée de 3 lignes, 4 colonnes et de 5 niveaux d'intensité entiers variant de 1 à 5. Les Figures 20 2B, 2D, et 2F présentent les valeurs contenues dans l'histogramme cumulé constitué de 3 lignes et 5 colonnes (correspondant aux 5 niveaux d'intensité) après remplissage respectivement de la première, de la deuxième et de la troisième ligne. On passe par exemple de la Figure 2A à la Figure 2B de la façon suivante : pour remplir l'élément TA(1,1), on compte le nombre de fois où apparaissent des valeurs d'intensité inférieures ou égales à 1 dans la première ligne de l' ; puis on 25 remplit l'élément TA(2,1) en comptant le nombre de fois où apparaissent des valeurs d'intensité inférieures ou égales à 2 dans la première ligne de l'image l', et on répète un tel mécanisme pour l'élément TA(3,1), TA(4,1) et TA (5,1). Ce pré-calcul est très rapide et est réalisé en un seul balayage de toutes les mailles de l'image l' obtenue après rotation de l'image I. De plus, l'histogramme cumulé résultant est organisé 30 en fonction de seuils possibles pour l'image l' puisque un élément TA(i,j) renseigne directement sur les mailles de la ligne j de l'image l' dont les valeurs sont au plus égales à la valeur de seuil i. De façon générale, on définit par TA(sn) l'ensemble des éléments de l'histogramme cumulé représentatifs du nombre de niveaux d'intensité de l'image l' inférieurs à une valeur de seuil sn. 3033068 12 Dans le cas d'une image 2D et dans le cas où ladite direction privilégiée est une droite parallèle à l'une des lignes j de l'image l', TA(sn) correspond à la snième colonne de l'histogramme cumulé TA. Dans le cas d'une image 3D, TA(sn) peut correspondre à un plan dans l'histogramme cumulé TA. 5 Selon un mode de réalisation de la présente invention pour lequel les projections ont été acquises via une géométrie d'acquisition non parallèle, on définit une matrice de projection fonction de l'angle de mesure A considéré et mettant en relation chaque point de l'image I avec un élément de l'histogramme cumulé TA. Puis on parcourt chaque point de l'image I pour remplir 10 l'histogramme cumulé TA de la même façon que dans le cas de la géométrie d'acquisition parallèle décrit précédemment. 2) Estimation des valeurs des seuils d'intensité 15 Au cours de cette étape, il s'agit de déterminer de façon automatique un nombre NS de valeurs de seuils d'intensité d'une image I, le nombre NS de seuils étant un entier supérieur ou égal à 1, prédéfini par le spécialiste. Selon l'invention, ce problème est résolu en cherchant les NS valeurs de seuils permettant de minimiser une fonction objectif F mesurant au moins la différence entre la projection P réelle selon l'angle A et la projection PA' simulée pour ledit angle A et calculée 20 à partir de la distribution représentative des niveaux d'intensités dans l'image I. Selon un mode de réalisation de la présente invention, on calcule une projection P,", d'angle A simulée à partir de l'histogramme cumulé TA d'angle A calculé dans un mode de réalisation de l'étape précédente selon une formule du type : PA = f (TA(si), TA(s.), TA(sNs)) où f est une combinaison linéaire, TA(sn) est l'ensemble des éléments de l'histogramme cumulé TA 25 d'angle A représentatifs du nombre de niveaux d'intensité de l'image l' au plus égales à un seuil sri, avec n variant de 1 à NS. Selon un mode de réalisation de la présente invention, on calcule la projection P,',1 d'angle A simulée à partir de l'histogramme cumulé TA d'angle A calculé dans un mode de réalisation de l'étape précédente de la façon suivante: 30 - dans le cas d'un seul seuil si (cas où NS=1) : /3:4(j) = ci . TA(j, sl) + cZ . (W - TAG Si)) où ci et cz sont les intensités représentatives des deux classes délimitées par le seuil s1; 3033068 13 dans le cas de deux seuils si et s2, avec si REVENDICATIONS1. Procédé pour la segmentation d'une image d'un objet représentée par un nombre fini de niveaux d'intensité, ladite image ayant été reconstruite à partir de mesures externes à l'objet et acquises selon différents angles pris par rapport à une référence, caractérisé en ce qu'on réalise, à partir d'un nombre prédéfini de seuils d'intensité de ladite image, les étapes suivantes pour au moins un desdits angles : a) on détermine, pour ledit angle, une distribution représentative de la répartition desdits niveaux d'intensités de ladite image selon une direction prédéterminée dans ladite image b) on détermine les valeurs desdits seuils d'intensité en minimisant une fonction objectif estimant la différence entre ladite mesure pour ledit angle et une mesure simulée pour ledit même angle, ladite mesure simulée étant fonction desdites valeurs desdits seuils et de ladite distribution ; c) on segmente ladite image par regroupement desdits niveaux d'intensité en classes d'intensité délimitées par lesdites valeurs desdits seuils. 2. Procédé selon la revendication précédente, dans lequel ladite direction prédéterminée est une droite parallèle à l'un des axes horizontaux de ladite image. 3. Procédé selon l'une des revendications précédentes, dans lequel ladite distribution est un histogramme cumulé desdits niveaux d'intensités rencontrés dans ladite image. 4. Procédé selon la revendication 3, dans lequel lesdites mesures ont été acquises dans le cadre d'une géométrie d'acquisition parallèle et dans lequel ledit histogramme cumulé est déterminé sur une rotation de ladite image selon ledit angle. 5. Procédé selon l'une des revendications 3 à 4, dans lequel on calcule ladite mesure simulée P pour un desdits angles A selon une formule du type : PÂ = f(TA(S1), TA(Sn), TA(sNS)) où f est une combinaison linéaire, NS est, ledit nombre prédéfini desdits seuils d'intensité, et TA(sn) est l'ensemble des éléments dudit histogramme cumulé représentatifs du nombre desdits niveaux d'intensité de ladite image inférieurs à une valeur de seuil sn, ledit histogramme étant calculé pour ledit angle A. 3033068 19 6. Procédé selon l'une des revendications précédentes, dans lequel on détermine une intensité représentative pour chacune desdites classes délimitées par lesdites valeurs desdits seuils. 7. Procédé la revendication 6, dans lequel ladite intensité représentative de ladite classe est la valeur d'un desdits seuils délimitant ladite classe, une moyenne desdites valeurs desdits seuils délimitant ladite classe, ou toute combinaison linéaire desdites valeurs desdits seuils délimitant ladite classe. 8. Procédé selon l'une des revendications précédentes, dans lequel on applique les étapes a) et b) pour chacun desdits angles de mesure et on détermine des valeurs de seuils 10 représentatives de l'ensemble desdits angles en choisissant soit la médiane, soit une combinaison desdites valeurs desdits seuils déterrninés pour chacun desdits angles de mesure. 9. Procédé selon l'une des revendications précédentes, dans lequel on détermine des valeurs de seuils représentatives de l'ensemble desdits angles de mesure en minimisant une 15 fonction objectif estimant les différences entre lesdites mesures et lesdites mesures simulées pour l'ensemble desdits angles de mesure. 10. Utilisation du procédé selon l'une des revendications précédentes pour la caractérisation des supports de catalyseurs dans le domaine du raffinage pétrochimique. 11. Utilisation du procédé selon l'une des revendications 1 à 9 pour la caractérisation des 20 milieux poreux dans le domaine de l'exploration et de l'exploitation pétrolière. 12. Utilisation selon l'une des revendications 1 à 9 pour la caractérisation de matériaux, pour le contrôle non destructif d'un objet, ou pour l'imagerie médicale. 13. Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, 25 comprenant des instructions de code de programme pour la mise en oeuvre du procédé selon l'une des revendications 1 à 9, lorsque ledit programme est exécuté sur un ordinateur.