Buongiorno a tutti, ho modificato lo script per calcolare la pendenza media dell'asta con la formula di Taylor-Schwarts. Quindi, al parametro mainchannel_slope ho aggiunto il parametro mainchannel_slope_taylor. Il progamma mi restituisce errore. Riporto la parte del codice con le tre righe aggiunte:
------------------------------------- for i in range(0,len(vertex)-1,2): dist = math.sqrt(math.fabs((nodi[i,0] - nodi[i+1,0]))**2 + math.fabs((nodi[i,1] - nodi[i+1,1]))**2) deltaz = math.fabs(nodi[i,2] - nodi[i+1,2]) # Control to prevent float division by zero (dist=0) try: pendenza = deltaz / dist pendenze.append(pendenza) mainchannel_slope = sum(pendenze) / len(pendenze) * 100 pendenza_taylor = dist /math.sqrt(pendenza) pendenze_taylor.append(pendenza_taylor) mainchannel_slope_taylor = (mainchannel / sum(pendenze_taylor)) ** 2 except : pass ----------------------------------- Il problema è che il alcuni punti ci sono pendenze negative; infatti ponendo la pendenza in valore assoluto il programma funziona, ma il risultato è una pendenza uguale a zero. _______________________________________________ GRASS-Italia mailing list GRASS-Italia@listserv.unipr.it http://listserv.unipr.it/mailman/listinfo/grass-italia