Warstwa shp z wysokościami budynków

[idzi]

Dziękuję za zainteresowanie tematem, ale:
1. Poprawianie geometrii niewiele daje, bo warstwa SHP z GML LoD zawiera nakładające się obiekty, masę zduplikowanych węzłów i niedomknięte poligony. Przy kilkudziesięciu i więcej obiektach to automatycznie się nie da.
2. I tu jest problem, bo nie znalazłem funkcji, która wyciąga z geometrii obiektu wartość "closest vertex z"

Znalazłem rozwiązanie trochę na około wykorzystując specyfikę plików LoD z Geoportalu więc za jakiś czas może nie działać.


1. Pobrać z Geoportalu BDOT10k
2. Wczytać warstwę budynki np. PL.PZGiK.333.0807__OT_BUBD_A
3. jeśli mamy potrzebę to Usunąć zbędne budynki
4. Pobrać budynki 3D z Geoportalu
5. Zaznaczuć te co nas interesują i skopiować
6. Wkleić na warstwę tymczasową
7. Do Wklejonej warstwy tymczasowej dodać pole H
8. Uzupełnić pole H formułą: z(end_point( $geometry ))-z(start_point($geometry))
z(start_point($geometry)) zwraca wartość z pierwszego punktu geometrii i jest to rzędna terenu
z(end_point( $geometry )) zwraca wartość z ostatniego punktu geometrii i tak się składa, że jest to wartość rzędnej dachu
9. Złączyć warstwę z budynkami z BDOT10k z wklejoną warstwą tymczasową polami identyfikującymi budynki zaznaczając do podłączenia pole z wyliczoną wysokością budynku H
10. Zapisać nową warstwę z budynkami i ich wysokością.