Conversión de raster de netCDF ORCA025 a proyección equirectangular

Question

Necesito utilizar los datos ráster de Copernicus Marine Service, que vienen en formato netCDF, y con la cuadrícula tripolar ORCA025 (1440 x 1021 puntos de cuadrícula, los tres polos están situados sobre la Antártida, Asia Central y el Norte de Canadá).

Pero lo que necesito son simples rásters usando latitudes y longitudes regulares (creo que se llama proyección equirectangular, WSG84...) No tengo un software GIS (e incluso, no tengo suficientes conocimientos para usarlos), sólo R.

Así que supongo que la conversión se podría hacer con R y quizás con el paquete gdal, pero no tengo ni idea de cómo hacerlo.

Answer 1

Puede ejecutar gdalinfo en el archivo netcdf para ver que tiene los siguientes subconjuntos de datos:

SUBDATASET_1_NAME=NETCDF:"BIOMER_FREEGLORYS2V3_20161215_R20170703_gridCHL.nc":nav_lon
SUBDATASET_1_DESC=[1021x1442] longitude (32-bit floating-point)
SUBDATASET_2_NAME=NETCDF:"BIOMER_FREEGLORYS2V3_20161215_R20170703_gridCHL.nc":nav_lat
SUBDATASET_2_DESC=[1021x1442] latitude (32-bit floating-point)
SUBDATASET_3_NAME=NETCDF:"BIOMER_FREEGLORYS2V3_20161215_R20170703_gridCHL.nc":CHL
SUBDATASET_3_DESC=[1x75x1021x1442] mass_concentration_of_chlorophyll_in_sea_water (32-bit floating-point)

Por lo tanto, gdalwarp puede tomar la información de lon y lat de los dos primeros subconjuntos de datos, y almacenar el resultado en un raster de +/-180/90 grados con una resolución de un cuarto de grado utilizando esta línea de comando:

gdalwarp -overwrite -geoloc -t_srs EPSG:4326 -tr 0.25 0.25 NETCDF:"BIOMER_FREEGLORYS2V3_20161215_R20170703_gridCHL.nc":CHL out.tif

El resultado se ve bastante bien, excepto por el borde del archivo de origen en el Océano Índico:

ACTUALIZACIÓN

Puede añadir -wo SOURCE_EXTRA=200 al comando gdalwarp para reducir el artefacto, pero eso lleva más tiempo de cálculo y no lo elimina por completo.

Answer 2

Bueno, finalmente lo hicimos usando CDO, como nos sugirió el soporte de Copernicus, porque no fuimos capaces de resolver los problemas detectados usando GDAL. Para Windows requiere instalar Cygwin primero. Esperando que pueda ayudar a otras personas, hicimos un PDF con los pasos que nos funcionaron, y un archivo de descripción de la rejilla de ejemplo utilizado con CDO. Se puede encontrar en el foro de usuarios de Copernicus: http://forum.marine.copernicus.eu/discussion/506/raster-conversion-from-orca025-netcdf-to-esri-ascii#Item_1 .

Esencialmente requiere:

• Instalación de Cygwin para Windows con los paquetes: netcd , udunits , libproj,gcc-core , gcc-g++
• Ir al sitio web de CDO https://code.mpimet.mpg.de/projects/cdo y descargar la última versión etiquetada como cygwin64.
• Añadir cdo.exe a cygwin \bin carpeta.
• utilice este comando en el terminal cywin:

cdo remapnn,grid.txt input_file.nc output_file.nc