Преобразование собственного NetCDF в Gridded NetCDF3

Question

Попытка получить доступ к необработанным данным NEXRAD в собственном формате NetCDF с серверов NOAA (например, thredds-aws, Amazon и т. Д.) И автоматически преобразовать их в файл NetCDF с сеткой для последующей обработки.

Инструментарий NOAA Weather and Climate может экспортировать файл NetCDF3 с сеткой (в папке раскрывающегося списка ниже).

Я обучил внешний код для обработки экспортированных файлов с сеткой. Но чтобы автоматизировать этот процесс и проанализировать несколько лет своего исследования, я должен выяснить, как преобразовать собственный файл с серверов NOAA в этот формат с сеткой. В частности, я хотел бы иметь массив MxN в конечном выходном файле для «отражательной способности», в дБз, самого низкого угла сканирования возвышения, ~ 0,5 градуса, с массивами 1D размера M и N для широты и долготы соответственно.

Чтобы помочь мне понять, что я ищу, здесь есть папка с двумя файлами. Один, который называется «KDOX_Gridded_V06_20170601_190648.nc», - это формат, который я ищу, а другой - RAW, собственный NetCDF, который легко доступен через онлайн-серверы: https://www.dropbox.com/sh/zgpb318ncpif894/AAAuBLsjC0_BRfB0T7peJk9ca?dl=0

Спасибо

Answer 1

Я согласен с DopplerShift, что вам, скорее всего, придется создавать свой собственный код, чтобы получить желаемый результат.

А именно, во входном файле у вас есть Reflectivity (scanR, radialR, gateR) с размерами Reflectivity (3,360,1332), но в выходной Reflectivity (время, широта, долгота) и размерами Reflectivity (1800,497). Я подозреваю, что ни один инструмент не может увидеть всего за один-два шага, что ему нужно для преобразования входных данных из начальных измерений в совершенно другой размер (и, конечно, имена).

По сути, вы могли бы просто начать с:

• вырезание нужной вам переменной

• переименование размеров / координатных переменных

• интерполяция от начальных координат до новых

Answer 2

Я думаю, что PyArt может иметь то, что вы хотите в его map_to_grid функциональности. Этот пример , в частности, кажется, имеет отношение к тому, что вы хотите сделать. После создания сетки вам, вероятно, придется написать некоторый код для ручного создания искомой структуры netCDF, если только она не соответствует модели данных PyArt с сеткой.