[QGIS-it-user] Digest di QGIS-it-user, Volume 66, Numero 18

[Andrea Giudiceandrea]

natan_b wrote
> Given DEM (raster) and Linea (SHP) calculate the slope on each segment of
> line1) Explode Linea--->Segment 2) Field calculator define two float field
> Z_end_point
> raster_value('Lake_100m',1,make_point(x(end_point($geometry)),y(end_point($geometry))))Z_start_point
> raster_value('Lake_100m',1,make_point(x(start_point($geometry)),y(start_point($geometry))))then
> slope= Z_end_point-Z_start_point/lenght($geometry)*100

Ciao natanb,
probabilmente l'ultima espressione è errata.

Infatti oltre ad esserci un refuso ('lenght' al posto di 'length'), mancano
anche delle parentesi.

A mio avviso, l'espressione corretta dovrebbe essere:
slope = (Z_end_point-Z_start_point)/length($geometry)*100

Inoltre bisogna specificare che per avere risultati corretti con queste
formule, sia il layer vettoriale lineare del percorso, sia il layer raster
del DEM devono avere un sistema di riferimento in coordinate piane con unità
di misura in metri (o comunque uguale a quella in cui è espressa
l'altitudine), possibilmente di tipo UTM (o comunque con una proiezione che
sia il più possibile equidistante).

A presto.

Andrea



--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html