[Investigación y cálculo de elementos de orientación interna y externa de fotografías basado en programación VC++6.0] Elementos de orientación externa de fotografías
Clasificación de bibliotecas chinas. :P23 Código de identificación del documento: A Número de documento: 1672-3791(2012)05(b)-0033-01 Durante el proceso de fotogrametría aérea, es necesario realizar un buen trabajo en levantamientos de campo y fotografía aérea.
El principio de cálculo de 1
Utilizando un cierto número de puntos de control terrestre y las coordenadas de los puntos correspondientes en el mapa, se pueden obtener los elementos direccionales.
La idea básica es: a partir de una única imagen y según la ecuación colineal, resolver los elementos de acimut de la imagen durante la fotografía aérea. Dado que la ecuación colineal del modelo matemático utilizado para la resección espacial es una función no lineal, la linealizamos. La ecuación colineal es:
Después de la linealización, quedan:
x =(x)+a 11 dxs+a 12 Dy+a 13 dzs+a 14dφ+a 15dω+a 16dκ+a 17x+a 18y+a 19f
y =(y)+a 21 dxs+a 22 dis+a 23 dzs+a24dφ+a25dω+a26dκ+a27x+a28y+a29f
Cuando utilizamos cuatro o más puntos de control terrestre y las coordenadas de los puntos de imagen correspondientes, utilizamos el método de ajuste de mínimos cuadrados para resolver el problema, más los errores accidentales Vx, Vy,
VX = a 11 dxs +a 12 dis+a 13 dzs+a 14dφ+a 15dω+a 16dκ+a 17dx+a 18dy+a 19df-Lx
vy = a 21 dxs+a 22 dis +a 23 dzs+a24dφ +a25dω+a26dκ+a27dx+a28dy+a29df-Ly
El coeficiente de cada término se puede obtener resolviendo la derivada parcial, donde:
a 17 = (x-x0) /f; a27 = (y-y0)/f; a 18 = 1; a28 = 0; a 19 = 0;
Entonces podemos obtener la ecuación de error:
V=AX-L
Las ecuaciones normales se enumeran de la siguiente manera: (A-T? PA)X=A-T? PL; La solución de la ecuación normal es: X=(A-TA)-1ATL.
Resolviendo así los nueve elementos azimutales.
2 Diagrama de flujo del algoritmo
Los resultados del cálculo de la pieza correcta se muestran en la Figura 1.
El contenido del bloc de notas se exporta de la siguiente manera.
(1) El elemento de dirección externa de la foto es:
xs =-1038.485358; ys =-612.3039098 zs = 533.895357.
(2) Los elementos de orientación interna de la foto son:
f = 1.129641785; Ángulo de inclinación lateral: -0.02851629385;
(3) La matriz de rotación r es:
0.9994657309 -0.0241878564
-0.02114962887
0.02551734651 0.9992674772
0.02851971209
0.02115172777 -0.0290441573
0.9993694608
(4) Los resultados de la evaluación de precisión son los siguientes:
El error por unidad de peso es (micras):
3.3878e+005.
Los resultados del cálculo del corte izquierdo se muestran en la Figura 2.
El contenido del bloc de notas se exporta de la siguiente manera.
(1) Los elementos de orientación externa de la foto son:
xs =-1038.430617; ys =-612.276523 Zs = 533.6477522.
(2) Los elementos de orientación interna de la foto son:
f =-0.7559497103; Ángulo de inclinación lateral: -0.01600371623;
(3) La matriz de rotación r es:
0.999605192 -0.007871353955
-0.02697224397
0.008304493335 0.9998377524
0.01598449904
0.02697478945 -0.01620217905
0.9995083761
(4) Los resultados de la evaluación de precisión son los siguientes:
El error por unidad de peso es (micras):
5.4711e+005
3 Conclusión
Según el análisis de coordenadas del punto de la imagen, el punto medio de la foto es el origen , el eje vertical es el eje Y y el eje horizontal es el eje X. Existe una gran brecha entre los resultados calculados y los reales, especialmente cuando el valor f es negativo. Después de utilizar el procedimiento de inspección de datos estándar, aunque hay una cierta brecha en los elementos de orientación externos, especialmente xs e ys, los valores f. son similares. Si utiliza sus propios datos para calcular, la desviación es demasiado grande y se convierte en un error. Intente realizar varias modificaciones en los datos de coordenadas del punto de control de entrada, pero hay algunos problemas. El número de iteraciones será 1 o 29 y la diferencia de tamaño es demasiado grande. Por lo tanto, los problemas con la precisión de los datos o las diferencias en los sistemas de coordenadas de los datos a menudo traen consecuencias erróneas para el trabajo posterior.