Punktidest aladeks

ehk kuidas kasvatada kringlit rosinate ümber



@tkardi
Tõnis Kärdi
LonLat OÜ
tonis.kardi@gmail.com
Palun, siin on hunnik punkte, meil on vaja atribuutika alusel loodud alasid!
(tegelikult nii ei olnud, aga no oletame)
"Hah! see on ju lihtne!"
Voronoid
ja neid saab genereerida otse PostGISis
Grupeerime tükid atribuudi (või mitme) alusel ja ..
via giphy
aaaaaaaga...
nurga taga on standby režiimis ootel ventilaator
jõudlus?
päris keeruline paralleliseerida arvutamist
andmete uuendused/muudatused?
sisuliselt peab alati terve andmestiku uuesti arvutama
ja tegelikult / peamiselt:
mida need alad üldse peaks näitama?
antud juhul räägime sihtnumbrialadest
vt nt sihtnumbrite reeglistik
Omniva kodulehelt csv-failina allalaetav
seega vaja natuke teistsugust lähenemist
alustame hoopis teisest otsast
teejooned
OSM , ETAK
raudteed
OSM , ETAK
vooluveed
OSM , ETAK , KKR avalik WFS
... ja asustusüksused
EHAK
seejuures
asustusüksuste piirjooned venitame välja veekogudele/riigipiirini
kui huvi, siis selle kohta kuidas, saab lugeda siit
lükkame kõik need andmed kokku
ja saame terve hunniku risti-rästi jooni
ehk "müra"
kuigi tegelikult
joonistuvad siit asulates välja päris ilusad kvartalid
ning tuues pildile tagasi selle loo peategelased
värvitud sihtnumbri järgi
seega:
oleks vaja joontest mingisugused mõistlikud polügonid saada
ja seejärel neisse jäävate punktide alusel alad lihtsalt klassifitseerida
appi tuleb PostGIS
- jeee! :) -
esiteks
purusta mürajooned kõigis omavahelistes ristumiskohtades.
pm sama, mis teeksid teejoontest navi-kõlbulike andmete loomisel
teiseks
kombineeri tekkinud lõigud omavahel kokku polügonideks (kutsume neid kvartaliteks)
vt lisaks st_polygonize, st_node, st_snap
kolmandaks
topokontrollid - auke ei tohiks olla, ülekatteid ei tohiks olla
nt GRASSGIS
ja siit edasi klassifitseerime kvartali
temasse jäävate punktide atribuutika (nt sihtnumber) järgi
rääkides sihtnumbritest
siis tegelikult on
väga suures enamuses asustusüksuse kohta kasutusel vaid üks (kandeala) sihtnumber
nii et tegelikult saame me enamjagu aladest paika juba eelnevalt
kuid edasi...
esiteks
üks kvartal - üks sihtnumber
värvides need sihtnumbri väärtuse järgi
teiseks
üks kvartal - mitu sihtnumbrit
(peame jagama kvartali osadeks)
Voronoid aadresspunktidele iga kvartali kaupa?
tänu katastriandmete avamisele
MAA-AMET, AITÄH! :)
saame täpsemal jagamisel kasutada katastriüksuseid.
siin samuti
üks kataster - üks sihtnumber *
ja pannes kaks kihti kokku saame
oot-oot-oot!!!
katastri järgi pole terve kvartal ju täidetud!
täpselt nii...
pole hullu, jagame kvartali tühja ruumi selles olevate alade vahel lihtsalt ära
...
(järgneb lõikamine, kleepimine, saagimine ja haamrikoputused PostGISis)
kolmandaks
üks kvartal - mitte ühtegi sihtnumbrit
alustame aukudest
kui klassifitseerimata kvartalit ümbritseb ainult sama sihtnumbriga kvartal / kvartalid, siis määrame samasse
samas ülejäänud on samuti augud
lahendame nad siis ühe käiguga
kaks võimalust
- jagame võrdselt piirnevate kvartalite vahel
- klassifitseerime selleks, millega pikim ühine piir
AITÄH!
UEC OÜ, Omniva
võimaluse eest selle ülesandega tegeleda
@tkardi
tonis.kardi@gmail.com