Ing. [Link]. Edgar Montenegro T.
UNIVERSIDAD MAYOR DE SAN
SIMON
Facultad de Ciencias y Tecnología
Ing. - Civil
PRÁCTICA 6
UNIVERSITARIO:
SORIA BECERRA NELSON
MATERIA: GEOMATICA
DOCENTE: Ing. EDGAR LUIS MONTENEGRO TERRAZAS
GRUPO: 2
FECHA: 30/05/2022
COCHABAMBA - BOLIVIA
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
GEOMATICA PRACTICA N° 6
Postproceso de datos DGPS, Asignación de atributos a polígonos mediante tablas, selección
de datos por cruce con zona de interés
1. Para un trabajo de estudio de inundaciones en el río Rocha, se utilizará topografía satelital,
imágenes satélite y un levantamiento topográfico tradicional, motivo por el que se desea
colocar toda la información en un sistema de coordenadas común.
La estrategia empleada para georeferenciar la imagen consiste en identificar puntos que se
puedan reconocer fácilmente en el terreno y la imagen, determinando sus coordenadas
correctas con el uso de un GPS diferencial (DGPS)y su postproceso.
Después de identificar los puntos GCP en el Google Earth, se obtiene las coordenadas
mostradas en el Cuadro 1.
Cuadro 1: Puntos de control de tierra (GCP Ground Control Points), proyecto Río Rocha
punto Latitud Longitud
GCP_1 -17.3968 -66.3295
GCP_2 -17.4405 -66.3112
GCP_3 -17.4417 -66.2687
GCP_4 -17.3934 -66.2468
PUNTO LATITUD LONGITUD
GCP1 -17.396806° -66.329439°
GCP2 -17.440444° -66.311344°
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
GCP3 -17.441708° -66.268751°
GCP4 -17.393430° -66.246841°
Debido a que cada versión de Google Earth que usted esta manejando puede variar en las
coordenadas, verifique que los GCP, evidentemente coinciden con los sitios previamente
seleccionados y que se encuentran ilustrados en la carpeta “ilustraciones GCP”, de ser
necesario muévalos a su posición correcta, realice capturas que permitan evaluar la
coincidencia mencionada, proporcione un Cuadro equivalente al Cuadro 1, con las
posiciones correctas de su imagen de Google Earth.
GCP1
GCP2
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
GCP3
GCP4
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
Después de levantar datos con el Rober del DGPS, durante 4 horas como mínimo en cada
punto y en vista de que no se tubo acceso a un punto geodésico confiable, los datos de cada
GCP en formato “rinex” fueron sometidos a un post proceso de la herramienta AUSPOS,
de la entidad denominada “Geociencia Australia” del gobierno de Australia, que tiene el
siguiente Link: [Link]
1.1) Usar el el servicio AUSPOS y procesar el “rinex” del punto GCP 4, que se encuentra
en la carpeta: GCP4-VICUSTARIZ y completar el Cuadro 2, en el que figuran las
coordenadas de otros puntos, que fueron procesados con la misma metodología.
1.2) Usar también el servicio de post proceso del IBGE, que tiene el siguiente link:
[Link]
geodesico/servicos-para-posicionamento-geodesico/16334-servico-online-para-pos-
[Link]?=&t=processar-os-dados . Debido a que el
punto GCP4, no se encuentra en territorio Brasilero el servicio del IBGE, no reporta la
altura geoidal, en consecuencia obtener esta altura a través del siguiente link
[Link] , usar para el efecto el modelo geoidal egm
2008 de 1´ de resolución y calcular la altura ortometrica del punto GCP 4.
´+ Altura elipsoidal Altura geoidal Altura ortometrica msnm
GCPS 2589,36 45,01828 2.544
1.3) Presentar un Lay out con las alturas geoidales para toda Bolivia.
Cuadro 2: Coordenadas obtenidas con DGPS y postproceso con AUSPOS
DGCPS LONGITUDE LATITUDE WGS-84 Ellipsoidal Height (m) MSL Elevation (m)
1 W 66° 19' 46.02551" S 17° 23' 48.47686" 2593.496 2547.749
2 W 66° 18' 40.51585" S 17° 26' 25.81157" 2555.15 2509.827
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
3 W 66° 16' 7.47236" S 17° 26' 30.09594" 2593.416 2548.452
4
Ing. [Link]. Edgar Montenegro T.
4 W66°14’48.6301” S 17°23’36.27740” 2589.354 2544.324
El procedimiento para el proceso con AUSPOS, se encuentra en el archivo denominado
“procedimiento”
Adicionalmente, responder lo siguiente:
a) Incorporar los DGCPS del Cuadro 2 en el Google Earth, elaborar una tabla con las
coordenadas UTM en WGS84 Z19 y determinar las diferencias en ΔX y ΔY de la
georeferencia de su Google Earth.
coordenadas GCP coordenadas Google earth
PUNTO X Y X Y ERROR X(m) ERROR Y(m)
GCP1 783736,65 8074570,08 783741,75 8074570,11 5,1 0,03
GCP2 785614,14 8069703.94 785608,19 806903,68 5,95 0,26
GCP3 790129,55 8069506,97 790126,25 8069506,57 3,3 0,4
GCP4 792534,32 8074822,06 792530,55 8074818,35 3,77 3,71
PROMEDIO 4,53 1,1
b) Cual el motivo de las diferencias reportadas en el ítem a.
c) Si usted desea georeferenciar con mucha precisión una imagen Google Earth de esta
zona, como lo haría, describir el procedimiento que aplicaría y que puntos de amarre
usaría ??.
R.-Para poder georeferenciar la imagen con una buena precisión se debería adquirir y usar
un GPS diferencial y realizar todos los trabajos cin este GPS para poder tener una mejor
precisión, el punto de amarre se usaría alrededor de cada punto buscado.
d) Describa los métodos que podría usar para determinar las coordenadas (x,y,z) de un
punto mediante el uso del GPS diferencial.
R.- La navegación autónoma con u n receptor simple, la precisión es de 20 metros para
usuarios militares y 100 metros para usuarios civiles. Utilizado en navegación marina
posicionamiento diferencial corregido (DGPS) con precisiones de 0.5 posicionamiento de
vehículos. Posicionamiento entre 0 diferencial de fase con precisiones.520, utilizado en
control de maquinaria y topografía. Para determinar las coordenadas (x,y,z) de un punto
utilizaríamos cualquiera de los tres métodos aprendidos en clases de geomática.
Primer método.- Para obtener datos precisos, debemos contar con un punto geodésico de
primer orden o antena base que se encuentra en cotapachi y un Robert para después
procesar los datos obtenidos en nuestro Rinex.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
Segundo método.- Utilizando el servidor AUSPOS el servidor auspos utiliza puntos
geodésicos regionales, procesa los datos simultáneos y obtenemos (x,y,z).
Tercer método.-Para este método debemos contar con una antena base y un Robert que se
comunican mediante una radio (RTK) y hacemos un postproceso con los dos archivos y asi
obtendremos nuevo punto (x´,y´,z´).
e) Cuál es la altura que proporciona un DGPS?
R.-Posicionamiento, Diferencial, Corregido. Más comúnmente conocido como DGPS,
el cual proporciona precisiones del orden 0,5-5 m.
f) Que es la altura ortométrica ?
R.- Es la diferencia entre la altura elipsoidal (GPS) y altura geoidal
g) La altura del punto proporcionada por postproceso del IBGE, que en su reporte es
denominada Alt. Geo.(m), que altura es?, por que hace esa afirmación?
R.-La altura proporcionada del IBGE Alt. Geo(m) es la Altura ELIPSOIDAL es 2.589,36.
h) Explicar los aspectos más destacados para la determinación de las coordenadas (x,y,z)
de un punto, mediante el uso de un GPS + RTK-NTRIP
R.-Es cuando tienes un dispositivo lleno de internet y conectas con una antena empírica
interceptándolo ya conoceríamos un punto. La única diferencia seria que el costo sería poco
más elevado.
i) Con relación a la pregunta (h) cual sería la diferencia del método conocido como GPS
RTK
R.-La única diferencia seria vía radio en la que la señal llegaría a perderse en cualquier
momento y adicional a ello es que el monto es también poco más elevado.
2. Sobre la imagen [Link], digitalizar polígonos para representar los lotes mostrados
en la Figura 1. Durante la digitalización usar la opción autocompletar a fin de no
tener espacios entre los límites de cada lote.
Figura 1: Identificación de lotes de interés
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
2.1. A cada uno de los lotes presentados anteriormente introducir en una tabla de ArcGis
las características mostradas en el Cuadro 3, elaborar para el efecto una tabla en Excel
e importar a la tabla de atributos de ArcGis, mostrar una captura de la tabla de atributos
donde se pueda visualizar la columna de vinculo Excel-ArcGis utilizada.
Además incorporar en la tabla una columna donde muestre el área de cada lote, presentar
la captura de pantalla correspondiente para verificar la conclusión de esta etapa.
Cuadro 3: Características de los lotes mostrados en la Figura 1
Lote Telefono pro_excel Tipo_ propiedad Tipo Suelo Precio U$/m2
1 71433490 Sofia Orellana Propio Arcilloso 150
2 67325215 Guido Nogales Herencia Arcillo arenoso 120
3 77424213 Octavio Quiroga Propio Arcillo siltoso 100
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
2.2. Con la opción idetify de ArcGis, pinchar el lote tres, mostrando sus
características, capturar la pantalla correspondiente.
2.3. Exporta la tabla de atributos completa de ArcGis a Excel e Indicar el presupuesto
requerido para adquirir cada lote y el monto total para los tres lotes.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
TIPO PRECIO
LOTE TELEFONO PRO EXCEL PROPIRDAD TIPO SUELO U$/m2 Área precio (U$)
1 71433490 Sofía Orellana Propio Arcilloso 150 9157,79 1373668,32
2 67325215 Guido Nogales Herencia Arcillo arenoso 120 5629,46 675534,785
3 77424213 Octavio Quiroga Propio Arcillo sitoso 100 10467,15 1046714,99
Monto
total 3095918,1
2.4. En función del error esperado en la estimación del área, cual considera que
podría ser el error en el presupuesto asignado para la compra de los tres lotes ?
Error en la Error en la Error en la
PRECIO estimación del estimación del estimación del
U$/m2 Área precio (U$) superficie (%) superficie (m2) presupuesto (U$)
150 9157,79 1373668,32 4,5 412,10 61815,07
120 5629,46 675534,785 4,5 253,33 30399,07
100 10467,15 1046714,99 4,5 471,02 47102,17
Monto
total 3095918,1
2.5. Que método y que instrumentos usaría para realizar una cuantificación precisa
de la superficie de los lotes ?
R.- Se puede usar el método de la estación total o aerofotometría para ser más preciso.
2.6. Elaborar un mapa en ArcGis, etiquetando cada lote con el nombre del
propietario.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
3. En un mapa que tenga solo los contornos de Bolivia mostrar la localización de la
cuenca de los ríos Ichilo y Mamore, mostrar la localización exclusivamente de las
estaciones “Hido-meteorologicas” que se encuentran al interior de esta cuenca,
incorporar la localización de las capitales de los nueve departamentos de Bolivia, así
como la red de drenaje principal, solo para la cuenca del Ichilo y Mamore.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
3.1. Indicar la superficie de la cuenca del río Ichilo y Memoré, expresando la misma
en Km2 y en m2.
3.2. Elaborar un Layout con todos los elementos que debe contener un mapa.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
3.3. Exportar la tabla de las estaciones a Excel y presentar un cuadro que contenga el
nombre de la estación, sus coordenadas en Longitud, Latitud y Altitud, así como
la institución a cargo de la colecta de datos.
4. Presentar en un solo mapa los centros de salud y vías secundarias del departamento que le
corresponda en función de la terminación de su número de Carnet de identidad, según lo
establecido en el Cuadro 4.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
Cuadro 4: Departamento que le corresponde en función de la terminación de su CI
Terminación CI Departamento Terminación CI Departamento
0, 1 Tarija 6 Cochabamba
2 Beni 7 Santa Cruz
3 Pando 8 La Paz
4 Potosí 9 Oruro
5 Chuquisaca
Agregar una foto en la que aparezca su nombre y su número de carnet
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
Elaborar el citado mapa a partir de las bases de datos existentes, proporcionar el link de
donde obtuvo el mapa mencionado.
Presentar su mapa en un layout con todos los elementos que debe tener un mapa
cartográfico.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
Ing. [Link]. Edgar Montenegro T.
5. Se desea realizar un estudio hidrológico en la cuenca delimitada por usted en
anteriores prácticas, indicar las estaciones meteorológicas del SENAMHI que usaría
para el efecto, presentar un mapa en el que se ilustre las estaciones por usted
adoptadas.
R.- Las estaciones próximas a nuestra cuenca de son las Estaciones: PUENTE ARCE,
PUENTE ARCE_GRANDE, TABOADA
6. Presentar sus conclusiones respecto a la presente practica.
R.- Las conclusiones se refieren tanto a la experiencia adquirida en la realización del
trabajo como a los resultados prácticos del mismo: este software llamado Arcgis tienen
demasiadas herramientas de las cuales en este proyecto no logramos utilizar todas,
solamente las necesarias como para sacar el área, etc. También descargamos archivos como
el AUSPOS, IBGE.
La finalidad de la práctica es aprender a aplicar los conocimientos transmitidos del docente
tanto como lo teórico y como lo práctico
D:\geomatica\lixo 6 2022\pra6_1_22\pra6_1_2022.doc
