Cómo visualizar un conjunto de datos NetCDF

Question

Supongamos que tengo un NetCDFarchivo bastante grande que abro como un conjunto de datos usando xarray.

import xarray

nc_in = r"....\D_P_WGS84_comp.nc"
nc_ds = xarray.open_dataset(nc_in)

Se compone de múltiples variables ráster, cada una con una sola banda, coordenadas x, y y un crs (4326):

nc_ds

Supongamos que quiero visualizar esto de manera muy simple en mi Python Notebook. Realmente no quiero agregar mapas base o leyendas ni nada por el estilo. Intenté simplemente ejecutar show(nc_ds) but throws an errors AttributeError: el objeto 'Dataset' no tiene atributo 'shape'`. También noté que las dimensiones están congeladas. ¿Cuál es la forma más sencilla de visualizar este conjunto de datos o una sola variable del mismo en Python? ¿Existe una modificación o extensión de paquete que haga esto?

Answer 1

Dado que un Datasetpuede contener más de una variable, no sería trivial que el software sepa qué desea que muestre. DataArraySin embargo, si selecciona una de sus variables , se puede mostrar fácilmente (si está en un formato apropiado).

Las variables de a Datasetestán disponibles a través de su data_varsatributo. Esto expone las variables como DataArrays en un diccionario con los nombres de las variables como claves.

Paso a paso:

variables = nc_ds.data_vars
f02_nbr_e0002 = variables['F02_NBR_E0002']
show(f02_nbr_e0002)

Oneliner eficiente:

show(nc_ds.data_vars['F02_NBR_E0002'])