Vektorizacia rastrovych map

[V12]

Tiez som sa niecim podobnym zaoberal, no filtrovanie a upravy robil vo fotosope a CorelTrace. Bohuzial bol vysledok tragicky. jedine rozumne riesenie sa mi javila vektorizacia cez stredovu ciaru a nasledny export do DXF. A prave tu bol kamen urazu - DXF je sice fajn vektorovy format, no nie je georeferencovatelny. A napisat konvertor som doteraz este nemal cas hoci premyslene to uz mam - v DXF najdem lavy dolny extrem, pravy horny extrem a z tychto bodov sa vypocitaju koeficienty, ktorymi upravim vsetky X a Y suradnice v DXF. Nasledne by sa to dalo vysypat napr. do GPX/PLT trackov po jednotlivych polylines, takze by to bolo spracovatelne priamo v OZI. Cely problem je ale v tom, ze je tam privela manualnej roboty, neviem to zautomatizovat. Okrem toho by bol urcite problem s dorovnanim projekcie - toto by slapalo iba na JTSK.

[bogi]

suttorad, to znie prilis zlozite.
V georeferencovani v JTSK suradniciach nevidim problem.

[lacop]

Ešte raz ďakujem za rady (hlavne Bogi). Tu je dôkaz, že má zmysel pokračovať.

[bogi]

Vidim tam .plt 

Ako si to georeferencoval?

EDIT:
Clovece, natiahol som si ten .PLT do Mapy les. ciest ... Fakt dobra robota! 

[lacop]

Mozny postup georeferencovania vo WinTopo Pro
1. Pixelove a WGS84 suradnice LH, PH, PD, LD bodu bitmapy precitas z Oziho kalibracneho suboru. Ide o sekcie MMPXY a MMPLL. Tieto udaje zadas vo WinTopo Pro v dialogu Georeference Data.
2. Vektorizaciou vzniknu georeferencovane utvary (shapes), ktore budu obsahovat WGS84 suradnice. Tie je vhodne vo WinTopo ulozit v .shp formate.
3. Shapefiles ide natiahnut do OziExplorera a v nom skonvertovat na GPX, alebo PLT.

No predsa podľa tvojich rád. Zatiaľ je to len pokus pre malú mapu (1536 x 2048 pixelov, 6,6 m na pixel). O nejakých konverziách v projekcii alebo dátume mapy som preto zatiaľ neuvažoval.

[bogi]

Lat/Lon nie je pravouhly grid - poludniky nie su v hornej a dolnej casti mapy od seba rovnako daleko, su to vlastne krivky. Ked aj zadas rohove body spravne, zavislost zmeny Lat/Lon na pixel nie je medzi rohovymi bodmi celkom linearna. JTSK grid je navyse pootoceny - na mape les. ciest zdola hore nie je presny sever. Pri vacsom mapovom liste by som preto ocakaval nejake posunutia vysledku voci origo mape.

Avsak ziaden posuv neocakavam, keby si robil georeferencovanie v JTSK suradniciach. To by fungovalo dobre aj s vacsou mapou. Dovod je ten, co som pisal: Mapa lesnych ciest je ulozena v pravidelnej (linearnej) pravouhlej mriezke JTSK suradnic. Cize idealne prostredie pre georeferencovanie vo WinTopo, v ktorom som nenasiel moznost zvolit projekciu.
Takymto sposobom by vsetky "trasove body" tvojho teraz uz PLT suboru najprv "nabrali" JTSK suradnice z tejto rovnomernej pravouhlej mriezky a potom konverziou v sk2sk ziskali zodpovedajuce WGS84 suradnice.

Proste vyskusaj s vacsou mapou a uvidis.

[lacop]

To si uvedomujem, že pri väčšej mape to nebude také jednoduché. Budem v tom pokračovať. Pre mňa sú samozrejme zaujímavé väčšie mapy, napr. 40 x 40 km. Určite tvoju radu využijem.

[lacop]

Rady uvedené vyššie som využil a dospel som k výsledku, s ktorým som spokojný. Ďakujem.

Podkladom mi bola mapa lesných ciest pre územie od Košíc po Krompachy, ktorú som si svojho času vyrobil. "Červené" lesné cesty som odfiltroval v bitmapovom grafickom editore Photofiltre Studio pomocou nástrojov na znižovanie farebnej hĺbky a na nahradzovanie farieb. Výsledný obrázok som vektorizoval v programe WinTopo Professional, kde som ho aj georeferencoval v JTSK súradniciach podľa BOGI-ho rád. Na prevod WGS84 -JTSK som použil on-line tool z adresy http://www.sopsr.sk/istb/index.php?p=2. Následne po vektorizácii som výsledok uložil v SHP formáte. Súbor SHP som potom pomocou programu TRANSDAT (nastavenie je na obrázku v prílohe) konvertoval z JTSK na WGS84. Tento som otvoril v OziExploreri a uložil v PLT formáte. Výsledok je na stiahnutie tu : KE-KROMPACHY.plt (2,84 MB)

... takže si to teraz skúsiť ešte raz so "zelenými" lesnými cestami.

[bogi]

laco, si sikovny! 

Konvertovat suradnice v shapefile sa dalo aj cez sk2sk, ale aspon som videl, ako vyzera Transdat. 

Keby sa tak chcelo niekomu urobit oblast BA - Pezinska Baba, alebo este lepsie aj par km dalej ...
Pocas dlhych zimnych vecerov sa mozno do toho pustim. Ale nebude mi vobec vadit, keby sa nasiel niekto iniciativnejsi. 
Keby bolo treba, budem zase kibycovat.

[lacop]

Povzbudený som narýchlo skúsil urobiť tie "zelené". Tu je výsledok: KE-KROMPACHY_zelene.plt (7,1 MB)

Ako vyzerajú zelené a červené lesné cesty v mape VKÚ:




Nejaké info o našich výsledkoch som dal aj na MTBiker.sk: http://forum.mtbiker.sk/viewtopic.php?f=12&p=224077#p224077

[bogi]

Tak ma Laco navnadil a zacal som sa aj ja hrat s vektorizaciou lesnych ciest.
Zistil som, ze vobec to neboli. 
Pouzil som free GIMP na extrakciu ciest podla farby & shareware WinTopo Pro na vektorizaciu a georeferencovanie.

Moj vysledok:



Posunut to o pixel doprava ci dole nie je problem ...
Vsimnite si, ze trasove body nie su vytvarane zbytocne husto. V priamej ciare nie su zbytocne body.

Zelene lesne cesty na ortofotomape s pridanymi vrstevnicami:

[lacop]

Tvoje "zelené" vyzerajú perfektne. Ja so svojimi nie som až tak spokojný. Skúsil by som ich extrahovať v GIMPe. Ktoré tooly v GIMPe si používal a by si ich aj odporučil?

[bogi]

Robil som to takto:
Vysoky zoom (trebarz 400%), volba Vyber | Podla farby, klik na pixeli so zelenou, otvorim dalsiu bitmapu rovnakej pixelovej velkosti, pastnem to tam a necham ukotvit vrstvu. Potom v origo bitmape kliknem na "inej zelenej" (lebo ja mam tmavsiu, aj svetlejsiu zelenu, neviem, ci to vzniklo predoslym spracovanim bitmapy, alebo to tak bolo vzdy) pastnem to do druhej bitmapy, cez Shift-sipky a sipky bez Shift posuniem pastnutu vrstvu, aby pasovala a ukotvim ju. Potom som Dalej robil uz len v tej druhej bitmape, co obsahuje len zelene cesty. Oznacil som zelene prvej farby, dal som Upravit | Vyplnenie farbou popredia - tym sa mi zmenila farba na ciernu. To som spravil aj s "druhou zelenou". Potom som este raz oznacil vsetko teraz uz cierne a pouzil som zaujimavu funkciu v GIMPe - obkolesenie selekcie. Zadal som 1 pixel obkolesenia. Je to pristupne cez Vyber | Editor vyberu a prava ikonka. Neporovnaval som vysledok po vektorizacii s a bez tejto upravy. Porovnaval som len pred vektorizaciou. Po tom obkoleseni vzniknu sice tucnejsie ciary, ale su menej zubate, su take spojitejsie a taktiez sa zaplnili take tesne medzery, ktore v zelenych niekedy su. Po tom obkoleseni opat nebolo vsetko cierne, tak som tie novovzniknute okraje pozmenil na farbu popredia (ciernu). Vzniknu pomerne hrube ciary, ale su celkom fajn spojite a nie su zbytocne zubate.

Dalej som pokracoval vo WinTopo. Potesilo ma, ze som nasiel nastroj na stencenie ciar. Experimentovanim som zistil, ze z Image | Thinning Methods je najlepsia Stentiford Method. Zvolil som Vector | One Touch Vectorisation (bez zmeny defaultov) a potom som georeferencoval ako Ty. Na zaver som ulozil .shp.

[lacop]

Díky, budem ešte tie zelené vylepšovať.

[bogi]

lacop, diki za obrazok s TransDat!

Dnes som sa totiz pustil do vektorizacie Malych Karpatov a Zahoria. Najprv som to len zo srandy skusil georeferencovat priamo v LatLon suradniciach. Vysledok na velkej mape nebol dobry. Niekde to vcelku sedelo, niekde inde bol posun aj niekolko desiatok px. Potom som to georeferencoval este raz v JTSK suradniciach. Pouzil som sk2sk na konverziu suradnic z WGS84 do S-JTSK. Potom som pustil konverziu georeferencovaneho .shp v sk2sk. Zial, sk2sk sa na tomto zrubal. Tak som pokracoval s TransDat podla Tvojho obrazku. Vysledok je velmi dobry, na celej mape je nepresnost na urovni 1 px.

Zatial mam georeferencovane cervene, teraz idem na zelene.