Detekce městského obvodu Plzně ze souřadnic

Webová služba, která na základě vstupních souřadnic ve formátu WGS 84 vrátí výstup s informací ve formátu XML či JSON, v jakém městském obvodu Plzně se daný bod nachází a v jaké jeho části.

Provoz

Webová služba běží na adrese http://tools.jasnapaka.com/mestske-obvody-plzen/service.php.

Instalace

Je vyžadován webový server s podporou PHP, databáze PostgreSQL s nadstavbou PostGIS. Otestována byla následující konfigurace, ale může být i jiná.

Apache 2.4
PHP 7.2 či vyšší (vyžaduje moduly "curl" a "pgsql")
PostgreSQL 9.3 či vyšší
PostGIS k odpovídající verzi PostgreSQL

Tento projekt si stačí stáhnout, a následně přejmenovat soubor config-default.php v rootu projektu na config.php. V něm pak nastavit přístupové údaje k databázi s daty hranic městských obvodů.

Příprava dat

Hranice městských obvodů v Plzni lze volně získat v Databázi otevřených dat města Plzně ve formátu shp. Na témže místě lze získat i části obce. Obojí lze následně naimportovat do databáze pomocí příkazu shp2pgsql, který je součástí PostGisu.

shp2pgsql -s 5514 mestskecasti_SHP.shp public.umo | psql -h localhost -d plzen -U postgres
shp2pgsql -s 5514 castiobce_SHP.shp public.castiobce | psql -h localhost -d plzen -U postgres

Výše uvedený příklad předpokládá, že výsledná databáze se bude jmenovat plzen, informace o městských obvodech se naimportuje do tabulky umo a části obce budou v tabulce castiobce. Tu první následně upravíme tak, aby jednotlivé sloupce byly pojmenovány následovně:

gid
id
nazev
geom
kod

Reálně přejmenujeme sloupec pro název městského obvodu a přidáme sloupec kod. Protože se chybně naimportovala diakritika, bude potřeba ještě upravit názvy ve sloupci nazev. Do nově přidaného sloupce kod je potřeba přidat zkratky městských obvodů. Např. umo3 pro Plzeň 3 a ekvivalentně pro další. U druhé tabulky castiobce pouze upravíme diakritiku ve sloupci nazev.

Pozor! Pokud budete vytvořenou databázi zálohovat a obnovovat jinde, ujistěte se, že v tabulce spatial_ref_sys, která je součástí databáze, bude po obnovení uveden souřadnicový systém S-JTSK EPSG:5514. Stačí si vyfiltrovat hodnotu 5514 ve sloupci srid. Pokud ji zde nenaleznete, stáhněte si následující skript a nad databází jej proveďte:

wget http://epsg.io/5514.sql
psql -U postgres -f 5514.sql plzen

Bez uvedeného souřadnicového systému v databázi nic nenaleznete.

Příklad použití

Výstup získáte zavoláním skriptu service.php, který jako parametr bere lat a long odpovídající hodnotám ze souřadnicového systému WGS 84.

http://mujvlastniserver.cz/service.php?lat=49.738065&long=13.382195

Pokud je vše v pořádku, vrátí se vám XML soubor s HTTP stavovým kódem 200. Jeho podoba bude následující:

<?xml version="1.0" encoding="UTF-8"?>
<area>
    <code>umo3</code>
    <umo>Plzeň 3</umo>
    <part>Jižní Předměstí</part>
</area>

Výstup obsahuje kód městské části ve formátu umoX, kde X odpovídá číslu městské části Plzně. Následuje název městské části.

Pokud nebyl na základě souřadnic nalezen žádný městský obvod či došlo při hledání k chybě, vrací se chybové XML. Stav, kdy nebylo nic nalezeno, je reprezentován HTTP stavovým kódem 404. Chybu pak značí stavový kód 500. V návratovém XML je pak kromě číselného kódu chyby i jeho popis. Jejich kompletní výčet lze nalézt v souboru src/Error.php.

<?xml version="1.0" encoding="UTF-8"?>
<error>
    <code>4</code>
    <msg>Nastala interní chyba služby. Databáze není dostupná.</msg>
</error>

Výstup lze vrátit i ve formátu json. K tomu lze využít parametr URL format, který akceptuje hodnotu jak xml (není třeba uvádět), tak json. Příklad navráceného JSONu:

{"code":"umo3","umo":"Plzeň 3","part":"Doudlevce"}

Případně chybový výstup:

{"code":2,"msg":"Vstupní parametry 'lat' a 'long' musí být reálná čísla."}

Hromadná žádost

Pokud se potřebujete hromadně dotázat na více bodů, je to možné. Postačí zavolat samotný skript a do POST Payloadu uvést JSON s žádostí. Výstupem bude JSON s hromadnou odpovědí. Příklad žádosti:

[{"lat":49.725,"long":13.37661},{"lat":49.761248}]

Odpověď bude:

{"count" : 2, "items" : [{"status" : 200,"code" : "umo3", "umo" : "Plzeň 3", "part" : "Doudlevce"},{"status" : 500, "code" : 2, "msg" : "Vstupní parametry 'lat' a 'long' musí být reálná čísla."}]}