from eoreader.reader import Reader
from eoreader.bands import RED, GREEN, BLUE, NDVI, VV, VH

reader = Reader()
# prod = reader.open("your_satellite_product_path")
# bands = prod.load([NDVI, GREEN, CLOUDS])
# stack = prod.stack([RED, GREEN, BLUE], stack_path="output.tif")

print("Available optical bands:", [RED, GREEN, BLUE, NDVI])
print("Available SAR bands:", [VV, VH])