település térkép (durvulás)

[Kádár Iván]

Sziasztok!

Páran már ismertek, de azt hiszem ez az első levelem a csoportba, úgyhogy 
kezdem egy kis bemutatkozással.
Fehérváron földmérőként, majd az ELTÉn térképészként végeztem, azóta pedig 
a KSH-ban vagyok térinformatikus, a térstatisztikai rendszer fejlesztésében 
veszek részt. OSM-et eddig nagyon koca módon szerkesztgettem: ha valami 
durva hibát találtam, belenyúltam 1-2 helyen.

Az aktuális problémám: szükségem lenne egy Magyarország közigazgatási 
határai (település) térképre Esri shape-ben (.shp).
*Lehet, hogy elbonyolítottam, várom az ötleteket.*
Kis keresgélés után az tűnt a legegyszerűbbnek, ha a kérésemet beírom az 
overpass-ba, az eredményt meg QGIS-ben átalakítom shp-be.

Kis tanulással erre jutottam:
// Magyarország települései

[bbox:45.7,16.1,48.6,23];
(
  node
      ["ksh_ref"]
     ["is_in:subregion"!="Budapesti"];
);
(
  rel(bn)
    ["boundary"="administrative"]
    ["admin_level"="8"]
);
(._;>;);
out;

illetve kerületekkel:
// Magyarország települései Bp. kerületeivel
[bbox:45.7,16.1,48.6,23];
(
  node
      ["ksh_ref"]
     ["name"!="Budapest"];
)->.kp;
(
  rel(bn.kp)
    ["boundary"="administrative"]
    ["admin_level"="8"];
  rel(bn.kp)
    ["boundary"="administrative"]
    ["admin_level"="9"]
);
(._;>;);
out;
Közben kipucoltam a ksh_ref-tag-et az adatbázisban, és pár településnek 
igyekeztem helyre rázni a középpontját (admin_center)

Tisztább lenne a letöltött OSM forrás, ha
(._;>;);
helyett használhatnám az
out geom;
parancsot, mert ekkor csak a poligonok jönnének le, nem az összes way, node 
minden adatukkal. *De ezt nem nyitotta meg a QGIS. Ezzel van valakinek 
tapasztalata?*

Sajnos a határoknak csak nevük van, a ksh_ref az admin_center node-on van.
Ezért letöltöttem egy táblázatot is név-ksh_ref kapcsolathoz:
[out:csv(::id,"name","ksh_ref")];
(
  node
      ["ksh_ref"]
);

out;

QGIS-be betöltöttem az osm-et, hozzákapcsoltam a csv-t, és az eredményt 
elmentettem shp-be.
Azonban még mindig pár problémám:

   1. Gadács bent van Somogyszilen belül, és valamiért kétszer van meg a 
   poligonja - azon túl, hogy Somogyszil szépen lyukas, ahogy kell. Persze 
   kézzel tudom törölni, de jó lenne már jól leszedni, ha lehetséges.
   2. Kőszegszerdahely területe negatív lett, amikor kiszámoltattam a 
   QGIS-szel.
   3. Ezen kívül van kb. 30 település, aminek a területe több mint 10%-kal 
   eltér a hivatalostól, de ezt majd egy másik posztba teszem.
   
Szóval jó lenne valami egyszerű mód, amivel elő tudok állítani poligon 
shp-et, ami tartalmazza az admin_level=8 (esetleg 9 - Bp kerületei) 
határokat, és a ksh_ref-eket is. Van ötletetek?

Tudom, hogy durva lett, remélem sikerül vele megszólítanom a célközönséget. 
:)

-- 
 Iván


--------- következő rész ---------
Egy csatolt HTML állomány át lett konvertálva...
URL: <http://lists.openstreetmap.org/pipermail/talk-hu/attachments/20150914/67048376/attachment.htm>