con lo que es posible hacer las conversiones de UTM a latitud y longitud y viceversa.<\/p>\n\n\n\n
La f\u00f3rmula de distancia en UTM es la tradicional para distancia entre dos puntos en el plano.<\/p>\n\n\n
\n# Datos desde GPS Diferencial a un archivo.txt\n# Incorpora distancias a los vertices (baliza),\n# a\u00f1ade grupo por sector, tipo de punto y LOS.\n# Girni 2020-10-07 edelros@espol.edu.ec\nimport numpy as np\nimport pandas as pd\nimport utm\n\n# INGRESO\n# Archivo de coordenadas\narch_gps = 'ubicaPuntos01.txt'\n\n# Archivo salida procesado con distancias\narch_gpsrsm = 'resumen_ubica03.txt'\nzona ='17 M'\n\n# Referencias\nbaliza = {'d1':'gtwRECT',\n 'd2':'gtwFIEC',\n 'd3':'gtwFCNM'}\ngrupo = ['FIEC','FCNM','RECT','CIRC']\ntipo = ['punto','1m','gtw','dispositivo']\n\n# digitos decimales en distancias\ndigitos = 3\n\n# PROCEDIMIENTO\n# leer coordenadas\nubica = pd.read_csv(arch_gps,index_col='etiqueta')\nubica.drop('item',inplace=True, axis=1)\nubica = pd.DataFrame(ubica)\nn = len(ubica)\nzonanum = int(zona[0:2])\nzonalet = zona[3]\n\n# vertices con balizas\nbaliza_key = list(baliza.keys())\nbaliza_val = list(baliza.values())\nvertices = ubica.loc[baliza_val]\n\n# distancias a vertices\nfor fila in vertices.index:\n x1 = vertices['c_este'][fila]\n y1 = vertices['c_norte'][fila]\n x2 = ubica['c_este']\n y2 = ubica['c_norte']\n dist = np.sqrt((x2-x1)**2 + (y2-y1)**2)\n dist = np.round(dist,digitos)\n donde = baliza_val.index(fila)\n cual = baliza_key[donde]\n etiq_col = 'dist_'+cual\n ubica[etiq_col] = dist\n\n# redondear a dos digitos\nubica['c_este'] = np.round(ubica['c_este'],digitos)\nubica['c_norte'] = np.round(ubica['c_norte'],digitos)\nubica['altitud'] = np.round(ubica['altitud'],digitos)\n\n# a\u00f1adir coordenadas en latitud y longitud\nubica['longitud'] = 0.0\nubica['latitud'] = 0.0\nfor fila in ubica.index:\n lon_utm = ubica['c_este'][fila]\n lat_utm = ubica['c_norte'][fila]\n coord_gra = utm.to_latlon(lon_utm,lat_utm,\n zonanum,zonalet)\n ubica['latitud'][fila] = coord_gra[0]\n ubica['longitud'][fila] = coord_gra[1]\n\n# grupo y tipo en cada punto\nubica.insert(0,'grupo','')\nubica.insert(1,'tipo','')\nfor cadauno in ubica.index:\n \n # etiqueta de grupo\n esgrupo = 'CIRC'\n for ungrupo in grupo:\n cond1 = cadauno.startswith(ungrupo)\n cond2 = cadauno.endswith(ungrupo)\n if cond1 or cond2:\n esgrupo = ungrupo\n \n # etiqueta de tipo\n estipo = tipo[0]\n for untipo in tipo:\n cond1 = cadauno.startswith(untipo)\n if cond1:\n estipo = untipo\n if esgrupo == 'CIRC':\n estipo = 'dispositivo'\n \n ubica.loc[cadauno,'grupo'] = esgrupo\n ubica.loc[cadauno,'tipo'] = estipo\n\n# generar columna de Linea de vista LOS\ncolumna = 2\nfor cadauno in baliza.keys():\n ubica.insert(columna,'LOS_'+cadauno,1)\n columna = columna + 1\n\n# SALIDA\nprint('archivo resumen:', arch_gpsrsm)\nprint('una muestra de archivo: ')\nprint(ubica.head())\nprint('vertices: \\n',vertices)\nubica.to_csv(arch_gpsrsm)\n<\/pre><\/div>","protected":false},"excerpt":{"rendered":"Procesa las coordenadas de cada punto medido, registradas en un archivo de texto usando un GPS Diferencial. Las coordenadas se encuentran en formato UTM ubicadas en la zona \"17 M\". Ejemplo de archivo de coordenadas del GPS Diferencial A los datos de cada punto se a\u00f1aden los c\u00e1lculos de distancia hacia cada \"baliza\". En el […]<\/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-2724","post","type-post","status-publish","format-standard","hentry","category-lorawan-localizacion-rssi","tag-python"],"_links":{"self":[{"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts\/2724","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=2724"}],"version-history":[{"count":14,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts\/2724\/revisions"}],"predecessor-version":[{"id":4836,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/posts\/2724\/revisions\/4836"}],"wp:attachment":[{"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/media?parent=2724"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/categories?post=2724"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/girni\/wp-json\/wp\/v2\/tags?post=2724"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}
![\"GPS](\"http:\/\/blog.espol.edu.ec\/girni\/files\/2020\/11\/GPSDiferencial01.png\")
<\/a><\/figure>\n\n\n\n