Brevet FR3089286A1 — publié le 20200605

Description Titre de l'invention : Procédé et système pour mesurer une surface d'un objet comprenant des structures différentes par interfé- rométrie à faible cohérence Domaine technique 100011 La présente invention concerne un procédé pour la mesure de la surface d'un objet comprenant des structures différentes par interférométrie à faible cohérence. 100021 Le domaine de l'invention comprend, de manière non limitative, celui des dispositifs de mesure et de contrôle pour l'industrie de la microélectronique ou de l'optique intégrée. Etat de la technique [0003] La profilométrie optique permet classiquement de déterminer le profil ou la to- pographie d'une face ou d'une surface d'un objet. Dans le cas de la profilométrie optique par interférométrie à faible cohérence, une source de lumière à large spectre est utilisée. En faisant varier les chemins optiques des rayonnements de référence ou d'inspection, un signal d'interférence (interférogramme) est obtenu. C'est un signal en intensité qui est fonction de la différence des chemins parcourus par les rayonnements. [0004] L'onde réfléchie sur, ou émergeant de, la surface de l'objet. à mesurer subit un déphasage qui dépend des caractéristiques physiques des matériaux présents que l'onde traverse et/ou sur lesquels elle se réfléchit. Ces caractéristiques sont les indices de réfraction complexes et/ou les épaisseurs des empilements de couches des matériaux transparents dans la gamme de longueur d'onde de la source de lumière. Ce déphasage à la réflexion se rajoute au déphasage dû à la topographie. [0005] Lorsque les motifs ou structures présents sur la surface de l'objet analysé sont réalisés avec des matériaux différents ou des empilements de matériaux différents, ou des empilements multicouches transparents, le déphasage à la réflexion est alors différent pour chacune de ces structures. [0006] De par ce fait, la hauteur relative obtenue à partir des signaux d'interférences pour des structures différentes peut être erronée. Dc même, le cas échéant, des valeurs d'épaisseurs de couches de motifs ou de structures calculées peuvent être erronées. [0007] Il est connu de calculer des phases théoriques à la réflexion en exploitant les in- formations que Von a sur les structures, et notamment les matériaux ou les empilements de couches de matériaux, présents. On construit ainsi des modèles théoriques de rétlectance que l'on utilise pour corriger les mesures. [0008] Il est également connu de calculer des modèles théoriques de réflectance pour des empilements de matériaux connus au préalable et d'appliquer des méthodes de 2 recalage de ces modèles sur des signaux provenant des interférogrammes pour en déduire des mesures d'épaisseurs de couches minces (inférieures à un micromètre). [0009] Pour corriger des mesures topographiques et/ou évaluer des épaisseurs de couches à partir de modèles théoriques, il est cependant nécessaire de connaître ou d'identifier les différentes structures de matériaux présentes dans le champ de vue afin de pouvoir appliquer le modèle approprié. Ceci peut être problématique pour des surfaces comportant des structures multiples. En effet, de manière générale, la nature dcs structures présentes dans le champ de vue et leur répartition spatiale (layout en langue anglaise) n'est pas connue. [0010] On tonnait par ailleurs le document « Interferometric microscope with Crue color imaging », J.L. Bevcrage et al., Proc. SPIE 9203, August 2014, qui décrit un procédé basé sur la génération d'images de couleurs en plus des acquisitions de signaux interférométriques dans des microscopes interférométriques pour visualiser la présence de matériaux différents dans le champ de vue à partir de leur réflectance spectrale, analysée par exemple sur trois longueurs d'onde (rouge, vert et bleu). Toutefois, ce procédé nécessite plusieurs acquisitions (une par couleur) et des modifications du système intcrférométrique avec l'ajout ou la modification de sources de lumière, filtres ou systèmes de détection, rendant ainsi cette approche très complexe et coûteuse en termes de temps de mesure et d'équipements. De plus, en général, on est limité à l'utilisation de trois longueurs d'ondes fixes qui ne sont pas nécessairement les plus discriminatives concernant la réflectance des différents empilements de matériaux de l'objet mesuré. [0011] L'invention a pour objectif de résoudre les problèmes décrits ci-dessus. Exposé de l'invention [0012] Un but de la présente invention est dc proposer un procédé et un système de mesure de la surface d'un objet pouvant comprendre au moins deux structures différentes qui soient simple à mettre en oeuvre et qui ne nécessitent notamment pas l'utilisation de plusieurs éclairages de couleur ou de capteurs. [0013] Un autre but dc la présente invention est de proposer un procédé et un système dc mesure dc la surface d'un objet pouvant comprendre au moins deux structures différentes permettant d'obtenir des mesures de profils de la surface d'un objet corrigées des effets dus aux matériaux ou aux structures présents à la surface de cet objet. [0014] Il est également un but de la présente invention de proposer un procédé et un système de mesure de la surface d'un objet pouvant comprendre au moins deux structures différentes permettant d'obtenir des mesures d'épaisseurs de couches transparentes présentes sur la surface de l'objet. [0015] Il est aussi un but de la présente invention de proposer un procédé et un système de 3 mesure de la surface d'un objet pouvant comprendre au moins deux structures différentes permettant de détecter et caractériser ces structures avec des appareils de mesure actuels avec très peu de, voire sans, modification de leur architecture actuelle. [0016] Au moins un de ces objectifs est atteint avec un procédé de mesure d'une surface d'un objet comprenant au moins une structure par interférométrie optique à faible cohérence, le procédé comprenant les étapes suivantes : acquisition d'un signal interférométrique en une pluralité de points, dits de mesure, de ladite surface dans un champ de vue ; [0017] pour au moins un point de mesure : attribution du signal interférométrique acquis à une classe de signaux interférométriques parmi une pluralité de classes, chaque classe étant associée à un signal interférométrique de référence représentatif d'une structure type ; et analyse du signal interférométrique pour en déduire une information sur la structure au point de mesure, en fonction de sa classe. [0018] L'objet peut être ou comprendre notamment un substrat pour l'industrie de la micro- électronique ou de l'optique intégrée sous ses formes, tel qu'une plaquette (ou wafcr) de semi-conducteur ou de tout autre matériau (verre), un substrat sur film (« dicing frarne » en langue anglaise), un substrat avec des circuits intégrés et des interconnexions (« electrical interconnects » en langue anglaise), etc. [0019] Dans le présent document, le ternie « structure » désigne tout agencement pouvant se trouver à la surface de l'objet, tel que des pistes, couches de matériaux ou empilements de couches de matériaux diélectriques ou conducteurs, composants optiques ou électroniques, etc. Une structure peut ainsi être formée d'au moins une couche de matériau. [0020] Un objet peut ainsi comprendre au moins deux, ou une pluralité de structures dif- férentes réparties sur sa surface. [0021] Les signaux interférométriques peuvent être obtenus par un système de profilométrie. Ces signaux interférométriques peuvent contenir les informations relatives à une altitude ou une hauteur en différents points de la surface à mesurer. Ils peuvent également contenir des informations relatives à la nature de la surface elle-même (matériaux, épaisseurs de couches, etc.). [0022] L'étape d'attribution des signaux interférométriques acquis à des classes de signaux interférométriques, proposé par le procédé selon l'invention, permet ainsi de classifier les signaux interférométriques individuellement selon le type de la structure de laquelle chaque signal a été obtenu. Le type de la structure (par exemple, une simple couche transparente ou un empilement de couches minces de matériaux transparents) est représenté par une structure type avec un signal interférométrique de référence associé. [0023] L'étape d'attribution correspond à une segmentation du champ de vue mesuré sur la surface de l'objet. Le procédé selon l'invention permet ainsi d'aboutir à une seg- mentation de l'image de la surface de l'objet à mesurer. [0024] Une fois attribué à une classe, chaque signal intcrférométrique est analysé selon un modèle, ou un modèle profilométriquc, correspondant ou propre à sa classe. Une classe peut être définie, par exemple, par un modèle de réflectance, ce modèle de réflectance prenant en compte le ou les matériaux et la géométrie (épaisseurs des couches) de la structure. [0025] Ainsi, grâce au procédé selon l'invention, il est possible d'obtenir une information sur la répartition de différentes structures sur la surface d'un objet, en exploitant un signal intcrférométrique mesuré par interférométrie optique à faible cohérence, sans avoir besoin d'une connaissance précise de la topologie ou de la localisation de ces structures. [0026] Le procédé selon l'invention permet également d'obtenir des mesures de profils ou de forme de la surface de l'objet corrigées des effets (en particulier de phase) dus aux matériaux et/ou aux structures présentes, et/ou des mesures d'épaisseurs de couches, sans avoir une connaissance préalable précise de la topologie de ces matériaux et/ou de ces structures sur la surface de l'objet. [0027] En outre, l'étape d'attribution peut être réalisée par traitement numérique ou analogique, et n'impacte pas, ou peu, les étapes de mesures des signaux interférométriques. Par conséquent, le procédé selon l'invention peut être mis en oeuvre par les dispositifs de mesure actuels sans, ou avec très peu de, modification architecturale, en particulier de la partie de détection de ces dispositifs. [0028] Selon un mode de réalisation, l'étape d'attribution du signal interférométrique à une classe peut comprendre une étape d'initialisation de classes. [0029] Cette étape d'initialisation de classes peut comprendre la fourniture d'au moins l'un des éléments suivants - une liste de classes ; - une estimation a priori d'un nombre de classes présentes dans le champ de vue. [0030] Cette étape d'initialisation peut résulter en un ensemble de classes susceptibles d'être représentées par les signaux interférométriques mesurés dans le champ de vue. Elle permet ainsi de disposer d'un ensemble de classes correspondant, ou susceptibles de correspondre, aux structures réellement présentes dans le champ de vue. Cette étape d'initialisation peut également résulter dans un nombre de classes défini a priori, permettant une répartition des signaux interférométriques dans ce nombre de classes. [0031] Avantageusement, la fourniture d'une liste de classes peut comprendre au moins l'une des étapes suivantes : - détermination de classes à partir de mesures de signaux interférométriques sur une surface de référence (du même objet ou d'un autre objet) ayant des structures connues ; - détermination de classes à partir de connaissances à priori (telles que des masques, dessins ou descriptions) sur les structures présentes sur la surface de l'objet. [0032] On dispose ainsi d'une liste de classes susceptibles d'être présentes sur l'objet à mesurer avec leurs signaux interférométriques de référence associés. Une liste de classes ainsi fournie est robuste et permet, par la suite, d'utiliser des méthodes de classification simples. [0033] Cette liste de classes peut être, par exemple, obtenue à partir de masques, dessins (« layout ») ou descriptions des structures présentes sur l'objet, lorsque ces informations sont disponibles. On peut alors en déduire par calcul des signaux d'interférence de référence théoriques et ainsi déterminer un signal interférométrique de référence pour chaque classe. [0034] La liste de classes peut également être obtenue expérimentalement, à partir, par exemple, de mesures effectuées sur un objet de référence avec des structures connues. [0035] Alternativement ou en plus, l'étape d'initialisation de classes peut comprendre, ou être réalisée par, les étapes suivantes : répartition des signaux interférométriques acquis selon des régions homogènes du champ de vue, où les régions sont homogènes selon un critère prédéfini ; détermination des classes à partir des signaux interférométriques dans lesdites régions homogènes. [0036] L'étape d' initialisation selon ce mode de réalisation peut également être qualifiée de « segmentation préliminaire » Elle réalise en effet une première segmentation du champ de vue en régions homogènes (selon le critère prédéfini) qui peut être utilisée pour rendre l'étape d'initialisation des classes plus robuste et plus précise en termes de choix ou de construction de classes, et améliorer la classification pixel par pixel des signaux interférométriques. Les classes ainsi déterminées peuvent être obtenues directement à partir des signaux interférométriques clans les zones homogènes, par exemple à partir de moyennes ou d'autres paramètres statistiques. Les signaux interférométriques clans les zones homogènes peuvent être également utilisés pour identifier des classes parmi une liste prédéfinie de classes. [0037] Les régions homogènes peuvent être notamment déterminées en exploitant une propriété des signaux interférométriques, telle qu'une intensité, une réflectance ou une réflectance spectrale de ces signaux interférométriques. Le critère prédéfini peut être par exemple lié à des variations d'intensité ou de réflectance maximales admissibles. [0038] Avantageusement, l'étape d'attribution du signal interférométrique à une classe peut comprendre, ou être réalisée par, la comparaison du signal interférométrique acquis au signal interférométrique de référence. [0039] De préférence, cette comparaison peut comprendre, ou être réalisée par, la déter- mination d'une distance (telle qu'une distance euclidienne) entre le signal interféro- 6 métrique acquis et le signal interférométrique de référence. [0040] Alternativement ou en plus, l'étape d'attribution du signal interférométrique à une classe peut comprendre, ou être réalisée par, une itération des étapes suivantes : association du signal interférométrique à une classe ; pour chaque classe, détermination d'un nouveau signal interférométrique de référence à partir des signaux interférométriques attribués à cette classe. L'itération est continuée jusqu'à ce qu'un critère de convergence prédéterminé soit satisfait. [0041] Ainsi, en effectuant une itération des étapes d'association du signal interférométrique à une classe et de détermination d'un nouveau signal de référence pour cette classe, la précision de l'identification des classes pour chaque signal interférométrique peut être améliorée. En outre, les classes peuvent être initialisées à partir des signaux interférométriques eux-mêmes, dans le cas notamment où l'étape d' initialisation ne fournit qu'un nombre de classes. [0042] Pour ce faire, chaque classe, ou plus particulièrement le signal interférométrique de référence, est déterminé en exploitant les signaux interférométriques attribués à cette classe. Un nouveau signal interférométrique de référence peut, par exemple, être calculé à partir de la moyenne des signaux interférométriques attribués à cette classe, et éventuellement du signal d'interférence de référence précédent. Les itérations sont poursuivies jusqu'à satisfaire un critère de convergence prédéterminé. Le critère de convergence peut comprendre, par exemple, une mesure de variation ou de distance entre des signaux interférométriques issus d'itérations successives. Il peut également être basé sur le fait qu'entre 2 itérations, plus aucun point de mesure ne change de classe. [0043] Avantageusement, l'étape d'analyse du signal interférométrique peut comprendre une identification d'une structure ou d'une propriété d'une structure en fonction de la classe du signal interférométrique. Cette identification peut notamment être réalisée à partir de la détermination de la densité spectrale optique du signal interférométrique, ou de sa réflectance spectrale. [0044] En effet, l'utilisation de la densité spectrale optique ou de la réflectance spectrale permet de différencier des structures ou matériaux présents sur la surface de l'objet, et ce même en utilisant un capteur d'intensité ou d'imagerie produisant des interférogranurtes en niveaux de gris. [0045] Grâce à la connaissance de la classe à laquelle chaque signal interférométrique a été attribué, il est ainsi possible d'attribuer ou de faire correspondre les signaux interférométriques à un type de structure et/ou de matériau que représente leur classe. [0046] Il est ainsi possible de déterminer précisément la géométrie de l'objet, c'est-à-dire, la répartition des différentes structures à la surface de l'objet. 7 [0047] Avantageusement, l'étape d'analyse du signal interférométrique peut comprendre la détermination de l'épaisseur d'une ou de plusieurs couches transparentes en appliquant un modèle de reflectance spectrale dépendant de la classe du signal interférométrique. [0048] On peut ainsi mettre en oeuvre des techniques dc réflectométric connues par recalage de modèle pour la mesure d'épaisseurs dc couche mince. L'attribution du signal interférométrique à une classe permet de choisir le bon modèle. [0049] Avantageusement, l'étape d'analyse du signal interférométrique peut comprendre une détermination de hauteur de la surface au point de mesure en prenant en compte un facteur de déphasage dudit signal interférométrique dépendant de sa classe. [0050] Cc facteur de déphasage permet de corriger le déphasage à la réflexion dû, par exemple, à des pertes d'absorption ou des réflexions dans des couches de matériau. [0051] Avantageusement, l'étape d'acquisition peut être réalisée par une mesure d'un signal interférométriquc pour chaque pixel d'un capteur effectuant une mesure plein champ, où chaque pixel du champ de vue du capteur correspond à un point de mesure. [0052] Avantageusement, le procédé selon l'invention peut être mis en oeuvre avec un objet de type substrat semi-conducteur. [0053] Selon un autre aspect de la même invention, il est proposé un système dc mesure par interférométrie optique à faible cohérence, pour mesurer une surface d'un objet comprenant au moins une structure, le système comprenant : un dispositif d'acquisition d'un signal interférométriquc en une pluralité de points, dits de mesure, de ladite surface dans un champ de vue, et un module de traitement des signaux interférométriques acquis, configuré pour : - attribuer chaque signal interférométrique acquis à une classe de signaux interférométriques parmi une pluralité de classes, chaque classe étant associée à un signal interférométrique de référence représentant une structure (particulière) ; et - réaliser une analyse du signal interférométrique pour en déduire une information sur la structure au point de mesure, en fonction de sa classe. [0054] Le système de profilométrie produit des mesures interférométriques qui exploitent les différences de chemin optique (ou temps de propagation) entre un faisceau de référence et un faisceau de mesure réfléchi par la surface à mesurer. En faisant varier la différence de chemin optique au cours du temps, on obtient un signal interférométrique dans le domaine temporel. [0055] Alternativement, d'autres systèmes de mesures interférométriques peuvent être utilisés, par exemple, permettant d'obtenir des signaux interférométriques dans le domaine spectral, directement fonction des longueurs d'onde utilisées. De tels 8 systèmes peuvent mettre en oeuvre, par exemple, un interféromètre illuminé par une source de lumière à large spectre, avec un dispositif de détection permettant une analyse spectrale du signal interférométrique, tel qu'un spectromètre placé devant des photodétecteurs. 11 est également possible d'utiliser des sources de lumière ajustables en longueur d'onde. [0056] Dans unc version avantageuse, le dispositif d'acquisition d'un signal interféro- métrique peut comprendre un capteur interférométrique plein champ. [0057] Dans ce cas, un point de mesure peut correspondre à un pixel du champ de vue dudit capteur. Un tel capteur matriciel peut être, par exemple, unc caméra CCD (pour Charged Coupled Device ») ou CMOS (pour « Cornplernentary Metal Oxide Semiconductor »). [0058] Le procédé et le système selon l'invention peuvent chacun être utilisé pour la mesure d'une surface d'un substrat semi-conducteur ou d'un substrat transparent dans le spectre visible sous toutes formes possibles (wafer, wafer on frame, panel, etc.). La surface à mesurer peut comporter différents types de structures telles que des couches de matériaux diélectriques ou conducteurs, et des composants optiques ou électroniques. En particulier, le procédé et le système selon l'invention peuvent chacun être utilisé pour la mesure de leurs profils corrigés et de la topologie des structures et/ou matériaux présents sur la surface, ou encore pour imager ladite surface. Un profil correspond à la hauteur relative, ou à la hauteur par rapport à une référence, d'un ensemble de points de la surface. La mesure du profil permet également de faire état de la rugosité de la surface de l'objet. [0059] Plus généralement, le procédé et le système selon l'invention peuvent chacun être utilisé pour la mesure d'une surface d'un objet, en particulier pour la caractérisation et/ ou l'imagerie d'au moins un motif de ladite surface. Description des dessins [0060] D'autres avantages et caractéristiques apparaitront à l'examen de la description détaillée d'exemples nullement limitatifs, et des dessins annexés sur lesquels : [0061] [fig.1] est une représentation schématique d'un exemple de réalisation non limitatif d'un procédé selon l'invention ; [0062] [fig.2] est une représentation schématique d'un exemple de mise en oeuvre du procédé de la Figure 1 ; [0063] [fig.3] est une représentation schématique d'un exemple de réalisation non limitatif d'un système selon l'invention ; [0064] [fig.4] est une représentation schématique d'un exemple de mise en oeuvre du procédé de la Figure 1 ; [0065] [fig.5] est une représentation schématique d'un exemple de mise en oeuvre du 9 procédé de la Figure 1 ; [0066] [fig.6] est une représentation schématique d'un exemple de mise en oeuvre du procédé de la Figure 1 ; [0067] [fig.7] est une représentation schématique d'un exemple de mise en oeuvre du procédé de la Figure 1 ; et [0068] [fig.8] est une représentation schématique d'étapes du procédé selon l'invention selon un exemple non limitatif. Description détaillée de l'invention [0069] 11 est bien entendu que les modes de réalisation qui seront décrits dans la suite ne sont nullement limitatifs. On pourra notamment imaginer des variantes (le l'invention ne comprenant qu'une sélection de caractéristiques décrites par la suite isolées des autres caractéristiques décrites, si cette sélection de caractéristiques est suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique antérieur. Cette sélection comprend au moins une caractéristique de préférence fonctionnelle sans détail structurel, ou avec seulement une partie des détails structurels si cette partie uniquement est suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique antérieur. [0070] En particulier toutes les variantes et tous les modes de réalisation décrits sont com- binables entre eux si rien ne s'oppose à cette combinaison sur le plan technique. [0071] Sur les figures, les éléments communs à plusieurs figures conservent la même référence. [0072] La Figure 1 est une représentation schématique d'un exemple de réalisation non limitatif d'un procédé selon l'invention. [0073] Lc procédé 100 comprend une étape 102 d'acquisition d'un signal interférométrique en une pluralité de points de mesure sur une surface d'un objet, par exemple à l'aide d'un capteur interférométrique plein champ comme détaillé ci-après. Dans ce cas, chaque pixel du capteur correspond à un point de mesure, et un signal interférométrique est mesuré pour chaque pixel dans le champ de vue du capteur. [0074] Un exemple de l'acquisition des signaux interférométriques sera décrit plus en détails en référence à la Figure 3 ci-dessous. [0075] Le procédé 100 comprend en outre une étape 104 d'attribution de chaque signal in- terférométrique mesuré à une classe parmi un ensemble de classes. Chaque classe est représentée par un signal interférométrique de référence. [0076] Dans certains modes de mise en oeuvre, ces classes de signaux interférométriques de référence correspondent, ou sont associées à, des structures ou des types de structures particulières pouvant être présentes sur la surface de l'objet à mesurer. Les signaux interférométriques de référence peuvent alors être choisis de telle sorte à être repré- sentatifs de signaux interférométriques particuliers acquis sur différents types de structures telles que : des structures en matériaux particuliers (silice, silicium, etc.), transparents ou opaques (métal, etc.) ; et/ou des structures constituées de matériaux transparents, ou d'empilements de matériaux transparents, selon des épaisseurs ou des gammes d'épaisseurs particulières. [0077] L'ensemble de classes est construit ou initialisé pendant une étape 106 d'initialisation de classes. [0078] De manière générale, l'initialisation permet de définir et/ou de sélectionner des classes appropriées, c'est-à-dire, des classes qui sont représentatives des types de structures réellement présents dans le champ de vue du capteur. [0079] On peut disposer d'une liste de classes susceptibles d'être présentes sur la surface de l'objet à mesurer avec leurs signaux interférométriques de référence, cette liste pouvant être obtenue à partir de masques, dessins (« layout ») ou descriptions des structures présentes sur l'objet. [0080] Des classes peuvent également être obtenues expérimentalement, à partir, par exemple, de mesures effectuées sur un objet de référence avec des structures connues. [0081] Selon une autre variante, on peut simplement définir un nombre de classes a priori, et utiliser un algorithme de classification, tel que l'algorithme connu des K-moyennes ( « K-Means » en langue anglaise), pour générer les classes à partir de l'ensemble des signaux interférométriques mesurés. Toutefois, cette méthode peut avoir l'inconvénient de générer des classes qui ne sont pas optimales. [0082] Bien entendu, toutes ces variantes peuvent être combinées, avec des classes définies théoriquement, et/ou définies expérimentalement, et/ou définies en nombre. [0083] Selon un mode de réalisation, l'étape d'initialisation 106 du procédé 100 selon l'invention peut comprendre une étape de segmentation préliminaire. Cette étape de segmentation préliminaire peut être réalisée par le regroupement des signaux interférométriques mesurés, dans le champ de vue, selon des régions considérées comme homogènes selon un critère. Ce critère peut être, par exemple, une valeur d'une réflectance globale ou à des longueurs d'ondes particulières. On considère alors qu'une région homogène correspond à une classe, ou en d'autres termes ne comprend que des mesures interférométriques appartenant à une même classe. La segmentation préliminaire permet ainsi plus efficacement de sélectionner ou d'identifier les classes effectivement présentes et appropriées pour effectuer la classification pixel par pixel de l'étape suivante. En effet, la classification pixel par pixel peut permettre une meilleure précision de classification sur des motifs fins (structures à l'échelle de quelques pixels). 11 [0084] Selon un premier exemple, il est possible de ne conserver que des classes qui ont pu être associées aux régions homogènes telles que définies par la segmentation préliminaire, choisies par exemple parmi un ensemble de classes prédéfinics. Il est également possible d'affiner les signaux interférométriques de référence des classes sélectionnées grâce aux données issues des mesures, en remplaçant, par exemple, le signal interférométrique de référence d'une classe par le signal moyen calculé sur l'ensemble des régions homogènes associées à cette classe. [0085] Selon un deuxième exemple, il est possible de définir les classes directement à partir de signaux d'interférences obtenus, par exemple en moyenne, sur des zones homogènes. [0086] Les méthodes de segmentation préliminaire applicables peuvent comprendre, par exemple : des méthodes basées sur une détection des contours, de type ligne de partage des eaux par exemple, ces méthodes visant à utiliser les contours ou les gradients dans l'image pour délimiter les différentes régions ; des méthodes de segmentation basées sur une détection des régions, en utilisant des champs de Markov par exemple, ces méthodes visant à séparer le champ de vue en différentes régions homogènes. [0087] Un exemple de cette étape de segmentation préliminaire sera détaillé ci-dessous. [0088] Lors d'une étape 108 du procédé 100, les signaux interférométriques mesurés I(x,y) sont associés ou attribués, par pixels, aux classes construites préalablement. [0089] De manière générale, cette étape 108 d'association met en oeuvre une comparaison entre chaque signal interférométrique, à chaque point mesuré, et le signal interférométrique de référence de chaque classe. Les signaux interférométriques sont respectivement attribués à la classe pour laquelle la meilleure similitude ou la meilleure proximité avec le signal interférométrique de référence a été trouvée. [0090] La comparaison des signaux interférométriques par pixels avec les signaux interféro- métriques de référence peut être réalisée en effectuant des calculs de distance tels que des calculs de distance Euclidienne, ou toute autre distance ou toute autre mesure mettant en évidence un écart entre deux signaux. [0091] De préférence, la comparaison peut être effectuée notamment en utilisant la ré- flectance spectrale des signaux interférométriques, obtenue par transformée de Fourier des signaux interférométriques temporels mesurés. [0092] Bien entendu, d'autres méthodes de comparaison peuvent être mises en oeuvre, qui exploitent des informations temporelles, de forme ou de phase des signaux interférométriques. [0093] La Figure 2 illustre un exemple d'attribution de signaux interférométriques à des classes prédéfinies, en utilisant les réflectances spectrales. En traits pleins sont re- 12 présentés des réflectances spectrales correspondant aux signaux interférométriques de référence représentatifs de trois classes. En pointillés sont représentés des réflectances spectrales correspondant aux signaux interférométriques mesurés en trois points de mesure ou pixels sur l'image, et attribués, respectivement à ces classes. Chaque paire 35, 36, 37 de signal mesuré - signal de référence représente ainsi un type de structure différente. [0094] Différents algorithmes de classification peuvent être mis en oeuvre pour attribuer les signaux interférométriques aux classes. [0095] Selon des modes de réalisation, les signaux sont attribués aux classes en utilisant des techniques par apprentissage. [0096] Selon un premier exemple, une telle technique d'apprentissage met en oeuvre des machines à vecteurs de support (« Support vector machine » ou SVM en anglais) qui nécessitent un ensemble de classes d'apprentissage. Cet ensemble peut être fourni, par exemple, lors de l'étape étape 106 d'initialisation du procédé 100, et notamment par une connaissance théorique ou expérimentale des caractéristiques de certains signaux interférométriques attendus ou mesurés sur la surface de l'objet. [0097] Selon un deuxième exemple, une technique d'apprentissage implique des réseaux de neurones. Les données d'apprentissage (correspondant aux signaux interférométriques de référence) peuvent ainsi être prédéfinis théoriquement ou expérimentalement ou être déduites lors des apprentissages en acceptant ou en rejetant certaines d'entre elles. [0098] Selon d'autres modes de réalisation, des méthodes de classification simple au plus proche voisin peuvent être utilisées. Pour tout pixel du champ de vue du capteur, son signal interférométrique ou le spectre de réflectance correspondant est comparé à tous les signaux interférométriques de référence ou aux spectres de réflectance de référence correspondant. Cette comparaison peut être effectuée, par exemple, en utilisant une distance euclidienne. Le signal interférométrique du pixel est ensuite attribué à la classe la plus « proche » ou, en d'autres termes, pour laquelle la distance ainsi déterminée est minimale. Cette méthode est notamment appropriée pour classifier des pixels dans un ensemble de classes prédéfinies. [0099] Selon encore d'autres modes de réalisation, on peut également utiliser des méthodes fondées sur une opération de classification, par exemple celle des K-moyennes évoquée précédemment. Les points de mesure associés aux pixels de l'image sont ainsi divisés en k groupes ou classes issues, par exemple, de l'étape 106 d'initialisation du procédé 100 selon l'invention. [0100] Pour l'algorithme des K-moyennes, les classes peuvent être issues notamment d'une information a priori qui peut être le nombre de classes (valeur de k) attendues. Alternativement, l'algorithme peut également être initialisé avec un ensemble de classes issues de connaissances complètes ou approximatives des caractéristiques des signaux 13 d'interférence ou des signaux de réflectance spectrale de référence associés à des structures mesurées ou théoriques présentes sur la surface de l'objet. Cela permet d'initialiser l'algorithme avec une meilleure précision et de minimiser les erreurs de classification. [0101] L'étape d'attribution ou association 108 peut ensuite être réalisée par la minimisation d'une fonction de distance qui peut être la somme des carrés des distances entre les signaux de réflectance spectrale de chaque pixel par rapport aux réflectances spectrales de référence de la classe à laquelle le pixel est attribué. Ces signaux de référence correspondent aux signaux moyens en chaque classe. [0102] Un exemple d'application de la technique des K-moyennes sera décrit ci-dessous. [0103] De manière avantageuse, le procédé 100 selon l'invention peut également comprendre une étape 110 de détermination, pour chaque classe à laquelle des signaux interférométriques ont été associés, d'un nouveau signal interférométrique de référence à partir des signaux interférométriques attribués à cette classe. Un nouvel ensemble de classes est ainsi formé. Les nouveaux signaux interférométriques de référence (ou les réflectances spectrales correspondantes) peuvent, par exemple, être calculés à partir de la moyenne des signaux interférométriques attribués à la classe concernée et, éventuellement, du signal d'interférence de référence initial. Ensuite, les signaux interférométriques mesurés sont de nouveau attribués aux nouvelles classes (étape 112 de la Figure 1). Les itérations sont poursuivies jusqu'à satisfaire, par exemple, un critère de convergence prédéterminé. Ce critère de convergence peut, par exemple, prendre en compte le déplacement des centres de classes : la convergence est atteinte à partir du moment où les centres ne bougent plus (ou moins qu'une certaine valeur). Ce critère de convergence peut également prendre en compte le fait que la convergence est atteinte lorsque plus aucun pixel (ou un nombre de pixels inférieur à une valeur) ne change de classe d'une itération à l'autre. [0104] L'itération des étapes d'attribution 108 et de détermination 110 d'un nouveau signal de référence est particulièrement adaptée pour la classification selon la méthode des K- moyennes, les classes pouvant être améliorées de manière itérative à partir des mesures. [0105] L'étape d'attribution 104 des signaux interférométriques à des classes se termine après la classification de tous les signaux. [0106] En référence à la Figure 1, le procédé 100 selon l'invention comprend en outre une étape 114 de traitement ou d'analyse des signaux interférométriques associés aux différents points de mesure. De manière avantageuse, cette étape 114 de traitement est réalisée selon des méthodes propres aux classes auxquelles ils sont associés. [0107] Le procédé de l'invention permet ainsi notamment d'appliquer des analyses parti- culières, ou des traitements particuliers, aux signaux interférométriques des différents 14 pixels, en fonction de la classe à laquelle ils appartiennent ou de la structure correspondante. Il est ainsi possible, par exemple, d'effectuer : un calcul de statistique différenciée, par exemple de rugosité, en fonction des différentes structures identifiées ; un calcul d'une cartographie d'estimation d'épaisseur de matériaux transparents, par exemple en appliquant un algorithme de réflectométrie ; des mesures de profils de topographie corrigées, et en particulier affranchies des contributions des caractéristiques physiques des matériaux, telles que leurs indices de réfraction complexes et les épaisseurs pour les empilements de couches dans le cas de matériaux transparents. [0108] Tl est notamment possible d'appliquer, pour obtenir des informations de hauteur corrigées : un facteur de déphasage simple, par exemple comme pour le cas d'une couche métallique ; un modèle de réflectance d'empilement de couches, pour obtenir une estimation de l'épaisseur d'une couche transparente. [0109] Par la suite, un exemple de mise en oeuvre du procédé selon l'invention sera décrit de manière détaillée. [0110] En référence à la Figure 3, un exemple de système profilométriquc ou profilomètre optique mis en oeuvre selon l'invention est décrit. [0111] Ce système profilométrique ou profilomètre optique 1 comprend un dispositif interfé- rométrique 13 qui peut être de type Michelson, Mirau ou autre. Dans le mode de réalisation présenté, ce dispositif interférométrique 13 est illustré sous la forme d'un objectif interférométrique 13 de type Mirau avec une lentille ou une optique d'imagerie 14, une lame séparatrice 15 semi-réfléchissante et un miroir de référence 16. Le profilomètre optique 100 comprend également un système d'éclairage 10 qui émet une lumière avec un spectre large et qui peut inclure, par exemple, une diode électroluminescente LED (pour « Light Einitting Diode ») ou une lampe halogène, et un capteur d'imagerie ou une caméra 11 avec par exemple un capteur d'imagerie matriciel de type CCD ou CMOS. Le profilomètre optique 1 comprend également un dispositif de balayage 18 actionné, par exemple, par un actuateur piézoélectrique qui permet de faire varier relativement les trajets optiques dans les deux bras de l'interféromètre, entre la lame séparatrice 15 et, respectivement, la surface de l'objet à mesurer 17 et le miroir de référence 16. Le profilomètre optique comprend enfin un module de traitement d'image 19 ou de signaux interféroméniques, agencé pour mettre en oeuvre le procédé de traitement des mesures selon l'invention. [0112] La Figure 4 illustre une image 30 en intensité d'un champ de vue d'un substrat semi- conducteur sur lequel des structures constituées d'empilements de matériaux différents ont été réalisées. La Figure 4 montre également des exemples de signaux intcrférométriques 31, 32, 33, 34 correspondants à des types de structures différents. [0113] On peut observer que les différentes zones distinctes sur l'image d'intensité 30 génèrent des signaux d'interférences 31-34 significativement différents. Ces signaux acquis aux différents points de l'image 30 dans le champ de vue peuvent être identifiés et exploités selon le procédé 100 de l'invention, afin de pouvoir les traiter de manière différentiée en fonction de la nature des structures présentes dans ce champ de vue. [0114] Dans le mode de mise en oeuvre décrit, les signaux d'interférence sont analysés et classifiés en exploitant leur réflectance spectrale. [0115] Un signal d'interférence en lumière blanche (ou polychromatiquc) peut être décomposé en une collection de signaux d'interférence en éclairage monochromatique. Cette décomposition du signal d'interférence peut se faire en fonction du nombre d'onde k, avec k = 2rsfiL,7,. étant la longueur d'onde. [0116] Dans un cas pour lequel la contribution liée à l'angle d'incidence est négligée (approximation paraxialc), le signal d'interférence à un pixel donné peut être modélisé par l'équation suivante : [0117] [Math.1] I( ) (k, e )1/ (k )01k [0118] I est l'intensité lumineuse mesurée sur le pixel donné au pas de déplacement du dispositif de balayage 18, V(k) est la distribution du spectre de la source lumineuse, et g(k,,§) peut s'écrire : [0119] [Math.2] g(k,e) = R(k) Z(k) 4- 2 -IR (k).Z(k).cos [k(h-e) 4 y(k)[ [0120] Z correspond à la rétlectance effective de l'objet à mesurer 17. Ce terme inclut la ré- flectance spectrale de l'objet au pixel donné et la transmittance du système de mesure comprenant la lame séparatrice 15 (ou le cube séparateur). [0121] 11 est possible d'écrire : [0122] [Math.3] k ) = (k ),p(k) 101231 T représente la transmittance du système de mesure, comprenant la lame séparatrice 15, et p la réllectance spectrale de l'objet à mesurer 17. [0124] La réflectance effective de référence, notée R, inclut les réflectances ou trans- mittances des éléments optiques du système 1, telles que la réflectance du miroir de référence 16 et la transmittance de la lame séparatrice 15. Lc paramètre h représente la hauteur de l'objet et le paramètre y est un terme de déphasage de l'onde lumineuse à la réflexion sur l'objet et à travers le système de mesure. [0125] La Figure 4 montre des exemples de signaux d'interférence mesurés ou 16 d' interférogrammes 1(e) indiqués par les signes 31-34. [0126] Pour plus de clarté dans les notations, on ne mentionne pas jusqu'ici et dans tout le reste du document la dépendance de la localisation (x,y) dans l'image (ou dans le plan conjugué par le système optique d'imagerie de l'objet) des ternies V, R, Z, h et y. Celle-ci n'est cependant pas négligée. [0127] L'amplitude des sinusoïdes formant des signaux d'interférence en lumière mono- chromatique est proportionnelle à : [0128] [Math.4] V(k)-\/R(k)Z(k) [0129] Elle est l'expression du module de la transformée de Fourier de l' interférogramme correspondant à la densité spectrale optique P(k). [0130] Une étape de calcul de la transformée de Fourier de l' interférogramme permet d'obtenir une mesure de la densité spectrale optique de la lumière réfléchie par l'objet, à travers le module de la transformée de Fourier. La densité spectrale permet d'obtenir une information sur la réflectance spectrale p(k) de l'objet. [0131] Lc champ de vue du système 1 ou microscope interférométrique tel qu'illustré à la Figure 3 est décomposé en sous-régions, en général définies par les pixels de la caméra de détection 11. Un interférogramme, ou signal interférométrique, est ainsi acquis pour chaque pixel. [0132] Selon une variante, il est aussi possible de considérer la phase de la transformée de Fourier de l' interférogramme, qui permet d'obtenir la réflectance complexe de l'empilement de matériaux. [0133] Dans la suite du document, nous ne décrivons que l'utilisation du module de la transformée de Fourier de l'interférogramme correspondant à la densité spectrale optique P(k), ce qui n'est nullement limitatif. [0134] La variable d'intérêt étant la réflectance spectrale p(k), il peut être utile de normaliser la densité spectrale optique P(k) afin d'améliorer les classifications et/ou calculs des épaisseurs et/ou topographies subséquents. [0135] Une étape de normalisation du module de la transformée de Fourier de l'interférograrnme (et donc de la densité spectrale optique) permet de déduire une réflectance spectrale en tout point (x,y) du champ de vue. [0136] Pour ce faire, plusieurs stratégies sont envisageables, notamment une normalisation à partir d'un signal de référence, ou une auto-normalisation, qui seront décrits ci-après. [0137] Le principe de la normalisation à partir d'un signal de référence consiste à réaliser une acquisition d'un signal interférométrique sur un objet ou échantillon de référence connu recouvert ou constitué d'un matériau homogène. Cet échantillon de référence peut être, par exemple, une plaque de silicium. A partir du signal acquis, le module de la transformée de Fourier sur cet échantillon de référence est calculé en chaque point 17 de l'image, ce qui fournit une densité spectrale optique de référence. [0138] Il est ensuite possible de déterminer, en chaque pixel de l'image, les valeurs de ré- flectance spectrale 0(k) en calculant le ratio entre la densité spectrale optique sur l'objet à mesurer et la densité spectrale optique de référence. Bien entendu, les valeurs p(k) ne sont calculées que dans la bande spectrale du système d'éclairage 10, pour des valeurs de densité spectrale optique de référence non nulles. [0139] Cette approche nécessite de calibrer le système 1 de mesure sur un objet de référence connu (voire de recalihrcr périodiquement le système). La qualité de la normalisation dépend grandement de la régularité de déplacement du dispositif de balayage 18 le long de sa course et de la stabilité de l'éclairage d'une acquisition à l'autre. [0140] Suivant une variante, la densité spectrale optique de référence peut être une valeur globale ou moyenne calculée sur toute la surface de l'échantillon de référence si les caractéristiques optiques du système de mesure sont suffisamment homogènes. [0141] Selon le principe de l'auto-normalisation, on pose comme approximation que la plupart des grandeurs ne sont pas dépendantes de leur position sur l'image (éclairage homogène, pas de variation spatiale des propriétés du miroir de référence 16 et de la lame séparatrice 15). Seules les valeurs de la réflectance spectrale p(k) sont considérées comme pouvant varier selon la localisation sur l'image en fonction des motifs de l'objet. On calcule le module moyen de la transformée de Fourier sur l'image pour obtenir ainsi une expression de la densité spectrale optique de référence. On réalise alors le rapport en chaque pixel du module de la transformée de Fourier sur le spectre moyen (ou de la densité spectrale optique sur la densité spectrale optique de référence). [0142] L'approche d'auto-normalisation a pour avantage de ne pas nécessiter d'étape de ca- libration au préalable. Cependant, on obtient alors une carte de pseudo-réflectance comme le ratio entre la réflectance du pixel donné à chaque nombre d'onde avec la moyenne des rétlectances à ce nombre d'onde. [0143] Un autre exemple de mesure est illustré sur la Figure 5. Les épaisseurs des couches des structures 41, 42 d'un objet 40 représenté sur la 6a doivent être déterminées. Il est connu que les structures à mesurer sont formées par des couches d'oxyde de différentes épaisseurs et déposées sur du silicium. Une valeur précise des épaisseurs des couches d'oxyde pour chaque point d'une image de mesure de l'objet 40 doit être déterminée. [0144] Des signaux d'interférence mesurés associés à un pixel localisé dans chacune des deux régions, montrés sur la Figure 5b, et définies par les structures sont représentés sur la Figure 5c. [0145] Sur la base de mesures de signaux d'interférence sur un objet de référence (par exemple identique à celui illustré sur la Figure 7a), ou de connaissances théoriques à 18 priori, on définit les empilements attendus sur la surface de l'objet. Dans le cas de l'exemple de la Figure 5, les épaisseurs théoriques ou attendues sont, respectivement, 1,2 pm de silice (SiO2) sur un substrat de silicium (Si) d'une part et 0,1 pm de SiO2 sur un substrat de Si d'autre part. [0146] Il est ainsi possible de calculer des réflectances spectrales théoriques, ou de référence, pour générer des classes de signaux interférométriques comme décrit précédemment. [0147] La Figure 6 illustre un exemple de réflectance spectrale théorique 51 pour une couche de 1,2 pm de SiO2 déposée sur un substrat de Si. [0148] Ces signaux de réflectance spectrale théorique constituent les signaux interféro- métriques de référence associés aux deux classes de référence utilisées pour la segmentation de l'objet, ces classes correspondant aux deux types de couches ou structures connues. [0149] A titre d'exemple, la réflectance spectrale théorique ou de référence correspondant au cas d'une couche transparente simple déposée sur un milieu réfléchissant, telle qu'illustrée sur la Figure 6, peut être construite mathématiquement, en incidence normale, grâce à l'équation suivante : [0150] [Math.5] [0151] Les coefficients roi et sont les coefficients de réflexion de Fresnel qui dépendent de la longueur d'onde. Le coefficient roi correspond à l'interface entre le milieu ambiant et le matériau transparent, et le coefficient r12 à l'interface entre la couche transparente et le substrat. [0152] La constante [3 est définie comme p= 2A, où d est l'épaisseur de la couche, et n est son indice de réfraction qui dépend également de la longueur d'onde. [0153] Sur la Figure 6 est également montrée une mesure de réflectance spectrale 52, attribuée à la structure de référence correspondant à une couche de SiO, du fait de sa plus proche ressemblance à cette référence en termes de réflectance spectrale. 101541 La Figure 7a illustre des résultats de la classification, ou segmentation, obtenue par l'étape d'attribution des signaux interférométriques à des classes du procédé selon l'invention. Cette classification aboutit à une répartition des pixels de l'image entre les deux structures de référence connues, à partir d'une attribution des signaux interférométriques mesurés à l'une des deux classes préalablement définies et correspondant à ces structures de référence. 101551 Une fois le pixel attribué à la classe de structure de référence correspondante, il est possible d'estimer la mesure (l'épaisseur de silice déposée sur la couche de silicium. 19 Un rccalage des paramètres de la courbe théorique de la réflectance spectrale (voir l'équation Math. 1) avec la mesure de réflectance du pixel permet ainsi d'obtenir une valeur estimée de l'épaisseur d pour chaque pixel mesuré. [0156] Les Figures 7h et 7c montrent la mise en correspondance des réflectances théoriques 73, 74 sur les réflectances mesurées 71, 72 après ajustement du paramètre d'épaisseur de couche pour un pixel de chacune des deux classes, respectivement, de la Figure 7a. Ce paramètre d'épaisseur donne ainsi une estimation de l'épaisseur de la couche transparente de silice en chaque pixel. Pour la Figure 7h, l'épaisseur est estimée à 1,185 pm, et pour la Figure 7c, l'épaisseur est estimée à 0,140 prn, ces deux valeurs étant cohérentes avec la forme de l'objet attendue. [0157] Ce principe de mesure peut être étendu à des empilements multiples de couches ou des indices de réfraction complexes. [0158] Par la suite, deux exemples de mise en oeuvre de la segmentation de la surface de l'objet seront décrits plus en détails : premièrement, une implémentation de la classification des signaux interférométriques en classes basée sur l'algorithme des K-moyennes, et deuxièmement, une implémentation de la classification des signaux interférométriques en classes basée sur une combinaison d'un algorithme des K-moyennes et d'un algorithme de ligne de partage des eaux sur des données multispectrales pour réaliser une segmentation préliminaire et initialiser l'algorithme des K-moyennes. L'algorithme des K-moyennes est connu par ailleurs, il n'est donc pas détaillé ici. [0159] Utilisation de l'algorithme des K-moyennes pour la classification des signaux interfé- rométriques [0160] Dans cet exemple, on considère qu'il est possible de définir des réflectances spectrales de référence (ou des signaux interférométriques de référence) à partir d'une information approximative et a priori sur les structures ou les empilements sur la surface d'un objet à mesurer. Ces réflectances spectrales de référence forment un ensemble de classes prédéfinies. [0161] L'architecture 60 de l'algorithme est illustrée à la Figure 8. [0162] Les réflectances spectrales de référence peuvent servir à initialiser un algorithme des K-moyennes en définissant des positions initiales des centres de classes (étape 62). Ensuite, à chaque pixel du champ de vue est assigné la classe dont le centre est plus proche en termes de signature spectrale, en utilisant, par exemple, une distance Euclidienne (étape 64). Puis pour chaque classe, le centre prédéfini est remplacé par le centre de gravité défini par le spectre moyen de tous les pixels assignés à cette classe (étape 66). Les étapes 64 et 66 sont itérées (flèche 68) jusqu'à ce que les centres des classes ne changent plus, au sens d'un critère de convergence basé, par exemple, sur des variations de réflectances spectrales de référence d'une itération à l'autre, ou de changements de classe des pixels entre deux itérations. [0163] L'utilisation de l'algorithme des K-moyennes parait particulièrement adaptée par rapport à une simple classification au plus proche voisin directement sur les réflectances spectrales de référence car elle permet une adaptation de ces réflectances spectrales de référence (déplacement des centres). [0164] Selon une variante, il peut être utile d'ajuster les réflectances spectrales de référence car on ne cannait pas exactement la composition des empilements, en particulier des épaisseurs exactes de couches transparentes. De plus, sur toute la surface de l'objet, les structures peuvent présenter des variations d'épaisseur. Les réflectances spectrales de référence ne se superposent donc en général pas exactement avec les réflectances spectrales des pixels attribués à la classe correspondante. [0165] S'il existe des zones dont on ne peut définir une réflectance spectrale de référence, en cas par exemple de présence de motifs sous-résolus ou de zones dont l'opérateur ignore complètement la composition, il est possible, pour cette étape d'attribution des pixels aux classes, de définir une initialisation partiellement ou complètement empirique. L'utilisateur peut, par exemple, désigner des pixels types sur un champ de vue de référence. Les signatures spectrales de ces pixels types servent alors à initialiser les centres de classes et l'algorithme s'exécute comme décrit précédemment. [0166] Après convergence des K-moyennes sur des données où toutes les classes sont re- présentées, il est possible d'utiliser les réflectances spectrales des centres, ou spectre moyen de chaque classe, pour initialiser la classification pour tout champ de vue sur l'objet. [0167] Lorsque toutes les classes ne sont pas nécessairement représentées sur chaque champ de vue, ne pas en tenir compte conduirait à un effet de sur-segmentation du champ de vue. Des régions homogènes peuvent, par exemple, être divisées en plusieurs zones et donc mal classifiées, les empilements de matériaux seraient donc mal identifiés pour une grande partie des pixels. Il est nécessaire de supprimer les classes inutiles pour chaque prise de vue. [0168] On peut établir différents critères pour sélectionner les classes à supprimer : [0169] - Avant convergence des K-moyennes, on peut comparer le spectre moyen de la classe à la première itération à son spectre de référence. Au-delà d'une valeur seuil, on estime que la position du centre a divergée de sa position initiale, et la classe est considérée comme invalide et est supprimée. [0170] - Avant convergence des K-moyennes, on peut considérer le nombre de pixels associés à chaque classe. En dessous d'un nombre seuil, la classe est considérée comme non représentative et est supprimée. [0171] - Après convergence des K-moyennes, on peut comparer la distance entre les centres de classes 2 à 2. En dessous d'un certain seuil, on considère que des pixels correspondant à un même empilement de matériaux. Les classes concernées sont alors fu- 21 sionnécs. [0172] Utilisation conjointe d'un algorithme de segmentation préliminaire: algorithme de ligne de partage des eaux multispectral [0173] Comme décrit précédemment, pour la classification pixel igue, et en particulier pour le mode d'implémentation choisi des K-moyennes, l'initialisation des classes est une étape importante. En effet, si un trop grand nombre de classes sont définies ou si ces classes ne correspondent pas aux empilements réellement présents dans le champ de vue, le résultat de classification est compromis. [0174] Il est ainsi possible d'effectuer une étape de segmentation préliminaire sans tenir compte de modèles de réflectance préétabli et de classifier ensuite les pixels région par région. De manière avantageuse, les résultats de cette étape de segmentation préliminaire peuvent être utilisés pour sélectionner les classes et ajuster les spectres de référence avant d'effectuer une classification pixel par pixel sur ce sous-ensemble de classes. [0175] Dans l'exemple illustré, on effectue cette segmentation préliminaire en utilisant la ré- flectance, ou la réflectance spectrale issue des signaux interférométriques. En effet, dans beaucoup de cas, comme par exemple illustré à la Figure 4, cette réflectance spectrale varie entre les différentes structures ou matériaux. [0176] Pour effectuer la segmentation préliminaire, un algorithme de ligne de partage des eaux sur données multispectrales tel que décrit par Li et al., « Muhispectral image segmentation by a multichannel watershed-based approach », International Journal of Remote Sensing, 2004, vol. 28, p. 4429-4452, est utilisé. Il s'agit d'une approche basée sur les contours, c'est-à-dire qu'on considère que les régions homogènes au sens de la réflectance et que l'on cherche à segmenter, sont séparées par des gradients multispectraux. [0177] La première étape de l'algorithme de segmentation préliminaire est le calcul de la carte de gradient multispectral. [0178] Le gradient multispectral est défini comme suit. Les données d'entrée de l'algorithme sont un cube de réflectance, chaque pixel est défini par un vecteur de dimension N, où N est le nombre de longueurs d'onde dont on dispose pour décrire la réflectance de ce pixel. L'objectif est de créer une image de gradient en niveau de gris à partir des données multispectrales pour pouvoir appliquer ensuite un algorithme classique de ligne de partage des eaux. [0179] La seconde étape de l'algorithme de segmentation préliminaire est l'application de l'algorithme de ligne de partage des eaux sur la carte de gradient calculée (carte en niveaux de gris). [0180] On obtient ainsi une répartition des signaux interférométriques en régions homogènes. 'n [0181] On peut alors déterminer des classes à partir des régions homogènes. On peut, par exemple, calculer pour chaque région un signal interférométrique de référence à partir des signaux interférométriques de cette région, et l'utiliser pour initialiser une classe. Bien entendu, on peut générer des classes à partir de signaux interférométriques issus de régions distinctes mais considérés comme similaires, etiou utiliser des connaissances a priori sur l'objet pour générer les classes. [0182] Une fois ces étapes de segmentation préliminaire effectuées, la classification des signaux interférométriques en classes pixel par pixel peut être réalisée, notamment en utilisant l'algorithme des K-moyennes comme expliqué précédemment. [0183] Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de invention. 23 [Revendication 1]