GDAL y reproyección de un archivo NetCDF

Question

Tengo un archivo .nc que puede descargar aquí. Esto son imágenes de satélite.

Quiero volver a proyectar esto y guardarlo en un nuevo archivo .nc, sin guardar en raster. Sólo NetCDF a NetCDF. ¿Es esto posible?

Estoy recibiendo un error que dice "ERROR 4: `test.nc' no se reconoce como un formato de archivo compatible".

El código que he probado es :

gdalwarp -t_srs EPSG:3857 test.nc test-projected.nc

No soy un profesional de esto, pero ¿debería incluso usar gdalwarp para hacer esto? (¿tal vez gdaltransform?) Sólo quiero cambiar la proyección - nada más en este momento.

Answer 1

Gdalinfo en sus datos devueltos a continuación:

Warning 1: Recode from UTF-8 to CP_ACP failed with the error: "Invalid argument".
Driver: netCDF/Network Common Data Format
Files: C:\Users\userr\Downloads\OR_ABI-L1b-RadC-M3C01_G16_s20190321802143_e20190321804516_c20190321804558.nc
Size is 5000, 3000
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["unknown",
        DATUM["unknown",
            SPHEROID["Spheroid",6378137,298.2572221]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Geostationary_Satellite"],
    PARAMETER["central_meridian",-75],
    PARAMETER["satellite_height",35786023],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    EXTENSION["PROJ4","+proj=geos +lon_0=-75 +h=35786023 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs  +sweep=x"]]
Origin = (-3627271.328111211769283,1583173.512170388596132)
Pixel Size = (1002.008657743770527,1002.008657743770527)

Geostationary_Satellite ... Sinceramente, no sé nada de esta proyección... Así que intenté gdalwarp dando la cadena proj4 explícitamente a su nc (capa Rad) como abajo:

gdalwarp -of netCDF 
-s_srs "+proj=geos +lon_0=-75 +h=35786023 +x_0=0 +y_0=0
+ellps=GRS80 +units=m +no_defs +sweep=x" 
-t_srs EPSG:3857  -tr 1002.008657743770527 1002.008657743770527 
NETCDF:"OR_ABI-L1b-RadC-M3C01_G16_s20190321802143_e20190321804516_c20190321804558.nc":Rad 
"OR_ABI-L1b-RadC-M3C01_G16_s20190321802143_e20190321804516_c20190321804558_3857.nc"

[Abajo] Ambas imágenes están subyacentes a la imagen de Bing (en QGIS... La imagen de la izquierda se muestra en la proyección sobre la marcha).

---

[Actualización]

1. A raíz de este post empezamos a trabajar en un nuevo ámbito para crear una imagen en color verdadero a partir de tres imágenes de GOES-16 (Banda1, 2 y 3).
2. La idea inicial era reproyectar los archivos NetCDF proporcionados junto con todas las variables (Rad para la radiancia y DQF como valor QC). Luego, con las tres bandas reproyectadas, el software de David debería ser capaz de convertirlas en una imagen Jpeg de color mejorado... Usando la herramienta GDAL pude reproyectar una variable (Rad) pero no dos variables juntas.
3. En su lugar, extraje Rad variable de cada banda y creó el compuesto RGB por SAGA. Había un buen artículo Bah et al. (2018) cómo crear la banda "verde" mezclando las bandas roja, azul y NIR. Siguiendo este documento, he mezclado el 45% de la Banda2 (rojo), el 45% de la Banda1 (azul) y el 10% de la Banda3 (Nir) mediante una calculadora de trama. La imagen compuesta de salida es la siguiente: