s2_test.py

import datacube
from datacube.storage import masking
import matplotlib.pyplot as plt
import numpy as np

dc = datacube.Datacube(app='load-data-example')

query = {
        'lat': (-35.27, -35.33),
        'lon': (149.07, 149.15),
        'output_crs': 'EPSG:3577',
        'resolution': (-20, 20),
        'time': ('2018-01-01', '2018-02-01')
        }

Canberra = dc.load(product='s2a_ard_granule', group_by='solar_day', **query)

Canberra = masking.mask_invalid_data(Canberra)

plt.imsave("red.png", Canberra['nbar_red'][0,:,:].data)
plt.imsave("mask.png", Canberra['fmask'][0,:,:].data)