Re: [gdal-dev] gdal_calc.py IF ELSE Condition

2016-07-01 Thread Gabriele N.
Ciao Antonio. Now everything looks perfect. Thank you very much Gabriele -- View this message in context: http://osgeo-org.1560.x6.nabble.com/gdal-calc-py-IF-ELSE-Condition-tp5274166p5274365.html Sent from the GDAL - Dev mailing list archive at Nabble.com. __

Re: [gdal-dev] gdal_calc.py IF ELSE Condition

2016-06-30 Thread Antonio Falciano
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 r

Re: [gdal-dev] gdal_calc.py IF ELSE Condition

2016-06-30 Thread Gabriele N.
Ciao Antonio. @arctan I think is ok. It works but there are many negative values...I must investigate :) Meanwhile thanks Gabrlele -- View this message in context: http://osgeo-org.1560.x6.nabble.com/gdal-calc-py-IF-ELSE-Condition-tp5274166p5274178.html Sent from the GDAL - Dev mailing

Re: [gdal-dev] gdal_calc.py IF ELSE Condition

2016-06-30 Thread Antonio Falciano
Il 30/06/2016 15:25, Gabriele N. ha scritto: Hi all, my goal is to apply the following condition: if map A > 0 then MAP C = arctan((A/B)+180 ELSE MAP C = 90-arctan(A/B) I tried in several ways, including this: gdal_calc.py -A u-wind$DATE.tif -B v-wind$DATE.tif --outfile=dir_wind_$DATE.tif -

[gdal-dev] gdal_calc.py IF ELSE Condition

2016-06-30 Thread Gabriele N.
Hi all, my goal is to apply the following condition: if map A > 0 then MAP C = arctan((A/B)+180 ELSE MAP C = 90-arctan(A/B) I tried in several ways, including this: gdal_calc.py -A u-wind$DATE.tif -B v-wind$DATE.tif --outfile=dir_wind_$DATE.tif --calc="arctan(A/B)+180*(A>0) w_dir=90-arctan(A/