Я пытаюсь преобразовать файлы netCDF в EPSG:3857 для использования с Mapbox с помощью GDAL. Это будет преобразование .nc в .nc. Не до растра. Я открыт для использования GDAL или других методов для этого. Эти данные должны быть перепроецированы, прежде чем они попадут в консольное приложение — и этот процесс занимает недели, чтобы найти решение для — я полагал, что это было просто.
Я работаю над раскрашиванием спутниковых данных. Есть 3 файла .nc (синий, красный и инфракрасный), которые при объединении и обработке создают цветное изображение. После того, как 3 файла загружены (из Amazon AWS), консольное приложение Python выполняет обработку и выгружает .jpg в ту же папку. Исходный код этого приложения находится здесь, чтобы вы могли проверить данные. (Это медленно, так как файлы имеют сверхвысокое разрешение).
Код, который я пробовал:
gdalwarp -t_srs EPSG:3857 test.nc test-projected.nc
Тем не менее, было испробовано несколько других вариантов, и ничего не работает.
Я не профессионал в этом, но должен ли я вообще использовать для этого gdalwarp? Я хочу только изменить проекцию — больше ничего, поэтому приложение Python все еще может работать с данными. Он должен иметь возможность создавать .jpg с использованием перепроецированных файлов.
Следующие ссылки являются примерами данных, которые необходимо преобразовать:
Цветной канал 1 (разрешение синего – 1 км)
Color Channel 2 (красный, более высокое разрешение 0,5 км и больший размер файла)
файл"rel=".nc" на AWS > Цветовой канал 3 (инфракрасный — зеленый)
Кроме того, кто-то еще в Интернете сделал это, используя аналогичную проекцию через модуль pyproj по адресу https://github.com/blaylockbk/pyBKB_v2/tree/master/BB_GOES16. (Мой должен быть EPSG: 3857 для использования с Mapbox). Если бы код Python был изменен, чтобы делать все это за один раз, это тоже было бы здорово. Я открываю награду как последнюю надежду.
Я не знаю Python, поэтому по большей части пытался использовать GDAL, однако работающий код Python, добавленный к моему исходному коду для достижения ожидаемого результата (или работающий скрипт GDAL), принесет награду.