It would be great if aggregation supported dimension labels on both dense and sparse arrays, as it does multi_index.
Currently, the below throws AttributeError: 'Aggregation' object has no attribute 'label_index'
with tiledb.open(array_uri, "r") as A:
q = A.query(attrs='', dims=['dim1'])
# res = q.agg("sum").multi_index[[0,1]]
res = q.agg("sum").label_index(["l"])[["l0", "l1"]]
print(res)