search

JuliaDataCubes / YAXArrays.jl

Yet Another XArray-like Julia package
https://juliadatacubes.github.io/YAXArrays.jl/
Other
89 stars 14 forks source link

How to select variables in a Dataset #316

Closed dpabon closed 11 months ago

dpabon commented 11 months ago

I'm trying to select certain variables of a Dataset:

using Pkg
cd("/Net/Groups/BGI/work_3/OEMC/oemc_towers/")
Pkg.activate("/Net/Groups/BGI/work_3/OEMC/oemc_towers/")

using YAXArrays, CairoMakie, YAXArraysToolbox, Zarr, Statistics,  DiskArrays, DiskArrayTools, Revise, DataFrames, CSV, Plots, ProgressMeter, DimensionalData

YAXArrays.YAXdir("/Net/Groups/BGI/scratch/dpabon/YAXA_tmp")

cubeloc = "/Net/Groups/BGI/scratch/zmhamdi/sen2cube"

sen2cube = open_dataset(cubeloc)

site_name = "ES-LMa"

site_cube = sen2cube[site = At(site_name)]

If I do:

site_10m = site_cube[Variable = At(["Sentinel2__L2A_10m__RED", "Sentinel2__L2A_10m__GREEN","Sentinel2__L2A_10m__NIR", "Sentinel2__L2A_10m__BLUE"])]

No variable is selected and site_10m is equal to site_cube

dpabon commented 11 months ago 
site_10m = sen2cube[site = At(site_name), Variable = At(["Sentinel2__L2A_10m__RED", "Sentinel2__L2A_10m__GREEN","Sentinel2__L2A_10m__NIR", "Sentinel2__L2A_10m__BLUE"])]

Also didn't select the variables

dpabon commented 11 months ago 
site_cube[["Sentinel2__L2A_10m__RED", "Sentinel2__L2A_10m__GREEN","Sentinel2__L2A_10m__NIR", "Sentinel2__L2A_10m__BLUE"]]