FR 3030847 A1 20160624 FR 1402929 A 20141219 PROCEDE DE DISCRIMINATION ET D'IDENTIFICATION PAR IMAGERIE 3D D'OBJETS D'UNE SCENE Le domaine de l'invention est celui de la discrimination et de l'identification d'objets quelconques qui peuvent être dissimulés dans une scène complexe comme illustré figure 4. Ces objets peuvent se masquer les uns les autres, ce qui rend l'identification et la discrimination des objets particulièrement difficiles. Différents types de visualisation tridimensionnelle existent. Certaines visualisations offrent une image en profondeur, ou image 2,5 D. Ces techniques ne permettent pas de récupérer l'ensemble des voxels (« volumetric pixels ») tridimensionnels, elles sont fondées sur les techniques stéréoscopiques liées à l'optimisation d'une différence de marche du rayon optique ou électromagnétique entre les différents points de deux images ; cette différence de marche est calculable si les objets se situent à courte portée du système d'imagerie. L'identification est donc limitée aux courtes distances. Certaines techniques d'imagerie nécessitent un balayage du faisceau laser ou de lumière sur la scène, et une séquence de mesures donne un nuage de points 3D. Plusieurs balayages à plusieurs angles de vue sont nécessaires pour obtenir une information tridimensionnelle mais le mouvement relatif du porteur contenant le système de balayage produit une distorsion de nuage de points à nuage de points, le recalage est complexe et dégrade la résolution et l'identification. Les techniques dites des silhouettes obtenues à partir de vues multiples donnent une enveloppe extérieure peu détaillée. Si les objets contiennent une zone d'ombre ou si un objet se trouve dans l'ombre d'un autre objet, une grande partie de l'enveloppe est perdue donc l'identification et la discrimination des deux objets n'est pas réalisable. Les technologies type Spin Image nécessitent des bases de 30 données et ne peuvent pas être appliquées pour des objets dissimulés dans des scènes complexes sans connaissance a priori. Les techniques de regroupement de voxels lié à une reconstruction de surface par isodensité permettent d'obtenir la surface externe des objets mais éliminent les informations internes contenues dans l'objet. Les techniques de reconstruction tridimensionnelle liées à une optimisation de bases de connaissances et à l'extraction de signaux faibles et dont un exemple est présenté dans le document « Procédé de reconstruction synthétique d'objets illuminés par une onde électromagnétique et/ou excités par une onde élastique » (EP 2 929 421 ou US 8 345 960), nécessitent des bases de données et que celles-ci soient optimisées afin de réaliser la meilleure reconstruction tridimensionnelle possible. Les bases de connaissance ne représentent souvent les objets que par un ensemble de surfaces tridimensionnelles externes qui peuvent être modélisées par des facettes. Les structures transparentes (vitres...) sont mal prises en compte, ce qui nuit fortement à la visualisation complète de l'objet en 3D. Ces techniques sont aussi tributaires de la disponibilité des bases de connaissance. La technique MIP (Maximum lntensity Projection) permet la visualisation 2D de données tridimensionnelles. Cette technique projette les voxels sur un plan de projection ; les voxels sont déterminés par les rayons joignant le plan de projection au point d'observation avec un seuil d'intensité imposé. Les résultats 2D ne permettent pas d'obtenir une valeur de la profondeur et des distances. Pour créer une illusion de rotation et une notion de profondeur et ainsi améliorer le rendu 3D, plusieurs plans de projection sont réalisés à des angles d'observation successifs. Les techniques de rendu de l'intensité des voxels permettent une visualisation tridimensionnelle bruitée, ce qui diminue la discrimination entre différents contigus de la scène observée. Les nuages de points 3D obtenus directement par un positionnement cartésien des voxels ne permettent qu'une faible discrimination et apportent des artefacts liés aux fausses alarmes. Les techniques de complétion surfacique et dont un exemple est présenté dans le document « Procédé de reconstruction 3d d'un objet d'une scène » (US2013/0100131), permettent de remplir les zones lacunaires par génération de surface tridimensionnelle pour obtenir un ensemble de données complétées de l'objet, sans recourir à une base de données extérieure. Elles donnent ainsi une surface externe complète de l'objet sans discontinuité surfacique ; les données internes situées derrière des transparents (vitres par exemple) ne sont pas reconstruites et se trouvent éliminées du processus d'identification. Un processus de discrimination de deux objets proches est plus difficilement atteignable, parce que la génération de surfaces peut être mathématiquement complexe. Le but de l'invention est de pallier ces inconvénients. Plus précisément l'invention a pour objet un procédé de discrimination et d'identification par imagerie 3D d'un objet dans une scène complexe. Il est principalement caractérisé en ce qu'il comporte les étapes suivantes : - générer une séquence d'images dites images 2D MIP de l'objet, à partir d'un volume 3D en voxels de la scène complexe, ce volume étant prédéterminé et visualisé par un opérateur en utilisant un processus itératif de type MIP à partir d'un plan de projection et d'un seuil d'intensité déterminés par l'opérateur à chaque itération, extraire automatiquement à partir de la séquence d'images 2D MIP, des coordonnées d'un volume réduit correspondant à la séquence d'images 2D MIP, - choisir un des seuils d'intensité utilisés au cours des itérations du processus de type MIP, - extraire automatiquement à partir du volume 3D en voxels de la scène complexe, des coordonnées et du seuil d'intensité choisi, un volume 3D en voxels réduit contenant l'objet, - générer automatiquement à partir du volume réduit, par optimisation en seuillage d'intensité, un seuil d'intensité optimisé et un volume en voxels optimisé, une couleur étant associée à chaque intensité, identifier l'objet par visualisation. Il comprend éventuellement une étape de visualisation du volume 3D en voxels réduit et/ou du volume en voxels optimisé.35 Selon une caractéristique de l'invention, il comprend en outre une étape pour générer à partir de la séquence d'images 2D MIP, des coordonnées et du seuil d'intensité choisi, un nuage de points 3D (8) de l'objet. Selon une autre caractéristique de l'invention, il comprend il comporte en outre les étapes suivantes : générer à partir du volume 3D de la scène complexe et du seuil d'intensité choisi, un nuage de points 3D brut de la scène complexe, générer à partir du volume 3D optimisé et du seuil d'intensité optimisé, un nuage de points 3D optimisé de l'objet, générer à partir d'une superposition du nuage de points 3D brut de la scène, du nuage de points 3D optimisé de l'objet, et éventuellement du nuage de points 3D de l'objet, un nuage de points 3D global optimisé de l'objet inclus dans la scène complexe, visualiser le nuage 3D global optimisé. Ces nuages peuvent également être visualisés par l'opérateur. Le procédé selon l'invention permet une discrimination, une identification et éventuellement un positionnement de l'objet. La précision des détails de l'objet est rendue possible par la densification du nuage de points, ce qui permet aussi d'améliorer l'identification de l'objet de la scène. Les nuages de points tridimensionnels permettent de séparer et de discriminer les différents objets de la scène et de discriminer l'objet à identifier des autres éléments de la scène. L'ensemble des nuages de points associés à l'objet est positionné dans un repère cartésien absolu, ce qui permet de calculer simplement les positions de l'objet et les différentes distances entre l'objet et d'autres éléments de la scène. L'objet est replacé dans son contexte opérationnel. La solution apportée comble les inconvénients individuels de chaque technique de visualisation existante. Il n'est pas nécessaire de compléter les données (notamment les données masquées par d'autres éléments de la scène) car elles sont implicitement présentes.35 D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels : la figure 1 représente schématiquement un exemple de volume 3D 5 avec le quadrillage correspondant aux différents voxels, obtenu par reconstruction 3D, la figure 2 montre un exemple d'enchaînement des différentes étapes du procédé selon l'invention, la figure 3 montre un exemple de visualisation 3D d'un volume 3D 10 en voxels d'une scène complexe, sans le quadrillage et obtenu par reconstruction 3D, la figure 4 montre un exemple d'image 2D d'une scène complexe, la figure 5 montre un exemple de séquence d'images 2D MIP d'un objet à identifier, 15 la figure 6 montre un exemple de visualisation 3D du volume 3D réduit de l'objet à identifier, la figure 7 montre un exemple de visualisation 3D du volume 3D optimisé de l'objet à identifier, la figure 8 montre un exemple de visualisation 3D d'un nuage de 20 points 3D MIP de l'objet à identifier, la figure 9 montre un exemple de visualisation 3D d'un nuage de points 3D brut de la scène complexe, la figure 10 montre un exemple de visualisation 3D d'un nuage de points 3D optimisé de l'objet à identifier, 25 les figures 11 montre un exemple de visualisation 3D d'un nuage de points 3D global optimisé de l'objet dans la scène complexe, avec comme objet une voiture (fig 11a), ou deux voitures (fig 11b). D'une figure à l'autre, les mêmes éléments sont repérés par les 30 mêmes références. Le procédé selon l'invention est appliqué à un volume 3D en voxels d'une scène complexe. Cette scène a donc été préalablement transformée par imagerie 3D en ce volume 3D (1) dont un exemple est 35 montré figure I. Ce volume tridimensionnel peut être obtenu grâce à un procédé de reconstruction par transmission ou par fluorescence (Optical Projection Tomography, imagerie nucléaire ou X-Ray Computed Tomography) ou par réflexion (réflexion d'une onde laser ou par réflexion solaire dans le cas de la bande visible (entre 0,4pm et 0,7pm) ou proche infrarouge (entre 0,7pm et 1pm ou SWIR entre 1pm et 3pm) ou en prenant en compte l'émission thermique de l'objet (imagerie thermique entre 3pm et 5pm et entre 8pm et 12pm), ce processus de reconstruction tridimensionnelle est décrit dans le brevet « Système optronique et procédé d'élaboration d'images en trois dimensions dédiés à l'identification » (US 8,836,762, EP 2 333 481). On utilise l'ensemble des voxels issus d'une reconstruction tridimensionnelle avec l'intensité associée, cette reconstruction ayant de préférence été obtenue par réflexion. Le procédé de discrimination et d'identification selon l'invention comporte les étapes suivantes décrites en relation avec la figure 2. On désigne par objet un ou plusieurs éléments de la scène : un objet peut être une voiture comme montré dans les exemples des figures 2 à 11a, ou deux voitures comme montré dans l'exemple de la figure 11b, ou encore un élément de voiture comme la jante. Parmi les étapes suivantes, certaines requièrent la participation d'un opérateur à partir d'une visualisation sur écran, d'autres sont mises en oeuvre automatiquement. - A) à partir du volume 3D en voxels (1) de la scène complexe (représenté sur les figures 2 et 3 sans le quadrillage des voxels), la technique MIP est utilisée pour que l'opérateur visualise en temps réel sur son écran par exemple la scène complexe suivant différents plans de projection (correspondant à différents angles d'observation) et selon un seuil d'intensité déterminés expérimentalement par l'opérateur, et pour qu'au cours de ces différentes visualisations dont un exemple est montré figure 4, il dégrossisse le volume à analyser incluant un objet qu'il souhaite discriminer et/ou identifier. A l'issue de ce processus itératif de type MIP, une séquence d'images 2D MIP (2) de l'objet à identifier, a été générée, qui correspond visuellement à un volume 3D réduit contenant l'objet à identifier, mais ce volume réduit n'est pas encore déterminé à ce stade. Un exemple d'une telle séquence de trois images 2D MIP (2) est montré figures 2 et 5. - B1) extraire automatiquement à partir de cette séquence 5 d'images 2D MIP (2), des coordonnées (3) de ce volume réduit contenant l'objet à identifier. - B2) déterminer un des seuils en intensité (4) qui a été utilisé dans le processus MIP de l'étape A) ; ce seuil est déterminé 10 expérimentalement par l'opérateur - C) extraire automatiquement à partir du volume 3D en voxels (1) de la scène complexe, et en utilisant les coordonnées (3) et le seuil en intensité (4), un volume en voxels (5) correspondant à l'objet à identifier : il 15 s'agit du volume réduit de l'étape A). L'opérateur peut visualiser ce volume réduit (5) s'il le souhaite, comme montré figures 2 et 6 ; cette visualisation peut lui offrir une possibilité d'identifier l'objet. - D) générer automatiquement à partir du volume en voxels (5), 20 par un processus d'optimisation en seuillage d'intensité, chaque intensité étant associée à une couleur : - un seuil optimisé (7), - puis en utilisant ce seuil optimisé, un volume 3D optimisé en voxels (6) de l'objet à identifier. 25 On obtient ainsi en quelque sorte le volume réduit, en couleur après seuillage, que l'opérateur peut également visualiser comme montré figures 2 et 7, et identifier. On obtient un volume 3D en voxels qui présente visuellement beaucoup plus de détails que le volume 3D d'origine. Au cours des étapes précédentes, A à C, l'opérateur a eu la 30 possibilité d'identifier l'objet avec une probabilité croissante. A l'issue de cette étape D, l'opérateur a atteint la plus forte probabilité d'identifier l'objet. L'objet à identifier peut être lui-même divisé en sous-objets. Les sous-objets d'une voiture sont par exemple les roues, le capot, ... . Les 35 étapes précédentes sont alors appliquées à chaque sous-objet. On obtient pour chaque sous-objet, une séquence d'images 2D MIP (2), des coordonnées (3) de volume 3D réduit, un seuil d'intensité choisi (4), un volume 3D en voxels réduit (5), un seuil optimisé (7) et un volume 3D optimisé (6). A l'issue de cette étape D, lorsque l'objet comporte plusieurs éléments ou sous-objets présentant un intérêt pour l'opérateur, ceux-ci ont pu être discriminés. Les étapes suivantes visent à remettre en situation l'objet identifié. - E) éventuellement générer automatiquement à partir de la séquence d'images 2D MIP (2) en utilisant les coordonnées (3) et le seuil en intensité (4), un nuage de points 3D (8) de l'objet à identifier, que l'opérateur 15 peut éventuellement visualiser, comme montré figures 2 et 8. - F1) générer automatiquement à partir de la séquence d'images 2D MIP (2) et en utilisant le seuil en intensité (4), un nuage de points 3D brut (9) de la scène complexe incluant donc l'objet à identifier, que l'opérateur 20 peut éventuellement visualiser, comme montré figures 2 et 9. - F2) générer automatiquement à partir du volume 3D optimisé en voxels (6) et en utilisant le seuil en intensité (7), un nuage de points 3D optimisé (10) de l'objet à identifier, que l'opérateur peut éventuellement 25 visualiser, comme montré figures 2 et 10. - G) générer automatiquement par superposition : - du nuage de points 3D brut (9) de la scène complexe, 30 du nuage de points 3D optimisé (10) de l'objet, qui présente une forte densité au niveau des détails de l'objet permettant son identification, et du nuage de points 3D (8) de l'objet s'il a été généré, 35 un nuage de points 3D (11) global optimisé de l'objet inclus dans la scène complexe. Ce nuage permet de positionner l'objet dans son contexte opérationnel en prenant en compte les coordonnées absolues, de le discriminer, d'aider l'opérateur à l'identifier avec une forte probabilité. En effet, l'opérateur visualise ce nuage (11) global optimisé, dont un exemple est montré figures 2 et 11 a pour un objet qui représente une voiture ou figure llb pour un objet qui représente deux voitures ; cela lui permet de visualiser l'objet identifié dans son contexte opérationnel. Les différentes extractions des nuages de points permettent une diminution des artéfacts et du bruit. L'ensemble des nuages de points associés à l'objet est positionné dans un repère cartésien absolu, ce qui permet de calculer simplement la position de l'objet et les différentes distances entre l'objet et d'autres éléments de la scène. L'objet est replacé dans son contexte opérationnel. Lorsque des sous-objets ont été identifiés, chacune de ces étapes de mise en situation est appliquée en superposant les nuages respectivement obtenus à partir de chaque sous-objet. Lorsque l'opérateur considère que le résultat n'est pas satisfaisant, ces étapes A à G peuvent être itérées afin d'obtenir la meilleure densification possible de ces détails. A l'issue de l'étape G, lorsque l'objet comporte plusieurs éléments ou sous-objets discriminés entre eux, ceux-ci sont replacés dans leur contexte et leur positionnement, les uns par rapport aux autres, peut être mesuré dans des coordonnées cartésiennes. Avec ce procédé, il est particulièrement aisé de visualiser en temps réel les différents éléments de la scène observée. Selon l'invention, l'identification des objets de la scène est améliorée grâce à une densification du nuage de points tridimensionnel. Ce procédé de discrimination et d'identification s'applique notamment au domaine de la sécurité pour l'identification et la localisation précise de cibles, et au domaine médical pour l'identification et la localisation de tumeurs. La présente invention peut s'implémenter à partir d'éléments 5 matériel et/ou logiciel. Elle peut être disponible en tant que produit programme d'ordinateur sur un support lisible par ordinateur, ce programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé de discrimination et d'identification. Le support peut être électronique, magnétique, optique, électromagnétique ou être un support de 10 diffusion de type infrarouge. De tels supports sont par exemple, des mémoires à semi-conducteur (Random Access Memory RAM, Read-Only Memory ROM), des bandes, des disquettes ou disques magnétiques ou optiques (Compact Disk - Read Only Memory (CD-ROM), Compact Disk - ReadANrite (CD-R/W) and DVD). 15 Bien que l'invention ait été décrite en liaison avec des modes de réalisation particuliers, il est bien évident qu'elle n'y est nullement limitée et qu'elle comprend tous les équivalents techniques des moyens décrits ainsi que leurs combinaisons si celles-ci entrent dans le cadre de l'invention. 20 REVENDICATIONS1. Procédé de discrimination et d'identification par imagerie 3D d'un objet dans une scène complexe, caractérisé en ce qu'il comporte les étapes suivantes : générer (A) une séquence d'images dites images 2D MIP de l'objet (2), à partir d'un volume 3D en voxels de la scène complexe, ce .10 volume étant prédéterminé (1) et visualisé par un opérateur en utilisant un processus itératif de type MIP à partir d'un plan de projection et d'un seuil d'intensité déterminés par l'opérateur à chaque itération, - extraire (B1) automatiquement à partir de la séquence d'images 15 2D MIP (2), des coordonnées (3) d'un volume réduit correspondant à la séquence d'images 2D MIP (2), choisir (B2) un (4) des seuils d'intensité utilisés au cours des itérations de l'étape A), - extraire (C) automatiquement à partir du volume 3D en voxels de 20 la scène complexe (1), des coordonnées (3) et du seuil d'intensité choisi (4), un volume 3D en voxels réduit (5) contenant l'objet, générer (D) automatiquement à partir du volume réduit (5), par optimisation en seuillage d'intensité, un seuil d'intensité optimisé (7) et un volume en voxels optimisé (6), une couleur étant 25 associée à chaque intensité, - identifier l'objet par visualisation. 2. Procédé de discrimination et d'identification d'un objet dans une scène complexe selon la revendication précédente, caractérisé en ce qu'il 30 comprend une étape de visualisation du volume 3D en voxels réduit (5). 3. Procédé de discrimination et d'identification d'un objet dans une scène complexe selon l'une des revendications précédentes, caractérisé en 35 ce qu'il comprend une étape de visualisation du volume en voxels optimisé (6). 4. Procédé de discrimination et d'identification d'un objet dans une scène complexe selon l'une des revendications précédentes, caractérisé en ce qu'il comprend en outre une étape pour générer (E) à partir de la séquence d'images 2D MIP (2), des coordonnées (3) et du seuil d'intensité choisi (4), un nuage de points 3D (8) de l'objet. 5. Procédé de discrimination et d'identification d'un objet dans une scène complexe selon l'une des revendications précédentes, caractérisé en ce qu'il comporte en outre les étapes suivantes : - générer (F1) à partir du volume 3D de la scène complexe (1) et du seuil d'intensité choisi (4), un nuage de points 3D brut (9) de la scène complexe, générer (F2) à partir du volume 3D optimisé (6) et du seuil d'intensité optimisé (7), un nuage de points 3D optimisé (10) de l'objet, générer (G) à partir d'une superposition du nuage de points 3D brut (9) de la scène et du nuage de points 3D optimisé (10) de l'objet, un nuage de points 3D global optimisé (11) de l'objet inclus dans la scène complexe, visualiser le nuage 3D global optimisé. 6. Procédé de discrimination et d'identification d'un objet dans une scène complexe selon les revendications 4 et 5 prises en combinaison, caractérisé en ce que le nuage de points global optimisé (11) est généré à partir de la superposition du nuage de points 3D brut (9) de la scène, du nuage de points 3D optimisé (10) de l'objet, et du nuage de points 3D (8) de l'objet. 7. Procédé selon la revendication 4 ou 5, caractérisé en ce que le nuage de points 3D (8) de l'objet et/ou le nuage de points 3D (9) de la scène et/ou le nuage de points 3D (10) optimisé de l'objet est visualisé par l'opérateur.35 8. Procédé selon l'une des revendications précédentes, caractérisé en ce que le volume 3D en voxels (1) est obtenu par réflexion d'une onde optique visible ou IR sur la scène complexe ou par émission thermique de l'objet. 9. Un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé selon l'une quelconque des revendications 1 à 8, lorsque ledit programme est exécuté sur un ordinateur. 10