OpenStreetMap

{{
- display_map:
50.1040, 14.3890 | width=307 | height=200 | zoom=14 | service=openlayers }}

OpenStreetMap je projekt, jehož cílem je tvorba volně dostupných geografických dat a následně jejich vizualizace do podoby silniční mapy, uliční mapy měst atd. Pro tvorbu geodat se využívá zejména záznamů z přijímačů GPS, která jsou následně kontrolována a editována. Je založen na kolektivní spolupráci a na koncepci Open source. Data jsou poskytována pod licencí Creative Commons Attribution-ShareAlike 2.0. OpenStreetMap byl inspirován projekty jako je například Wikipedia, umožňuje editaci dat a uchovává kompletní historii provedených změn.
Převzato z české wikipedie.
Výuka
Cvičná databáze pgis_osm
Na serveru geo102 je umístěna databáze pgis_osm s importovanými OSM daty ČR. Databáze slouží pro výuku je určena pro experimentování.
pgis_osm=> \d List of relations Schema | Name | Type | Owner --------+------------------+-------+------- public | czech_line | table | landa public | czech_point | table | landa public | czech_polygon | table | landa public | czech_roads | table | landa public | geometry_columns | table | landa public | spatial_ref_sys | table | landa (6 rows)
Příklad vytvoření tématické vrstvy 'lesni_porosty' (další příklady).
CREATE TABLE lesni_porosty AS
SELECT osm_id,name,way FROM czech_polygon WHERE landuse = 'forest';
ALTER TABLE lesni_porosty ADD PRIMARY KEY (osm_id);
SELECT Populate_geometry_columns('lesni_porosty'::regclass);
CREATE INDEX lesni_porosty_gist ON lesni_porosty USING GIST (way);
-- přidělení práv pro QGIS
GRANT SELECT ON lesni_porosty TO postgis;
Příklad přidání sloupce do tabulky.
ALTER TABLE lesni_porosty ADD COLUMN wood text;
UPDATE lesni_porosty SET wood = czech_polygon.wood FROM czech_polygon
WHERE lesni_porosty.osm_id = czech_polygon.osm_id;
Vzdálená vizualizace v QGIS.
![]() |
![]() |
Příklad vizualizace dat pomocí UMN MapServer.
MAP NAME "Cvicna databaze OSM" STATUS ON IMAGETYPE PNG EXTENT 1347456.13 6198566.23 2099386.47 6630407.12 SIZE 800 600 IMAGECOLOR 255 255 255 UNITS meters WEB IMAGEPATH "/var/tmp/ms_tmp/" IMAGEURL "/tmp/" END LAYER NAME lesy CONNECTIONTYPE postgis CONNECTION "user=postgis password=XXXXXXXX dbname=pgis_osm host=localhost port=5432" DATA "way FROM czech_polygon" FILTER "landuse = 'forest'" STATUS DEFAULT TYPE POLYGON CLASS NAME "Lesy" STYLE COLOR 89 167 89 OUTLINECOLOR 32 32 32 END END END END
Import
Import dat do PostGIS
Nejprve stáhneme zkomprimovaná data ve formátu OSM, viz [1]. Např.
wget http://download.geofabrik.de/osm/europe/czech_republic.osm.bz2
Projekce Google Mercator ve standardní distribuci PostGISu může chybět, potom je potřeba přidat do tabulky spatial_ref_sys příslušný záznam.
INSERT INTO spatial_ref_sys VALUES (900913, 'local', 900913,
'PROJCS["Google Mercator",GEOGCS["WGS 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",
6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,
AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],
AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4326"]],
PROJECTION["Mercator_1SP"],PARAMETER["semi_minor",6378137.0],
PARAMETER["latitude_of_origin",0.0],PARAMETER["central_meridian",0.0],PARAMETER["scale_factor",1.0],
PARAMETER["false_easting",0.0],PARAMETER["false_northing",0.0],UNIT["m",1.0],AXIS["Easting",EAST],
AXIS["Northing",NORTH],AUTHORITY["EPSG","900913"]]',
'+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs');
Pomocí konzolové aplikace osm2pgsql data naimportujeme do existující databáze PostGIS (Google Mercator).
osm2pgsql -d pgis_osm -p czech -s czech_republic.osm.bz2
Příklad pro EPSG 4326 (WGS-84):
osm2pgsql -d pgis_osm -p czech -E 4326 -s czech_republic.osm.bz2
Pro import do PostGISu můžete kromě osm2pgsql použít také Imposm.
Import dat do GRASS GIS
Poznámka: Návod pro GRASS 7
Nejprve stáhneme zkomprimovaná data ve formátu ESRI Shapefile, viz [2]. Např.
mkdir ~/osm cd ~/osm wget http://download.geofabrik.de/osm/europe/czech_republic.shp.zip
Data jsou referencována v souřadnicovém systému WGS-84 (EPSG:4326), před importem data dekomprimujeme, např.
unzip czech_republic.shp.zip
Spustíme GRASS
grass70
Novou GRASS lokaci můžeme vytvořit pomocí wxGUI nástroje Location Wizard nebo z běžící seance pomocí příkazu g.proj.
g.proj location=cr_osm wkt=~/osm/buildings.prj
GRASS opustíme a spustíme znovu se zvolenou lokací 'cr_osm'.
Archiv obsahuje několik tématických vrstev
v.in.ogr -l dsn=~/osm Data source contains 7 layers: roads railways places waterways points buildings natural
Import dané vrstvy provedeme příkazem v.in.ogr.
v.in.ogr dsn=~/osm layer=buildings
Import může trvat poměrně dlouho (navíc je náročný na paměť), modul při importu opravuje topologické chyby, které se ve vstupních datech mohou vyskytnout.
Příklad Python skriptu pro import data:
#!/usr/bin/env python
#%module
#% description: Import OSM data provided in ESRI Shapefile format into GRASS
#%end
#%option
#% key: directory
#% description: Name of input directory
#% type: string
#% required: yes
#%end
import sys
import grass.script as grass
def main():
layers = grass.read_command('v.in.ogr',
quiet = True,
flags = 'l',
dsn = options['directory']).splitlines()
if not layers:
grass.fatal("No layers available")
for l in layers:
grass.message("Importing <%s>" % l)
grass.run_command('v.in.ogr',
dsn = options['directory'],
layer = l)
return 0
if __name__ == "__main__":
options, flags = grass.parser()
sys.exit(main())
Odkazy
- Přednáška z předmětu Free software GIS
- Projekt OpenStreetMap z pohledu geoinformatika
- Wiki projektu OpenStreetMap
- Článek na abclinuxu.cz
- Rozšířené možnosti editace OSM pomocí editoru JOSM
- Blog o OSM
- Promo materiály
- OpenStreetMap Basics (videos)
- switch2osm - Make the switch to OpenStreetMap
- Project OpenStreetMap and Open and Free Source of Geodata and Maps
- Rendering OpenStreetMap Data with MapServer
- http://yournavigation.org/
- http://sourceforge.net/projects/osm3d/
- http://www.maposmatic.org
- GeoFabrik Tools
- Česká e-mailová konference
- Humanitarian OpenStreetMap Team
- Open Street Map community responds to Haiti crisis
- Příspěvky z konference SotMEU11
- OSM History
- How did you contribute to OpenStreetMap?
- http://openrouteservice.org/
- http://map.project-osrm.org/
- Vespucci - online editor for OSM data