Brevet FR3047830A1 (20170818)
REVENDICATIONS1. Procédé d'estimation de la direction moyenne de déplacements d'objets en mouvement dans un flux d'images numériques {MO d'une scène comprenant du ciel, les images numériques dudit flux étant munies d'un axe de référence Are f, le procédé comprenant : - l'acquisition du flux d'images numériques {MO ; - pour chaque image IN?, du flux d'images numériques {MO o le calcul d'une image segmentée INe, identifiant les objets en mouvement dans l'image INn; o le calcul d'un ensemble de zones d'image {Zka,} respectivement associées à des directions d'analyse {Dk} formant des angles {ak} différents prédéterminés avec l'axe de référence A' f o pour chaque zone Z km, dudit ensemble de zones {Zka,} , le calcul d'une distribution DDk,i, des valeurs des pixels de ladite image segmentée INe, compris dans ladite zone Z kat, la distribution DDka, étant calculée le long d'une direction normale à la direction d'analyse Dk de ladite zone Z km; - pour chaque direction d'analyse Dk, le calcul d'un critère de variabilité CVka, de distributions {DDka,_p, DDka,_p+i, , DDk,72_1, DDka,} calculées pour des zones Z km, associées à ladite direction Dk; et - la détermination de la direction moyenne Dkoptal de déplacement des objets dans le flux d'images numériques {/Nn} comme étant égale à la direction d'analyse pour laquelle le critère de variabilité CV/c,' est optimum.\n<CLM>\1\tProcédé selon la revendication 1, caractérisé en ce que l'acquisition du flux d'image numériques {/Nn} comprend : - l'acquisition de la scène par une dispositif de prise de vue ayant un angle de champ supérieur ou égal à 100° selon la diagonale d'un capteur au format 35mm ou en équivalent au format 35mm, en particulier un angle de champ égal à 180° ; - et préalablement au calcul des distributions DDk,', le redressement des images les images numériques {/Nn} .\n<CLM>\1\tProcédé selon la revendication 1 ou 2, caractérisé en ce que la scène est acquise à l'aide d'un dispositif d'imagerie comprenant un système optique apte à former une image de la scène sur un plan focal et un capteur numérique plan disposé sur le plan focal du système optique pour produire le flux d'images numérique {/Nn} , et en 3047830 34 5\n<CLM>\1\t10\n<CLM>\1\t15\n<CLM>\1\t20\n<CLM>\1\t25 30 ce que le procédé comprend le calcul de l'ensemble de zones d'image {Zka,} pour chaque image I A In du flux d'images numériques {MO , ledit calcul comprenant : - le calcul d'un ensemble de rectangles {Rk,'} ayant chacun un petit côté centré sur un même point prédéterminé du plan focal, et ayant leurs grands côtés respectivement parallèles aux directions d'analyse {Dk} ; et - le calcul des zones {Zka,} comme étant égales respectivement à l'intersection des rectangles {Rka,} avec l'image I Ni,. Procédé selon la revendication 3, caractérisé en ce que la scène comprend du ciel, et en ce que le point prédéterminé du plan focal est la position du centre du disque solaire sur le plan focal du système optique. Procédé selon la revendication 4, caractérisé en ce le petit côté de chaque rectangle Rical de l'ensemble de rectangles {Rk,72} est au moins égal au diamètre du halo Di,' du disque solaire. Procédé selon la revendication 4 ou 5, caractérisé en ce que la position du centre du disque solaire est calculée en fonction d'un modèle astronomique de l'azimut et de l'élévation du soleil dans le ciel définissant la position dudit centre sur le plan focal. Procédé selon la revendication 6, caractérisé en ce que le modèle astronomique est réglé pour une position prédéterminée (Of , 0 , Azr) du dispositif d'imagerie, définie par la position (Of ,\n<CLM>\1\tdu capteur numérique par rapport à un plan horizontal et par la position d'un axe du capteur par rapport à un azimut (Azr), et en ce que le procédé comprend une étape de détermination de la position prédéterminée (Of , 0 , Azr) du dispositif d'imagerie, comportant : - l'acquisition d'un flux d'images numériques de calibration horodatées {/N79 comprenant chacune le disque solaire ; - pour chaque image In dudit flux {/N79 , le calcul d'une position Prf du centre du disque solaire dans l'image I n en analysant l'image In ; - le calcul de la position du dispositif d'imagerie (Of , 0 , Azr) comme solution d'un problème d'optimisation selon a relation : /Pi 1 (Pf (Bi, 02, AZ)) (0f, 0 , Azr) = ar gmin (01,02,Az)612 vsnc:Pn (Bi, 02, AZ) 35 où D. est un ensemble de recherche, et II II est une fonction de coût quantifiant un écart entre deux vecteurs de 3 3047830 35\n<CLM>\1\tProcédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le critère de de variabilité CVka, pour la direction Dk est calculé en fonction des distributions IDDk,n_p, DDk,n_p+i, , DDk,n1 associées aux p > 2 dernières images numériques acquises.\n<CLM>\1\tProcédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le critère de variabilité CVk,' pour la direction Dk est calculé selon la relation : CVkal = Dpk,j11 j=n-p+1 5 10 15\n<CLM>\1\t20\n<CLM>\1\t25 30\n<CLM>\1\texpression dans laquelle j référence la position d'une image numérique dans le flux d'images numériques {/Nn} , n et n-p référence respectivement la nième image et la (n-p)lème image dans le flux d'images numériques {/Nn} , avec p strictement inférieur à n, k référençant la klème direction Dk, et 11 V II est une norme du vecteur V de l'espace l'Zi , I étant la dimension des distributions DDk,n. Procédé selon l'une quelconque des revendications, caractérisé en ce qu'il comporte le calcul d'un masque d'exclusion des éléments immobiles de la scène, et en ce que préalablement au calcul des distributions DDk,n, ledit masque est appliqué aux images du flux d'images numériques {/Nn} et/ou aux images segmentées {Min . Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre : - le calcul d'une vitesse moyenne Vn des objets en mouvement dans le flux d'images numériques {/Nn} selon ladite direction moyenne de déplacement calculée ; - et la prédiction de la position d'un objet en mouvement en fonction de ladite vitesse moyenne Vn et de ladite direction moyenne Dkopt,n. Procédé selon la revendication 11, caractérisé en ce qu'il comprend la prévision de l'ensoleillement au-dessus d'un parc photovoltaïque, en particulier la chute d'ensoleillement. 3047830 36\n<CLM>\1\tProcédé selon la revendication 11 ou 12, caractérisé en ce que le calcul de ladite vitesse moyenne Vn comprend : pour chaque image segmentée INn, le calcul d'une distribution DVkopt,n des valeurs des pixels de l'image segmentée INe, appartenant à la zone Zkoptal 5 associée à la direction moyenne déterminée Dkopt,n, la distribution DVkopt,n étant calculée le long de la direction moyenne Dkoptal déterminée ; le calcul de l'inclinaison d'au moins une structure affine comprise dans une matrice DVn = [DVkopt,n_qDVkopt,n_q+1 DVkopt,n_1DVkopt,n1 issue de la concaténation des distributions DVkopt,n des images entre la (n-q)lème image du 10 flux jusqu'à la nième du flux d'images numériques {INn} , où q est un entier supérieur ou égal à 1; et le calcul de la vitesse moyenne Vn en fonction de l'inclinaison calculée.\n<CLM>\1\tProcédé selon la revendication 13, caractérisé en ce que le calcul de l'inclinaison 15 de la au moins une structure affine comprise dans une matrice DV, comprend : l'application de différentes rotations d'angles leu prédéterminés à la matrice DV, ; le calcul d'un histogramme H istnx pour chaque matrice DV,,' obtenue par rotation et le calcul du maximum Mn,' de cet histogramme ; 20 la détermination d'inclinaison comme étant égal à l'angle &opta,. correspondant au maximum des maxima Mn,' calculés.\n<CLM>\1\tSystème d'estimation de la direction moyenne de déplacements d'objets en mouvement dans un flux d'images numériques {INn} d'une scène comprenant du 25 ciel, les images numériques dudit flux étant munies d'un axe de référence Are f, le système comprenant : - une caméra pour l'acquisition du flux d'images numériques {INn} de ladite scène; une unité de traitement d'information configurée pour : 30 o pour chaque image INn du flux d'images numériques {INn} ^ le calcul d'une image segmentée INe, identifiant les objets en mouvement dans l'image INn; ^ le calcul d'un ensemble de zones d'image {Zka2} respectivement associées à des directions d'analyse {Dk} formant des angles {ak} 35 différents prédéterminés avec l'axe de référence A' f 3047830 37 ^ pour chaque zone Z kat dudit ensemble de zones {Zka,} , le calcul d'une distribution DDka, des valeurs des pixels de ladite image segmentée INn compris dans ladite zone Z kat, la distribution DDka, étant calculée le long d'une direction normale à la direction d'analyse 5 Dk de ladite zone Z km; o pour chaque direction d'analyse Dk, le calcul d'un critère de variabilité CVk,i, de distributions {DDka,_p, DDka,_p+i, , DDk,72_1, DDka,} calculées pour des zones Z kat associées à ladite direction Dk; et o la détermination de la direction moyenne Dkoptal de déplacement des 10 objets dans le flux d'images numériques {/Nn} comme étant égale à la direction d'analyse pour laquelle le critère de variabilité CVka, est optimum.\n<CLM>\1\tSystème selon la revendication 15 configuré pour mettre en oeuvre un procédé selon 15 l'une quelconque des revendications 2 à 14.