Početkom decembra 2021. smo pokrenuli akciju mapiranja objekata (zgrada) u Mladenovcu. Pisali smo o tome tada. Ovo je prva akcija mapiranja neke teritorije koju smo napravili i nismo znali šta da očekujemo (da budemo iskreni, imali smo pre toga akciju unošenja pošti, ali razmera problema nije ni blizu bila uporediva).

Sa zadovoljstvom možemo da vas obavestimo da je ceo projekat već gotov! Sve smo očekivali, ali da bude gotov ovako brzo – nismo😀. Projekat je trajao dva cela meseca (ceo decembar i ceo januar) i obrađeno je 100% teritorije opštine Mladenovac. Na projektu je učestvovalo 8 ljudi i uneto je čak 37.000 objekata! Ovim putem treba da pomenemo i posebno pohvalimo člana „soliMM” koji je uradio najviše i praktično izneo ceo projekat! Pošto je ceo Mladenovac oko 0,5% teritorije Srbije, završili bismo ovim tempom celu Srbiju za 33 godine (deluje dosta sporo😄, taman da odmah nakon toga uradimo još jedan prolaz😀), ali treba uzeti u obzir da broj ljudi koji mogu da doprinesu može da ide samo na gore.

Pre i posle

Ako ima jedna slika da obeleži ovaj projekat, to je ova (klik za HD verziju od 6 MB!):

Na sreću (i optimistično) smo pre početka projekta uzeli slike Mladenovca, pa možemo da prikažemo ovde sad par uporednih fotki. Zato smo odradili par GIF-ova koji treba najbolje da dočaraju ove pre-posle momente, nalik ovako jednom (klik na veću rezoluciju):

Pošto su slike ogromne i nismo želeli da ih uključimo u ovaj članak (da ne opterete Vaš browser i protok), ovde možete naći linkove na ono što smo smatrali najzanimljivijim. Ne zaboravite da svaku sliku možete da zumirate i da se najbolje vide na 100% uvećanju:

• Zum 15
  • Opština Mladenovac, 5 MB, 6836 x 3270px
  • Amerić i Vlaška, 500 KB, 2227 x 1103px
  • Crkvine i Markovac, 400 KB, 1785 x 1107px
  • Velika Ivanča, 300 KB, 1355 x 1102px
• Zum 16
  • Mladenovac grad, 1,5 MB, 3078 x 2048px
  • Koraćica, 600 KB, 1846 x 1862px
  • Kovačevac, 800 KB, 1728 x 2958px
  • Pružatovac, 300 KB, 1604 x 1482px
  • Velika Krsna, 650 KB, 2088 x 2020px
• Zum 17
  • Mladenovac grad, 3,1 MB, 5092 x 3584px
  • Kovačevac, 1,8 MB, 3492 x 6105 px
  • Pružatovac, 400 KB, 2223 x 2238 px
  • Velika Krsna, 1,4 MB, 3980 x 3980 px
• Zum 18
  • Mladenovac grad, 6,9 MB, 10616 x 6912px

Uživajte!

Analize

Kada svi objekti ucrtani na mapu postoje, i to ucrtani sa velikim kvalitetom, otvaraju se mogućnosti i za razne analize. Kao što je već rečeno, ukupan broj zgrada u Mladenovcu je sada oko 37.000. To se može videti sledećim Overpass upitom:

[out:json][timeout:15];
area["name"="Gradska opština Mladenovac"]->.searchArea;
(
  way[building](area.searchArea);
);
(._;>;);
out;

Broj objekata po naselju

Pošto smo u prethodne dve godine dobili od RGZ-a i granice svih naselja, o čemu možete da čitate opširno na forumu, takođe se lako vidi i broj zgrada po naseljima u opštini Mladenovac, i on izgleda ovako nekako:

Da biste dobili ovu računicu, možete iskoristiti sledeći Overpass upit:

[out:csv( name, ways)];
area["name"="Gradska opština Mladenovac"]->.searchArea;
(
  relation(area.searchArea)["admin_level"=9]["ref:RS:naselje"];
  map_to_area;
);
foreach->.regio (
  (
    way(area.regio)[building];
  );
  make count
             name = regio.set(t["name"]),
             ways = count(ways);
  out;
);

Naravno, znajući površine naselja (ili broj ljudi u naselju) možemo dalje izračunati i gustinu objekata po naselju, ali to neka bude ostavljeno kao „vežba čitaocu”😀.

Brzinu unosa zgrada na projektu

Za potrebe našeg projekta, a koristeći upit sličan ovom:

[out:json][timeout:15];
area["name"="Gradska opština Mladenovac"]->.searchArea;
(
	way[building](area.searchArea)(newer:"2021-11-01T00:00:00Z");
);
(._;>;);
out;

(samo menjajući datum u upitu), možemo dobiti kako smo unosili objekte na dnevnom nivou:

Ono što se lako primeti je kako smo krenuli polako, vidi se i kada se „soliMM” pridružio (oko 15.12.), vidi se i da nas Novogodišnji i Božićni praznici nisu nešto usporili, a vidi se i jedan period kada smo odlučili da „objavimo primirje” (21.-27.02.) jer smo razmišljali da li da ostavimo nešto zgrada da ih uradimo (i završimo) na nekom mapatonu. Pošto se mapaton nije desio (omikron varijanta je tad dominirala), nastavili smo i dovršili projekat.

Obim i površina svih objekata

Jedna zanimljiva analiza koja može da se uradi je da se izračuna površina svih ovih objekata (tj. površina tlocrta)! Naravno, to umnogome zavisi od kvaliteta ucrtanih zgrada, ali nas ne zanima ovde ništa precizno, već eto za neke okvirne analize. Npr. ono što možemo da izvučemo sa teritorije Mladenovca je da je prosečan obim objekata 42m, medijana obima 40m, a što se tiče površina – prosečna površina je 115m², a medijana 91m². Ovo deluje malo više nego što bi se očekivalo, moguće je da maperi uvek zahvate malo više od objekta dok ga crtaju nego što je realno, pa zato ispadne ovako visoko. Interesantno je videti i histogram površine (x-osa su opsezi površina, a y-osa je broj objekata u tom opsegu):

Orijentacija kuća

Inspirisani projektom orijentacije ulica, bilo je zanimljivo proveriti da li ima ikakve pravilnosti kod orijentacije kuća. Tj. da li se kuće zidaju tako da prozori gledaju na neku određenu stranu sveta. Za ovo je iz QGIS-a iskorišćena metoda „Oriented minimum bounding box” koja za neki poligon izvuče njegov glavni ugao. Ti uglovi idu od 0° do 180°, pa je za potrebe grafičkog prikaza napravljena simetrija (na grafiku su od 180° do 360° iste vrednosti). Kada se iscrta na polarnom grafiku, to izgleda nekako ovako:

Sad, definitivno postoji neka pravilnost… Dominiraju istok-zapad orijentacije (i u manjoj meri ostali pravci). Ovo mi nema baš mnogo smisla, jer ako je duža strana kuće u pravcu istok-zapad, onda prozori većinom gledaju na sever ili jug? Ova druga orijentacija (severozapad-jugoistok) mi ima više smisla, pa čak i ova sever-jug ima najviše smisla jer je valjda najbolje da prozori budu istok i zapad? Severoistok-jugozapad opet i ne baš. Vrlo moguće objašnjenje je da su kuće većinski kvadratne, pa ovo nije simetrija 0°-180° na 180° na 360°, već se 0°-90° preslikava u sva četiri kvadranta (što, gledajući grafik ima smisla donekle). Definitivno se primećuje da dominiraju glavni smerovi (sever, istok, jug, zapad). Bilo bi zanimljivo odraditi istu ovu analizu i za druge opštine, da li bi smo dobili iste rezultate? I ovo odlažemo kao „vežbu čitaocu” da nastavi😀.

Najusamljenija kuća

Jedna interesantnija analiza, u ovim vremenima kada svi gledaju da pobegnu van grada, je naći najizolovaniju kuću u Mladenovcu. Npr. želite vikendicu, želite da je blizu Beograda, ali ne želite ljude oko Vas. Kako naći takvu kuću? Pa lako sada kada imamo ucrtane sve objekte! Prvo se sve kuće učitaju u programu QGIS i onda se uradi analiza „Join attributes by nearest”. Ovim se za svaki objekat nađe objekat najbliži njemu. Pošto van opštine Mladenovac zgrade mahom nisu ucrtane, odlučili smo da i sama granica opštine bude fiktivna „zgrada” za potrebe ove analize (inače ne bi bilo fer pošto bi se u suprotnom kuće najbliže granici opštine vodile kao najudaljenije). Kada se to izračuna, evo top 10 objekata koji se nalaze najudaljeniji od drugih objekata (klik za veću rezoluciju):

Ljubičasto su ucrtani objekti. Analiza može da se izvede verovatno bolje jer npr. dva objekta jedan do drugoga (recimo štala i kuća) poništavaju jedan drugog i, čak iako su udaljeni 2km od sledeće najbliže kuće, neće se videti na ovoj analizi. Da li možete da smislite pametniji pristup? No, neke stvari se naziru i ovde. Npr. ako hoćete da se izolujete na Kosmaju – srećno Vam bilo, to izgleda nemoguće😀. Takođe potez od Mladenovca ka selu Dubona ima najviše ovakvih izolovanih objekata. I ovaj deo severno od sela Jagnjilo isto. Ali apsolutni pobednik je neki objekat udaljen 1174m od prve sledeće kuće – evo je veza ka objektu u OpenStreetMap-ama. Hajde da ga vidimo iz satelita😀 (klik za veću rezoluciju):

Vidite ovu tačkicu u žutom krugu? E, to je možda Vaša nova vikendica😀. Evo malo bliže (opet klik za veću rezoluciju):

Da, definitivno je neki objekat za život, a ne štala, plastenik ili sl. Svaka čast, ko god da je, može da kaže da je najudaljeniji od drugih u Mladenovcu😀.

Najudaljenija tačka od drugih kuća

Sada kada znamo koja je najudaljenija kuća u Mladenovcu, i ako pretpostavimo da ne možemo da je kupimo, sledeće logično pitanje je – a gde da kupimo plac da sagradimo novu takvu kuću? Tj. koja je tačka u Mladenovcu koja je najudaljenija od prve sledeće kuće. Da bismo uradili ovakvu analizu, moramo da koristimo QGIS. Prvo treba da napravimo „grid“ tačaka po Mladenovcu (razmak od 20m od tačke do tačke će biti dovoljan za naše potrebe preciznosti), a kada imamo izgenerisanu mrežu, onda treba da izvršimo istu operaciju „Join attributes by nearest” (samo sad treba da poredimo sloj tačaka sa novonapravljene mreže sa slojem zgrada). Onda ovakvu mrežu treba da rasterizujemo („Rasterize (vector to raster)”. Kada to iscrtamo dobijemo nešto ovako:

Opet dominira deo ka Duboni, deo severno od Jagnjila i deo severno od Vlaške. Pošto ovakvo crtanje raster sloja preko OSM mape ne govori mnogo, možemo nad ovim rasterom da uradimo operaciju „Contour” i da dobijemo linije iste udaljenosti od najbliže kuće. Tako nešto izgleda ovako:

Evo još jedna slika uveličana na ovu zonu sa najvećom udaljenošću:

Kao što se vidi, ima čak jedan deo gde je udaljenost čak preko 1700m!

Sve u svemu, ovi podaci neprocenjivo obogaćuju OSM mapu. Ako ne računamo besmislene analize za traženje vikendica😄, oni imaju široke primene, od orijentacije planinara do doktorskih radova na temu lokacijske procene štete zemljotresa na osnovu broja i tipa objekata (da, i to se radi na osnovu OSM podataka!). I zato treba da poboljšavamo OSM mapu, zgradu po zgradu, opštinu po opštinu!

Sledeća lokacija putujućeg cirkusa?

Kao što smo rekli prošli put, nagrada za onog ko uradi najviše je da bira sledeću lokaciju😄. Da li će Miloš da bira ili ćemo da napravimo neki proporcionalni sistem glasanja će odlučiti Miloš, ali sad ide zaslužena pauza. Na proleće možemo da biramo sledeću lokaciju, do tad – hvala svima na učestvovanju i iskren poziv da se pridružite ovom putujućem cirkusu u sledećoj avanturi mapiranja.