Daha çox

GRASS GIS-də maili torpaq üzərində Həcmi hesablayın

GRASS GIS-də maili torpaq üzərində Həcmi hesablayın


GRASS GIS ilə bir DEM üzərində bəzi həcm hesablamaları aparmağa çalışıram.

Ancaq inanıram ki, r.volume hər şeyi 0-dan ölçür.

Ölçdüyüm xovlu altındakı torpaq meyllidir.

Bunu əks etdirmək üçün həcmi hesablamaq istərdim, amma bunun üçün ən yaxşı yolu bilmirəm.


Torpaq adlanan bir DEM ilə başladığınızı düşünürəmsə:

I. Yığın ətrafındakı haşiyəni rəqəmləşdirin. Ya:

  1. Çoxbucaqlı təpələrlə mətn faylı hazırlayın. Biçim 'standart' adlanır və https://grass.osgeo.org/grass64/manuals/v.in.ascii.html-də təsvir edilmişdir. Bənzər bir şey görünəcək:
B N x1 y1 x2 y2… xN yN x1 y1

burada N - bir sərhəddəki zirvələrin sayı. Quruluş birdən çox sərhəd üçün təkrarlanır. Sonra mətn sənədini idxal edin:v.in.ascii -n formatı = standart giriş = textFile.txt çıxışı = yığın_ sərhəd

  1. Ot6-dan istifadə edirsinizsə, yalnız sərhədləri rəqəmləşdirə bilərsiniz: v.digit -n map = pile_boundary bgcmd = "d.rast map = terrain_surface"
  2. Ot7 üçün v.edit istifadə etməməyimə baxmayaraq kömək edə bilər.

II. Yığını zımbalayaraq və xovlu altındakı zəmini təxmin edərək yeni bir ərazi səthi düzəldin

# rəqəmsallaşdırılmış vektor sərhədini raster sahəyə çevirmək v.to.rast növü = sahə istifadəsi = pişik girişi = qalaq_ sərhəd çıxışı = qalaq_area # yığını yumruqlayın və altındakı reinterpolatı yığın altındakı ərazi səthini təxmin edin. # Məlumat, ot 6 və 7 r.mapcalc sintaksis biraz fərqlidir. r.mapcalc ifadə = "terrain_wo_pile = if (isnull (pile_area), arazi, null ())" r.fillnulls input = arazi_wo_pile output = terrain_below_pile

III. Xovun sahəsini hesablayın

r.mapcalc ifadə = "qalaq_height_above_terrain = arazi-arazi_below_pile" r.volume data = pile_height_above_terrain clump = pile_area

Videoya baxın: GRASS GIS Tutorial - 3. How Add Multispectral Data in GRASS. How to Add raster data in GRASS