Tag hospital: POI: Overpass-Turbo.eu vs osmfilter vs osmconvert

tagtheworld · 2020-08-07T12:44:53.000Z

hallo Community,

ich arbeite derzeit daran, die beste Methode herauszufinden, mit der ich zum Thema “Krankenhaus” weiterkomme. Beabsichtigt ist, eine Übersicht auf europäische Krankenhäuser herzustellen. M.a.W.: ich brauche einen Datensatz der ggf. einfach so aussieht:

 @id    @type    name    addr:postcode    addr:city    addr:street    addr:housenumber    website     contact:email=*
2656877    relation    Unfallkrankenhaus Lorenz Böhler    1200    Wien    Donaueschingenstraße    13    [https://www.auva.at/ukhboehler/](https://www.auva.at/ukhboehler/)

dies wurde mittels Overpass-Turbo.eu und einer dortigen Abfrage gewonnen (vgl. unten). Die Frage ist,
**a. **kann ich diese Overpass-Turbo-Abfrage auf kpl. Europa ausbauen - m.a.W. die Abfrage so formulieren dass ich einen Datensatz erhalte - den lokal speichern kann und dann immer wieder Aktualisierungen in die DB einpflegen kann. oder ist ggf. ein anderes Verfahren -

b. z.B. jenes mit osmfilter zu arbeiten oder mit osmconvert noch bessser und zielführender?!

zu Beginn habe ich mit dem Herunterladen von Auszügen von anderen Diensten wie Geofabrik gespielt. Ich habe selbst keine Daten abgerufen.

Dabei ist die Frage entstanden: Was ist der beste Weg, um dieses Thema anzugehen?

Als Randnotiz; Es besteht eine hohe Wahrscheinlichkeit, dass beispielsweise viele nicht-krankenhausmedizinische Einrichtungen kartiert wurden. Ein weiteres Problem, das ebenfalls häufig auftritt, sind mehrere Krankenhausobjekte, die an einem Krankenhausstandort vorhanden sind:

https://wiki.openstreetmap.org/wiki/Downloading_data#Huge_amounts_of_data

**einige Ansätze: **

osmconvert: https://wiki.openstreetmap.org/wiki/Osmconvert
osmfilter: https://wiki.openstreetmap.org/wiki/Osmfilter
XAPI, Overpass API: https://wiki.openstreetmap.org/wiki/Xapi
overpass-turbo: http://overpass-turbo.eu/
osmium: https://osmcode.org/osmium-tool/

“XAPI, Overpass API and the Export Tool allow to download custom data sets like arbitrary bounding boxes, elements with specific tags, public transport networks or other features.” might take a bit of experimenting with parameters and validation because you’re relying on consistent tagging by other users

Beispiele: osmium https://osmcode.org/osmium-tool/ osmfilter https://wiki.openstreetmap.org/wiki/Osmfilter

healthsites.io: hier könnte ich mit den Tags arbeiten: hospital? amenity=hospital
Weitere Daten & Anleitungen - auch zur API: healthsites.io : etc https://healthsites.io/map

mit overpass turbo http://overpass-turbo.eu/ habe ich insgesamt sehr gute Resultate:

[out:csv(::id,::type,“name”,“addr:postcode”,“addr:city”,“addr:street”,“addr:housenumber”,“website”,“contact:email=*”)][timeout:600];
area[“ISO3166-1”=“AT”]->.a;
( node(area.a)[amenity=hospital];
way(area.a)[amenity=hospital];
rel(area.a)[amenity=hospital];);
out;

Resultate:

 @id    @type    name    addr:postcode    addr:city    addr:street    addr:housenumber    website     contact:email=*
2656877    relation    Unfallkrankenhaus Lorenz Böhler    1200    Wien    Donaueschingenstraße    13    [https://www.auva.at/ukhboehler/](https://www.auva.at/ukhboehler/)   
2685003    relation    Landeskrankenhaus Hörgas-Enzenbach                    [http://www.lkh-hoergas.at](http://www.lkh-hoergas.at)   
2685004    relation    Landeskrankenhaus Judenburg-Knittelfeld                    [http://www.lkh-judenburg.at](http://www.lkh-judenburg.at)   
2685005    relation    Landeskrankenhaus Hochsteiermark                    [http://www.lkh-hochsteiermark.at](http://www.lkh-hochsteiermark.at)   
2685006    relation    Landeskrankenhaus Mürzzuschlag-Mariazell                    [http://www.lkh-muerzzuschlag.at](http://www.lkh-muerzzuschlag.at)   
2685007    relation    Landeskrankenhaus Rottenmann-Bad Aussee                    [http://www.lkh-rottenmann.at](http://www.lkh-rottenmann.at)   
2755244    relation    Landesklinikum Scheibbs    3270    Scheibbs    Feldgasse    26       
2764083    relation    Allgemeines Krankenhaus Wien    1090    Wien    Währinger Gürtel    18-20    [https://www.akhwien.at](https://www.akhwien.at)   
2882269    relation    Unfallkrankenhaus Salzburg                    [http://www.ukh-salzburg.at](http://www.ukh-salzburg.at)   
3333391    relation    Krankenhaus St. Josef    5280    Braunau am Inn    Ringstraße    60    [https://www.khbr.at/](https://www.khbr.at/)

weiteres Beispiel:

30850910	way	Hôpital Nord Guillaume et René Laennec	44093	Nantes Cedex 1	Bd professeur Jacques-Monod - Saint-Herblain		[https://www.chu-nantes.fr/](https://www.chu-nantes.fr/)	
30979413	way	Hôpital Albert Chenevier	94010	Créteil	rue de Mesly	40	[http://chu-mondor.aphp.fr/hopital-albert-chenevier1/](http://chu-mondor.aphp.fr/hopital-albert-chenevier1/)	
31150222	way	Centre de Psychiatrie Les Tournesols						
31150223	way	Clinique du Saint-Coeur					[http://www.saint-coeur.groupe-elsan.com/](http://www.saint-coeur.groupe-elsan.com/)	
31252310	way	Hôpital Avicenne					[http://chu93.aphp.fr/hopital-avicenne-bobigny/](http://chu93.aphp.fr/hopital-avicenne-bobigny/)	
32229408	way	Centre hospitalier du Rouvray	76300	Sotteville-lès-Rouen	Rue Paul Éluard	4	[www.ch-lerouvray.fr](http://www.ch-lerouvray.fr)	
32354044	way	Centre hospitalier intercommunal Fréjus Saint-Raphaël					[https://www.chi-fsr.fr/](https://www.chi-fsr.fr/)	
32388820	way	Clinique des Cèdres					[https://clinique-cedres.capio.fr/](https://clinique-cedres.capio.fr/)	
32471086	way	Hôpital Saint-Jacques						
32522247	way	Centre hospitalier de Dieppe	76200	Dieppe	Avenue Pasteur		[http://www.ch-dieppe.fr](http://www.ch-dieppe.fr)	
32617022	way	Hôpital Désandrouins	55100	Verdun	Avenue d'Étain			
32679701	way	Hôpital Belle-Isle					[https://www.hpmetz.fr/fr/etablissement/hopital-belle-isle](https://www.hpmetz.fr/fr/etablissement/hopital-belle-isle)	
32851894	way	Clinique d'Occitanie	31600	Muret	Avenue Bernard IV	20	[http://www.clinique-occitanie.com/](http://www.clinique-occitanie.com/)

Die Frage ist, ist ggf. der Ansatz mit overpass-turbo.eu zu arbeiten viel vielversprechender als andere Tools und Ansätze wie
osmconvert & osmfilter. Denn mit overpass-turbo scheine ich sehr klar festlegen zu können, was ich brauche.

wie seht ihr das!?

Wenn die Daten einmal heruntergeladen sind, dann sollen sie alle in einer Datenbank (MySQL oder PostgreSQL db) gepflegt werden.

Freue mich auf einen Tipp,

vg Tagtheworld

toc-rox · 2020-08-07T13:08:36.000Z

Ich würde es mit osmconvert und osmfilter versuchen:

osmconvert mit der Option --all-to-nodes um alle Objekte in POIs umzuwandeln
osmfilter um die gewünschten Objekte herauszufiltern

tagtheworld · 2020-08-07T20:43:45.000Z

hallo und guten Abend Klaus,

vielen Dank für Deine rasche Antwort und die Tipps. Das werde ich machen. ich werde mit

1 . osmconvert anfangen - und “mit der Option --all-to-nodes um alle Objekte in POIs umzuwandeln” und dann
2 . mit osmfilter weiterarbeiten u. die gewünschten Objekte herausfiltern…

dann werde ich die Daten alle in einer db. speichern

vielen Dank noch einmal!!!
#tagtheworld

http://manpages.ubuntu.com/manpages/trusty/man1/osmconvert.1.html

   --all-to-nodes

          Some applications do not have the ability to process ways or relations,  they  just
          accept  nodes as input. However, more and more complex object are mapped as ways or
          even relations in order to get all their details into  the  database.   Apply  this
          option  if  you  want  to convert ways and relations to nodes and thereby make them
          available to applications which can only deal with nodes.  For each way a  node  is
          created.  The  way's id is increased by 10^15 and taken as id for the new node. The
          node's longitude and latitude are  set  to  the  way's  geographical  center.  Same
          applies to relations, however they get 2*10^15 as id offset.

Anm; ich hatte mal damit - also mit diesen Themen schon zu tun - vor einigen Jahren - im Jahre 2014. vgl.: overpass-api or osmconvert: which method is more appropiate to use for files from 10 to 400 mb?

https://gis.stackexchange.com/questions/95635/overpass-api-or-osmconvert-which-method-is-more-appropiate-to-use-for-files-fro

**Zu den Optionen: **
Der XML-output der Overpass Api hat im Grunde die gleiche Struktur wie alle .osm files, man kann deshalb auch gut osmconvert einsezten.
bzgl. des Einsatzes einer mysql DB Export XML Data to a MySQL DB ( Link: https://gis.stackexchange.com/questions/94534/export-xml-data-to-a-mysql-db )

Alternativ kann man eine Postgis-DB mit osm2pgsql einsetzen oder eine sog. spatialite database.
CSV files sind eher unpraktisch.

GDAL für den Import .osm data bzw. für den Export in andere Formate (incl CSV).

**Anm: **eine Postgis DB mit osmosis hat einen anderen Setup als eine osm2pgsql. Die osm2pgsql version ist für GIS applicationen gemacht, wie etwa QGIS und rendering tiles mit mapnik. Ich werde mal sehen wie ich weiterkomme - und ob ich PostGIS tatsächl. brauche.

See this Q&A for some guidance:

https://stackoverflow.com/questions/20125242/import-openstreetmap-data-osm-file-to-a-mysql-database-sql-file

http://forum.openstreetmap.org/viewtopic.php?pid=152695

http://mysqlserverteam.com/mysql-5-7-and-gis-an-example/

MySQL allows only the xml formats as stated in the linked document. You have to modify the osm xml into one of those. Alternatively, you can use osmosis or GDAL to do the import.

Zusatz: …um das ganze anzugehen, hab ich noch ein paar Links gesammelt - bzgl. OSM im Einsatz mit Python:

Using Osmconvert with Python
https://stackoverflow.com/questions/26513716/using-osmconvert-with-python

cmd = r"““c:\temp\osmconvert.exe c:\temp\770.osc.gz -b=1,1,3,3 -o=c:\temp\extract.o5m””"

https://github.com/ElcoK/osm_clipper

OpenStreetMap clipper: Small python tool to clip pre-defined areas from large osm.pbf files, using osmconvert or osmosis. The tool is at the moment specifically written to split OpenStreetMap Protocolbuffer Binary Format (PBF) files into a set of countries and/or regions, based on the GADM36 classification. Future versions will allow for user-defined regional classifications as well. Please note: This package is still in development phase. In case of any problems, or if you have any suggestions for improvements, please raise an issue.

**OSM and Python: **To be able to extract the data from OpenStreetMap (OSM), a few steps and downloads are required.
https://gmtra.readthedocs.io/en/latest/osmtutorial.html

Loading Data from OpenStreetMap with Python and the Overpass API
https://towardsdatascience.com/loading-data-from-openstreetmap-with-python-and-the-overpass-api-513882a27fd0