Radio de Curvatura y Formulas

IX CONGRESO INTERNACIONAL DE CIENCIAS DE LA TIERRA
( SANTIAGO DE CHILE, 6 10 NOVIEMBRE, 2006 )
PONENCIA: REDUCCIN Y AJUSTE POR MINIMOS CUADRADOS DE

DISTANCIAS APLICANDO EL ANALISIS VECTORIAL.
AUTOR:
FIS. MARIO ALBERTO CRUZ DAZ
DEPENDENCIA:
INSTITUTO DE CATASTRO DEL ESTADO DE PUEBLA
DOMICILIO:
AVE. 11 ORIENTE NUM. 2003

COL. AZCARATE, C.P. 72501
PUEBLA , PUE.
MXICO
TELFONO:
222 2340948
FECHA:
OCT-2006
Ponencia: Reduccin y ajuste por mnimos cuadrados de distancias aplicando el anlisis vectorial / Fis. Mario Alberto Cruz Daz
I . INTRODUCCIN
Entre las atribuciones del Instituto de Catastro del Estado de Puebla, se encuentra
la de realizar levantamientos topogrficos y geodsicos de predios rsticos y urbanos en el
Estado de Puebla, con el propsito de ligar las observaciones a la Red Geodsica Estatal para
alcanzar la precisin necesaria en la georreferenciacin de los predios medidos .
Se ha observado que en algunos casos, se tienen variaciones sobre la cartografa
al ubicar los predios medidos en campo, y que pueden ser resultado de observaciones de
campo deficientes de inconsistencias en la misma cartografa. Esta problemtica dio origen a
un estudio para analizar el status de la cartografa. Este consiste bsicamente en comparar las
coordenadas de los puntos determinados por mtodos satelitales (GPS), con las coordenadas
de los correspondientes puntos de la cartografa .
Debido a que no es posible colocar un equipo con tecnologa satelital en los
vrtices que conforman los predios y que las construcciones impiden la recepcin de las seales
electromagnticas de los satlites, se propone como solucin el mtodo de reseccin de
distancias.
El Teorema de Taylor lineariza las funciones no lineales que se presentan en el
mtodo de reseccin, induciendo a un modelo matemtico de mnimos cuadrados, que para
su solucin requiere de la posicin del punto en primera aproximacin .
Con base en lo anterior, la Tesis que se presenta en este documento, consiste en
aplicar el producto escalar de vectores para establecer esta aproximacin, lo cual nos ha
permitido la convergencia de las distancias, en una sola iteracin. El principio de este mtodo
tiene su fundamento terico en el producto escalar de vectores unitarios a lo largo de las
distancias entre los vrtices con posicin conocida y los vrtices cuyas coordenadas son
desconocidas.
Es importante mencionar que el ajuste se realiza sobre la Proyeccin
Cartogrfica UTM, lo cul implica, que las distancias medidas sobre la superficie terrestre, deben
ser reducidas a esta proyeccin, antes de aplicarse el modelo matemtico de Mnimos
Cuadrados. En Geodesia se ha adoptado al elipsoide como la figura geomtrica que describe
la forma de la Tierra desde un punto de vista matemtico.
II . PLANTEAMIENTO CONCEPTUAL
El radio de curvatura de una funcin en una vecindad de algn punto contenido
en ella, se define como el radio del crculo osculador que toca al punto; y debido a las
caractersticas geomtricas del elipsoide, para cada punto en la superficie de esta figura, se
tiene un nmero infinito de radios de curvatura, sin embargo, para propsitos Geodsicos se
consideran nicamente los radios de curvatura mximo y mnimo, y que se conocen como
radio de curvatura en el primer vertical ( N ) y en el meridiano ( M ) respectivamente.
El radio de curvatura de una funcin en una vecindad de algn punto contenido
en ella, se define como el radio del crculo osculador que toca al punto; y debido a las
caractersticas geomtricas del elipsoide, para cada punto en la superficie de esta figura, se
tiene un nmero infinito de radios de curvatura, sin embargo, para propsitos Geodsicos se
consideran nicamente los radios de curvatura mximo y mnimo, y que se conocen como
radio de curvatura en el primer vertical ( N ) y en el meridiano ( M ) respectivamente. Estos
radios de curvatura nos permiten reducir las distancias medidas en campo, a las distancias de
cuerda en el elipsoide, y para lo cual se re -quieren las alturas geodsicas ( elipsoidales ) de
ambos puntos, y un radio terrestre promedio que este en funcin de ambos radios, para ste
propsito se considera al Radio Medio Gaussiano.
Una proyeccin cilndrica transversa, es aqulla en la que la colocacin del
cilindro secante al elipsoide es girada 90 y representa la figura donde se proyectan los puntos
del elipsoide, para posteriormente desarrollar el cilindro, dando origen as a la proyeccin UT.M.
Entonces, existe un paso intermedio entre la distancia de cuerda y la distancia
sobre la proyeccin U.T.M., que consiste en transformar la distancia de cuerda a la distancia en
el elipsoide para lo cual se hace una aproximacin esfrica con el radio medio Gaussiano
como el radio de la esfera.
Una vez que las distancias de cuerda han sido transformadas en distancias
elipsoidales, se procede a proyectarlas sobre el cilindro, en donde se presenta una distorsin en
las magnitudes lineales, que aumenta a medida que se aleja del meridiano central, por este
motivo, se aplica un factor de escala en el meridiano central ( k0 = 0.9996 ), que se introduce

en la ecuacin del factor de escala que es funcin de la latitud del punto.
Al desarrollar las ecuaciones de distancias en el modelo de Reseccin, se llega a
un modelo matemtico de Mnimos Cuadrados, cuya solucin esta en funcin de la matriz de
diseo ( A ) y la matriz de valores observados ( L ).
III. RADIOS DE CURVATURA

El radio de curvatura en el meridiano ( M ) se deduce a partir de la definicin
clsica de radio de curvatura de una curva.
Donde es el radio de curvatura que esta dado por :
1+
( )]
2
dz
3
2
dx
2
d z
dx
Donde ( dz / dx ) es la derivada de la lnea tangente al elipsoide en el punto, y (

d2z
dx2 )
representa la segunda derivada.

( X , Z ) son las coordenadas cartesianas del punto situado en la superficie del
elipsoide y que estn en funcin de la latitud geodsica del punto. El desarrollo de estas
derivadas induce la expresin matemtica para el radio de curvatura M .
Esquemticamente el radio de curvatura en el meridiano ( M ) se puede visualizar
de la siguiente manera :
La expresin matemtica para M es:
M=
a(1e2)
( 1 e 2 sen 2 )3/2
El radio de curvatura en el primer vertical ( N ) se deduce a travs del radio del

paralelo:
De la figura anterior se tiene que:
a
N=
( 1 e 2 sen 2 )1/2
II .1. Radio Medio Gaussiano.
Este radio se define como el valor integral medio de R tomado sobre el acimut de
0 a 360:
R=
1
2
1
R d =
2
0
M
M sen 2 + N cos 2
Cambiando el lmite de integracin a / 2 se tiene que:
R=
/2
N cos 2
M sen 2
N cos 2
Finalmente tenemos que el Radio de Curvatura Medio Gaussiano esta dado por:
R=
a [ 1 e 2 ] 1/2
( 1 e 2 sen 2 )
III . FACTOR DE ELEVACIN .
El radio de curvatura medio Gaussiano y las alturas geodsicas son los

parmetros que definen el factor de elevacin y aplicando la ley de cosenos en la siguiente
figura:
l
h1
h2
l0
Se tiene que:
l0 =
l2 - h2
1 +
h1
R
1 +
h2
R
1
2
Donde l0 es la distancia medida en campo que ha sido reducida a la

distancia de cuerda en el elipsoide, ahora se requiere transformar la distancia de cuerda a la
distancia en la superficie del elipsoide S ; para este propsito se utiliza una aproximacin
esfrica cuyos parmetros son, el Radio Medio Gaussiano y el ngulo entre las normales al
elipsoide en cada uno de los puntos, as tenemos que:
a [ 1 e 2 ] 1/2
S=
( 1 e sen )
cos -1
1-
l0 2 ( 1 e 2 sen 2 ) 2
2a [ 1 e 2 ]
IV . FACTOR DE ESCALA
Las distancias en el elipsoide son proyectadas sobre un cilindro que es secante al
elipsoide de referencia, y que posteriormente se desarrolla dando origen a la proyeccin U.T.M.,
en la cual, el cilindro se coloca transversalmente, es decir, el eje del cilindro, se encuentra en el
Ecuador.
La imagen del meridiano central es una lnea recta; todos los dems meridianos
son lneas curvas que convergen en el polo y son perpendiculares a la imagen del Ecuador, el
cual coincide con el eje X. Por supuesto que los paralelos mapeados son perpendiculares a los
meridianos mapeados.
El factor de escala k0 es un elemento de diseo importante, ya que si por
ejemplo se analiza la magnitud de un elemento de mapeo x, se encuentra que la magnitud se
incrementa a medida que se aleja del meridiano central, siendo despreciable en el meridiano
central. Al seleccionar un valor de k0 < 1, se permite una pequea distorsin en el meridiano
central, para as, tener menos distorsin al alejarse del meridiano central, de esta manera, el
cubrimiento longitudinal del rea de un mapa puede extenderse para un nivel de distorsin
aceptable.
Debido a las caractersticas de esta Proyeccin, es necesario calcular un factor
de escala que es funcin de la latitud, y que esta dado por:
k = k0 [ 1 + ( 2 / 2 ) cos2 ( 1 + 2 ) + 4 ( cos4 / 24 ) ( 5 4t2 +14 2 +

13 4 28 t2 2 + 4 6 48 t2 4 24 t2 6 ) + 6 ( cos6 / 720 )
( 61 148 t2 + 16 t4 ) + ..........]
Donde;
t = tan
2 =
e2 cos2
1 e2
Entonces, la distancia U.T.M. se calcula con la siguiente expresin:

dUTM = k . s
V . TEORIA DE MINIMOS CUADRADOS .

Sea el siguiente sistema de ecuaciones de observacin:
a11 x1 + a12 x2 + .................+ a1n xn = L1 + V1

a21 x1 + a22 x2 + .................+ a2n xn = L2 + V2
.
.
am1 x1 + am2 x2 + ...............+ amn xn = Lm + Vm
Este sistema de ecuaciones se representa en notacin matricial como:
a11 a12 ...........a1n

a21 a22 ...........a2n
.
.
.
.
.
am1 am2...........amn
x1
x2
.
.
.
.
.
xn
L1
L2
.
.
.
.
.
Lm
V1
V2
.
.
.
.
.
Vm
La solucin de mnimos cuadrados minimiza la suma de los residuos elevados al

cuadrado:
V2 = min
10
En lgebra matricial se minimiza V T V .La teora del Clculo Diferencial establece

que el mnimo de una funcin se obtiene al derivar dicha funcin con respecto a las variables e
igualando el resultado a cero, entonces:
( VTV)
X
donde
2 VT
V
X
V=AXL
V
X
=>
(AXL)
X
= A
=>
2V A=0
Transponiendo las matrices y multiplicando por un medio, se tiene que:
AT V = 0
pero
V = AX L
AT ( AX L ) = 0
( )( 2 ) ( VT A )T = 0
AT A X AT L = 0
AT A X = AT L
( AT A )1 ( AT A ) X = ( AT A )1 AT L
X = ( AT A )1 AT L
V I . PRODUCTO ESCALAR .
Por definicin, el producto escalar de dos vectores esta dado por:
a . b = ( a1 , .............., an ) . ( b1 , ............., bn )
11
= a1 b1 + ..............................+ an bn
Ahora consideremos la proyeccin de un vector
vectoriales unitarias
sobre sus componentes
b, y b :
Sea s la componente del vector a sobre b, es decir:
s = Comp b a
Sea t la componente del vector a sobre b, es decir:
t = Comp b a
=>
a = s b + t b = ( Comp b a ) b + ( Comp b a ) b
= Proy b a + Proy b a
Realizando el producto escalar de la expresin anterior con el vector b, se tiene

que:
s=
a.b
b.b
a.b
b
Donde II b II es la magnitud del vector b, y de la figura anterior:
a . b = II a II II b II cos
12
V I I . METODO DE RESECCIN VECTORIAL .

Para dar solucin al modelo matemtico de mnimos cuadrados, es necesario
conocer las coordenadas aproximadas del punto cuya posicin se va a calcular; dicha
posicin se determina con el producto escalar de vectores unitarios a lo largo de las direcciones
entre los puntos fijos que han sido georreferenciados por posicionamiento G.P.S.:
Y
P
d0
d1
e
A
d2
e2
e 1
b2
b1
B
X
Donde; e, e1, y e2, son los vectores unitarios en la direccin de ( A a P ) , ( A a B )
y ( A a C ) respectivamente, y cuyas componentes son:
e = ( ex , ey )
e1 = ( ex1 , ey1 )
e2 = ( ex2 , ey2 )
d0, d1, y d2 son las distancias medidas en campo, que han sido reducidas a la
proyeccin cartogrfica U.T.M., y las distancias b1 y b2, se calculan en la misma proyeccin:
13
b1 = [ ( xB xA ) + ( yB yA ) ]
b2 = [ ( xC xA ) + ( yC yA ) ]
De la figura anterior, se tiene por ley de cosenos que:
d02 + b12 d12
e . e1 =
2 d0 b1
d02 +b22 d22
e . e2 =
2 d0 b2
Por lo tanto:
ey1
ey2
( ex , ey ) . ( ex ,
( ex , ey ) . ( ex ,
)=
d02 + b12 d12

2 d0 b1
)=
d02 +b22 d22

2 d0 b2
)=
d02 + b12 d12

2 d0 b1
=>
1
( ex ex +
ey ey1
( ex ex2 + ey ey2 ) =
d02 +b22 d22

2 d0 b2
En las ecuaciones anteriores, las incgnitas son las componentes ( ex , ey ) del

vector unitario e, donde las componentes de los vectores unitarios estn dadas por:
e1 = ( ex1 , ey1 ) =
( xB xA )
b1
( yB yA )
b1
14
e2 = ( ex ,
ey2
)=
( xC xA )
( yC yA )
b2
b2
En notacin matricial, se tiene que:
( xB xA )
( yB yA )
b1
b1
( xC xA )
( yC yA )
b2
b2
d02 + b12 d12
ex
2 d0 b1
=
ey
d02 +b22 d22

2 d0 b2
Que es el modelo matemtico de mnimos cuadrados en notacin matricial, y

cuya solucin esta dada por:
X=(A A)
*A L
Entonces, si es el ngulo de direccin del vector e , se tiene que:
tan =
ex
ey
Y las coordenadas aproximadas para el punto P se calculan con las siguientes

expresiones:
15
XP = XA + d0 * sen ( )
YP = YA + d0 * cos ( )
Con las coordenadas aproximadas del punto P ( que se obtienen con el
producto escalar de vectores unitarios ) se calculan las coordenadas precisas de P ,
considerando el Teorema de Taylor.
En la siguiente figura se muestran las distancias U.T.M. que se utilizan en el ajuste
por mnimos cuadrados por el mtodo de reseccin:
Y
P ( XP,YP )
d0
A ( X0,Y0 )
d2
d1
B ( X1,Y1 )
C ( X2,Y2 )
X
Las distancias ( do , d1 , y d2 ) son funciones de las coordenadas de los vrtices

conocidos y libres de error ( A , B , y C ), as como de las coordenadas del punto P ( XP , YP )
desconocidas:
16
f0 ( XP0 , YP0 , X0 , Y0 ) = [ ( XP0 X0 ) 2 + ( YP0 Y0 ) 2 ] = d0

f1 ( XP0 , YP0 , X1 , Y1 ) = [ ( XP0 X1 ) 2 + ( YP0 Y1 ) 2 ] = d1
f2 ( XP0 , YP0 , X2 , Y2 ) = [ ( XP0 X2 ) 2 + ( YP0 Y2 ) 2 ] = d2
El superndice ( 0 )en las coordenadas de P indica que son aproximadas, y las
cuales sern calculadas con precisin en el modelo matemtico de mnimos cuadrados.
Las ecuaciones de distancia anteriores son funciones no lineales, y para
linearizarlas se aplica el Teorema de Taylor: Toda funcin no lineal, que admite las primeras
derivadas, se aproxima por la funcin evaluada en valores aproximados ( f 00 ), ms las primeras
derivadas de la funcin con respecto a las variables :
f0
0
f0 = f0 +
f0
d XP +
YP
XP
f1 = f10 +
f2 = f20 +
f1
d XP
XP
f2
d YP
f1
d YP
YP
0
d XP
XP
f2
d YP
YP
0
Desarrollando las derivadas parciales de las funciones anteriores, y rearreglando

los resultados en notacin matricial se tiene que:
17
XP0 X 0
d00
XP0 X1
d10
XP0 X2
d20
YP0 Y0
d00
f0 f00
d XP
YP0 Y1
f1 f10
d10
YP0 Y2
d20
d YP
f2 f20
Estas matrices presentan el diseo del modelo matemtico de mnimos

cuadrados, que tiene como solucin:
X = ( AT A ) 1 AT L
Donde la matriz solucin representada por X , tiene como componentes los
elementos diferenciales que se suman algebraicamente a las coordenadas aproximadas del
punto P :
Es decir,
XP = XP0 + d XP
YP = YP0 + d YP
Este proceso contina hasta que la solucin converge.
18
CONCLUSIONES.
El mtodo vectorial propuesto para la georreferenciacin de puntos del terreno,
tiene como fundamento terico, el producto escalar de vectores para la determinacin de las
coordenadas de los puntos en primera aproximacin, que son utilizados en el ajuste por mnimos
cuadrados.
Este mtodo se aplica en las distancias medidas que han sido reducidas a la
proyeccin U.T.M., por lo que para cada punto se tiene que las incgnitas son las coordenadas
( X , Y ), y ya que la Teora de mnimos cuadrados establece que debe existir redundancia en la
solucin ( siendo la redundancia igual al nmero de ecuaciones menos el nmero de incgnitas
), entonces se deben medir al menos tres distancias a partir de tres puntos diferentes ( puntos
G.P.S. ) hacia los puntos en consideracin.
Se realizaron estudios de campo para probar este mtodo encontrndose que:

*Las coordenadas obtenidas tienen una aproximacin del orden de 2 a 3 cms., y en el
mejor de los casos de 2 a 3 mm. con respecto a las coordenadas ajustadas por el
mtodo de mnimos cuadrados.
*Al introducir las coordenadas obtenidas por este mtodo como primera aproximacin
en el ajuste por mnimos cuadrados, la solucin converge en la primera iteracin.
*Se reduce el tiempo de observacin, ya que solo se miden distancias y no es necesario
medir adicionalmente ngulos, como se lleva a cabo en un levantamiento geodsico
tradicional.
*Ya que resulta imposible el colocar un equipo de posicionamiento satelital clsico en las
esquinas de predios urbanos, y debido a que se tienen cortinas que impiden la
recepcin de las seales satelitales, este mtodo puede resultar una opcin adecuada
para la georreferenciacin precisa de estos elementos de construccin.
19
*Ya que en nuestro Pas ( Mxico ), y en especial en el Estado de Puebla, la superficie

terrestre sufre de accidentes topogrficos de consideracin, nos encontramos
frecuentemente con sitios que resultan inaccesibles para las mediciones de campo, y
esta resultara otra de las aplicaciones para este mtodo.
*Adicionalmente, se tiene que ya que la solucin se da por un modelo matricial
bastante simple, resulta obvio que este mtodo reduce el tiempo de proceso de
computo, reflejandose en un ahorro en los recursos econmicos asignados para esta
actividad.
BIBLIOGRAFIA .
*FOWLES,G.R.,ANALYTICAL MECHANICS,SAUNDERS COLLEGE PUBLISHING,USA,1993.

*TORGE,W.,GEODESY,WALTER DE GRUYTER INC.,U.S.A.,1991.
*ZAKATOV,P.S.,CURSO DE GEODESIA SUPERIOR,EDITORIAL MIR,MOSCU,1997
*LIPSCHUTZ,S.,LINEAR ALGEBRA,Mc GRAW-HILL,U.S.A.,2000
*GANTMACHER,F.R.,APPLICATIONS OF THE THEORY OF MATRICES, DOVER
PUBLICATIONS,U.S.A.,2005.
*RAPP,R,GEODESIA TERICA AVANZADA, NOTAS DE CLASE, FORT
CLAYTON,PANAMA,1985.
20

Radio de Curvatura y Formulas

Cargado por

Copyright:

Formatos disponibles

Radio de Curvatura y Formulas

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Radio de Curvatura y Formulas

Cargado por

Copyright:

Formatos disponibles

IX CONGRESO INTERNACIONAL DE CIENCIAS DE LA TIERRA

( SANTIAGO DE CHILE, 6 10 NOVIEMBRE, 2006 )

PONENCIA: REDUCCIN Y AJUSTE POR MINIMOS CUADRADOS DE

FIS. MARIO ALBERTO CRUZ DAZ

INSTITUTO DE CATASTRO DEL ESTADO DE PUEBLA

AVE. 11 ORIENTE NUM. 2003

motivo, se aplica un factor de escala en el meridiano central ( k0 = 0.9996 ), que se introduce

III. RADIOS DE CURVATURA

Donde es el radio de curvatura que esta dado por :

Donde ( dz / dx ) es la derivada de la lnea tangente al elipsoide en el punto, y (

representa la segunda derivada.

La expresin matemtica para M es:

El radio de curvatura en el primer vertical ( N ) se deduce a travs del radio del

De la figura anterior se tiene que:

Cambiando el lmite de integracin a / 2 se tiene que:

III . FACTOR DE ELEVACIN .

El radio de curvatura medio Gaussiano y las alturas geodsicas son los

Donde l0 es la distancia medida en campo que ha sido reducida a la

k = k0 [ 1 + ( 2 / 2 ) cos2 ( 1 + 2 ) + 4 ( cos4 / 24 ) ( 5 4t2 +14 2 +

Entonces, la distancia U.T.M. se calcula con la siguiente expresin:

V . TEORIA DE MINIMOS CUADRADOS .

a11 x1 + a12 x2 + .................+ a1n xn = L1 + V1

a11 a12 ...........a1n

La solucin de mnimos cuadrados minimiza la suma de los residuos elevados al

En lgebra matricial se minimiza V T V .La teora del Clculo Diferencial establece

sobre sus componentes

Sea s la componente del vector a sobre b, es decir:

Realizando el producto escalar de la expresin anterior con el vector b, se tiene

Donde II b II es la magnitud del vector b, y de la figura anterior:

V I I . METODO DE RESECCIN VECTORIAL .

De la figura anterior, se tiene por ley de cosenos que:

d02 + b12 d12

d02 + b12 d12

d02 +b22 d22

d02 + b12 d12

d02 +b22 d22

En las ecuaciones anteriores, las incgnitas son las componentes ( ex , ey ) del

En notacin matricial, se tiene que:

d02 + b12 d12

d02 +b22 d22

Que es el modelo matemtico de mnimos cuadrados en notacin matricial, y

Entonces, si es el ngulo de direccin del vector e , se tiene que:

Y las coordenadas aproximadas para el punto P se calculan con las siguientes

Las distancias ( do , d1 , y d2 ) son funciones de las coordenadas de los vrtices

f0 ( XP0 , YP0 , X0 , Y0 ) = [ ( XP0 X0 ) 2 + ( YP0 Y0 ) 2 ] = d0

Desarrollando las derivadas parciales de las funciones anteriores, y rearreglando

Estas matrices presentan el diseo del modelo matemtico de mnimos

Se realizaron estudios de campo para probar este mtodo encontrndose que:

*Ya que en nuestro Pas ( Mxico ), y en especial en el Estado de Puebla, la superficie

*FOWLES,G.R.,ANALYTICAL MECHANICS,SAUNDERS COLLEGE PUBLISHING,USA,1993.

También podría gustarte