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.
__
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
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/
