FR 3155342 A1 20250516 FR 2312406 A 20231113 Un procédé d’imagerie par visualisation de variations d’un contraste de tavelure permet de révéler de manière plus fiable des mouvements de parties d’une cible diffusante à analyser (100). Le procédé comprend d’utiliser une cible de référence (101) pour vérifier qu’une source laser d’éclairage qui est utilisée ne subit pas de fluctuations de fonctionnement pendant une exécution du procédé. L’utilisation de la cible de référence permet en outre de connaître un contraste de tavelure maximum que peut générer la source laser d’éclairage. Le procédé peut aussi être poursuivi en construisant une image de variation de contraste de tavelure, en soustrayant chaque valeur de contraste calculée pour la cible à analyser à une valeur moyenne de contraste de tavelure calculée pour la cible de référence. Ainsi, des images de variation de contraste de tavelure qui sont construites en utilisant des sources laser différentes, notamment avec des largeurs spectrales différentes, des optiques différentes de formation d’image et des capteurs d’images différents, et/ou lors de répétitions du procédé qui sont séparées dans le temps, peuvent être comparées entre elles valablement. Pour cela, chaque image de contraste de tavelure de la cible à analyser est corrigée d’une contribution spécifique à la source laser d’éclairage et éventuellement en plus d’une autre contribution due à un bruit de détection d’image. Figure d’abrégé : Figure 1 PROCEDE D’IMAGERIE PAR VISUALISATION DE VARIATIONS D’UN CONTRASTE DE TAVELURE La présente description concerne un procédé d’imagerie par visualisation de variations d’un contraste de tavelure, ou «speckle contrast» en anglais. L’imagerie par variations de contraste de tavelure est utile pour révéler des mouvements, ou des micromouvements, qui se produisent dans un milieu («medium» en anglais). Pour cela, le milieu doit être semi-transparent ou rugueux, et en plus diffusant pour la longueur d’onde d’éclairage qui est utilisée. Cet éclairage est effectué en utilisant une source laser, pour satisfaire des exigences de cohérence lumineuse qui sont nécessaires pour cette technique d’imagerie. Du fait de sa nature semi-transparente ou rugueuse et en plus diffusante à la longueur d’onde considérée, le milieu éclairé par la source laser produit un état d’interférence lumineuse en chaque point d’un plan d’objet d’un instrument imageur qui est focalisé sur ce milieu. Le plan d’objet contient alors un motif de tavelure, qui est constitué de grains de tavelure juxtaposés, plus clairs ou plus sombres selon l’état d’interférence qui existe en chaque point du plan d’objet. Ce motif de tavelure peut alors être saisi en image en utilisant l’instrument imageur qui conjugue optiquement le milieu à analyser avec la surface photosensible d’un capteur d’images matriciel. La taille et le contraste des grains de tavelure dépendent de plusieurs paramètres qui sont relatifs au milieu analysé, mais aussi de certains paramètres de l’instrument imageur et de la source laser d’éclairage qui sont utilisés. Un procédé d’imagerie par visualisation de variations de contraste de tavelure tel que connu avant la présente invention consiste à construire une image de contraste de tavelure en attribuant à chaque point de l’image telle que saisie par l’instrument imageur une valeur de contraste de tavelure qui existe en ce point dans l’image ou les images saisie(s). Cette valeur de contraste de tavelure peut être calculée selon l’un des trois modes suivants de construction de l’image de contraste de tavelure : premier mode : séparément pour chaque image du milieu qui est saisie en utilisant l’instrument imageur, la valeur de contraste de tavelure en chaque point de la matrice d’image, ou en une sélection de certains d’entre eux, est calculée localement à partir du motif de tavelure tel que saisi par le capteur d’images matriciel, en limitant ce calcul de contraste à un voisinage du point. Puis les valeurs de contraste qui sont obtenues de cette façon pour chacune d’une pluralité d’images saisies successivement, et qui sont relatives à un même point de la matrice d’image, sont moyennées entre elles afin d’obtenir une image au contraste renforcé. Ce premier mode concerne donc un contraste spatial, et fournit une image de contraste spatial ; deuxième mode : une série d’images successives du milieu est saisie en utilisant l’instrument imageur, par exemple plusieurs dizaines ou plusieurs centaines d’images, et la valeur de contraste de tavelure est calculée en chaque point de la matrice d’image à partir des intensités qui ont été saisies pour ce point dans toutes les images successives. Il s’agit alors d’un contraste temporel ; et troisième mode : de nouveau lorsque plusieurs images du milieu sont saisies successivement, mais le contraste de tavelure qui est calculé pour chaque point de la matrice d’image est hybride, en utilisant une formule qui combine les valeurs d’intensité saisies à l’intérieur d’un voisinage de ce point dans toutes les images saisies. Un tel contraste est dit spatio-temporel. Lorsque des micromouvements du milieu sont suffisamment rapides pour se produire en partie pendant la durée d’intégration qui est mise en œuvre par le capteur d’images matriciel afin de saisir chaque image, l’état d’interférence local varie pendant cette durée d’intégration. Il en résulte une réduction de la valeur de contraste spatial de tavelure qui est calculée dans le premier mode de construction de l’image de contraste de tavelure, aux endroits où se produisent les micromouvements. Lorsque le milieu est immobile, les images qui sont saisies successivement sont toutes quasi-identiques, aux bruits de détection près, et le contraste temporel de tavelure qui est calculé dans le deuxième mode de construction de l’image de contraste de tavelure est quasi-nul pour tous les points de la matrice d’image. Mais lorsque le milieu subit des micromouvements locaux entre deux images qui sont saisies successivement, le motif de tavelure varie entre ces images, si bien que la valeur de contraste temporel de tavelure devient non-nulle pour les endroits où se produisent les micromouvements. Si les micromouvements s’accélèrent alors que la durée d’intégration qui est utilisée pour saisir chaque image est maintenue constante, les tavelures vont avoir le temps de beaucoup varier pendant cette durée d’intégration et chaque photodétecteur du capteur d’images va intégrer ces variations. Les valeurs de lecture du capteur d’images pour les images saisies successivement, en un même point d’image, tendent alors vers une même valeur intermédiaire, aboutissant à une réduction de leur valeur d’écart-type, et donc de leur valeur de contraste temporel. En l’absence de diffuseurs fixes, ces valeurs de contraste temporel tendent vers zéro si les micromouvements sont extrêmement rapides par rapport à la durée d’intégration. Le régime opératoire utilisé, correspondant soit à une réduction soit une augmentation du contraste temporel en fonction de la rapidité des micromouvements, dépend des applications. Ainsi, l’imagerie par visualisation de variations d’un contraste de tavelure est efficace pour repérer des micromouvements qui se produisent à la surface d’un milieu rugueux diffusant ou dans un milieu semi-transparent diffusant. Pour cette raison, elle est utile pour de nombreuses applications. Toutefois, il est nécessaire de pouvoir attribuer les variations de contraste de tavelure qui sont mises en évidence à des mouvements du milieu analysé. Pour cela, l’instrument imageur est gardé sans modification entre des images qui sont saisies successivement. Mais certaines caractéristiques de la source laser qui est utilisée pour éclairer le milieu peuvent varier involontairement, en provoquant des contributions aux variations des valeurs calculées de contraste de tavelure, indépendamment de l’existence de mouvements dans le milieu analysé. Il est alors important de pouvoir attribuer au milieu analysé les variations de contraste de tavelure qui sont observées, en excluant que des variations non-contrôlées se soient produites dans le fonctionnement de la source laser. En outre, il peut être nécessaire de connaître précisément certaines des caractéristiques de la source laser qui est utilisée pour éclairer le milieu analysé par contraste de tavelure. Certaines causes de variations non-contrôlées qui affectent le fonctionnement de la source laser utilisée pour éclairer le milieu analysé, sont : /1/ une partie de la lumière d’éclairage qui est émise par la source laser peut se trouver réinjectée dans la source laser elle-même après avoir été rétrodiffusée par le milieu analysé ; /2/ une variation de bruit d’alimentation électrique peut affecter la source laser ; /3/ la source laser peut subir des variations de température, qui altèrent les caractéristiques spectrales du rayonnement d’éclairage ; /4/ la largeur spectrale d’émission de la source laser peut fluctuer, et son fonctionnement monomode peut être instable. Notamment, la source laser peut basculer intempestivement vers un fonctionnement multimode non-désiré ; et /5/ d’autres causes existent encore, selon la nature et la technologie de la source laser utilisée. Il existe des remèdes qui permettent d’éviter ces variations dans le fonctionnement de la source laser, notamment : pour la cause de variations /1/ : utiliser des isolateurs optiques entre la source laser et le milieu analysé ; pour la cause de variations /2/ : utiliser une alimentation électrique pour la source laser qui soit très stabilisée ; pour la cause de variations /3/ : mettre en œuvre une stabilisation thermique pour certains composants de la source laser ; et pour la cause de variations /4/ : il peut être nécessaire de mettre en œuvre un asservissement du fonctionnement de la source laser par rapport à une référence atomique. Toutefois ces remèdes sont complexes et onéreux, voire très onéreux. En outre, ils ne sont pas compatibles avec des applications où l’instrument imageur à utiliser doit être léger, de faible encombrement, mobile et éventuellement à appliquer manuellement sur le milieu à analyser. En outre, des variations non-contrôlées du fonctionnement de la source laser peuvent encore survenir, qui sont autres que celles pour lesquelles des remèdes ont été mis en œuvre. Problème technique A partir de cette situation, un but de la présente invention est de permettre d’attribuer avec certitude des variations de contraste de tavelure à des mouvements qui se produisent dans le milieu analysé. Des buts annexes de l’invention sont d’exclure des résultats de mesure qui ont été obtenus lors de variations du fonctionnement de la source laser d’éclairage, sans engendrer de surcoût significatif pour le système d’imagerie utilisé, ni augmenter son encombrement, et ni réduire sa capacité à être déplacé ou manipulé facilement. Un autre but de l’invention est d’être compatible avec des systèmes d’imagerie par variations de contraste de tavelure qui existent déjà, en effectuant une adaptation en rétrofit de ces systèmes. Un autre but de l’invention est encore de fournir des valeurs et des images de contraste de tavelure qui soient calibrées et dépourvues de contributions spécifiques aux sources laser utilisées pour éclairer le milieu à analyser, ou dans lesquelles de telles contributions dues à la source laser sont réduites. Autrement dit, il est recherché de disposer d’images de contraste de tavelure qui soient dépourvues au moins en partie d’artéfacts causés par la source laser, afin de pouvoir comparer significativement des images de contraste de tavelure saisies en utilisant des sources laser différentes. Un autre but de l’invention est encore de fournir des valeurs et des images de variation de contraste de tavelure en minimisant l’effet sur ces valeurs, des bruits de détection qui se produisent dans le capteur d’images matriciel utilisé. Enfin, un autre but de l’invention est de permettre de caractériser la source laser d’éclairage en temps réel pendant son fonctionnement. Pour atteindre l’un au moins de ces buts ou un autre, un premier aspect de l’invention propose un procédé d’imagerie par visualisation de variations d’un contraste de tavelure, suivant lequel des mouvements de parties diffusantes d’une cible à analyser provoquent des variations d’un contraste de tavelure dans une ou plusieurs images de la cible à analyser, le procédé comprenant les étapes suivantes : /1/ fournir une cible de référence qui possède un motif diffusant statique et constant ; /2/ en utilisant un instrument imageur et une source laser, saisir au moins deux images de la cible à analyser à deux instants différents, en disposant la cible de référence, en plus de la cible à analyser, à l’intérieur d’un champ de vue de l’instrument imageur qui est effectif pour chaque image, de sorte que la source laser éclaire simultanément la cible à analyser et la cible de référence pendant que chaque image est saisie, et que l’image saisie soit formée par du rayonnement produit par la source laser puis diffusé par la cible à analyser et par la cible de référence ; /3/ séparément pour chaque image saisie à l’étape /2/, calculer une valeur de contraste de tavelure, notée C static , à partir d’une partie de l’image qui correspond à la cible de référence ; puis /4/ valider les images saisies à l’étape /2/ uniquement si les valeurs C static calculées respectivement pour ces images ont des variations de l’une quelconque des images à une autre qui sont inférieures à un seuil, sinon recommencer les étapes /2/ à /4/ ; et ensuite /5/ si les images ont été validées, déduire d’au moins une d’entre elles des informations sur des mouvements de parties de la cible à analyser, à partir de variations du contraste de tavelure qui sont relatives à ces parties de la cible à analyser. Au sens de l’invention, on entend par cible à analyser une portion d’un milieu à laquelle le procédé de l’invention est appliqué pour visualiser des variations de contraste de tavelure, quelle que soit la nature du milieu. Pour cela, la cible à analyser est au moins partiellement diffusante pour le rayonnement laser utilisé. En particulier, la cible à analyser peut être constituée par une portion de milieu dans laquelle des mouvements ou micromouvements sont susceptibles de se produire localement ou globalement. On entend par motif statique un motif dont aucune partie n’est en mouvement à un instant considéré ou pendant une durée considérée, par exemple pendant la durée d’intégration qui est mise en œuvre pour saisir une image, couramment aussi appelée durée d’accumulation de l’image. On entend par motif constant un motif qui reste identique entre deux instants séparés, par exemple entre les instants auxquels sont saisies des images successives. On entend enfin par instrument imageur un système qui est adapté pour former et saisir des images, tel que couramment appelé caméra. Un tel instrument imageur comprend une optique de formation d’image, un capteur d’images matriciel, une électronique de lecture et de traitement des signaux de saisie d’image qui sont délivrés par le capteur d’images matriciel, et optionnellement un diaphragme pupillaire et un diaphragme de champ. De préférence, la cible de référence et la source laser sont sélectionnées pour produire une tavelure circulaire Gaussienne, ou tavelure complètement développée («fully developped speckle» en anglais) dans la partie de chaque image saisie à l’étape /2/ qui correspond à la cible de référence. Pour cela, la cible de référence peut posséder une valeur d’écart-type de hauteur de rugosité qui est supérieure à un tiers d’une valeur de longueur d’onde de la source laser, de préférence supérieure à une fois cette valeur de longueur d’onde. En outre, et aussi pour obtenir la tavelure circulaire gaussienne, la cible de référence et l’instrument imageur peuvent être sélectionnés de sorte que la partie de chaque image saisie à l’étape /2/ qui correspond à la cible de référence intègre du rayonnement diffusé par au moins dix, de préférence au moins cinquante, reliefs de rugosité séparés par des creux de rugosité intermédiaires. Selon l’invention, la cible à analyser et la cible de référence apparaissent simultanément dans chaque image qui est saisie. La partie de cette image qui correspond à la cible de référence, par conjugaison optique à travers l’optique de formation d’image, est constituée par un motif de tavelure que produit le rayonnement de la source laser en interagissant avec cette cible de référence. Ce motif de tavelure dans l’image de la cible de référence résulte d’une combinaison de caractéristiques radiatives de la source laser avec des caractéristiques du motif diffusant de la cible de référence. Comme ce dernier est statique et constant, une constance ou quasi-constance du motif de tavelure qui apparaît dans les parties d’images successives qui correspondent à la cible de référence garantit que les caractéristiques du rayonnement de la source laser ont été identiques ou quasi-identiques entre les instants respectifs de saisie des images. Un autre motif de tavelure qui apparaît dans chaque image à des endroits de celle-ci correspondant à la cible à analyser n’est alors pas affecté par des fluctuations de fonctionnement de la source laser. Cet autre motif de tavelure peut donc être utilisé valablement pour en déduire des informations sur la cible à analyser. Plus précisément, lorsqu’une seule image saisie est utilisée pour analyser la cible d’intérêt, le procédé de l’invention permet d’affirmer avec un niveau élevé de certitude que le fonctionnement de la source laser n’a pas varié pendant la durée d’intégration de cette image, l’autre image qui a été saisie dans le procédé de l’invention servant à vérifier la constance du fonctionnement de la source laser. Lorsque plusieurs images saisies sont utilisées pour analyser la cible d’intérêt, le procédé de l’invention permet d’assurer que la source laser a eu un fonctionnement identique pour toutes ces images. Le procédé de l’invention peut être combiné avec chacun des modes de construction de l’image de contraste de tavelure qui ont été rappelés plus haut. Lorsque le premier mode de construction d’image de contraste de tavelure est utilisé, une valeur de contraste spatial de tavelure est calculée localement dans l’une au moins des images qui ont été validées à l’étape /4/, pour certains au moins des points d’image de celle-ci qui correspondent à la cible à analyser. Ensuite, une image de variation de contraste de tavelure peut être construite en attribuant à chacun de ces points d’image une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies conformément à l’étape /2/, et d’autre part la valeur de contraste spatial de tavelure qui a été calculée pour ce point d’image. Lorsque le deuxième mode de construction d’image de contraste de tavelure est utilisé, une série d’images successives de la cible à analyser est saisie conformément à l’étape /2/. Puis, si la série d’images est validée à l’étape /4/, une valeur de contraste temporel de tavelure est calculée séparément pour chacun d’un ensemble de points d’image qui correspondent à la cible à analyser, à partir de valeurs d’intensité respectives des images de la série à ce point d’image. Ensuite, une image de variation de contraste de tavelure peut être construite en attribuant à chacun de ces points d’image une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies conformément à l’étape /2/, et d’autre part la valeur de contraste temporel de tavelure qui a été calculée pour ce point d’image. Lorsque le troisième mode de construction d’image de contraste de tavelure est utilisé, une série d’images successives de la cible à analyser est encore saisie conformément à l’étape /2/. Si cette série d’images est validée à l’étape /4/, une valeur de contraste spatio-temporel de tavelure est alors calculée séparément pour chacun d’un ensemble de points d’image qui correspondent à la cible à analyser, à partir de valeurs d’intensité respectives des images successives dans un voisinage du point d’image. Ensuite, une image de variation de contraste de tavelure peut être construite en attribuant à chacun des points d’image une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies conformément à l’étape /2/, et d’autre part la valeur de contraste spatio-temporel de tavelure qui a été calculée pour ce point d’image. Grâce à l’opération de soustraction introduite selon l’invention dans chacun des trois modes de construction de l’image de contraste de tavelure, l’image de variation de contraste de tavelure qui est obtenue pour la cible à analyser peut être indépendante des caractéristiques radiatives de la source laser, ou dépend de celles-ci dans une mesure résiduelle qui est très petite. Autrement dit, l’invention fournit des valeurs et des images de variation de contraste de tavelure en minimisant l’effet sur ces valeurs et images de la sonde optique qui est utilisée et des conditions de saisie d’images. Ceci est obtenu par l’opération de soustraction entre la valeur de contraste de tavelure calculée pour la cible de référence et chaque valeur de contraste de tavelure relative à la cible à analyser. Ainsi, l’image de variation de contraste de tavelure qui est obtenue pour la cible à analyser met en évidence directement l’amplitude des mouvements ou micromouvements qui se produisent dans la cible à analyser, indépendamment de la source laser qui est utilisée. Autrement dit, l’invention fournit un résultat d’analyse absolu, ou “universel”, qui permet de comparer entre elles des analyses effectuées avec des sondes optiques différentes, et à des moments séparés. Il est possible en outre de réduire ou supprimer en plus les effets du bruit de détection d’image, notamment du bruit de grenaille («shot noise» en anglais), en corrigeant préalablement chacune des valeurs de contraste de tavelure relatives à la cible à analyser et chaque valeur de contraste de tavelure C static qui est calculée pour la cible de référence, d’une contribution des fluctuations de bruit de détection d’image. La valeur du contraste de tavelure C static qui est calculée à partir de la partie de chaque image saisie correspondant à la cible de référence, et qui est corrigée ainsi, est essentiellement dépourvue d’effet de bruit de détection d’image. Elle permet alors une caractérisation exacte de la source laser d’éclairage. De façon générale pour l’invention, la cible de référence peut être un morceau de papier à valeur identifiée de granulométrie. Elle peut être maintenue par un support qui est proche du plan d’objet du système imageur, aussi appelé plan de focalisation de ce système, à l’intérieur du champ de vue, de préférence à proximité d’un bord de ce champ de vue. La cible de référence est ainsi systématiquement jointe à la cible à analyser pour chaque image qui est saisie. Aussi de façon générale pour l’invention, à l’étape /1/, la cible de référence peut être sélectionnée préférablement de sorte qu’elle possède un niveau d’intensité de diffusion pour le rayonnement de la source laser, tel que chaque grain de tavelure dans la partie de chaque image saisie qui correspond à cette cible de référence, soit à l’intérieur d’un intervalle de détection linéaire de l’instrument imageur. Chaque valeur de contraste de tavelure C static qui est calculée à l’étape /3/ est ainsi plus significative. Encore de façon générale pour l’invention, l’instrument imageur peut comprendre : - une optique de formation d’image ; - un capteur d’images matriciel, comprenant des photodétecteurs qui sont disposés à des intersections de lignes et de colonnes ; et - un diaphragme pupillaire. Alors, une ouverture du diaphragme pupillaire peut être ajustée préalablement à la saisie des images à l’étape /2/, de sorte que chaque grain de tavelure recouvre au moins deux photodétecteurs dans la partie de chaque image saisie qui correspond à la cible de référence. Il est ainsi possible de bien discrétiser spatialement toutes les teintes des grains de tavelure dans chaque image qui est saisie. Encore de façon générale pour l’invention, la cible de référence peut être sélectionnée à l’étape /1/ préférablement pour posséder une taille et des caractéristiques de rugosité telles que la partie de chaque image saisie qui correspond à cette cible de référence contienne au moins cinquante grains de tavelure. Dans ces conditions, chaque valeur de contraste de tavelure C static qui est calculée à l’étape /3/ peut posséder un niveau de fiabilité élevé. Encore de façon générale pour l’invention, le procédé peut comprendre additionnellement de caractériser une répartition spectrale du rayonnement qui est produit par la source laser, telle que cette répartition spectrale est effective pour au moins une des images saisies. Cette caractérisation de répartition spectrale est déduite du motif de tavelure dans la partie de l’image saisie qui correspond à la cible de référence. Une telle caractérisation de répartition spectrale peut notamment comprendre de déterminer une largeur spectrale de la source laser utilisée. Le procédé de l’invention peut être utilisé pour de nombreuses applications. En particulier, la cible à analyser peut être une portion d’un fluide qui est partiellement transparent pour le rayonnement de la source laser et qui contient des particules diffusantes mobiles. Pour de telles applications, qui peuvent notamment concerner des unités de purification d’eau, le procédé peut comprendre de caractériser une combinaison d’une concentration des particules dans le fluide avec une vitesse de déplacement de ces particules. Alternativement, la cible à analyser peut être une portion d’un solide diffusant qui vibre. Alors, le procédé peut comprendre de caractériser des amplitudes locales de vibration qui existent en des endroits distincts dans la portion de solide. Pour d’autres applications qui sont aussi possibles, la cible à analyser peut être une portion d’un tissu biologique. Dans ce cas, le procédé peut comprendre de caractériser des niveaux de mouvements qui existent en des endroits distincts dans le tissu biologique. Ces mouvements peuvent indiquer des niveaux de vascularisation qui varient entre différentes parties du tissu biologique, ou entre des dates différentes auxquelles le procédé de l’invention est répété pour un même endroit du tissu biologique. Enfin, un second aspect de l’invention propose un système d’imagerie par visualisation de variations d’un contraste de tavelure, ce système comprenant : - un instrument imageur, qui est adapté pour saisir des images d’une cible à analyser ; - une source laser, qui est disposée pour éclairer la cible à analyser pendant que chaque image est saisie, de sorte que l’image soit formée par du rayonnement produit par la source laser puis diffusé par la cible à analyser ; - une unité de calcul, qui est configurée pour calculer des valeurs d’un contraste de tavelure relatives à la cible à analyser, à partir d’une ou plusieurs images de celle-ci ; - un système d’affichage, connecté pour montrer une image qui a été construite par l’unité de calcul, - une cible de référence, qui possède un motif diffusant statique et constant ; et - des moyens de support, adaptés pour maintenir la cible de référence à l’intérieur d’un champ de vue de l’instrument imageur qui est effectif pour chaque image, et de sorte que la source laser éclaire simultanément la cible à analyser et la cible de référence pendant que chaque image est saisie. Selon l’invention, l’unité de calcul est configurée en outre pour calculer, séparément pour chaque image saisie, une valeur de contraste de tavelure, notée C static , à partir d’une partie de l’image qui correspond à la cible de référence. Elle est configurée en plus pour valider plusieurs images qui ont été saisies successivement si les valeurs C static qui ont été calculées respectivement pour ces images ont des variations de l’une quelconque des images à une autre qui sont inférieures à un seuil. L’unité de calcul est en plus configurée pour construire une image de variation de contraste de tavelure en attribuant à chacun d’un ensemble de points d’image qui correspondent à la cible à analyser, une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies, et d’autre part une valeur de contraste de tavelure qui a été calculée pour le point d’image à partir d’une ou plusieurs images saisies puis validées. Un tel système d’imagerie est adapté pour mettre en œuvre un procédé conforme au premier aspect de l’invention, y compris ses perfectionnements et/ou conditions préférées de mise en œuvre. Optionnellement, le système d’imagerie par visualisation de variations de contraste de tavelure de l’invention peut comprendre en outre un dispositif supplémentaire de saisie d’images de la cible à analyser, qui est adapté pour fournir des images en couleurs. Avantageusement, ce dispositif supplémentaire peut permettre d’éclairer la cible à analyser successivement en trois couleurs séparées, et peut utiliser le même instrument imageur pour saisir des images séparées de chaque couleur, que celui utilisé pour les images saisies avec éclairage par la source laser. Optionnellement aussi, l’instrument imageur peut comprendre en outre un polariseur linéaire qui est orienté perpendiculairement à une direction de polarisation de la source laser. Un tel polariseur linéaire, orienté de cette façon, améliore les images de variation de contraste de tavelure qui sont obtenues en supprimant la partie du rayonnement d’éclairage qui est diffusée par la surface de la cible à analyser. Brève description des figures Les caractéristiques et avantages de la présente invention apparaîtront plus clairement dans la description détaillée ci-après d’exemples de mise en œuvre non-limitatifs, en référence aux figures annexées parmi lesquelles : FIG. 1 est un schéma optique d’un système d’imagerie qui est conforme à l’invention ; FIG. 2 montre schématiquement un contenu d’une image telle que saisie par le système d’imagerie de FIG. 1 ; FIG. 3 montre des étapes d’un procédé d’imagerie qui est conforme à l’invention ; FIG. 4 montre un exemple d’image de tissu organique éclairé par une source laser, telle que saisie par le système d’imagerie de FIG. 1 ; FIG. 4 montre une image de variation de contraste spatial de tavelure qui a été construite selon le procédé de l’invention, à partir d’une série d’images successives chacune saisie comme celle de FIG. 4 ; FIG. 5 montre un exemple d’image en couleurs de tissu organique telle que saisie par le système d’imagerie de FIG. 1 , en utilisant un mode d’obtention d’images en couleurs tel que connu de l’art antérieur ; et FIG. 5 montre une image de variation de contraste temporel de tavelure qui a été construite selon le procédé de l’invention, à partir d’une série d’images saisies successivement avec éclairage laser. Description détaillée de l’invention Pour raison de clarté, les dimensions des éléments qui sont représentés dans FIG. 1 et FIG. 2 ne correspondent ni à des dimensions réelles, ni à des rapports de dimensions réels. En outre, certains des éléments représentés dans FIG. 1 ne le sont que symboliquement. Conformément à FIG. 1 , un système d’imagerie comprend une sonde optique 10, une unité de calcul 20 et un écran d’affichage 30. L’unité de calcul 20 est notée COMPUT. et peut être un ordinateur portable. L’écran d’affichage 30, noté DISPLAY, est de préférence d’un type à haute définition. La sonde optique 10 peut être d’un modèle portatif, avec un écran latéral («baffle» en anglais) 11 et une poignée de préhension 12. Elle regroupe une voie d’éclairage et une voie de détection, qui sont agencées pour que la voie d’éclairage éclaire tout un champ de vue de la voie de détection. Pour cela, une direction d’éclairage de la voie d’éclairage peut être oblique par rapport à un axe optique de la voie de détection, noté A-A, afin de juxtaposer les deux voies d’éclairage et de détection à l’intérieur de la sonde optique 10. Le champ de vue est entouré par l’écran latéral 11 pour supprimer d’éventuels rayons de lumière parasite qui ne proviendraient pas de la voie d’éclairage ni du champ de vue. La voie d’éclairage comprend une source laser 1 et une optique divergente 2 par laquelle du rayonnement qui est émis par la source laser 1 éclaire l’ensemble du champ de vue, ou une partie essentielle de ce champ de vue qui est utile pour mettre en œuvre l’invention. La source laser 1 peut être d’un type qui produit du rayonnement à une longueur d’onde adaptée à l’application envisagée. Par exemple, pour observer des déplacements de globules rouges dans un tissu organique, la longueur d’onde de la source laser 1 peut être dans l’intervalle de 785 nm (nanomètre) à 810 nm. Cet intervalle est pertinent car le rayonnement de la source laser 1 pénètre alors à travers de la peau qui peut être présente à la surface du tissu organique. Pour la détection de vibrations de métaux, des valeurs plus faibles de longueur d’onde dans le domaine visible, par exemple 512 nm, 630 nm..., sont aussi adaptées. Une largeur spectrale de la source laser 1 est de préférence inférieure à 1 nm, d’une part pour avoir une meilleure sensibilité de mesure, et d’autre part pour que l’invention fournisse des images de variation de contraste de tavelure qui soient indépendantes de la source laser utilisée dans une mesure suffisante. La source laser 1 peut être dépourvue de moyens dédiés à stabiliser son fonctionnement d’émission, si bien qu’elle peut être peu onéreuse, de faible encombrement et légère, notamment pour être incorporée dans la sonde optique portative 10. La voie de détection comprend un instrument imageur 4, qui est constitué d’une optique de formation d’image 41, d’un capteur d’images matriciel 42, d’un diaphragme pupillaire 43, aussi appelé diaphragme d’ouverture, d’un diaphragme de champ optionnel 44 et d’un polariseur linéaire optionnel 45. Le capteur d’images matriciel 42, noté DETECT., comporte de préférence une densité élevée de photodétecteurs séparés dans sa surface photosensible S. Une sortie de signaux de lecture du capteur d’images matriciel 42 est connectée à une entrée de l’unité de calcul 20, et une sortie vidéo de cette dernière est connectée à l’écran d’affichage 30. Un intervalle de sensibilité du capteur d’images matriciel 42 contient la valeur de longueur d’onde de la source laser 1. Le cas échéant, le diaphragme de champ 44 peut être situé contre la surface photosensible S du capteur d’images matriciel 42. De façon connue, la surface photosensible S est formée par les photodétecteurs qui sont disposés à des intersections de lignes et de colonnes d’un agencement matriciel. La durée d’intégration qui est mise en œuvre par le capteur 42 pour saisir chaque image produite en utilisant un éclairage par la source laser 1, peut être comprise entre 1 ms (milliseconde) et 100 ms, à titre d’exemple. Etant donné que le rayonnement d’éclairage qui est produit par la source laser 1 est polarisé rectilignement, le polariseur linéaire 45 est avantageusement orienté perpendiculairement à la direction de polarisation du rayonnement de la source laser 1. De cette façon, le polariseur 45 supprime une partie du rayonnement laser qui est diffusée par le contenu du champ de vue. Il favorise en conséquence la détection d’une autre partie du rayonnement laser qui subit des diffusions multiples produites par un contenu diffusant du champ de vue, comme décrit dans FR 3 062 542 A1. Selon un perfectionnement facultatif de l’invention, la sonde optique 10 peut être complétée pour saisir des images en couleurs du contenu du champ de vue. Pour cela, le capteur d’images matriciel 42 peut être d’un type monochrome sensible sur tout le domaine visible, et des sources de lumière rouge, verte et bleue peuvent être ajoutées dans la sonde optique 10. Par exemple, des diodes électroluminescentes rouges, vertes et bleues 3 peuvent être juxtaposées de façon répétée le long d’un anneau qui entoure le champ de vue au niveau du diaphragme pupillaire 43, autour de l’ouverture de ce diaphragme comme cela apparaît dans FIG. 1 . Les diodes électroluminescentes 3 produisent un éclairage non-cohérent du contenu du champ de vue pour chacune des couleurs bleue, verte et rouge. Alors, trois images qui sont saisies successivement avec des durées intermédiaires minimales par rapport à la durée de cycle de fonctionnement du capteur d’images 42, l’une des images en utilisant un éclairage exclusivement par les diodes électroluminescentes bleues, une autre image avec uniquement les diodes électroluminescentes vertes, et la dernière image avec uniquement les diodes électroluminescentes rouges, permet de reconstruire une image en couleurs RVB en n’utilisant que le capteur d’images 42 de type monochrome. Ainsi, l’éclairage par la source laser 1 permet de saisir des images de tavelure pour mettre en œuvre l’invention, et la façon qui vient d’être décrite pour acquérir des images en couleurs RVB fournit de telles images en couleurs qui peuvent être superposées aux images de tavelure. La référence 13 désigne une gâchette de déclenchement de saisie d’images. Un appui court sur la gâchette 13 par un opérateur déclenche la saisie d’images successives avec éclairage par la source laser 1 uniquement, suivi du calcul et de l’affichage en temps réel d’images de variation de contraste de tavelure successives, par exemple avec une fréquence de rafraichissement qui est comprise entre 5 Hz (hertz) environ et 25 Hz environ en fonction d’une valeur voulue pour un rapport signal-sur-bruit de chaque image de variation de contraste de tavelure. Les images de variation de contraste de tavelure qui sont affichées ainsi sont conformes à l’art antérieur tel qu’existant avant la présente invention, tout en bénéficiant cependant du filtrage polarimétrique qui procure une meilleure sensibilité à une partie du rayonnement qui a pénétré plus profondément dans la cible à analyser 100. Un appui long qui est effectué alternativement par l’opérateur sur la gâchette 13, par exemple avec une durée d’appui supérieure à 1 s (seconde), déclenche la saisie d’images successives dont les trois premières sont saisies avec des éclairages respectifs par les diodes électroluminescences 3, séparément et successivement pour les couleurs bleue, verte et rouge, puis les images suivantes sont saisies en utilisant pour chacune un éclairage par la source laser 1. Ces images suivantes saisies avec l’éclairage laser sont utilisées pour calculer une seule image de variation de contraste de tavelure. Celle-ci possède une valeur de rapport signal-sur-bruit qui est supérieure à celle de chaque image de variation contraste de tavelure d’un appui court sur la gâchette 13, grâce à un nombre plus grand d’images saisies qui sont combinées pour obtenir chaque image de variation de contraste de tavelure. Chaque appui long sur la gâchette 13 aboutit ainsi à une image en couleurs RVB, telle que connue usuellement, et à une image de variation de contraste de tavelure à valeur élevée de rapport signal-sur-bruit, qui sont toutes les deux relatives au même contenu de champ de vue et donc superposables. L’invention qui est décrite en détail dans la suite peut être appliquée aux images de variation de contraste de tavelure d’un appui court ou d’un appui long sur la gâchette 13. Une cible à analyser est désignée par la référence 100. Elle est semi-transparente ou rugueuse, et diffusante pour le rayonnement de la source laser 1. Ce peut être un milieu biologique, tel qu’une partie de tissu d’organe ou autre. Lorsque la sonde optique 10 est utilisée pour saisir des images de la cible à analyser 100, le bord antérieur de l’écran latéral 11 peut être appliqué sur une surface de la cible à analyser 100, autour de la partie de celle-ci qui est à analyser. Le rayonnement laser de la source 1 pénètre dans la cible à analyser 100 sur une profondeur qui dépend de la longueur d’onde de la source laser 1 et de la nature de la cible à analyser 100. Cette profondeur de pénétration du rayonnement laser dans la cible à analyser peut être de l’ordre de quelques dizaines de nanomètres, notamment pour des métaux rugueux qui constituent la cible à analyser, à quelques millimètres, notamment pour des milieux biologiques, tels que des tissus cutanés. Le rayonnement laser est rétrodiffusé par des centres de diffusion qui sont présents dans la cible 100, et des parties rétrodiffusées du rayonnement laser interfèrent pour créer un motif de tavelure dans le plan de focalisation du système imageur 4. L’optique de formation d’image 41 est telle que la surface de la cible à analyser 100 soit sensiblement conjuguée optiquement avec la surface photosensible S du capteur d’images matriciel 42. Dans ces conditions d’éclairage et de saisie d’image, l’image qui est saisie par ce système imageur 4 reproduit alors le motif de tavelure. Ce motif est constitué de grains de tavelure qui sont juxtaposés, en étant alternativement plus clairs et plus foncés. Pour mettre en œuvre l’invention, une cible de référence 101 est ajoutée dans le champ de vue de l’instrument imageur 4, de façon à être éclairée par la source laser 1 en même temps que la cible à analyser 100, et à apparaître avec cette dernière dans chaque image qui est saisie. De préférence, la cible de référence 101 est placée dans le plan de focalisation de l’instrument imageur 4, ou à proximité de ce plan de focalisation. Elle est donc ainsi sensiblement au même niveau que la surface de la cible à analyser 100 le long de l’axe optique A-A de la voie de détection de l’instrument imageur 4. Dans ces conditions, le contenu de chaque image qui est saisie est tel que représenté dans FIG. 2 . Il comprend une partie d’image 101’ qui est conjuguée optiquement avec la cible de référence 101 par l’optique de formation d’image 41, et le reste de l’étendue de chaque image, désignée par la référence 100’, est conjugué optiquement avec la partie visée de la cible à analyser 100. La cible de référence 101 peut être maintenue fixement à la position voulue dans le champ de vue en étant accrochée à l’écran latéral 11. Dans ce cas, l’écran latéral 11 réalise la fonction supplémentaire de moyens de support tels qu’introduits dans la partie générale de la présente description. La cible de référence 101 possède un motif qui est diffusant pour le rayonnement de la source laser 1, et qui est statique et constant. Notamment, la cible de référence 101 peut être constituée par une portion d’un matériau solide et opaque qui présente une rugosité à sa surface. Par exemple, ce peut être un morceau de feuille de papier commercialisé sous la marque CANSON ® . La rugosité de cette cible de référence 101, telle que caractérisée par une valeur d’écart-type de la distribution de hauteur de surface de la cible 101, est de préférence supérieure à la valeur de longueur d’onde du rayonnement laser qui est utilisé, et une largeur moyenne de relief et creux de la rugosité est de préférence telle qu’au moins dix reliefs séparés de rugosité, de préférence cinquante, soient éclairés par la source laser 1. Dans ces conditions, le rayonnement qui est rétrodiffusé par la cible de référence 101 produit des tavelures qui sont conformes à des statistiques circulaires Gaussiennes, et qui présentent un contraste de tavelure maximum. Alors, la partie d’image 101’ est aussi constituée par un motif de tavelure, qui est utilisé selon l’invention d’une part pour calculer une valeur de contraste statique relative à la cible de référence 101, et d’autre part pour vérifier a posteriori la stabilité du fonctionnement de la source laser 1. Le motif de rugosité de la cible de référence 101 étant constant et dépourvu de mouvements internes, des variations du motif de tavelure dans la partie d’image 101’, entre des images qui ont été saisies successivement, sont nécessairement causées soit par des mouvements de torsion excessifs de la sonde optique lorsqu’elle est en contact avec la cible à analyser, soit par des variations de fonctionnement de la source laser 1, ce qui dans tous les cas invalide la mesure. Lorsque plusieurs images sont saisies successivement avec éclairage par la source laser 1, et que ces images saisies présentent entre elles des variations de contraste spatial de tavelure dans la partie d’image 101’ qui sont inférieures à un seuil, en valeurs absolues, alors le procédé de l’invention propose de valider ces images saisies en considérant que le fonctionnement de la source laser 1 a été constant pendant la durée d’intégration de chacune d’elles et a été identique pour toutes ces images saisies. Par fonctionnement constant de la source laser 1, on entend une émission laser qui est notamment : dépourvue d’effets de réinjection de rayonnement rétrodiffusé variable dans la source laser, dépourvue d’effets de variations de bruit d’alimentation électrique, en amplitude et en phase, de la source laser, dépourvue d’effets de variations thermiques qui pourraient affecter certains composants de la source laser, et dépourvue de fluctuations temporelles dans la répartition spectrale du rayonnement qui est émis par la source laser, qui pourraient être dues à des transitions non-contrôlées entre une émission monomode voulue et une émission multimode non-voulue, ou qui pourraient apparaître sous forme de variations temporelles de la largeur spectrale du rayonnement laser émis. Pour mettre en œuvre l’invention de manière optimale, les conditions opératoires suivantes peuvent être adoptées : - la cible de référence 101 est sélectionnée pour posséder un pouvoir de diffusion à la longueur d’onde de la source laser 101, tel que les variations spatiales d’intensité lumineuse du motif de tavelure soient à l’intérieur d’un intervalle de détection linéaire du capteur d’images matriciel 42. Lorsque la cible de référence 101 est un morceau de feuille de papier CANSON ® , la couleur de papier peut notamment être choisie pour ne pas correspondre à la longueur d’onde de la source laser 101 ; - l’ouverture du diaphragme pupillaire 43 est ajustée pour que chaque grain de tavelure dans la partie d’image 101’ recouvre au moins deux photodétecteurs voisins dans la surface photosensible S du capteur d’images matriciel 42 ; et - la cible de référence 101 est sélectionnée pour posséder une taille et des caractéristiques de rugosité telles que la partie d’image 101’ contienne au moins cinquante grains de tavelure. Lorsque les conditions mentionnées plus haut sont remplies, pour que le motif de tavelure qui est saisi dans la partie d’image 101’ soit proche d’un motif circulaire Gaussien, tel que défini dans le livre “Speckle Phenomena in Optics, Theory and Applications”, de J.W. Goodman, et dans l’article intitulé «Study of a circular Gaussian transition in an optical speckle field», de I. Bergoënd, X. Orlik et Eric Lacot, Journal of the European Optical Society, Rapid Publications 3, 08028 (2008), il est possible de déduire des valeurs numériques pour certaines caractéristiques de la répartition spectrale du rayonnement qui est émis par la source laser 1, à partir du motif de tavelure qui est saisi dans la partie d’image 101’. La possibilité d’une telle analyse est montrée notamment dans l’article intitué «Choosing a laser for laser speckle contrast imaging» de D.D. Postnov, X. Cheng, S.E. Erdeneret D.A.Boas, Scientific Reports, 22 février 2019,9(1):2542. doi: 10.1038/s41598-019-39137-x. Lorsque l’émission laser est monomode, des valeurs respectives de la longueur d’onde centrale d’émission et de la largeur spectrale d’émission peuvent être obtenues de cette façon. Lorsqu’une telle déduction de valeurs numériques qui caractérisent l’émission laser est répétée séparément pour plusieurs images qui ont été saisies successivement, à partir des parties d’images 101’, les résultats obtenus caractérisent des fluctuations du fonctionnement de la source laser 1. Ces fluctuations sont acceptables lorsque les valeurs numériques qui sont obtenues pour chaque caractéristique considérée de la répartition spectrale du rayonnement de la source laser 1 présentent des variations inférieures à un seuil fixé. Dans ce cas, les images sont validées pour analyser la cible 100. Dans le contraire, c’est-à-dire lorsque la source laser 1 a subi des fluctuations de fonctionnement qui apparaissent trop importantes d’après les valeurs numériques obtenues successivement pour la répartition spectrale de son rayonnement, les images sont déclarées inacceptables pour analyser valablement la cible 100. Selon l’invention, il est suffisant de vérifier qu’un contraste de tavelure à l’intérieur de la partie d’image 101’ varie dans une mesure qui est suffisamment faible entre toutes les images de la série, pour accepter ou rejeter ces images aux fins d’analyse de la cible 100. Ce contraste de tavelure, noté C static , est calculé séparément pour chaque image saisie selon la formule : C static = σ / , où et σ désignent respectivement la valeur moyenne et l’écart-type des valeurs d’intensité saisies dans la partie d’images 101’ de chaque image saisie. En référence à FIG. 3 , le procédé d’imagerie de l’invention pour analyser la cible 100 comprend les étapes suivantes : à l’étape S1 : une série de plusieurs images successives de la cible à analyser 100, de deux images à quelques centaines d’images, est saisie en utilisant la sonde optique 10 munie de la cible de référence 101, avec éclairage par la source laser 1. Cette cible de référence 101 est assemblée dans la sonde optique 10 identiquement pour toutes les images de la série. De préférence, la cible de référence 101 est laissée dans la sonde optique 10 sans être démontée ni déplacée entre deux images successives ; à l’étape S2 : les valeurs du contraste de tavelure C static sont calculées séparément pour chaque image de la série, à partir de la partie 101’ de cette image ; à l’étape S3 : l’existence de fluctuations dans le fonctionnement de la sonde laser 1, qui se seraient produites entre certaines des images de la série, est recherchée. Pour cela, des différences de valeurs du contraste de tavelure C static , notées ΔC static , sont calculées entre deux des images de la série, et comparées à un seuil ΔC max , pour certaines ou toutes les paires d’images dans la série. Si l’une au moins de ces différences excède le seuil ΔC max , la série d’images est rejetée et le procédé d’imagerie est repris à l’étape S1. Sinon, la série d’images est validée et le procédé d’imagerie est poursuivi avec l’étape S4 ; à l’étape S4 : pour chaque point d’image de la partie d’image 100’, qui était conjuguée optiquement avec la cible à analyser 100 pendant que chaque image de la série a été saisie, une valeur de contraste C cible (i, j) est calculée pour ce point d’image conformément au mode choisi pour la construction d’une image de variation de contraste de tavelure, où i et j sont les coordonnées du point d’image dans la matrice de la surface photosensible S du capteur d’images matriciel 2 ; à l’étape S5 : construction d’une nouvelle image, qui est l’image de variation de contraste de tavelure obtenue conformément à l’invention, en attribuant au point d’image de coordonnées i et j dans la partie d’image 100’ une nouvelle valeur de point d’image notée ΔC(i, j) et calculée selon la formule suivante : ΔC(i, j) = static > - C cible (i, j), où static > est la moyenne des valeurs respectives du contraste de tavelure C static obtenues à l’étape S2 pour toutes les images de la série, pour les trois modes de construction d’image de variation de contraste de tavelure indiqués au début de la présente description. Les nouvelles valeurs de points d’image qui sont calculées ainsi, ΔC(i, j), sont toujours positives ; et à l’étape S6 : affichage de l’image de variation de contraste de tavelure sur l’écran 30. Cette image de variation de contraste de tavelure met en évidence les parties de la cible à analyser qui ont subi des mouvements, notamment des mouvements à l’échelle submicronique. Elle possède des valeurs de points d’image qui dépendent de l’amplitude et/ou de la vitesse de ces mouvements. En outre, grâce à l’opération de soustraction dans le calcul de l’étape S5 pour obtenir chaque valeur de point d’image ΔC(i, j), des images de variation de contraste de tavelure qui sont élaborées à des dates différentes, ou avec des sondes optiques 10 différentes, peuvent être comparées entre elles valablement. La valeur de contraste de tavelure C static qui est calculée à l’étape S2 comprend une contribution de tavelure qui résulte du bruit de détection survenant dans le capteur d’images matriciel 42. Ce bruit de détection est une combinaison de plusieurs types de bruit, dont les principaux sont le bruit de grenaille et le bruit d’obscurité. Le bruit d’obscurité est le plus souvent négligeable pour un fonctionnement usuel du système d’imagerie, en raison du nombre suffisant de photons qui sont émis par la source laser 1 puis rétrodiffusés par la cible de référence 101 vers le capteur d’images 42. Par contre, le bruit de grenaille, ou «shot noise» an anglais, qui est proportionnel à la racine carrée de l’intensité saisie I, peut contribuer notablement à la valeur d’écart-type σ, et intervient donc dans chaque valeur de contraste de tavelure. La valeur d’écart-type d’intensité qui résulte du bruit de grenaille peut donc avantageusement être soustraite à la valeur d’écart-type d’intensité qui est calculée pour la succession des images saisies à l’étape S1 dans les parties d’images 101’ pour chaque point d’image dans le cas du contraste temporel ou pour chaque sous-ensemble de points d’image dans le cas du contraste spatial ou le cas du contraste spatio-temporel. Cette opération de soustraction, pour supprimer la contribution du bruit de grenaille, est de préférence appliquée au contraste de tavelure C static pour la cible de référence 101 et pour chaque valeur de contraste C cible (i, j) pour la cible à analyser 101. L’autre opération de soustraction qui est réalisée à l’étape S5 pour calculer chaque valeur de point d’image ΔC(i, j) permet d’attribuer directement, sous réserve des conditions mentionnées plus haut, les variations des valeurs ΔC(i, j) aux mouvements qui surviennent dans la cible à analyser 100, malgré l’utilisation de sources laser différentes et de systèmes optiques différents entre des exécutions séparées du procédé de l’invention. Il est à noter qu’un même micromouvement dans la cible à analyser 100 se traduira par des valeurs de variation de contraste de tavelure ΔC(i, j) qui sont différentes si la valeur de contraste de tavelure C static calculée pour la cible de référence 101 n’est pas la même pour les sources laser respectives des deux sondes optiques. En effet, la baisse de contraste ne varie pas linéairement avec la vitesse des micromouvements, et dépend donc de la valeur du contraste de tavelure qui existerait en l’absence des micromouvements. Cette non-linéarité est cependant connue et publiée, notamment dans l’article «Review of laser speckle contrast techniques for visualizing tissue perfusion», de Matthijs Draijer et al., Lasers Med. Sci. (2009) 24:639-651, FIG. 2 , p. 641, ce qui permet une comparaison rigoureuse des résultats malgré la différence de valeur de contraste de tavelure C static entre les deux sources laser utilisées, ou plus généralement entre les deux sondes optiques utilisées. Les étapes S2 à S5 sont exécutées par l’unité de calcul 20, de préférence d’une façon automatisée. Avec le premier mode de construction d’image de contraste de tavelure, des valeurs de contraste spatial sont calculées pour chaque point d’image (i, j) dans la partie d’image 100’, séparément pour chaque image qui a été saisie à l’étape S1. Elles sont calculées chacune sur un voisinage du point d’image de coordonnées i et j dans la partie d’image 100’, ce voisinage pouvant être un carré de 5 x 5, ou 7 x 7, ou 9 x 9 points d’image centré sur le point d’image i, j. Chaque valeur de contraste spatial est alors égale au quotient de l’écart-type des valeurs d’intensité saisies à l’intérieur de ce carré de voisinage, par la moyenne de ces mêmes valeurs d’intensité saisies encore à l’intérieur du carré de voisinage. Ces valeurs de contraste spatial sont ensuite moyennées sur toutes les images saisies, séparément pour chaque point d’image, afin de renforcer le signal obtenu. Les valeurs moyennes qui sont calculées ainsi sont celles de l’étape /4/ notées C cible (i, j). Pour ce premier mode de construction d’image de contraste de tavelure, l’image qui est affichée à l’étape S6 est donc une image de variation de contraste spatial de tavelure. L’existence de mouvements à l’intérieur de la cible à analyser 100 provoque des réductions locales des valeurs C cible (i, j), et par conséquent des augmentations locales des valeurs ΔC(i, j) dans l’image de variation de contraste de tavelure qui est affichée à l’étape S6, puisque C static correspond à la valeur maximale possible pour C cible (i, j). Avec le deuxième mode de construction d’image de contraste de tavelure, les valeurs de contraste C cible (i, j) qui sont calculées à l’étape S4 sont des valeurs de contraste temporel qui sont calculées respectivement pour les points d’image de la partie d’image 100’. Chaque valeur de contraste C cible (i, j) est alors égale au quotient de l’écart-type des valeurs d’intensité de toutes les images de la série au point d’image de coordonnées i et j, par une moyenne de ces mêmes valeurs d’intensité. L’image qui est affichée à l’étape S6 est alors une image de variation de contraste de tavelure entre le contraste spatial calculé sur la cible de référence 101 et le contraste temporel de tavelure calculé sur la cible à analyser 100. L’existence de mouvements à l’intérieur de la cible à analyser 100 qui se produisent entre des images saisies successivement provoque alors, selon la durée d’intégration du capteur d’images matriciel 42 par rapport au temps caractéristique qui résulte de la vitesse des mouvements dans la cible à analyser 100, des baisses ou des augmentations locales du contraste temporel de tavelure qui vont donc provoquer respectivement des augmentations ou des baisses des valeurs ΔC(i, j), visibles dans l’image de variation de contraste de tavelure affichée à l’étape S6. Avec le troisième mode de construction d’image de contraste de tavelure, les valeurs de contraste C cible (i, j) qui sont calculées à l’étape S4 sont des valeurs de contraste spatio-temporel qui sont calculées respectivement pour les points d’image de la partie d’image 100’. Par exemple, chaque valeur de contraste C cible (i, j) peut être calculée selon la même formule que précédemment pour le contraste temporel, mais où l’intensité de chaque point d’image dans la partie d’image 100’ est remplacée par une valeur moyenne d’intensité calculée sur un environnement de ce point d’image, par exemple un carré de dimensions 3 x 3, 5 x 5,…. L’image qui est affichée à l’étape S6 est alors une image de variation de contraste spatio-temporel de tavelure, calculée comme la différence entre le contraste spatial de tavelure de la cible de référence 101 et le contraste spatio-temporel de tavelure calculé sur la cible à analyser 100. L’existence de mouvements à l’intérieur de la cible à analyser 100 qui se produisent entre des images saisies successivement provoque alors, selon la durée d’intégration du capteur d’images matriciel 42 par rapport au temps caractéristique qui résulte de la vitesse des mouvements dans la cible à analyser 100, des baisses ou des augmentations locales de contraste spatio-temporel qui vont donc provoquer respectivement des augmentations ou des baisses des valeurs ΔC(i, j), comme pour le cas du contraste temporel. De façon générale pour tous les modes de construction de l’image de variation de contraste de tavelure, c’est la variation ΔC(i, j) à chaque point d’image qui est affichée à l’étape S6. Pour les cas de contraste temporel et de contraste spatio-temporel, des valeurs de paramètres opératoires pour la sonde optique 10 peuvent être adoptés de préférence, notamment pour la durée d’intégration du capteur d’images matriciel 42, qui maximisent l’écart entre la valeur de contraste spatial de tavelure static > obtenue pour la cible de référence 101 et les valeurs de contraste temporel ou spatio-temporel de tavelure C cible (i, j) obtenues pour la cible à analyser 100. Lorsque la cible à analyser 100 est une portion d’un fluide partiellement transparent, les mouvements dans celui-ci qui sont mis en évidence par l’image de variation de contraste de tavelure correspondent à des endroits où une vitesse et/ou une concentration de particules diffusantes mobiles dans le fluide est (sont) plus importante(s). Lorsque la cible à analyser 100 est une portion d’un solide diffusant qui vibre, les mouvements dans celui-ci qui sont mis en évidence par l’image de variation de contraste de tavelure correspondent à des endroits de maximum d’amplitude de vibration, couramment appelés ventres de vibration ou «vibration antinodes» en anglais. Lorsque la cible à analyser 100 est un tissu biologique vascularisé, que la longueur d’onde de la source laser 1 se situe dans le domaine spectral proche infrarouge, et que le temps d’intégration du capteur d’images matriciel 42 se situe entre 1 ms et quelques dizaines de millisecondes, les mouvements dans celui-ci qui sont mis en évidence par l’image de variation de contraste de tavelure sont majoritairement produits par des inhomogénéités de vascularisation. Ces inhomogénéités concernent la vitesse de déplacement des globules rouges et/ou la concentration des globules rouges mobiles. FIG. 4 reproduit une image d’un tissu organique, telle que saisie par le système d’imagerie de FIG. 1 , en utilisant l’éclairage par la source laser 1 et une durée d’intégration de 10 ms dans le capteur d’images matriciel 42. Cette image n’a subi aucun traitement, si bien qu’elle correspond à la désignation usuelle d’image brute. Les parties d’image 100’ et 101’ qui correspondent respectivement à la cible à analyser 100 et à la cible de référence 101 sont indiquées. Chacune de ces parties d’image est constituée par un motif de tavelure, dont les origines respectives ont été expliquées plus haut. FIG. 4 montre l’image de variation de contraste spatial de tavelure qui est construite selon l’invention à partir d’une série d’images successives saisies de la même façon que l’image de FIG. 4 , et en appliquant le premier mode de construction d’image de variation de contraste de tavelure. Cette image de variation de contraste spatial de tavelure révèle des micromouvements qui se produisent dans un réseau de micro-vascularisation du tissu organique. FIG. 5 reproduit une image en couleurs RVB de la pulpe d’extrémité d’un doigt humain, telle que fournie par le système d’imagerie de FIG. 1 en utilisant successivement un éclairage par les diodes électroluminescentes 3 bleues, puis par celles vertes, puis celles rouges. L’image en couleurs de FIG. 5 a été reconstruite à partir des trois images saisies séparément pour chaque couleur. Elle fait principalement apparaître la surface de la peau avec le motif d’empreinte digitale. Elle est destinée à être superposée à une image de variation de contraste de tavelure obtenue selon la présente invention. Pour cela, une série d’une centaine d’images successives est alors saisie en utilisant l’éclairage par la source laser 1, sans déplacer la sonde optique 10 par rapport au doigt, par exemple à raison d’une centaine d’images par seconde avec une durée d’intégration de 10 ms pour chaque image. FIG. 5 montre l’image de variation de contraste temporel de tavelure qui est construite selon l’invention à partir de cette série d’images avec éclairage laser, en appliquant le deuxième mode de construction d’image de variation de contraste de tavelure. Elle révèle les micromouvements du réseau microvasculaire sous-cutané et/ou transcutané. L’image de FIG. 5 est la superposition de l’image de variation de contraste temporel avec l’image en couleurs de FIG. 5 . Enfin, et comme de façon générale en imagerie de tavelure dynamique, il est possible d’utiliser l’indice de flux sanguin, ou BFI pour «Blood Flux Index» en anglais, qui est égal à 1/C cible (i, j) 2 pour quantifier les flux. Il est aussi possible d’utiliser alternativement le paramètres VMAI qui est égal à 1/(τ∙C cible (i, j) 2 ), où τ est la durée d’intégration du capteur d’images matriciel 42, et C cible (i, j) le contraste calculé en chaque point d’image i, j, que ce contraste soit spatial, temporel ou spatio-temporel. D’autres équations utilisant le contraste de tavelure pour déduire des informations sur la cible à analyser sont aussi proposées dans la littérature. De manière générale, l’invention proposée implique d’adapter ces différentes équations pour prendre en compte non plus chaque valeur C cible (i, j) comme mesure physique mais ΔC(i, j) = C static (i, j) - cible >, afin de calibrer le résultat. Il est entendu que l’invention peut être reproduite en modifiant des aspects secondaires des modes de mise en œuvre qui ont été décrits en détail ci-dessus, tout en conservant certains au moins des avantages cités. Notamment, la sonde optique peut avoir une constitution et une apparence qui sont différentes de celles qui ont été décrites. Procédé d’imagerie par visualisation de variations d’un contraste de tavelure, suivant lequel des mouvements de parties diffusantes d’une cible à analyser (100) provoquent des variations d’un contraste de tavelure dans une ou plusieurs images de la cible à analyser, le procédé comprenant les étapes suivantes : /1/ fournir une cible de référence (101) qui possède un motif diffusant statique et constant ; /2/ en utilisant un instrument imageur (4) et une source laser (1), saisir au moins deux images de la cible à analyser (100) à deux instants différents, en disposant la cible de référence (101), en plus de la cible à analyser, à l’intérieur d’un champ de vue de l’instrument imageur qui est effectif pour chaque image, de sorte que la source laser éclaire simultanément la cible à analyser et la cible de référence pendant que chaque image est saisie, et que l’image saisie soit formée par du rayonnement produit par la source laser puis diffusé par la cible à analyser et par la cible de référence ; /3/ séparément pour chaque image saisie à l’étape /2/, calculer une valeur de contraste de tavelure, notée C static , à partir d’une partie de l’image qui correspond à la cible de référence (101) ; puis /4/ valider les images saisies à l’étape /2/ uniquement si les valeurs C static calculées respectivement pour lesdites images ont des variations d’une quelconque des images à une autre qui sont inférieures à un seuil, sinon recommencer les étapes /2/ à /4/ ; et ensuite /5/ si les images ont été validées, déduire d’au moins une desdites images des informations sur des mouvements de parties de la cible à analyser (100), à partir de variations du contraste de tavelure relatives auxdites parties de la cible à analyser. Procédé selon la revendication 1, suivant lequel la cible de référence (101) et la source laser (1) sont sélectionnées de sorte que ladite cible de référence possède une valeur d’écart-type de hauteur de rugosité qui est supérieure à un tiers d’une valeur de longueur d’onde de la source laser, de préférence supérieure à une fois ladite valeur de longueur d’onde, et suivant lequel la cible de référence et l’instrument imageur (4) sont sélectionnés en outre de sorte que la partie de chaque image saisie à l’étape /2/ qui correspond à la cible de référence intègre du rayonnement diffusé par au moins dix, de préférence au moins cinquante, reliefs de rugosité séparés par des creux de rugosité intermédiaires. Procédé selon la revendication 1 ou 2, suivant lequel, pour une au moins des images qui ont été validées à l’étape /4/, une valeur de contraste spatial de tavelure est calculée localement dans l’image pour certains au moins des points d’image qui correspondent à la cible à analyser (100), puis une image de variation de contraste de tavelure est construite en attribuant à chacun desdits points d’image une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies conformément à l’étape /2/, et d’autre part la valeur de contraste spatial de tavelure calculée pour le point d’image. Procédé selon la revendication 1 ou 2, suivant lequel une série d’images successives de la cible à analyser(100) est saisie conformément à l’étape /2/ puis, si la série d’images est validée à l’étape /4/, une valeur de contraste temporel de tavelure est calculée séparément pour chacun d’un ensemble de points d’image qui correspondent à la cible à analyser, à partir de valeurs d’intensité respectives des images de la série audit point d’image, puis une image de variation de contraste de tavelure est construite en attribuant à chacun desdits points d’image une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies conformément à l’étape /2/, et d’autre part la valeur de contraste temporel de tavelure calculée pour le point d’image. Procédé selon la revendication 1 ou 2, suivant lequel une série d’images successives de la cible à analyser (100) est saisie conformément à l’étape /2/ puis, si la série d’images est validée à l’étape /4/, une valeur de contraste spatio-temporel de tavelure est calculée séparément pour chacun d’un ensemble de points d’image qui correspondent à la cible à analyser, à partir de valeurs d’intensité respectives des images successives dans un voisinage dudit point d’image, puis une image de variation de contraste de tavelure est construite en attribuant à chacun desdits points d’image une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies conformément à l’étape /2/, et d’autre part la valeur de contraste spatio-temporel de tavelure calculée pour le point d’image. Procédé selon l’une des revendications précédentes, suivant lequel la cible de référence (101) est sélectionnée à l’étape /1/ de sorte que ladite cible de référence possède un niveau d’intensité de diffusion pour le rayonnement de la source laser (1), tel que chaque grain de tavelure dans la partie de chaque image saisie qui correspond à ladite cible de référence, soit à l’intérieur d’un intervalle de détection linéaire de l’instrument imageur (4). Procédé selon l’une des revendications précédentes, suivant lequel l’instrument imageur (4) comprend : - une optique de formation d’image (41) ; - un capteur d’images matriciel (42), comprenant des photodétecteurs disposés à des intersections de lignes et de colonnes ; et - un diaphragme pupillaire (43), et suivant lequel une ouverture du diaphragme pupillaire (43) est ajustée préalablement à la saisie des images à l’étape /2/, de sorte que chaque grain de tavelure recouvre au moins deux photodétecteurs dans la partie de chaque image saisie qui correspond à la cible de référence (101). Procédé selon l’une des revendications précédentes, suivant lequel la cible de référence (101) est sélectionnée à l’étape /1/ pour posséder une taille et des caractéristiques de rugosité telles que la partie de chaque image saisie qui correspond à ladite cible de référence contienne au moins cinquante grains de tavelure. Procédé selon l’une des revendications précédentes, comprenant en outre de caractériser une répartition spectrale d’un rayonnement qui est produit par la source laser (1), telle que ladite répartition spectrale est effective pour au moins une des images saisies, à partir du motif de tavelure dans la partie de ladite image saisie qui correspond à la cible de référence (101). Procédé selon l’une des revendications précédentes, suivant lequel la cible à analyser (100) est une portion d’un fluide qui est partiellement transparent pour le rayonnement de la source laser (1) et qui contient des particules diffusantes mobiles, et le procédé comprend de caractériser une combinaison d’une concentration des particules dans le fluide avec une vitesse de déplacement desdites particules ; ou la cible à analyser (100) est une portion d’un solide diffusant qui vibre, et le procédé comprend de caractériser des amplitudes locales de vibration qui existent en des endroits distincts dans la portion de solide ; ou la cible à analyser (100) est une portion d’un tissu biologique, et le procédé comprend de caractériser des niveaux de mouvements qui existent en des endroits distincts dans le tissu biologique. Système d’imagerie par visualisation de variations d’un contraste de tavelure, comprenant : - un instrument imageur (4), adapté pour saisir des images d’une cible à analyser (100) ; - une source laser (1), disposée pour éclairer la cible à analyser (100) pendant que chaque image est saisie, de sorte que l’image soit formée par du rayonnement produit par la source laser puis diffusé par la cible à analyser ; - une unité de calcul (20), configurée pour calculer des valeurs d’un contraste de tavelure relatives à la cible à analyser (100), à partir d’une ou plusieurs images de ladite cible à analyser ; et - un système d’affichage (30), connecté pour montrer une image qui a été construite par l’unité de calcul (20), caractérisé en ce qu’il comprend en outre : - une cible de référence (101), qui possède un motif diffusant statique et constant ; et - des moyens de support (11), adaptés pour maintenir la cible de référence (101) à l’intérieur d’un champ de vue de l’instrument imageur (4) qui est effectif pour chaque image, et de sorte que la source laser (1) éclaire simultanément la cible à analyser (100) et la cible de référence pendant que chaque image est saisie, et en ce que l’unité de calcul (20) est configurée en outre pour calculer, séparément pour chaque image saisie, une valeur de contraste de tavelure, notée C static , à partir d’une partie de l’image qui correspond à la cible de référence (101), et configurée en outre pour valider plusieurs images saisies successivement si les valeurs C static qui ont été calculées respectivement pour lesdites images ont des variations d’une quelconque des images à une autre qui sont inférieures à un seuil, et pour construire une image de variation de contraste de tavelure en attribuant à chacun d’un ensemble de points d’image qui correspondent à la cible à analyser (100), une valeur de point d’image qui dépend d’un résultat de différence entre, d’une part un résultat de moyenne static > calculé pour les valeurs de contraste de tavelure C static sur toutes les images saisies, et d’autre part une valeur de contraste de tavelure calculée pour le point d’image à partir d’une ou plusieurs images saisies puis validées.