Brevet FR2943145A1 — publié le 20100917

APPAREIL DE TOMOGRAPHIE OPTIQUE COHERENTE PLEIN CHAMP EN LUMIERE INCOHERENTE A SPECTRE LARGE AVEC OPTIQUE ADAPTATIVE. La présente invention concerne un système de tomographie optique cohérente plein champ en lumière incohérente à spectre large. Elle concerne également un procédé de calibration d'un tel appareil et un produit programme d'ordinateur pour mettre en oeuvre le procédé. La technologie de tomographie optique cohérente (OCT) plein champ utilisant une source de lumière incohérente à spectre large est décrite en particulier dans la demande de brevet FR2817030. Grâce à l'emploi d'une source lumineuse à spectre large, supérieur à 300nm, donc à très faible cohérence temporelle, cette technologie permet d'obtenir des résolutions spatiales d'environ 1 micron dans les trois dimensions pour une profondeur atteignant 1 millimètre dans les tissus vivants. Cette technologie permet, en pratique, d'obtenir une résolution anatomopathologique permettant de � voir � les détails d'une coupe de tissu vivant en trois dimensions. Cette technologie utilise des optiques de focalisation à faible ouverture donc moins onéreuses que les optiques utilisées, par exemple, en microscopie confocale. En particulier dans le cadre d'étude de tissus vivants, on utilise des optiques à immersion permettant de s'affranchir des différences d'indice optique entre l'air et l'échantillon étudié et de mieux adapter les différences de marche dans les deux bras de l'interféromètre. En effet, ces optiques sont telles qu'un liquide de transmission est placé entre l'optique et l'échantillon, liquide ayant un indice optique proche de celui de l'échantillon. Classiquement, l'eau, avec un indice de 1,33, est utilisée car proche de l'indice du tissu vivant, lui-même composé en majorité d'eau. La demande de brevet FR2889584 décrit l'utilisation de ce type d'optique avec un interféromètre OCT. Cependant, l'observation montre que la qualité de l'image décroît avec l'augmentation de profondeur des couches visualisées de l'échantillon. Cela est dû à ce que les couches supérieures de l'échantillon, que doit traverser le signal optique de mesure, ont des indices optiques légèrement différents de celui du milieu d'immersion. Ce qui modifie la longueur du chemin de mesure et rend donc difficile la conservation de l'égalité du chemin optique entre le bras de référence et le bras de mesure, en déformant le plan d'onde. Il est en effet important de se rappeler qu'en OCT plein champ, on ne procède pas à un balayage de la surface à mesurer mais on observe toute la surface à un instant donné. En particulier, sur la surface de l'échantillon visualisée, le chemin optique de chaque élément de cette surface varie en fonction des éléments rencontrés sur son parcours. Autrement dit, pour chaque voxel le chemin de mesure varie. La demande de brevet FR 2 791 548 propose d'utiliser les principes de l'optique adaptative pour corriger ces aberrations optiques qui peuvent se produire dans la zone de propagation libre du faisceau lumineux, par exemple les aberrations introduites par l'oeil lors de l'observation rétinienne. Pour cela, une mesure est effectuée avec l'aide d'un faisceau laser envoyé sur une zone proche de la zone à étudier. Le front d'onde de ce faisceau laser reçu à travers l'instrument d'optique est déformé par des aberrations du plan d'onde de mesure et un dispositif d'optique adaptative, par exemple un miroir déformable, est piloté pour minimiser les déformations du front d'onde qui sont mesurées avec un analyseur de front d'onde laser. Cette correction est alors appliquée pour la visualisation de la zone à étudier. Cette adaptation des techniques d'optique adaptative à l'observation microscopique a permis d'améliorer significativement la qualité et la précision des images obtenues. Cependant ces techniques sont relativement complexes à mettre en oeuvre et augmentent significativement les temps d'observation puisque qu'il faut d'abord illuminer la zone d'intérêt avec le faisceau laser, en déduire les corrections à appliquer puis faire l'observation proprement dite. Lorsque, comme dans l'OCT, il est nécessaire de multiplier les images pour obtenir une vue de l'échantillon en 3 dimensions, cette augmentation de temps peut devenir prohibitive. Cela nécessite de plus une parfaite synchronisation entre le dispositif de mesure et le dispositif d'observation. Et, lorsque l'image à corriger se situe en profondeur, le front d'onde laser à analyser n'est pas accessible car la lumière est diffusée par la traversée de l'échantillon. Il sera donc avantageux d'obtenir un système d'OCT permettant d'obtenir des images en 3 dimensions avec une résolution de l'ordre du micron qui s'affranchisse aisément des aberrations créées par des variations de chemin optique, sans sonde laser, sans manipulation particulière ni augmentation excessive des temps d'observation. Pour résoudre un ou plusieurs des inconvénients cités précédemment, un appareil de tomographie optique cohérente plein champ en lumière incohérente à spectre large comporte : • une source lumineuse ayant un spectre supérieur à 300 nm, • deux chemins d'onde lumineuse définissant un chemin de mesure et un chemin de référence, issus de la source lumineuse, • une platine support d'un objet à mesurer sur le chemin de mesure ; • une optique de focalisation à faible ou moyenne ouverture disposée sur le chemin de mesure de façon à focaliser le chemin de mesure sur une zone de mesure de l'objet à mesurer, définissant ainsi un front d'onde de mesure, • un capteur 2D des signaux lumineux issus des interférences des ondes lumineuses issues du chemin de mesure et du chemin de référence, ledit capteur étant connecté à un calculateur adapté pour générer une image à partir des signaux d'interférence captés, et • une optique adaptative située sur l'un des chemins et paramétrable pour modifier la longueur dudit chemin de façon différentiée pour au moins deux parties différentes de la zone de mesure de façon à ce que le chemin de référence et le chemin de mesure soient équivalents pour lesdites au moins deux parties. L'optique adaptative comporte un analyseur de contraste de l'image obtenue connectés à un contrôleur de l'optique adaptative adapté pour maximiser le contraste de ladite image pour lesdites au moins deux parties, le chemin de mesure étant alors équivalent au chemin de référence pour lesdites au moins deux parties. Ainsi la correction est réalisable automatiquement sans nécessité d'appareillage complexe complémentaire tel qu'un laser. Des caractéristiques ou des modes de réalisation particuliers, utilisables seuls ou en combinaison, sont : • le contrôleur est adapté pour paramétrer l'optique adaptative afin d'obtenir un contraste maximal par utilisation de la méthode des polynômes orthogonaux. Ce qui permet avantageusement de régler localement l'optique adaptative pour une aberration donnée sans interférence entre aberrations séparées spatialement. • pour maximiser le contraste, le contrôleur est adapté pour faire varier la longueur du chemin de façon à fournir à l'analyseur au moins 3 mesures correspondant à 3 longueurs différentes pour chaque au moins deux parties de la zone de mesure, et ledit contrôleur étant adapté pour, en outre, déduire le paramétrage de ladite optique adaptative maximisant ledit contraste à partir desdites mesures. Ce qui permet avantageusement de limiter le temps de calibration. • l'optique adaptative est disposée sur le chemin de mesure ; ce qui permet avantageusement de corriger directement les variations du chemin de mesure liées aux aberrations. • l'optique adaptative est disposée au conjugué du plan focal image de l'optique de focalisation ; ce qui permet avantageusement d'obtenir une modification optimale du chemin lié à une aberration. • la conjugaison de l'optique adaptative avec l'optique de focalisation est réalisée par un dispositif afocal ; ce qui permet avantageusement d'éviter toute déformation parasite de l'onde lumineuse par le système optique. Dans un deuxième aspect de l'invention, un procédé de calibration d'un appareil de tomographie optique cohérente plein champ en lumière incohérente à spectre large tel que décrit ci-dessus comporte : • une pluralité de mesures de contraste d'une partie de l'image, chaque mesure correspondant à un paramétrage différent de l'optique adaptive tel que la longueur du chemin pour ladite partie varie ; • détermination du paramétrage de l'optique adaptative pour ladite partie tel que le contraste de ladite partie d'image est maximal, ledit paramétrage correspondant à ce que le chemin de mesure et le chemin de référence sont équivalents pour ladite partie d'image. Dans des modes particuliers de réalisation du procédé : • la détermination du paramétrage est réalisée par la méthode des polynômes orthogonaux ; • il ne comporte que 3 mesures de contraste pour ladite partie d'image. Dans un troisième aspect de l'invention, un procédé de calibration d'un appareil de tomographie optique cohérente plein champ en lumière incohérente à spectre large tel que décrit ci-dessus, mais sans obligatoirement une optique adaptative, comporte : • une pluralité de mesures de contraste de l'image, chaque mesure correspondant à une longueur différente du chemin de référence ; • détermination de la longueur du chemin de référence correspondant à un contraste maximal de l'image, ladite longueur correspondant également à la longueur équivalente dudit chemin de mesure. Cela permet avantageusement de faire une mise au point automatique d'un appareil OCT. Dans un quatrième aspect de l'invention, un produit programme d'ordinateur comportant des instructions de code de programme pour l'exécution des étapes des procédés ci-dessus lorsque ledit programme est exécuté sur un ordinateur. L'invention sera mieux comprise à la lecture de la description qui suit, faite uniquement à titre d'exemple, et en référence aux figures en annexe dans lesquelles : - la figure 1 est une vue schématique d'un interféromètre de Linnik dans un mode de réalisation de l'invention ; - la figure 2 est une vue schématique des aberrations autour d'un voxel ; et - la figure 3 est un ordinogramme d'un procédé de calibration selon un mode de réalisation de l'invention. Dans la description qui va suivre, on appelle surface isoplanétique une surface d'observation de forme sensiblement plane, qui remplit la condition que les fronts d'onde issus d'un point quelconque lui appartenant sont identiquement déformés. Une telle surface peut, par exemple, être constituée matériellement par une partie de plan de coupe dans un milieu diffusant. La zone d'isoplanétisme est donc la zone sur laquelle la qualité de l'image observée est stationnaire quand on corrige les aberrations : par exemple pour un défaut de mise au point toute l'image est corrigée si on remet au point. Par contre, pour une aberration qui dépend de la position dans le champ, toute l'image ne peut pas être corrigée en même temps. En référence à la figure 1, un interféromètre 1 utilisable pour réaliser une tomographie optique cohérente (OCT) plein champ en lumière incohérente à spectre large. est représenté sous forme d'un microscope de Linnik. Il comprend quatre organes principaux que sont un bloc source 3 incohérente de lumière polychromatique, un dispositif d'optique asservi 5, un bloc support d'échantillon 7 et un dispositif de captation et de réalisation d'image 9. Comme indiqué dans le préambule, le bloc source 3 génère une lumière à spectre large, typiquement supérieur à 300nm, donc à très faible cohérence temporelle. Le dispositif d'optique asservi 5 comporte un miroir séparateur 11 de la lumière issue du bloc source 3 permettant de définir deux chemins optiques : un chemin de mesure 13 et un chemin de référence 15. Sur le chemin de référence 15 est disposé un miroir de référence 17 disposé sur une platine oscillante 19 capable de faire osciller à faible fréquence, typiquement 50 Hz, et sur une très faible amplitude, inférieure à 1 micron, le miroir dans l'axe du chemin optique. Ainsi la longueur du chemin optique de référence 15 varie selon une oscillation déterminée par la platine oscillante 19. Un dispositif de focalisation 21 à faible ou moyenne ouverture sous forme d'une optique à immersion est disposée sur le chemin de mesure 13. Un dispositif optique d'ajustement 23 est associé au miroir de référence 17 de telle sorte que le chemin de référence 15 soit égal au chemin de mesure 13 à moins de la longueur de cohérence de la source lumineuse 3 près. Dans un souci de simplification du dessin, la conjugaison du miroir de référence 17 avec la pupille du dispositif de focalisation 21 n'a pas été représentée. Le dispositif de captation et de réalisation d'image 9 comporte une caméra CCD à 2 dimensions 25 ayant une fréquence d'échantillonnage multiple entier de la fréquence d'oscillation du miroir de référence 17 et synchronisée avec celle-ci. Un calculateur 27 récupère les images de la caméra CCD et les intègre pour fournir une image avec une résolution en x, y et z de l'ordre du micron. Ainsi, pour une mesure donnée, l'image représente une zone, ou un plan de mesure, de l'objet à mesurer, correspondant au front d'onde de mesure. Une description plus complète d'un système d'OCT plein champ en 20 lumière incohérente à spectre large est contenue dans le document FR 2 817 030 auquel l'homme du métier est invité à se référer. L'interféromètre 1 comporte en outre, sur le chemin de mesure 13, une optique adaptative 28, par exemple un miroir déformable, capable de modifier la longueur du chemin de référence pour chaque élément du chemin 25 lumineux correspondant à une partie de la zone de mesure. Cette optique adaptative est placée de telle sorte qu'elle est conjuguée de la pupille de l'objectif 21 ou plus exactement du plan image de l'objectif ou plan Fourier. Pour ne pas déformer l'onde plane, par exemple en l'absence d'aberrations, qui irait converger au foyer de l'objectif, on réalise 30 cette conjugaison à l'aide d'un système afocal. Cette optique adaptative est pilotée par des moyens de pilotage 29 connectés au calculateur 27, et plus particulièrement à un analyseur de contraste 30. Le fonctionnement du système d'OCT étant connue, seul le fonctionnement du dispositif d'optique adaptative sera décrit ci-après. La figure 2 représente un faisceau incident 31 focalisé sur un voxel de l'image 33 et ayant un front d'onde plan symbolisé par des flèches d'égale longueur. Par voxel, il faut comprendre ici le volume élémentaire au niveau du front d'onde. Le voxel a une section correspondant à la tâche de diffraction en l'absence d'aberrations et la tranche 39 d'où revient le signal rétrodiffusé à une épaisseur égale à la longueur de cohérence temporelle de la source divisé par deux fois l'indice de réfraction n du milieu. Sur le chemin lumineux se trouve des structures 35 jouant le rôle de diffuseur. Les diffuseurs ont une position aléatoire. Le signal 37 rétrodiffusé résulte donc de la somme des ondelettes rétrodiffusées pour les diffuseurs 35 du volume du voxel considéré. Ainsi, d'un voxel à l'autre, l'onde rétrodiffusée 37 a une amplitude et une phase qui varient aléatoirement. Cette onde, d'amplitude et de phase aléatoires, reste une onde plane si le milieu traversé n'introduit pas de distorsion du front d'onde. L'onde 37 correspond à un front d'onde distordu symbolisé par des flèches de longueurs inégales. Le signal optique ainsi générée est connu des utilisateurs de l'OCT sous le terme anglais de �speckle � ou en français de � tavelures �. Initialement, ce terme � speckle � désigne la granularité laser qui arrive quand un objet diffusant est éclairé par un laser très cohérent. Dans le cas de l'OCT, le phénomène physique est de nature différente même si, mathématiquement, il se représente dans les deux cas par une marche aléatoire dans le plan complexe : partie réelle - partie imaginaire. Traditionnellement, les utilisateurs d'OCT minimisent les fluctuations spatiales de � speckle � dans les images par du filtrage numérique des données reçues en utilisant une diffusion grossissante de région (en anglais � region growing diffusion �) ou une décomposition en ondelette. Un � speckle � possède une propriété particulière dans la distribution des intensités : le contraste, qui est égal au rapport de l'écart type à la valeur moyenne pour la distribution des intensités, est égal à l'unité. De plus, l'addition de N figures de � speckle � non corrélées divise le contraste par la racine carrée de N. Par exemple, la superposition de 16, respectivement 100, � speckles � non corrélés divise le contraste par 4, respectivement 10. Cette propriété est applicable aux aberrations géométriques introduites par les hétérogénéités macroscopiques des structures traversées de tailles largement supérieures à celles microscopiques des diffuseurs. Ces aberrations ont pour effet un élargissement de la tache focale. L'image obtenue est donc le produit de convolution de l'image idéale, sans aberration, par la tache de diffusion qui résulte de la présence d'aberrations. Cette tache de diffusion peut être considérée comme générée par N �speckles � non corrélés. La convolution consiste donc à faire la superposition de N speckles non corrélés, avec la partie de contraste en rapport. Ainsi, on constate, que l'image obtenue contient une information, le contraste, qui est directement liée à la présence des aberrations qui doivent être corrigées par le dispositif d'optique adaptative, ce dispositif devant introduire une correction variable du front d'onde en fonction de l'aberration constatée pour un voxel donné. Chaque aberration contribuant à la tâche de diffusion peut être associée à des déformations du front d'onde du chemin de mesure. Or ces déformations peuvent s'exprimer sur une base de polynômes orthogonaux permettant de modifier le niveau d'une déformation présente sans toucher aux niveaux des autres déformations. La méthode des polynômes orthogonaux parfois également nommés polynômes de Lukosz est décrite par exemple dans le document de D. Debarre, D. M. J. Booth et T. Wilson � Image based adaptative optics through optimization of low spatial frequences � Optics Express (2007), vol. 15 (13) pp 8176-8190. Ainsi, il est possible de piloter l'optique adaptative pour corriger une aberration sans toucher aux niveaux des autres aberrations. Et, pour chaque aberration, le dispositif de pilotage va, en coopération avec le calculateur d'images, rechercher le paramétrage du dispositif d'optique adaptative qui va maximiser le contraste de l'image à la hauteur de cette aberration. La propriété des polynômes orthogonaux choisis est telle que trois mesures différentes du niveau de l'aberration à corriger suffisent pour déterminer les paramètres de la correction à appliquer. Ainsi, le procédé de calibration et de correction des aberrations d'un appareil d'OCT plein champ en lumière incohérente à spectre large comporte, figure 3, les étapes de : • faire, étape 41, une pluralité de mesures de contraste pour la partie d'image à corriger tout en faisant varier la longueur du chemin optique par l'optique adaptative ; • déterminer, étape 43, le paramétrage de l'optique adaptative qui permet la maximisation du contraste. En corrigeant ainsi l'ensemble des aberrations constatées, on obtient ainsi en tout point une correction du chemin qui permet le maintien de l'équivalence des chemins de mesure et de référence et donc un fonctionnement optimal de l'interférométrie OCT. Il est à noter que la corrélation d'un défaut de mise au point lié à la différence d'indice est réalisable par la même méthode en se basant sur la différence de marche, donc sans miroir déformable ou autre dispositif d'optique adaptative. L'analyse de contraste se base alors sur l'image entière ou sur une partie représentative de celle-ci. Avec une optique adaptative, il est alors aisé de faire une correction complète en effectuant dans un premier temps une correction globale de mise au point puis en utilisant l'optique adaptative pour réaliser les corrections locales, pour chaque groupe de voxels. Il est à noter que ce procédé de calibration peut être mise en oeuvre en utilisant un produit programme d'ordinateur comportant des instructions de programme pour exécuter ce procédé sur l'ordinateur. Cependant cet ordinateur commande alors des actionneurs par logiciel pour appliquer les modifications de paramétrage sur l'optique adaptative, et plus généralement pour régler l'interféromètre OCT. L'invention a été illustrée et décrite en détail dans les dessins et la description précédente. Celle-ci doit être considérée comme illustrative et donnée à titre d'exemple et non comme limitant l'invention a cette seule description. De nombreuses variantes de réalisation sont possibles. Par exemple, l'optique adaptative peut être installée sur le chemin de référence. Il est constaté cependant que les corrections sont plus aisément réalisées lorsque celle-ci est sur le chemin de mesure comme décrite ci-dessus. Dans les revendications, le mot � comprenant � n'exclue pas d'autres éléments et l'article indéfini � un/une � n'exclue pas une pluralité. REVENDICATIONS1. Appareil de tomographie optique cohérente plein champ en lumière incohérente à spectre large comportant : • une source lumineuse (3) ayant un spectre supérieur à 300 nm, • deux chemins d'onde lumineuse définissant un chemin de mesure (13) et un chemin de référence (15), issus de la source lumineuse, • une platine (7) support d'un objet à mesurer sur le chemin de mesure ; • une optique de focalisation (21) à faible ou moyenne ouverture disposée sur le chemin de mesure de façon à focaliser le chemin de mesure sur une zone de mesure de l'objet à mesurer, définissant ainsi un front d'onde de mesure, • un capteur 2D (25) des signaux lumineux issus des interférences des ondes lumineuses issues du chemin de mesure et du chemin de référence, ledit capteur étant connecté à un calculateur (27) adapté pour générer une image à partir des signaux d'interférence captés, et • une optique adaptative (28) située sur l'un des chemins et paramétrable pour modifier la longueur dudit chemin de façon différentiée pour au moins deux parties différentes de la zone de mesure de façon à ce que le chemin de référence et le chemin de mesure soient équivalents pour lesdites au moins deux parties, caractérisé en ce que l'optique adaptative (28) comporte un analyseur de contraste (30) de l'image obtenue connectés à un contrôleur (29) de l'optique adaptative adapté pour maximiser le contraste de ladite image pour lesdites au moins deux parties, le chemin de mesure étant alors équivalent au chemin de référence pour lesdites au moins deux parties. 2. Appareil selon la revendication 1, caractérisé en ce que le contrôleur est adapté pour paramétrer l'optique adaptative afin d'obtenir un contraste maximal par utilisation de la méthode des polynômes orthogonaux. 3. Appareil selon la revendication 2, caractérisé en ce que, pour maximiser le contraste, le contrôleur (29) est adapté pour faire varier la longueur duchemin de façon à fournir à l'analyseur au moins 3 mesures correspondant à 3 longueurs différentes pour chaque au moins deux parties de la zone de mesure, et ledit contrôleur étant adapté pour, en outre, déduire le paramétrage de ladite optique adaptative maximisant ledit contraste à partir desdites mesures. 4. Appareil selon l'une quelconque des revendications précédentes, caractérisé en ce que l'optique adaptative est disposée sur le chemin de mesure. 5. Appareil selon la revendication 4, caractérisé en ce que l'optique adaptative est disposée au conjugué du plan focal image de l'optique de focalisation (21). 15 6. Appareil selon la revendication 4, caractérisé en ce que la conjugaison de l'optique adaptative avec l'optique de focalisation est réalisée par un dispositif afocal. 7. Procédé de calibration d'un appareil de tomographie optique cohérente 20 plein champ en lumière incohérente à spectre large comportant : • une source lumineuse (3) ayant un spectre supérieur à 300 nm, • deux chemins d'onde lumineuse définissant un chemin de mesure (13) et un chemin de référence (15), issus de la source lumineuse ; • une platine (7) support d'un objet à mesurer sur le chemin de mesure ; 25 • une optique de focalisation (21) à faible ou moyenne ouverture disposée sur le chemin de mesure de façon à focaliser le chemin de mesure sur une zone de mesure de l'objet à mesurer, définissant ainsi un front d'onde de mesure, • un capteur 2D (25) des signaux lumineux issus des interférences des 30 ondes lumineuses issues du chemin de mesure et du chemin de référence, ledit capteur étant connecté à un calculateur (27) adapté pour générer une image à partir des signaux d'interférence captés, et10• une optique adaptative (28) située sur l'un des chemins et paramétrable pour modifier la longueur dudit chemin de façon différentiée pour au moins deux parties différentes de la zone de mesure de façon à ce que le chemin de référence et le chemin de mesure soient équivalents pour lesdites au moins deux parties ; ledit procédé comportant : - une pluralité (41) de mesures de contraste d'une partie de l'image, chaque mesure correspondant à un paramétrage différent de l'optique adaptive tel que la longueur du chemin pour ladite partie varie ; - détermination (43) du paramétrage de l'optique adaptative pour ladite partie tel que le contraste de ladite partie d'image est maximal, ledit paramétrage correspondant à ce que le chemin de mesure et le chemin de référence sont équivalents pour ladite partie d'image. 8 Procédé de calibration selon la revendication 7, caractérisé en ce que la détermination du paramétrage est réalisée par la méthode des polynômes orthogonaux. 9. Procédé de calibration selon la revendication 8, caractérisé en ce qu'il ne comporte que 3 mesures de contraste pour ladite partie d'image. 1O.Procédé de calibration d'un appareil de tomographie optique cohérente plein champ en lumière incohérente à spectre large comportant : • une source lumineuse (3) ayant un spectre supérieur à 300 nm, • deux chemins d'onde lumineuse définissant un chemin de mesure (13) et un chemin de référence (15), issus de la source lumineuse ; • une platine (7) support d'un objet à mesurer sur le chemin de mesure ; • une optique de focalisation (21) à faible ou moyenne ouverture disposée sur le chemin de mesure de façon à focaliser le chemin de mesure sur une zone de mesure de l'objet à mesurer, définissant ainsi un front d'onde de mesure, • un capteur 2D (25) des signaux lumineux issus des interférences des ondes lumineuses issues du chemin de mesure et du chemin deréférence, ledit capteur étant connecté à un calculateur (27) adapté pour générer une image à partir des signaux d'interférence captés, ledit procédé comportant : - une pluralité de mesures (41) de contraste de l'image, chaque mesure correspondant à une longueur différente du chemin de référence ; - détermination (43) de la longueur du chemin de référence correspondant à un contraste maximal de l'image, ladite longueur correspondant également à la longueur équivalente dudit chemin de mesure 11. Produit programme d'ordinateur comportant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une quelconque des revendications 7 à 10 lorsque ledit programme est exécuté sur un ordinateur.