[Talk-es] Coordenadas UTM extrañérrimas
Hola Mirando por el portal de opendata de bcn me he encontrado un listado de todas las fuentes, tanto las de beber como las ornamentales, de la ciudad de Barcelona. Están con licencia CC BY 3.0 que por lo que he visto es compatible por lo que he pensando en importarlas, pero ahora es cuando me encuentro el problema. Según ellos el listado de fuentes usa coordenadas UTM sistema municipal He estado probando varias cosas para ver si conseguía transformar esto a algo manejable pero todos los conversores que he probado me dicen que las coordenadas están fuera de rango. Ejemplo: X-Y : 29621467,82587692 debería estar más o menos por aqui (41,3898301, 2,1571738) http://www.openstreetmap.org/node/2964058469 Si alguien me echa un cable con algo que pueda scriptear le estaré altamente agradecido Un saludo -- Celso González @PerroVerd ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
Yo en su momento me escribí este script en python para arreglar las
coordenadas en UTM y pasarlo a osm, espero que te pueda servir de ayuda.
Hola
Mirando por el portal de opendata de bcn me he encontrado un listado de
todas las fuentes, tanto las de beber como las ornamentales, de la
ciudad de Barcelona.
Están con licencia CC BY 3.0 que por lo que he visto es compatible por lo
que he pensando en importarlas, pero ahora es cuando me encuentro el
problema.
Según ellos el listado de fuentes usa coordenadas UTM sistema municipal
He estado probando varias cosas para ver si conseguía transformar esto a
algo manejable pero todos los conversores que he probado me dicen que
las coordenadas están fuera de rango.
Ejemplo:
X-Y : 29621467,82587692
debería estar más o menos por aqui
(41,3898301, 2,1571738) http://www.openstreetmap.org/node/2964058469
Si alguien me echa un cable con algo que pueda scriptear le estaré
altamente agradecido
Un saludo
--
Cruz Enrique Borges Hernández
Email: cruz.bor...@deusto.es
DeustoTech Energy
Telefono: 944139000 ext.2052
Avda. Universidades, 24
48007 Bilbao, Spain#! /usr/bin/python
# -*- coding: iso-8859-15 -*-
# Arregla los datos de los suministros y genera un osm con su geolocalizacion.
import csv
from math import pi, sin, cos, tan, sqrt, radians, degrees
from pyproj import Proj, transform
wgs84 = Proj(proj='latlong',datum='WGS84')
utm30 = Proj(proj='utm',zone='30',init='epsg:23030',nadgrids=/home/cruzki/src/cat2osm/peninsula.gsb)
#utm30 = Proj(proj='utm',zone='30',init='epsg:25830')
fields = [OBJECTID,IDEST,ESTACION,DIRECC,POBLACION,LINK_ESTAC,ZONA,LINK_ZONA,PROVCOD,NO2,NO_,NOX,O3,CO,SO2,PST,PM10,PM2_5,BEN,ETB,MPX,O_X,TOL,NH3,SH2,HC,MSV,CL2,HCL,HCMET,HCNOMET,TEMPERATUR,HUMEDAD,PRESION,RADIACION,RAD_UV,LLUVIA,LL_ACIDA,DIR_VIENTO,V_VIENTO,XUTM,YUTM]
f = csv.DictReader(open('estaciones.csv'),delimiter=';',quotechar='')
g = open('Bizkaia-corregidos.txt','w')
h = open('localizacion.osm','w')
h.write(?xml version='1.0' encoding='UTF-8'?\n)
h.write(osm version='0.6' generator='JOSM'\n)
for row in f:
# Transformamos de utm a latlong
lon, lat, z = transform(utm30,wgs84,float(row[XUTM]),float(row[YUTM]),0)
# Transformamos de latlong a utm para corregir los datos de los clientes
easting, northing, z = transform(wgs84,utm30,lon,lat,0)
# Guardamos OSM
h.write(node id=\ + row[OBJECTID] + \ version=\6\ timestamp=\2011-10-27 08:26:29.565\ uid=\292702\ user=\AnderPijoan\ changeset=\5407370\ lat=\+str(lat)+\ lon=\+str(lon)+\\n/node\n)
# Guardamos clientes
g.write(row[OBJECTID]+,+str(lat)+,+str(lon)+\n)
h.write(/osm\n)
g.close()
h.close()___
Talk-es mailing list
Talk-es@openstreetmap.org
https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
On Mon, Sep 08, 2014 at 02:42:58AM -0700, Cruz Enrique Borges wrote: Yo en su momento me escribí este script en python para arreglar las coordenadas en UTM y pasarlo a osm, espero que te pueda servir de ayuda. Mi script viene a ser igual pero cambio utm30 por utm31 utm31 = Proj(proj='utm',zone='31',init='epsg:23031',nadgrids=/home/celso/openstreetmap/peninsula.gsb) y los datos me vienen como enteros por los que lo divido por 100 lon, lat, z = transform(utm31,wgs84,float(row[INVENTARI_COORDENADA_X])/100,float(row[INVENTARI_COORDENADA_Y])/100,0) Por lo demás es lo mismo, Lo bueno El mapa de puntos que me da tiene la forma d e la ciudad Lo malo Me la situa en África Un saludo y gracias ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
Sin mirar nada del resultado, simplemente al oir lo de que te lo sitúa en África. ¿Es por el este de África? Esto nos ha pasado alguna vez que nos colamos y metimos las coordenadas del revés. Hay librerías q primero reciben LAT y luego LON y otras que al revés. El 8 de septiembre de 2014, 12:39, Celso González ce...@mitago.net escribió: On Mon, Sep 08, 2014 at 02:42:58AM -0700, Cruz Enrique Borges wrote: Yo en su momento me escribí este script en python para arreglar las coordenadas en UTM y pasarlo a osm, espero que te pueda servir de ayuda. Mi script viene a ser igual pero cambio utm30 por utm31 utm31 = Proj(proj='utm',zone='31',init='epsg:23031',nadgrids=/home/celso/openstreetmap/peninsula.gsb) y los datos me vienen como enteros por los que lo divido por 100 lon, lat, z = transform(utm31,wgs84,float(row[INVENTARI_COORDENADA_X])/100,float(row[INVENTARI_COORDENADA_Y])/100,0) Por lo demás es lo mismo, Lo bueno El mapa de puntos que me da tiene la forma d e la ciudad Lo malo Me la situa en África Un saludo y gracias ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es -- Ander Pijoan Lamas Research Assistant, Deustotech Computer Science Engineer University of Deusto E-mail: ander.pij...@deusto.es Phone: +34 664471228 in: http://www.linkedin.com/profile/view?id=162888312 ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
Y lo de dividir por 100, por qué lo haces? Me extraña mucho que lo tengas que hacer. Por otro lado, creo que al usar floats alguna vez tuvimos problemas. On Mon, Sep 08, 2014 at 02:42:58AM -0700, Cruz Enrique Borges wrote: Yo en su momento me escribí este script en python para arreglar las coordenadas en UTM y pasarlo a osm, espero que te pueda servir de ayuda. Mi script viene a ser igual pero cambio utm30 por utm31 utm31 = Proj(proj='utm',zone='31',init='epsg:23031',nadgrids=/home/celso/openst reetmap/peninsula.gsb) y los datos me vienen como enteros por los que lo divido por 100 lon, lat, z = transform(utm31,wgs84,float(row[INVENTARI_COORDENADA_X])/100,float(row [INVENTARI_COORDENADA_Y])/100,0) Por lo demás es lo mismo, Lo bueno El mapa de puntos que me da tiene la forma d e la ciudad Lo malo Me la situa en África Un saludo y gracias ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es -- Cruz Enrique Borges Hernández Email: cruz.bor...@deusto.es DeustoTech Energy Telefono: 944139000 ext.2052 Avda. Universidades, 24 48007 Bilbao, Spain ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
On Mon, Sep 08, 2014 at 12:43:09PM +0200, Ander Pijoan wrote: Sin mirar nada del resultado, simplemente al oir lo de que te lo sitúa en África. ¿Es por el este de África? Esto nos ha pasado alguna vez que nos colamos y metimos las coordenadas del revés. Hay librerías q primero reciben LAT y luego LON y otras que al revés. Sí, ya me ha pasado alguna vez, probé a dar la vuelta a las coords pero tampoco es este caso ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
On Mon, Sep 08, 2014 at 12:54:15PM +0200, Cruz Enrique Borges wrote: Y lo de dividir por 100, por qué lo haces? Me extraña mucho que lo tengas que hacer. Por otro lado, creo que al usar floats alguna vez tuvimos problemas. Realmente lo de dividir por 100, por 1000 etc lo he ido probando a ojimetro ya que los enteros que proporcionan estarían fuera de rango -- Celso González @PerroVerd ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas
Celso, Yo te diría que a la X le falta un 4 y a la Y un 45 y la división por 1000. Que en goolzoom caería por donde tu quieres 41.3898957869,2.15707415667 Aunque en osm no, pero eso es un problema de proyecciones que aunque a mi se me escapa ahora mismo, seguro que se puede solucionar. On 8Sep, 2014, at 11:37 AM, Celso González ce...@mitago.net wrote: X-Y : 29621467,82587692 429621467, 4582587692 ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas (Celso González)
Hola Celso, Después de investigar un poco parece que la formula para obtener las coordenadas es la siguiente: Pasar coordenadas en UTM sistema municipal a ED50 31N (EPSG:23031) En ambas coordenadas los 3 últimos digitos son los decimales En las coordenadas X hay que agregar un 4 al inicio En las coordenadas Y hay que agregar un 45 al inicio. Ejemplo: X,Y BCN 32424958, 82336236 == X,Y ED50 432424.958, 4582336.236 En InstaMaps puedes ver un mapa de las fuentes de BCN (te puedes descargar las fuentes en varios formatos y en WGS84 para que las puedas agregar en OSM) http://www.instamaps.cat/geocatweb/visor.html?businessid=9b0b91ff6ee6a5b5969a4dffbbec754fid=47949 -- Saludos, Bolo www.geoinquiets.cat ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
Re: [Talk-es] Coordenadas UTM extrañérrimas (Celso González)
On Mon, Sep 08, 2014 at 02:58:15PM +0200, Wladimir Szczerban wrote: Hola Celso, Después de investigar un poco parece que la formula para obtener las coordenadas es la siguiente: Pasar coordenadas en UTM sistema municipal a ED50 31N (EPSG:23031) En ambas coordenadas los 3 últimos digitos son los decimales En las coordenadas X hay que agregar un 4 al inicio En las coordenadas Y hay que agregar un 45 al inicio. Ejemplo: X,Y BCN 32424958, 82336236 == X,Y ED50 432424.958, 4582336.236 En InstaMaps puedes ver un mapa de las fuentes de BCN (te puedes descargar las fuentes en varios formatos y en WGS84 para que las puedas agregar en OSM) http://www.instamaps.cat/geocatweb/visor.html?businessid=9b0b91ff6ee6a5b5969a4dffbbec754fid=47949 Efectivamente, con añadiendo el 4 y el 45 y haciendo la conversión a ED50 31N (EPSG:23031) sale tal como era de esperar, pero... Veo que la mayoría de fuentes me caen dentro de edificios o en medio de calles, hay una ligera desviación de pocos metros pero no se si es debido al sistema de conversión o si este se puede apurar un poco más o bien que los datos originales no tienen la precisión correcta. Gracias a los geoinquiets por la ayuda -- Celso ___ Talk-es mailing list Talk-es@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-es
