"Masowa" edycja danych OSM

kocio · 2018-08-22T19:42:49.000Z

Wybrać budynki z bazy wraz z położeniem punktów wchodzących w ich skład jest łatwo:

Warto jeszcze sprawdzić czy nie ma przesunięcia, a potem można np. w centroidzie każdego budynku sprawdzać wartości i pokombinować na ile to dokładnie działa na małej próbce.

Samo wrzucanie potem tagów z wysokością powinno być technicznie proste, np.:

https://help.openstreetmap.org/questions/21524/how-to-mass-edit-nodes-with-specific-tags-in-josm?page=1&focusedAnswerId=38850#38850

Tu są ogólne uwagi dotyczące masowych edycji:

https://wiki.openstreetmap.org/wiki/Automated_edits

kocio · 2018-08-23T10:32:49.000Z

Tu jest ciekawy artykuł Tordanika sprzed 5 lat o generowaniu mapy powierzchni na podstawie próbek wysokości - akurat dla celów renderowania 3D, ale tu się też może przydać do wyliczania wysokości gdy teren nie jest idealnie płaski:

http://tobias-knerr.de/publications/thesis/Master%20Thesis%20Knerr.pdf

Lynx86 · 2018-08-23T11:24:58.000Z

Ok, dzięki zapoznam się z treściami i dam znać
Bo etap przypisania wysokości (centroida / mediana / moda / średnia poligonu) jest faktycznie bardzo prosty.
Do napisania później i dzięki za odpowiedzi!

Dajcie tylko proszę znać jaka waszym zdaniem wysokość obiektu byłaby optymalna (komin lub kalenica (najwyższy punkt), średnia z obrysu, czy najczęstsza wartość z obrysu).

kocio · 2018-08-23T12:29:27.000Z

Zdecydowanie ta najwyższa:

https://wiki.openstreetmap.org/wiki/Key:height#Height_of_buildings

możesz od razu dodawać tag source:height=*, żeby było łatwiej to weryfikować (np. ktoś kiedyś przyjdzie z miernikiem laserowym i poprawi poszczególne budynki).

Daj znać jak poszły próby na niewielkiej populacji - dopiero z tego można wyciągnąć jakieś wnioski jak to sensownie zrobić hurtem.

Lynx86 · 2018-08-23T13:14:14.000Z

Profil przez teren, atrybuty z pobranych danych + tabela atrybutów z wartością wysokości maksymalnej

Tu nie w kolorach, ale B&W - widoczne nadbudówki, kominy

z podkładem WMTS ortofoto geoportal / transparency 50%

Te same obiekty widoczne na widoku 3D google maps.
https://www.google.com/maps/@50.0575594,19.9812399,123a,35y,313.1h,47.76t/data=!3m1!1e3

Wygenerowanie centroidy dla każdego budynku też mam zrobione więc jej uzyskanie nie jest problemem.

kocio · 2018-08-23T13:34:06.000Z

Fajnie to wygląda i sensownie. Jak już dodasz to podrzuć link, popatrzę sobie w Kendzi3D, a może także w jakimś stylu 3D, jeśli się szybko zaktualizuje.

Lynx86 · 2018-08-23T13:56:35.000Z

hmm… jeśli nie przegryzę się przez overpass’a i to “technicznie proste wrzucenie” które podlinkowaliście, to będę prosił po weekendzie o pomoc z przejściem tej ścieżki, lub podeślę chętnemu SHP’a z centroidami budynków w wysokościami, które ktoś może spróbować załadować.
Pozdrawiam!

Gdyby ktoś jednak wcześniej chciał: paczka ZIP z centoidami dla kilkudziesięciu budynków
https://drive.google.com/open?id=1ctnNVK71u0_1a24xeLZzfK2GD1GMPTYh

kocio · 2018-08-23T14:06:39.000Z

Ten link do Overpassa zawiera w komentarzu kod, który się wrzuca do kreatora - on już sobie to rozwija. Zasadniczo ustawiasz okienko na interesujący cię obszar, w kreatorze podajesz tylko “building=*” (czyli wybierz dowolne budynki w tym obszarze), a następnie uruchamiasz i dostajesz dane do obróbki. Ja akurat ograniczyłem do Złotna, ale to tylko żeby pokazać że da się wybrać precyzyjnie.

Dawaj znać jak ci idzie.

RicoElectrico · 2018-08-23T14:12:36.000Z

Lynx86:

Wygenerowanie centroidy dla każdego budynku też mam zrobione więc jej uzyskanie nie jest problemem.

Jak sobie radzisz z budynkami, których centroid nie leży w środku (np. w kształcie L lub O)?

Lynx86 · 2018-08-23T14:48:00.000Z

Hm… w zamyśle myślałem że da się uniknąć centroid - realizowałem analizę w QGIS poprzez zonal statistics →
https://gis.stackexchange.com/questions/254974/how-does-zonal-statistics-work-exactly z wartością max
i na podstawie tej funkcji wyciągałem wysokość, gdzie wykorzystywałem cały obrys budynku.

Natomiast (o czym nie pomyślałem przy eksporcie dla Was teraz) faktycznie przy nietypowym kształcie jak L czy S może punkt uciekać poza obrys.

W moim zamyśle - nie wiem czy słusznym - kluczem do połączenia danych powinno być osm_id (takie skrzywienie z systemów bazodanowych) - gdzie jest ono o ile rozumiem unikalne i można na jego podstawie dodać ‘height’ obiektu do obiektu (jakiś join / spatial-join). Jak dla mnie nic nie stoi na przeszkodzie eksport poligonów https://drive.google.com/open?id=1JQcQgWc5dLXjLvB7JtroMxQJ-wZzOPm_ które jedynie trzeba byłoby wyczyścić z miejsc gdzie powstały nowe obiekty (nowe osiedla których nie ma w danych z pokrycia terenu - ale to jest nie aż taka ilość pracy.

kocio · 2018-08-24T10:24:57.000Z

Faktycznie osm_id jest unikalne i też bym po tym identyfikował obiekty, z tym, że nie masz bezpośredniej komunikacji z bazą, tylko przez API, które przyjmuje dane XML, stąd konieczność używania narzędzi do edycji XML zamiast SQL:

https://wiki.openstreetmap.org/wiki/API_v0.6

Ludwik · 2018-08-27T16:49:15.000Z

Zerknijcie jeszcze na to: https://www.geoportal.gov.pl/o-geoportalu/aktualnosci/-/asset_publisher/HCHq0YGNRszn/content/27-08-2018-modele-budynkow-3d-w-geoportalu-dostepne-bezplatnie-do-powszechnego-wykorzystania W komunikacie napisali: “Modele budynków są dostępne bezpłatnie do powszechnego wykorzystania.” - sądzicie, że można traktować to jako zgodę na wykorzystanie w OSM?

szydzio · 2018-08-27T20:24:35.000Z

Ludwik:

sądzicie, że można traktować to jako zgodę na wykorzystanie w OSM?

Ja bym powiedział, że tak, ale jeśli chcesz ostatecznie rozwiać wątpliwości to :

“W przypadku jakichkolwiek pytań prosimy o kontakt z Wydziałem Fotogrametrii i Teledetekcji GUGiK odpowiedzialnym za produkcję przedmiotowych modeli 3D budynków”

RicoElectrico · 2018-08-27T20:34:32.000Z

szydzio:

“W przypadku jakichkolwiek pytań prosimy o kontakt z Wydziałem Fotogrametrii i Teledetekcji GUGiK odpowiedzialnym za produkcję przedmiotowych modeli 3D budynków”

Trzeba to pytanie odpowiednio sformułować, bo jeśli coś im powiesz o ODbL, to pewnie zgłupieją.

kocio · 2018-08-28T10:39:53.000Z

O, super! Ja też rozumiem, że to jest jasna zgoda także na wykorzystanie w OSM.

Mateusz_Konieczny · 2018-09-08T21:02:33.000Z

Lynx86:

Czy jest możliwość edycji “automatycznej” lub “pół-automatycznej” kilkudziesięciu tysięcy obiektów (a może nawet kilkuset tysięcy?).
Już tłumaczę: poprzez analizy teledetekcyjne udało mi się dokonać pomiaru licznych obiektów budowlanych (np. budynki w Krakowie).

Jest taka możliwość jednak na pewno należy zacżąć od czegoś na dużo mniejszą skalę by mieć pewność że nie pojawaiają się jakieś niespodziewane problemy.

Zmianowy · 2018-09-11T13:42:20.000Z

Jakby kto mógł, to hurtowo przydałoby się pozmieniać nazwę Bank Zachodni WBK na Santander.

szydzio · 2018-09-11T19:33:48.000Z

Zmianowy:

Jakby kto mógł, to hurtowo przydałoby się pozmieniać nazwę Bank Zachodni WBK na Santander.

Samą podmianę tagu name dla obiektów z wartościami np. BZ WBK (186 obiektów) lub Bank Zachodni WBK (105 obiektów) mogę zrobić od ręki, pytanie tylko na jaką wartość ?
Oficjalny komunikat na stronie banku głosi, iż „7 września 2018 roku Bank Zachodni WBK S.A. zmienił nazwę na Santander Bank Polska S.A.“ a baner na głównej stronie że „Bank Zachodni WBK to teraz Santander Bank Polska“.
Zatem chyba raczej „Santander Bank Polska“ a nie samo „Santander“ ?
Pytanie też co z obiektami, które mają w nazwie jednocześnie też oddział - mamy pewność, że po zmianie marki oddziały pozostaną tak samo numerowane ?
Czy po prostu wszystko jak leci zmienić tylko na Santander Bank Polska ?

Zmianowy · 2018-09-11T21:51:31.000Z

Jak szukam banku, wpisuję miasto i “Santander” i znajduję. Mnie to wystarczy i tak też już raz zmieniłem nazwę, a jak komuś mało, niech sobie dopisze więcej.

rmikke · 2018-09-12T09:54:16.000Z

IMHO name=“Santander Bank Polska”, w short_name=“Santander” (żeby się wyszukiwało po samym “Santander”), a z numerami oddziałów to bym założył, że reset. Lepiej żeby się dawało łatwo znaleźć po nazwie, niż się wstrzymywać, aż ktoś posprawdza na miejscu.

W drugiej kolejności ręcznie posprawdzać na stronie Santandera który oddział gdzie i pouzupełniać name=* żeby zawierało oddział.