Procédé de génération d’un modèle de représentation de relief à partir d’une pluralité de N paires d’images stéréo, chacune des N paires d’images stéréo étant associée à une carte de disparités photogrammétriques qui est obtenue par calcul de corrélation à partir de la paire d’images stéréo en question, le procédé étant caractérisé en ce qu’il comporte les étapes itératives suivantes : - calcul de cartes de disparités géométriques à partir d’un modèle de représentation de relief prédéterminé, par projection du modèle de représentation de relief prédéterminé dans une géométrie épipolaire des paires d’images stéréo ; - calcul d’une fonction de coût, représentative d’une différence entre les disparités géométriques et photogrammétriques, associée au modèle de représentation de relief prédéterminé ; et - mise à jour du modèle de représentation de relief prédéterminé par optimisation de type descente de gradient de la fonction de coût, et réitération jusqu’à atteinte d’un critère d’arrêt prédéfini. Procédé selon la revendication 1, dans lequel la fonction de coût FC est exprimée ainsi : où k est un indice de parcours des N paires d’images stéréo, N m est une fonction de normalisation, D k (i,j) est la valeur de disparité pour le pixel (i,j) dans la carte de disparités géométriques associée à la paire d’images stéréo pointée par la valeur de l’indice k , et D k,0 (i,j) est la valeur de disparité pour le pixel (i,j) dans la carte de disparités photogrammétriques associée à la paire d’images stéréo pointée par la valeur de l’indice k , et W k,0 (i,j) est la valeur correspondante pour le pixel (i,j) dans une carte de scores de corrélation associée, à l’issue du calcul de corrélation, à ladite carte de disparités photogrammétriques associée à la paire d’images stéréo pointée par la valeur de l’indice k . Procédé selon la revendication 2, dans lequel la fonction de normalisation N m est la pseudo-norme de Huber ainsi appliquée : Procédé selon la revendication 2 ou 3, dans lequel la fonction de coût FC est ajustée avec un terme R de régularisation du modèle de représentation de relief, de telle sorte que : où est un facteur de pondération de régularisation du modèle de représentation de relief et FCA est la fonction de coût ajustée. Procédé selon la revendication 4, dans lequel le terme R de régularisation du modèle de représentation de relief est la variation totale du modèle de représentation de relief, ou sa variante Huber-TV, laquelle est normalisée par un pas de grille du modèle de représentation de relief. Procédé selon l’une quelconque des revendications 1 à 5, dans lequel la projection du modèle de représentation de relief prédéterminé dans la géométrie épipolaire des paires d’images stéréo, appelée ici géométrie épipolaire initiale, est effectuée dans une géométrie épipolaire sur-échantillonnée, et le résultat de la projection est ensuite sous-échantillonné dans la géométrie épipolaire initiale. Procédé selon la revendication 6, dans lequel le résultat de la projection est sous-échantillonné dans la géométrie épipolaire initiale à l’aide d’un noyau de convolution bilinéaire réalisant un filtrage passe-bas. Procédé selon l’une quelconque des revendications 1 à 7, dans lequel le modèle de représentation de relief est un modèle numérique d’élévation. Produit programme d’ordinateur comportant des instructions pour implémenter le procédé selon l’une quelconque des revendications 1 à 8, lorsque lesdites instructions sont exécutées par un processeur. Support de stockage d’informations stockant un programme d’ordinateur comportant des instructions pour implémenter le procédé selon l’une quelconque des revendications 1 à 8, lorsque lesdites instructions sont lues depuis le support de stockage d’informations et exécutées par un processeur. Système de génération d’un modèle de représentation de relief à partir d’une pluralité de N paires d’images stéréo, chacune des N paires d’images stéréo étant associée à une carte de disparités photogrammétriques qui est obtenue par calcul de corrélation à partir de la paire d’images stéréo en question, le système comportant de la circuiterie électronique configurée pour implémenter les étapes itératives suivantes : - calcul de cartes de disparités géométriques à partir d’un modèle de représentation de relief prédéterminé, par projection du modèle de représentation de relief prédéterminé dans une géométrie épipolaire des paires d’images stéréo ; - calcul d’une fonction de coût, représentative d’une différence entre les disparités géométriques et photogrammétriques, associée au modèle de représentation de relief prédéterminé ; et - mise à jour du modèle de représentation de relief prédéterminé par optimisation de type descente de gradient de la fonction de coût, et réitération jusqu’à atteinte d’un critère d’arrêt prédéfini.