Forum QGIS

Czy istnieje możliwość liczenia objętości na pdst. NMT albo warstwic (poligonów) w Qgis ?

Nie wiem jak bezpośrednio w QGISie, ale może rzuć okiem na operacje na 3d rastrach, czyli "voxelach", modułami r3.*, https://grass.osgeo.org/grass70/manuals/...intro.html, w GRASSie. r3.univar, r3.mapcalc powinny policzyć co trzeba. Sposobów na import na pewno jest parę, zależnie od rodzaju i formatu danych wejściowych. Jak to rastrowy NMT to np. r3.in.xyz albo r.in.gdal; jak TIN to też się, ale może nie wprost. Przepraszam za dygresję. Ewidentnie mi się nudzi.

Dzię☺ki,

Saga Gis daje radę.

W QGISie w kalkulatorze pól warstwy poligonowej wystarczy stworzyć wirtualne pole nazwane np. powierzchnia i w polu 'wyrażenie' wpisać $area.
Dzięki temu QGIS będzie na 'żywo' wyliczał (i wyświetlał np. jako etykiety) pole powierzchni danego poligonu.
Mam nadzieje, że o to chodziło.
Pozdrawiam!

(05-08-2016, 07:20)sebaq napisał(a): [ -> ]Mam nadzieje, że o to chodziło

Raczej nie, chodziło mu o objętość: oprócz Sagi, to się robi natychmiast za pomocą R.VOLUME od menu Processing...

Ehh ale zonk
Sorrki za niedopatrzenie...

(05-08-2016, 07:49)Italiano napisał(a): [ -> ]Raczej nie, chodziło mu o objętość: oprócz Sagi, to się robi natychmiast za pomocą R.VOLUME od menu Processing...

W jaki sposób liczysz objętość za pomoca r.volume?

@Mazur: Właściwie objętość czego chcesz policzyć?

Proponuję jeszcze następujący "włoski trik" by obliczyć objętości DEM-a.

W Kalkulatorze Rastra, wystarczy wprowadzić coś takiego:

Nazwa_NMT * powierzchnia_pixel

niczym "całki" całego modelu.

W celu oceniania dokładności mojej metody, ściągałem obszar Tatr z zasobów SRTM (rozdzielczość 3"), i wycinałem go w promieniu 10 km. wokół Giewontu.

Oto wyniki końcowe:

R.VOLUME => 405717325197 m³

PROCEDURA ITALIANA => 405717628573 m³

czyli błąd procentowy wynosi 0.9999992522484146, całkiem nieźle...

Załączam NMT (my nazwamy go DEM = Digital Elevation Model) od w/w Giewontu.

(05-08-2016, 11:26)msieczka napisał(a): [ -> ]W jaki sposób liczysz objętość za pomoca r.volume?

Przeczytam ją w raporcie, pod "Total volume = ..."