MANUAL DE PROCESAMIENTO DE PUNTOS GPS CON SOFTWARE GAMIT/GLOBK
Balvín Velásquez, Alvaro M. Grijalva Calzada, Carlos A. Mendoza del Aguila, Mario C.
Dirección de Geodesia Dirección General de Cartografía Instituto Geográfico Nacional
3 de Junio de 2011
CONTENIDO
Introducción Procedimiento
INTRODUCCIÓN
GAMIT Es un software que procesa datos GPS para obtener estimaciones de coordenadas de las estaciones, parámetros orbitales, parámetros de rotación terrestre. El software está diseñado para ejecutarse en cualquier sistema operativo LINUX.
GLOBK Combina soluciones de GPS, SLR y VLBI (son sistemas de almacenamiento de datos y observaciones) obteniendo repetitividad de las coordenadas, ajuste de coordenadas a un Marco de referencia establecido, velocidades de puntos, ajuste de las órbitas de satélites, ajuste de parámetros de rotación terrestre,
El software GAMIT debe ser instalado en el sistema operativo Linux-Ubuntu. Todo procedimiento y pasos a seguir serán realizados mediante la consola del sistema operativo. Entorno del Sistema Linux (Desktop, PC)
Consola del Sistema Operativo Linux
Recomendamos al usuario organizar su información (data de los días a procesar, información de las estaciones IGS y convertidores de archivos) en una carpeta adicional en el escritorio para facilitar el trabajo.
PROCEDIMIENTO 1.1.
Abriremos la consola del sistema LINUX haciendo clic en recuadro negro que aparece en la parte superior del PC. Clic aquí
1.2.
Digitaremos el comando “ls” y presionaremos “ENTER” para la visualización de nuestro desktop.
1.3.
Digitaremos el comando “mkdir” para crear una carpeta donde guardaremos los datos de nuestro proyecto. Ej. “mkdir neww”, presionamos la tecla “ENTER”
1.4.
Digitaremos el comando “mkdir” para crear la carpeta “2011”. Ej. “mkdir 2011” y presionamos la tecla “ENTER”
1.5.
Digitaremos el comando “mkdir” para crear la carpeta “rinex”. Ej. “mkdir rinex” y presionamos la tecla “ENTER”.
1.6.
Digitamos el comando “cd” para entrar a la carpeta “rinex” y le damos “ENTER”. Ej. “cd rinex”.
1.7.
Digitamos el comando “pwd” para ubicar la ruta de la carpeta antes creada y le damos “ENTER”.
1.8.
Abrimos una segunda consola presionando el botón negro que mostramos inicialmente.
1.9.
Digitamos el comando “cd Escritorio”, para entrar al escritorio.
1.10. Escribimos el comando “cd” para entrar a la carpeta donde tenemos los datos a procesar. Ej. “cd ica”
1.11. Digitamos el comando “ls” y le damos “ENTER” para visualizar la data de la carpeta “ica”
1.12. Digitamos el comando “cp *”, para copiar la ruta de la carpeta “rinex” mostrada en la ventana anterior. Ej. “cp * /home/ign/ing11/2011/rinex”.
1.13. Regresamos a la ventana inicial y digitamos el comando “ls” y le damos “ENTER” para verificar que los datos hayan sido copiados.
1.14. Digitamos el comando “cd ..” y le damos “ENTER”.
1.15. Digitamos el comando “sh_setup –yr 2011” y le damos “ENTER” para crear la carpeta “tables”.
1.16. Digitamos el comando “cd” y le damos “ENTER”, para ingresar a la carpeta “tables”.
1.17. Digitamos el comando “ls” y le damos “ENTER”, para visualizar los datos de la carpeta creada.
1.18. Digitamos el archivo “process.defaults” con el comando “gedit process.defaults”.
1.19. Le damos “ENTER” y se deberá abrir la siguiente ventana.
1.20. Buscamos en la ventana o consola la ruta de la carpeta “rinex”, copiamos y pegamos en la fila “15” de la ventana de texto que se ha abierto.
1.21. Buscamos la fila “65” y digitamos el nombre del “ITRF” (MARCO INTERNACIONAL DE REFERENCIA TERRESTRE) con el que vamos a procesar.
1.22. Guardamos y cerramos la ventana.
1.23. En la segunda ventana crearemos un programa con el cual editaremos los “IGS” (Servicio Internacional del GPS), el comando a usar es “gedit prog.bash” y le damos “ENTER”.
1.24. Se abrirá la siguiente venta sin dato alguno.
1.25. En esta ventana digitaremos los siguientes comandos: Guardamos y cerramos la ventana.
neww
1.26. Digitaremos el comando “chmod +x” para volver ejecutable el “prog.bash” y poder editar sus archivos. Ej. “chmod +x prog.bash” y le damos “ENTER”
1.27. Volvemos a la primera ventana (consola), y editaremos el archivo “sites.defaults”. ej. “gedit sites.defaults” y le damos “ENTER”.
1.28. Se abrirá la siguiente ventana.
1.29. En esta ventana comentaremos desde la fila “26” hasta la “32” colocando el signo numeral “#” de la siguiente manera:
1.30. Luego desde la fila “34” hasta la “37”escribiremos el comando “glreps” de la siguiente manera: (solo colocamos este comando en la estaciones IGS mas no en la estación a determinar sus coordenadas, que en este caso es “ica01”).
1.31. Guardamos y cerramos el archivo. 1.32. Digitaremos los siguientes comando en la primera consola en el mismo orden (uno por uno) en el que se encuentran.
Damos “ENTER”
Damos “ENTER”
Damos “ENTER”
Damos “ENTER”
Con este último paso vemos en la consola que empieza a generar varios archivos e información, esperamos unos cuantos minutos que termine de cargar los datos y nos daremos cuenta que terminara cuando regresa a la carpeta “tables”; tal y como se muestra en la imagen.
1.33. Digitamos el comando “cd ..” para regresar a la carpeta “2011” y le damos “ENTER”.
1.34. Una vez en la carpeta “2011” digitaremos el comando “sh_gamit –expt neww –s 2011 038 042 –orbit IGSF cbpt x k p –dopt c” y le damos “ENTER” para dar inicio al procesamiento.
Luego de apretar “ENTER… empezará a correr el programa…
1.35. Luego de correr todo. Digitamos el comando “ls” y luego “ENTER” y aparecerán los días procesados (días establecidos en un principio). Ver figura:
1.36. Luego el siguiente paso consiste en copiar algunos archivos que están dentro del “gg” del escritorio. Para entrar ponemos “ls” y luego “ENTER”.
1.37. Luego entramos al “gg” con el comando “cd” es decir: “cd gg”.
1.38. Luego entramos al tables de la siguiente manera: “cd tables”
1.39. Luego entramos con “ls” y aparecen los archivos.
1.40. Luego el siguiente paso consiste en copiar los archivos “globk y glorg” ejecutando el comando así como muestra en la figura.
1.41. Luego de ejecutar ese comando inmediatamente se copiará en el tables del experimento neww (no cerrar esta consola). Comprobando sale lo siguiente.
Entonces vemos que se copio. (Esta resaltado en rojo).
1.42. Luego editamos “globk_comb.cmd”
1.43. Luego modificamos lo siguiente:
1.44. Luego editamos el “glorg_comb.cmd”:
1.45. Editamos la linea 6 para cambiar el marco de referencia que se va a utilizar, en la linea 25 escribimos “stab_site (nombres de los IGS a utilizar)” y comentamos la linea 28.
1.46 Luego guardamos los cambios que se hicieron en el archivo
1.47 Después cerramos este archivo y luego en la consola salimos de la carpeta de donde estamos, digitando “cd ..” y luego damos ENTER:
1.48 Ahora dentro de la carpeta “2011” digitamos el comando: “ sh_glred -s 2011 038 2011 042 -expt neww -opt H G E ”, luego damos ENTER aquí 038 y 042 son los días (día inicial hasta el día final) que se van a utilizar y neww es el nombre de la carpeta experimento.
1.49. Después que haya corrido totalmente se va a poder digitar “ls” y damos ENTER, esto es para verificar si se creó el archivo “disk.used” verificar:
1.50 Después entramos a la carpeta tables digitando “cd tables”
1.51
Una vez entrado a tables digitamos ls damos ENTER para enlistar los archivos:
1.52
1.53
Luego en la otra consola como estamos en la carpeta tables salimos de esta carpeta digitando “cd ..”, damos ENTER, y después creamos dos carpetas digitando “mkdir tables vsoln” y damos ENTER:
Una vez creado estas 2 carpetas entramos a tables digitando “cd tables” y luego damos ENTER:
1.54
1.55
Luego en la otra consola salimos del tables que está dentro de la carpeta 2011 digitando “cd ..”, damos ENTER y después digitamos “ls” y también damos ENTER:
Entramos a la carpeta del día inicial digitando “cd 038” y damos ENTER:
1.56
1.57
Ahora digitamos todo este comando: “sh_oneway –sites IC01 -autsum autcln.post.sum –prefix DPH –scale 0.5 – page L” y damos ENTER (aqui IC01 es la estación local):
Despues de haber corrido, otra vez digitamos otro comando parecido: “sh_oneway –sites IC01 -autsum autcln.post.sum –prefix DPH –scale 0.5 – page L -elev”, y luego damos ENTER:
1.58
Luego de haber corrido salimos de esta carpeta y nos vamos a la carpeta usuario digitando “cd” y damos ENTER:
1.59
Digitamos “cd gg” y damos ENTER para entrar a la carpeta gg:
1.60
Después entramos a la carpeta tables de “gg” digitando “cd tables”:
1.61 Enlistamos los archivos que se van a copiar a la carpeta tables que se creó últimos, para hacer esto digitamos “ls” y damos ENTER:
1.62
Copiamos algunos archivos digitando “cp IGS08_101117.eq itrf08.apr pmu.usno /home/ign/neww/tables”
1.63
También copiamos otros archivos a la carpeta vsoln digitando: “cp globk_vel.cmd glorg_vel.cmd /home/ign/neww/vsoln”
1.64
En la otra consola digitamos “cd /home/ign/neww/vsoln” damos ENTER y luego editamos los archivos copiados digitando “gedit globk_vel.cmd ”
1.65
Digitamos algunas líneas (escribimos igual de lo que está en círculos):
1.66
1.67.
Luego guardamos los cambios:
Cerramos el archivo, después editamos el otro archivo para eso digitamos “gedit glorg_vel.cmd” y damos ENTER:
1.68. Editamos algunas líneas mostradas
1.69. Luego guardamos los cambios
1.70. Entramos a la carpeta vsoln y digitamos los siguientes comandos: “ ls ../2011/glbf/*.glx > bien_2011.gdl ” y le damos “ENTER”
“ globk 6 PTR LOG bien 2011.gdl globk_vel.cmd ” y le damos “ENTER”
“ sh_plotvel ” y le damos “ENTER” “ sh_plotvel –f globk_vel.org –R244.5/251/40.5/46.5 ” y le damos “ENTER”
“ sh_plotvel –file PRT.org –ps test_vel_itrf –R-80/-30/-30/5 ” y le damos “ENTER”
“ gs test_vel_itrf.ps ” y le damos “ENTER”
“ sh_org2vel –file PRT.org –plate SOAM ” y le damos “ENTER”
“ ps test_vel_soam.ps ” y le damos “ENTER”
“ sh_plotvel –file PRT_SOAM.org –ps test_vel_soam –R-80/-30/-30/5” y le damos “ENTER”
“ ps test_vel_soam.ps ” y le damos “ENTER”
Luego de digitar todos estos comandos se muestra la siguiente pantalla:
1.71. Para observar las coordenadas UTM empecemos ingresando a la carpeta de experimento “neww” y abrir la carpeta “ vsoln”.
1.72. Dentro de la carpeta “vsoln” tenemos lo siguiente
1.73. Abrimos el archivo “PRT_SOAM.apr”
1.74. Entonces obtenemos las coordenadas UTM de la estación de ICA (ICA01)
1.75. Para observar las coordenadas Geográficas abrimos el archivo “PRT_SOAM.vel” de la carpeta “vsoln” abierto anteriormente.
1.76. Entonces obtenemos las coordenadas Geográficas de la estación de ICA (ICA01)
Magnitud de los errores de posición 1.77. Regresamos a la carpeta “neww” y abrimos la carpeta “2011”.
1.78. Ahora abrimos la carpeta “gsoln”
1.79. Dentro de la carpeta “gsoln” se encuentran todos los archivos que se muestra
1.80. Ahora hacemos clic sobre el archivo “psbase_neww.IC01”
1.81. El archivo muestra la magnitud de error en las coordenadas
BIBLIOGRAFIA 1. Introduction to GAMIT/GLOBK; T.A.Herring, R.W.King, S.C. McClusky. 2. Sistema y Marco de Referencia Terrestre; Angel Martín Furones. 3. Introduccion al sistema de Posicionamiento Global; Pedro Tipula T. and Osorio. 4. GPS & GLONASS “Descripcion y aplicaciones”; Juan Carlos Bermejo Ortega 5. Manual de LINUX; Javier García de Jaló, Ikern Aquinaga. 6. Manual Basico de Ubuntu
