Il 30/06/2016 15:47, Gabriele N. ha scritto:
Ciao Antonio.
@arctan I think is ok.
It works but there are many negative values...I must investigate :)
Meanwhile .... thanks
Gabriele,
the previous expression has some errors yet: two parenthesis and a
conversion are missing. This should be the right one:
--calc="(degrees(arctan(A/B))+180)*(A>0)+(90-degrees(arctan(A/B)))*(A<=0)"
where the numpy "degrees" function is used to convert from radians to
degrees more elegantly than before. Tested on GDAL 2.1.0.
ciao
Antonio
--
Antonio Falciano
http://www.linkedin.com/in/antoniofalciano
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
