Извлечение данных об осадках CRU из netcdf в отдельные файлы GeoTIFF

Я хотел бы извлечь данные об осадках CRU в формате netcdf в отдельные файлы GeoTIFF. Обычно, если файл netcdf содержит только переменные: lon, lat, time и pre, мне удается извлечь его, используя следующий скрипт:

for t in `cdo showdate input.nc`; do
  cdo seldate,$t input.nc dummy.nc
  gdal_translate -of GTiff -a_ullr <top_left_lon> <top_left_lat> <bottom_right_lon> <bottom_right_lat> -a_srs EPSG:4326 dummy.nc $t.tif
done

Данные об осадках CRU имеют переменные: lon, lat, time, pre и stn.

Я не могу использовать приведенный выше сценарий, потому что он содержит 2 поднабора данных, получено сообщение от CDO: Input file contains subdatasets. Please, select one of them for reading.

Как выбрать pre переменные в CDO и применить в приведенном выше скрипте?


cdo-climate
person user97103    schedule 03.09.2020    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Если вы имеете в виду, что файлы имеют более одной переменной, вы можете выбрать переменную pre с помощью команды selvar, которую затем можно передать по каналу seldate:

cdo seldate,$t -selvar,pre input.nc dummy.nc
person Adrian Tompkins    schedule 15.09.2020
comment
Спасибо! Это сработало. Я воспользуюсь вашим предложением для следующего анализа. Некоторое время назад, прежде чем найти ваш путь, я использовал cdo delname,stn input.nc output.nc, удалив переменную stn, а затем выполнил, используя мой сценарий выше, работал, но неэффективно :) - person user97103; 16.09.2020