Przybliżanie traku gpx do najbliższej drogi

kamilek · 2012-08-16T19:52:42.000Z

Gotowy skrypt do pobierania danych z osrm (zapytania na bazie tracku gpx), i przetwarzania odpowiedzi do gotowego tracku gpx

<?php $xml = simplexml_load_file('track.gpx'); file_put_contents('tmp.txt', ''); $tmp = fopen("tmp.txt","a"); foreach($xml -> trk -> trkseg -> trkpt as $trkpt) { $atr = $trkpt->attributes(); $adr = "[http://0.0.0.0:5000/nearest?loc=](http://0.0.0.0:5000/nearest?loc=)".$atr['lat'].",".$atr['lon']; $odp_serv = file_get_contents($adr); fputs($tmp,$odp_serv."\n"); } fclose($tmp); $xml_gpx = simplexml_load_file('track_approx.gpx'); foreach($xml -> trk -> trkseg -> trkpt as $trkpt) { unset($xml_gpx->trk->trkseg->trkpt); } $tablica = file("tmp.txt"); for ($i=0; $itrk->trkseg->addChild('trkpt'); $trkpt->addAttribute('lat',$lat); $trkpt->addAttribute('lon',$lon); } $str = $xml_gpx -> asXML(); $q = fopen('track_approx.gpx', 'w+'); fwrite($q, $str); fclose($q); ?>

kamilek · 2012-08-20T13:29:42.000Z

Może się przyda komuś, dzięki jeszcze raz za wcześniejszą pomoc

kamilek · 2012-08-21T10:10:37.000Z

serekmedia:

Ja to robię w Trafficu:biorę punkt drogi z gps i szukam najbliższych dróg w sensownym bbox porównuję azymut śladu z azymutem drogi - jak jest w odpowiednim zakresie stopni to trafiłem w drogę i mam też od razu w znaku kierunek względem drogi - czyli na dobrą sprawę w całej kontynentalnej EU stronę jezdni (należy w to jeszcze wmieszać oneway, bo tam pod prąd minus nie może wyjść).

Tak, bez bazy to robota jest w ogóle bez sensu. Jeśli to plik osm.xml to on zawiera węzły, a z nich dopiero powstają drogi, drogi są ciągami węzłów - Ciebie interesuje odcinek z tego. Postgis ma do tego odpowiednie narzędzia - bez tego będziesz wyważał otwarte drzwi - chyba że tego dotyczy Twoja praca - więc zapytaj czy masz do tego użyć postgisa czy masz właśnie opracować to na gołym pliku - co też jest możliwe, ale wymagałoby 2 przebiegów przez plik - traktując go jako strumień - wybrać bboxowe nody (węzły) w zakresie traku GPX i wybrać linie lub lepiej odcinki składające się z tych nodów - dostaniesz 2 tablice już użytecznych informacji.Jednak sądzę że ktoś były niepoważny dając Ci robotę tego typu czyli zrobić nie używając bazy.

Muszę jednak aprosymować ślad bez osrm, choć już zrobiłem to za pomocą tego serwera. Mógłbyś przybliżyć mi trochę dokładniej metodę którą zaproponowałeś wcześniej? Co to jest właściewie ten bbox i traffic?

serekmedia · 2012-08-21T11:22:43.000Z

Traffic to jest warstwa mapy ze średnimi prędkościami, po nałożeniu na inną mapę widać gdzie nie ma płynności jazdy itp.
BBOX jest to po prostu “kwadrat” ograniczony 4 współrzędnymi. Nie jest to kwadrat geometryczny bo południki się zwężają ku biegunom, ale na nich łatwo (i szybko) się ogranicza obszar poszukiwań. Np jeśli szukasz drogi w odległości 1km od punktu to musiałbyś sprawdzać każdą drogę czy nie leży 1km od zadanego punktu. Jeżeli do zadanych współrzędnych dodasz sobie ±0.1 stopnia, to uzyskasz “kwadrat” który szybko wybierze te drogi które przecinają bbox i tylko im sprawdzę dokładny 1km.

Miszka · 2012-08-21T11:56:39.000Z

kamilek:

Muszę jednak aprosymować ślad bez osrm, choć już zrobiłem to za pomocą tego serwera. Mógłbyś przybliżyć mi trochę dokładniej metodę którą zaproponowałeś wcześniej? Co to jest właściewie ten bbox i traffic?

Napisz co było źle z OSRMem. Podpowiemy Ci coś znowu.
Podstawowe pytanie to czy możesz używać bazy danych.

kamilek · 2012-08-21T12:18:56.000Z

Hmm co było źle z OSRMem? Wg mnie było dobrze i to działało, przybliżało track do najbliższej drogi i go wyświetlało ale mój opiekun stwierdził, że to ma działać bez tego serwera. O bazie danych nic nie mówił, twierdził ze mam to zrobić za pomocą opellayers i operować na pliku mapowym (moge stworzyć sobie nawet swój własny z jedną trasą) - open layers ma chyba wbudowane funkcje do obsługi plików xml. O ile zadeklarowanie tego bboxa rozumiem i wydaje się dosyć logiczne ograniczenie obszaru poszukiwań do minimum, to nie rozumiem założeń tego traffica, może mogę prosić o jakiś przykład tego rozwiązania. Jeszcze jedno małe pytanie, czy pasy ruchu w normalnej drodze są rozgraniczone w OSM?

serekmedia · 2012-08-21T13:04:43.000Z

Traffic masz tu: http://osmapa.pl/beta/
Włącz sobie warstwę Traffic (korki)
Ja mam do dyspozycji (podobnie jak Ty) ślad po samochodzie, przy czym on ma rozrzut przez niedokładność gps, rzadką rejestrację śladu itp. Aby wyrysować traffic po drogach muszę (tak jak Ty) znaleźć drogę (i później w niej odpowiedni odcinek) najbardziej pasujący do kursu pojazdu.
Więc biorę punkt ze śladu GPS, robię wokoło niego bbox i szukam dróg przecinających ten bbox (w postgis robi to odpowiednia funkcja np:
select way from planet_roads where ST_Intersects(way,ST_GeomFromText(‘POLYGON((14 49, 24 49, 24 55, 14 55, 14 49))’,4326))
wyjmie wszystkie drogi w zgrubnym BBOX Polski. Oczywiście szukając dróg do traffica używam malutkiego bboxa który obejmuje mniej więcej kwadrat 50x50m. Takie użycie spowoduje że zadziała index gis postgresa - ST_distance nie użyje indexu.
Następnym warunkiem w where jest ST_Distance które wymaga już policzenia dystansu i on jest też wspomagany przez ST_ClosestPoint szukający najbliższego punktu z tej drogi i szukam wtedy podejrzanych odcinków dróg po których mogę potencjalnie jechać, sortuję je po odległości. Jak mam wycięte podejrzane odcinki to już aplikacyjnie szukam tego który ma odpowiedni azymut do odczytu z GPS i ten odcinek notuję - sumując prędkość i inkremetując ilość pomiarów.
Ta baza pozwala mi wyrysować wyżej widoczne kafelki - tu tez się bbox przedaje - tym razem bbox kafelka wyliczony z XYZ kafelka.

kamilek · 2012-08-22T11:55:12.000Z

Troche mnie już to przerasta ale coś zaczełem robić:
Zainstalowałem PostreQSL i PostGis wg instrukcji http://wiki.openstreetmap.org/wiki/Mapnik/PostGIS

for Ubuntu >=12.04:

sudo apt-get install postgresql-9.1-postgis postgresql-contrib-9.1
sudo -u postgres -i -H
createuser -SdR gisuser
createdb -E UTF8 -O gisuser gis
createlang plpgsql gis
psql -d gis -f /usr/share/postgresql/9.1/contrib/postgis-1.5/postgis.sql

to populate the table spatial_ref_sys (mandatory for use with osm2pgsql):

psql -d gis -f /usr/share/postgresql/9.1/contrib/postgis-1.5/spatial_ref_sys.sql

common for all three

psql gis -c “ALTER TABLE geometry_columns OWNER TO gisuser”
psql gis -c “ALTER TABLE spatial_ref_sys OWNER TO gisuser”
exit

Rozumiem że teraz muszę wysłać odpowiednie zapytanie:
coś w stylu:
psql gis “select way from planet_roads where ST_Intersects(way,ST_GeomFromText(‘POLYGON((14 49, 24 49, 24 55, 14 55, 14 49))’,4326))” ??

czy planet_roads ma być ścieżką dostępu do pliku OSMowego? Ok, sprawdziłem te punkty biorą w kwadrat całą powierzchnię Polski. A czym jest liczba 4326 na końcu zapytania?
Co właściwie zwróci mi te zapytanie?

Dotevo · 2012-08-22T12:02:32.000Z

Najpierw musisz zaimportować dane z OSM do bazy np. poczytaj o osmosis lub osm2pgsql.

Ta liczba to odwzorowanie geograficzne (projekcja).

pbabik · 2012-08-22T12:33:08.000Z

kamilek:

Co właściwie zwróci mi te zapytanie?

Coś takiego:

                        way                         
----------------------------------------------------
 0101000020840800000BE37D338C7B1D41603AAF3641670F41
 01010000208408000086ACE512ABB91D4160AA3493B9BF0F41
 0101000020840800002641477D08AD1D41806E2A457F011041

....

są to zakodowane geometrie dróg, które przecinają (intersects) zadany obszar. Same z siebie na wiele się nie zdadzą, ale można je wsadzić w kolejne zapytanie.

kamilek · 2012-08-22T13:09:55.000Z

Za pomocą narzędzia osmosis wyciągam z poland.osm wszystko co jest wewnątrz granicy miasta (wewnątrz polygonu):
bzcat poland.osm.bz2 | ~/Pobieranie/osmosis-0.40.1/bin/osmosis --read-xml enableDateParsing=no file=- --bounding-polygon file=tmp/out.poly --write-xml file=tmp/out.osm
I tym sposobem w tmp/out.osm mamy całą zawartość miejscowości w formacie OSM ( http://wiki.openstreetmap.org/wiki/.osm )"

Czyli biorę pojedynczy punkt z tracka gpx, robię z niego bboxa załóżmy 50 metrów czyli około 0.004504 stopnia - bbox (lon lat, lon+0.004504 lat, lon+0.004504 lat+0.004504, lon lat+lat+0.004504 , lon lat) importuje te dane z ososis do bazy i potem dalej wybieram way z tego wyniku?

kamilek · 2012-08-23T13:07:46.000Z

witam,
udało mi się za pomocą osmosis’a, za pomocą polecenia
osmosis --read-xml planet.osm --bounding-box top=52.4058 left=20.895 bottom=52.405 right=20.896 --write-xml plik.xml
wyodrębnić jakąś tam drogę.

czy po zaimportowaniu pliku plenet.osm do bazy postgresql
osmosis --read-xml file=“planet.osm” --write-apidb host=“x” database=“x” user=“x” password=“x”
polecenie te zrobi to samo co to:
select way from planet_roads where ST_Intersects(way,ST_GeomFromText(‘POLYGON((14 49, 24 49, 24 55, 14 55, 14 49))’,4326)) ???

jako planet.osm jest u mnie tylko polska część osm, polecenie wykonuje się bardzo długo, jedno zapytanie i musze czekać ponad minute czasu na odpowiedz.

Jaki teraz powinien być mój algorytm wykonywanych działań?

Czy każdy punkt tracku powinienem przepuszczać przez te polecenie?

Bazę PostgreSQL i Postgisa tworzyłem na podstawie tej instrukcji http://www.earthisflat.net/zdp/ ćw 1 i 2

serekmedia · 2012-08-23T15:48:17.000Z

Osmosis zassa bazę inaczej niż osm2pgsql, ale dla Twoich potrzeb nie ma to aż tak dużego znaczenia. Zassaj PL do bazy.
przykład z tabelą planet_roads był po zassaniu przez osm2pgsql - tam jest osobna tabela na ulice.

W przypadku zassania osmosis metodą --write-pgsql drogi trafią do tabeli ways - ale wszystkie drogi, będziesz musiał wyjąc te których ways.id=way_tags.way_id i way_tags.k=‘highway’ i ewentualnie jeszcze typ way_tags.v = ‘motorway’, ‘primary’ itd. jakie tam potrzebujesz (budynek w osm też może mieć way, ale będzie miał k=‘building’ v=‘yes’. Co do obszaru bbox to będziesz musiał użyć ways.bbox - najszybciej - osmosis trzyma pudełko w bazie (chociaż nie wiem po co, łatwiej trzymać chyba index na st_envelope(linesting) i tak się do tego index na bbox sprowadzi, a użyteczność pudełka poza wyszukiwaniem jest IMO niewielka), a do dokładnego porównania ways.linestring.

Osmosis bardziej uniwersalnie zasysa dane, do takich operacji wygodniej niż osm2pgsql.

No co do traku to tak, powinieneś każdy punkt spróbować przyciągnąć do jakiejś drogi, zapewne kilka punktów się do 1 też może przyciągnąć.
Co gorsza, Ciebie interesuje FRAGMENT drogi, nie cała ale co i jak to już pisałem we wątku, myśl, bo nie zaliczysz praktyk

kamilek · 2012-08-24T07:57:21.000Z

Mam problem z zaimportowaniem planet.osm do postgisa, wywołuje w terminalu komendę:

osmosis --read-xml file=/home/kamil/osmosis-0.40.1/planet.osm --write-apidb host=“localhost” database=“osm” user=“kamil” password=“haslo”;

postgres@kamil-MS-1675:/home/kamil/osmosis-0.40.1/script$ osmosis --read-xml file=/home/kamil/osmosis-0.40.1/planet.osm --write-pgsimp host=“localhost” database=“osm” user=“kamil” password=“******”
Aug 24, 2012 10:51:45 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Osmosis Version 0.40.1
Aug 24, 2012 10:51:46 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Preparing pipeline.
Aug 24, 2012 10:51:46 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Launching pipeline execution.
Aug 24, 2012 10:51:46 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Pipeline executing, waiting for completion.
Aug 24, 2012 10:51:46 AM org.openstreetmap.osmosis.core.pipeline.common.ActiveTaskManager waitForCompletion
SEVERE: Thread for task 1-read-xml failed
org.openstreetmap.osmosis.core.OsmosisRuntimeException: Unable to read the schema version from the schema info table.
at org.openstreetmap.osmosis.pgsimple.common.SchemaVersionValidator.validateDBVersion(SchemaVersionValidator.java:90)
at org.openstreetmap.osmosis.pgsimple.common.SchemaVersionValidator.validateVersion(SchemaVersionValidator.java:50)
at org.openstreetmap.osmosis.pgsimple.v0_6.PostgreSqlWriter.initialize(PostgreSqlWriter.java:185)
at org.openstreetmap.osmosis.pgsimple.v0_6.PostgreSqlWriter.process(PostgreSqlWriter.java:779)
at org.openstreetmap.osmosis.xml.v0_6.impl.NodeElementProcessor.end(NodeElementProcessor.java:118)
at org.openstreetmap.osmosis.xml.v0_6.impl.OsmHandler.endElement(OsmHandler.java:107)
at org.apache.xerces.parsers.AbstractSAXParser.endElement(Unknown Source)
at org.apache.xerces.impl.XMLDocumentFragmentScannerImpl.scanEndElement(Unknown Source)
at org.apache.xerces.impl.XMLDocumentFragmentScannerImpl$FragmentContentDispatcher.dispatch(Unknown Source)
at org.apache.xerces.impl.XMLDocumentFragmentScannerImpl.scanDocument(Unknown Source)
at org.apache.xerces.parsers.XML11Configuration.parse(Unknown Source)
at org.apache.xerces.parsers.XML11Configuration.parse(Unknown Source)
at org.apache.xerces.parsers.XMLParser.parse(Unknown Source)
at org.apache.xerces.parsers.AbstractSAXParser.parse(Unknown Source)
at org.apache.xerces.jaxp.SAXParserImpl$JAXPSAXParser.parse(Unknown Source)
at org.apache.xerces.jaxp.SAXParserImpl.parse(Unknown Source)
at javax.xml.parsers.SAXParser.parse(SAXParser.java:195)
at org.openstreetmap.osmosis.xml.v0_6.XmlReader.run(XmlReader.java:108)
at java.lang.Thread.run(Thread.java:679)
Caused by: org.postgresql.util.PSQLException: ERROR: permission denied for relation schema_info
at org.postgresql.core.v3.QueryExecutorImpl.receiveErrorResponse(QueryExecutorImpl.java:2102)
at org.postgresql.core.v3.QueryExecutorImpl.processResults(QueryExecutorImpl.java:1835)
at org.postgresql.core.v3.QueryExecutorImpl.execute(QueryExecutorImpl.java:257)
at org.postgresql.jdbc2.AbstractJdbc2Statement.execute(AbstractJdbc2Statement.java:500)
at org.postgresql.jdbc2.AbstractJdbc2Statement.executeWithFlags(AbstractJdbc2Statement.java:374)
at org.postgresql.jdbc2.AbstractJdbc2Statement.executeQuery(AbstractJdbc2Statement.java:254)
at org.openstreetmap.osmosis.pgsimple.common.SchemaVersionValidator.validateDBVersion(SchemaVersionValidator.java:71)
… 18 more
Aug 24, 2012 10:51:46 AM org.openstreetmap.osmosis.core.Osmosis main
SEVERE: Execution aborted.
org.openstreetmap.osmosis.core.OsmosisRuntimeException: One or more tasks failed.
at org.openstreetmap.osmosis.core.pipeline.common.Pipeline.waitForCompletion(Pipeline.java:146)
at org.openstreetmap.osmosis.core.Osmosis.run(Osmosis.java:92)
at org.openstreetmap.osmosis.core.Osmosis.main(Osmosis.java:37)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:57)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
at java.lang.reflect.Method.invoke(Method.java:616)
at org.codehaus.plexus.classworlds.launcher.Launcher.launchStandard(Launcher.java:329)
at org.codehaus.plexus.classworlds.launcher.Launcher.launch(Launcher.java:239)
at org.codehaus.plexus.classworlds.launcher.Launcher.mainWithExitCode(Launcher.java:409)
at org.codehaus.plexus.classworlds.launcher.Launcher.main(Launcher.java:352)
at org.codehaus.classworlds.Launcher.main(Launcher.java:47)
postgres@kamil-MS-1675:/home/kamil/osmosis-0.40.1/script$

Proszę o pomoc

serekmedia · 2012-08-24T21:22:44.000Z

Użyj --write-psql bo taki masz załadowany schemat bazy.

kamilek · 2012-08-27T09:21:53.000Z

Dalej jest coś nie tak. Wszystkie czynności które wykonuje:

createuser db_user
postgres=# ALTER USER db_user WITH ENCRYPTED PASSWORD ‘password’;
createdb db_osm
createlang plpgsql db_osm
psql -d db_osm -f postgis.sql
psql -d db_osm -f spatial_ref_sys.sql
psql db_osm -f pgsimple_schema_0.6.sql

postgres@kamil-MS-1675:/home/kamil/osmosis-0.40.1$ osmosis --read-xml file=/home/kamil/osmosis-0.40.1/planet.osm --write-psql host=“localhost” database=“db_osm” user=“db_user” password=“haslo”;
Aug 27, 2012 11:18:51 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Osmosis Version 0.40.1
Aug 27, 2012 11:18:51 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Preparing pipeline.
Aug 27, 2012 11:18:51 AM org.openstreetmap.osmosis.core.Osmosis main
SEVERE: Execution aborted.
org.openstreetmap.osmosis.core.OsmosisRuntimeException: Task type write-psql doesn’t exist.
at org.openstreetmap.osmosis.core.pipeline.common.TaskManagerFactoryRegister.getInstance(TaskManagerFactoryRegister.java:60)
at org.openstreetmap.osmosis.core.pipeline.common.Pipeline.buildTasks(Pipeline.java:50)
at org.openstreetmap.osmosis.core.pipeline.common.Pipeline.prepare(Pipeline.java:112)
at org.openstreetmap.osmosis.core.Osmosis.run(Osmosis.java:86)
at org.openstreetmap.osmosis.core.Osmosis.main(Osmosis.java:37)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:57)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
at java.lang.reflect.Method.invoke(Method.java:616)
at org.codehaus.plexus.classworlds.launcher.Launcher.launchStandard(Launcher.java:329)
at org.codehaus.plexus.classworlds.launcher.Launcher.launch(Launcher.java:239)
at org.codehaus.plexus.classworlds.launcher.Launcher.mainWithExitCode(Launcher.java:409)
at org.codehaus.plexus.classworlds.launcher.Launcher.main(Launcher.java:352)
at org.codehaus.classworlds.Launcher.main(Launcher.java:47)

Dotevo · 2012-08-27T09:31:32.000Z

moze być schemat --write-pgsimp chociaż łatwiej byłoby gdybyś napisał jakiego schematu użyłeś.

To ten? http://wiki.openstreetmap.org/wiki/Osmosis/Detailed_Usage_0.40#PostGIS_Tasks_.28Simple_Schema.29

kamilek · 2012-08-27T09:36:51.000Z

tego użyłem psql db_osm -f pgsimple_schema_0.6.sql

ok, zmieniłem na pgsimp i ruszyło

osmosis --read-xml file=/home/kamil/osmosis-0.40.1/planet.osm --write-pgsimp host=“localhost” database=“db_osm” user=“db_user” password=“haslo”;
Aug 27, 2012 11:37:27 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Osmosis Version 0.40.1
Aug 27, 2012 11:37:27 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Preparing pipeline.
Aug 27, 2012 11:37:28 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Launching pipeline execution.
Aug 27, 2012 11:37:28 AM org.openstreetmap.osmosis.core.Osmosis run
INFO: Pipeline executing, waiting for completion.

kamilek · 2012-08-27T11:19:56.000Z

Ok tak wygląda moja baza.

           List of relations

Ok, rzeczą którą teraz chce zrobić to wyjąć jakąś tam drogę z bazy. Załóżmy że to będzie bbox ograniczony współrzędnymi: 52.679504,22.824512 52.67581,22.848458 52.668003,22.841249 52.669304,22.825713

(by Serekmedia)
biorę punkt ze śladu GPS, robię wokoło niego bbox i szukam dróg przecinających ten bbox (w postgis robi to odpowiednia funkcja np:
select way from planet_roads where ST_Intersects(way,ST_GeomFromText(‘POLYGON((14 49, 24 49, 24 55, 14 55, 14 49))’,4326))
wyjmie wszystkie drogi w zgrubnym BBOX Polski. Oczywiście szukając dróg do traffica używam malutkiego bboxa który obejmuje mniej więcej kwadrat 50x50m.

Rozumiem, że zapytanie powinno wyglądać mniejwięcej tak, ale mam wątpliwości bo w mojej bazie nie ma takiej tabeli jak planet_roads
select way from planet_roads where ST_Intersects(way,ST_GeomFromText(‘POLYGON((22.824512 52.679504, 22.848458 52.67581, 22.841249 52.668003, 22.825713 52.669304, 22.824512 52.679504))’,4326))
Proszę o jakąś podpowiedz

Dotevo · 2012-08-27T11:36:24.000Z

Masz nieco inny schemat. Proponuję abyś sobie dodatkowo wrzucił:


    pgsimple_schema_0.6_bbox.sql - Adds the optional bbox column to the way table.
    pgsimple_schema_0.6_linestring.sql - Adds the optional linestring column to the way table.

Najważniejszy dla Ciebie w tym przypadku to linestring, abyś miał odpowiednią kolumnę w tabeli ways.