Modèles de distorsion

Les différents logiciels de photogrammétrie implémentent des modèles mathématiques variés pour corriger la distorsion. Un étalonnage réalisé avec un logiciel ne peut pas facilement être utilisé dans un autre logiciel. Si les modèles reposent sur des principes différents (corrections par grilles, non radiales, etc), il sera nécessaire de prévoir un module de transformation spécifique, à condition de disposer des documentations adéquates.

Mais, dans certains cas, la transformation entre modèles suffisamment proches peut être simple, par exemple entre Redresseur et Photoscan, qui implémentent tous deux un modèle radial polynômial. Le modèle de Photoscan est un peu moins général, il ne sépare pas le point principal du centre de symétrie de la distorsion. .Par contre, Photoscan prévoit des pixels rectangulaires, ce qui amène la présence de deux distances principales en X et Y (mais je les ai toujours vues identiques). De plus, l’axe Y image est orienté différemment. Et les coefficients du polynôme ne se réfèrent pas à la même unité, pixel pour Redresseur distance principale pour Photoscan.

Redresseur Photoscan
xc Cx = xc
yc Cy = H-yc (H : nb de lignes image)
xpps absent
ypps absent
p Fx,fy = p
a K1 = a.p^2
b K2 = b.p^4
c K3 = c.p^6

A titre de test, j’ai calculé le même étalonnage (sur le chevet de Salagon) avec Redresseur et Photoscan. Voici les résultats :

Redresseur Photoscan Redresseur transformé
Xc =2391.6 Cx = 2408.8 2391.6
Yc=1606.8 Cy = 1651.2 1657.2
Xpps=2364.3  —-
Ypps=1607.0  —-
P=3770.8 Fx,fy = 3769.9 3770.8
a=-6.340229e-9 K1 = -8.93962e-2 -9.01513e-2
b=8.605516e-17 K2 = 1.70934e-2 1.73998e-2
c=2.308141e-24 K3 = 6.84810e-3 6.63534e-3

Comme souvent, les valeurs de centrages sont un peu différentes (ce sont les plus sensibles à la répartition des points). Focale et distorsion sont très proches.

Pour passer le modèle de distorsion de Redresseur à Apero-Micmac, rien de plus simple: Les modèles RadialExtended et Fraser (avec LibAff et LibDec à false) incluent le modèle de Redresseur (en mettant à 0 les coefficients supplémentaires)

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *