If there are many geotiff files need to be converted to netCDF files, how can we do?
The easy way is to use GDAL.
We can write a simple script to implement this task.
suffix=".tif"
for i in $(ls *.tif); do
# remove suffix:
foo=${i%$suffix}
gdal_translate -of netCDF -co "FORMAT=NC4" $i $foo.nc
done
When I want to merge a lot of daily data files to a single one, sometimes cdo will tell me：
cdo: Argument list too long
Because there are so many files if ten years requested.A solution is to combine files for each year, then combine all years.We need to use a bash for loop to implement:
#!/bin/bash
rm -rf flx.nc
rm -rf slv.nc
year_start=2001
year_end=2009
for year in $(eval echo "{$year_start..$year_end}")
do
echo $year
cdo cat flx/*$year*.nc4 flx$year.nc
cdo cat slv/*$year*.nc4 slv$year.nc
done
cdo cat flx*.nc flx.nc
cdo cat slv*.nc slv.nc
#rm -rf flx*.nc
#rm -rf slv*.nc