<!DOCTYPE html><html><head><title></title><style type="text/css">p.MsoNormal,p.MsoNoSpacing{margin:0}</style></head><body><div style="font-family:Arial;">Hello,<br></div><div style="font-family:Arial;"><br></div><div style="font-family:Arial;">I've never tried it, but I believe you can describe some computations in a VRT. But gdal_calc can't do this, so you'd need to write the VRT directly as XML. That's probably not much easier than writing some Python or C++ to perform the same computation.<br></div><div style="font-family:Arial;"><br></div><div style="font-family:Arial;">Laurentiu</div><div style="font-family:Arial;"><br></div><div>On Thu, Oct 20, 2022, at 21:06, Jonas Ardö via gdal-dev wrote:<br></div><blockquote type="cite" id="qt" style=""><div><span class="font" style="font-family:"Courier New", Courier, monospace;">Dear gdal-gurus,<br> <br> Trying to use<b> gdal_calc</b> to do the same calculation for tow
366 band (daily data) files and store it in a *.nc file with the
same dimensions as the input files. <br> <br> This works OK (output to 366 tiff files):<br> <span class="colour" style="color:#0000ff;">for i in {1..366} # One band per day output
to tif<br> do<br> gdal_calc.py --calc="exp( (17.625*D) /(243.04+D) ) / exp(
(17.625*T)/(243.04+T) ) * 100" --outfile=Rh_2020_day$i.tif -D
DPTa_dayavg_C_2020.nc -T Ta_dayavg_C_2020.nc --D_band=$i
--T_band=$i<br> done</span><br> <br> <br> This do not work (output to one netcdf file):<br> <span class="colour" style="color:#0000ff;">for i in {1..366} # One band per day output
to tif<br> do<br> gdal_calc.py --calc="exp( (17.625*D) /(243.04+D) ) / exp(
(17.625*T)/(243.04+T) ) * 100" --outfile=Rh_2020.nc -D
DPTa_dayavg_C_2020.nc -T Ta_dayavg_C_2020.nc --D_band=$i
--T_band=$i<br> done<br> <br> <span class="colour" style="color:#000000;">yields:<br> Error! Output exists, but is the wrong geotransform. Use the
--overwrite option to automatically overwrite the existing
file<br> Error! Output exists, but is the wrong geotransform. Use the
--overwrite option to automatically overwrite the existing
file<br> </span><br> <span class="colour" style="color:#000000;">Any suggestions how to make the output
into a 366-band netcdf file?</span><br> </span></span> </div><pre class="qt-moz-signature" cols="72"><span class="font" style="font-family:"Courier New", Courier, monospace;">using on Cygwin:
gdalinfo.exe --version
GDAL 3.5.1, released 2022/06/30
Regards
/Jonas, Lund University, Sweden
</span><br></pre><div>_______________________________________________<br></div><div>gdal-dev mailing list<br></div><div><a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br></div><div><a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></div><div><br></div></blockquote></body></html>