Visualizar ativos de cubo de dados do Computador Planetário Pro da Microsoft

O Computador Planetário Pro da Microsoft inclui uma visualização em blocos que pode ser usada para visualizar alguns ativos do NetCDF.

Verificar a visualização do NetCDF

Nem todos os conjuntos de dados do NetCDF que podem ser ingeridos no Computador Planetário da Microsoft são compatíveis com a visualização em bloco do Computador Planetário Pro. Um conjunto de dados deve ter eixos X e Y, coordenadas de latitude e longitude e dimensões espaciais e limites a serem visualizados. Por exemplo, um conjunto de dados no qual latitude e longitude são variáveis, mas não coordenadas, não é compatível com o renderizador de mapas do Planetary Computer Pro.

Antes de tentar visualizar o conjunto de dados do NetCDF, você pode usar o seguinte para verificar se ele atende aos requisitos.

1. Instalar as dependências necessárias

   pip install xarray[io] rioxarray cf_xarray
   

2. Execute a seguinte função:

   import xarray as xr
   import cf_xarray
   import rioxarray
   
   def is_dataset_visualizable(ds: xr.Dataset):
       """
       Test if the dataset is compatible with the Planetary Computer tiler API.
       Raises an informative error if the dataset is not compatible.
       """
       if not ds.cf.axes:
           raise ValueError("Dataset does not have CF axes")
       if not ds.cf.coordinates:
           raise ValueError("Dataset does not have CF coordinates")
       if not {"X", "Y"} <= ds.cf.axes.keys():
           raise ValueError(f"Dataset must have CF X and Y axes, found: {ds.cf.axes.keys()}")
   
       if not {"latitude", "longitude"} <= ds.cf.coordinates.keys():
           raise ValueError("Dataset must have CF latitude and longitude coordinates, "
                            f"actual: {ds.cf.coordinates.keys()}")
   
       if ds.rio.x_dim is None or ds.rio.y_dim is None:
           raise ValueError("Dataset does not have rioxarray spatial dimensions")
   
       if ds.rio.bounds() is None:
           raise ValueError("Dataset does not have rioxarray bounds")
   
       left, bottom, right, top = ds.rio.bounds()
       if left < -180 or right > 180 or bottom < -90 or top > 90:
           raise ValueError("Dataset bounds are not valid; they must be within [-180, 180] and [-90, 90]")
   
       if ds.rio.resolution() is None:
           raise ValueError("Dataset does not have rioxarray resolution")
   
       if ds.rio.transform() is None:
           raise ValueError("Dataset does not have rioxarray transform")
   
       print("✅ Dataset is compatible with the Planetary Computer tiler API.")