{"id":2730,"date":"2020-11-21T16:31:39","date_gmt":"2020-11-21T21:31:39","guid":{"rendered":"http:\/\/blog.espol.edu.ec\/girni\/?p=2730"},"modified":"2022-02-25T21:21:07","modified_gmt":"2022-02-26T02:21:07","slug":"procesa-datos-ubica-los-puntos-en-grafica-2d-o-3d","status":"publish","type":"post","link":"https:\/\/blog.espol.edu.ec\/girni\/procesa-datos-ubica-los-puntos-en-grafica-2d-o-3d\/","title":{"rendered":"3.4 Procesa datos. Ubica los puntos en gr\u00e1fica 2D o 3D"},"content":{"rendered":"

Para revisar los datos de coordenadas y ubicaciones obtenidos en el proceso anterior, se grafica la ubicaci\u00f3n de los puntos usando las coordenadas este, norte y altura registradas con el GPS diferencial.<\/p>\n

Las gr\u00e1ficas se pueden realizar en 2D y en 3D, se adjunta los resultados:<\/p>\n

\"\"<\/a>Como referencia se usan las etiquetas de los gateways que permiten observar la posici\u00f3n relativa de cada punto.<\/p>\n

Los grupos de puntos a mostrar se pueden seleccionar en el bloque de ingreso dentro de los par\u00e1metros de la gr\u00e1fica, permitiendo obervar con mayor detalle la ubicaci\u00f3n de cada punto<\/p>\n

Una observaci\u00f3n en 3D de los puntos permite revisar que lso valores de altura presentan menor precisi\u00f3n que la de posici\u00f3n, principalmente en puntos ubicados entre vegetaci\u00f3n.<\/p>\n

\"\"<\/a><\/p>\n

Se observa que las medidas de alturas tienen valores muy variables, por lo que el uso de los valores de altura se descarta inicialmente de los modelos. Se podr\u00eda realizar posteriormente un an\u00e1lisis m\u00e1s detallado de lo presentado con los valores de las alturas.<\/p>\n

archivo de resumen<\/strong><\/em> de coordenadas y distancias : resumen_ubica01.txt<\/a><\/p>\n

\n

Algoritmo en Python<\/h2>\n

En el caso de gr\u00e1ficas 3D, para crear el archivo.gif animado es necesario instalar imagemagic. El nombre del archivo animado es 'rotando3D.gif' que se guarda en el mismo directorio del algoritmo.py.<\/p>\n

# Datos desde GPS Diferencial archivo.txt<\/span>\r\n# Graficas de coordenadas 2D o 3D<\/span>\r\n# Girni 2020-10-07 edelros@espol.edu.ec<\/span>\r\nimport<\/span> numpy as<\/span> np\r\nimport<\/span> pandas as<\/span> pd\r\nimport<\/span> matplotlib.pyplot as<\/span> plt\r\nfrom<\/span> mpl_toolkits.mplot3d import<\/span> Axes3D\r\nimport<\/span> matplotlib.animation as<\/span> animation\r\n\r\n# INGRESO<\/span>\r\n# Archivo procesado con distancias<\/span>\r\narch_gpsrsm = 'resumen_ubica01.txt'<\/span>\r\n\r\n# Parametros de gr\u00e1fica<\/span>\r\nmostrargrp = ['FIEC'<\/span>,'FCNM'<\/span>,'RECT'<\/span>]\r\n# ['FIEC','FCNM','RECT','CIRC']<\/span>\r\n\r\nmostrartip = ['punto'<\/span>,'1m'<\/span>,'dispositivo'<\/span>,'gtw'<\/span>]\r\n# ['punto','1m','dispositivo','gtw']<\/span>\r\n\r\ntipograf = '3D'<\/span> # '2D','3D'<\/span>\r\narch_rotacion = 'rotando3D.gif'<\/span>\r\n\r\n# Referencias<\/span>\r\nbaliza  = {'d1'<\/span>:'gtwRECT'<\/span>,\r\n           'd2'<\/span>:'gtwFIEC'<\/span>,\r\n           'd3'<\/span>:'gtwFCNM'<\/span>}\r\ngrupo   = ['FIEC'<\/span> ,'FCNM'<\/span>  ,'RECT'<\/span>,'CIRC'<\/span>]\r\ncolores = ['green'<\/span>,'orange'<\/span>,'grey'<\/span>,'magenta'<\/span>]\r\ntipo    = ['punto'<\/span>,'1m'<\/span> ,'gtw'<\/span>,'dispositivo'<\/span>]\r\nmarcas  = [    'o'<\/span>,'D'<\/span>  ,'D'<\/span>  ,'*'<\/span> ]\r\n\r\n# PROCEDIMIENTO<\/span>\r\n# leer coordenadas<\/span>\r\nubica = pd.read_csv(arch_gpsrsm, index_col='etiqueta'<\/span>)\r\nubica = pd.DataFrame(ubica)\r\nn = len<\/span>(ubica)\r\n\r\n# puntos de vertices<\/span>\r\nvertices = ubica.loc[baliza.values()]\r\n\r\n# segmento de grupo\/tipo<\/span>\r\nubica['usar'<\/span>]  = False<\/span>\r\nubica['color'<\/span>] = 'yellow'<\/span>\r\nubica['marca'<\/span>] = 'o'<\/span>\r\nfor<\/span> fila in<\/span> ubica.index:\r\n    unaeti  = str<\/span>(fila)\r\n    ungrupo = ubica['grupo'<\/span>][fila]\r\n    untipo  = ubica['tipo'<\/span>][fila]\r\n    cond1  = ungrupo in<\/span> mostrargrp\r\n    cond2  = untipo in<\/span> mostrartip\r\n    if<\/span> cond1:\r\n        cual = grupo.index(ungrupo)\r\n        ubica.loc[fila,'color'<\/span>] = colores[cual]\r\n    if<\/span> cond2:\r\n        cual = tipo.index(untipo)\r\n        ubica.loc[fila,'marca'<\/span>] = marcas[cual]\r\n    if<\/span> (cond1 and<\/span> cond2):\r\n        ubica.loc[fila,'usar'<\/span>]  = True<\/span>\r\n\r\n# SALIDA<\/span>\r\nprint<\/span>(ubica.head())\r\n\r\n# Grafica 2D --------------------------<\/span>\r\nif<\/span> tipograf == '2D'<\/span>:\r\n    figura, grafica = plt.subplots()\r\n    for<\/span> fila in<\/span> ubica.index:\r\n        unalon = ubica['c_este'<\/span>][fila]\r\n        unalat = ubica['c_norte'<\/span>][fila]\r\n        unaalt = ubica['altitud'<\/span>][fila]\r\n        usar   = ubica['usar'<\/span>][fila]\r\n        uncolor = ubica['color'<\/span>][fila] \r\n        unamarca = ubica['marca'<\/span>][fila]\r\n        untipo = ubica['tipo'<\/span>][fila]\r\n        if<\/span> usar:\r\n            grafica.scatter(unalon,unalat,\r\n                            color = uncolor,\r\n                            marker = unamarca,\r\n                            label = fila)\r\n        if<\/span> usar and<\/span> untipo=='gtw'<\/span>:\r\n            grafica.annotate(fila,\r\n                             (unalon,unalat))\r\n    plt.xlabel('UTM_este'<\/span>)\r\n    plt.ylabel('UTM_norte'<\/span>)\r\n    plt.title('Ubicacion UTM'<\/span>)\r\n    plt.show()\r\n\r\n# Grafica 3D --------------------------<\/span>\r\nif<\/span> tipograf == '3D'<\/span>:\r\n    figura = plt.figure()\r\n    grafica = Axes3D(figura)\r\n    for<\/span> fila in<\/span> ubica.index:\r\n        unalon  = ubica['c_este'<\/span>][fila]\r\n        unalat  = ubica['c_norte'<\/span>][fila]\r\n        unaalt  = ubica['altitud'<\/span>][fila]\r\n        usar    = ubica['usar'<\/span>][fila]\r\n        uncolor = ubica['color'<\/span>][fila] \r\n        unamarca = ubica['marca'<\/span>][fila]\r\n        untipo  = ubica['tipo'<\/span>][fila]\r\n        if<\/span> usar:\r\n            grafica.scatter(unalon,unalat,\r\n                            unaalt,\r\n                            marker = unamarca,\r\n                            color = uncolor,\r\n                            label = fila)\r\n        if<\/span> (usar and<\/span> (untipo=='gtw'<\/span>)):\r\n            grafica.text(unalon,unalat,\r\n                         unaalt,fila)\r\n    grafica.set_xlabel('UTM_este'<\/span>)\r\n    grafica.set_ylabel('UTM_norte'<\/span>)\r\n    grafica.set_zlabel('altitud'<\/span>)\r\n    grafica.set_title('Ubicacion UTM'<\/span>)\r\n\r\n    def<\/span> rotate<\/span>(angle):\r\n        grafica.view_init(azim=angle)\r\n\r\n    print<\/span>(\"realizando animation\"<\/span>)\r\n    rot_animation = animation.FuncAnimation(figura,\r\n                    rotate,\r\n                    frames = np.arange(45,360+45,10),\r\n                    interval = 200)\r\n    rot_animation.save(arch_rotacion, dpi=80,\r\n                    writer = animation.PillowWriter(fps=5))\r\n    plt.show()\r\n<\/pre>\n
Referencia<\/strong>: Graficas 3D puntos dispersos-scatter<\/a> en m\u00e9todos num\u00e9ricos<\/p>\n","protected":false},"excerpt":{"rendered":"
Para revisar los datos de coordenadas y ubicaciones obtenidos en el proceso anterior, se grafica la ubicaci\u00f3n de los puntos usando las coordenadas este, norte y altura registradas con el GPS diferencial. Las gr\u00e1ficas se pueden realizar en 2D y en 3D, se adjunta los resultados: Como referencia se usan las etiquetas de los gateways … <\/p>\n
Continuar leyendo \"3.4 Procesa datos. Ubica los puntos en gr\u00e1fica 2D o 3D\"<\/span><\/a><\/p>\n","protected":false},"author":8043,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1463489],"tags":[20777],"class_list":["post-2730","post","type-post","status-publish","format-standard","hentry","category-localizacion-rssi-lora","tag-python"],"_links":{"self":[{"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts\/2730","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/users\/8043"}],"replies":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/comments?post=2730"}],"version-history":[{"count":15,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts\/2730\/revisions"}],"predecessor-version":[{"id":3974,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts\/2730\/revisions\/3974"}],"wp:attachment":[{"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/media?parent=2730"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/categories?post=2730"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/tags?post=2730"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}